58 research outputs found
<i>In Silico</i>Tools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review
No full textOzbek, Pemra/0000-0002-3043-0015; Yazar, Metin/0000-0002-2657-3072Single-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines thein silicotools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination ofin silicoapproaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions madein silico.[3454]TUSEB project number 3454 is kindly acknowledged
Dynamic characterization of HLA-B*44 Alleles: A comparative molecular dynamics simulation study
No full textHuman Leukocyte Antigens (HLA) are highly polymorphic proteins that play a key role in the immune system. HLA molecule is present on the cell membrane of antigen-presenting cells of the immune system and presents short peptides, originating from the proteins of invading pathogens or self-proteins, to the T-cell Receptor (TCR) molecule of the T-cells. In this study, peptide-binding characteristics of HLAB*44:02, 44:03, 44:05 alleles bound to three nonameric peptides were studied using molecular dynamics simulations. Polymorphisms among these alleles (Asp116Tyr and Asp156Leu) result in major differences in the allele characteristics. While HLA-B*44:02 (Asp116, Asp156) and HLA-B*44:03 (Asp116, Leu156) depend on tapasin for efficient peptide loading, HLA-B*44:05 (Tyr116, Asp156) is tapasin independent. On the other hand, HLA-B*44:02 and HLA-B*44:03 mismatch is closely related to transplant rejection and acute-graft-versus-host disease. In order to understand the dynamic characteristics, the simulation trajectories were analyzed by applying Root Mean Square Deviation (RMSD) and Root Mean Square Fluctuation (RMSF) calculations and hydrogen bonding analysis. Binding dynamics of the three HLAB*44 alleles and peptide sequences are comparatively discussed. In general, peptide binding stability is found to depend on the peptide rather than the allele type for HLA-B*44 alleles. (C) 2016 Elsevier Ltd. All rights reserved
Revisiting allostery in CREB-binding protein (CBP) using residue-based interaction energy
No full textYazar, Metin/0000-0002-2657-3072; Ozbek, Pemra/0000-0002-3043-0015CREB-binding protein (CBP) is a multi-subunit scaffold protein complex in transcription regulation process, binding and interacting with ligands such as mixed-lineage leukemia (MLL) and c-Myb allosterically. Here in this study, we have revisited the concept of allostery in CBP via residue-based interaction energy calculation based on molecular dynamics (MD) simulations. To this end, we conducted MD simulations of KIX:MLL:c-Myb ternary complex, its binary components and kinase-inducible domain (KID) interacting domain (KIX) backbone. Interaction energy profiles and cross correlation analysis were performed and the results indicated that KIX:MLL and KIX:c-Myb:MLL complexes demonstrate significant similarities according to both analysis methods. Two regions in the KIX backbone were apparent from the interaction energy and cross correlation maps that hold a key to allostery phenomena observed in CBP. While one of these regions are related to the ligand binding residues, the other comprises of L-12-G(2) loop and alpha(3) helix regions that have been found to have a significant role in allosteric signal propagation. All in all, residue-based interaction energy calculation method is demonstrated to be a valuable calculation technique for the detection of allosteric signal propagation and ligand interaction regions.TUBITAK [218M861]TUBITAK 218M861 is acknowledged. Fruitful discussions and coding support from Onur Sercinoglu and Elif Naz Bingol are greatly acknowledged
Assessment of 13 <i>in silico</i> pathogenicity methods on cancer-related variants
No full textYazar, Metin/0000-0002-2657-3072; Ozbek, Pemra/0000-0002-3043-0015Single nucleotide variants (SNVs) are single base substitutions that could influence many biological functions in the cell including gene expression, protein folding, and protein-protein interactions among many others. Thus, predictions of functional effects of cancer-related variants are crucial for drug responses and treatment options in clinical oncology. Experimental identification of these effects could be slow, inefficient, and inconvenient, hence in silico methods are gaining popularity in predicting the variants' effects. There are many studies on the cancer variants, however, up to date, none of these have been aimed to assess the performance metrics of in silico pathogenicity methods on functional relevance of cancer variants obtained from ClinVar. To this end, we examined the pathogenicity predictions of cancer-related variant datasets of 8 cancer types (bladder, breast, colon, colorectal, kidney, liver, lung, and pancreas cancer) retrieved from ClinVar using 13 different in silico methods including SIFT, CADD, FATHMM-weighted, FATHMM-unweighted, GERP(++), MetaSVM, Mutation Assessor, MutationTaster, MutPred, PolyPhen-2, Provean, Revel and VEST4. A combination of statistical performance metric analysis, prediction distribution frequency data and ROC curve analysis results have suggested that; among all in silico prediction tools, top three tools with the highest discriminatory power were found to be MutPred (AUC = 0.677), MetaSVM (AUC = 0.645) and Revel (AUC = 0.637).[3454]PO acknowledges TUSEB Project number 3454
Sequence-structure-function relationships in class I MHC: A local frustration perspective
No full textClass I Major Histocompatibility Complex (MHC) binds short antigenic peptides with the help of Peptide Loading Complex (PLC), and presents them to T-cell Receptors (TCRs) of cytotoxic T-cells and Killer-cell Immunglobulin-like Receptors (KIRs) of Natural Killer (NK) cells. With more than 10000 alleles, human MHC (Human Leukocyte Antigen, HLA) is the most polymorphic protein in humans. This allelic diversity provides a wide coverage of peptide sequence space, yet does not affect the three-dimensional structure of the complex. Moreover, TCRs mostly interact with HLA in a common diagonal binding mode, and KIR-HLA interaction is allele-dependent. With the aim of establishing a framework for understanding the relationships between polymorphism (sequence), structure (conserved fold) and function (protein interactions) of the human MHC, we performed here a local frustration analysis on pMHC homology models covering 1436 HLA I alleles. An analysis of local frustration profiles indicated that (1) variations in MHC fold are unlikely due to minimally-frustrated and relatively conserved residues within the HLA peptide-binding groove, (2) high frustration patches on HLA helices are either involved in or near interaction sites of MHC with the TCR, KIR, or tapasin of the PLC, and (3) peptide ligands mainly stabilize the F-pocket of HLA binding groove
gRINN: a tool for calculation of residue interaction energies and protein energy network analysis of molecular dynamics simulations
No full textAtomistic molecular dynamics (MD) simulations generate a wealth of information related to the dynamics of proteins. If properly analyzed, this information can lead to new insights regarding protein function and assist wet-lab experiments. Aiming to identify interactions between individual amino acid residues and the role played by each in the context of MD simulations, we present a stand-alone software called gRINN (get Residue Interaction eNergies and Networks). gRINN features graphical user interfaces (GUIs) and a command-line interface for generating and analyzing pairwise residue interaction energies and energy correlations from protein MD simulation trajectories. gRINN utilizes the features of NAMD or GROMACS MD simulation packages and automatizes the steps necessary to extract residue-residue interaction energies from user-supplied simulation trajectories, greatly simplifying the analysis for the end-user. A GUI, including an embedded molecular viewer, is provided for visualization of interaction energy time-series, distributions, an interaction energy matrix, interaction energy correlations and a residue correlation matrix. gRINN additionally offers construction and analysis of Protein Energy Networks, providing residue-based metrics such as degrees, betweenness-centralities, closeness centralities as well as shortest path analysis. gRINN is free and open to all users without login requirement at http://grinn.readthedocs.io
Revisiting allostery in CREB-binding protein (CBP) using residue-based interaction energy
No full textCREB-binding protein (CBP) is a multi-subunit scaffold protein complex in transcription regulation process, binding and interacting with ligands such as mixed-lineage leukemia (MLL) and c-Myb allosterically. Here in this study, we have revisited the concept of allostery in CBP via residue-based interaction energy calculation based on molecular dynamics (MD) simulations. To this end, we conducted MD simulations of KIX:MLL:c-Myb ternary complex, its binary components and kinase-inducible domain (KID) interacting domain (KIX) backbone. Interaction energy profiles and cross correlation analysis were performed and the results indicated that KIX:MLL and KIX:c-Myb:MLL complexes demonstrate significant similarities according to both analysis methods. Two regions in the KIX backbone were apparent from the interaction energy and cross correlation maps that hold a key to allostery phenomena observed in CBP. While one of these regions are related to the ligand binding residues, the other comprises of L-12-G(2) loop and alpha(3) helix regions that have been found to have a significant role in allosteric signal propagation. All in all, residue-based interaction energy calculation method is demonstrated to be a valuable calculation technique for the detection of allosteric signal propagation and ligand interaction regions
In SilicoTools and Approaches for the Prediction of Functional and Structural Effects of Single-Nucleotide Polymorphisms on Proteins: An Expert Review
No full textSingle-nucleotide polymorphisms (SNPs) are single-base variants that contribute to human biological variation and pathogenesis of many human diseases. Among all SNP types, nonsynonymous single-nucleotide polymorphisms (nsSNPs) can alter many structural, biochemical, and functional features of a protein such as folding characteristics, charge distribution, stability, dynamics, and interactions with other proteins/nucleotides. These modifications in the protein structure can lead nsSNPs to be closely associated with many multifactorial diseases such as cancer, diabetes, and neurodegenerative diseases. Predicting structural and functional effects of nsSNPs with experimental approaches can be time-consuming and costly; hence, computational prediction tools and algorithms are being widely and increasingly utilized in biology and medical research. This expert review examines thein silicotools and algorithms for the prediction of functional or structural effects of SNP variants, in addition to the description of the phenotypic effects of nsSNPs on protein structure, association between pathogenicity of variants, and functional or structural features of disease-associated variants. Finally, case studies investigating the functional and structural effects of nsSNPs on selected protein structures are highlighted. We conclude that creating a consistent workflow with a combination ofin silicoapproaches or tools should be considered to increase the performance, accuracy, and precision of the biological and clinical predictions madein silico
Computational characterization of residue couplings and micropolymorphism-induced changes in the dynamics of two differentially disease-associated human MHC class-I alleles
No full textHuman major histocompatibility complex class I (MHC I) - or human leukocyte antigen (HLA) - proteins present intracellularly processed peptides to cytotoxic T lymphocytes in the adaptive immune response to pathogens. A high level of polymorphism in human MHC I proteins defines the peptide-binding specificity of thousands of different MHC alleles. However, polymorphism as well as the peptide ligand can also affect the global dynamics of the complex. In this study, we conducted classical molecular dynamics simulations of two HLA alleles, the ankylosing spondylitis (AS) associated/tapasin-dependent HLA-B*27:05 and nondisease-associated/tapasin-independent HLA-B*27:09, both in peptide-free forms as well as complex with four different peptides ligands. Our results indicate that in peptide-free form, the single amino acid substitution distinguishing the two alleles (D116H), leads to a weaker dynamic coupling of residues in the tapasin-dependent HLA-B*27:05. In peptide-bound form, several residues of the binding-groove, mostly in A and B pockets, show hinge-like behavior in the global motion of the MHC. Moreover, allele-dependent changes are shown in residue interactions, affecting the B-pocket as well as the beta-2-microglobulin (2m)-facing residues of the HLA chain
- …
