Hierarchy of simulation models in predicting molecular recognition mechanisms from the binding energy landscapes: structural analysis of the peptide complexes with SH2 domains

DbPTM 3.0: an informative resource for investigating substrate site specificity and functional association of protein post-translational modifications

Protein modification is an extremely important post-translational regulation that adjusts the physical and chemical properties, conformation, stability and activity of a protein; thus altering protein function. Due to the high throughput of mass spectrometry (MS)-based methods in identifying site-specific post-translational modifications (PTMs), dbPTM (http://dbPTM.mbc.nctu.edu.tw/) is updated to integrate experimental PTMs obtained from public resources as well as manually curated MS/MS peptides associated with PTMs from research articles. Version 3.0 of dbPTM aims to be an informative resource for investigating the substrate specificity of PTM sites and functional association of PTMs between substrates and their interacting proteins. In order to investigate the substrate specificity for modification sites, a newly developed statistical method has been applied to identify the significant substrate motifs for each type of PTMs containing sufficient experimental data. According to the data statistics in dbPTM, >60% of PTM sites are located in the functional domains of proteins. It is known that most PTMs can create binding sites for specific protein-interaction domains that work together for cellular function. Thus, this update integrates protein–protein interaction and domain–domain interaction to determine the functional association of PTM sites located in protein-interacting domains. Additionally, the information of structural topologies on transmembrane (TM) proteins is integrated in dbPTM in order to delineate the structural correlation between the reported PTM sites and TM topologies. To facilitate the investigation of PTMs on TM proteins, the PTM substrate sites and the structural topology are graphically represented. Also, literature information related to PTMs, orthologous conservations and substrate motifs of PTMs are also provided in the resource. Finally, this version features an improved web interface to facilitate convenient access to the resource.

Genome-wide prediction of SH2 domain targets using structural information and the FoldX algorithm

Current experiments likely cover only a fraction of all protein-protein interactions. Here, we developed a method to predict SH2-mediated protein-protein interactions using the structure of SH2-phosphopeptide complexes and the FoldX algorithm. We show that our approach performs similarly to experimentally derived consensus sequences and substitution matrices at predicting known in vitro and in vivo targets of SH2 domains. We use our method to provide a set of high-confidence interactions for human SH2 domains with known structure filtered on secondary structure and phosphorylation state. We validated the predictions using literature-derived SH2 interactions and a probabilistic score obtained from a naive Bayes integration of information on coexpression, conservation of the interaction in other species, shared interaction partners, and functions. We show how our predictions lead to a new hypothesis for the role of SH2 domains in signaling.

Understanding the functional role of every protein in the cell is a long-standing goal of cellular biology. An important step in this direction is to discover how and when proteins interact inside the cell to accomplish their tasks. Many of the cellular functions depend on reversible protein modifications like phosphorylation. To sense these modifications, cells have protein domains capable of binding phosphorylated proteins such as the SH2 domain. In this work, we show that it is possible to use the three-dimensional structure of protein domains to predict its binding preferences. Using a computational tool called FoldX, we have predicted the binding specificity of several human SH2 domains. These predictions, based on the computational analysis of the 3-D structure, were shown to be of similar accuracy as those obtained from experimental binding assays. We show here that it is also possible to understand how a mutation changes the binding preference of protein binding domains, opening the way for better understanding of some disease causing mutations. The combination of this novel computational approach with other sources of information allowed us to provide a set of high-confidence novel interactions for the proteins here studied.

Computational analysis of ligand binding dynamics at the intermolecular hot spots with the aid of simulated tempering and binding free energy calculations

Equilibrium binding dynamics is studied for a panel of benzimidazole-containing compounds at the remodeled interface between human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp), engineered by mutating to glycine hot spot residues T175 from the hormone and W104 from the receptor. Binding energetics is predicted in a good agreement with the experimental data for a panel of these small molecules that complement the engineered defect and restore the binding affinity of the wild-type hGH-hGHbp complex. The results of simulated tempering ligand dynamics at the protein-protein interface reveals a diversity of ligand binding modes that is consistent with the structural orientation of the benzimidazole ring which closely mimics the position of the mutated W104 hot spot residue in the wild-type hGH-hGHbp complex. This structural positioning of the benzimidazole core motif is shown to be a critical feature of the low-energy ligand conformations binding in the engineered cavity. The binding free energy analysis provides a plausible rationale behind the experimental dissociation constants and suggests a key role of ligand-protein van der Waals interactions in restoring binding affinity.

Computational detection of the binding-site hot spot at the remodeled human growth hormone-receptor interface

A hierarchical computational approach is used to identify the engineered binding-site cavity at the remodeled intermolecular interface between the mutants of human growth hormone (hGH) and the extracellular domain of its receptor (hGHbp). Multiple docking simulations are conducted with the remodeled hGH-hGHbp complex for a panel of potent benzimidazole-containing inhibitors that can restore the binding affinity of the wild-type complex, and for a set of known nonactive small molecules that contain different heterocyclic motifs. Structural clustering of ligand-bound conformations and binding free-energy calculations, using the AMBER force field and a continuum solvation model, can rapidly locate and screen numerous ligand-binding modes on the protein surface and detect the binding-site hot spot at the intermolecular interface. Structural orientation of the benzimidazole motif in the binding-site cavity closely mimics the position of the hot spot residue W104 in the crystal structure of the wild-type complex, which is recognized as an important structural requirement for restoring binding affinity. Despite numerous pockets on the protein surface of the mutant hGH-hGHbp complex, the binding-site cavity presents the energetically favorable hot spot for the benzimidazole-containing inhibitors, whereas for a set of nonactive molecules, the lowest energy ligand conformations do not necessarily bind in the engineered cavity. The results reveal a dominant role of the intermolecular van der Waals interactions in providing favorable ligand-protein energetics in the redesigned interface, in agreement with the experimental and computational alanine scanning of the hGH-hGHbp complex.

Simulating disorder-order transitions in molecular recognition of unstructured proteins: where folding meets binding

A microscopic study of functional disorder-order folding transitions coupled to binding is performed for the p27 protein, which derives a kinetic advantage from the intrinsically disordered unbound form on binding with the phosphorylated cyclin A-cyclin-dependent kinase 2 (Cdk2) complex. Hierarchy of structural loss during p27 coupled unfolding and unbinding is simulated by using high-temperature Monte Carlo simulations initiated from the crystal structure of the tertiary complex. Subsequent determination of the transition-state ensemble and the proposed atomic picture of the folding mechanism coupled to binding provide a microscopic rationale that reconciles the initiation recruitment of p27 at the cyclin A docking site with the kinetic benefit for a disordered alpha-helix in the unbound form of p27. The emerging structural polarization in the ensemble of unfolding/unbinding trajectories and in the computationally determined transition-state ensemble is not determined by the intrinsic folding preferences of p27 but rather is attributed to the topological requirements of the native intermolecular interface to order beta-hairpin and beta-strand of p27 that could be critical for nucleating rapid folding transition coupled to binding. In agreement with the experimental data, the disorder-order folding transition for p27 is largely determined by the functional requirement to form a specific intermolecular interface that ultimately dictates the folding mechanism and overwhelms any local folding preferences for creating a stable alpha-helix in the p27 structure before overcoming the major free energy barrier.

Monte Carlo simulations of the peptide recognition at the consensus binding site of the constant fragment of human immunoglobulin G: the energy landscape analysis of a hot spot at the intermolecular interface

Monte Carlo simulations of molecular recognition at the consensus binding site of the constant fragment (Fc) of human immunoglobulin G (Ig) protein have been performed to analyze structural and thermodynamic aspects of binding for the 13-residue cyclic peptide DCAWHLGELVWCT. The energy landscape analysis of a hot spot at the intermolecular interface using alanine scanning and equilibrium-simulated tempering dynamics with the simplified, knowledge-based energy function has enabled the role of the protein hot spot residues in providing the thermodynamic stability of the native structure to be determined. We have found that hydrophobic interactions between the peptide and the Met-252, Ile-253, His-433, and His-435 protein residues are critical to guarantee the thermodynamic stability of the crystallographic binding mode of the complex. Binding free energy calculations, using a molecular mechanics force field and a solvation energy model, combined with alanine scanning have been conducted to determine the energetic contribution of the protein hot spot residues in binding affinity. The conserved Asn-434, Ser-254, and Tyr-436 protein residues contribute significantly to the binding affinity of the peptide-protein complex, serving as an energetic hot spot at the intermolecular interface. The results suggest that evolutionary conserved hot spot protein residues at the intermolecular interface may be partitioned in fulfilling thermodynamic stability of the native binding mode and contributing to the binding affinity of the complex.

Complexity and simplicity of ligand-macromolecule interactions: the energy landscape perspective

The energy landscape approach has contributed to recent progress in understanding the complexity and simplicity of ligand-macromolecule interactions. Significant advances in computational structure prediction of ligand-protein complexes have been made using approaches that include the effects of protein flexibility and incorporate a hierarchy of energy functions. The results suggest that the complexity of structure prediction in molecular recognition may be determined by low-resolution properties of the underlying binding energy landscapes and by the nature of the energy funnels near the native structures of the complexes.