1
|
Cancela F, Rendon-Marin S, Quintero-Gil C, Houston DR, Gumbis G, Panzera Y, Pérez R, Arbiza J, Mirazo S. Modelling of Hepatitis E virus RNA-dependent RNA polymerase genotype 3 from a chronic patient and in silico interaction analysis by molecular docking with Ribavirin. J Biomol Struct Dyn 2023; 41:705-721. [PMID: 34861797 DOI: 10.1080/07391102.2021.2011416] [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] [Indexed: 01/04/2023]
Abstract
Hepatitis E Virus (HEV) infection is an emergent zoonotic disease, where chronic hepatitis E associated to solid organ transplant (SOT) recipients, related to genotype 3, is the clinical manifestation of major concern. In this setting, ribavirin (RBV) treatment is the only available therapy, though drug-resistant variants could emerge leading to a therapeutic failure. Crystallographic structures have not been reported for most of the HEV proteins, including the RNA-polymerase (RdRp). Therefore, the mechanism of action of RBV against HEV and the molecular interactions between this drug and RdRp are largely unknown. In this work, we aimed to model in silico the 3 D structure of a novel HEV3 RdRp (HEV_C1_Uy) from a chronically HEV infected-SOT recipient treated with RBV and to perform a molecular docking simulation between RBV triphosphate (RBVT), 7-methyl-guanosine-5'-triphosphate and the modelled protein. The models were generated using I-TASSER server and validated with multiple bioinformatics tools. The docking analysis were carried out with AutoDock Vina and LeDock software. We obtained a suitable model for HEV_C1_Uy (C-Score=-1.33, RMSD = 10.4 ± 4.6 Å). RBVT displayed a binding affinity of -7.6 ± 0.2 Kcal/mol by molecular docking, mediated by 6 hydrogen-bonds (Q195-O14, S198-O11, E257-O13, S260-O2, O3, S311-O11) between the finger's-palm-domains and a free binding energy of 31.26 ± 16.81 kcal/mol by molecular dynamics simulations. We identified the possible HEV RdRp interacting region for incoming nucleotides or analogs and provide novel insights that will contribute to better understand the molecular interactions of RBV and the enzyme and the mechanism of action of this antiviral drug.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Florencia Cancela
- Sección Virología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Santiago Rendon-Marin
- Grupo de Investigación en Ciencias Animales - GRICA, Facultad de Medicina Veterinaria y Zootecnia, Universidad Cooperativa de Colombia, sede Bucaramanga, Bucaramanga, Colombia
| | - Carolina Quintero-Gil
- Grupo de Investigación en Ciencias Animales - GRICA, Facultad de Medicina Veterinaria y Zootecnia, Universidad Cooperativa de Colombia, sede Bucaramanga, Bucaramanga, Colombia
| | - Douglas R Houston
- Institute of Quantitative Biology, Biochemistry and Biotechnology, The University of Edinburgh, Edinburgh, UK
| | - Gediminas Gumbis
- Institute of Quantitative Biology, Biochemistry and Biotechnology, The University of Edinburgh, Edinburgh, UK
| | - Yanina Panzera
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Ruben Pérez
- Sección Genética Evolutiva, Instituto de Biología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Juan Arbiza
- Sección Virología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay
| | - Santiago Mirazo
- Sección Virología, Facultad de Ciencias, Universidad de la República, Montevideo, Uruguay.,Departamento de Bacteriología y Virología, Instituto de Higiene, Facultad de Medicina, Universidad de la República, Montevideo, Uruguay
| |
Collapse
|
2
|
Ammar A, Cavill R, Evelo C, Willighagen E. PSnpBind: a database of mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow. J Cheminform 2022; 14:8. [PMID: 35227289 PMCID: PMC8886843 DOI: 10.1186/s13321-021-00573-5] [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: 08/11/2021] [Accepted: 11/18/2021] [Indexed: 11/15/2022] Open
Abstract
A key concept in drug design is how natural variants, especially the ones occurring in the binding site of drug targets, affect the inter-individual drug response and efficacy by altering binding affinity. These effects have been studied on very limited and small datasets while, ideally, a large dataset of binding affinity changes due to binding site single-nucleotide polymorphisms (SNPs) is needed for evaluation. However, to the best of our knowledge, such a dataset does not exist. Thus, a reference dataset of ligands binding affinities to proteins with all their reported binding sites' variants was constructed using a molecular docking approach. Having a large database of protein-ligand complexes covering a wide range of binding pocket mutations and a large small molecules' landscape is of great importance for several types of studies. For example, developing machine learning algorithms to predict protein-ligand affinity or a SNP effect on it requires an extensive amount of data. In this work, we present PSnpBind: A large database of 0.6 million mutated binding site protein-ligand complexes constructed using a multithreaded virtual screening workflow. It provides a web interface to explore and visualize the protein-ligand complexes and a REST API to programmatically access the different aspects of the database contents. PSnpBind is open source and freely available at https://psnpbind.org .
Collapse
Affiliation(s)
- Ammar Ammar
- Department of Bioinformatics—BiGCaT, NUTRIM, Maastricht University, Maastricht, The Netherlands
| | - Rachel Cavill
- Department of Data Science and Knowledge Engineering, Maastricht University, Maastricht, The Netherlands
| | - Chris Evelo
- Department of Bioinformatics—BiGCaT, NUTRIM, Maastricht University, Maastricht, The Netherlands
| | - Egon Willighagen
- Department of Bioinformatics—BiGCaT, NUTRIM, Maastricht University, Maastricht, The Netherlands
| |
Collapse
|
3
|
Tarapara B, Shah F. An in-silico analysis to identify structural, functional and regulatory role of SNPs in hMRE11. J Biomol Struct Dyn 2022; 41:2160-2174. [PMID: 35048780 DOI: 10.1080/07391102.2022.2028678] [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/19/2022]
Abstract
Meiotic recombination 11 (MRE11) is a component of the tri-molecular MRE11-RAD50-NBS1 (MRN) complex, which functions as an exonuclease and endonuclease which is involved in identifying, signalling, protecting and repairing double-strand breaks in DNA (DSBs). Ataxia-telangiectasia-like disorder (ATLD) 1 and Nijmegen breakage syndrome (NBS)-like disorder are MRE11 associated diseases. In the present study, we used an integrated computational approach to identify the most deleterious SNPs and their structural and functional impact on human MRE11. Five of the 68 observed non-synonymous SNP (nsSNPs; I162T, S273C, W210C, D311Y and R364L) should be worked on due to their strong possible pathogenicity and the risk of changing protein properties. All the nsSNPs were highly conserved and decrease the protein stability located in the MRE11 nuclease and MRE11 DNA binding presumed domain. R364L and I162T were predicted to be involved in post-translational modification (PTM) sites. Furthermore, we also analysed the regulatory effect of noncoding SNPs on MRE11 gene regulation in which 6 SNPs were found to affect gene regulation. All six noncoding SNPs predicted chromatin interactive site whereas only one SNP was noted its association with miRNA binding site which disrupts 5 miRNA conserved site. These findings help future studies to get more insights into the role of these variants in the alteration of the MRE11 function. Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Bhoomi Tarapara
- Department of Cancer Biology, Stem Cell Biology Lab, The Gujarat Cancer and Research Institute, Ahmedabad, India
| | - Franky Shah
- Department of Cancer Biology, Stem Cell Biology Lab, The Gujarat Cancer and Research Institute, Ahmedabad, India
| |
Collapse
|
4
|
Lam SD, Babu MM, Lees J, Orengo CA. Biological impact of mutually exclusive exon switching. PLoS Comput Biol 2021; 17:e1008708. [PMID: 33651795 PMCID: PMC7954323 DOI: 10.1371/journal.pcbi.1008708] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 03/12/2021] [Accepted: 01/14/2021] [Indexed: 12/27/2022] Open
Abstract
Alternative splicing can expand the diversity of proteomes. Homologous mutually exclusive exons (MXEs) originate from the same ancestral exon and result in polypeptides with similar structural properties but altered sequence. Why would some genes switch homologous exons and what are their biological impact? Here, we analyse the extent of sequence, structural and functional variability in MXEs and report the first large scale, structure-based analysis of the biological impact of MXE events from different genomes. MXE-specific residues tend to map to single domains, are highly enriched in surface exposed residues and cluster at or near protein functional sites. Thus, MXE events are likely to maintain the protein fold, but alter specificity and selectivity of protein function. This comprehensive resource of MXE events and their annotations is available at: http://gene3d.biochem.ucl.ac.uk/mxemod/. These findings highlight how small, but significant changes at critical positions on a protein surface are exploited in evolution to alter function.
Collapse
Affiliation(s)
- Su Datt Lam
- Institute of Structural and Molecular Biology, University College London, Darwin Building, Gower Street, London, United Kingdom
- Department of Applied Physics, Faculty of Science and Technology, Universiti Kebangsaan Malaysia, Bangi, Malaysia
- * E-mail: (SDL); (JL); (CO)
| | - M. Madan Babu
- MRC Laboratory of Molecular Biology, Francis Crick Avenue, Cambridge Biomedical Campus, Cambridge, United Kingdom
- Department of Structural Biology and Center for Data Driven Discovery, St Jude Children’s Research Hospital, Memphis, Tennessee, United States of America
| | - Jonathan Lees
- Faculty of Health and Life Sciences, Oxford Brookes University, Oxford, United Kingdom
- * E-mail: (SDL); (JL); (CO)
| | - Christine A. Orengo
- Institute of Structural and Molecular Biology, University College London, Darwin Building, Gower Street, London, United Kingdom
- * E-mail: (SDL); (JL); (CO)
| |
Collapse
|
5
|
Mohamadian M, Ghandil P, Naseri M, Bahrami A, Momen AA. A novel homozygous variant in an Iranian pedigree with cerebellar ataxia, mental retardation, and dysequilibrium syndrome type 4. J Clin Lab Anal 2020; 34:e23484. [PMID: 33079427 PMCID: PMC7676196 DOI: 10.1002/jcla.23484] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2020] [Revised: 06/12/2020] [Accepted: 06/26/2020] [Indexed: 01/20/2023] Open
Abstract
BACKGROUND Cerebellar ataxia, mental retardation, and dysequilibrium (CAMRQ) syndrome is a rare and early-onset neurodevelopmental disorder. Four subtypes of this syndrome have been identified, which are clinically and genetically different. To date, altogether 32 patients have been described with ATP8A2 mutations and phenotypic features assigned to CAMRQ type 4. Herein, three additional patients in an Iranian consanguineous family with non-progressive cerebellar ataxia, severe hypotonia, intellectual disability, dysarthria, and cerebellar atrophy have been identified. METHODS Following the thorough clinical examination, consecutive detections including chromosome karyotyping, chromosomal microarray analysis, and whole exome sequencing (WES) were performed on the proband. The sequence variants derived from WES interpreted by a standard bioinformatics pipeline. Pathogenicity assessment of candidate variant was done by in silico analysis. The familial cosegregation of the WES finding was carried out by PCR-based Sanger sequencing. RESULTS A novel homozygous missense variant (c.1339G > A, p.Gly447Arg) in the ATP8A2 gene was identified and completely segregated with the phenotype in the family. In silico analysis and structural modeling revealed that the p.G477R substitution is deleterious and induced undesired effects on the protein stability and residue distribution in the ligand-binding pocket. The novel sequence variant occurred within an extremely conserved subregion of the ATP-binding domain. CONCLUSION Our findings expand the spectrum of ATP8A2 mutations and confirm the reported genotype-phenotype correlation. These results could improve genetic counseling and prenatal diagnosis in families with clinical presentations related to CAMRQ4 syndrome.
Collapse
Affiliation(s)
- Malihe Mohamadian
- Department of Molecular Medicine, Birjand University of Medical Sciences, Birjand, Iran
| | - Pegah Ghandil
- Diabetes Research Center, Health Research Institute, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran.,Department of Medical Genetics, School of Medicine, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| | - Mohsen Naseri
- Cellular and Molecular Research Center, Birjand University of Medical Sciences, Birjand, Iran
| | - Afsane Bahrami
- Cellular and Molecular Research Center, Birjand University of Medical Sciences, Birjand, Iran
| | - Ali Akbar Momen
- Department of Paediatric Neurology, Golestan Medical, Educational, and Research Center, Ahvaz Jundishapur University of Medical Sciences, Ahvaz, Iran
| |
Collapse
|
6
|
Yura K. Preface of Special Issue "Protein-Ligand Interactions". Biophys Physicobiol 2016; 13:85-86. [PMID: 27924261 PMCID: PMC5042156 DOI: 10.2142/biophysico.13.0_85] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/04/2022] Open
Affiliation(s)
- Kei Yura
- Graduate School of Humanities and Sciences, Ochanomizu University, 2-1-1 Otsuka, Bunkyo, Tokyo 112-8610, Japan; Center for Informational Biology, Ochanomizu University, 2-1-1 Otsuka, Bunkyo, Tokyo 112-8610, Japan; National Institute of Genetics, 1111 Yata, Mishima, Shizuoka 411-8540, Japan
| |
Collapse
|