201
|
Huang YT, Chen JM, Ho BC, Wu ZY, Kuo RC, Liu PY. Genome Sequencing and Comparative Analysis of Stenotrophomonas acidaminiphila Reveal Evolutionary Insights Into Sulfamethoxazole Resistance. Front Microbiol 2018; 9:1013. [PMID: 29867899 PMCID: PMC5966563 DOI: 10.3389/fmicb.2018.01013] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2018] [Accepted: 04/30/2018] [Indexed: 11/23/2022] Open
Abstract
Stenotrophomonas acidaminiphila is an aerobic, glucose non-fermentative, Gram-negative bacterium that been isolated from various environmental sources, particularly aquatic ecosystems. Although resistance to multiple antimicrobial agents has been reported in S. acidaminiphila, the mechanisms are largely unknown. Here, for the first time, we report the complete genome and antimicrobial resistome analysis of a clinical isolate S. acidaminiphila SUNEO which is resistant to sulfamethoxazole. Comparative analysis among closely related strains identified common and strain-specific genes. In particular, comparison with a sulfamethoxazole-sensitive strain identified a mutation within the sulfonamide-binding site of folP in SUNEO, which may reduce the binding affinity of sulfamethoxazole. Selection pressure analysis indicated folP in SUNEO is under purifying selection, which may be owing to long-term administration of sulfonamide against Stenotrophomonas.
Collapse
Affiliation(s)
- Yao-Ting Huang
- Department of Computer Science and Information Engineering, National Chung Cheng University, Chiayi, Taiwan
| | - Jia-Min Chen
- Department of Computer Science and Information Engineering, National Chung Cheng University, Chiayi, Taiwan
| | - Bing-Ching Ho
- Department of Clinical Laboratory Sciences and Medical Biotechnology, National Taiwan University Hospital, Taipei, Taiwan
| | - Zong-Yen Wu
- DOE Joint Genome Institute, Walnut Creek, CA, United States.,Department of Veterinary Medicine, National Chung Hsing University, Taichung, Taiwan
| | - Rita C Kuo
- DOE Joint Genome Institute, Walnut Creek, CA, United States
| | - Po-Yu Liu
- The Department of Nursing, Shu-Zen Junior College of Medicine and Management, Kaohsiung, Taiwan.,Rong Hsing Research Center for Translational Medicine, College of Life Sciences, National Chung Hsing University, Taichung, Taiwan.,Division of Infectious Diseases, Department of Internal Medicine, Taichung Veterans General Hospital, Taichung, Taiwan
| |
Collapse
|
202
|
Setiawan D, Brender J, Zhang Y. Recent advances in automated protein design and its future challenges. Expert Opin Drug Discov 2018; 13:587-604. [PMID: 29695210 DOI: 10.1080/17460441.2018.1465922] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Abstract
INTRODUCTION Protein function is determined by protein structure which is in turn determined by the corresponding protein sequence. If the rules that cause a protein to adopt a particular structure are understood, it should be possible to refine or even redefine the function of a protein by working backwards from the desired structure to the sequence. Automated protein design attempts to calculate the effects of mutations computationally with the goal of more radical or complex transformations than are accessible by experimental techniques. Areas covered: The authors give a brief overview of the recent methodological advances in computer-aided protein design, showing how methodological choices affect final design and how automated protein design can be used to address problems considered beyond traditional protein engineering, including the creation of novel protein scaffolds for drug development. Also, the authors address specifically the future challenges in the development of automated protein design. Expert opinion: Automated protein design holds potential as a protein engineering technique, particularly in cases where screening by combinatorial mutagenesis is problematic. Considering solubility and immunogenicity issues, automated protein design is initially more likely to make an impact as a research tool for exploring basic biology in drug discovery than in the design of protein biologics.
Collapse
Affiliation(s)
- Dani Setiawan
- a Department of Computational Medicine and Bioinformatics , University of Michigan , Ann Arbor , MI , USA
| | - Jeffrey Brender
- b Radiation Biology Branch , Center for Cancer Research, National Cancer Institute - NIH , Bethesda , MD , USA
| | - Yang Zhang
- a Department of Computational Medicine and Bioinformatics , University of Michigan , Ann Arbor , MI , USA.,c Department of Biological Chemistry , University of Michigan , Ann Arbor , MI , USA
| |
Collapse
|
203
|
AIMOES: Archive information assisted multi-objective evolutionary strategy for ab initio protein structure prediction. Knowl Based Syst 2018. [DOI: 10.1016/j.knosys.2018.01.028] [Citation(s) in RCA: 40] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/02/2023]
|
204
|
Nicoludis JM, Gaudet R. Applications of sequence coevolution in membrane protein biochemistry. BIOCHIMICA ET BIOPHYSICA ACTA. BIOMEMBRANES 2018; 1860:895-908. [PMID: 28993150 PMCID: PMC5807202 DOI: 10.1016/j.bbamem.2017.10.004] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/19/2017] [Revised: 09/28/2017] [Accepted: 10/02/2017] [Indexed: 12/22/2022]
Abstract
Recently, protein sequence coevolution analysis has matured into a predictive powerhouse for protein structure and function. Direct methods, which use global statistical models of sequence coevolution, have enabled the prediction of membrane and disordered protein structures, protein complex architectures, and the functional effects of mutations in proteins. The field of membrane protein biochemistry and structural biology has embraced these computational techniques, which provide functional and structural information in an otherwise experimentally-challenging field. Here we review recent applications of protein sequence coevolution analysis to membrane protein structure and function and highlight the promising directions and future obstacles in these fields. We provide insights and guidelines for membrane protein biochemists who wish to apply sequence coevolution analysis to a given experimental system.
Collapse
Affiliation(s)
- John M Nicoludis
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, MA 02138, United States
| | - Rachelle Gaudet
- Department of Molecular and Cellular Biology, Harvard University, Cambridge, MA, 02138, United States.
| |
Collapse
|
205
|
Manavalan B, Lee J. SVMQA: support-vector-machine-based protein single-model quality assessment. Bioinformatics 2018; 33:2496-2503. [PMID: 28419290 DOI: 10.1093/bioinformatics/btx222] [Citation(s) in RCA: 130] [Impact Index Per Article: 18.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/15/2016] [Accepted: 04/12/2017] [Indexed: 01/03/2023] Open
Abstract
Motivation The accurate ranking of predicted structural models and selecting the best model from a given candidate pool remain as open problems in the field of structural bioinformatics. The quality assessment (QA) methods used to address these problems can be grouped into two categories: consensus methods and single-model methods. Consensus methods in general perform better and attain higher correlation between predicted and true quality measures. However, these methods frequently fail to generate proper quality scores for native-like structures which are distinct from the rest of the pool. Conversely, single-model methods do not suffer from this drawback and are better suited for real-life applications where many models from various sources may not be readily available. Results In this study, we developed a support-vector-machine-based single-model global quality assessment (SVMQA) method. For a given protein model, the SVMQA method predicts TM-score and GDT_TS score based on a feature vector containing statistical potential energy terms and consistency-based terms between the actual structural features (extracted from the three-dimensional coordinates) and predicted values (from primary sequence). We trained SVMQA using CASP8, CASP9 and CASP10 targets and determined the machine parameters by 10-fold cross-validation. We evaluated the performance of our SVMQA method on various benchmarking datasets. Results show that SVMQA outperformed the existing best single-model QA methods both in ranking provided protein models and in selecting the best model from the pool. According to the CASP12 assessment, SVMQA was the best method in selecting good-quality models from decoys in terms of GDTloss. Availability and implementation SVMQA method can be freely downloaded from http://lee.kias.re.kr/SVMQA/SVMQA_eval.tar.gz. Contact jlee@kias.re.kr. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Balachandran Manavalan
- Center for In Silico Protein Science and School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea
| | - Jooyoung Lee
- Center for In Silico Protein Science and School of Computational Sciences, Korea Institute for Advanced Study, Seoul 130-722, Korea
| |
Collapse
|
206
|
Fanuel S, Tabesh S, Mokhtarian K, Saroddiny E, Fazlollahi MR, Pourpak Z, Falak R, Kardar GA. Construction of a recombinant B-cell epitope vaccine based on a Der p1-derived hypoallergen: a bioinformatics approach. Immunotherapy 2018; 10:537-553. [PMID: 29569512 DOI: 10.2217/imt-2017-0163] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/21/2022] Open
Abstract
AIM House dust mite (HDM) allergens are important elicitors of IgE-mediated allergies. This study was aimed at constructing and characterizing a recombinant fusion protein, DpTTDp, which was based on carrier-bound Der p 1-derived peptides for HDM allergen immunotherapy. METHODS Using the Immune Epitope Database (IEDB), we identified from Der p 1, a 34-mer hypoallergenic peptide. Two copies of the hypoallergen were then fused to a partial fragment of a tetanus toxoid molecule's N-and C terminus and expressed in Escherichia coli. After purification to homogeneity, the protein was evaluated for allergenicity and its ability to induce blocking antibodies upon immunization. RESULTS Upon immunization of mice, DpTTDp induced high levels of protective IgG-antibodies that blocked allergic patients' IgE reactivity to HDM. In addition, DpTTDp lacked relevant IgE-reactivity, induced low T-cell proliferation and IFN-γ in peripheral blood mononuclear cells of HDM-allergic patients' sera. CONCLUSION The protein represents a promising HDM-allergy immunotherapy candidate vaccine.
Collapse
Affiliation(s)
- Songwe Fanuel
- Department of Medical Biotechnology, School of Advanced Technologies in Medicine, Tehran University of Medical Sciences - International Campus (IC-TUMS) Tehran, Iran.,Immunology, Asthma & Allergy Research Institute (IAARI), Tehran University of Medical Science, Tehran, Iran
| | - Saeideh Tabesh
- Department of Immunology, School of Public Health, Tehran University of Medical Sciences, Tehran, Iran
| | - Kobra Mokhtarian
- Medicinal Plant Research Center, Shahrekord University of Medical Sciences, Shahrekord, Iran
| | - Esmaeil Saroddiny
- Department of Medical Biotechnology, School of Advanced Technologies in Medicine, Tehran University of Medical Sciences - International Campus (IC-TUMS) Tehran, Iran
| | - Mohammad Reza Fazlollahi
- Immunology, Asthma & Allergy Research Institute (IAARI), Tehran University of Medical Science, Tehran, Iran
| | - Zahra Pourpak
- Immunology, Asthma & Allergy Research Institute (IAARI), Tehran University of Medical Science, Tehran, Iran
| | - Reza Falak
- Department of Immunology, School of Medicine, Iran University of Medical Sciences, Tehran, Iran
| | - Gholam Ali Kardar
- Department of Medical Biotechnology, School of Advanced Technologies in Medicine, Tehran University of Medical Sciences - International Campus (IC-TUMS) Tehran, Iran.,Immunology, Asthma & Allergy Research Institute (IAARI), Tehran University of Medical Science, Tehran, Iran
| |
Collapse
|
207
|
Kato Y, Kihara H, Fukui K, Kojima M. A ternary complex model of Sirtuin4-NAD +-Glutamate dehydrogenase. Comput Biol Chem 2018; 74:94-104. [PMID: 29571013 DOI: 10.1016/j.compbiolchem.2018.03.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/16/2016] [Revised: 11/09/2017] [Accepted: 03/08/2018] [Indexed: 10/17/2022]
Abstract
Sirtuin4 (Sirt4) is one of the mammalian homologues of Silent information regulator 2 (Sir2), which promotes the longevity of yeast, C. elegans, fruit flies and mice. Sirt4 is localized in the mitochondria, where it contributes to preventing the development of cancers and ischemic heart disease through regulating energy metabolism. The ADP-ribosylation of glutamate dehydrogenase (GDH), which is catalyzed by Sirt4, downregulates the TCA cycle. However, this reaction mechanism is obscure, because the structure of Sirt4 is unknown. We here constructed structural models of Sirt4 by homology modeling and threading, and docked nicotinamide adenine dinucleotide+ (NAD+) to Sirt4. In addition, a partial GDH structure was docked to the Sirt4-NAD+ complex model. In the ternary complex model of Sirt4-NAD+-GDH, the acetylated lysine 171 of GDH is located close to NAD+. This suggests a possible mechanism underlying the ADP-ribosylation at cysteine 172, which may occur through a transient intermediate with ADP-ribosylation at the acetylated lysine 171. These results may be useful in designing drugs for the treatment of cancers and ischemic heart disease.
Collapse
Affiliation(s)
- Yusuke Kato
- School of Life Sciences, Tokyo University of Pharmacy and Life Sciences, 1432-1 Horinouchi, Hachioji 192-0392, Japan; Himeji Hinomoto College, 890 Koro, Himeji 679-2151, Japan; Institute for Enzyme Research, Tokushima University, 3-18-15 Kuramoto, Tokushima 770-8503, Japan.
| | - Hiroshi Kihara
- Himeji Hinomoto College, 890 Koro, Himeji 679-2151, Japan
| | - Kiyoshi Fukui
- Institute for Enzyme Research, Tokushima University, 3-18-15 Kuramoto, Tokushima 770-8503, Japan
| | - Masaki Kojima
- School of Life Sciences, Tokyo University of Pharmacy and Life Sciences, 1432-1 Horinouchi, Hachioji 192-0392, Japan
| |
Collapse
|
208
|
Correnti CE, Gewe MM, Mehlin C, Bandaranayake AD, Johnsen WA, Rupert PB, Brusniak MY, Clarke M, Burke SE, De Van Der Schueren W, Pilat K, Turnbaugh SM, May D, Watson A, Chan MK, Bahl CD, Olson JM, Strong RK. Screening, large-scale production and structure-based classification of cystine-dense peptides. Nat Struct Mol Biol 2018; 25:270-278. [PMID: 29483648 PMCID: PMC5840021 DOI: 10.1038/s41594-018-0033-9] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2017] [Accepted: 01/23/2018] [Indexed: 12/04/2022]
Abstract
Peptides folded through interwoven disulfides display extreme biochemical properties and unique medicinal potential. However, their exploitation has been hampered by the limited amounts isolatable from natural sources and the expense of chemical synthesis. We developed reliable biological methods for high-throughput expression, screening and large-scale production of these peptides: 46 were successfully produced in multimilligram quantities, and >600 more were deemed expressible through stringent screening criteria. Many showed extreme resistance to temperature, proteolysis and/or reduction, and all displayed inhibitory activity against at least 1 of 20 ion channels tested, thus confirming their biological functionality. Crystal structures of 12 confirmed proper cystine topology and the utility of crystallography to study these molecules but also highlighted the need for rational classification. Previous categorization attempts have focused on limited subsets featuring distinct motifs. Here we present a global definition, classification and analysis of >700 structures of cystine-dense peptides, providing a unifying framework for these molecules.
Collapse
Affiliation(s)
- Colin E Correnti
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Mesfin M Gewe
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Christopher Mehlin
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Ashok D Bandaranayake
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - William A Johnsen
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Peter B Rupert
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Mi-Youn Brusniak
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Midori Clarke
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Skyler E Burke
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | | | - Kristina Pilat
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Shanon M Turnbaugh
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Damon May
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
- Department of Genome Sciences, University of Washington, Seattle, WA, USA
| | - Alex Watson
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | - Man Kid Chan
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA
| | | | - James M Olson
- Clinical Research Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA.
| | - Roland K Strong
- Basic Sciences Division, Fred Hutchinson Cancer Research Center, Seattle, WA, USA.
| |
Collapse
|
209
|
Nielsen JT, Mulder FAA. POTENCI: prediction of temperature, neighbor and pH-corrected chemical shifts for intrinsically disordered proteins. JOURNAL OF BIOMOLECULAR NMR 2018; 70:141-165. [PMID: 29399725 DOI: 10.1007/s10858-018-0166-5] [Citation(s) in RCA: 113] [Impact Index Per Article: 16.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/29/2017] [Accepted: 01/25/2018] [Indexed: 05/04/2023]
Abstract
Chemical shifts contain important site-specific information on the structure and dynamics of proteins. Deviations from statistical average values, known as random coil chemical shifts (RCCSs), are extensively used to infer these relationships. Unfortunately, the use of imprecise reference RCCSs leads to biased inference and obstructs the detection of subtle structural features. Here we present a new method, POTENCI, for the prediction of RCCSs that outperforms the currently most authoritative methods. POTENCI is parametrized using a large curated database of chemical shifts for protein segments with validated disorder; It takes pH and temperature explicitly into account, and includes sequence-dependent nearest and next-nearest neighbor corrections as well as second-order corrections. RCCS predictions with POTENCI show root-mean-square values that are lower by 25-78%, with the largest improvements observed for 1Hα and 13C'. It is demonstrated how POTENCI can be applied to analyze subtle deviations from RCCSs to detect small populations of residual structure in intrinsically disorder proteins that were not discernible before. POTENCI source code is available for download, or can be deployed from the URL http://www.protein-nmr.org .
Collapse
Affiliation(s)
- Jakob Toudahl Nielsen
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000, Aarhus C, Denmark.
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000, Aarhus C, Denmark.
| | - Frans A A Mulder
- Interdisciplinary Nanoscience Center (iNANO), Aarhus University, Gustav Wieds Vej 14, 8000, Aarhus C, Denmark.
- Department of Chemistry, Aarhus University, Langelandsgade 140, 8000, Aarhus C, Denmark.
| |
Collapse
|
210
|
Moult J, Fidelis K, Kryshtafovych A, Schwede T, Tramontano A. Critical assessment of methods of protein structure prediction (CASP)-Round XII. Proteins 2018; 86 Suppl 1:7-15. [PMID: 29082672 PMCID: PMC5897042 DOI: 10.1002/prot.25415] [Citation(s) in RCA: 219] [Impact Index Per Article: 31.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/21/2017] [Revised: 10/25/2017] [Accepted: 10/27/2017] [Indexed: 12/24/2022]
Abstract
This article reports the outcome of the 12th round of Critical Assessment of Structure Prediction (CASP12), held in 2016. CASP is a community experiment to determine the state of the art in modeling protein structure from amino acid sequence. Participants are provided sequence information and in turn provide protein structure models and related information. Analysis of the submitted structures by independent assessors provides a comprehensive picture of the capabilities of current methods, and allows progress to be identified. This was again an exciting round of CASP, with significant advances in 4 areas: (i) The use of new methods for predicting three-dimensional contacts led to a two-fold improvement in contact accuracy. (ii) As a consequence, model accuracy for proteins where no template was available improved dramatically. (iii) Models based on a structural template showed overall improvement in accuracy. (iv) Methods for estimating the accuracy of a model continued to improve. CASP continued to develop new areas: (i) Assessing methods for building quaternary structure models, including an expansion of the collaboration between CASP and CAPRI. (ii) Modeling with the aid of experimental data was extended to include SAXS data, as well as again using chemical cross-linking information. (iii) A team of assessors evaluated the suitability of models for a range of applications, including mutation interpretation, analysis of ligand binding properties, and identification of interfaces. This article describes the experiment and summarizes the results. The rest of this special issue of PROTEINS contains papers describing CASP12 results and assessments in more detail.
Collapse
Affiliation(s)
- John Moult
- Institute for Bioscience and Biotechnology Research and Department of Cell Biology and Molecular Genetics, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
| | - Krzysztof Fidelis
- Genome Center, University of California, Davis, 451 Health Sciences Drive, Davis, CA 95616, USA
| | - Andriy Kryshtafovych
- Genome Center, University of California, Davis, 451 Health Sciences Drive, Davis, CA 95616, USA
| | - Torsten Schwede
- University of Basel, Biozentrum & SIB Swiss Institute of Bioinformatics, Basel, Switzerland
| | - Anna Tramontano
- Department of Physics and Istituto Pasteur - Fondazione Cenci Bolognetti, Sapienza University of Rome, P.le Aldo Moro, 5, 00185 Rome, Italy
| |
Collapse
|
211
|
Álvarez Ó, Fernández-Martínez JL, Fernández-Brillet C, Cernea A, Fernández-Muñiz Z, Kloczkowski A. Principal component analysis in protein tertiary structure prediction. J Bioinform Comput Biol 2018; 16:1850005. [PMID: 29566640 DOI: 10.1142/s0219720018500051] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
We discuss applicability of principal component analysis (PCA) for protein tertiary structure prediction from amino acid sequence. The algorithm presented in this paper belongs to the category of protein refinement models and involves establishing a low-dimensional space where the sampling (and optimization) is carried out via particle swarm optimizer (PSO). The reduced space is found via PCA performed for a set of low-energy protein models previously found using different optimization techniques. A high frequency term is added into this expansion by projecting the best decoy into the PCA basis set and calculating the residual model. This term is aimed at providing high frequency details in the energy optimization. The goal of this research is to analyze how the dimensionality reduction affects the prediction capability of the PSO procedure. For that purpose, different proteins from the Critical Assessment of Techniques for Protein Structure Prediction experiments were modeled. In all the cases, both the energy of the best decoy and the distance to the native structure have decreased. Our analysis also shows how the predicted backbone structure of native conformation and of alternative low energy states varies with respect to the PCA dimensionality. Generally speaking, the reconstruction can be successfully achieved with 10 principal components and the high frequency term. We also provide a computational analysis of protein energy landscape for the inverse problem of reconstructing structure from the reduced number of principal components, showing that the dimensionality reduction alleviates the ill-posed character of this high-dimensional energy optimization problem. The procedure explained in this paper is very fast and allows testing different PCA expansions. Our results show that PSO improves the energy of the best decoy used in the PCA when the adequate number of PCA terms is considered.
Collapse
Affiliation(s)
- Óscar Álvarez
- * Group of Inverse Problems, Optimization and Machine Learning, Department of Mathematics, University of Oviedo, C. Federico García Lorca, 18, 33007 Oviedo, Spain
| | - Juan Luis Fernández-Martínez
- * Group of Inverse Problems, Optimization and Machine Learning, Department of Mathematics, University of Oviedo, C. Federico García Lorca, 18, 33007 Oviedo, Spain
| | - Celia Fernández-Brillet
- * Group of Inverse Problems, Optimization and Machine Learning, Department of Mathematics, University of Oviedo, C. Federico García Lorca, 18, 33007 Oviedo, Spain
| | - Ana Cernea
- * Group of Inverse Problems, Optimization and Machine Learning, Department of Mathematics, University of Oviedo, C. Federico García Lorca, 18, 33007 Oviedo, Spain
| | - Zulima Fernández-Muñiz
- * Group of Inverse Problems, Optimization and Machine Learning, Department of Mathematics, University of Oviedo, C. Federico García Lorca, 18, 33007 Oviedo, Spain
| | - Andrzej Kloczkowski
- † Batelle Center for Mathematical Medicine, Nationwide Children's Hospital, Columbus, OH, USA.,‡ Department of Pediatrics, The Ohio State University, Columbus, OH, USA
| |
Collapse
|
212
|
Wang X, Zhang D, Huang SY. New Knowledge-Based Scoring Function with Inclusion of Backbone Conformational Entropies from Protein Structures. J Chem Inf Model 2018; 58:724-732. [DOI: 10.1021/acs.jcim.7b00601] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Affiliation(s)
- Xinxiang Wang
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, P. R. China
| | - Di Zhang
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, P. R. China
| | - Sheng-You Huang
- School of Physics, Huazhong University of Science and Technology, Wuhan, Hubei 430074, P. R. China
| |
Collapse
|
213
|
dos Santos RN, Ferrari AJR, de Jesus HCR, Gozzo FC, Morcos F, Martínez L. Enhancing protein fold determination by exploring the complementary information of chemical cross-linking and coevolutionary signals. Bioinformatics 2018; 34:2201-2208. [DOI: 10.1093/bioinformatics/bty074] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/19/2017] [Accepted: 02/10/2018] [Indexed: 11/13/2022] Open
Affiliation(s)
- Ricardo N dos Santos
- Institute of Chemistry, University of Campinas, Campinas, Brazil
- Center for Computational Engineering and Sciences, University of Campinas, Campinas, Brazil
| | | | | | - Fábio C Gozzo
- Institute of Chemistry, University of Campinas, Campinas, Brazil
| | - Faruck Morcos
- Department of Biological Sciences, University of Texas at Dallas, Richardson, USA
| | - Leandro Martínez
- Institute of Chemistry, University of Campinas, Campinas, Brazil
- Center for Computational Engineering and Sciences, University of Campinas, Campinas, Brazil
| |
Collapse
|
214
|
Guardiani C, Fedorenko OA, Roberts SK, Khovanov IA. On the selectivity of the NaChBac channel: an integrated computational and experimental analysis of sodium and calcium permeation. Phys Chem Chem Phys 2018; 19:29840-29854. [PMID: 29090695 DOI: 10.1039/c7cp05928k] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Ion channel selectivity is essential for their function, yet the molecular basis of a channel's ability to select between ions is still rather controversial. In this work, using a combination of molecular dynamics simulations and electrophysiological current measurements we analyze the ability of the NaChBac channel to discriminate between calcium and sodium. Our simulations show that a single calcium ion can access the Selectivity Filter (SF) interacting so strongly with the glutamate ring so as to remain blocked inside. This is consistent with the tiny calcium currents recorded in our patch-clamp experiments. Two reasons explain this scenario. The first is the higher free energy of ion/SF binding of Ca2+ with respect to Na+. The second is the strong electrostatic repulsion exerted by the resident ion that turns back a second potentially incoming Ca2+, preventing the knock-on permeation mechanism. Finally, we analyzed the possibility of the Anomalous Mole Fraction Effect (AMFE), i.e. the ability of micromolar Ca2+ concentrations to block Na+ currents. Current measurements in Na+/Ca2+ mixed solutions excluded the AMFE, in agreement with metadynamics simulations showing the ability of a sodium ion to by-pass and partially displace the resident calcium. Our work supports a new scenario for Na+/Ca2+ selectivity in the bacterial sodium channel, challenging the traditional notion of an exclusion mechanism strictly confining Ca2+ ions outside the channel.
Collapse
Affiliation(s)
- Carlo Guardiani
- School of Engineering, University of Warwick, Coventry CV4 7AL, UK.
| | | | | | | |
Collapse
|
215
|
Schneider J, Korshunova K, Musiani F, Alfonso-Prieto M, Giorgetti A, Carloni P. Predicting ligand binding poses for low-resolution membrane protein models: Perspectives from multiscale simulations. Biochem Biophys Res Commun 2018; 498:366-374. [PMID: 29409902 DOI: 10.1016/j.bbrc.2018.01.160] [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: 12/24/2017] [Revised: 01/24/2018] [Accepted: 01/25/2018] [Indexed: 12/21/2022]
Abstract
Membrane receptors constitute major targets for pharmaceutical intervention. Drug design efforts rely on the identification of ligand binding poses. However, the limited experimental structural information available may make this extremely challenging, especially when only low-resolution homology models are accessible. In these cases, the predictions may be improved by molecular dynamics simulation approaches. Here we review recent developments of multiscale, hybrid molecular mechanics/coarse-grained (MM/CG) methods applied to membrane proteins. In particular, we focus on our in-house MM/CG approach. It is especially tailored for G-protein coupled receptors, the largest membrane receptor family in humans. We show that our MM/CG approach is able to capture the atomistic details of the receptor/ligand binding interactions, while keeping the computational cost low by representing the protein frame and the membrane environment in a highly simplified manner. We close this review by discussing ongoing improvements and challenges of the current implementation of our MM/CG code.
Collapse
Affiliation(s)
- Jakob Schneider
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany; Department of Physics, Rheinisch-Westfälische Technische Hochschule Aachen, Aachen, Germany; JARA Institute Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich GmbH, Jülich, Germany
| | - Ksenia Korshunova
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany; Department of Physics, Rheinisch-Westfälische Technische Hochschule Aachen, Aachen, Germany
| | - Francesco Musiani
- Laboratory of Bioinorganic Chemistry, Department of Pharmacy and Biotechnology, University of Bologna, Bologna, Italy
| | - Mercedes Alfonso-Prieto
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany; Cécile and Oskar Vogt Institute for Brain Research, Medical Faculty, Heinrich Heine University Düsseldorf, Düsseldorf, Germany.
| | - Alejandro Giorgetti
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany; Department of Biotechnology, University of Verona, Verona, Italy
| | - Paolo Carloni
- Computational Biomedicine, Institute for Advanced Simulations IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich GmbH, Jülich, Germany; Department of Physics, Rheinisch-Westfälische Technische Hochschule Aachen, Aachen, Germany; JARA Institute Molecular Neuroscience and Neuroimaging (INM-11), Forschungszentrum Jülich GmbH, Jülich, Germany; VNU Key Laboratory "Multiscale Simulation of Complex Systems", VNU University of Science, Vietnam National University, Hanoi, Viet Nam.
| |
Collapse
|
216
|
Li B, Fooksa M, Heinze S, Meiler J. Finding the needle in the haystack: towards solving the protein-folding problem computationally. Crit Rev Biochem Mol Biol 2018; 53:1-28. [PMID: 28976219 PMCID: PMC6790072 DOI: 10.1080/10409238.2017.1380596] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/16/2017] [Revised: 08/22/2017] [Accepted: 09/13/2017] [Indexed: 12/22/2022]
Abstract
Prediction of protein tertiary structures from amino acid sequence and understanding the mechanisms of how proteins fold, collectively known as "the protein folding problem," has been a grand challenge in molecular biology for over half a century. Theories have been developed that provide us with an unprecedented understanding of protein folding mechanisms. However, computational simulation of protein folding is still difficult, and prediction of protein tertiary structure from amino acid sequence is an unsolved problem. Progress toward a satisfying solution has been slow due to challenges in sampling the vast conformational space and deriving sufficiently accurate energy functions. Nevertheless, several techniques and algorithms have been adopted to overcome these challenges, and the last two decades have seen exciting advances in enhanced sampling algorithms, computational power and tertiary structure prediction methodologies. This review aims at summarizing these computational techniques, specifically conformational sampling algorithms and energy approximations that have been frequently used to study protein-folding mechanisms or to de novo predict protein tertiary structures. We hope that this review can serve as an overview on how the protein-folding problem can be studied computationally and, in cases where experimental approaches are prohibitive, help the researcher choose the most relevant computational approach for the problem at hand. We conclude with a summary of current challenges faced and an outlook on potential future directions.
Collapse
Affiliation(s)
- Bian Li
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Michaela Fooksa
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
- Chemical and Physical Biology Graduate Program, Vanderbilt University, Nashville, TN, USA
| | - Sten Heinze
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| | - Jens Meiler
- Department of Chemistry, Vanderbilt University, Nashville, TN, USA
- Center for Structural Biology, Vanderbilt University, Nashville, TN, USA
| |
Collapse
|
217
|
Antunes D, Jorge NAN, Caffarena ER, Passetti F. Using RNA Sequence and Structure for the Prediction of Riboswitch Aptamer: A Comprehensive Review of Available Software and Tools. Front Genet 2018; 8:231. [PMID: 29403526 PMCID: PMC5780412 DOI: 10.3389/fgene.2017.00231] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Accepted: 12/21/2017] [Indexed: 12/14/2022] Open
Abstract
RNA molecules are essential players in many fundamental biological processes. Prokaryotes and eukaryotes have distinct RNA classes with specific structural features and functional roles. Computational prediction of protein structures is a research field in which high confidence three-dimensional protein models can be proposed based on the sequence alignment between target and templates. However, to date, only a few approaches have been developed for the computational prediction of RNA structures. Similar to proteins, RNA structures may be altered due to the interaction with various ligands, including proteins, other RNAs, and metabolites. A riboswitch is a molecular mechanism, found in the three kingdoms of life, in which the RNA structure is modified by the binding of a metabolite. It can regulate multiple gene expression mechanisms, such as transcription, translation initiation, and mRNA splicing and processing. Due to their nature, these entities also act on the regulation of gene expression and detection of small metabolites and have the potential to helping in the discovery of new classes of antimicrobial agents. In this review, we describe software and web servers currently available for riboswitch aptamer identification and secondary and tertiary structure prediction, including applications.
Collapse
Affiliation(s)
- Deborah Antunes
- Scientific Computing Program (PROCC), Computational Biophysics and Molecular Modeling Group, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Natasha A N Jorge
- Laboratory of Functional Genomics and Bioinformatics, Oswaldo Cruz Institute, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratory of Gene Expression Regulation, Carlos Chagas Institute, Fundação Oswaldo Cruz, Curitiba, Brazil
| | - Ernesto R Caffarena
- Scientific Computing Program (PROCC), Computational Biophysics and Molecular Modeling Group, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil
| | - Fabio Passetti
- Laboratory of Functional Genomics and Bioinformatics, Oswaldo Cruz Institute, Fundação Oswaldo Cruz, Rio de Janeiro, Brazil.,Laboratory of Gene Expression Regulation, Carlos Chagas Institute, Fundação Oswaldo Cruz, Curitiba, Brazil
| |
Collapse
|
218
|
Marks C, Nowak J, Klostermann S, Georges G, Dunbar J, Shi J, Kelm S, Deane CM. Sphinx: merging knowledge-based and ab initio approaches to improve protein loop prediction. Bioinformatics 2018; 33:1346-1353. [PMID: 28453681 PMCID: PMC5408792 DOI: 10.1093/bioinformatics/btw823] [Citation(s) in RCA: 28] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/25/2016] [Accepted: 01/09/2017] [Indexed: 01/31/2023] Open
Abstract
Motivation Loops are often vital for protein function, however, their irregular structures make them difficult to model accurately. Current loop modelling algorithms can mostly be divided into two categories: knowledge-based, where databases of fragments are searched to find suitable conformations and ab initio, where conformations are generated computationally. Existing knowledge-based methods only use fragments that are the same length as the target, even though loops of slightly different lengths may adopt similar conformations. Here, we present a novel method, Sphinx, which combines ab initio techniques with the potential extra structural information contained within loops of a different length to improve structure prediction. Results We show that Sphinx is able to generate high-accuracy predictions and decoy sets enriched with near-native loop conformations, performing better than the ab initio algorithm on which it is based. In addition, it is able to provide predictions for every target, unlike some knowledge-based methods. Sphinx can be used successfully for the difficult problem of antibody H3 prediction, outperforming RosettaAntibody, one of the leading H3-specific ab initio methods, both in accuracy and speed. Availability and Implementation Sphinx is available at http://opig.stats.ox.ac.uk/webapps/sphinx. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Claire Marks
- Department of Statistics, University of Oxford, Oxford, UK
| | - Jaroslaw Nowak
- Department of Statistics, University of Oxford, Oxford, UK
| | | | - Guy Georges
- Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, DE, Germany
| | - James Dunbar
- Pharma Research and Early Development, Large Molecule Research, Roche Innovation Center Munich, Penzberg, DE, Germany
| | - Jiye Shi
- Department of Informatics, UCB Pharma, Slough, UK
| | | | | |
Collapse
|
219
|
Assessing Exhaustiveness of Stochastic Sampling for Integrative Modeling of Macromolecular Structures. Biophys J 2018; 113:2344-2353. [PMID: 29211988 DOI: 10.1016/j.bpj.2017.10.005] [Citation(s) in RCA: 66] [Impact Index Per Article: 9.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2017] [Revised: 09/22/2017] [Accepted: 10/02/2017] [Indexed: 12/22/2022] Open
Abstract
Modeling of macromolecular structures involves structural sampling guided by a scoring function, resulting in an ensemble of good-scoring models. By necessity, the sampling is often stochastic, and must be exhaustive at a precision sufficient for accurate modeling and assessment of model uncertainty. Therefore, the very first step in analyzing the ensemble is an estimation of the highest precision at which the sampling is exhaustive. Here, we present an objective and automated method for this task. As a proxy for sampling exhaustiveness, we evaluate whether two independently and stochastically generated sets of models are sufficiently similar. The protocol includes testing 1) convergence of the model score, 2) whether model scores for the two samples were drawn from the same parent distribution, 3) whether each structural cluster includes models from each sample proportionally to its size, and 4) whether there is sufficient structural similarity between the two model samples in each cluster. The evaluation also provides the sampling precision, defined as the smallest clustering threshold that satisfies the third, most stringent test. We validate the protocol with the aid of enumerated good-scoring models for five illustrative cases of binary protein complexes. Passing the proposed four tests is necessary, but not sufficient for thorough sampling. The protocol is general in nature and can be applied to the stochastic sampling of any set of models, not just structural models. In addition, the tests can be used to stop stochastic sampling as soon as exhaustiveness at desired precision is reached, thereby improving sampling efficiency; they may also help in selecting a model representation that is sufficiently detailed to be informative, yet also sufficiently coarse for sampling to be exhaustive.
Collapse
|
220
|
Slater PG, Gutierrez-Maldonado SE, Gysling K, Lagos CF. Molecular Modeling of Structures and Interaction of Human Corticotropin-Releasing Factor (CRF) Binding Protein and CRF Type-2 Receptor. Front Endocrinol (Lausanne) 2018; 9:43. [PMID: 29515519 PMCID: PMC5826306 DOI: 10.3389/fendo.2018.00043] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The corticotropin-releasing factor (CRF) system is a key mediator of the stress response and addictive behavior. The CRF system includes four peptides: The CRF system includes four peptides: CRF, urocortins I-III, CRF binding protein (CRF-BP) that binds CRF with high affinity, and two class B G-protein coupled receptors CRF1R and CRF2R. CRF-BP is a secreted protein without significant sequence homology to CRF receptors or to any other known class of protein. Recently, it has been described a potentiation role of CRF-BP over CRF signaling through CRF2R in addictive-related neuronal plasticity and behavior. In addition, it has been described that CRF-BP is capable to physically interact specifically with the α isoform of CRF2R and acts like an escort protein increasing the amount of the receptor in the plasma membrane. At present, there are no available structures for CRF-BP or for full-length CRFR. Knowing and studying the structure of these proteins could be beneficial in order to characterize the CRF-BP/CRF2αR interaction. In this work, we report the modeling of CRF-BP and of full-length CRF2αR and CRF2βR based on the recently solved crystal structures of the transmembrane domains of the human glucagon receptor and human CRF1R, in addition with the resolved N-terminal extracellular domain of CRFRs. These models were further studied using molecular dynamics simulations and protein-protein docking. The results predicted a higher possibility of interaction of CRF-BP with CRF2αR than CRF2βR and yielded the possible residues conforming the interacting interface. Thus, the present study provides a framework for further investigation of the CRF-BP/CRF2αR interaction.
Collapse
Affiliation(s)
- Paula G. Slater
- Department of Cellular and Molecular Biology, Faculty of Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
| | | | - Katia Gysling
- Department of Cellular and Molecular Biology, Faculty of Biological Sciences, Pontificia Universidad Católica de Chile, Santiago, Chile
- *Correspondence: Katia Gysling, ; Carlos F. Lagos,
| | - Carlos F. Lagos
- Department of Endocrinology, School of Medicine, Pontificia Universidad Católica de Chile, Santiago, Chile
- *Correspondence: Katia Gysling, ; Carlos F. Lagos,
| |
Collapse
|
221
|
Abstract
Functional genomics encompasses diverse disciplines in molecular biology and bioinformatics to comprehend the blueprint, regulation, and expression of genetic elements that define the physiology of an organism. The deluge of sequencing data in the postgenomics era has demanded the involvement of computer scientists and mathematicians to create algorithms, analytical software, and databases for the storage, curation, and analysis of biological big data. In this chapter, we discuss on the concept of functional genomics in the context of systems biology and provide examples of its application in human genetic disease studies, molecular crop improvement, and metagenomics for antibiotic discovery. An overview of transcriptomics workflow and experimental considerations is also introduced. Lastly, we present an in-house case study of transcriptomics analysis of an aromatic herbal plant to understand the effect of elicitation on the biosynthesis of volatile organic compounds.
Collapse
Affiliation(s)
- Hoe-Han Goh
- Institute of Systems Biology, Universiti Kebangsaan Malaysia (UKM), Bangi, Malaysia.
| | - Chyan Leong Ng
- Institute of Systems Biology, Universiti Kebangsaan Malaysia (UKM), Bangi, Malaysia
| | - Kok-Keong Loke
- Institute of Systems Biology, Universiti Kebangsaan Malaysia (UKM), Bangi, Malaysia
| |
Collapse
|
222
|
Loo JSE, Emtage AL, Ng KW, Yong ASJ, Doughty SW. Assessing GPCR homology models constructed from templates of various transmembrane sequence identities: Binding mode prediction and docking enrichment. J Mol Graph Model 2017; 80:38-47. [PMID: 29306746 DOI: 10.1016/j.jmgm.2017.12.017] [Citation(s) in RCA: 18] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2017] [Revised: 11/27/2017] [Accepted: 12/26/2017] [Indexed: 11/15/2022]
Abstract
GPCR crystal structures have become more readily accessible in recent years. However, homology models of GPCRs continue to play an important role as many GPCR structures remain unsolved. The new crystal structures now available provide not only additional templates for homology modelling but also the opportunity to assess the performance of homology models against their respective crystal structures and gain insight into the performance of such models. In this study we have constructed homology models from templates of various transmembrane sequence identities for eight GPCR targets to better understand the relationship between transmembrane sequence identity and model quality. Model quality was assessed relative to the crystal structure in terms of structural accuracy as well as performance in two typical structure-based drug design applications: ligand binding pose prediction and docking enrichment in virtual screening. Crystal structures significantly outperformed homology models in both assessments. Accurate ligand binding pose prediction was possible but difficult to achieve using homology models, even with the use of induced fit docking. In virtual screening using homology models still conferred significant enrichment compared to random selection, with a clear benefit also observed in using models optimized through induced fit docking. Our results indicate that while homology models that are reasonably accurate structurally can be constructed, without significant refinement homology models will be outperformed by crystal structures in ligand binding pose prediction and docking enrichment regardless of the template used, primarily due to the extremely high level of structural accuracy needed for such applications.
Collapse
Affiliation(s)
- Jason S E Loo
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia.
| | - Abigail L Emtage
- School of Pharmacy, The University of Nottingham Malaysia Campus, Jalan Broga, 43500 Semenyih, Selangor, Malaysia
| | - Kar Weng Ng
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia
| | - Alene S J Yong
- School of Pharmacy, Taylor's University, No.1 Jalan Taylor's, 47500 Subang Jaya, Selangor, Malaysia
| | - Stephen W Doughty
- Penang Medical College, No. 4 Jalan Sepoy Lines, 10450 George Town, Penang, Malaysia
| |
Collapse
|
223
|
Suplatov D, Sharapova Y, Timonina D, Kopylov K, Švedas V. The visualCMAT: A web-server to select and interpret correlated mutations/co-evolving residues in protein families. J Bioinform Comput Biol 2017; 16:1840005. [PMID: 29361894 DOI: 10.1142/s021972001840005x] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
The visualCMAT web-server was designed to assist experimental research in the fields of protein/enzyme biochemistry, protein engineering, and drug discovery by providing an intuitive and easy-to-use interface to the analysis of correlated mutations/co-evolving residues. Sequence and structural information describing homologous proteins are used to predict correlated substitutions by the Mutual information-based CMAT approach, classify them into spatially close co-evolving pairs, which either form a direct physical contact or interact with the same ligand (e.g. a substrate or a crystallographic water molecule), and long-range correlations, annotate and rank binding sites on the protein surface by the presence of statistically significant co-evolving positions. The results of the visualCMAT are organized for a convenient visual analysis and can be downloaded to a local computer as a content-rich all-in-one PyMol session file with multiple layers of annotation corresponding to bioinformatic, statistical and structural analyses of the predicted co-evolution, or further studied online using the built-in interactive analysis tools. The online interactivity is implemented in HTML5 and therefore neither plugins nor Java are required. The visualCMAT web-server is integrated with the Mustguseal web-server capable of constructing large structure-guided sequence alignments of protein families and superfamilies using all available information about their structures and sequences in public databases. The visualCMAT web-server can be used to understand the relationship between structure and function in proteins, implemented at selecting hotspots and compensatory mutations for rational design and directed evolution experiments to produce novel enzymes with improved properties, and employed at studying the mechanism of selective ligand's binding and allosteric communication between topologically independent sites in protein structures. The web-server is freely available at https://biokinet.belozersky.msu.ru/visualcmat and there are no login requirements.
Collapse
Affiliation(s)
- Dmitry Suplatov
- 1 Belozersky Institute of Physicochemical Biology, Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Leninskiye Gory 1-73, Moscow 119991, Russia
| | - Yana Sharapova
- 1 Belozersky Institute of Physicochemical Biology, Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Leninskiye Gory 1-73, Moscow 119991, Russia
| | - Daria Timonina
- 1 Belozersky Institute of Physicochemical Biology, Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Leninskiye Gory 1-73, Moscow 119991, Russia
| | - Kirill Kopylov
- 1 Belozersky Institute of Physicochemical Biology, Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Leninskiye Gory 1-73, Moscow 119991, Russia
| | - Vytas Švedas
- 1 Belozersky Institute of Physicochemical Biology, Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Leninskiye Gory 1-73, Moscow 119991, Russia
| |
Collapse
|
224
|
Deng L, Fan C, Zeng Z. A sparse autoencoder-based deep neural network for protein solvent accessibility and contact number prediction. BMC Bioinformatics 2017; 18:569. [PMID: 29297299 PMCID: PMC5751690 DOI: 10.1186/s12859-017-1971-7] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/08/2023] Open
Abstract
BACKGROUND Direct prediction of the three-dimensional (3D) structures of proteins from one-dimensional (1D) sequences is a challenging problem. Significant structural characteristics such as solvent accessibility and contact number are essential for deriving restrains in modeling protein folding and protein 3D structure. Thus, accurately predicting these features is a critical step for 3D protein structure building. RESULTS In this study, we present DeepSacon, a computational method that can effectively predict protein solvent accessibility and contact number by using a deep neural network, which is built based on stacked autoencoder and a dropout method. The results demonstrate that our proposed DeepSacon achieves a significant improvement in the prediction quality compared with the state-of-the-art methods. We obtain 0.70 three-state accuracy for solvent accessibility, 0.33 15-state accuracy and 0.74 Pearson Correlation Coefficient (PCC) for the contact number on the 5729 monomeric soluble globular protein dataset. We also evaluate the performance on the CASP11 benchmark dataset, DeepSacon achieves 0.68 three-state accuracy and 0.69 PCC for solvent accessibility and contact number, respectively. CONCLUSIONS We have shown that DeepSacon can reliably predict solvent accessibility and contact number with stacked sparse autoencoder and a dropout approach.
Collapse
Affiliation(s)
- Lei Deng
- School of Software, Central South University, No.22 Shaoshan South Road, Changsha, 410075 China
| | - Chao Fan
- School of Software, Central South University, No.22 Shaoshan South Road, Changsha, 410075 China
| | - Zhiwen Zeng
- School of Information Science and Engineering, Central South University, No.932 South Lushan Road, Changsha, 410083 China
| |
Collapse
|
225
|
Affiliation(s)
- Rahul Kaushik
- Kusuma
School of Biological Sciences, Indian Institute of Technology, Delhi, India
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology, Delhi, India
| | - Ankita Singh
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology, Delhi, India
- Department
of Bioinformatics, Banasthali Vidyapith, Banasthali, India
| | - B. Jayaram
- Kusuma
School of Biological Sciences, Indian Institute of Technology, Delhi, India
- Supercomputing Facility for Bioinformatics & Computational Biology, Indian Institute of Technology, Delhi, India
- Department
of Chemistry, Indian Institute of Technology, Delhi, India
| |
Collapse
|
226
|
Glusman G, Rose PW, Prlić A, Dougherty J, Duarte JM, Hoffman AS, Barton GJ, Bendixen E, Bergquist T, Bock C, Brunk E, Buljan M, Burley SK, Cai B, Carter H, Gao J, Godzik A, Heuer M, Hicks M, Hrabe T, Karchin R, Leman JK, Lane L, Masica DL, Mooney SD, Moult J, Omenn GS, Pearl F, Pejaver V, Reynolds SM, Rokem A, Schwede T, Song S, Tilgner H, Valasatava Y, Zhang Y, Deutsch EW. Mapping genetic variations to three-dimensional protein structures to enhance variant interpretation: a proposed framework. Genome Med 2017; 9:113. [PMID: 29254494 PMCID: PMC5735928 DOI: 10.1186/s13073-017-0509-y] [Citation(s) in RCA: 44] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/10/2022] Open
Abstract
The translation of personal genomics to precision medicine depends on the accurate interpretation of the multitude of genetic variants observed for each individual. However, even when genetic variants are predicted to modify a protein, their functional implications may be unclear. Many diseases are caused by genetic variants affecting important protein features, such as enzyme active sites or interaction interfaces. The scientific community has catalogued millions of genetic variants in genomic databases and thousands of protein structures in the Protein Data Bank. Mapping mutations onto three-dimensional (3D) structures enables atomic-level analyses of protein positions that may be important for the stability or formation of interactions; these may explain the effect of mutations and in some cases even open a path for targeted drug development. To accelerate progress in the integration of these data types, we held a two-day Gene Variation to 3D (GVto3D) workshop to report on the latest advances and to discuss unmet needs. The overarching goal of the workshop was to address the question: what can be done together as a community to advance the integration of genetic variants and 3D protein structures that could not be done by a single investigator or laboratory? Here we describe the workshop outcomes, review the state of the field, and propose the development of a framework with which to promote progress in this arena. The framework will include a set of standard formats, common ontologies, a common application programming interface to enable interoperation of the resources, and a Tool Registry to make it easy to find and apply the tools to specific analysis problems. Interoperability will enable integration of diverse data sources and tools and collaborative development of variant effect prediction methods.
Collapse
Affiliation(s)
| | - Peter W Rose
- San Diego Supercomputer Center, University of California San Diego, La Jolla, CA, 98093, USA
| | - Andreas Prlić
- San Diego Supercomputer Center, University of California San Diego, La Jolla, CA, 98093, USA.,RCSB Protein Data Bank, University of California San Diego, La Jolla, CA, 98093, USA
| | | | - José M Duarte
- RCSB Protein Data Bank, University of California San Diego, La Jolla, CA, 98093, USA
| | - Andrew S Hoffman
- Human Centered Design & Engineering, University of Washington, Seattle, WA, 98195, USA
| | - Geoffrey J Barton
- Division of Computational Biology, School of Life Sciences, University of Dundee, Dundee, DD1 5EH, UK
| | - Emøke Bendixen
- Department of Molecular Biology and Genetics, Aarhus University, 8000, Aarhus, Denmark
| | - Timothy Bergquist
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA
| | - Christian Bock
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA
| | - Elizabeth Brunk
- University of California San Diego, La Jolla, CA, 92093, USA
| | - Marija Buljan
- Institute of Molecular Systems Biology, ETH Zurich, CH-8093, Zurich, Switzerland
| | - Stephen K Burley
- San Diego Supercomputer Center, University of California San Diego, La Jolla, CA, 98093, USA.,RCSB Protein Data Bank, University of California San Diego, La Jolla, CA, 98093, USA.,Institute for Quantitative Biomedicine, Rutgers, The State University of New Jersey, Piscataway, NJ, 08854, USA
| | - Binghuang Cai
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA
| | - Hannah Carter
- University of California San Diego, La Jolla, CA, 92093, USA
| | - JianJiong Gao
- Kravis Center for Molecular Oncology, Memorial Sloan Kettering Cancer Center, New York, NY, 10065, USA
| | - Adam Godzik
- SBP Medical Discovery Institute, La Jolla, CA, 92037, USA
| | - Michael Heuer
- AMPLab, University of California, Berkeley, CA, 94720, USA
| | | | - Thomas Hrabe
- SBP Medical Discovery Institute, La Jolla, CA, 92037, USA
| | - Rachel Karchin
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, 21218, USA.,Department of Oncology, Johns Hopkins Medicine, Baltimore, MD, 21287, USA
| | - Julia Koehler Leman
- Flatiron Institute, Center for Computational Biology, Simons Foundation, New York, NY, 10010, USA.,Department of Biology and Center for Genomics and Systems Biology, New York University, New York, NY, 10003, USA
| | - Lydie Lane
- SIB Swiss Institute of Bioinformatics and University of Geneva, CH-1211, Geneva, Switzerland
| | - David L Masica
- Department of Biomedical Engineering, Institute for Computational Medicine, Johns Hopkins University, Baltimore, MD, 21218, USA
| | - Sean D Mooney
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA
| | - John Moult
- Institute for Bioscience and Biotechnology Research, University of Maryland, Rockville, MD, 20850, USA.,Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, MD, 20742, USA
| | - Gilbert S Omenn
- Institute for Systems Biology, Seattle, WA, 98109, USA.,Department of Computational Medicine & Bioinformatics, University of Michigan, Ann Arbor, MI, 48109-2218, USA
| | - Frances Pearl
- School of Life Sciences, University of Sussex, Brighton, BN1 9QG, UK
| | - Vikas Pejaver
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA.,The University of Washington eScience Institute, Seattle, WA, 98195, USA
| | | | - Ariel Rokem
- The University of Washington eScience Institute, Seattle, WA, 98195, USA
| | - Torsten Schwede
- SIB Swiss Institute of Bioinformatics and Biozentrum University of Basel, CH-4056, Basel, Switzerland
| | - Sicheng Song
- Department of Biomedical Informatics and Medical Education, University of Washington, Seattle, WA, 98109, USA
| | - Hagen Tilgner
- Brain and Mind Research Institute, Weill Cornell Medicine, New York City, NY, 10021, USA
| | - Yana Valasatava
- RCSB Protein Data Bank, University of California San Diego, La Jolla, CA, 98093, USA
| | - Yang Zhang
- Department of Computational Medicine & Bioinformatics, University of Michigan, Ann Arbor, MI, 48109-2218, USA
| | | |
Collapse
|
227
|
Dai T, Yan X, Li Q, Li T, Liu C, McClements DJ, Chen J. Characterization of binding interaction between rice glutelin and gallic acid: Multi-spectroscopic analyses and computational docking simulation. Food Res Int 2017; 102:274-281. [DOI: 10.1016/j.foodres.2017.09.020] [Citation(s) in RCA: 77] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2017] [Revised: 09/02/2017] [Accepted: 09/08/2017] [Indexed: 10/18/2022]
|
228
|
Gludovacz E, Maresch D, Lopes de Carvalho L, Puxbaum V, Baier LJ, Sützl L, Guédez G, Grünwald-Gruber C, Ulm B, Pils S, Ristl R, Altmann F, Jilma B, Salminen TA, Borth N, Boehm T. Oligomannosidic glycans at Asn-110 are essential for secretion of human diamine oxidase. J Biol Chem 2017; 293:1070-1087. [PMID: 29187599 DOI: 10.1074/jbc.m117.814244] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/25/2017] [Revised: 11/14/2017] [Indexed: 01/28/2023] Open
Abstract
N-Glycosylation plays a fundamental role in many biological processes. Human diamine oxidase (hDAO), required for histamine catabolism, has multiple N-glycosylation sites, but their roles, for example in DAO secretion, are unclear. We recently reported that the N-glycosylation sites Asn-168, Asn-538, and Asn-745 in recombinant hDAO (rhDAO) carry complex-type glycans, whereas Asn-110 carries only mammalian-atypical oligomannosidic glycans. Here, we show that Asn-110 in native hDAO from amniotic fluid and Caco-2 cells, DAO from porcine kidneys, and rhDAO produced in two different HEK293 cell lines is also consistently occupied by oligomannosidic glycans. Glycans at Asn-168 were predominantly sialylated with bi- to tetra-antennary branches, and Asn-538 and Asn-745 had similar complex-type glycans with some tissue- and cell line-specific variations. The related copper-containing amine oxidase human vascular adhesion protein-1 also exclusively displayed high-mannose glycosylation at Asn-137. X-ray structures revealed that the residues adjacent to Asn-110 and Asn-137 form a highly conserved hydrophobic cleft interacting with the core trisaccharide. Asn-110 replacement with Gln completely abrogated rhDAO secretion and caused retention in the endoplasmic reticulum. Mutations of Asn-168, Asn-538, and Asn-745 reduced rhDAO secretion by 13, 71, and 32%, respectively. Asn-538/745 double and Asn-168/538/745 triple substitutions reduced rhDAO secretion by 85 and 94%. Because of their locations in the DAO structure, Asn-538 and Asn-745 glycosylations might be important for efficient DAO dimer formation. These functional results are reflected in the high evolutionary conservation of all four glycosylation sites. Human DAO is abundant only in the gastrointestinal tract, kidney, and placenta, and glycosylation seems essential for reaching high enzyme expression levels in these tissues.
Collapse
Affiliation(s)
- Elisabeth Gludovacz
- From the Departments of Biotechnology.,the Departments of Clinical Pharmacology and
| | | | - Leonor Lopes de Carvalho
- the Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, Tykistökatu 6A, 20520 Turku, Finland
| | | | | | - Leander Sützl
- Food Science and Technology, University of Natural Resources and Life Sciences, Muthgasse 18, 1190 Vienna, Austria
| | - Gabriela Guédez
- the Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, Tykistökatu 6A, 20520 Turku, Finland
| | | | | | | | - Robin Ristl
- the Section for Medical Statistics (IMS), Center of Medical Statistics, Informatics and Intelligent Systems, Medical University of Vienna, Waehringer Guertel 18-20, 1090 Vienna, Austria, and
| | | | - Bernd Jilma
- the Departments of Clinical Pharmacology and
| | - Tiina A Salminen
- the Structural Bioinformatics Laboratory, Biochemistry, Faculty of Science and Engineering, Åbo Akademi University, Tykistökatu 6A, 20520 Turku, Finland
| | | | | |
Collapse
|
229
|
Hydrophobicity of proteins and nanostructured solutes is governed by topographical and chemical context. Proc Natl Acad Sci U S A 2017; 114:13345-13350. [PMID: 29158409 DOI: 10.1073/pnas.1700092114] [Citation(s) in RCA: 79] [Impact Index Per Article: 9.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Hydrophobic interactions drive many important biomolecular self-assembly phenomena. However, characterizing hydrophobicity at the nanoscale has remained a challenge due to its nontrivial dependence on the chemistry and topography of biomolecular surfaces. Here we use molecular simulations coupled with enhanced sampling methods to systematically displace water molecules from the hydration shells of nanostructured solutes and calculate the free energetics of interfacial water density fluctuations, which quantify the extent of solute-water adhesion, and therefore solute hydrophobicity. In particular, we characterize the hydrophobicity of curved graphene sheets, self-assembled monolayers (SAMs) with chemical patterns, and mutants of the protein hydrophobin-II. We find that water density fluctuations are enhanced near concave nonpolar surfaces compared with those near flat or convex ones, suggesting that concave surfaces are more hydrophobic. We also find that patterned SAMs and protein mutants, having the same number of nonpolar and polar sites but different geometrical arrangements, can display significantly different strengths of adhesion with water. Specifically, hydroxyl groups reduce the hydrophobicity of methyl-terminated SAMs most effectively not when they are clustered together but when they are separated by one methyl group. Hydrophobin-II mutants show that a charged amino acid reduces the hydrophobicity of a large nonpolar patch when placed at its center, rather than at its edge. Our results highlight the power of water density fluctuations-based measures to characterize the hydrophobicity of nanoscale surfaces and caution against the use of additive approximations, such as the commonly used surface area models or hydropathy scales for characterizing biomolecular hydrophobicity and the associated driving forces of assembly.
Collapse
|
230
|
Martínez-Archundia M, Colín-Astudillo B, Moreno-Vargas LM, Ramírez-Galicia G, Garduño-Juárez R, Deeb O, Contreras-Romo MC, Quintanar-Stephano A, Abarca-Rojano E, Correa-Basurto J. Ligand recognition properties of the vasopressin V2 receptor studied under QSAR and molecular modeling strategies. Chem Biol Drug Des 2017; 90:840-853. [DOI: 10.1111/cbdd.13005] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2016] [Revised: 04/06/2017] [Accepted: 04/08/2017] [Indexed: 11/28/2022]
Affiliation(s)
- Marlet Martínez-Archundia
- Laboratorio de Modelado Molecular y Diseño de Fármacos; Escuela Superior de Medicina-Instituto Politécnico Nacional; México City Mexico
| | - Brenda Colín-Astudillo
- Laboratorio de Modelado Molecular y Diseño de Fármacos; Escuela Superior de Medicina-Instituto Politécnico Nacional; México City Mexico
| | - Liliana M. Moreno-Vargas
- Unidad de Investigación en Enfermedades Oncológicas; Hospital Infantil de México; Mexico City México
| | | | - Ramón Garduño-Juárez
- Instituto de Ciencias Físicas; Universidad Nacional Autónoma de México; Cuernavaca Morelos Mexico
| | - Omar Deeb
- Faculty of Pharmacy; Al-Quds University; Jerusalem Palestine
| | - Martha Citlalli Contreras-Romo
- Departamento de Fisiología y Farmacología; Centro de Ciencias Básicas; Universidad Autónoma de Aguascalientes; Aguascalientes Mexico
| | - Andres Quintanar-Stephano
- Departamento de Fisiología y Farmacología; Centro de Ciencias Básicas; Universidad Autónoma de Aguascalientes; Aguascalientes Mexico
| | - Edgar Abarca-Rojano
- Sección de Estudios de Posgrado e Investigación, Escuela Superior de Medicina; Instituto Politécnico Nacional; Mexico DF Mexico
| | - José Correa-Basurto
- Laboratorio de Modelado Molecular y Diseño de Fármacos; Escuela Superior de Medicina-Instituto Politécnico Nacional; México City Mexico
| |
Collapse
|
231
|
Xie H, Zeng D, Chen X, Huo D, Liu L, Zhang D, Jin Q, Ke K, Hu M. Prediction on the risk population of idiosyncratic adverse reactions based on molecular docking with mutant proteins. Oncotarget 2017; 8:95568-95576. [PMID: 29221149 PMCID: PMC5707043 DOI: 10.18632/oncotarget.21509] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2017] [Accepted: 09/20/2017] [Indexed: 01/11/2023] Open
Abstract
Idiosyncratic adverse drug reactions are drug reactions that occur rarely and unpredictably among the population. These reactions often occur after a drug is marketed, which means that they are strongly related to the genotype of the population. The prediction of such adverse reactions is a major challenge because of the lack of appropriate test models during the drug development process. In this study, we chose withdrawn drugs because the reasons why they were withdrawn and from which countries or regions is easily obtained. We selected Dilevalol and its chiral drug (Labetalol) as the investigatory drugs, as they have been withdrawn from a European market (Britain) because of serious hepatotoxicity. First, we searched for and obtained the Dilevalol-induced- liver-injury related protein, multidrug resistance protein 1 (MDR1), from the Comparative Toxicogenomics Database (CTD). Then, we searched and extracted 477 non-synonymous single nucleotide polymorphisms (nsSNP) on MDR1 in the dbSNP database. Second, we used the VarMod tool to predict the functional changes of MDR1 induced by these nsSNPs, from which we extracted the nsSNPs that significantly change the functions of this protein. Third, we built the three-dimensional structures of those variant proteins and used AutoDock to perform a docking study, choosing the best model to determine the sites of nsSNPs. Finally, we used the data from the 1000 Genomes Project to verify the dominant population distribution of the risk SNP. We applied the same strategy to the post-marketing drug-induced liver injury drugs to further test the feasibility of our method.
Collapse
Affiliation(s)
- Hongbo Xie
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Diheng Zeng
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Xiujie Chen
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Diwei Huo
- The 2nd Affiliated Hospital of Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Lei Liu
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Denan Zhang
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Qing Jin
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Kehui Ke
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| | - Ming Hu
- Department of Pharmacogenomics, College of Bioinformatics Science and Technology, Harbin Medical University, Harbin, Heilongjiang 150081, PR China
| |
Collapse
|
232
|
Kihara D, Yang YD, Hawkins T. Bioinformatics Resources for Cancer Research with an Emphasis on Gene Function and Structure Prediction Tools. Cancer Inform 2017. [DOI: 10.1177/117693510600200020] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/19/2023] Open
Abstract
The immensely popular fields of cancer research and bioinformatics overlap in many different areas, e.g. large data repositories that allow for users to analyze data from many experiments (data handling, databases), pattern mining, microarray data analysis, and interpretation of proteomics data. There are many newly available resources in these areas that may be unfamiliar to most cancer researchers wanting to incorporate bioinformatics tools and analyses into their work, and also to bioinformaticians looking for real data to develop and test algorithms. This review reveals the interdependence of cancer research and bioinformatics, and highlight the most appropriate and useful resources available to cancer researchers. These include not only public databases, but general and specific bioinformatics tools which can be useful to the cancer researcher. The primary foci are function and structure prediction tools of protein genes. The result is a useful reference to cancer researchers and bioinformaticians studying cancer alike.
Collapse
Affiliation(s)
- Daisuke Kihara
- Department of Biological Sciences; College of Science, Purdue University, West Lafayette, IN, 47907, USA
- Department of Computer Science; College of Science, Purdue University, West Lafayette, IN, 47907, USA
- Markey Center for Structural Biology; College of Science, Purdue University, West Lafayette, IN, 47907, USA
- The Bindley Bioscience Center, College of Science, Purdue University, West Lafayette, IN, 47907, USA
| | - Yifeng David Yang
- Department of Biological Sciences; College of Science, Purdue University, West Lafayette, IN, 47907, USA
| | - Troy Hawkins
- Department of Biological Sciences; College of Science, Purdue University, West Lafayette, IN, 47907, USA
| |
Collapse
|
233
|
Poondla V, Chikati R, Kallubai M, Chennupati V, Subramanyam R, Obulam VSR. Characterization and molecular modeling of polygalacturonase isoforms from Saccharomyces cerevisiae. 3 Biotech 2017; 7:285. [PMID: 28828292 DOI: 10.1007/s13205-017-0912-5] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2017] [Accepted: 08/02/2017] [Indexed: 10/19/2022] Open
Abstract
Earlier, low-temperature-active polygalacturonase isoforms from Saccharomyces cerevisiae PVK4 were isolated and purified. Substrate specificity of polygalacturonase isoforms indicated high affinity for pectins and very low enzyme activity towards non-pectic polysaccharides. To characterize the polygalacturonase isoforms, biochemical, spectral, and in silico approaches were used. The apparent Km and Vmax values for hydrolysis of pectin and galacturonic acid were 0.31 mg/ml and 3.15 mmol min/mg, respectively. Interestingly, the polygalacturonase isoforms were found to be more stable at optimal pH and temperature of 4.5 and 40 °C, respectively. These isoforms were reacted with different metal ions; Cd2+ and Ni2+ severely inhibited the enzyme activity, while Mg2+, Zn2+, Cd2+, Fe2+ Cu2+, and Ni2+ inhibited to a lesser extent, which clearly demonstrated that variations in enzyme activity were due to their differential binding affinity of metal ions. Furthermore, decrease in the viscosity of polygalacturonic acid and citrus pectin by these isoforms was approximately four and six times higher than the rate of release of reducing sugars. This indicates that polygalacturonase isoforms have an endo-mode of action. In addition to the above, thermostability of purified polygalacturonase isoforms was studied by circular dichroism and fluorescence spectroscopy. Circular dichroism showed 18% alpha helix and 57% beta sheets at pH 5, while at pH 7, 8, and 9 there was an increase of random coil. Fluorescence studies revealed small conformational changes, which were observed at 30-50 °C, while unfolding transition region was noticed between 60 and 70 °C. The purified enzyme fractions were analyzed by MALDI-TOF MS. Finally, 3D model structures for isoenzymes of polygalacturonase of S. cerevisiae were generated and validated as good quality models, which are also suitable for molecular interaction studies.
Collapse
|
234
|
Zeng J, Zhu Q, Wu Y, Chen H, Lin X. Characterization of a polycyclic aromatic ring-hydroxylation dioxygenase from Mycobacterium sp. NJS-P. CHEMOSPHERE 2017; 185:67-74. [PMID: 28686888 DOI: 10.1016/j.chemosphere.2017.07.001] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2016] [Revised: 06/01/2017] [Accepted: 07/01/2017] [Indexed: 06/07/2023]
Abstract
Ring-hydroxylating dioxygenases (RHDs) play a critical role in the biodegradation of polycyclic aromatic hydrocarbons (PAHs). In this study, genes pdoAB encoding a dioxygenase capable of oxidizing various PAHs with up to five-ring benzo[a]pyrene were cloned from Mycobacterium sp. NJS-P. The α-subunit of the PdoAB showed 99% and 93% identity to that from Mycobacterium sp. S65 and Mycobacterium sp. py136, respectively. An Escherichia coli expression experiment revealed that the enzyme is able to oxidize anthracene, phenanthrene, pyrene and benzo[a]pyrene, but not to fluoranthene and benzo[a]anthracene. Furthermore, the results of in silico analysis showed that PdoAB has a large substrate-binding pocket satisfying for accommodation of HMW PAHs, and suggested that the binding energy of intermolecular interaction may predict the substrate conversion of RHDs towards HMW PAHs, especially those may have steric constraints on the substrate-binding pocket, such as benzo[a]pyrene and benzo[a]anthracene.
Collapse
Affiliation(s)
- Jun Zeng
- Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Beijing East Road, 71, Nanjing 210008, PR China; Joint Open Laboratory of Soil and the Environment, Hong Kong University and Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, PR China
| | - Qinghe Zhu
- Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Beijing East Road, 71, Nanjing 210008, PR China; Joint Open Laboratory of Soil and the Environment, Hong Kong University and Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, PR China
| | - Yucheng Wu
- Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Beijing East Road, 71, Nanjing 210008, PR China; Joint Open Laboratory of Soil and the Environment, Hong Kong University and Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, PR China
| | - Hong Chen
- Soil and Environment Analysis Center, Institute of Soil Science, Chinese Academy of Science, PR China
| | - Xiangui Lin
- Key Laboratory of Soil Environment and Pollution Remediation, Institute of Soil Science, Chinese Academy of Sciences, Beijing East Road, 71, Nanjing 210008, PR China; Joint Open Laboratory of Soil and the Environment, Hong Kong University and Institute of Soil Science, Chinese Academy of Sciences, Nanjing 210008, PR China.
| |
Collapse
|
235
|
Monzon AM, Zea DJ, Marino-Buslje C, Parisi G. Homology modeling in a dynamical world. Protein Sci 2017; 26:2195-2206. [PMID: 28815769 DOI: 10.1002/pro.3274] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2017] [Revised: 08/09/2017] [Accepted: 08/09/2017] [Indexed: 12/31/2022]
Abstract
A key concept in template-based modeling (TBM) is the high correlation between sequence and structural divergence, with the practical consequence that homologous proteins that are similar at the sequence level will also be similar at the structural level. However, conformational diversity of the native state will reduce the correlation between structural and sequence divergence, because structural variation can appear without sequence diversity. In this work, we explore the impact that conformational diversity has on the relationship between structural and sequence divergence. We find that the extent of conformational diversity can be as high as the maximum structural divergence among families. Also, as expected, conformational diversity impairs the well-established correlation between sequence and structural divergence, which is nosier than previously suggested. However, we found that this noise can be resolved using a priori information coming from the structure-function relationship. We show that protein families with low conformational diversity show a well-correlated relationship between sequence and structural divergence, which is severely reduced in proteins with larger conformational diversity. This lack of correlation could impair TBM results in highly dynamical proteins. Finally, we also find that the presence of order/disorder can provide useful beforehand information for better TBM performance.
Collapse
Affiliation(s)
- Alexander Miguel Monzon
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, CONICET, B1876BXD, Bernal, Argentina
| | - Diego Javier Zea
- Structural Bioinformatics Unit, Fundación Instituto Leloir, CONICET, C1405BWE Ciudad Autónoma de Buenos Aires, Argentina
| | - Cristina Marino-Buslje
- Structural Bioinformatics Unit, Fundación Instituto Leloir, CONICET, C1405BWE Ciudad Autónoma de Buenos Aires, Argentina
| | - Gustavo Parisi
- Departamento de Ciencia y Tecnología, Universidad Nacional de Quilmes, CONICET, B1876BXD, Bernal, Argentina
| |
Collapse
|
236
|
Fierro F, Suku E, Alfonso-Prieto M, Giorgetti A, Cichon S, Carloni P. Agonist Binding to Chemosensory Receptors: A Systematic Bioinformatics Analysis. Front Mol Biosci 2017; 4:63. [PMID: 28932739 PMCID: PMC5592726 DOI: 10.3389/fmolb.2017.00063] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2017] [Accepted: 08/22/2017] [Indexed: 12/17/2022] Open
Abstract
Human G-protein coupled receptors (hGPCRs) constitute a large and highly pharmaceutically relevant membrane receptor superfamily. About half of the hGPCRs' family members are chemosensory receptors, involved in bitter taste and olfaction, along with a variety of other physiological processes. Hence these receptors constitute promising targets for pharmaceutical intervention. Molecular modeling has been so far the most important tool to get insights on agonist binding and receptor activation. Here we investigate both aspects by bioinformatics-based predictions across all bitter taste and odorant receptors for which site-directed mutagenesis data are available. First, we observe that state-of-the-art homology modeling combined with previously used docking procedures turned out to reproduce only a limited fraction of ligand/receptor interactions inferred by experiments. This is most probably caused by the low sequence identity with available structural templates, which limits the accuracy of the protein model and in particular of the side-chains' orientations. Methods which transcend the limited sampling of the conformational space of docking may improve the predictions. As an example corroborating this, we review here multi-scale simulations from our lab and show that, for the three complexes studied so far, they significantly enhance the predictive power of the computational approach. Second, our bioinformatics analysis provides support to previous claims that several residues, including those at positions 1.50, 2.50, and 7.52, are involved in receptor activation.
Collapse
Affiliation(s)
- Fabrizio Fierro
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany
| | - Eda Suku
- Department of Biotechnology, University of VeronaVerona, Italy
| | - Mercedes Alfonso-Prieto
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Cécile and Oskar Vogt Institute for Brain Research, Medical Faculty, Heinrich Heine University DüsseldorfDüsseldorf, Germany
| | - Alejandro Giorgetti
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Department of Biotechnology, University of VeronaVerona, Italy
| | - Sven Cichon
- Institute of Neuroscience and Medicine INM-1, Forschungszentrum JülichJülich, Germany.,Institute for Human Genetics, Department of Genomics, Life&Brain Center, University of BonnBonn, Germany.,Division of Medical Genetics, Department of Biomedicine, University of BaselBasel, Switzerland
| | - Paolo Carloni
- Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum JülichJülich, Germany.,Department of Physics, Rheinisch-Westfälische Technische Hochschule AachenAachen, Germany.,VNU Key Laboratory "Multiscale Simulation of Complex Systems", VNU University of Science, Vietnam National UniversityHanoi, Vietnam
| |
Collapse
|
237
|
Terashi G, Kihara D. Protein structure model refinement in CASP12 using short and long molecular dynamics simulations in implicit solvent. Proteins 2017; 86 Suppl 1:189-201. [PMID: 28833585 DOI: 10.1002/prot.25373] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/21/2017] [Revised: 08/01/2017] [Accepted: 08/18/2017] [Indexed: 12/21/2022]
Abstract
Protein structure prediction has matured over years, particularly those which use structure templates for building a model. It can build a model with correct overall conformation in cases where appropriate templates are available. Models with the correct topology can be practically useful for limited purposes that need residue-level accuracy, but further improvement of the models can allow the models to be used in tasks that need detailed structures, such as molecular replacement in X-ray crystallography or structure-based drug screening. Thus, model refinement is an important final step in protein structure prediction to bridge predictions to real-life applications. Model refinement is one of the categories in recent rounds of critical assessment of techniques in protein structure prediction (CASP) and has recently been drawing more attention due to its realized importance. Here we report our group's performance in the refinement category in CASP12. Our method is based on inexpensive short molecular dynamics (MD) simulations in implicit solvent. Our performance in CASP12 was among the top, which was consistent with the previous round, CASP11. Our method with short MD runs achieved comparable performance with other methods that used longer simulations. Detailed analyses found that improvements typically occurred in entire regions of a structure rather than only in flexible loop regions. The remaining challenge in the structure refinement includes large conformational refinement which involves substantial motions of secondary structure elements or domains.
Collapse
Affiliation(s)
- Genki Terashi
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907
| | - Daisuke Kihara
- Department of Biological Sciences, Purdue University, West Lafayette, Indiana, 47907.,Department of Computer Science, Purdue University, West Lafayette, Indiana, 47907
| |
Collapse
|
238
|
Sirovetz BJ, Schafer NP, Wolynes PG. Protein structure prediction: making AWSEM AWSEM-ER by adding evolutionary restraints. Proteins 2017; 85:2127-2142. [PMID: 28799172 DOI: 10.1002/prot.25367] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/29/2017] [Revised: 07/29/2017] [Accepted: 08/08/2017] [Indexed: 11/07/2022]
Abstract
Protein sequences have evolved to fold into functional structures, resulting in families of diverse protein sequences that all share the same overall fold. One can harness protein family sequence data to infer likely contacts between pairs of residues. In the current study, we combine this kind of inference from coevolutionary information with a coarse-grained protein force field ordinarily used with single sequence input, the Associative memory, Water mediated, Structure and Energy Model (AWSEM), to achieve improved structure prediction. The resulting Associative memory, Water mediated, Structure and Energy Model with Evolutionary Restraints (AWSEM-ER) yields a significant improvement in the quality of protein structure prediction over the single sequence prediction from AWSEM when a sufficiently large number of homologous sequences are available. Free energy landscape analysis shows that the addition of the evolutionary term shifts the free energy minimum to more native-like structures, which explains the improvement in the quality of structures when performing predictions using simulated annealing. Simulations using AWSEM without coevolutionary information have proved useful in elucidating not only protein folding behavior, but also mechanisms of protein function. The success of AWSEM-ER in de novo structure prediction suggests that the enhanced model opens the door to functional studies of proteins even when no experimentally solved structures are available.
Collapse
Affiliation(s)
- Brian J Sirovetz
- Center for Theoretical Biological Physics, Rice University, Houston, Texas.,Department of Chemistry, Rice University, Houston, Texas
| | - Nicholas P Schafer
- Center for Theoretical Biological Physics, Rice University, Houston, Texas
| | - Peter G Wolynes
- Center for Theoretical Biological Physics, Rice University, Houston, Texas.,Department of Chemistry, Rice University, Houston, Texas.,Department of Physics, Rice University, Houston, Texas.,Department of Biosciences, Rice University, Houston, Texas
| |
Collapse
|
239
|
Xiong D, Mao W, Gong H. Predicting the helix-helix interactions from correlated residue mutations. Proteins 2017; 85:2162-2169. [PMID: 28833538 DOI: 10.1002/prot.25370] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/15/2017] [Revised: 08/03/2017] [Accepted: 08/13/2017] [Indexed: 12/30/2022]
Abstract
Helix-helix interactions are crucial in the structure assembly, stability and function of helix-rich proteins including many membrane proteins. In spite of remarkable progresses over the past decades, the accuracy of predicting protein structures from their amino acid sequences is still far from satisfaction. In this work, we focused on a simpler problem, the prediction of helix-helix interactions, the results of which could facilitate practical protein structure prediction by constraining the sampling space. Specifically, we started from the noisy 2D residue contact maps derived from correlated residue mutations, and utilized ridge detection to identify the characteristic residue contact patterns for helix-helix interactions. The ridge information as well as a few additional features were then fed into a machine learning model HHConPred to predict interactions between helix pairs. In an independent test, our method achieved an F-measure of ∼60% for predicting helix-helix interactions. Moreover, although the model was trained mainly using soluble proteins, it could be extended to membrane proteins with at least comparable performance relatively to previous approaches that were generated purely using membrane proteins. All data and source codes are available at http://166.111.152.91/Downloads.html or https://github.com/dpxiong/HHConPred.
Collapse
Affiliation(s)
- Dapeng Xiong
- MOE Key Laboratory of Bioinformatics, School of Life Sciences, Tsinghua University, Beijing, China
- Beijing Innovation Center of Structural Biology, Tsinghua University, Beijing, China
| | - Wenzhi Mao
- MOE Key Laboratory of Bioinformatics, School of Life Sciences, Tsinghua University, Beijing, China
- Beijing Innovation Center of Structural Biology, Tsinghua University, Beijing, China
| | - Haipeng Gong
- MOE Key Laboratory of Bioinformatics, School of Life Sciences, Tsinghua University, Beijing, China
- Beijing Innovation Center of Structural Biology, Tsinghua University, Beijing, China
| |
Collapse
|
240
|
Colas C, Schlessinger A, Pajor AM. Mapping Functionally Important Residues in the Na +/Dicarboxylate Cotransporter, NaDC1. Biochemistry 2017; 56:4432-4441. [PMID: 28731330 DOI: 10.1021/acs.biochem.7b00503] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Transporters from the SLC13 family couple the transport of two to four Na+ ions with a di- or tricarboxylate, such as succinate or citrate. We have previously modeled mammalian members of the SLC13 family, including the Na+/dicarboxylate cotransporter NaDC1 (SLC13A2), based on a structure of the bacterial homologue VcINDY in an inward-facing conformation with one sodium ion bound at the Na1 site. In the study presented here, we modeled the outward-facing conformation of rabbit and human NaDC1 (rbNaDC1 and hNaDC1, respectively) using an outward-facing model of VcINDY as a template and identified residues in or near the putative Na2 and Na3 cation binding sites. Guided by the structural models in both conformations, we performed site-directed mutagenesis in rbNaDC1 for residues proposed to be in the Na+ or substrate binding sites. Cysteine substitution of T474 in the predicted Na2 binding site results in an inactive protein. The M539C mutant has a low apparent affinity for both sodium and lithium cations, suggesting that M539 may form part of the putative Na3 binding site. The Y432C and T86C mutants have increased Km values for succinate, supporting their proposed location in the outward-facing substrate binding site. In addition, cysteine labeling by MTSEA-biotin shows that Y432C is accessible from the outside of the cell, and the accessibility changes in the presence or absence of Na+. The results of this study improve our understanding of substrate and ion recognition in the mammalian members of the SLC13 family and provide a framework for developing conformationally specific inhibitors against these transporters.
Collapse
Affiliation(s)
- Claire Colas
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai , New York, New York 10029, United States
| | - Avner Schlessinger
- Department of Pharmacological Sciences, Icahn School of Medicine at Mount Sinai , New York, New York 10029, United States
| | - Ana M Pajor
- Skaggs School of Pharmacy and Pharmaceutical Sciences, University of California-San Diego , La Jolla, California 92130-0714, United States
| |
Collapse
|
241
|
Lam SD, Das S, Sillitoe I, Orengo C. An overview of comparative modelling and resources dedicated to large-scale modelling of genome sequences. Acta Crystallogr D Struct Biol 2017; 73:628-640. [PMID: 28777078 PMCID: PMC5571743 DOI: 10.1107/s2059798317008920] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/17/2016] [Accepted: 06/14/2017] [Indexed: 12/02/2022] Open
Abstract
Computational modelling of proteins has been a major catalyst in structural biology. Bioinformatics groups have exploited the repositories of known structures to predict high-quality structural models with high efficiency at low cost. This article provides an overview of comparative modelling, reviews recent developments and describes resources dedicated to large-scale comparative modelling of genome sequences. The value of subclustering protein domain superfamilies to guide the template-selection process is investigated. Some recent cases in which structural modelling has aided experimental work to determine very large macromolecular complexes are also cited.
Collapse
Affiliation(s)
- Su Datt Lam
- Institute of Structural and Molecular Biology, UCL, Darwin Building, Gower Street, London WC1E 6BT, England
- School of Biosciences and Biotechnology, Faculty of Science and Technology, University Kebangsaan Malaysia, 43600 Bangi, Selangor, Malaysia
| | - Sayoni Das
- Institute of Structural and Molecular Biology, UCL, Darwin Building, Gower Street, London WC1E 6BT, England
| | - Ian Sillitoe
- Institute of Structural and Molecular Biology, UCL, Darwin Building, Gower Street, London WC1E 6BT, England
| | - Christine Orengo
- Institute of Structural and Molecular Biology, UCL, Darwin Building, Gower Street, London WC1E 6BT, England
| |
Collapse
|
242
|
Kim HS, Asmis R. Mitogen-activated protein kinase phosphatase 1 (MKP-1) in macrophage biology and cardiovascular disease. A redox-regulated master controller of monocyte function and macrophage phenotype. Free Radic Biol Med 2017; 109:75-83. [PMID: 28330703 PMCID: PMC5462841 DOI: 10.1016/j.freeradbiomed.2017.03.020] [Citation(s) in RCA: 37] [Impact Index Per Article: 4.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 01/13/2017] [Revised: 03/03/2017] [Accepted: 03/17/2017] [Indexed: 12/21/2022]
Abstract
MAPK pathways play a critical role in the activation of monocytes and macrophages by pathogens, signaling molecules and environmental cues and in the regulation of macrophage function and plasticity. MAPK phosphatase 1 (MKP-1) has emerged as the main counter-regulator of MAPK signaling in monocytes and macrophages. Loss of MKP-1 in monocytes and macrophages in response to metabolic stress leads to dysregulation of monocyte adhesion and migration, and gives rise to dysfunctional, proatherogenic monocyte-derived macrophages. Here we review the properties of this redox-regulated dual-specificity MAPK phosphatase and the role of MKP-1 in monocyte and macrophage biology and cardiovascular diseases.
Collapse
Affiliation(s)
- Hong Seok Kim
- Department of Molecular Medicine, College of Medicine, Inha University, Incheon 22212, Republic of Korea; Hypoxia-related Disease Research Center, College of Medicine, Inha University, Incheon 22212, Republic of Korea
| | - Reto Asmis
- Department of Clinical Laboratory Sciences, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA; Department of Biochemistry, University of Texas Health Science Center at San Antonio, San Antonio, TX 78229, USA.
| |
Collapse
|
243
|
Gaded V, Anand R. Selective Deamination of Mutagens by a Mycobacterial Enzyme. J Am Chem Soc 2017; 139:10762-10768. [DOI: 10.1021/jacs.7b04967] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/07/2023]
Affiliation(s)
- Vandana Gaded
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| | - Ruchi Anand
- Department of Chemistry, Indian Institute of Technology Bombay, Powai, Mumbai 400076, India
| |
Collapse
|
244
|
Raghuraman P, Jesu Jaya Sudan R, Lesitha Jeeva Kumari J, Sudandiradoss C. Systematic prioritization of functional hotspot in RIG-1 domains using pattern based conventional molecular dynamic simulation. Life Sci 2017; 184:58-70. [PMID: 28705469 DOI: 10.1016/j.lfs.2017.07.011] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/07/2017] [Revised: 06/29/2017] [Accepted: 07/09/2017] [Indexed: 10/19/2022]
Abstract
BACKGROUND Retinoic acid inducible gene 1 (RIG-1), multi-domain protein has a role-play in detecting viral nucleic acids and stimulates the antiviral response. Dysfunction of this protein due to mutations makes the route vulnerable to viral diseases. AIM Identification of functional hotspots that maintains conformational stability in RIG-1 domains. METHODS In this study, we employed a systematic in silico strategy on RIG-1 protein to understand the mechanism of structural changes upon mutation. We computationally investigated the protein sequence signature for all the three domains of RIG-1 protein that encloses the mutation within the motif. Further, we carried out a structural comparison between RIG-1 domains with their respective distant orthologs which revealed the minimal number of interactions required to maintain its structural fold. This intra-protein network paved the way to infer hotspot residues crucial for the maintenance of the structural architecture and folding pattern. KEY FINDINGS Our analysis revealed about 40 hotspot residues that determine the folding pattern of the RIG-1 domains. Also, conventional molecular dynamic simulation coupled with essential dynamics provides conformational transitions of hot spot residues among native and mutant structures. Structural variations owing to hotspot residues in mutants again confirm the significance of these residues in structural characterization of RIG-1 domains. We believe our results will help the researchers to better comprehend towards regulatory regions and target-binding sites for therapeutic design within the pattern recognition receptor proteins. SIGNIFICANCE Our protocol employed in this work describes a novel approach in identifying signature residues that would provide structural insights in protein folding.
Collapse
Affiliation(s)
- P Raghuraman
- Department of Biotechnology, School of Biosciences and Technology, VIT University, Vellore 632014, India
| | - R Jesu Jaya Sudan
- Department of Biotechnology, School of Biosciences and Technology, VIT University, Vellore 632014, India
| | - J Lesitha Jeeva Kumari
- Department of Biotechnology, School of Biosciences and Technology, VIT University, Vellore 632014, India
| | - C Sudandiradoss
- Department of Biotechnology, School of Biosciences and Technology, VIT University, Vellore 632014, India.
| |
Collapse
|
245
|
Sharifi M. Computational approaches to understand the adverse drug effect on potassium, sodium and calcium channels for predicting TdP cardiac arrhythmias. J Mol Graph Model 2017; 76:152-160. [PMID: 28756335 DOI: 10.1016/j.jmgm.2017.06.012] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2017] [Revised: 06/08/2017] [Accepted: 06/10/2017] [Indexed: 02/08/2023]
Abstract
Ion channels play a crucial role in the cardiovascular system. Our understanding of cardiac ion channel function has improved since their first discoveries. The flow of potassium, sodium and calcium ions across cardiomyocytes is vital for regular cardiac rhythm. Blockage of these channels, delays cardiac repolarization or tend to shorten repolarization and may induce arrhythmia. Detection of drug risk by channel blockade is considered essential for drug regulators. Advanced computational models can be used as an early screen for torsadogenic potential in drug candidates. New drug candidates that are determined to not cause blockage are more likely to pass successfully through preclinical trials and not be withdrawn later from the marketplace by manufacturer. Several different approved drugs, however, can cause a distinctive polymorphic ventricular arrhythmia known as torsade de pointes (TdP), which may lead to sudden death. The objective of the present study is to review the mechanisms and computational models used to assess the risk that a drug may TdP. KEY POINTS There is strong evidence from multiple studies that blockage of the L-type calcium current reduces risk of TdP. Blockage of sodium channels slows cardiac action potential conduction, however, not all sodium channel blocking antiarrhythmic drugs produce a significant effect, while late sodium channel block reduces TdP. Interestingly, there are some drugs that block the hERG potassium channel and therefore cause QT prolongation, but they are not associated with TdP. Recent studies confirmed the necessity of studying multiple distinctionic ion channels which are responsible for cardiac related diseases or TdP, to obtain an improved clinical TdP risk prediction of compound interactions and also for designing drugs.
Collapse
Affiliation(s)
- Mohsen Sharifi
- Division of Systems Biology, National Center for Toxicological Research, US Food and Drug Administration, Jefferson, AR 72079, USA.
| |
Collapse
|
246
|
Brodie NI, Popov KI, Petrotchenko EV, Dokholyan NV, Borchers CH. Solving protein structures using short-distance cross-linking constraints as a guide for discrete molecular dynamics simulations. SCIENCE ADVANCES 2017; 3:e1700479. [PMID: 28695211 PMCID: PMC5501500 DOI: 10.1126/sciadv.1700479] [Citation(s) in RCA: 72] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/13/2017] [Accepted: 05/19/2017] [Indexed: 05/21/2023]
Abstract
We present an integrated experimental and computational approach for de novo protein structure determination in which short-distance cross-linking data are incorporated into rapid discrete molecular dynamics (DMD) simulations as constraints, reducing the conformational space and achieving the correct protein folding on practical time scales. We tested our approach on myoglobin and FK506 binding protein-models for α helix-rich and β sheet-rich proteins, respectively-and found that the lowest-energy structures obtained were in agreement with the crystal structure, hydrogen-deuterium exchange, surface modification, and long-distance cross-linking validation data. Our approach is readily applicable to other proteins with unknown structures.
Collapse
Affiliation(s)
- Nicholas I. Brodie
- University of Victoria–Genome British Columbia Proteomics Centre, Vancouver Island Technology Park, #3101-4464 Markham Street, Victoria, British Columbia V8Z7X8, Canada
| | - Konstantin I. Popov
- Department of Biochemistry and Biophysics, University of North Carolina, Genetic Medicine Building, 120 Mason Farm Road, Chapel Hill, NC 27599, USA
| | - Evgeniy V. Petrotchenko
- University of Victoria–Genome British Columbia Proteomics Centre, Vancouver Island Technology Park, #3101-4464 Markham Street, Victoria, British Columbia V8Z7X8, Canada
| | - Nikolay V. Dokholyan
- Department of Biochemistry and Biophysics, University of North Carolina, Genetic Medicine Building, 120 Mason Farm Road, Chapel Hill, NC 27599, USA
| | - Christoph H. Borchers
- University of Victoria–Genome British Columbia Proteomics Centre, Vancouver Island Technology Park, #3101-4464 Markham Street, Victoria, British Columbia V8Z7X8, Canada
- Department of Biochemistry and Microbiology, University of Victoria, Room 270d, Petch Building, 3800 Finnerty Road, Victoria, British Columbia V8P 5C2, Canada
- Gerald Bronfman Department of Oncology, Jewish General Hospital, Suite 720, 5100 de Maisonneuve Boulevard West, Montreal, Quebec H4A 3T2, Canada
- Proteomics Centre, Segal Cancer Centre, Lady Davis Institute, Jewish General Hospital, McGill University, 3755 Côte-Sainte-Catherine Road, Montreal, Quebec H3T 1E2, Canada
| |
Collapse
|
247
|
|
248
|
Shams H, Soheilypour M, Peyro M, Moussavi-Baygi R, Mofrad MRK. Looking "Under the Hood" of Cellular Mechanotransduction with Computational Tools: A Systems Biomechanics Approach across Multiple Scales. ACS Biomater Sci Eng 2017; 3:2712-2726. [PMID: 33418698 DOI: 10.1021/acsbiomaterials.7b00117] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Signal modulation has been developed in living cells throughout evolution to promote utilizing the same machinery for multiple cellular functions. Chemical and mechanical modules of signal transmission and transduction are interconnected and necessary for organ development and growth. However, due to the high complexity of the intercommunication of physical intracellular connections with biochemical pathways, there are many missing details in our overall understanding of mechanotransduction processes, i.e., the process by which mechanical signals are converted to biochemical cascades. Cell-matrix adhesions are mechanically coupled to the nucleus through the cytoskeleton. This modulated and tightly integrated network mediates the transmission of mechanochemical signals from the extracellular matrix to the nucleus. Various experimental and computational techniques have been utilized to understand the basic mechanisms of mechanotransduction, yet many aspects have remained elusive. Recently, in silico experiments have made important contributions to the field of mechanobiology. Herein, computational modeling efforts devoted to understanding integrin-mediated mechanotransduction pathways are reviewed, and an outlook is presented for future directions toward using suitable computational approaches and developing novel techniques for addressing important questions in the field of mechanotransduction.
Collapse
Affiliation(s)
- Hengameh Shams
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and Mechanical Engineering, University of California, Berkeley, California 94720-1762, United States
| | - Mohammad Soheilypour
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and Mechanical Engineering, University of California, Berkeley, California 94720-1762, United States
| | - Mohaddeseh Peyro
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and Mechanical Engineering, University of California, Berkeley, California 94720-1762, United States
| | - Ruhollah Moussavi-Baygi
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and Mechanical Engineering, University of California, Berkeley, California 94720-1762, United States
| | - Mohammad R K Mofrad
- Molecular Cell Biomechanics Laboratory, Departments of Bioengineering and Mechanical Engineering, University of California, Berkeley, California 94720-1762, United States
| |
Collapse
|
249
|
HÖft TA, Alpert BK. FAST UPDATING MULTIPOLE COULOMBIC POTENTIAL CALCULATION. SIAM JOURNAL ON SCIENTIFIC COMPUTING : A PUBLICATION OF THE SOCIETY FOR INDUSTRIAL AND APPLIED MATHEMATICS 2017; 39:https://doi.org/10.1137/16M1096189. [PMID: 33088167 PMCID: PMC7574401] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Subscribe] [Scholar Register] [Indexed: 06/11/2023]
Abstract
We present a numerical method to efficiently and accurately re-compute the Coulomb potential of a large ensemble of charged particles after a subset of the particles undergoes a change of position. Errors are bounded even after a large number of such shifts, making it practical for use in Monte Carlo Markov chain methods in molecular dynamics, computational astrophysics, computational chemistry, and other applications. The method uses truncated multipole expansions of the potential energy functional and a tree decomposition of the computational domain to reduce the computational complexity. Computational costs scale logarithmically in the size of the problem. Scaling, accuracy, and efficiency are confirmed with numerical experiments. The new method outperforms a direct calculation for moderate problem sizes.
Collapse
Affiliation(s)
- Thomas A HÖft
- Department of Mathematics, University of St. Thomas, Saint Paul, MN 55105
| | - Bradley K Alpert
- National Institute of Standards and Technology, Boulder, Colorado 80305
| |
Collapse
|
250
|
Abstract
Antibiotic resistance is a prevalent problem in public health worldwide. In general, the carbapenem β-lactam antibiotics are considered a final resort against lethal infections by multidrug-resistant bacteria. Colistin is a cationic polypeptide antibiotic and acts as the last line of defense for treatment of carbapenem-resistant bacteria. Very recently, a new plasmid-borne colistin resistance gene, mcr-2, was revealed soon after the discovery of the paradigm gene mcr-1, which has disseminated globally. However, the molecular mechanisms for MCR-2 colistin resistance are poorly understood. Here we show a unique transposon unit that facilitates the acquisition and transfer of mcr-2 Evolutionary analyses suggested that both MCR-2 and MCR-1 might be traced to their cousin phosphoethanolamine (PEA) lipid A transferase from a known polymyxin producer, Paenibacillus Transcriptional analyses showed that the level of mcr-2 transcripts is relatively higher than that of mcr-1 Genetic deletions revealed that the transmembrane regions (TM1 and TM2) of both MCR-1 and MCR-2 are critical for their location and function in bacterial periplasm, and domain swapping indicated that the TM2 is more efficient than TM1. Matrix-assisted laser desorption ionization-time of flight mass spectrometry (MALDI-TOF MS) confirmed that all four MCR proteins (MCR-1, MCR-2, and two chimeric versions [TM1-MCR-2 and TM2-MCR-1]) can catalyze chemical modification of lipid A moiety anchored on lipopolysaccharide (LPS) with the addition of phosphoethanolamine to the phosphate group at the 4' position of the sugar. Structure-guided site-directed mutagenesis defined an essential 6-residue-requiring zinc-binding/catalytic motif for MCR-2 colistin resistance. The results further our mechanistic understanding of transferable colistin resistance, providing clues to improve clinical therapeutics targeting severe infections by MCR-2-containing pathogens.IMPORTANCE Carbapenem and colistin are the last line of refuge in fighting multidrug-resistant Gram-negative pathogens. MCR-2 is a newly emerging variant of the mobilized colistin resistance protein MCR-1, posing a potential challenge to public health. Here we report transfer of the mcr-2 gene by a unique transposal event and its possible origin. Distribution of MCR-2 in bacterial periplasm is proposed to be a prerequisite for its role in the context of biochemistry and the colistin resistance. We also define the genetic requirement of a zinc-binding/catalytic motif for MCR-2 colistin resistance. This represents a glimpse of transferable colistin resistance by MCR-2.
Collapse
|