Dataset Information

Assessment of the utility of contact-based restraints in accelerating the prediction of protein structure using molecular dynamics simulations.

ABSTRACT: Molecular dynamics (MD) simulation is a well-established tool for the computational study of protein structure and dynamics, but its application to the important problem of protein structure prediction remains challenging, in part because extremely long timescales can be required to reach the native structure. Here, we examine the extent to which the use of low-resolution information in the form of residue-residue contacts, which can often be inferred from bioinformatics or experimental studies, can accelerate the determination of protein structure in simulation. We incorporated sets of 62, 31, or 15 contact-based restraints in MD simulations of ubiquitin, a benchmark system known to fold to the native state on the millisecond timescale in unrestrained simulations. One-third of the restrained simulations folded to the native state within a few tens of microseconds-a speedup of over an order of magnitude compared with unrestrained simulations and a demonstration of the potential for limited amounts of structural information to accelerate structure determination. Almost all of the remaining ubiquitin simulations reached near-native conformations within a few tens of microseconds, but remained trapped there, apparently due to the restraints. We discuss potential methodological improvements that would facilitate escape from these near-native traps and allow more simulations to quickly reach the native state. Finally, using a target from the Critical Assessment of protein Structure Prediction (CASP) experiment, we show that distance restraints can improve simulation accuracy: In our simulations, restraints stabilized the native state of the protein, enabling a reasonable structural model to be inferred.

SUBMITTER: Raval A

PROVIDER: S-EPMC4815320 | biostudies-literature | 2016 Jan

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Assessment of the utility of contact-based restraints in accelerating the prediction of protein structure using molecular dynamics simulations.

Raval Alpan A Piana Stefano S Eastwood Michael P MP Shaw David E DE

Protein science : a publication of the Protein Society 20150830 1

Molecular dynamics (MD) simulation is a well-established tool for the computational study of protein structure and dynamics, but its application to the important problem of protein structure prediction remains challenging, in part because extremely long timescales can be required to reach the native structure. Here, we examine the extent to which the use of low-resolution information in the form of residue-residue contacts, which can often be inferred from bioinformatics or experimental studies, ...[more]

PMID: 26266489

Similar Datasets

Project description:The G-protein coupled receptor, GPR120, has ubiquitous expression and multifaceted roles in modulating metabolic and anti-inflammatory processes. Recent implications of its role in cancer progression have presented GPR120 as an attractive oncogenic drug target. GPR120 gene knockdown in breast cancer studies revealed a role of GPR120-induced chemoresistance in epirubicin and cisplatin-induced DNA damage in tumour cells. Higher expression and activation levels of GPR120 is also reported to promote tumour angiogenesis and cell migration in colorectal cancer. Some agonists targeting GPR120 have been reported, such as TUG891 and Compound39, but to date development of small-molecule inhibitors of GPR120 is limited. Herein, following homology modelling of the receptor a pharmacophore hypothesis was derived from 300 ns all-atomic molecular dynamics (MD) simulations on apo, TUG891-bound and Compound39-bound GPR120S (short isoform) receptor models embedded in a water solvated lipid bilayer system. We performed comparative MD analysis on protein-ligand interactions between the two agonist and apo simulations on the stability of the "ionic lock" - a Class A GPCRs characteristic of receptor activation and inactivation. The detailed analysis predicted that ligand interactions with W277 and N313 are critical to conserve the "ionic-lock" conformation (R136 of Helix 3) and prevent GPR120S receptor activation. The results led to generation of a W277 and N313 focused pharmacophore hypothesis and the screening of the ZINC15 database using ZINCPharmer through the structure-based pharmacophore. 100 ns all-atomic molecular dynamics (MD) simulations were performed on 9 small molecules identified and Cpd 9, (2-hydroxy-N-{4-[(6-hydroxy-2-methylpyrimidin-4-yl) amino] phenyl} benzamide) was predicted to be a small-molecule GPR120S antagonist. The conformational results from the collective all-atomic MD analysis provided structural information for further identification and optimisation of novel druggable inhibitors of GPR120S using this rational design approach, which could have future potential for anti-cancer drug development studies.

Dataset Information

Assessment of the utility of contact-based restraints in accelerating the prediction of protein structure using molecular dynamics simulations.

Publications

Assessment of the utility of contact-based restraints in accelerating the prediction of protein structure using molecular dynamics simulations.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets