Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Calculation of the free energy of association for protein complexes.

ABSTRACT: We have developed a method for calculating the association energy of quaternary complexes starting from their atomic coordinates. The association energy is described as the sum of two solvation terms and an energy term to account for the loss of translational and rotational entropy. The calculated solvation energy, using atomic solvation parameters and the solvent accessible surface areas, has a correlation of 96% with experimentally determined values. We have applied this methodology to examine intermediates in viral assembly and to assess the contribution isomerization makes to the association energy of molecular complexes. In addition, we have shown that the calculated association can be used as a predictive tool for analyzing modeled molecular complexes.

SUBMITTER: Horton N

PROVIDER: S-EPMC2142085 | biostudies-other | 1992 Jan

REPOSITORIES: biostudies-other

ACCESS DATA

Json Xml

Publications

Calculation of the free energy of association for protein complexes.

Horton N N Lewis M M

Protein science : a publication of the Protein Society 19920101 1

We have developed a method for calculating the association energy of quaternary complexes starting from their atomic coordinates. The association energy is described as the sum of two solvation terms and an energy term to account for the loss of translational and rotational entropy. The calculated solvation energy, using atomic solvation parameters and the solvent accessible surface areas, has a correlation of 96% with experimentally determined values. We have applied this methodology to examine ...[more]

PMID: 1339024

Similar Datasets

Free energy landscapes of encounter complexes in protein-protein association.

Project description:We report the computer generation of a high-density map of the thermodynamic properties of the diffusion-accessible encounter conformations of four receptor-ligand protein pairs, and use it to study the electrostatic and desolvation components of the free energy of association. Encounter complex conformations are generated by sampling the translational/rotational space of the ligand around the receptor, both at 5-A and zero surface-to-surface separations. We find that partial desolvation is always an important effect, and it becomes dominant for complexes in which one of the reactants is neutral or weakly charged. The interaction provides a slowly varying attractive force over a small but significant region of the molecular surface. In complexes with no strong charge complementarity this region surrounds the binding site, and the orientation of the ligand in the encounter conformation with the lowest desolvation free energy is similar to the one observed in the fully formed complex. Complexes with strong opposite charges exhibit two types of behavior. In the first group, represented by barnase/barstar, electrostatics exerts strong orientational steering toward the binding site, and desolvation provides some added adhesion within the local region of low electrostatic energy. In the second group, represented by the complex of kallikrein and pancreatic trypsin inhibitor, the overall stability results from the rather nonspecific electrostatic attraction, whereas the affinity toward the binding region is determined by desolvation interactions.

| S-EPMC1300098 | biostudies-other

Calculation of the free energy and cooperativity of protein folding.

Project description:Calculation of the free energy of protein folding and delineation of its pre-organization are of foremost importance for understanding, predicting and designing biological macromolecules. Here, we introduce an energy smoothing variant of parallel tempering replica exchange Monte Carlo (REMS) that allows for efficient configurational sampling of flexible solutes under the conditions of molecular hydration. Its usage to calculate the thermal stability of a model globular protein, Trp cage TC5b, achieves excellent agreement with experimental measurements. We find that the stability of TC5b is attained through the coupled formation of local and non-local interactions. Remarkably, many of these structures persist at high temperature, concomitant with the origin of native-like configurations and mesostates in an otherwise macroscopically disordered unfolded state. Graph manifold learning reveals that the conversion of these mesostates to the native state is structurally heterogeneous, and that the cooperativity of their formation is encoded largely by the unfolded state ensemble. In all, these studies establish the extent of thermodynamic and structural pre-organization of folding of this model globular protein, and achieve the calculation of macromolecular stability ab initio, as required for ab initio structure prediction, genome annotation, and drug design.

| S-EPMC1865387 | biostudies-literature

Prediction of protein solubility from calculation of transfer free energy.

Project description:Solubility plays a major role in protein purification, and has serious implications in many diseases. We studied the effects of pH and mutations on protein solubility by calculating the transfer free energy from the condensed phase to the solution phase. The condensed phase was modeled as an implicit solvent, with a dielectric constant lower than that of water. To account for the effects of pH, the protonation states of titratable side chains were sampled by running constant-pH molecular dynamics simulations. Conformations were then selected for calculations of the electrostatic solvation energy: once for the condensed phase, and once for the solution phase. The average transfer free energy from the condensed phase to the solution phase was found to predict reasonably well the variations in solubility of ribonuclease Sa and insulin with pH. This treatment of electrostatic contributions combined with a similar approach for nonelectrostatic contributions led to a quantitative rationalization of the effects of point mutations on the solubility of ribonuclease Sa. This study provides valuable insights into the physical basis of protein solubility.

| S-EPMC2527251 | biostudies-literature

Revisiting free energy calculations: a theoretical connection to MM/PBSA and direct calculation of the association free energy.

Project description:The prediction of absolute ligand-receptor binding affinities is essential in a wide range of biophysical queries, from the study of protein-protein interactions to structure-based drug design. End-point free energy methods, such as the Molecular Mechanics Poisson-Boltzmann Surface Area (MM/PBSA) model, have received much attention and widespread application in recent literature. These methods benefit from computational efficiency as only the initial and final states of the system are evaluated, yet there remains a need for strengthening their theoretical foundation. Here a clear connection between statistical thermodynamics and end-point free energy models is presented. The importance of the association free energy, arising from one molecule's loss of translational and rotational freedom from the standard state concentration, is addressed. A novel method for calculating this quantity directly from a molecular dynamics simulation is described. The challenges of accounting for changes in the protein conformation and its fluctuations from separate simulations are discussed. A simple first-order approximation of the configuration integral is presented to lay the groundwork for future efforts. This model has been applied to FKBP12, a small immunophilin that has been widely studied in the drug industry for its potential immunosuppressive and neuroregenerative effects.

| S-EPMC1303837 | biostudies-literature

Calculation of absolute protein-ligand binding free energy from computer simulations.

Project description:A general methodology for calculating the equilibrium binding constant of a flexible ligand to a protein receptor is formulated on the basis of potentials of mean force. The overall process is decomposed into several stages that can be computed separately: the free ligand in the bulk is first restrained into the conformation it adopts in the bound state, position, and orientation by applying biasing potentials, then it is translated into the binding site, where it is released completely. The conformational restraining potential is based on the root-mean-square deviation of the peptide coordinates relative to its average conformation in the bound complex. Free energy contributions from each stage are calculated by means of free energy perturbation potential of mean force techniques by using appropriate order parameters. The present approach avoids the need to decouple the ligand from its surrounding (bulk solvent and receptor protein) as is traditionally performed in the double-decoupling scheme. It is believed that the present formulation will be particularly useful when the solvation free energy of the ligand is very large. As an application, the equilibrium binding constant of the phosphotyrosine peptide pYEEI to the Src homology 2 domain of human Lck has been calculated. The results are in good agreement with experimental values.

| S-EPMC1100764 | biostudies-literature

Calculation of mutational free energy changes in transition states for protein folding.

Project description:Recent advances in experimental and computational methods have made it possible to determine with considerable accuracy the structures whose formation is rate limiting for the folding of some small proteins-the transition state ensemble, or TSE. We present a method to analyze and validate all-atom models of such structures. The method is based on the comparison of experimental data with the computation of the change in free energy of the TSE resulting from specific mutations. Each mutation is modeled individually in all members of an ensemble of transition state structures using a method originally developed to predict mutational changes in the stability of native proteins. We first apply this method to six proteins for which we have determined the TSEs with a technique that uses experimental mutational data (Phi-values) as restraints in the structure determination and find a highly significant correlation between the calculated free energy changes and those derived from experimental kinetic data. We then use the procedure to analyze transition state structures determined by molecular dynamics simulations of unfolding, again finding a high correlation. Finally, we use the method to estimate changes in folding rates of several hydrophobic core mutants of Fyn SH3. Taken together, these results show that the procedure developed here is a tool of general validity for analyzing, assessing, and improving the quality of the structures of transition states for protein folding.

| S-EPMC1303238 | biostudies-literature

Saturation Mutagenesis by Efficient Free-Energy Calculation.

Project description:Single-point mutations in proteins can greatly influence protein stability, binding affinity, protein function or its expression per se. Here, we present accurate and efficient predictions of the free energy of mutation of amino acids. We divided the complete mutational free energy into an uncharging step, which we approximate by a third-power fitting (TPF) approach, and an annihilation step, which we approximate using the one-step perturbation (OSP) method. As a diverse set of test systems, we computed the solvation free energy of all amino acid side chain analogues and obtained an excellent agreement with thermodynamic integration (TI) data. Moreover, we calculated mutational free energies in model tripeptides and established an efficient protocol involving a single reference state. Again, the approximate methods agreed excellently with the TI references, with a root-mean-square error of only 3.6 kJ/mol over 17 mutations. Our combined TPF+OSP approach does show not only a very good agreement but also a 2-fold higher efficiency than full blown TI calculations.

| S-EPMC5813279 | biostudies-literature

Trypsin-ligand binding free energy calculation with AMOEBA.

Project description:The binding free energies of several benzamidine -like inhibitors to trypsin were examined using a polarizable molecular mechanics potential. All the computed binding free energies are in good agreement with the experimental data. From free energy decomposition, electrostatic interaction was indicated to be the driving force for the binding. MD simulations show that the ligands form hydrogen bonds with trypsin and water molecules nearby in a competitive fashion. While the binding free energy is independent of the ligand dipole moment, it shows a strong correlation with the ligand molecular polarizability.

| S-EPMC2819397 | biostudies-literature

Calculation of protein-ligand binding free energy by using a polarizable potential.

Project description:The binding of charged ligands benzamidine and diazamidine to trypsin was investigated by using a polarizable potential energy function and explicit-water molecular dynamics simulations. The binding free energies were computed from the difference between the free energies of decoupling the ligand from water and protein environments. Both the absolute and the relative free energies from the perturbation simulations agree with experimental measurements to within 0.5 kcal.mol(-1). Comparison of free-energy components sampled from different thermodynamic paths indicates that electrostatics is the main driving force behind benzamidine recognition of trypsin. The contribution of electronic polarization to binding appears to be crucial. By computing the free-energy contribution caused by the polarization between the ligand and its surroundings, we found that polarization has the opposite effect in dissimilar environments. Although polarization favors ligand solvation in water, it weakens the protein-ligand attraction by screening the electrostatic interaction between trypsin and benzamidine. We also examined the relative binding free energies of a benzamidine analog diazamidine to trypsin. The changes in free energy on benzamidine-diazamidine substitution were tens of kilocalories in both water and trypsin environments; however, the change in the total binding free energy is <2 kcal.mol(-1) because of cancellation, consistent with the experimental results. Overall, our results suggest that the use of a polarizable force field, given adequate sampling, is capable of achieving chemical accuracy in molecular simulations of protein-ligand recognition.

| S-EPMC2359813 | biostudies-other

Empirical free energy calculation: comparison to calorimetric data.

Project description:An effective free energy potential, developed originally for binding free energy calculation, is compared to calorimetric data on protein unfolding, described by a linear combination of changes in polar and nonpolar surface areas. The potential consists of a molecular mechanics energy term calculated for a reference medium (vapor or nonpolar liquid), and empirical terms representing solvation and entropic effects. It is shown that, under suitable conditions, the free energy function agrees well with the calorimetric expression. An additional result of the comparison is an independent estimate of the side-chain entropy loss, which is shown to agree with a structure-based entropy scale. These findings confirm that simple functions can be used to estimate the free energy change in complex systems, and that a binding free energy evaluation model can describe the thermodynamics of protein unfolding correctly. Furthermore, it is shown that folding and binding leave the sum of solute-solute and solute-solvent van der Waals interactions nearly invariant and, due to this invariance, it may be advantageous to use a nonpolar liquid rather than vacuum as the reference medium.

| S-EPMC2143797 | biostudies-other

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