Browse
Submit Data
Databases
API
Help

Dataset Information

20 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Funneled energy landscape unifies principles of protein binding and evolution.

ABSTRACT: Most proteins have evolved to spontaneously fold into native structure and specifically bind with their partners for the purpose of fulfilling biological functions. According to Darwin, protein sequences evolve through random mutations, and only the fittest survives. The understanding of how the evolutionary selection sculpts the interaction patterns for both biomolecular folding and binding is still challenging. In this study, we incorporated the constraint of functional binding into the selection fitness based on the principle of minimal frustration for the underlying biomolecular interactions. Thermodynamic stability and kinetic accessibility were derived and quantified from a global funneled energy landscape that satisfies the requirements of both the folding into the stable structure and binding with the specific partner. The evolution proceeds via a bowl-like evolution energy landscape in the sequence space with a closed-ring attractor at the bottom. The sequence space is increasingly reduced until this ring attractor is reached. The molecular-interaction patterns responsible for folding and binding are identified from the evolved sequences, respectively. The residual positions participating in the interactions responsible for folding are highly conserved and maintain the hydrophobic core under additional evolutionary constraints of functional binding. The positions responsible for binding constitute a distributed network via coupling conservations that determine the specificity of binding with the partner. This work unifies the principles of protein binding and evolution under minimal frustration and sheds light on the evolutionary design of proteins for functions.

SUBMITTER: Yan Z

PROVIDER: S-EPMC7959555 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Similar Datasets

Discrete kinetic models from funneled energy landscape simulations.

Project description:A general method for facilitating the interpretation of computer simulations of protein folding with minimally frustrated energy landscapes is detailed and applied to a designed ankyrin repeat protein (4ANK). In the method, groups of residues are assigned to foldons and these foldons are used to map the conformational space of the protein onto a set of discrete macrobasins. The free energies of the individual macrobasins are then calculated, informing practical kinetic analysis. Two simple assumptions about the universality of the rate for downhill transitions between macrobasins and the natural local connectivity between macrobasins lead to a scheme for predicting overall folding and unfolding rates, generating chevron plots under varying thermodynamic conditions, and inferring dominant kinetic folding pathways. To illustrate the approach, free energies of macrobasins were calculated from biased simulations of a non-additive structure-based model using two structurally motivated foldon definitions at the full and half ankyrin repeat resolutions. The calculated chevrons have features consistent with those measured in stopped flow chemical denaturation experiments. The dominant inferred folding pathway has an "inside-out", nucleation-propagation like character.

| S-EPMC3520928 | biostudies-literature

How well does a funneled energy landscape capture the folding mechanism of spectrin domains?

Project description:Three structurally similar domains from α-spectrin have been shown to fold very differently. First, there is a contrast in the folding mechanism, as probed by Φ-value analysis, between the R15 domain and the R16 and R17 domains. Second, there are very different contributions from internal friction to folding: the folding rate of the R15 domain was found to be inversely proportional to solvent viscosity, showing no apparent frictional contribution from the protein, but in the other two domains, a large internal friction component was evident. Non-native misdocking of helices has been suggested to be responsible for this phenomenon. Here, I study the folding of these three proteins with minimalist coarse-grained models based on a funneled energy landscape. Remarkably, I find that, despite the absence of non-native interactions, the differences in folding mechanism of the domains are well captured by the model, and the agreement of the Φ-values with experiment is fairly good. On the other hand, within the context of this model, there are no significant differences in diffusion coefficient along the chosen folding coordinate, and the model cannot explain the large differences in folding rates between the proteins found experimentally. These results are nonetheless consistent with the expectations from the energy landscape perspective of protein folding, namely, that the folding mechanism is primarily determined by the native-like interactions present in the Gō-like model, with missing non-native interactions being required to explain the differences in "internal friction" seen in experiment.

| S-EPMC3808457 | biostudies-literature

Equilibrium and kinetic folding pathways of a TIM barrel with a funneled energy landscape.

Project description:The role of native contact topology in the folding of a TIM barrel model based on the alpha-subunit of tryptophan synthase (alphaTS) from Salmonella typhimurium (Protein Data Bank structure 1BKS) was studied using both equilibrium and kinetic simulations. Equilibrium simulations of alphaTS reveal the population of two intermediate ensembles, I1 and I2, during unfolding/refolding at the folding temperature, Tf = 335 K. Equilibrium intermediate I1 demonstrates discrete structure in regions alpha0-beta6 whereas intermediate I2 is a loose ensemble of states with N-terminal structure varying from at least beta1-beta3 (denoted I2A) to alpha0-beta4 at most (denoted I2B). The structures of I1 and I2 match well with the two intermediate states detected in equilibrium folding experiments of Escherichia coli alphaTS. Kinetic folding simulations of alphaTS reveal the sequential population of four intermediate ensembles, I120Q, I200Q, I300Q, and I360Q, during refolding. Kinetic intermediates I120Q, I200Q, and I300Q are highly similar to equilibrium alphaTS intermediates I2A, I2B, and I1, respectively, consistent with kinetic experiments on alphaTS from E. coli. A small population (approximately 10%) of kinetic trajectories are trapped in the I120Q intermediate ensemble and require a slow and complete unfolding step to properly refold. Both the on-pathway and off-pathway I120Q intermediates show structure in beta1-beta3, which is also strikingly consistent with kinetic folding experiments of alphaTS. In the off-pathway intermediate I(120Q), helix alpha2 is wrapped in a nonnative chiral arrangement around strand beta3, sterically preventing the subsequent folding step between beta3 and beta4. These results demonstrate the success of combining kinetic and equilibrium simulations of minimalist protein models to explore TIM barrel folding and the folding of other large proteins.

| S-EPMC1366549 | biostudies-literature

On the importance of a funneled energy landscape for the assembly and regulation of multidomain Src tyrosine kinases.

Project description:Regulation of signaling pathways in the cell often involves multidomain allosteric enzymes that are able to adopt alternate active or inactive conformations in response to specific stimuli. It is therefore of great interest to elucidate the energetic and structural determinants that govern the conformational plasticity of these proteins. In this study, free-energy computations have been used to address this fundamental question, focusing on one important family of signaling enzymes, the Src tyrosine kinases. Inactivation of these enzymes depends on the formation of an assembly comprising a tandem of SH3 and SH2 modules alongside a catalytic domain. Activation results from the release of the SH3 and SH2 domains, which are then believed to be structurally uncoupled by virtue of a flexible peptide link. In contrast to this view, this analysis shows that inactivation depends critically on the intrinsic propensity of the SH3-SH2 tandem to adopt conformations that are conducive to the assembled inactive state, even when no interactions with the rest of the kinase are possible. This funneling of the available conformational space is encoded within the SH3-SH2 connector, which appears to have evolved to modulate the flexibility of the tandem in solution. To further substantiate this notion, we show how constitutively activating mutations in the SH3-SH2 connector shift the assembly equilibrium toward the disassembled, active state. Based on a similar analysis of several constructs of the kinase complex, we propose that assembly is characterized by the progressive optimization of the protein's conformational energy, with little or no energetic frustration.

| S-EPMC1959435 | biostudies-literature

Effects of Ligand Binding on the Energy Landscape of Acyl-CoA-Binding Protein.

Project description:Binding of ligands is often crucial for function yet the effects of ligand binding on the mechanical stability and energy landscape of proteins are incompletely understood. Here, we use a combination of single-molecule optical tweezers and MD simulations to investigate the effect of ligand binding on the energy landscape of acyl-coenzyme A (CoA)-binding protein (ACBP). ACBP is a topologically simple and highly conserved four-α-helix bundle protein that acts as an intracellular transporter and buffer for fatty-acyl-CoA and is active in membrane assembly. We have previously described the behavior of ACBP under tension, revealing a highly extended transition state (TS) located almost halfway between the unfolded and native states. Here, we performed force-ramp and force-jump experiments, in combination with advanced statistical analysis, to show that octanoyl-CoA binding increases the activation free energy for the unfolding reaction of ACBP without affecting the position of the transition state along the reaction coordinate. It follows that ligand binding enhances the mechanical resistance and thermodynamic stability of the protein, without changing its mechanical compliance. Steered molecular dynamics simulations allowed us to rationalize the results in terms of key interactions that octanoyl-CoA establishes with the four α-helices of ACBP and showed that the unfolding pathway is marginally affected by the ligand. The results show that ligand-induced mechanical stabilization effects can be complex and may prove useful for the rational design of stabilizing ligands.

| S-EPMC7677128 | biostudies-literature

Imaging G protein-coupled receptors while quantifying their ligand-binding free-energy landscape.

Project description:Imaging native membrane receptors and testing how they interact with ligands is of fundamental interest in the life sciences but has proven remarkably difficult to accomplish. Here, we introduce an approach that uses force-distance curve-based atomic force microscopy to simultaneously image single native G protein-coupled receptors in membranes and quantify their dynamic binding strength to native and synthetic ligands. We measured kinetic and thermodynamic parameters for individual protease-activated receptor-1 (PAR1) molecules in the absence and presence of antagonists, and these measurements enabled us to describe PAR1's ligand-binding free-energy landscape with high accuracy. Our nanoscopic method opens an avenue to directly image and characterize ligand binding of native membrane receptors.

| S-EPMC5087271 | biostudies-literature

Prediction of native-state hydrogen exchange from perfectly funneled energy landscapes.

Project description:Simulations based on perfectly funneled energy landscapes often capture many of the kinetic features of protein folding. We examined whether simulations based on funneled energy functions can also describe fluctuations in native-state protein ensembles. We quantitatively compared the site-specific local stability determined from structure-based folding simulations, with hydrogen exchange protection factors measured experimentally for ubiquitin, chymotrypsin inhibitor 2, and staphylococcal nuclease. Different structural definitions for the open and closed states based on the number of native contacts for each residue, as well as the hydrogen-bonding state, or a combination of both criteria were evaluated. The predicted exchange patterns agree with the experiments under native conditions, indicating that protein topology indeed has a dominant effect on the exchange kinetics. Insights into the simplest mechanistic interpretation of the amide exchange process were thus obtained.

| S-EPMC3203634 | biostudies-literature

Protein-Ligand Binding Free-Energy Calculations with ARROW─A Purely First-Principles Parameterized Polarizable Force Field.

Project description:Protein-ligand binding free-energy calculations using molecular dynamics (MD) simulations have emerged as a powerful tool for in silico drug design. Here, we present results obtained with the ARROW force field (FF)─a multipolar polarizable and physics-based model with all parameters fitted entirely to high-level ab initio quantum mechanical (QM) calculations. ARROW has already proven its ability to determine solvation free energy of arbitrary neutral compounds with unprecedented accuracy. The ARROW FF parameterization is now extended to include coverage of all amino acids including charged groups, allowing molecular simulations of a series of protein-ligand systems and prediction of their relative binding free energies. We ensure adequate sampling by applying a novel technique that is based on coupling the Hamiltonian Replica exchange (HREX) with a conformation reservoir generated via potential softening and nonequilibrium MD. ARROW provides predictions with near chemical accuracy (mean absolute error of ∼0.5 kcal/mol) for two of the three protein systems studied here (MCL1 and Thrombin). The third protein system (CDK2) reveals the difficulty in accurately describing dimer interaction energies involving polar and charged species. Overall, for all of the three protein systems studied here, ARROW FF predicts relative binding free energies of ligands with a similar accuracy level as leading nonpolarizable force fields.

| S-EPMC9753910 | biostudies-literature

Mapping the energy landscape of protein-ligand binding <i>via</i> linear free energy relationships determined by protein NMR relaxation dispersion.

Project description:Biochemical signaling is mediated by complexes between macromolecular receptors and their ligands, with the duration of the signal being directly related to the lifetime of the ligand-receptor complex. In the field of drug design, the recognition that drug efficacy in vivo depends on the lifetime of the drug-protein complex has spawned the concept of designing drugs with particular binding kinetics. To advance this field it is critical to investigate how the molecular details of designed ligands might affect the binding kinetics, as well as the equilibrium binding constant. Here we use protein NMR relaxation dispersion to determine linear free energy relationships involving the on- and off-rates and the affinity for a series of congeneric ligands targeting the carbohydrate recognition domain of galectin-3. Using this approach we determine the energy landscape and the position of the transition state along the reaction coordinate of protein-ligand binding. The results show that ligands exhibiting reduced off-rates achieve this by primarily stabilizing the bound state, but do not affect the transition state to any greater extent. The transition state forms early, that is, it is located significantly closer to the free state than to the bound state, suggesting a critical role of desolvation. Furthermore, the data suggest that different subclasses of ligands show different behavior with respect to these characteristics.

| S-EPMC8341105 | biostudies-literature

Role of Electrostatics in Protein-RNA Binding: The Global vs the Local Energy Landscape.

Project description:U1A protein-stem loop 2 RNA association is a basic step in the assembly of the spliceosomal U1 small nuclear ribonucleoprotein. Long-range electrostatic interactions due to the positive charge of U1A are thought to provide high binding affinity for the negatively charged RNA. Short range interactions, such as hydrogen bonds and contacts between RNA bases and protein side chains, favor a specific binding site. Here, we propose that electrostatic interactions are as important as local contacts in biasing the protein-RNA energy landscape toward a specific binding site. We show by using molecular dynamics simulations that deletion of two long-range electrostatic interactions (K22Q and K50Q) leads to mutant-specific alternative RNA bound states. One of these states preserves short-range interactions with aromatic residues in the original binding site, while the other one does not. We test the computational prediction with experimental temperature-jump kinetics using a tryptophan probe in the U1A-RNA binding site. The two mutants show the distinct predicted kinetic behaviors. Thus, the stem loop 2 RNA has multiple binding sites on a rough RNA-protein binding landscape. We speculate that the rough protein-RNA binding landscape, when biased to different local minima by electrostatics, could be one way that protein-RNA interactions evolve toward new binding sites and novel function.

| S-EPMC6044910 | biostudies-literature

OmicsDI is part of the ELIXIR infrastructure

OmicsDI is an Elixir interoperability service. Learn more ›

Tweets

OmicsDI Databases

PRIDE
PeptideAtlas
MassIVE
JPOST Repository
Physiome Model Repository

EGA
EVA
ENA
LINCS
PAXDB
Cell Collective

MetaboLights
Metabolomics Workbench
MetabolomeExpress
GNPS
BioModels
FAIRDOMHub

ArrayExpress
dbGaP
ExpressionAtlas
GEO
NODE

Information

Databases
Help
API
Contact us
Code on GitHub
Terms of use
Submit Data