Search PubMed⌕ Search
Home/Archive/PubMed

PubMed · 9796822

Optimizing potentials for the inverse protein folding problem.

T L Chiu·R A Goldstein
View original source ↗Export metadata ↓

Abstract

Inverse protein folding, which seeks to identify sequences that fold into a given structure, has been approached by threading candidate sequences onto the structure and scoring them with database-derived potentials. The sequences with the lowest energies are predicted to fold into that structure. It has been argued that the limited success of this type of approach is not due to the discrepancy between the scoring potential and the true potential but is rather due to the fact that sequences choose their lowest-energy structure rather than structures choosing the lowest-energy sequences. Here we develop a non-physical potential scheme optimized for the inverse folding problem. We maximize the average probability of success for a set of lattice proteins to obtain the optimal potential energy function, and show that the potential obtained by our method is more likely to produce successful predictions than the true potential.

Explore related subjects

Protein FoldingProteinsJournal ArticleResearch Support, Non-U.S. Gov'tResearch Support, U.S. Gov't, Non-P.H.S.Research Support, U.S. Gov't, P.H.S.

Keep this discovery

Explore connections, maps & timelines

BibTeXRIS

T L Chiu, R A Goldstein. 1998. Optimizing potentials for the inverse protein folding problem.. https://doi.org/10.1093/protein%2F11.9.749

Cite the original work for its findings. Save a collection to share your selection of sources.

KEEP EXPLORING

Related citations

PubMed2025

NumSimEX: A method using EXX hydrogen exchange mass spectrometry to map the energetics of protein folding landscapes.

Hydrogen exchange mass spectrometry (HXMS) is a powerful tool to understand protein folding pathways and energetics. However, HXMS experiments to date have used exchange conditions termed EX1 or EX2 which limit the information that can be gained compared to the more general EXX exchange regime. If EXX behavior could be understood and analyzed, a single HXMS timecourse on an intact protein could fully map its folding landscape without requiring denaturation. To address this challenge, we developed a numerical simulation method called NumSimEX that models EXX exchange for arbitrarily complex folding pathways. NumSimEx fits protein folding dynamics to experimental HXMS data by iteratively comparing the simulated and experimental timecourses, allowing for determination of both kinetic and thermodynamic protein folding parameters. After analytically verifying NumSimEX's accuracy, we demonstrated its power on HXMS data from beta-2 microglobulin (β2M), a protein involved in dialysis-related amyloidosis. In particular, using NumSimEX, we identified three-state kinetics that near-perfectly matched experimental observation. This proof-of-principle application of NumSimEX sets the stage for harnessing HXMS to expand our understanding of proteins currently excluded from traditional protein folding methods. NumSimEX is freely available at https://github.com/JaswalLab/NumSimEX_Public.

Protein Folding↗
PubMed2006

A PDZ domain recapitulates a unifying mechanism for protein folding.

A unifying view has been recently proposed according to which the classical diffusion-collision and nucleation-condensation models may represent two extreme manifestations of an underlying common mechanism for the folding of small globular proteins. We report here the characterization of the folding process of the PDZ domain, a protein that recapitulates the three canonical steps involved in this unifying mechanism, namely: (i) the early formation of a weak nucleus that determines the native-like topology of a large portion of the structure, (ii) a global collapse of the entire polypeptide chain, and (iii) the consolidation of the remaining partially structured regions to achieve the native state conformation. These steps, which are clearly detectable in the PDZ domain investigated here, may be difficult to distinguish experimentally in other proteins, which would thus appear to follow one of the two limiting mechanisms. The analysis of the (un)folding kinetics for other three-state proteins (when available) appears consistent with the predictions ensuing from this unifying mechanism, thus providing a powerful validation of its general nature.

Protein Folding↗
PubMed2006

Exploring the energy landscape of protein folding using replica-exchange and conventional molecular dynamics simulations.

Two independent replica-exchange molecular dynamics (REMD) simulations with an explicit water model were performed of the Trp-cage mini-protein. In the first REMD simulation, the replicas started from the native conformation, while in the second they started from a nonnative conformation. Initially, the first simulation yielded results qualitatively similar to those of two previously published REMD simulations: the protein appeared to be over-stabilized, with the predicted melting temperature 50-150K higher than the experimental value of 315K. However, as the first REMD simulation progressed, the protein unfolded at all temperatures. In our second REMD simulation, which starts from a nonnative conformation, there was no evidence of significant folding. Transitions from the unfolded to the folded state did not occur on the timescale of these simulations, despite the expected improvement in sampling of REMD over conventional molecular dynamics (MD) simulations. The combined 1.42 micros of simulation time was insufficient for REMD simulations with different starting structures to converge. Conventional MD simulations at a range of temperatures were also performed. In contrast to REMD, the conventional MD simulations provide an estimate of Tm in good agreement with experiment. Furthermore, the conventional MD is a fraction of the cost of REMD and continuous, realistic pathways of the unfolding process at atomic resolution are obtained.

Protein Folding↗