Conformational contribution to thermodynamics of binding in protein-peptide complexes through microscopic simulation.
Ontology highlight
ABSTRACT:
SUBMITTER: Das A
PROVIDER: S-EPMC3602779 | biostudies-literature | 2013 Mar
REPOSITORIES: biostudies-literature
Publications
Conformational contribution to thermodynamics of binding in protein-peptide complexes through microscopic simulation.
Das Amit A Chakrabarti J J Ghosh Mahua M
Biophysical journal 20130319 6
We extract the thermodynamics of conformational changes in biomacromolecular complexes from the distributions of the dihedral angles of the macromolecules. These distributions are obtained from the equilibrium configurations generated via all-atom molecular dynamics simulations. The conformational thermodynamics data we obtained for calmodulin-peptide complexes using our methodology corroborate well with the experimentally observed conformational and binding entropies. The conformational free-en ...[more]
Similar Datasets
Project description:Structure-based drug design has often been restricted by the rather static picture of protein-ligand complexes presented by crystal structures, despite the widely accepted importance of protein flexibility in biomolecular recognition. Here we report a detailed experimental and computational study of the drug target, human heat shock protein 90, to explore the contribution of protein dynamics to the binding thermodynamics and kinetics of drug-like compounds. We observe that their binding properties depend on whether the protein has a loop or a helical conformation in the binding site of the ligand-bound state. Compounds bound to the helical conformation display slow association and dissociation rates, high-affinity and high cellular efficacy, and predominantly entropically driven binding. An important entropic contribution comes from the greater flexibility of the helical relative to the loop conformation in the ligand-bound state. This unusual mechanism suggests increasing target flexibility in the bound state by ligand design as a new strategy for drug discovery.
| S-EPMC5741624 | biostudies-literature
Project description:The oligopeptide-binding protein OppA provides a useful model system for studying the physical chemistry underlying noncovalent interactions since it binds a variety of readily synthesized ligands. We have studied the binding of eight closely related tripeptides of the type Lysine-X-Lysine, where X is an abnormal amino acid, by isothermal titration calorimetry (ITC) and X-ray crystallography. The tripeptides fall into three series of ligands, which have been designed to examine the effects of small changes to the central side chain. Three ligands have a primary amine as the second side chain, two have a straight alkane chain, and three have ring systems. The results have revealed a definite preference for the binding of hydrophobic residues over the positively charged side chains, the latter binding only weakly due to unfavorable enthalpic effects. Within the series of positively charged groups, a point of lowest affinity has been identified and this is proposed to arise from unfavorable electrostatic interactions in the pocket, including the disruption of a key salt bridge. Marked entropy-enthalpy compensation is found across the series, and some of the difficulties in designing tightly binding ligands have been highlighted.
| S-EPMC2144387 | biostudies-other
Project description:The behaviour of biomolecular systems is governed by their thermodynamic and kinetic properties. It is thus important to be able to calculate, for example, both the affinity and rate of binding and dissociation of a protein-ligand complex, or the populations and exchange rates between distinct conformational states. Because these are typically rare events, calculating these properties from long molecular dynamics simulations remains extremely difficult. Instead, one often adopts a divide-and-conquer strategy in which equilibrium free-energy differences and the fastest state-to-state transition (e.g. ligand association or minor-to-major state conversion) are combined to estimate the slow rate (e.g. ligand dissociation) using a two-state assumption. Here we instead address these problems by using a previously developed method to calculate both the forward and backward rates directly from simulations. We then estimate the thermodynamics from the rates, and validate these values by independent means. We applied the approach to three systems of increasing complexity, including the association and dissociation of benzene to a fully buried cavity inside the L99A mutant variant of T4 lysozyme. In particular, we were able to determine both millisecond association and dissociation rates, and the affinity, of the protein-ligand system by directly observing dozens of rare events in atomic detail. Our approach both sheds light on the precision of methods for calculating kinetics and further provides a generally useful test for the internal consistency of kinetics and thermodynamics. We also expect our route to be useful for obtaining both the kinetics and thermodynamics at the same time in more challenging cases.
| S-EPMC5859887 | biostudies-literature
Project description:Advances in simulation techniques and computing hardware have created a substantial overlap between the timescales accessible to atomic-level simulations and those on which the fastest-folding proteins fold. Here we demonstrate, using simulations of four variants of the human villin headpiece, how simulations of spontaneous folding and unfolding can provide direct access to thermodynamic and kinetic quantities such as folding rates, free energies, folding enthalpies, heat capacities, ?-values, and temperature-jump relaxation profiles. The quantitative comparison of simulation results with various forms of experimental data probing different aspects of the folding process can facilitate robust assessment of the accuracy of the calculations while providing a detailed structural interpretation for the experimental observations. In the example studied here, the analysis of folding rates, ?-values, and folding pathways provides support for the notion that a norleucine double mutant of villin folds five times faster than the wild-type sequence, but following a slightly different pathway. This work showcases how computer simulation has now developed into a mature tool for the quantitative computational study of protein folding and dynamics that can provide a valuable complement to experimental techniques.
| S-EPMC3497772 | biostudies-literature
Project description:The binding of short disordered peptide stretches to globular protein domains is important for a wide range of cellular processes, including signal transduction, protein transport, and immune response. The often promiscuous nature of these interactions and the conformational flexibility of the peptide chain, sometimes even when bound, make the binding specificity of this type of protein interaction a challenge to understand. Here we develop and test a Monte Carlo-based procedure for calculating protein-peptide binding thermodynamics for many sequences in a single run. The method explores both peptide sequence and conformational space simultaneously by simulating a joint probability distribution which, in particular, makes searching through peptide sequence space computationally efficient. To test our method, we apply it to 3 different peptide-binding protein domains and test its ability to capture the experimentally determined specificity profiles. Insight into the molecular underpinnings of the observed specificities is obtained by analyzing the peptide conformational ensembles of a large number of binding-competent sequences. We also explore the possibility of using our method to discover new peptide-binding pockets on protein structures.
| S-EPMC3812049 | biostudies-literature
Project description:PixelDB, the Peptide Exosite Location Database, compiles 1966 non-redundant, high-resolution structures of protein-peptide complexes filtered to minimize the impact of crystal packing on peptide conformation. The database is organized to facilitate study of structurally conserved versus non-conserved elements of protein-peptide engagement. PixelDB clusters complexes based on the structural similarity of the peptide-binding protein, and by comparing complexes within a cluster highlights examples of domains that engage peptides using more than one binding mode. PixelDB also identifies conserved peptide core structural motifs characteristic of each binding mode. Peptide regions that flank core motifs often make non-structurally conserved interactions with the protein surface in regions we call exosites. Many examples establish that exosite contacts can be important for enhancing protein binding and interaction specificity. PixelDB provides a resource for computational and structural biologists to study, model, and predict core-motif and exosite-contacting peptide interactions. PixelDB is available to the community without restriction in a convenient flat-file format with accompanying visualization tools.
| S-EPMC5734312 | biostudies-literature
Project description:Conformational entropy is expected to contribute significantly to the thermodynamics of structural transitions in intrinsically disordered proteins or regions in response to protein/ligand binding, posttranslational modifications, and environmental changes. We calculated the backbone (dihedral) conformational entropy of oligoglycine (GlyN), a protein backbone mimic and model intrinsically disordered region, as a function of chain length (N=3, 4, 5, 10, and 15) from simulations using three different approaches. The backbone conformational entropy scales linearly with chain length with a slope consistent with the entropy of folding of well-structured proteins. The entropic contributions of second-order dihedral correlations are predominantly through intraresidue ?-? pairs, suggesting that oligoglycine may be thermodynamically modeled as a system of independent glycine residues. We find the backbone conformational entropy to be largely independent of global structural parameters, like the end-to-end distance and radius of gyration. We introduce a framework referred to herein as "ensemble confinement" to estimate the loss (gain) of conformational free energy and its entropic component when individual residues are constrained to (released from) particular regions of the ?-? map. Quantitatively, we show that our protein backbone model resists ordering/folding with a significant, unfavorable ensemble confinement free energy because of the loss of a substantial portion of the absolute backbone entropy. Proteins can couple this free-energy reservoir to distal binding events as a regulatory mechanism to promote or suppress binding.
| S-EPMC6026333 | biostudies-literature
Project description:S100P is a member of the S100 subfamily of calcium-binding proteins that are believed to be associated with various diseases, and in particular deregulation of S100P expression has been documented for prostate and breast cancer. Previously, we characterized the effects of metal binding on the conformational properties of S100P and proposed that S100P could function as a Ca2+ conformational switch. In this study we used fluorescence and CD spectroscopies and isothermal titration calorimetry to characterize the target-recognition properties of S100P using a model peptide, melittin. Based on these experimental data we show that S100P and melittin can interact in a Ca2+-dependent and -independent manner. Ca2+-independent binding occurs with low affinity (Kd approximately 0.2 mM), has a stoichiometry of four melittin molecules per S100P dimer and is presumably driven by favorable electrostatic interactions between the acidic protein and the basic peptide. In contrast, Ca2+-dependent binding of melittin to S100P occurs with high affinity (Kd approximately 5 microM) has a stoichiometry of two molecules of melittin per S100P dimer, appears to have positive cooperativity, and is driven by hydrophobic interactions. Furthermore, Ca2+-dependent S100P-melittin complex formation is accompanied by significant conformational changes: Melittin, otherwise unstructured in solution, adopts a helical conformation upon interaction with Ca2+-S100P. These results support a model for the Ca2+-dependent conformational switch in S100P for functional target recognition.
| S-EPMC2373636 | biostudies-literature
Project description:Peptide toxins isolated from venomous creatures, long prized as research tools due to their innate potency for ion channels, are emerging as drugs as well. However, it remains challenging to understand why peptide toxins bind with high potency to ion channels, to identify residues that are key for activity, and to improve their affinities via mutagenesis. We use WaterMap, a molecular dynamics simulation-based method, to gain computational insight into these three questions by calculating the locations and thermodynamic properties of water molecules in the peptide toxin binding sites of five ion channels. These include an acid-sensing ion channel, voltage-gated potassium channel, sodium channel in activated and deactivated states, transient-receptor potential channel, and a nicotinic receptor whose structures were recently determined by crystallography and cryo-electron microscopy (cryo-EM). All channels had water sites in the peptide toxin binding site, and an average of 75% of these sites were stable (low-energy), and 25% were unstable (medium or high energy). For the sodium channel, more unstable water sites were present in the deactivated state structure than the activated. Additionally, for each channel, unstable water sites coincided with the positions of peptide toxin residues that previous mutagenesis experiments had shown were important for activity. Finally, for the sodium channel in the deactivated state, unstable water sites were present in the peptide toxin binding pocket but did not overlap with the peptide toxin, suggesting that future experimental efforts could focus on targeting these sites to optimize potency.
| S-EPMC7599640 | biostudies-literature
Project description:One hypothesis accounting for major histocompatibility complex (MHC) restriction by T cell receptors (TCRs) holds that there are several evolutionary conserved residues in TCR variable regions that contact MHC. While this "germline codon" hypothesis is supported by various lines of evidence, it has been difficult to test. The difficulty stems in part from the fact that TCRs exhibit low affinities for pep/MHC, thus limiting the range of binding energies that can be assigned to these key interactions using mutational analyses. To measure the magnitude of binding energies involved, here we used high-affinity TCRs engineered by mutagenesis of CDR3. The TCRs included a high-affinity, MART-1/HLA-A2-specific single-chain TCR and two other high-affinity TCRs that all contain the same Vα region and recognize the same MHC allele (HLA-A2), with different peptides and Vβ regions. Mutational analysis of residues in CDR1 and CDR2 of the three Vα2 regions showed the importance of the key germline codon residue Y51. However, two other proposed key residues showed significant differences among the TCRs in their relative contributions to binding. With the use of single-position, yeast-display libraries in two of the key residues, MART-1/HLA-A2 selections also revealed strong preferences for wild-type germline codon residues, but several alternative residues could also accommodate binding and, hence, MHC restriction. Thus, although a single residue (Y51) could account for a proportion of the energy associated with positive selection (i.e., MHC restriction), there is significant plasticity in requirements for particular side chains in CDR1 and CDR2 and in their relative binding contributions among different TCRs.
| S-EPMC3844621 | biostudies-literature