Browse
Submit Data
Databases
API
Help

Dataset Information

0 Views

0 Connections

0 Citations

0 Reanalyses

0 Downloads

Omics score: 0

Determination of the interfacial water content in protein-protein complexes from free energy simulations.

ABSTRACT: The question as to how many tightly or weakly bound water molecules are located in interfaces between protein-protein complex constituents is addressed from a phase equilibrium point of view by developing a theory in the canonical ensemble. A fast method based on free energy simulations is described for computing the number of water molecules in the interface regions. Results are given for 211 interfacial cavities of 26 antigen-antibody complexes for which experimentally determined structures are found in the Protein Data Bank. The accuracy of the method is assessed and the computational water content is compared with experimental data, revealing the amount of water molecules not resolved by experimental approaches.

SUBMITTER: Monecke P

PROVIDER: S-EPMC1367109 | biostudies-literature | 2006 Feb

REPOSITORIES: biostudies-literature

ACCESS DATA

Json Xml

Publications

Determination of the interfacial water content in protein-protein complexes from free energy simulations.

Monecke Peter P Borosch Thorsten T Brickmann Jürgen J Kast Stefan M SM

Biophysical journal 20051111 3

The question as to how many tightly or weakly bound water molecules are located in interfaces between protein-protein complex constituents is addressed from a phase equilibrium point of view by developing a theory in the canonical ensemble. A fast method based on free energy simulations is described for computing the number of water molecules in the interface regions. Results are given for 211 interfacial cavities of 26 antigen-antibody complexes for which experimentally determined structures ar ...[more]

PMID: 16284258

Similar Datasets

Accounting for the Central Role of Interfacial Water in Protein-Ligand Binding Free Energy Calculations.

Project description:Rigorous binding free energy methods in drug discovery are growing in popularity because of a combination of methodological advances, improvements in computer hardware, and workflow automation. These calculations typically use molecular dynamics (MD) to sample from the Boltzmann distribution of conformational states. However, when part or all of the binding sites is inaccessible to the bulk solvent, the time needed for water molecules to equilibrate between bulk solvent and the binding site can be well beyond what is practical with standard MD. This sampling limitation is problematic in relative binding free energy calculations, which compute the reversible work of converting ligand 1 to ligand 2 within the binding site. Thus, if ligand 1 is smaller and/or more polar than ligand 2, the perturbation may allow additional water molecules to occupy a region of the binding site. However, this change in hydration may not be captured by standard MD simulations and may therefore lead to errors in the computed free energy. We recently developed a hybrid Monte Carlo/MD (MC/MD) method, which speeds up the equilibration of water between bulk solvent and buried cavities, while sampling from the intended distribution of states. Here, we report on the use of this approach in the context of alchemical binding free energy calculations. We find that using MC/MD markedly improves the accuracy of the calculations and also reduces hysteresis between the forward and reverse perturbations, relative to matched calculations using only MD with or without the crystallographic water molecules. The present method is available for use in AMBER simulation software.

| S-EPMC7725968 | biostudies-literature

Protein thermostability calculations using alchemical free energy simulations.

Project description:Thermal stability of proteins is crucial for both biotechnological and therapeutic applications. Rational protein engineering therefore frequently aims at increasing thermal stability by introducing stabilizing mutations. The accurate prediction of the thermodynamic consequences caused by mutations, however, is highly challenging as thermal stability changes are caused by alterations in the free energy of folding. Growing computational power, however, increasingly allows us to use alchemical free energy simulations, such as free energy perturbation or thermodynamic integration, to calculate free energy differences with relatively high accuracy. In this article, we present an automated protocol for setting up alchemical free energy calculations for mutations of naturally occurring amino acids (except for proline) that allows an unprecedented, automated screening of large mutant libraries. To validate the developed protocol, we calculated thermodynamic stability differences for 109 mutations in the microbial Ribonuclease Barnase. The obtained quantitative agreement with experimental data illustrates the potential of the approach in protein engineering and design.

| S-EPMC2872215 | biostudies-literature

Free-energy landscape of protein oligomerization from atomistic simulations.

Project description:In the realm of protein-protein interactions, the assembly process of homooligomers plays a fundamental role because the majority of proteins fall into this category. A comprehensive understanding of this multistep process requires the characterization of the driving molecular interactions and the transient intermediate species. The latter are often short-lived and thus remain elusive to most experimental investigations. Molecular simulations provide a unique tool to shed light onto these complex processes complementing experimental data. Here we combine advanced sampling techniques, such as metadynamics and parallel tempering, to characterize the oligomerization landscape of fibritin foldon domain. This system is an evolutionarily optimized trimerization motif that represents an ideal model for experimental and computational mechanistic studies. Our results are fully consistent with previous experimental nuclear magnetic resonance and kinetic data, but they provide a unique insight into fibritin foldon assembly. In particular, our simulations unveil the role of nonspecific interactions and suggest that an interplay between thermodynamic bias toward native structure and residual conformational disorder may provide a kinetic advantage.

| S-EPMC3856838 | biostudies-literature

Blind protein structure prediction using accelerated free-energy simulations.

Project description:We report a key proof of principle of a new acceleration method [Modeling Employing Limited Data (MELD)] for predicting protein structures by molecular dynamics simulation. It shows that such Boltzmann-satisfying techniques are now sufficiently fast and accurate to predict native protein structures in a limited test within the Critical Assessment of Structure Prediction (CASP) community-wide blind competition.

| S-EPMC5106196 | biostudies-literature

Relative energies of binding for antibody-carbohydrate-antigen complexes computed from free-energy simulations.

Project description:Free-energy perturbation (FEP) simulations have been applied to a series of analogues of the natural trisaccharide epitope of Salmonella serotype B bound to a fragment of the monoclonal anti-Salmonella antibody Se155-4. This system was selected in order to assess the ability of free-energy perturbation (FEP) simulations to predict carbohydrate-protein interaction energies. The ultimate goal is to use FEP simulations to aid in the design of synthetic high affinity ligands for carbohydrate-binding proteins. The molecular dynamics (MD) simulations were performed in the explicit presence of water molecules, at room temperature. The AMBER force field, with the GLYCAM parameter set for oligosaccharides, was employed. In contrast to many modeling protocols, FEP simulations are capable of including the effects of entropy, arising from differential ligand flexibilities and solvation properties. The experimental binding affinities are all close in value, resulting in small relative free energies of binding. Many of the DeltaDeltaG values are on the order of 0-1 kcal mol(-1), making their accurate calculation particularly challenging. The simulations were shown to reasonably reproduce the known geometries of the ligands and the ligand-protein complexes. A model for the conformational behavior of the unbound antigen is proposed that is consistent with the reported NMR data. The best agreement with experiment was obtained when histidine 97H was treated as fully protonated, for which the relative binding energies were predicted to well within 1 kcal mol(-1). To our knowledge this is the first report of FEP simulations applied to an oligosaccharide-protein complex.

| S-EPMC1764640 | biostudies-literature

Mechanism-Dependent Modulation of Ultrafast Interfacial Water Dynamics in Intrinsically Disordered Protein Complexes.

Project description:The recognition of intrinsically disordered proteins (IDPs) is highly dependent on dynamics owing to the lack of structure. Here we studied the interplay between dynamics and molecular recognition in IDPs with a combination of time-resolving tools on timescales ranging from femtoseconds to nanoseconds. We interrogated conformational dynamics and surface water dynamics and its attenuation upon partner binding using two IDPs, IBB and Nup153FG, both of central relevance to the nucleocytoplasmic transport machinery. These proteins bind the same nuclear transport receptor (Importinβ) with drastically different binding mechanisms, coupled folding-binding and fuzzy complex formation, respectively. Solvent fluctuations in the dynamic interface of the Nup153FG-Importinβ fuzzy complex were largely unperturbed and slightly accelerated relative to the unbound state. In the IBB-Importinβ complex, on the other hand, substantial relative slowdown of water dynamics was seen in a more rigid interface. These results show a correlation between interfacial water dynamics and the plasticity of IDP complexes, implicating functional relevance for such differential modulation in cellular processes, including nuclear transport.

| S-EPMC6563697 | 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

Foaming properties and air-water interfacial behavior of corn protein hydrolyzate-tannic acid complexes.

Project description:The complexation of corn protein hydrolyzate (CPH) with tannic acid (TA) was utilized to improve the foaming properties of CPH itself, and the air-water interfacial behavior of CPH-TA complex was also investigated. The results showed that the surface hydrophobicity of pure CPH was significantly decreased in bulk solution after the complexation with TA. Compared with pure CPH, the foams stabilized by CPH-TA complex showed higher interfacial thickness between the bubbles, which well explained the better long term stability of the corresponding foams. Therefore, the complexation maintained the good foaming capacity of CPH itself, but considerably increased its foam stability. Moreover, the air-water interfacial behavior study demonstrated that the complexation slightly decreased the interfacial activity of CPH itself, but considerably increased its interfacial viscoelasticity, suggesting more stable of the air-water interface stabilized by CPH-TA complex compared with that stabilized by CPH alone. These findings indicated that foaming properties of the surface active components were closely related with its air-water interfacial behavior. The study suggested that CPH-TA complex could be used as a stabilizer in constructing the peptides-based foams.

| S-EPMC6400733 | biostudies-literature

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

Extensive free-energy simulations identify water as the base in nucleotide addition by DNA polymerase.

Project description:Transphosphorylation of nucleotide triphosphates is the central reaction in DNA replication by DNA polymerase as well as many other biological processes. Despite its importance, the microscopic chemical mechanism of transphosphorylation of nucleotide triphosphates is, in most cases, unknown. Here we use extensive simulations of DNA polymerase ? to test mechanistic hypotheses. We systematically survey the reactive space by calculating 2D free-energy surfaces for 10 different plausible mechanisms that have been proposed. We supplement these free-energy surfaces with calculations of pKa for a number of potentially acidic protons in different states relevant to the catalytic cycle. We find that among all of the conditions that we test, the smallest activation barrier occurs for a reaction where a Mg2+-coordinated water deprotonates the nucleophilic 3'-OH, and this deprotonation is concerted with the phosphoryl transfer. The presence of a third Mg2+ in the active site lowers the activation barrier for the water-as-base mechanism, as does protonation of the pyrophosphate leaving group, which is consistent with general acid catalysis. The results demonstrate the value of simulations, when used in conjunction with experimental data, to help establish a microscopic chemical mechanism in a complex environment.

| S-EPMC6911213 | 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