1
|
Basit A, Choudhury D, Bandyopadhyay P. Prediction of Ca 2+ Binding Site in Proteins With a Fast and Accurate Method Based on Statistical Mechanics and Analysis of Crystal Structures. Proteins 2024. [PMID: 39258438 DOI: 10.1002/prot.26743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/21/2024] [Revised: 08/20/2024] [Accepted: 08/26/2024] [Indexed: 09/12/2024]
Abstract
Predicting the precise locations of metal binding sites within metalloproteins is a crucial challenge in biophysics. A fast, accurate, and interpretable computational prediction method can complement the experimental studies. In the current work, we have developed a method to predict the location of Ca2+ ions in calcium-binding proteins using a physics-based method with an all-atom description of the proteins, which is substantially faster than the molecular dynamics simulation-based methods with accuracy as good as data-driven approaches. Our methodology uses the three-dimensional reference interaction site model (3D-RISM), a statistical mechanical theory, to calculate Ca2+ ion density around protein structures, and the locations of the Ca2+ ions are obtained from the density. We have taken previously used datasets to assess the efficacy of our method as compared to previous works. Our accuracy is 88%, comparable with the FEATURE program, one of the well-known data-driven methods. Moreover, our method is physical, and the reasons for failures can be ascertained in most cases. We have thoroughly examined the failed cases using different structural and crystallographic measures, such as B-factor, R-factor, electron density map, and geometry at the binding site. It has been found that x-ray structures have issues in many of the failed cases, such as geometric irregularities and dubious assignment of ion positions. Our algorithm, along with the checks for structural accuracy, is a major step in predicting calcium ion positions in metalloproteins.
Collapse
Affiliation(s)
- Abdul Basit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| | | | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi, India
| |
Collapse
|
2
|
Zamanos A, Ioannakis G, Emiris IZ. HydraProt: A New Deep Learning Tool for Fast and Accurate Prediction of Water Molecule Positions for Protein Structures. J Chem Inf Model 2024; 64:2594-2611. [PMID: 38552195 PMCID: PMC11005053 DOI: 10.1021/acs.jcim.3c01559] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2023] [Revised: 02/13/2024] [Accepted: 02/13/2024] [Indexed: 04/09/2024]
Abstract
Water molecules are integral to the structural stability of proteins and vital for facilitating molecular interactions. However, accurately predicting their precise position around protein structures remains a significant challenge, making it a vibrant research area. In this paper, we introduce HydraProt (deep Hydration of Proteins), a novel methodology for predicting precise positions of water molecule oxygen atoms around protein structures, leveraging two interconnected deep learning architectures: a 3D U-net and a Multi-Layer Perceptron (MLP). Our approach starts by introducing a coarse voxel-based representation of the protein, which allows for rapid sampling of candidate water positions via the 3D U-net. These water positions are then assessed by embedding the water-protein relationship in the Euclidean space by means of an MLP. Finally, a postprocessing step is applied to further refine the MLP predictions. HydraProt surpasses existing state-of-the-art approaches in terms of precision and recall and has been validated on large data sets of protein structures. Notably, our method offers rapid inference runtime and should constitute the method of choice for protein structure studies and drug discovery applications. Our pretrained models, data, and the source code required to reproduce these results are accessible at https://github.com/azamanos/HydraProt.
Collapse
Affiliation(s)
- Andreas Zamanos
- Archimedes, Athena Research Center, Marousi 15125, Greece
- Department
of Informatics and Telecommunications, National
and Kapodistrian University of Athens, Athens 16122, Greece
| | - George Ioannakis
- Institute
for Language and Speech Processing, Athena
Research Center, Xanthi 67100, Greece
| | - Ioannis Z. Emiris
- Department
of Informatics and Telecommunications, National
and Kapodistrian University of Athens, Athens 16122, Greece
- Athena
Research Center, Marousi 15125, Greece
| |
Collapse
|
3
|
Magyar ZÉ, Bauer J, Bauerová-Hlinková V, Jóna I, Gaburjakova J, Gaburjakova M, Almássy J. Eu 3+ detects two functionally distinct luminal Ca 2+ binding sites in ryanodine receptors. Biophys J 2023; 122:3516-3531. [PMID: 37533257 PMCID: PMC10502479 DOI: 10.1016/j.bpj.2023.07.029] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Revised: 06/26/2023] [Accepted: 07/31/2023] [Indexed: 08/04/2023] Open
Abstract
Ryanodine receptors (RyRs) are Ca2+ release channels, gated by Ca2+ in the cytosol and the sarcoplasmic reticulum lumen. Their regulation is impaired in certain cardiac and muscle diseases. Although a lot of data is available on the luminal Ca2+ regulation of RyR, its interpretation is complicated by the possibility that the divalent ions used to probe the luminal binding sites may contaminate the cytoplasmic sites by crossing the channel pore. In this study, we used Eu3+, an impermeable agonist of Ca2+ binding sites, as a probe to avoid this complication and to gain more specific information about the function of the luminal Ca2+ sensor. Single-channel currents were measured from skeletal muscle and cardiac RyRs (RyR1 and RyR2) using the lipid bilayer technique. We show that RyR2 is activated by the luminal addition of Ca2+, whereas RyR1 is inhibited. These results were qualitatively reproducible using Eu3+. The luminal regulation of RyR1 carrying a mutation associated with malignant hyperthermia was not different from that of the wild-type. RyR1 inhibition by Eu3+ was extremely voltage dependent, whereas RyR2 activation did not depend on the membrane potential. These results suggest that the RyR1 inhibition site is in the membrane's electric field (channel pore), whereas the RyR2 activation site is outside. Using in silico analysis and previous results, we predicted putative Ca2+ binding site sequences. We propose that RyR2 bears an activation site, which is missing in RyR1, but both isoforms share the same inhibitory Ca2+ binding site near the channel gate.
Collapse
Affiliation(s)
- Zsuzsanna É Magyar
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Jacob Bauer
- Institute of Molecular Biology, Slovak Academy of Sciences, Bratislava, Slovakia
| | | | - István Jóna
- Department of Physiology, Faculty of Medicine, University of Debrecen, Debrecen, Hungary
| | - Jana Gaburjakova
- Institute of Molecular Physiology and Genetics, Centre of Biosciences, Slovak Academy of Sciences, Bratislava, Slovakia
| | - Marta Gaburjakova
- Institute of Molecular Physiology and Genetics, Centre of Biosciences, Slovak Academy of Sciences, Bratislava, Slovakia
| | - János Almássy
- Department of Physiology, Semmelweis University, Budapest, Hungary.
| |
Collapse
|
4
|
Wang Y, Zhang D, Huang L, Zhang Z, Shi Q, Hu J, He G, Guo X, Shi H, Liang L. Uncovering the interactions between PME and PMEI at the gene and protein levels: Implications for the design of specific PMEI. J Mol Model 2023; 29:286. [PMID: 37610510 DOI: 10.1007/s00894-023-05644-y] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2023] [Accepted: 06/30/2023] [Indexed: 08/24/2023]
Abstract
CONTEXT Pectin methylesterase inhibitor (PMEI) can specifically bind and inhibit the activity of pectin methylesterase (PME), which has been widely used in fruit and vegetable juice processing. However, the limited three-dimensional structure, unclear action mechanism, low thermal stability and biological activity of PMEI severely limited its application. In this work, molecular recognition and conformational changes of PME and PMEI were analyzed by various molecular simulation methods. Then suggestions were proposed for improving thermal stability and affinity maturation of PMEI through semi-rational design. METHODS Phylogenetic trees of PME and PMEI were established using the Maximum likelihood (ML) method. The results show that PME and PMEI have good sequence and structure conservation in various plants, and the simulated data can be widely adopted. In this work, MD simulations were performed using AMBER20 package and ff14SB force field. Protein interaction analysis indicates that H-bonds, van der Waals forces, and the salt bridge formed of K224 with ID116 are the main driving forces for mutual molecular recognition of PME and PMEI. According to the analyses of free energy landscape (FEL), conformational cluster, and motion, the association with PMEI greatly disrupts PME's dispersed functional motion mode and biological function. By monitoring the changes of residue contact number and binding free energy, IG35M/ IG35R: IT93F and IT113W/ IT113W: ID116W mutations contribute to thermal stability and affinity maturation of the PME-PMEI complex system, respectively. This work reveals the interaction between PME and PMEI at the gene and protein levels and provides options for modifying specific PMEI.
Collapse
Affiliation(s)
- Yueteng Wang
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Derong Zhang
- School of Marxism, Chengdu Vocational & Technical College of Industry, Chengdu, 610081, China
| | - Lifen Huang
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Zelan Zhang
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Quanshan Shi
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Jianping Hu
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Gang He
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Xiaoqiang Guo
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China
| | - Hang Shi
- Institute of Bioinformatics and Medical Engineering, School of Electrical and Information Engineering, Jiangsu University of Technology, Changzhou, 213001, China.
| | - Li Liang
- Key Laboratory of Medicinal and Edible Plants Resources Development of Sichuan Education Department, School of Pharmacy, Chengdu University, Chengdu, 610106, China.
| |
Collapse
|
5
|
Naz S, Liu P, Farooq U, Ma H. Insight into de-regulation of amino acid feedback inhibition: a focus on structure analysis method. Microb Cell Fact 2023; 22:161. [PMID: 37612753 PMCID: PMC10464499 DOI: 10.1186/s12934-023-02178-z] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2023] [Accepted: 08/13/2023] [Indexed: 08/25/2023] Open
Abstract
Regulation of amino acid's biosynthetic pathway is of significant importance to maintain homeostasis and cell functions. Amino acids regulate their biosynthetic pathway by end-product feedback inhibition of enzymes catalyzing committed steps of a pathway. Discovery of new feedback resistant enzyme variants to enhance industrial production of amino acids is a key objective in industrial biotechnology. Deregulation of feedback inhibition has been achieved for various enzymes using in vitro and in silico mutagenesis techniques. As enzyme's function, its substrate binding capacity, catalysis activity, regulation and stability are dependent on its structural characteristics, here, we provide detailed structural analysis of all feedback sensitive enzyme targets in amino acid biosynthetic pathways. Current review summarizes information regarding structural characteristics of various enzyme targets and effect of mutations on their structures and functions especially in terms of deregulation of feedback inhibition. Furthermore, applicability of various experimental as well as computational mutagenesis techniques to accomplish feedback resistance has also been discussed in detail to have an insight into various aspects of research work reported in this particular field of study.
Collapse
Affiliation(s)
- Sadia Naz
- Biodesign Center, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
| | - Pi Liu
- Biodesign Center, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
| | - Umar Farooq
- Department of Chemistry, COMSATS University Islamabad, Abbottabad Campus, Islamabad, 22060, Pakistan
| | - Hongwu Ma
- Biodesign Center, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China.
| |
Collapse
|
6
|
Zsidó BZ, Bayarsaikhan B, Börzsei R, Szél V, Mohos V, Hetényi C. The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering. Int J Mol Sci 2023; 24:11784. [PMID: 37511543 PMCID: PMC10381018 DOI: 10.3390/ijms241411784] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2023] [Revised: 07/18/2023] [Accepted: 07/20/2023] [Indexed: 07/30/2023] Open
Abstract
Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.
Collapse
Affiliation(s)
- Balázs Zoltán Zsidó
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Bayartsetseg Bayarsaikhan
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Rita Börzsei
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Viktor Szél
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Violetta Mohos
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| | - Csaba Hetényi
- Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, 7624 Pécs, Hungary
| |
Collapse
|
7
|
Zan Q, Long M, Zheng N, Zhang Z, Zhou H, Xu X, Osire T, Xia X. Improving ethanol tolerance of ethyl carbamate hydrolase by diphasic high pressure molecular dynamic simulations. AMB Express 2023; 13:32. [PMID: 36920541 PMCID: PMC10017909 DOI: 10.1186/s13568-023-01538-7] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/01/2023] [Accepted: 03/04/2023] [Indexed: 03/16/2023] Open
Abstract
Ethyl carbamate (EC) is mainly found in fermented foods and fermented alcoholic beverages, which could cause carcinogenic potential to humans. Reducing EC is one of the key research priorities to address security of fermented foods. Enzymatic degradation of EC with EC hydrolase in food is the most reliable and efficient method. However, poor tolerance to ethanol severely hinders application of EC hydrolase. In this study, the mutants of EC hydrolase were screened by diphasic high pressure molecular dynamic simulations (dHP-MD). The best variant with remarkable improvement in specific activity and was H68A/K70R/S325N, whose specific activity was approximately 3.42-fold higher than WT, and relative enzyme activity under 20% (v/v) was 5.02-fold higher than WT. Moreover, the triple mutant increased its stability by acquiring more hydration shell and forming extra hydrogen bonds. Furthermore, the ability of degrading EC of the immobilized triple mutant was both detected in mock wine and under certain reaction conditions. The stability of immobilized triple mutant and WT were both improved, and immobilized triple mutant degraded nearly twice as much EC as that of immobilized WT. Overall, dHP-MD was proved to effectively improve enzyme activity and ethanol tolerance for extent application at industrial scale.
Collapse
Affiliation(s)
- Qijia Zan
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Mengfei Long
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Nan Zheng
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Zehua Zhang
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Huimin Zhou
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Xinjie Xu
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China
| | - Tolbert Osire
- Faculty of Biology, Shenzhen MSU-BIT University, Shenzhen, 518172, Guangdong, China
| | - Xiaole Xia
- Key Laboratory of Industrial Biotechnology, Ministry of Education, School of Biotechnology, Jiangnan University, Wuxi, 214122, Jiangsu, China.
| |
Collapse
|
8
|
Zhang J, Zhou F, Liang X, Yang G. SCAMPER: Accurate Type-Specific Prediction of Calcium-Binding Residues Using Sequence-Derived Features. IEEE/ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS 2023; 20:1406-1416. [PMID: 35536812 DOI: 10.1109/tcbb.2022.3173437] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/04/2023]
Abstract
Understanding molecular mechanisms involved in calcium-protein interactions and modeling corresponding docking rely on the accurate identification of calcium-binding residues (CaBRs). The defects of experimentally annotating protein functions enhances the development of computational approaches that correctly identify calcium-binding interactions. Studies have reported that current methods severely cross-predict residues that interact with other types of molecules (e.g., nucleic acids, proteins, and small ligands) as CaBRs. In this study, a novel predictor named SCAMPER (Selective CAlciuM-binding PrEdictoR) is proposed for the accurate and specific prediction of CaBRs. SCAMPER is designed using newly compiled dataset with complete UniProt sequences and annotations, which include calcium-binding, nucleic acid-binding, protein-binding, and small ligand-binding residues. We use a novel designed two-layer scheme to perform predictions as well as penalize cross-predictions. Empirical tests on an independent test dataset reveals that the proposed method significantly outperforms state-of-the-art predictors. SCAMPER is proved to be capable of distinguishing CaBRs from different types of metal-ion binding residues. We further perform CaBRs predictions on the whole human proteome, and use the results to hypothesize calcium-binding proteins (CaBPs). The latest experimental verified CaBPs and GO analysis prove the accuracy of our predictions. We implement the proposed method and share the data at http://www.inforstation.com/webservers/SCAMPER/.
Collapse
|
9
|
Kim TH, Zhou Z, Choi YS, Costanza V, Wang L, Bahng JH, Higdon NJ, Yun Y, Kang H, Kim S, Daraio C. Flexible biomimetic block copolymer composite for temperature and long-wave infrared sensing. SCIENCE ADVANCES 2023; 9:eade0423. [PMID: 36763652 PMCID: PMC9916982 DOI: 10.1126/sciadv.ade0423] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/22/2022] [Accepted: 01/09/2023] [Indexed: 06/18/2023]
Abstract
Biological compounds often provide clues to advance material designs. Replicating their molecular structure and functional motifs in artificial materials offers a blueprint for unprecedented functionalities. Here, we report a flexible biomimetic thermal sensing (BTS) polymer that is designed to emulate the ion transport dynamics of a plant cell wall component, pectin. Using a simple yet versatile synthetic procedure, we engineer the physicochemical properties of the polymer by inserting elastic fragments in a block copolymer architecture, making it flexible and stretchable. The thermal response of our flexible polymer outperforms current state-of-the-art temperature sensing materials, including vanadium oxide, by up to two orders of magnitude. Thermal sensors fabricated from these composites exhibit a sensitivity that exceeds 10 mK and operate stably between 15° and 55°C, even under repeated mechanical deformations. We demonstrate the use of our flexible BTS polymer in two-dimensional arrays for spatiotemporal temperature mapping and broadband infrared photodetection.
Collapse
Affiliation(s)
- Tae Hyun Kim
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| | - Zhun Zhou
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| | - Yeong Suk Choi
- Samsung Advanced Institute of Technology (SAIT), Samsung Electronics, Suwon 16678, South Korea
| | - Vincenzo Costanza
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| | - Linghui Wang
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| | - Joong Hwan Bahng
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| | - Nicholas J. Higdon
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
- Division of Chemistry and Chemical Engineering, California Institute of Technology, Pasadena, CA 91125, USA
| | - Youngjun Yun
- Samsung Advanced Institute of Technology (SAIT), Samsung Electronics, Suwon 16678, South Korea
| | - Hyunbum Kang
- Samsung Advanced Institute of Technology (SAIT), Samsung Electronics, Suwon 16678, South Korea
| | - Sunghan Kim
- Samsung Advanced Institute of Technology (SAIT), Samsung Electronics, Suwon 16678, South Korea
| | - Chiara Daraio
- Division of Engineering and Applied Science, California Institute of Technology, Pasadena, CA 91125, USA
| |
Collapse
|
10
|
Wang K. GPDOCK: highly accurate docking strategy for metalloproteins based on geometric probability. Brief Bioinform 2023; 24:6987821. [PMID: 36642411 DOI: 10.1093/bib/bbac620] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2022] [Revised: 12/09/2022] [Accepted: 12/17/2022] [Indexed: 01/17/2023] Open
Abstract
Accurately predicting the interaction modes for metalloproteins remains extremely challenging in structure-based drug design and mechanism analysis of enzymatic catalysis due to the complexity of metal coordination in metalloproteins. Here, we report a docking method for metalloproteins based on geometric probability (GPDOCK) with unprecedented accuracy. The docking tests of 10 common metal ions with 9360 metalloprotein-ligand complexes demonstrate that GPDOCK has an accuracy of 94.3% in predicting binding pose. What is more, it can accurately realize the docking of metalloproteins with ligand when one or two water molecules are engaged in the metal ion coordination. Since GPDOCK only depends on the three-dimensional structure of metalloprotein and ligand, structure-based machine learning model is employed for the scoring of binding poses, which significantly improves computational efficiency. The proposed docking strategy can be an effective and efficient tool for drug design and further study of binding mechanism of metalloproteins. The manual of GPDOCK and the code for the logistical regression model used to re-rank the docking results are available at https://github.com/wangkai-zhku/GPDOCK.git.
Collapse
Affiliation(s)
- Kai Wang
- School of Agriculture and Biology, Zhongkai University of Agriculture and Engineering, Guangzhou 510225, P. R. China.,Abinitio Technology Company, Ltd, Guangzhou 510640, P. R. China
| |
Collapse
|
11
|
Antibody CDR amino acids underlying the functionality of antibody repertoires in recognizing diverse protein antigens. Sci Rep 2022; 12:12555. [PMID: 35869245 PMCID: PMC9307644 DOI: 10.1038/s41598-022-16841-9] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/14/2022] [Accepted: 07/18/2022] [Indexed: 11/21/2022] Open
Abstract
Antibodies recognize protein antigens with exquisite specificity in a complex aqueous environment, where interfacial waters are an integral part of the antibody–protein complex interfaces. In this work, we elucidate, with computational analyses, the principles governing the antibodies’ specificity and affinity towards their cognate protein antigens in the presence of explicit interfacial waters. Experimentally, in four model antibody–protein complexes, we compared the contributions of the interaction types in antibody–protein antigen complex interfaces with the antibody variants selected from phage-displayed synthetic antibody libraries. Evidently, the specific interactions involving a subset of aromatic CDR (complementarity determining region) residues largely form the predominant determinant underlying the specificity of the antibody–protein complexes in nature. The interfacial direct/water-mediated hydrogen bonds accompanying the CDR aromatic interactions are optimized locally but contribute little in determining the epitope location. The results provide insights into the phenomenon that natural antibodies with limited sequence and structural variations in an antibody repertoire can recognize seemingly unlimited protein antigens. Our work suggests guidelines in designing functional artificial antibody repertoires with practical applications in developing novel antibody-based therapeutics and diagnostics for treating and preventing human diseases.
Collapse
|
12
|
Xi Y, Chalvon V, Padilla A, Cesari S, Kroj T. The activity of the RGA5 sensor NLR from rice requires binding of its integrated HMA domain to effectors but not HMA domain self-interaction. MOLECULAR PLANT PATHOLOGY 2022; 23:1320-1330. [PMID: 35766176 PMCID: PMC9366066 DOI: 10.1111/mpp.13236] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2022] [Revised: 05/08/2022] [Accepted: 05/09/2022] [Indexed: 05/25/2023]
Abstract
The rice nucleotide-binding (NB) and leucine-rich repeat (LRR) domain immune receptors (NLRs) RGA4 and RGA5 form a helper NLR/sensor NLR (hNLR/sNLR) pair that specifically recognizes the effectors AVR-Pia and AVR1-CO39 from the blast fungus Magnaporthe oryzae. While RGA4 contains only canonical NLR domains, RGA5 has an additional unconventional heavy metal-associated (HMA) domain integrated after its LRR domain. This RGA5HMA domain binds the effectors and is crucial for their recognition. Investigation of the three-dimensional structure of the AVR1-CO39/RGA5HMA complex by X-ray crystallography identified a candidate surface for effector binding in the HMA domain and showed that the HMA domain self-interacts in the absence of effector through the same surface. Here, we investigated the relevance of this HMA homodimerization for RGA5 function and the role of the RGA5HMA effector-binding and self-interaction surface in effector recognition. By analysing structure-informed point mutations in the RGA5HMA -binding surface in protein interaction studies and in Nicotiana benthamiana cell death assays, we found that HMA self-interaction does not contribute to RGA5 function. However, the effector-binding surface of RGA5HMA identified by X-ray crystallography is crucial for both in vitro and in vivo effector binding as well as effector recognition. These results support the current hypothesis that noncanonical integrated domains of NLRs act primarily as effector traps and deepen our understanding of the sNLRs' function within NLR pairs.
Collapse
Affiliation(s)
- Yuxuan Xi
- PHIM Plant Health Institute, Univ Montpellier, INRAE, CIRAD, Institut Agro, IRDMontpellierFrance
| | - Véronique Chalvon
- PHIM Plant Health Institute, Univ Montpellier, INRAE, CIRAD, Institut Agro, IRDMontpellierFrance
| | - André Padilla
- CBS, Univ Montpellier, CNRS, INSERMMontpellierFrance
| | - Stella Cesari
- PHIM Plant Health Institute, Univ Montpellier, INRAE, CIRAD, Institut Agro, IRDMontpellierFrance
| | - Thomas Kroj
- PHIM Plant Health Institute, Univ Montpellier, INRAE, CIRAD, Institut Agro, IRDMontpellierFrance
| |
Collapse
|
13
|
Heterotropic roles of divalent cations in the establishment of allostery and affinity maturation of integrin αXβ2. Cell Rep 2022; 40:111254. [PMID: 36001965 PMCID: PMC9440770 DOI: 10.1016/j.celrep.2022.111254] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/17/2021] [Revised: 05/23/2022] [Accepted: 08/02/2022] [Indexed: 11/21/2022] Open
Abstract
Allosteric activation and silencing of leukocyte β2-integrins transpire through cation-dependent structural changes, which mediate integrin biosynthesis and recycling, and are essential to designing leukocyte-specific drugs. Stepwise addition of Mg2+ reveals two mutually coupled events for the αXβ2 ligand-binding domain-the αX I-domain-corresponding to allostery establishment and affinity maturation. Electrostatic alterations in the Mg2+-binding site establish long-range couplings, leading to both pH- and Mg2+-occupancy-dependent biphasic stability change in the αX I-domain fold. The ligand-binding sensorgrams show composite affinity events for the αX I-domain accounting for the multiplicity of the αX I-domain conformational states existing in the solution. On cell surfaces, increasing Mg2+ concentration enhanced adhesiveness of αXβ2. This work highlights how intrinsically flexible pH- and cation-sensitive architecture endows a unique dynamic continuum to the αI-domain structure on the intact integrin, thereby revealing the importance of allostery establishment and affinity maturation in both extracellular and intracellular integrin events.
Collapse
|
14
|
Selvam P, Jain A, Abbott J, Ahuja AS, Cheema A, Bruno KA, Atwal H, Forghani I, Caulfield T, Atwal PS. Molecular Modeling and Phenotypic Description of a Patient with a Novel Exonic Deletion of GALNS with Resultant Morquio Syndrome with Two Successful Pregnancies. Mol Syndromol 2022; 13:282-289. [PMID: 36158053 PMCID: PMC9421684 DOI: 10.1159/000519326] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/07/2020] [Accepted: 08/30/2021] [Indexed: 01/03/2023] Open
Abstract
In this report, we describe phenotypic features of a patient with mucopolysaccharidosis type IVA (Morquio syndrome) harboring a novel exon 1 deletion in GALNS with enzymatic confirmation consistent with Morquio syndrome. To our knowledge, this is the first reported case of this variant. Additionally, we protein modelled wild-type GALNS and the pathogenic variant with an exon 1 deletion for comparative analysis using statistical mechanics methods described herein. We demonstrate that, even when the protein is translated, the mutation would affect protein stability and function via homodimer interaction modifications. Lastly, given the patient's 2 successful pregnancies, data about the management of pregnancies in mucopolysaccharidoses are reviewed, and we discuss the management of pregnancy in patients with Morquio syndrome.
Collapse
Affiliation(s)
| | | | | | - Abhimanyu S. Ahuja
- Charles E. Schmidt College of Medicine, Florida Atlantic University, Boca Raton, Florida, USA
| | | | | | | | - Irman Forghani
- Dr. John T. Macdonald Foundation Department of Human Genetics, Miller School of Medicine, University of Miami, Miami, Florida, USA
| | | | | |
Collapse
|
15
|
Roy P, Bhattacharyya D. MetBP: A Software Tool for Detection of Interaction between Metal Ion-RNA Base Pairs. Bioinformatics 2022; 38:3833-3834. [PMID: 35695777 DOI: 10.1093/bioinformatics/btac392] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/17/2022] [Revised: 05/09/2022] [Accepted: 06/10/2022] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION The role of metals in shaping and functioning of RNA is a well established fact and the understanding of that through the analysis of structural data has biological relevance. Often metal ions bind to one or more atoms of the nucleobase of an RNA. This fact becomes more interesting when such bases form a base pair with any other base. Furthermore, when metal ions bind to any residue of an RNA, the secondary structural features of the residue (helix, loop, unpaired etc) are also biologically important. The available metal binding related software tools cannot address such type specific queries. RESULTS To fill this limitation, we have designed a software tool, called MetBP, that meets the goal. This tool is a stand-alone command line based tool and has no dependency on the other existing software. It accepts a structure file in mmCIF or PDB format and computes the base pairs and thereafter reports all metals that bind to one or more nucleotides that form pairs with another. It reports binding distance, angles along with base pair stability. It also reports several other important aspects, e.g. secondary structure of the residue in the RNA. MetBP can be used as a generalized metal binding site detection tool for Proteins and DNA as well. AVAILABILITY https://github.com/computational-biology/metbp. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Parthajit Roy
- The Department of Computer Science, The University of Burdwan, Burdwan 713104, West Bengal, India
| | | |
Collapse
|
16
|
Smith ST, Shub L, Meiler J. PlaceWaters: Real-time, explicit interface water sampling during Rosetta ligand docking. PLoS One 2022; 17:e0269072. [PMID: 35639743 PMCID: PMC9154094 DOI: 10.1371/journal.pone.0269072] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/02/2022] [Accepted: 05/13/2022] [Indexed: 01/29/2023] Open
Abstract
Water molecules at the protein-small molecule interface often form hydrogen bonds with both the small molecule ligand and the protein, affecting the structural integrity and energetics of a binding event. The inclusion of these 'bridging waters' has been shown to improve the accuracy of predicted docked structures; however, due to increased computational costs, this step is typically omitted in ligand docking simulations. In this study, we introduce a resource-efficient, Rosetta-based protocol named "PlaceWaters" to predict the location of explicit interface bridging waters during a ligand docking simulation. In contrast to other explicit water methods, this protocol is independent of knowledge of number and location of crystallographic waters in homologous structures. We test this method on a diverse protein-small molecule benchmark set in comparison to other Rosetta-based protocols. Our results suggest that this coarse-grained, structure-based approach quickly and accurately predicts the location of bridging waters, improving our ability to computationally screen drug candidates.
Collapse
Affiliation(s)
- Shannon T. Smith
- Chemical and Physical Biology Program, Vanderbilt University, Nashville, Tennessee, United States of America
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
| | - Laura Shub
- Biomedical Informatics Program, University of California San Francisco, San Francisco, California, United States of America
- Institute for Neurodegenerative Diseases, University of California San Francisco, San Francisco, California, United States of America
| | - Jens Meiler
- Center for Structural Biology, Vanderbilt University, Nashville, Tennessee, United States of America
- Departments of Chemistry, Pharmacology, and Biomedical Informatics, Center for Structural Biology and Institute of Chemical Biology, Nashville, Tennessee, United States of America
- Institute for Drug Discovery, Leipzig University Medical School, SAC, Leipzig, Germany
- * E-mail:
| |
Collapse
|
17
|
Basit A, Yadav AK, Bandyopadhyay P. Calcium Ion Binding to the Mutants of Calmodulin: A Structure-Based Computational Predictive Model of Binding Affinity Using a Charge Scaling Approach in Molecular Dynamics Simulation. J Chem Inf Model 2022; 62:2821-2834. [PMID: 35608259 DOI: 10.1021/acs.jcim.2c00428] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The binding of calcium ions (Ca2+) to the calcium-binding proteins (CBPs) controls a plethora of regulatory processes. Among the roles played by CBPs in several diseases, the onset and progress of some cardiovascular diseases are caused by mutations in calmodulin (CaM), an important member of CBPs. Rationalization and prediction of the binding affinity of Ca2+ ions to the CaM can play important roles in understanding the origin of cardiovascular diseases. However, there is no robust structure-based computational method for predicting the binding affinity of Ca2+ ions to the different forms of CBPs in general and CaM in particular. In the current work, we have devised a fast yet accurate computational technique to accurately calculate the binding affinity of Ca2+ to the different forms of CaM. This method combines the well-known molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method and a charge scaling approach developed by previous authors that takes care of the polarization of CaM and Ca2+ ions. Our detailed analysis of the different components of binding free energy shows that subtle changes in electrostatics and van der Waals contribute to the difference in the binding affinity of mutants from that of the wild type (WT), and the charge scaling approach is superior in calculating these subtle changes in electrostatics as compared to the nonpolarizable force field used in this work. A statistically significant regression model made from our binding free energy calculations gives a correlation coefficient close to 0.8 to the experimental results. This structure-based predictive model can open up a new strategy to understand and predict the binding of Ca2+ to the mutants of CBPs, in general.
Collapse
Affiliation(s)
- Abdul Basit
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Ajeet Kumar Yadav
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| | - Pradipta Bandyopadhyay
- School of Computational and Integrative Sciences, Jawaharlal Nehru University, New Delhi 110067, India
| |
Collapse
|
18
|
Structural insights into collagen binding by platelet receptor glycoprotein VI. Blood 2022; 139:3087-3098. [PMID: 35245360 DOI: 10.1182/blood.2021013614] [Citation(s) in RCA: 19] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/04/2021] [Accepted: 02/08/2022] [Indexed: 02/02/2023] Open
Abstract
Glycoprotein VI (GPVI) mediates collagen-induced platelet activation after vascular damage and is an important contributor to the onset of thrombosis, heart attack, and stroke. Animal models of thrombosis have identified GPVI as a promising target for antithrombotic therapy. Although for many years the crystal structure of GPVI has been known, the essential details of its interaction with collagen have remained elusive. Here, we present crystal structures of the GPVI ectodomain bound to triple-helical collagen peptides, which reveal a collagen-binding site across the β-sheet of the D1 domain. Mutagenesis and binding studies confirm the observed binding site and identify Trp76, Arg38, and Glu40 as essential residues for binding to fibrillar collagens and collagen-related peptides (CRPs). GPVI binds a site on collagen comprising two collagen chains with the core formed by the sequence motif OGPOGP. Potent GPVI-binding peptides from Toolkit-III all contain OGPOGP; weaker binding peptides frequently contain a partial motif varying at either terminus. Alanine-scanning of peptide III-30 also identified two AGPOGP motifs that contribute to GPVI binding, but steric hindrance between GPVI molecules restricts the maximum binding capacity. We further show that no cooperative interactions could occur between two GPVI monomers binding to a stretch of (GPO)5 and that binding of ≥2 GPVI molecules to a fibril-embedded helix requires non-overlapping OGPOGP motifs. Our structure confirms the previously suggested similarity in collagen binding between GPVI and leukocyte-associated immunoglobulin-like receptor 1 (LAIR-1) but also indicates significant differences that may be exploited for the development of receptor-specific therapeutics.
Collapse
|
19
|
van der Kant R, Louros N, Schymkowitz J, Rousseau F. Thermodynamic analysis of amyloid fibril structures reveals a common framework for stability in amyloid polymorphs. Structure 2022; 30:1178-1189.e3. [DOI: 10.1016/j.str.2022.05.002] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/20/2021] [Revised: 03/30/2022] [Accepted: 04/29/2022] [Indexed: 11/28/2022]
|
20
|
Humanizing plant-derived snakins and their encrypted antimicrobial peptides. Biochimie 2022; 199:92-111. [PMID: 35472564 DOI: 10.1016/j.biochi.2022.04.011] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/13/2021] [Revised: 04/16/2022] [Accepted: 04/20/2022] [Indexed: 12/11/2022]
Abstract
Due to safety restrictions, plant-derived antimicrobial peptides (AMPs) need optimization to be consumed beyond preservatives. Herein, 175 GASA-domain-containing snakins were analyzed. Factors including charge, hydrophobicity, helicity, hydrophobic moment (μH), folding enthalpy, folding heat capacity, folding free energy, therapeutic index, allergenicity, and bitterness were considered. The most optimal snakins for oral consumption as preservatives were from Cajanus cajan, Cucumis melo, Durio zibethinus, Glycine soja, Herrania umbratica, and Ziziphus jujuba. Virtual digestion of snakins predicted ACE1 and DPPIV inhibitory as dominant effects upon oral use with antihypertensive and antidiabetic properties. To be applied as a therapeutic in parenteral administration, snakins were browsed for short 20-mer encrypted fragments that were non-toxic or with eliminated toxicity using directed mutagenesis yet retaining the AMP property. The most promising 20-mer AMPs were Mr-SNK2-1a in Morella rubra with BBB permeation, Na-SNK2-2a(C18W), and Na-SNK2-2b(C16F) from Nicotiana attenuata. These AMPs were cell-penetrating peptides (CPPs), with a charge of +6, a μH of about 0.40, and a Boman-index higher than 2.48 Kcalmol-1. Na-SNK2-2a(C18W) had putative activity against gram-negative bacteria with MIC lower than 25 μgml-1, and Na-SNK2-2b(C16F) was a potential anti-HIV with an IC50 of 3.04 μM. Other 20-mer AMPs, such as Cc-SNK1-2a from Cajanus cajan displayed an anti-HCV property with an IC50 of 13.91 μM. While Si-SNK2-3a(C17P) from Sesamum indicum was a cationic anti-angiogenic CPP targeting the acidic microenvironment of tumors, Cme-SNK2-1a(C11F) from Cucumis melo was an immunomodulator CPP applicable as a vaccine adjuvant. Because of combined mechanisms, investigating cysteine-rich peptides can nominate effective biotherapeutics.
Collapse
|
21
|
Vila JA. Proteins' Evolution upon Point Mutations. ACS OMEGA 2022; 7:14371-14376. [PMID: 35573218 PMCID: PMC9089682 DOI: 10.1021/acsomega.2c01407] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/09/2022] [Accepted: 04/05/2022] [Indexed: 05/03/2023]
Abstract
As the reader must be already aware, state-of-the-art protein folding prediction methods have reached a smashing success in their goal of accurately determining the three-dimensional structures of proteins. Yet, a solution to simple problems such as the effects of protein point mutations on their (i) native conformation; (ii) marginal stability; (iii) ensemble of high-energy nativelike conformations; and (iv) metamorphism propensity and, hence, their evolvability, remains as an unsolved problem. As a plausible solution to the latter, some properties of the amide hydrogen-deuterium exchange, a highly sensitive probe of the structure, stability, and folding of proteins, are assessed from a new perspective. The preliminary results indicate that the protein marginal stability change upon point mutations provides the necessary and sufficient information to estimate, through a Boltzmann factor, the evolution of the amide hydrogen exchange protection factors and, consequently, that of the ensemble of folded conformations coexisting with the native state. This work contributes to our general understanding of the effects of point mutations on proteins and may spur significant progress in our efforts to develop methods to determine the appearance of new folds and functions accurately.
Collapse
|
22
|
Schmalz S, Mayr V, Shosherova A, Gepp B, Ackerbauer D, Sturm G, Bohle B, Breiteneder H, Radauer C. Isotype-specific binding patterns of serum antibodies to multiple conformational epitopes of Bet v 1. J Allergy Clin Immunol 2021; 149:1786-1794.e12. [PMID: 34740603 DOI: 10.1016/j.jaci.2021.10.026] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/12/2021] [Revised: 09/27/2021] [Accepted: 10/08/2021] [Indexed: 12/26/2022]
Abstract
BACKGROUND Birch pollen is an important elicitor of respiratory allergy. The major allergen, Bet v 1, binds IgE exclusively via conformational epitopes. OBJECTIVE To identify Bet v 1-specific epitope repertoires of IgE and IgG from birch pollen-allergic and non-allergic subjects. METHODS Chimeric proteins were created by grafting individual epitope-sized, contiguous surface patches of Bet v 1 onto a non-allergenic structural homologue and expressed in Escherichia coli. Binding of IgE, IgG1 and IgG4 from sera of 30 birch pollen-allergic and 11 non-allergic subjects to Bet v 1, 13 chimeric proteins and four bacterial Bet v 1 homologues were measured by ELISA. The proportion of epitope-specific in total Bet v 1-specific IgE and the cross-reactivity of Bet v 1-specific IgE with bacterial homologues were determined by competitive ELISA. RESULTS Thirteen soluble, correctly folded chimeric proteins were produced. IgE from 27/30 birch pollen-allergic patients bound to 1-12 chimeric proteins (median 4.0) with patient-specific patterns. Three chimeras binding IgE from the majority of sera were identified, whose pgrafted patches overlapped with previously published epitopes. Patterns of IgG1 and IgG4 binding to the chimeric proteins did not correspond to the binding patterns of IgE. Sera of 19/30 birch pollen-allergic patients contained low amounts of IgE to bacterial homologues. Bacterial proteins were able to partially inhibit IgE binding to Bet v 1. CONCLUSION Epitopes recognized by Bet v 1-specific antibodies from birch pollen-allergic patients are specific to each patient and differ between IgE, IgG1 and IgG4.
Collapse
Affiliation(s)
- Stefanie Schmalz
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Vanessa Mayr
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Alexandra Shosherova
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Barbara Gepp
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria; Department Life Science Engineering, University of Applied Sciences Technikum Wien, Vienna Austria
| | - Daniela Ackerbauer
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Gunter Sturm
- Allergy Outpatient Clinic Reumannplatz, Vienna, Austria; Department of Dermatology, Medical University of Graz, Graz, Austria
| | - Barbara Bohle
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Heimo Breiteneder
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria
| | - Christian Radauer
- Department of Pathophysiology and Allergy Research, Center for Pathophysiology, Infectiology and Immunology, Medical University of Vienna, Vienna, Austria.
| |
Collapse
|
23
|
Bonet LFS, Loureiro JP, Pereira GRC, Da Silva ANR, De Mesquita JF. Molecular dynamics and protein frustration analysis of human fused in Sarcoma protein variants in Amyotrophic Lateral Sclerosis type 6: An In Silico approach. PLoS One 2021; 16:e0258061. [PMID: 34587215 PMCID: PMC8480726 DOI: 10.1371/journal.pone.0258061] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/14/2021] [Accepted: 09/16/2021] [Indexed: 11/18/2022] Open
Abstract
Amyotrophic lateral sclerosis (ALS) is the most frequent adult-onset motor neuron disorder. The disease is characterized by degeneration of upper and lower motor neurons, leading to death usually within five years after the onset of symptoms. While most cases are sporadic, 5%-10% of cases can be associated with familial inheritance, including ALS type 6, which is associated with mutations in the Fused in Sarcoma (FUS) gene. This work aimed to evaluate how the most frequent ALS-related mutations in FUS, R521C, R521H, and P525L affect the protein structure and function. We used prediction algorithms to analyze the effects of the non-synonymous single nucleotide polymorphisms and performed evolutionary conservation analysis, protein frustration analysis, and molecular dynamics simulations. Most of the prediction algorithms classified the three mutations as deleterious. All three mutations were predicted to reduce protein stability, especially the mutation R521C, which was also predicted to increase chaperone binding tendency. The protein frustration analysis showed an increase in frustration in the interactions involving the mutated residue 521C. Evolutionary conservation analysis showed that residues 521 and 525 of human FUS are highly conserved sites. The molecular dynamics results indicate that protein stability could be compromised in all three mutations. They also affected the exposed surface area and protein compactness. The analyzed mutations also displayed high flexibility in most residues in all variants, most notably in the interaction site with the nuclear import protein of FUS.
Collapse
Affiliation(s)
- L. F. S. Bonet
- Department of Genetics and Molecular Biology, Laboratory of Bioinformatics and Computational Biology, Federal University of the State of Rio de Janeiro, Rio de Janeiro, Brazil
| | - J. P. Loureiro
- Department of Genetics and Molecular Biology, Laboratory of Bioinformatics and Computational Biology, Federal University of the State of Rio de Janeiro, Rio de Janeiro, Brazil
| | - G. R. C. Pereira
- Department of Genetics and Molecular Biology, Laboratory of Bioinformatics and Computational Biology, Federal University of the State of Rio de Janeiro, Rio de Janeiro, Brazil
| | - A. N. R. Da Silva
- Department of Genetics and Molecular Biology, Laboratory of Bioinformatics and Computational Biology, Federal University of the State of Rio de Janeiro, Rio de Janeiro, Brazil
| | - J. F. De Mesquita
- Department of Genetics and Molecular Biology, Laboratory of Bioinformatics and Computational Biology, Federal University of the State of Rio de Janeiro, Rio de Janeiro, Brazil
| |
Collapse
|
24
|
Samways ML, Taylor RD, Bruce Macdonald HE, Essex JW. Water molecules at protein-drug interfaces: computational prediction and analysis methods. Chem Soc Rev 2021; 50:9104-9120. [PMID: 34184009 DOI: 10.1039/d0cs00151a] [Citation(s) in RCA: 29] [Impact Index Per Article: 9.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/05/2023]
Abstract
The fundamental importance of water molecules at drug-protein interfaces is now widely recognised and a significant feature in structure-based drug design. Experimental methods for analysing the role of water in drug binding have many challenges, including the accurate location of bound water molecules in crystal structures, and problems in resolving specific water contributions to binding thermodynamics. Computational analyses of binding site water molecules provide an alternative, and in principle complete, structural and thermodynamic picture, and their use is now commonplace in the pharmaceutical industry. In this review, we describe the computational methodologies that are available and discuss their strengths and weaknesses. Additionally, we provide a critical analysis of the experimental data used to validate the methods, regarding the type and quality of experimental structural data. We also discuss some of the fundamental difficulties of each method and suggest directions for future study.
Collapse
Affiliation(s)
- Marley L Samways
- School of Chemistry, University of Southampton, Highfield, Southampton SO17 1BJ, UK.
| | | | | | | |
Collapse
|
25
|
Sadeghian I, Hemmati S. Characterization of a Stable Form of Carboxypeptidase G2 (Glucarpidase), a Potential Biobetter Variant, From Acinetobacter sp. 263903-1. Mol Biotechnol 2021; 63:1155-1168. [PMID: 34268672 DOI: 10.1007/s12033-021-00370-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Accepted: 07/08/2021] [Indexed: 01/14/2023]
Abstract
Carboxypeptidase G2 (CPG2) is a bacterial enzyme widely used to detoxify methotrexate (MTX) and in enzyme/prodrug therapy for cancer treatment. However, several drawbacks, such as instability, have limited its efficiency. Herein, we have evaluated the properties of a putative CPG2 from Acinetobacter sp. 263903-1 (AcCPG2). AcCPG2 is compared with a CPG2 derived from Pseudomonas sp. strain RS-16 (PsCPG2), available as an FDA-approved medication called glucarpidase. After modeling AcCPG2 using the I-TASSER program, the refined model was validated by PROCHECK, VERIFY 3D and according to the Z score of the model. Using computational analyses, AcCPG2 displayed higher thermodynamic stability and a lower aggregation propensity than PsCPG2. AcCPG2 showed an optimum pH of 7.5 against MTX and was stable over a pH range of 5-10. AcCPG2 exhibited optimum activity at 50 °C and higher thermal stability at a temperature range of 20-70 °C compared to PsCPG2. The Km value of the purified AcCPG2 toward folate and MTX was 31.36 µM and 44.99 µM, respectively. The Vmax value of AcCPG2 for folate and MTX was 125.80 µmol/min/mg and 48.90 µmol/min/mg, respectively. Accordingly, thermostability and pH versatility makes AcCPG2 a potential biobetter variant for therapeutic applications.
Collapse
Affiliation(s)
- Issa Sadeghian
- Biotechnology Research Center, Shiraz University of Medical Sciences, Shiraz, Iran
- Pharmaceutical Sciences Research Center, Shiraz University of Medical Sciences, Shiraz, Iran
- Department of Pharmaceutical Biotechnology, School of Pharmacy, Shiraz University of Medical Sciences, Shiraz, Iran
| | - Shiva Hemmati
- Biotechnology Research Center, Shiraz University of Medical Sciences, Shiraz, Iran.
- Pharmaceutical Sciences Research Center, Shiraz University of Medical Sciences, Shiraz, Iran.
- Department of Pharmaceutical Biotechnology, School of Pharmacy, Shiraz University of Medical Sciences, Shiraz, Iran.
| |
Collapse
|
26
|
A single residue deletion in the barley HKT1;5 P189 variant restores plasma membrane localisation but not Na + conductance. BIOCHIMICA ET BIOPHYSICA ACTA-BIOMEMBRANES 2021; 1863:183669. [PMID: 34139196 DOI: 10.1016/j.bbamem.2021.183669] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Subscribe] [Scholar Register] [Received: 02/25/2021] [Revised: 06/01/2021] [Accepted: 06/02/2021] [Indexed: 12/12/2022]
Abstract
Leaf Na+ exclusion, mediated by plasma membrane-localised Class 1 High-affinity potassium (K+) Transporters (HKTs), is a key mechanism contributing to salinity tolerance of several major crop plants. We determined previously that the leucine to proline residue substitution at position 189 (L189P) in barley HvHKT1;5 disrupts its characteristic plasma membrane localisation and Na+ conductance. Here, we focus on a surprising observation that a single residue deletion of methionine at position 372 (M372del) within the conserved VMMYL motif in plant HKTs, restores plasma membrane localisation but not Na+ conductance in HvHKT1;5 P189. To clarify why the singular M372 deletion regains plasma membrane localisation, we built 3D models and defined α-helical assembly pathways of the P189 M372del mutant, and compared these findings to the wild-type protein, and the HvHKT1;5 L189 variant and its M372del mutant. We find that α-helical association and assembly pathways in HvHKT1;5 proteins fall in two contrasting categories. Inspections of structural flexibility through molecular dynamics simulations revealed that the conformational states of HvHKT1;5 P189 diverge from those of the L189 variant and M372del mutants. We propose that M372del in HvHKT1;5 P189 instigates structural rearrangements allowing routing to the plasma membrane, while the restoration of conductance would require further interventions. We integrate the microscopy, electrophysiology, and biocomputational data and discuss how a profound structural change in HvHKT1;5 P189 M372del impacts its α-helical protein association pathway and flexibility, and how these features underlie a delicate balance leading to restoring plasma membrane localisation but not Na+ conductance.
Collapse
|
27
|
Vera AM, Galera-Prat A, Wojciechowski M, Różycki B, Laurents DV, Carrión-Vázquez M, Cieplak M, Tinnefeld P. Cohesin-dockerin code in cellulosomal dual binding modes and its allosteric regulation by proline isomerization. Structure 2021; 29:587-597.e8. [PMID: 33561387 DOI: 10.1016/j.str.2021.01.006] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/24/2020] [Revised: 11/25/2020] [Accepted: 01/11/2021] [Indexed: 12/20/2022]
Abstract
Cellulose is the most abundant organic molecule on Earth and represents a renewable and practically everlasting feedstock for the production of biofuels and chemicals. Self-assembled owing to the high-affinity cohesin-dockerin interaction, cellulosomes are huge multi-enzyme complexes with unmatched efficiency in the degradation of recalcitrant lignocellulosic substrates. The recruitment of diverse dockerin-borne enzymes into a multicohesin protein scaffold dictates the three-dimensional layout of the complex, and interestingly two alternative binding modes have been proposed. Using single-molecule fluorescence resonance energy transfer and molecular simulations on a range of cohesin-dockerin pairs, we directly detect varying distributions between these binding modes that follow a built-in cohesin-dockerin code. Surprisingly, we uncover a prolyl isomerase-modulated allosteric control mechanism, mediated by the isomerization state of a single proline residue, which regulates the distribution and kinetics of binding modes. Overall, our data provide a novel mechanistic understanding of the structural plasticity and dynamics of cellulosomes.
Collapse
Affiliation(s)
- Andrés Manuel Vera
- Department of Chemistry and Center for NanoScience, Ludwig-Maximilians-Universität München, Butenandtstr. 5-13 Haus E, 81377 München, Germany.
| | - Albert Galera-Prat
- Biocenter Oulu and Faculty of Biochemistry and Molecular Medicine, University of Oulu, 90014 Oulu, Finland
| | - Michał Wojciechowski
- Institute of Physics, Polish Academy of Sciences, Al. Lotników, 32/46, 02-668 Warsaw, Poland
| | - Bartosz Różycki
- Institute of Physics, Polish Academy of Sciences, Al. Lotników, 32/46, 02-668 Warsaw, Poland
| | - Douglas V Laurents
- Instituto de Química Física "Rocasolano", CSIC, C/ Serrano 119, 28006 Madrid, Spain
| | | | - Marek Cieplak
- Institute of Physics, Polish Academy of Sciences, Al. Lotników, 32/46, 02-668 Warsaw, Poland
| | - Philip Tinnefeld
- Department of Chemistry and Center for NanoScience, Ludwig-Maximilians-Universität München, Butenandtstr. 5-13 Haus E, 81377 München, Germany
| |
Collapse
|
28
|
Greisch JF, den Boer MA, Beurskens F, Schuurman J, Tamara S, Bondt A, Heck AJR. Generating Informative Sequence Tags from Antigen-Binding Regions of Heavily Glycosylated IgA1 Antibodies by Native Top-Down Electron Capture Dissociation. JOURNAL OF THE AMERICAN SOCIETY FOR MASS SPECTROMETRY 2021; 32:1326-1335. [PMID: 33570406 PMCID: PMC8176452 DOI: 10.1021/jasms.0c00461] [Citation(s) in RCA: 13] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/15/2023]
Abstract
Immunoglobulins A (IgA) include some of the most abundant human antibodies and play an important role in defending mucosal surfaces against pathogens. The unique structural features of the heavy chain of IgA subclasses (called IgA1 and IgA2) enable them to polymerize via the joining J-chain, resulting in IgA dimers but also higher oligomers. While secretory sIgA oligomers are dominant in milk and saliva, IgAs exist primarily as monomers in serum. No method currently allows disentangling the millions of unique IgAs potentially present in the human antibody repertoire. Obtaining unambiguous sequence reads of their hypervariable antigen-binding regions is a prerequisite for IgA identification. We here report a mass spectrometric method that uses electron capture dissociation (ECD) to produce straightforward-to-read sequence ladders of the variable parts of both the light and heavy chains of IgA1s, in particular, of the functionally critical CDR3 regions. We directly compare the native top-down ECD spectra of a heavily and heterogeneously N- and O-glycosylated anti-CD20 IgA1, the corresponding N-glycosylated anti-CD20 IgG1, and their Fab parts. We show that while featuring very different MS1 spectra, the native top-down ECD MS2 spectra of all four species are nearly identical, with cleavages occurring specifically within the CDR3 and FR4 regions of both the heavy and light chain. From the sequence-informative ECD data of an intact glycosylated IgA1, we foresee that native top-down ECD will become a valuable complementary tool for the de novo sequencing of IgA1s from milk, saliva, or serum.
Collapse
Affiliation(s)
- Jean-Francois Greisch
- Biomolecular
Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular
Research and Utrecht Institute of Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584
CH Utrecht, The Netherlands
- Netherlands
Proteomics Center, Padualaan
8, 3584 CH Utrecht, The Netherlands
| | - Maurits A. den Boer
- Biomolecular
Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular
Research and Utrecht Institute of Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584
CH Utrecht, The Netherlands
- Netherlands
Proteomics Center, Padualaan
8, 3584 CH Utrecht, The Netherlands
| | - Frank Beurskens
- Genmab,
Utrecht, Uppsalalaan
15, 3584 CT Utrecht, The Netherlands
| | - Janine Schuurman
- Genmab,
Utrecht, Uppsalalaan
15, 3584 CT Utrecht, The Netherlands
| | - Sem Tamara
- Biomolecular
Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular
Research and Utrecht Institute of Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584
CH Utrecht, The Netherlands
- Netherlands
Proteomics Center, Padualaan
8, 3584 CH Utrecht, The Netherlands
| | - Albert Bondt
- Biomolecular
Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular
Research and Utrecht Institute of Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584
CH Utrecht, The Netherlands
- Netherlands
Proteomics Center, Padualaan
8, 3584 CH Utrecht, The Netherlands
| | - Albert J. R. Heck
- Biomolecular
Mass Spectrometry and Proteomics, Bijvoet Center for Biomolecular
Research and Utrecht Institute of Pharmaceutical Sciences, Utrecht University, Padualaan 8, 3584
CH Utrecht, The Netherlands
- Netherlands
Proteomics Center, Padualaan
8, 3584 CH Utrecht, The Netherlands
| |
Collapse
|
29
|
Coban MA, Morrison J, Maharjan S, Hernandez Medina DH, Li W, Zhang YS, Freeman WD, Radisky ES, Le Roch KG, Weisend CM, Ebihara H, Caulfield TR. Attacking COVID-19 Progression Using Multi-Drug Therapy for Synergetic Target Engagement. Biomolecules 2021; 11:biom11060787. [PMID: 34071060 PMCID: PMC8224684 DOI: 10.3390/biom11060787] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/15/2021] [Revised: 05/19/2021] [Accepted: 05/20/2021] [Indexed: 12/23/2022] Open
Abstract
COVID-19 is a devastating respiratory and inflammatory illness caused by a new coronavirus that is rapidly spreading throughout the human population. Over the past 12 months, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the virus responsible for COVID-19, has already infected over 160 million (>20% located in United States) and killed more than 3.3 million people around the world (>20% deaths in USA). As we face one of the most challenging times in our recent history, there is an urgent need to identify drug candidates that can attack SARS-CoV-2 on multiple fronts. We have therefore initiated a computational dynamics drug pipeline using molecular modeling, structure simulation, docking and machine learning models to predict the inhibitory activity of several million compounds against two essential SARS-CoV-2 viral proteins and their host protein interactors-S/Ace2, Tmprss2, Cathepsins L and K, and Mpro-to prevent binding, membrane fusion and replication of the virus, respectively. All together, we generated an ensemble of structural conformations that increase high-quality docking outcomes to screen over >6 million compounds including all FDA-approved drugs, drugs under clinical trial (>3000) and an additional >30 million selected chemotypes from fragment libraries. Our results yielded an initial set of 350 high-value compounds from both new and FDA-approved compounds that can now be tested experimentally in appropriate biological model systems. We anticipate that our results will initiate screening campaigns and accelerate the discovery of COVID-19 treatments.
Collapse
Affiliation(s)
- Mathew A. Coban
- Department of Cancer Biology, Mayo Clinic, 4500 San Pablo Road South, Jacksonville, FL 32224, USA; (M.A.C.); (E.S.R.)
| | - Juliet Morrison
- Department of Microbiology and Plant Pathology, University of California, 900 University, Riverside, CA 92521, USA;
| | - Sushila Maharjan
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, 65 Landsdowne St, Cambridge, MA 02139, USA; (S.M.); (D.H.H.M.); (W.L.); (Y.S.Z.)
| | - David Hyram Hernandez Medina
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, 65 Landsdowne St, Cambridge, MA 02139, USA; (S.M.); (D.H.H.M.); (W.L.); (Y.S.Z.)
| | - Wanlu Li
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, 65 Landsdowne St, Cambridge, MA 02139, USA; (S.M.); (D.H.H.M.); (W.L.); (Y.S.Z.)
| | - Yu Shrike Zhang
- Division of Engineering in Medicine, Department of Medicine, Brigham and Women’s Hospital, Harvard Medical School, 65 Landsdowne St, Cambridge, MA 02139, USA; (S.M.); (D.H.H.M.); (W.L.); (Y.S.Z.)
| | - William D. Freeman
- Department of Neurology, Mayo Clinic, 4500 San Pablo South, Jacksonville, FL 32224, USA;
| | - Evette S. Radisky
- Department of Cancer Biology, Mayo Clinic, 4500 San Pablo Road South, Jacksonville, FL 32224, USA; (M.A.C.); (E.S.R.)
| | - Karine G. Le Roch
- Department of Molecular, Cell and Systems Biology, University of California, 900 University, Riverside, CA 92521, USA;
| | - Carla M. Weisend
- Department of Molecular Medicine, Mayo Clinic, Rochester, MN 55905, USA; (C.M.W.); (H.E.)
| | - Hideki Ebihara
- Department of Molecular Medicine, Mayo Clinic, Rochester, MN 55905, USA; (C.M.W.); (H.E.)
| | - Thomas R. Caulfield
- Department of Cancer Biology, Mayo Clinic, 4500 San Pablo Road South, Jacksonville, FL 32224, USA; (M.A.C.); (E.S.R.)
- Department of Neuroscience, Mayo Clinic, Jacksonville, FL 32224, USA
- Department of Quantitative Health Science, Division of Computational Biology, Mayo Clinic, Jacksonville, FL 32224, USA
- Department of Clinical Genomics, Mayo Clinic, Rochester, MN 55905, USA
- Department of Biochemistry & Molecular Biology, Mayo Clinic, Rochester, MN 55905, USA
- Correspondence: ; Tel.: +1-904-953-6072
| |
Collapse
|
30
|
Variants in MHY7 Gene Cause Arrhythmogenic Cardiomyopathy. Genes (Basel) 2021; 12:genes12060793. [PMID: 34067482 PMCID: PMC8224781 DOI: 10.3390/genes12060793] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2021] [Revised: 05/14/2021] [Accepted: 05/20/2021] [Indexed: 02/08/2023] Open
Abstract
BACKGROUND Arrhythmogenic Cardiomyopathy (ACM) is a disease of the cardiac muscle, characterized by frequent ventricular arrhythmias and functional/ structural abnormalities, mainly of the right ventricle. To date, 20 different genes have been associated with ACM and the majority of them encode for desmosomal proteins. In this study, we describe the characterization of two novel variants in MHY7 gene, segregating in two ACM families. MYH7 encodes for myosin heavy chain β (MHC-β) isoform, involved in cardiac muscle contractility. METHOD AND RESULTS In family A, the autopsy revealed ACM with biventricular involvement in both the proband and his father. In family B, the proband had been diagnosed as affected by ACM and implanted with implantable cardioverter defibrillator (ICD), due to ECG evidence of monomorphic ventricular tachycardia after syncope. After clinical evaluation, a molecular diagnosis was performed using a NGS custom panel. The two novel variants identified predicted damaging, located in a highly conserved domain: c. 2630T>C is not described while c.2609G>A has a frequency of 0.00000398. In silico analyses evaluated the docking characteristics between proteins using the Haddock2.2 webserver. CONCLUSIONS Our results reveal two variants in sarcomeric genes to be the molecular cause of ACM, further increasing the genetic heterogeneity of the disease; in fact, sarcomeric variants are usually associated with HCM phenotype. Studies on the role of sarcomere genes in the pathogenesis of ACM are surely recommended in those ACM patients negative for desmosomal mutation screening.
Collapse
|
31
|
Li RX, Zhang NN, Wu B, OuYang B, Shen HB. Multiobjective heuristic algorithm for de novo protein design in a quantified continuous sequence space. Comput Struct Biotechnol J 2021; 19:2575-2587. [PMID: 34025944 PMCID: PMC8114120 DOI: 10.1016/j.csbj.2021.04.046] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Revised: 04/19/2021] [Accepted: 04/22/2021] [Indexed: 11/12/2022] Open
Abstract
Protein design usually involves sequence search process and evaluation criteria. Commonly used methods primarily implement the Monte Carlo or simulated annealing algorithm with a single-energy function to obtain ideal solutions, which is often highly time-consuming and limited by the accuracy of the energy function. In this report, we introduce a multiobjective algorithm named Hydra for protein design, which employs two different energy functions to optimize solutions simultaneously and makes use of the latent quantitative relationship between different amino acid types to facilitate the search process. The framework uses two kinds of prior information to transform the original disordered discrete sequence space into a relatively ordered space, and decoy sequences are searched in this ordered space through a multiobjective swarm intelligence algorithm. This algorithm features high accuracy and a high-speed search process. Our method was tested on 40 targets covering different fold classes, which were computationally verified to be well folded, and it experimentally solved the 1UBQ fold by NMR in excellent agreement with the native structure with a backbone RMSD deviation of 1.074 Å. The Hydra software package can be downloaded from: http://www.csbio.sjtu.edu.cn/bioinf/HYDRA/ for academic use.
Collapse
Affiliation(s)
- Rui-Xiang Li
- Institute of Image Processing and Pattern Recognition, Shanghai Jiao Tong University, and Key Laboratory of System Control and Information Processing, Ministry of Education of China, Shanghai 200240, China
| | - Ning-Ning Zhang
- State Key Laboratory of Molecular Biology, CAS Center for Excellence in Molecular Cell Science, Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, University of Chinese Academy of Sciences, Shanghai 201203, China
| | - Bin Wu
- National Facility for Protein Science in Shanghai, ZhangJiang Lab, Shanghai Advanced Research Institute, Chinese Academy of Sciences, Shanghai 201210, China
| | - Bo OuYang
- State Key Laboratory of Molecular Biology, CAS Center for Excellence in Molecular Cell Science, Shanghai Institute of Biochemistry and Cell Biology, Chinese Academy of Sciences, University of Chinese Academy of Sciences, Shanghai 201203, China
| | - Hong-Bin Shen
- Institute of Image Processing and Pattern Recognition, Shanghai Jiao Tong University, and Key Laboratory of System Control and Information Processing, Ministry of Education of China, Shanghai 200240, China.,Department of Computer Science, Shanghai Jiao Tong University, Shanghai 200240, China
| |
Collapse
|
32
|
Garg A, Pal D. Inferring metal binding sites in flexible regions of proteins. Proteins 2021; 89:1125-1133. [PMID: 33864411 DOI: 10.1002/prot.26085] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2021] [Revised: 03/15/2021] [Accepted: 04/12/2021] [Indexed: 12/20/2022]
Abstract
Metal ions are central to the molecular function of many proteins. Thus their knowledge in experimentally determined structure is important; however, such structures often lose bound metal ions during sample preparation. Identification of these metal-binding site(s) becomes difficult when the receptor is novel and/or their conformations differ in the bound/unbound states. Locating such sites in theoretical models also poses a challenge due to the uncertainties with side-chain modeling. We address the problem by employing the Geometric Hashing algorithm to create a template library of functionally important binding sites and match query structures with the available templates. The matching is done on the structure ensemble obtained from coarse-grained molecular dynamics simulation, where metal-specific amino acids are screened to infer the true site. Test on 1347 non-redundant monomer protein structures show that Ca2+ , Zn2+ , Mg2+ , Cu2+ , and Fe3+ binding site residues can be classified at 0.92, 0.95, 0.80, 0.90, and 0.92 aggregate performance (out of 1) across all possible thresholds. The performance for Ca2+ and Zn2+ is notably superior in comparison to state-of-the-art methods like IonCom and MIB. Specific case studies show that additionally predicted metal-binding site residues in proteins have features necessary for ion binding. These include new sites not predicted by other methods. The use of coarse-grained dynamics thus provides a generalized approach to improve metal-binding site prediction. The work is expected to contribute to improving our ability to correctly predict protein molecular function where knowledge of metal binding is a key requirement.
Collapse
Affiliation(s)
- Aditi Garg
- Department of Computational and Data Sciences, Indian Institute of Science, Bengaluru, India
| | - Debnath Pal
- Department of Computational and Data Sciences, Indian Institute of Science, Bengaluru, India
| |
Collapse
|
33
|
Patel D, Patel JS, Ytreberg FM. Implementing and Assessing an Alchemical Method for Calculating Protein-Protein Binding Free Energy. J Chem Theory Comput 2021; 17:2457-2464. [PMID: 33709712 PMCID: PMC8044032 DOI: 10.1021/acs.jctc.0c01045] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/18/2023]
Abstract
Protein-protein binding is fundamental to most biological processes. It is important to be able to use computation to accurately estimate the change in protein-protein binding free energy due to mutations in order to answer biological questions that would be experimentally challenging, laborious, or time-consuming. Although nonrigorous free-energy methods are faster, rigorous alchemical molecular dynamics-based methods are considerably more accurate and are becoming more feasible with the advancement of computer hardware and molecular simulation software. Even with sufficient computational resources, there are still major challenges to using alchemical free-energy methods for protein-protein complexes, such as generating hybrid structures and topologies, maintaining a neutral net charge of the system when there is a charge-changing mutation, and setting up the simulation. In the current study, we have used the pmx package to generate hybrid structures and topologies, and a double-system/single-box approach to maintain the net charge of the system. To test the approach, we predicted relative binding affinities for two protein-protein complexes using a nonequilibrium alchemical method based on the Crooks fluctuation theorem and compared the results with experimental values. The method correctly identified stabilizing from destabilizing mutations for a small protein-protein complex, and a larger, more challenging antibody complex. Strong correlations were obtained between predicted and experimental relative binding affinities for both protein-protein systems.
Collapse
Affiliation(s)
- Dharmeshkumar Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
| | - Jagdish Suresh Patel
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
- Department of Biological Sciences, University of Idaho, Moscow, Idaho 83844, United States
| | - F Marty Ytreberg
- Institute for Modeling Collaboration and Innovation, University of Idaho, Moscow, Idaho 83844, United States
- Department of Physics, University of Idaho, Moscow, Idaho 83844, United States
| |
Collapse
|
34
|
Pereira GRC, Vieira BDAA, De Mesquita JF. Comprehensive in silico analysis and molecular dynamics of the superoxide dismutase 1 (SOD1) variants related to amyotrophic lateral sclerosis. PLoS One 2021; 16:e0247841. [PMID: 33630959 PMCID: PMC7906464 DOI: 10.1371/journal.pone.0247841] [Citation(s) in RCA: 19] [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: 10/02/2020] [Accepted: 02/15/2021] [Indexed: 12/29/2022] Open
Abstract
Amyotrophic Lateral Sclerosis (ALS) is the most frequent motor neuron disorder, with a significant social and economic burden. ALS remains incurable, and the only drugs approved for its treatments confers a survival benefit of a few months for the patients. Missense mutations in superoxide dismutase 1 (SOD1), a major cytoplasmic antioxidant enzyme, has been associated with ALS development, accounting for 23% of its familial cases and 7% of all sporadic cases. This work aims to characterize in silico the structural and functional effects of SOD1 protein variants. Missense mutations in SOD1 were compiled from the literature and databases. Twelve algorithms were used to predict the functional and stability effects of these mutations. ConSurf was used to estimate the evolutionary conservation of SOD1 amino-acids. GROMACS was used to perform molecular dynamics (MD) simulations of SOD1 wild-type and variants A4V, D90A, H46R, and I113T, which account for approximately half of all ALS-SOD1 cases in the United States, Europe, Japan, and United Kingdom, respectively. 233 missense mutations in SOD1 protein were compiled from the databases and literature consulted. The predictive analyses pointed to an elevated rate of deleterious and destabilizing predictions for the analyzed variants, indicating their harmful effects. The ConSurf analysis suggested that mutations in SOD1 mainly affect conserved and possibly functionally essential amino acids. The MD analyses pointed to flexibility and essential dynamics alterations at the electrostatic and metal-binding loops of variants A4V, D90A, H46R, and I113T that could lead to aberrant interactions triggering toxic protein aggregation. These alterations may have harmful implications for SOD1 and explain their association with ALS. Understanding the effects of SOD1 mutations on protein structure and function facilitates the design of further experiments and provides relevant information on the molecular mechanism of pathology, which may contribute to improvements in existing treatments for ALS.
Collapse
Affiliation(s)
- Gabriel Rodrigues Coutinho Pereira
- Department of Genetics and Molecular Biology, Bioinformatics and Computational Biology Laboratory, Federal University of the State of Rio de Janeiro (UNIRIO), Rio de Janeiro, Rio de Janeiro, Brazil
| | | | - Joelma Freire De Mesquita
- Department of Genetics and Molecular Biology, Bioinformatics and Computational Biology Laboratory, Federal University of the State of Rio de Janeiro (UNIRIO), Rio de Janeiro, Rio de Janeiro, Brazil
- * E-mail:
| |
Collapse
|
35
|
Norris GA, Tsai ACH, Schneider KW, Wu YH, Caulfield T, Green AL. A novel, germline, deactivating CBL variant p.L493F alters domain orientation and is associated with multiple childhood cancers. Cancer Genet 2021; 254-255:18-24. [PMID: 33550024 DOI: 10.1016/j.cancergen.2021.01.008] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2020] [Revised: 11/26/2020] [Accepted: 01/18/2021] [Indexed: 10/22/2022]
Abstract
CBL is a mammalian gene encoding the protein CBL, which is an E3 ubiquitin-protein ligase involved in cell signaling and protein ubiquitination. Pathogenic variants in this gene have been implicated in a number of human cancers, particularly acute myeloid leukemia (AML). Here, we present a 5-year-old male patient with a history of AML, diffuse midline glioma, and left brain lesion with histiocytic features. A variant of uncertain significance (VUS): p.L493F was detected in his CBL gene via clinical evaluation. Protein modeling predicts this variant to be pathogenic. Details of the clinical evaluation and modeling assay are discussed.
Collapse
Affiliation(s)
- Gregory A Norris
- Department of Pediatrics, University of Colorado Anschutz Medical Campus, Children's Hospital Colorado, Aurora, CO USA
| | - Anne Chun-Hui Tsai
- Department of Pediatrics, University of Colorado Anschutz Medical Campus, Children's Hospital Colorado, Aurora, CO USA
| | - Kami Wolfe Schneider
- Department of Pediatrics, University of Colorado Anschutz Medical Campus, Children's Hospital Colorado, Aurora, CO USA
| | - Yuan-Haw Wu
- Department of Pediatrics, University of Colorado Anschutz Medical Campus, Children's Hospital Colorado, Aurora, CO USA
| | - Thomas Caulfield
- Department of Neuroscience, Mayo Clinic, Jacksonville, FL, 32224, USA; Mayo Graduate School, Neurobiology of Disease, Mayo Clinic, Jacksonville, FL USA
| | - Adam L Green
- Department of Pediatrics, University of Colorado Anschutz Medical Campus, Children's Hospital Colorado, Aurora, CO USA; Morgan Adams Foundation Pediatric Brain Tumor Research Program, Aurora, CO USA.
| |
Collapse
|
36
|
Richter JE, Hines S, Selvam P, Atwal H, Farres H, Caulfield TR, Atwal PS. Clinical description & molecular modeling of novel MAX pathogenic variant causing pheochromocytoma in family, supports paternal parent-of-origin effect. Cancer Genet 2021; 252-253:107-110. [PMID: 33493868 DOI: 10.1016/j.cancergen.2021.01.004] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2020] [Revised: 12/01/2020] [Accepted: 01/15/2021] [Indexed: 10/22/2022]
Abstract
The titular member of the MAX network of proteins, MYC-associated factor X (MAX), serves an important regulatory function in transcription of E-box genes associated with cell proliferation, differentiation, and apoptosis. Wild type MAX dimerizes with both MYC and MAD, both of which are members of the MAX network, and can promote or repress cell functions as needed. However, pathogenic variants in MAX are known to upset this balance, leading to uncontrolled oncogenic activity and disease phenotypes such as paragangliomas and pheochromocytomas. We report a 58-year-old male and his 32-year-old daughter, both of which have a history of pheochromocytoma and the unique nonsense MAX variant c.271C>T (p.Q91X). These individuals were diagnosed with pheochromocytomas in their early twenties that were later removed through corrective surgery. The father now presents with recurrent symptoms of hypertension, hyperhidrosis, and headaches, which accompany new pheochromocytomas of his remaining adrenal gland. Pathogenicity of this MAX variant is proven through molecular modeling. The case of this father-daughter pair supports both heritability of pheochromocytoma and the paternal parent-of-origin effect for MAX pathogenic variants.
Collapse
Affiliation(s)
- John E Richter
- Department of Clinical Genomics, Mayo Clinic, Jacksonville, FL 32224, USA
| | - S Hines
- Department of Medicine, Mayo Clinic, Jacksonville, FL, USA
| | - Pavalan Selvam
- Atwal Clinic: Genomic & Personalized Medicine, Palm Beach, FL 33480, USA
| | - Herjot Atwal
- Atwal Clinic: Genomic & Personalized Medicine, Palm Beach, FL 33480, USA
| | - Houssam Farres
- Department of Surgery, Mayo Clinic, Jacksonville, FL, USA
| | - Thomas R Caulfield
- Department of Neuroscience, Mayo Clinic, Jacksonville, FL 32224, USA; Mayo Graduate School, Neurobiology of Disease, Mayo Clinic, Jacksonville, FL, USA
| | - Paldeep S Atwal
- Atwal Clinic: Genomic & Personalized Medicine, Palm Beach, FL 33480, USA.
| |
Collapse
|
37
|
Abstract
![]()
We
present a novel web server, named gridSolvate, dedicated to
the prediction of biomolecular hydration properties. Given a solute
in atomic representation, such as a protein or protein–ligand
complex, the server determines positions and excess chemical potential
of buried and first hydration shell water molecules. Calculations
are based on our semiexplicit hydration model that provides computational
efficiency close to implicit solvent approaches, yet captures a number
of physical effects unique to explicit solvent representation. The
model was introduced and validated before in the context of bulk hydration
of drug-like solutes and determination of protein hydration sites.
Current methodological developments merge those two avenues into a
single, easily accessible tool. Here, we focus on the server’s
ability to predict water distribution and affinity within protein–ligand
interfaces. We demonstrate that with possibly minimal user intervention
the server correctly predicts the locations of 77% of interface water
molecules in an external set of test structures. The server is freely
available at https://gsolvate.biomod.cent.uw.edu.pl.
Collapse
Affiliation(s)
- Piotr Setny
- Centre of New Technologies, University of Warsaw, Banacha 2c, 02-097 Warsaw, Poland
| |
Collapse
|
38
|
Zheng Y, Wang J, Liang X, Huang H, Ma Y, Lin L, Wang C, Zhan X, Yang L, Zha G, Yang P, Zou X, Chen Z, Chen X, Chen W, Liu X, Lin M. Epidemiology, evolutionary origin, and malaria-induced positive selection effects of G6PD-deficient alleles in Chinese populations. Mol Genet Genomic Med 2020; 8:e1540. [PMID: 33128437 PMCID: PMC7767544 DOI: 10.1002/mgg3.1540] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/04/2020] [Revised: 09/02/2020] [Accepted: 09/17/2020] [Indexed: 02/05/2023] Open
Abstract
BACKGROUND Although glucose-6-phosphate dehydrogenase (G6PD) deficiency is the most common inherited disorder in the Chinese population, there is scarce evidence regarding the epidemiology, evolutionary origin, and malaria-induced positive selection effects of G6PD-deficient alleles in various Chinese ethnic populations. METHODS We performed a large population-based screening (n = 15,690) to examine the impact of selection on human nucleotide diversity and to infer the evolutionary history of the most common deficiency alleles in Chinese populations. RESULTS The frequencies of G6PD deficiency ranged from 0% to 11.6% in 12 Chinese ethnic populations. A frequency map based on geographic information showed that G6PD deficiency was highly correlated with historical malaria prevalence in China and was affected by altitude and latitude. The five most frequently occurring G6PD gene variants were NM_001042351.3:c.1376G>T, NM_001042351.3:c.1388G>A, NM_001042351.3:c.95A>G, NM_001042351.3:c.1311T>C, and NM_001042351.3:c.1024C>T, which were distributed with ethnic features. A pathogenic but rarely reported variant site (NM_001042351.3:c.448G>A) was identified in this study. Bioinformatic analysis revealed a strong and recent positive selection targeting the NM_001042351.3:c.1376G>T allele that originated in the past 3125 to 3750 years and another selection targeting the NM_001042351.3:c.1388G>A allele that originated in the past 5000 to 6000 years. Additionally, both alleles originated from a single ancestor. CONCLUSION These results indicate that malaria has had a major impact on the Chinese genome since the introduction of rice agriculture.
Collapse
Affiliation(s)
- Yuzhong Zheng
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Junli Wang
- Reproductive Medicine CenterThe Affiliated Hospital of Youjiang Medical University for NationalitiesBaiseChina
| | - Xueyan Liang
- Department of Medical GeneticsShantou University Medical CollegeShantouGuangdongChina
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| | - Huiying Huang
- Department of Medical GeneticsShantou University Medical CollegeShantouGuangdongChina
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| | - Yanbo Ma
- School of Mathematics and StatisticsHanshan Normal UniversityChaozhouGuangdongChina
| | - Liyun Lin
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Chunfang Wang
- Reproductive Medicine CenterThe Affiliated Hospital of Youjiang Medical University for NationalitiesBaiseChina
| | - Xiaofen Zhan
- Department of Medical LaboratoryChaozhou Central Hospital Affiliated to Southern Medical UniversityChaozhouGuangdongChina
| | - Liye Yang
- Department of Medical LaboratoryChaozhou Central Hospital Affiliated to Southern Medical UniversityChaozhouGuangdongChina
| | - Guangcai Zha
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Peikui Yang
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Xianghui Zou
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Zikai Chen
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
| | - Xinyao Chen
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| | - Weizhong Chen
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| | - Xiangzhi Liu
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| | - Min Lin
- School of Food Engineering and BiotechnologyHanshan Normal UniversityChaozhouGuangdong ProvinceChina
- Department of Medical LaboratoryChaozhou People’s Hospital Affiliated to Shantou University Medical CollegeChaozhouGuangdongChina
| |
Collapse
|
39
|
Delgado Blanco J, Hernandez-Alias X, Cianferoni D, Serrano L. In silico mutagenesis of human ACE2 with S protein and translational efficiency explain SARS-CoV-2 infectivity in different species. PLoS Comput Biol 2020; 16:e1008450. [PMID: 33284795 PMCID: PMC7746295 DOI: 10.1371/journal.pcbi.1008450] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 12/17/2020] [Accepted: 10/19/2020] [Indexed: 01/08/2023] Open
Abstract
The coronavirus disease COVID-19 constitutes the most severe pandemic of the last decades having caused more than 1 million deaths worldwide. The SARS-CoV-2 virus recognizes the angiotensin converting enzyme 2 (ACE2) on the surface of human cells through its spike protein. It has been reported that the coronavirus can mildly infect cats, and ferrets, and perhaps dogs while not pigs, mice, chicken and ducks. Differences in viral infectivity among different species or individuals could be due to amino acid differences at key positions of the host proteins that interact with the virus, the immune response, expression levels of host proteins and translation efficiency of the viral proteins among other factors. Here, first we have addressed the importance that sequence variants of different animal species, human individuals and virus isolates have on the interaction between the RBD domain of the SARS-CoV-2 spike S protein and human angiotensin converting enzyme 2 (ACE2). Second, we have looked at viral translation efficiency by using the tRNA adaptation index. We find that integration of both interaction energy with ACE2 and translational efficiency explains animal infectivity. Humans are the top species in which SARS-CoV-2 is both efficiently translated as well as optimally interacting with ACE2. We have found some viral mutations that increase affinity for hACE and some hACE2 variants affecting ACE2 stability and virus binding. These variants suggest that different sensitivities to coronavirus infection in humans could arise in some cases from allelic variability affecting ACE2 stability and virus binding.
Collapse
Affiliation(s)
- Javier Delgado Blanco
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, Spain
| | - Xavier Hernandez-Alias
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, Spain
| | - Damiano Cianferoni
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, Spain
| | - Luis Serrano
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology, Dr. Aiguader 88, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- ICREA, Pg. Lluís Companys 23, Barcelona, Spain
| |
Collapse
|
40
|
Ahuja AS, Selvam P, Vadlamudi C, Chopra H, Richter JE, Macklin SK, Samreen A, Helmi H, Mohammaad AN, Hines S, Davila MC, Atwal PS, Caulfield TR. Genomics combined with a protein informatics platform to assess a novel pathogenic variant c.1024 A>G (p.K342E) in OPA1 in a patient with autosomal dominant optic atrophy. Ophthalmic Genet 2020; 41:563-569. [PMID: 32940104 DOI: 10.1080/13816810.2020.1814344] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 10/23/2022]
Abstract
BACKGROUND Autosomal Dominant Optic Atrophy (ADOA) is caused by mutations in the Optic Atrophy 1 Gene which disrupts the OPA1 protein. This disruption affects the normal function of the protein; impairs fusion of the mitochondrial inner membrane; and prevents normal OPA1 protein degradation. These events cause damage in retinal ganglion cells that could affect the patients with symptoms ranging from none to legally blind. MATERIALS AND METHODS Our study identifies a missense variant mutation, c.1024 A > G (p.K342E), in OPA1 gene causing ADOA. Diagnosed clinically in three family members and the presence of this mutation was confirmed in two members by genetic testing. Pathogenic variants in OPA1 impact the secondary protein structure and function by causing non-conservative amino acid substitutions. We also modeled this mutation and compared it to the wild type using statistical mechanics. RESULTS AND CONCLUSIONS The proband's pathogenic variant, c.1024 A > G (p.K342E), is located in the GTPase domain of OPA1 and causes changes in the protein structure by affecting the oligomerization pattern thus resulting in ADOA. Identifying the pathogenic potential of the missense mutations in the OPA1 gene using neoteric protein modeling techniques would help in the early detection of ADOA in patients who have family history of blindness. This action would help in providing early follow up, possible treatment in the future, and genetic counseling. Abbreviations: ADOA: Autosomal Dominant Optic Atrophy; CYCS: Caspase Activator Cytochrome C; OPA1: Optic Atrophy Gene 1; RGC: Retinal Ganglion Cells; VUS: Variant of Uncertain Significance.
Collapse
Affiliation(s)
- Abhimanyu S Ahuja
- Charles E. Schmidt College of Medicine, Florida Atlantic University , Boca Raton, Florida, USA
| | - Pavalan Selvam
- Department of Clinical Genomics, Atwal Clinic , Jacksonville, Florida, USA
| | - Charitha Vadlamudi
- Department of Clinical Genomics, Atwal Clinic , Jacksonville, Florida, USA
| | - Hayley Chopra
- Department of Molecular and Human Genetics, Baylor College of Medicine , Houston, Texas, USA
| | | | | | - Ayesha Samreen
- Department of Molecular and Human Genetics, Baylor College of Medicine , Houston, Texas, USA
| | - Haytham Helmi
- Department of Endocrinology, Mayo Clinic , Jacksonville, Florida, USA
| | - Ahmed N Mohammaad
- Department of Endocrinology, Mayo Clinic , Jacksonville, Florida, USA
| | - Stephanie Hines
- Department of Medicine, Division of Diagnostic & Consultative Medicine, Mayo Clinic , Jacksonville, Florida, USA
| | | | - Paldeep S Atwal
- Department of Clinical Genomics, Atwal Clinic , Jacksonville, Florida, USA
| | - Thomas R Caulfield
- Department of Neuroscience, Mayo Clinic , Jacksonville, Florida, USA.,Department of Neurosurgery, Mayo Clinic , Jacksonville, Florida, USA.,Department of Cancer Biology, Mayo Clinic , Jacksonville, Florida, USA.,Department of Health Sciences Research, Mayo Clinic , Jacksonville, Florida, USA
| |
Collapse
|
41
|
Borjigin C, Schilling RK, Bose J, Hrmova M, Qiu J, Wege S, Situmorang A, Byrt C, Brien C, Berger B, Gilliham M, Pearson AS, Roy SJ. A single nucleotide substitution in TaHKT1;5-D controls shoot Na + accumulation in bread wheat. PLANT, CELL & ENVIRONMENT 2020; 43:2158-2171. [PMID: 32652543 PMCID: PMC7540593 DOI: 10.1111/pce.13841] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2020] [Revised: 05/20/2020] [Accepted: 05/24/2020] [Indexed: 05/22/2023]
Abstract
Improving salinity tolerance in the most widely cultivated cereal, bread wheat (Triticum aestivum L.), is essential to increase grain yields on saline agricultural lands. A Portuguese landrace, Mocho de Espiga Branca accumulates up to sixfold greater leaf and sheath sodium (Na+ ) than two Australian cultivars, Gladius and Scout, under salt stress in hydroponics. Despite high leaf and sheath Na+ concentrations, Mocho de Espiga Branca maintained similar salinity tolerance compared to Gladius and Scout. A naturally occurring single nucleotide substitution was identified in the gene encoding a major Na+ transporter TaHKT1;5-D in Mocho de Espiga Branca, which resulted in a L190P amino acid residue variation. This variant prevents Mocho de Espiga Branca from retrieving Na+ from the root xylem leading to a high shoot Na+ concentration. The identification of the tissue-tolerant Mocho de Espiga Branca will accelerate the development of more elite salt-tolerant bread wheat cultivars.
Collapse
Affiliation(s)
- Chana Borjigin
- Australian Centre for Plant Functional Genomics, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Rhiannon K. Schilling
- Australian Centre for Plant Functional Genomics, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Jayakumar Bose
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Centre of Excellence in Plant Energy BiologyThe University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Maria Hrmova
- Australian Centre for Plant Functional Genomics, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Life Sciences, Huaiyin Normal UniversityHuai'anChina
| | - Jiaen Qiu
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Centre of Excellence in Plant Energy BiologyThe University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Stefanie Wege
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Centre of Excellence in Plant Energy BiologyThe University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Apriadi Situmorang
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Caitlin Byrt
- Division of Plant SciencesResearch School of Biology, Australian National UniversityActonAustralian Capital TerritoryAustralia
| | - Chris Brien
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- Australian Plant Phenomics FacilityThe Plant Accelerator, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Bettina Berger
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- Australian Plant Phenomics FacilityThe Plant Accelerator, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Matthew Gilliham
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Centre of Excellence in Plant Energy BiologyThe University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Allison S. Pearson
- Australian Centre for Plant Functional Genomics, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Centre of Excellence in Plant Energy BiologyThe University of AdelaideGlen OsmondSouth AustraliaAustralia
| | - Stuart J. Roy
- Australian Centre for Plant Functional Genomics, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- School of Agriculture, Food and Wine, The University of AdelaideGlen OsmondSouth AustraliaAustralia
- ARC Industrial Transformation Research Hub for Wheat in a Hot Dry Climate, The University of AdelaideGlen OsmondSouth AustraliaAustralia
| |
Collapse
|
42
|
Borjigin C, Schilling RK, Bose J, Hrmova M, Qiu J, Wege S, Situmorang A, Byrt C, Brien C, Berger B, Gilliham M, Pearson AS, Roy SJ. A single nucleotide substitution in TaHKT1;5-D controls shoot Na + accumulation in bread wheat. PLANT, CELL & ENVIRONMENT 2020; 43:2158-2171. [PMID: 32652543 DOI: 10.1101/2020.01.21.909887] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/16/2020] [Revised: 05/20/2020] [Accepted: 05/24/2020] [Indexed: 05/28/2023]
Abstract
Improving salinity tolerance in the most widely cultivated cereal, bread wheat (Triticum aestivum L.), is essential to increase grain yields on saline agricultural lands. A Portuguese landrace, Mocho de Espiga Branca accumulates up to sixfold greater leaf and sheath sodium (Na+ ) than two Australian cultivars, Gladius and Scout, under salt stress in hydroponics. Despite high leaf and sheath Na+ concentrations, Mocho de Espiga Branca maintained similar salinity tolerance compared to Gladius and Scout. A naturally occurring single nucleotide substitution was identified in the gene encoding a major Na+ transporter TaHKT1;5-D in Mocho de Espiga Branca, which resulted in a L190P amino acid residue variation. This variant prevents Mocho de Espiga Branca from retrieving Na+ from the root xylem leading to a high shoot Na+ concentration. The identification of the tissue-tolerant Mocho de Espiga Branca will accelerate the development of more elite salt-tolerant bread wheat cultivars.
Collapse
Affiliation(s)
- Chana Borjigin
- Australian Centre for Plant Functional Genomics, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Rhiannon K Schilling
- Australian Centre for Plant Functional Genomics, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Jayakumar Bose
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Centre of Excellence in Plant Energy Biology, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Maria Hrmova
- Australian Centre for Plant Functional Genomics, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Life Sciences, Huaiyin Normal University, Huai'an, China
| | - Jiaen Qiu
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Centre of Excellence in Plant Energy Biology, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Stefanie Wege
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Centre of Excellence in Plant Energy Biology, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Apriadi Situmorang
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Caitlin Byrt
- Division of Plant Sciences, Research School of Biology, Australian National University, Acton, Australian Capital Territory, Australia
| | - Chris Brien
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- Australian Plant Phenomics Facility, The Plant Accelerator, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Bettina Berger
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- Australian Plant Phenomics Facility, The Plant Accelerator, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Matthew Gilliham
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Centre of Excellence in Plant Energy Biology, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Allison S Pearson
- Australian Centre for Plant Functional Genomics, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Centre of Excellence in Plant Energy Biology, The University of Adelaide, Glen Osmond, South Australia, Australia
| | - Stuart J Roy
- Australian Centre for Plant Functional Genomics, The University of Adelaide, Glen Osmond, South Australia, Australia
- School of Agriculture, Food and Wine, The University of Adelaide, Glen Osmond, South Australia, Australia
- ARC Industrial Transformation Research Hub for Wheat in a Hot Dry Climate, The University of Adelaide, Glen Osmond, South Australia, Australia
| |
Collapse
|
43
|
Examination of Molecular Effects of MYLK Deletion in a Patient with Extensive Aortic, Carotid, and Abdominal Dissections That Underlie the Genetic Dysfunction. Case Rep Med 2020; 2020:5108052. [PMID: 32655646 PMCID: PMC7327611 DOI: 10.1155/2020/5108052] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2019] [Revised: 04/12/2020] [Accepted: 05/27/2020] [Indexed: 11/17/2022] Open
Abstract
We describe the phenotype of a patient with extensive aortic, carotid, and abdominal dissections. The proband was found to have a heterozygous deletion of exons 21-34 in MYLK, which is a rare finding, as deletions in this gene have been infrequently reported. We describe this finding following detection in a proband with an extensive history of aortic, carotid, and abdominal dissections. Neoteric molecular modeling techniques to help determine the impact of this deletion on protein function indicated loss of function due to lack of any kinase domain. We also provide the electrostatics calculations from the wild type and mutant variant. Through a combined multiomic approach of clinical, functional, and protein informatics, we arrive at a data fusion for determination of pathogenicity embedded within the genetic code for this particular genetic variant, which, as a platform, continues to broaden its scope across the field of variants of uncertain significance classification.
Collapse
|
44
|
Bhat AS, Dustin Schaeffer R, Kinch L, Medvedev KE, Grishin NV. Recent advances suggest increased influence of selective pressure in allostery. Curr Opin Struct Biol 2020; 62:183-188. [PMID: 32302874 DOI: 10.1016/j.sbi.2020.02.004] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/01/2019] [Revised: 02/21/2020] [Accepted: 02/25/2020] [Indexed: 12/19/2022]
Abstract
Allosteric regulation of protein functions is ubiquitous in organismal biology, but the principles governing its evolution are not well understood. Here we discuss recent studies supporting the large-scale existence of latent allostery in ancestor proteins of superfamilies. As suggested, the evolution of allostery could be driven by the need for specificity in paralogs of slow evolving protein complexes with conserved active sites. The same slow evolution is displayed by purifying selection exhibited in allosteric proteins with somatic mutations involved in cancer, where disease-associated mutations are enriched in both orthosteric and allosteric sites. Consequently, disease-associated variants can be used to identify druggable allosteric sites that are specific for paralogs in protein superfamilies with otherwise similar functions.
Collapse
Affiliation(s)
- Archana S Bhat
- Department of Biophysics, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States
| | - Richard Dustin Schaeffer
- Department of Biophysics, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States
| | - Lisa Kinch
- Howard Hughes Medical Institute, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States
| | - Kirill E Medvedev
- Department of Biophysics, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States
| | - Nick V Grishin
- Department of Biophysics, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States; Howard Hughes Medical Institute, University of Texas Southwestern Medical Center, Dallas, TX 75390-9050, United States.
| |
Collapse
|
45
|
The Ca 2+ permeation mechanism of the ryanodine receptor revealed by a multi-site ion model. Nat Commun 2020; 11:922. [PMID: 32066742 PMCID: PMC7026163 DOI: 10.1038/s41467-020-14573-w] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2019] [Accepted: 01/16/2020] [Indexed: 12/22/2022] Open
Abstract
Ryanodine receptors (RyR) are ion channels responsible for the release of Ca2+ from the sarco/endoplasmic reticulum and play a crucial role in the precise control of Ca2+ concentration in the cytosol. The detailed permeation mechanism of Ca2+ through RyR is still elusive. By using molecular dynamics simulations with a specially designed Ca2+ model, we show that multiple Ca2+ ions accumulate in the upper selectivity filter of RyR1, but only one Ca2+ can occupy and translocate in the narrow pore at a time, assisted by electrostatic repulsion from the Ca2+ within the upper selectivity filter. The Ca2+ is nearly fully hydrated with the first solvation shell intact during the whole permeation process. These results suggest a remote knock-on permeation mechanism and one-at-a-time occupation pattern for the hydrated Ca2+ within the narrow pore, uncovering the basis underlying the high permeability and low selectivity of the RyR channels.
Collapse
|
46
|
Characterization of a Pathogenic Variant in the ABCD1 Gene Through Protein Molecular Modeling. Case Rep Genet 2020; 2020:3256539. [PMID: 32047678 PMCID: PMC7003284 DOI: 10.1155/2020/3256539] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/15/2019] [Accepted: 11/19/2019] [Indexed: 11/18/2022] Open
Abstract
Background The ATP-binding cassette, subfamily D, member 1 (ABCD1) protein is a peroxisomal half-transporter that allows for very long chain fatty acid (VLCFA) degradation. Pathogenic variants of ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death. Methods A newly characterized and suspected pathogenic variant in ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death. Results A case of adult onset adrenomyeloneuropathy (AMN) and a novel ABCD1 cause VLCFAs to build up in various tissues and bodily fluids, resulting in a disorder called X-linked adrenoleukodystrophy (X-ALD). This disorder is most commonly marked by adrenocortical insufficiency and high VLCFA concentration, and has varying levels of neurological involvement depending on phenotype. For example, the Addison-only form of X-ALD has no neurological impact, while the cerebral form of X-ALD often causes severe sensory loss, motor function impairment, cognitive decline, and death. Conclusions Data fusion from multiple sources was combined in a comprehensive approach yielding an enriched assessment of the patient's disease and prognosis. Molecular modeling was performed on the variant to better characterize its clinical significance and confirm pathogenicity.
Collapse
|
47
|
Martins IC, Santos NC. Intrinsically disordered protein domains in flavivirus infection. Arch Biochem Biophys 2020; 683:108298. [PMID: 32045581 DOI: 10.1016/j.abb.2020.108298] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2019] [Revised: 02/03/2020] [Accepted: 02/05/2020] [Indexed: 12/30/2022]
Abstract
Intrinsically disordered protein regions are at the core of biological processes and involved in key protein-ligand interactions. The Flavivirus proteins, of viruses of great biomedical importance such as Zika and dengue viruses, exemplify this. Several proteins of these viruses have disordered regions that are of the utmost importance for biological activity. Disordered proteins can adopt several conformations, each able to interact with and/or bind to different ligands. In fact, such interactions can help stabilize a particular fold. Moreover, by being promiscuous in the number of target molecules they can bind to, these protein regions increase the number of functions that their small proteome (10 proteins) can achieve. A folding energy waterfall better describes the protein folding landscape of these proteins. A disordered protein can be thought as rolling down the folding energy cascade, in order "to fall, fold and function". This is the case of many viral protein regions, as seen in the flaviviruses proteome. Given their small size, flaviviruses are a good model system for understanding the role of intrinsically disordered protein regions in viral function. Finally, studying these viruses disordered protein regions will certainly contribute to the development of therapeutic approaches against such promising (yet challenging) targets.
Collapse
Affiliation(s)
- Ivo C Martins
- Instituto de Medicina Molecular, Faculdade de Medicina, Universidade de Lisboa, Lisbon, Portugal.
| | - Nuno C Santos
- Instituto de Medicina Molecular, Faculdade de Medicina, Universidade de Lisboa, Lisbon, Portugal.
| |
Collapse
|
48
|
Pal LR, Kundu K, Yin Y, Moult J. Matching whole genomes to rare genetic disorders: Identification of potential causative variants using phenotype-weighted knowledge in the CAGI SickKids5 clinical genomes challenge. Hum Mutat 2020; 41:347-362. [PMID: 31680375 PMCID: PMC7182498 DOI: 10.1002/humu.23933] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/09/2019] [Revised: 09/26/2019] [Accepted: 10/13/2019] [Indexed: 02/06/2023]
Abstract
Precise identification of causative variants from whole-genome sequencing data, including both coding and noncoding variants, is challenging. The Critical Assessment of Genome Interpretation 5 SickKids clinical genome challenge provided an opportunity to assess our ability to extract such information. Participants in the challenge were required to match each of the 24 whole-genome sequences to the correct phenotypic profile and to identify the disease class of each genome. These are all rare disease cases that have resisted genetic diagnosis in a state-of-the-art pipeline. The patients have a range of eye, neurological, and connective-tissue disorders. We used a gene-centric approach to address this problem, assigning each gene a multiphenotype-matching score. Mutations in the top-scoring genes for each phenotype profile were ranked on a 6-point scale of pathogenicity probability, resulting in an approximately equal number of top-ranked coding and noncoding candidate variants overall. We were able to assign the correct disease class for 12 cases and the correct genome to a clinical profile for five cases. The challenge assessor found genes in three of these five cases as likely appropriate. In the postsubmission phase, after careful screening of the genes in the correct genome, we identified additional potential diagnostic variants, a high proportion of which are noncoding.
Collapse
Affiliation(s)
- Lipika R. Pal
- Institute for Bioscience and Biotechnology Research, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
| | - Kunal Kundu
- Institute for Bioscience and Biotechnology Research, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
- Computational Biology, Bioinformatics and Genomics, Biological Sciences Graduate Program, University of Maryland, College Park, MD 20742, USA
| | - Yizhou Yin
- Institute for Bioscience and Biotechnology Research, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
| | - John Moult
- Institute for Bioscience and Biotechnology Research, University of Maryland, 9600 Gudelsky Drive, Rockville, MD 20850, USA
- Department of Cell Biology and Molecular Genetics, University of Maryland, College Park, MD 20742, USA
| |
Collapse
|
49
|
Goonetilleke SN, Croxford AE, March TJ, Wirthensohn MG, Hrmova M, Mather DE. Variation among S-locus haplotypes and among stylar RNases in almond. Sci Rep 2020; 10:583. [PMID: 31953457 PMCID: PMC6969032 DOI: 10.1038/s41598-020-57498-6] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/22/2019] [Accepted: 12/31/2019] [Indexed: 11/23/2022] Open
Abstract
In many plant species, self-incompatibility systems limit self-pollination and mating among relatives. This helps maintain genetic diversity in natural populations but imposes constraints in agriculture and plant breeding. In almond [Prunus dulcis (Mill.) D.A. Webb], the specificity of self-incompatibility is mainly determined by stylar ribonuclease (S-RNase) and S-haplotype-specific F-box (SFB) proteins, both encoded within a complex locus, S. Prior to this research, a nearly complete sequence was available for one S-locus haplotype. Here, we report complete sequences for four haplotypes and partial sequences for 11 haplotypes. Haplotypes vary in sequences of genes (particularly S-RNase and SFB), distances between genes and numbers and positions of long terminal repeat transposons. Haplotype variation outside of the S-RNase and SFB genes may help maintain functionally important associations between S-RNase and SFB alleles. Fluorescence-based assays were developed to distinguish among some S-RNase alleles. With three-dimensional modelling of five S-RNase proteins, conserved active sites were identified and variation was observed in electrostatic potential and in the numbers, characteristics and positions of secondary structural elements, loop anchoring points and glycosylation sites. A hypervariable region on the protein surface and differences in the number, location and types of glycosylation sites may contribute to determining S-RNase specificity.
Collapse
Affiliation(s)
- Shashi N Goonetilleke
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia
| | - Adam E Croxford
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia
| | - Timothy J March
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia
| | - Michelle G Wirthensohn
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia
| | - Maria Hrmova
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia.,School of Life Sciences, Huaiyin Normal University, Huai'an, 223300, China
| | - Diane E Mather
- School of Agriculture, Food and Wine, Waite Research Institute, The University of Adelaide, PMB 1, Glen Osmond, SA, 5064, Australia.
| |
Collapse
|
50
|
Yu TG, Kim HS, Choi Y. B-SIDER: Computational Algorithm for the Design of Complementary β-Sheet Sequences. J Chem Inf Model 2019; 59:4504-4511. [PMID: 31512871 DOI: 10.1021/acs.jcim.9b00548] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The β-sheet is an element of protein secondary structure, and intra-/intermolecular β-sheet interactions play pivotal roles in biological regulatory processes including scaffolding, transporting, and oligomerization. In nature, a β-sheet formation is tightly regulated because dysregulated β-stacking often leads to severe diseases such as Alzheimer's, Parkinson's, systemic amyloidosis, or diabetes. Thus, the identification of intrinsic β-sheet-forming propensities can provide valuable insight into protein designs for the development of novel therapeutics. However, structure-based design methods may not be generally applicable to such amyloidogenic peptides mainly owing to high structural plasticity and complexity. Therefore, an alternative design strategy based on complementary sequence information is of significant importance. Herein, we developed a database search method called β-Stacking Interaction DEsign for Reciprocity (B-SIDER) for the design of complementary β-strands. This method makes use of the structural database information and generates target-specific score matrices. The discriminatory power of the B-SIDER score function was tested on representative amyloidogenic peptide substructures against a sequence-based score matrix (PASTA 2.0) and two popular ab initio protein design score functions (Rosetta and FoldX). B-SIDER is able to distinguish wild-type amyloidogenic β-strands as favored interactions in a more consistent manner than other methods. B-SIDER was prospectively applied to the design of complementary β-strands for a splitGFP scaffold. Three variants were identified to have stronger interactions than the original sequence selected through a directed evolution, emitting higher fluorescence intensities. Our results indicate that B-SIDER can be applicable to the design of other β-strands, assisting in the development of therapeutics against disease-related amyloidogenic peptides.
Collapse
Affiliation(s)
- Tae-Geun Yu
- Department of Biological Sciences , Korea Advanced Institute of Science and Technology (KAIST) , Daejeon 34141 , Republic of Korea
| | - Hak-Sung Kim
- Department of Biological Sciences , Korea Advanced Institute of Science and Technology (KAIST) , Daejeon 34141 , Republic of Korea
| | - Yoonjoo Choi
- Department of Biological Sciences , Korea Advanced Institute of Science and Technology (KAIST) , Daejeon 34141 , Republic of Korea
| |
Collapse
|