Dataset Information

Novel knowledge-based mean force potential at the profile level.

ABSTRACT: BACKGROUND: The development and testing of functions for the modeling of protein energetics is an important part of current research aimed at understanding protein structure and function. Knowledge-based mean force potentials are derived from statistical analyses of interacting groups in experimentally determined protein structures. Current knowledge-based mean force potentials are developed at the atom or amino acid level. The evolutionary information contained in the profiles is not investigated. Based on these observations, a class of novel knowledge-based mean force potentials at the profile level has been presented, which uses the evolutionary information of profiles for developing more powerful statistical potentials. RESULTS: The frequency profiles are directly calculated from the multiple sequence alignments outputted by PSI-BLAST and converted into binary profiles with a probability threshold. As a result, the protein sequences are represented as sequences of binary profiles rather than sequences of amino acids. Similar to the knowledge-based potentials at the residue level, a class of novel potentials at the profile level is introduced. We develop four types of profile-level statistical potentials including distance-dependent, contact, Phi/Psi dihedral angle and accessible surface statistical potentials. These potentials are first evaluated by the fold assessment between the correct and incorrect models generated by comparative modeling from our own and other groups. They are then used to recognize the native structures from well-constructed decoy sets. Experimental results show that all the knowledge-base mean force potentials at the profile level outperform those at the residue level. Significant improvements are obtained for the distance-dependent and accessible surface potentials (5-6%). The contact and Phi/Psi dihedral angle potential only get a slight improvement (1-2%). Decoy set evaluation results show that the distance-dependent profile-level potentials even outperform other atom-level potentials. We also demonstrate that profile-level statistical potentials can improve the performance of threading. CONCLUSION: The knowledge-base mean force potentials at the profile level can provide better discriminatory ability than those at the residue level, so they will be useful for protein structure prediction and model refinement.

SUBMITTER: Dong Q

PROVIDER: S-EPMC1534065 | biostudies-literature | 2006

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Novel knowledge-based mean force potential at the profile level.

Dong Qiwen Q Wang Xiaolong X Lin Lei L

BMC bioinformatics 20060627

<h4>Background</h4>The development and testing of functions for the modeling of protein energetics is an important part of current research aimed at understanding protein structure and function. Knowledge-based mean force potentials are derived from statistical analyses of interacting groups in experimentally determined protein structures. Current knowledge-based mean force potentials are developed at the atom or amino acid level. The evolutionary information contained in the profiles is not inv ...[more]

PMID: 16803615

Similar Datasets

Project description:Murine double minute 2 (MDM2) proteins are found to be overproduced by many human tumors in order to inhibit the functioning of p53 molecules, a tumor suppressor protein. Thus, reactivating p53 functioning in cancer cells by disrupting p53-MDM2 interactions may offer a significant approach in cancer treatment. However, the structural characterization of the p53-MDM2 complex at the atomistic level and the mechanism of binding/unbinding of the p53-MDM2 complex still remain unclear. Therefore, we demonstrate here the probable binding (unbinding) pathway of transactivation domain 1 of p53 during the formation (dissociation) of the p53-MDM2 complex in terms of free energy as a function of reaction coordinate from the potential of mean force (PMF) study using two different force fields: ff99SB and ff99SB-ILDN. From the PMF plot, we noticed the PMF to have a minimum value at a p53-MDM2 separation of 12 Å, with a dissociation energy of 30 kcal mol-1. We also analyzed the conformational dynamics and stability of p53 as a function of its distance of separation from MDM2. The secondary structure content (helix and turns) in p53 was found to vary with its distance of separation from MDM2. The p53-MDM2 complex structure with lowest potential energy was isolated from the ensemble at the reaction coordinate corresponding to the minimum PMF value and subjected to molecular dynamics simulation to identify the interface surface area, interacting residues at the interface, and the stability of the complex. The simulation results highlight the importance of hydrogen bonds and the salt bridge between Lys94 of MDM2 and Glu17 of p53 in the stability of the p53-MDM2 complex. We also carried out the binding free energy calculations and the per residue energy decomposition analyses of the interface residues of the p53-MDM2 complex. We found that the binding affinity between MDM2 and p53 is indeed high [ΔG bind = -7.29 kcal mol-1 from molecular mechanics/Poisson-Boltzmann surface area (MM/PBSA) and ΔG bind = -53.29 kcal mol-1 from molecular mechanics/generalized borne surface area]. The total binding energy obtained using the MM/PBSA method was noticed to be closer to the experimental values (-6.4 to -9.0 kcal mol-1). The p53-MDM2 complex binding profile was observed to follow the same trend even in the duplicate simulation run and also in the simulation carried out with different force fields. We found that Lys51, Leu54, Tyr100, and Tyr104 from MDM2 and the residues Phe19, Trp23, and Leu26 from p53 provide the highest energy contributions for the p53-MDM2 interaction. Our findings highlight the prominent structural and binding characteristics of the p53-MDM2 complex that may be useful in designing potential inhibitors to disrupt the p53-MDM2 interactions.

Project description:Changes in free energy provide valuable information for molecular recognition, including both ligand-receptor binding thermodynamics and kinetics. Umbrella sampling (US), a widely used free energy calculation method, has long been used to explore the dissociation process of ligand-receptor systems and compute binding free energy. In existing publications, the binding free energy computed from the potential of mean force (PMF) with US simulation mostly yielded "ball park" values with experimental data. However, the computed PMF values are highly influenced by factors such as initial conformations and/or trajectories provided, the reaction coordinate, and sampling of conformational space in each US window. These critical factors have rarely been carefully studied. Here we used US to study the guest aspirin and 1-butanol dissociation processes of β-cyclodextrin (β-CD) and an inhibitor SB2 dissociation from a p38α mitogen-activated protein kinase (MAPK) complex. For β-CD, we used three different β-CD conformations to generate the dissociation path with US windows. For p38α, we generated the dissociation pathway by using accelerated molecular dynamics followed by conformational relaxing with short conventional MD, steered MD, and manual pulling. We found that, even for small β-CD complexes, different β-CD conformations altered the height of the PMF, but the pattern of PMF was not affected if the MD sampling in each US window was well-converged. Because changing the macrocyclic ring conformation needs to rotate dihedral angles in the ring, a bound ligand largely restrains the motion of cyclodextrin. Therefore, once a guest is in the binding site, cyclodextrin cannot freely change its initial conformation, resulting in different absolute heights of the PMF, which cannot be overcome by running excessively long MD simulations for each US window. Moreover, if the US simulations were not converged, the important barrier and minimum were missed. For ligand-protein systems, our studies also suggest that the dissociation trajectories modeled by an enhanced sampling method must maintain a natural molecular movement to avoid biased PMF plots when using US simulations.

Dataset Information

Novel knowledge-based mean force potential at the profile level.

Publications

Novel knowledge-based mean force potential at the profile level.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets