Dataset Information

Exploring the use of molecular dynamics in assessing protein variants for phenotypic alterations.

ABSTRACT: With the advent of rapid sequencing technologies, making sense of all the genomic variations that we see among us has been a major challenge. A plethora of algorithms and methods exist that try to address genome interpretation through genotype-phenotype linkage analysis or evaluating the loss of function/stability mutations in protein. Critical Assessment of Genome Interpretation (CAGI) offers an exceptional platform to blind-test all such algorithms and methods to assess their true ability. We take advantage of this opportunity to explore the use of molecular dynamics simulation as a tool to assess alteration of phenotype, loss of protein function, interaction, and stability. The results show that coarse-grained dynamics based protein flexibility analysis on 34 CHEK2 and 1719 CALM1 single mutants perform reasonably well for class-based predictions for phenotype alteration and two-thirds of the predicted scores return a correlation coefficient of 0.6 or more. When all-atom dynamics is used to predict altered stability due to mutations for Frataxin protein (8 cases), the predictions are comparable to the state-of-the-art methods. The competitive performance of our straightforward approach to phenotype interpretation contrasts with heavily trained machine learning approaches, and open new avenues to rationally improve genome interpretation.

SUBMITTER: Garg A

PROVIDER: S-EPMC7318789 | biostudies-literature | 2019 Sep

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Exploring the use of molecular dynamics in assessing protein variants for phenotypic alterations.

Garg Aditi A Pal Debnath D

Human mutation 20190712 9

With the advent of rapid sequencing technologies, making sense of all the genomic variations that we see among us has been a major challenge. A plethora of algorithms and methods exist that try to address genome interpretation through genotype-phenotype linkage analysis or evaluating the loss of function/stability mutations in protein. Critical Assessment of Genome Interpretation (CAGI) offers an exceptional platform to blind-test all such algorithms and methods to assess their true ability. We ...[more]

PMID: 31106920

Similar Datasets

Project description:BackgroundSG2NA is a member of the striatin sub-family of WD-40 repeat proteins. Striatin family members have been associated with diverse physiological functions. SG2NA has also been shown to have roles in cell cycle progression, signal transduction etc. They have been known to interact with a number of proteins including Caveolin and Calmodulin and also propagate the formation of a multimeric protein unit called striatin-interacting phosphatase and kinase. As a pre-requisite for such interaction ability, these proteins are known to be unstable and primarily disordered in their arrangement. Earlier we had identified that it has multiple isoforms (namely 35, 78, 87 kDa based on its molecular weight) which are generated by alternative splicing. However, detailed structural information of SG2NA is still eluding the researchers.ResultsThis study was aimed towards three-dimensional molecular modeling and characterization of SG2NA protein and its isoforms. One structure out of five was selected for each variant having the least value for C score. Out of these, m35 kDa with a C score value of -3.21 was the most poorly determined structure in comparison to m78 kDa and m87 kDa variants with C scores of -1.16 and -1.97 respectively. Further evaluation resulted in about 61.6% residues of m35 kDa, 76.6% residues of m78 kDa and 72.1% residues of m87 kDa falling in the favorable regions of Ramchandran Plot. Molecular dynamics simulations were also carried out to obtain biologically relevant structural models and compared with previous atomic coordinates. N-terminal region of all variants was found to be highly disordered.ConclusionThis study provides first-hand detailed information to understand the structural conformation of SG2NA protein variants (m35 kDa, m78 kDa and m87 kDa). The WD-40 repeat domain was found to constitute antiparallel strands of β-sheets arranged circularly. This study elucidates the crucial structural features of SG2NA proteins which are involved in various protein-protein interactions and also reveals the extent of disorder present in the SG2NA structure crucial for excessive interaction and multimeric protein complexes. The study also potentiates the role of computational approaches for preliminary examination of unknown proteins in the absence of experimental information.

Project description:The significant work that has been invested toward understanding protein-protein interaction has not translated into significant advances in structure-based predictions. In particular redesigning protein surfaces to bind to unrelated receptors remains a challenge, partly due to receptor flexibility, which is often neglected in these efforts. In this work, we computationally graft the binding epitope of various small proteins obtained from the RCSB database to bind to barnase, lysozyme, and trypsin using a previously derived and validated algorithm. In an effort to probe the protein complexes in a realistic environment, all native and designer complexes were subjected to a total of nearly 400 ns of explicit-solvent molecular dynamics (MD) simulation. The MD data led to an unexpected observation: some of the designer complexes were highly unstable and decomposed during the trajectories. In contrast, the native and a number of designer complexes remained consistently stable. The unstable conformers provided us with a unique opportunity to define the structural and energetic factors that lead to unproductive protein-protein complexes. To that end we used free energy calculations following the MM-PBSA approach to determine the role of nonpolar effects, electrostatics and entropy in binding. Remarkably, we found that a majority of unstable complexes exhibited more favorable electrostatics than native or stable designer complexes, suggesting that favorable electrostatic interactions are not prerequisite for complex formation between proteins. However, nonpolar effects remained consistently more favorable in native and stable designer complexes reinforcing the importance of hydrophobic effects in protein-protein binding. While entropy systematically opposed binding in all cases, there was no observed trend in the entropy difference between native and designer complexes. A series of alanine scanning mutations of hot-spot residues at the interface of native and designer complexes showed less than optimal contacts of hot-spot residues with their surroundings in the unstable conformers, resulting in more favorable entropy for these complexes. Finally, disorder predictions revealed that secondary structures at the interface of unstable complexes exhibited greater disorder than the stable complexes.

Dataset Information

Exploring the use of molecular dynamics in assessing protein variants for phenotypic alterations.

Publications

Exploring the use of molecular dynamics in assessing protein variants for phenotypic alterations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets