The binding of a fully flexible peptide, loop, or terminus of a protein to a relatively structured site on the surface of its molecular partner is a prevalent type of a transient macromolecular interaction that mediates the majority of recognition and signaling cascades. While spanning a wide range of affinities, such interactions uniformly involve a structured and an unstructured partner that may exist as independent entities in the absence of one another and only associate for a limited period of time resulting in conformational changes, recruitment or dissociation of other proteins to/from the complex, and activation or inhibition of the downstream pathways. Synthetic peptide modulators of protein-peptide interactions are frequently promising therapeutic candidates in cancer, inflammation and endocrine disorders. Despite the critical importance of transient protein-peptide interactions for biomedical research and therapeutic discovery, only a small fraction of these complexes are amenable to experimental structure determination. Therefore, only accurate peptide docking may lead to a breakthrough in structural understanding of protein-peptide interactions. Yet for peptides longer than 6-8 amino acids, the overwhelming size and complexity of the conformational search space, the inevitable inaccuracies in the binding site representation due to induced fit, inefficient or insufficiently thorough sampling, and accumulation of scoring function errors prohibit the accurate determination of peptide binding poses and interactions by computational methods. The present proposal intends to dramatically expand the range of peptide and protein sizes for which accurate complex geometry prediction can be achieved by global conformational optimization. This advance will be made by pursuing two Specific Aims:
(Aim 1) Development and optimization of a reliable peptide cross-docking procedure using chemical field-enhanced binding site representations and improved entropy calculations in the stochastic global conformational search in internal coordinates;
and (Aim 2) Extension of the developed protocol to specific biological projects with experimental validation of the predicted geometries and peptide variants. The innovative strategies proposed in Aim 1 include using a new force field, the enrichment of binding site representations with chemical fields, and the optimal conformational expansion of the pocket for induced fit. The targets of Aim 2 include complexes of Class A and B GPCRs with their protein and peptide modulators and G-protein ?i interactions with GEF peptides. The attainment of the aims of this proposal will lead to breakthrough advances in flexible macromolecular docking. It will result in valuable software tools, protocols and shared resources for the biological community. It will also lead to discovery of new peptide modulators of important therapeutic, immune, and diagnostic targets.

Public Health Relevance

Interactions between fully flexible peptides and relatively structured surfaces of their molecular partners are of paramount importance in biology, but extremely challenging for experimental structure determination. Until now, computational prediction of these interactions was practically impossible due to the complexity of the search space and the inevitable inaccuracies in representation of molecular partners. The present proposal seeks dramatic expansion of the range of molecular sizes for which accurate complex geometry prediction can be achieved by global conformational optimization.

National Institute of Health (NIH)
Research Project (R01)
Project #
Application #
Study Section
Macromolecular Structure and Function D Study Section (MSFD)
Program Officer
Preusch, Peter
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
University of California San Diego
Schools of Pharmacy
La Jolla
United States
Zip Code
Lin, Changsheng; Ear, Jason; Midde, Krishna et al. (2014) Structural basis for activation of trimeric Gi proteins by multiple growth factor receptors via GIV/Girdin. Mol Biol Cell 25:3654-71
McRobb, Fiona M; Kufareva, Irina; Abagyan, Ruben (2014) In silico identification and pharmacological evaluation of novel endocrine disrupting chemicals that act via the ligand-binding domain of the estrogen receptor ?. Toxicol Sci 141:188-97
Kufareva, Irina; Katritch, Vsevolod; Participants of GPCR Dock 2013 et al. (2014) Advances in GPCR modeling evaluated by the GPCR Dock 2013 assessment: meeting new challenges. Structure 22:1120-39
Chen, Yu-Chen; Totrov, Max; Abagyan, Ruben (2014) Docking to multiple pockets or ligand fields for screening, activity prediction and scaffold hopping. Future Med Chem 6:1741-55
Kennedy, Dylan P; McRobb, Fiona M; Leonhardt, Susan A et al. (2014) The second extracellular loop of the adenosine A1 receptor mediates activity of allosteric enhancers. Mol Pharmacol 85:301-9
Gabrielsen, Mari; Kurczab, Rafa?; Siwek, Agata et al. (2014) Identification of novel serotonin transporter compounds by virtual screening. J Chem Inf Model 54:933-43
Acharya, Chayan; Kufareva, Irina; Ilatovskiy, Andrey V et al. (2014) PeptiSite: a structural database of peptide binding sites in 4D. Biochem Biophys Res Commun 445:717-23
McRobb, Fiona M; SahagĂșn, Virginia; Kufareva, Irina et al. (2014) In silico analysis of the conservation of human toxicity and endocrine disruption targets in aquatic species. Environ Sci Technol 48:1964-72
Kufareva, Irina; Lenoir, Marc; Dancea, Felician et al. (2014) Discovery of novel membrane binding structures and functions. Biochem Cell Biol 92:555-63
Tan, Qiuxiang; Zhu, Ya; Li, Jian et al. (2013) Structure of the CCR5 chemokine receptor-HIV entry inhibitor maraviroc complex. Science 341:1387-90

Showing the most recent 10 out of 52 publications