Biological processes often depend on ligand binding to receptors yet accurate calculation of the associated binding free energy remains a significant challenge of central importance to structure-based drug design. Water molecules participate in biological processes, including ligand binding to proteins, and the mechanical stability and aggregation of proteins. We propose to use molecular simulations and theory to provide deeper insights into water's role in the above processes. We propose to develop methods for the accurate calculation of protein-ligand binding affinities based on our recent progress towards developing a displaced solvent functional methodology (DSM). We propose to develop our "WaterMap" method into a practical tool for medicinal chemistry, by extending it to heterogeneous protein receptor interfaces with mixtures of hydrophilic, hydrophobic and neutral residues and to apply it to libraries of receptors and ligands. This work is part of an ongoing fruitful collaboration with the Friesner group at Columbia. To expedite this research we propose to combine our new colored noise multiple time step MD algorithm (CN-RESPA) with our replica exchange with solute tempering (REST) algorithm and the existing;-hopping free energy perturbation algorithm (;-hopping FEP) to provide an extremely powerful methodology for calculating binding efficiencies and sampling conformational states in the above research. Huntington's disease is caused by a genetic expansion of a region containing the CAG codon that leads to an increase in the length of the protein huntingtin's polyglutamine tract and a subsequent triggering of the formation of cytotoxic polyglutamine amyloid fibers and plaques. It is thought that a single rare misfolded polyglutamine polypeptide serves as a nucleus for the rapid formation of aggregated oligomers and finally the characteristic amyloid, but structural biologists have been unable to determine structures of these misfolded conformations. In preliminary simulations we found examples of highly collapsed long lived polyQ double-back structures that are mechanically resilient against large pulling forces (as is seen in AFM experiments). We propose to perform long enough molecular simulations using powerful computers like ANTON to pin down the structures of the mechanically resilient misfolded conformations of polyQ, to explore the balance between residue-water and residue-residue interactions to better understand why polyQ with its highly polar side-chain groups forms such compact structure, and, then, to get insights into conformational transitions of polyQ during aggregation.

Public Health Relevance

This proposal will focus on: developing novel methods for including water in the computation of protein-ligand binding affinities and protein active site modeling;
aim to understand the molecular biophysics of the precursors to aggregation and the structural, dynamical, and mechanical properties of polyglutamine repeats in Huntington's disease;and develop new sampling methods required for simulations of the relevant systems.

National Institute of Health (NIH)
National Institute of General Medical Sciences (NIGMS)
Research Project (R01)
Project #
Application #
Study Section
Macromolecular Structure and Function D Study Section (MSFD)
Program Officer
Preusch, Peter C
Project Start
Project End
Budget Start
Budget End
Support Year
Fiscal Year
Total Cost
Indirect Cost
Columbia University (N.Y.)
Other Domestic Higher Education
New York
United States
Zip Code
Stirnemann, Guillaume; Kang, Seung-gu; Zhou, Ruhong et al. (2014) How force unfolding differs from chemical denaturation. Proc Natl Acad Sci U S A 111:3413-8
Mondal, Jagannath; Morrone, Joseph A; Berne, B J (2013) How hydrophobic drying forces impact the kinetics of molecular recognition. Proc Natl Acad Sci U S A 110:13277-82
Wang, Lingle; Berne, B J; Friesner, Richard A (2012) On achieving high accuracy and reliability in the calculation of relative protein-ligand binding affinities. Proc Natl Acad Sci U S A 109:1937-42
Wang, Lingle; Berne, B J; Friesner, R A (2011) Ligand binding to protein-binding pockets with wet and dry regions. Proc Natl Acad Sci U S A 108:1326-30
Zhou, Ruhong; Li, Jingyuan; Hua, Lan et al. (2011) Comment on "urea-mediated protein denaturation: a consensus view". J Phys Chem B 115:1323-6; discussion 1327-8
Wang, Lingle; Friesner, Richard A; Berne, B J (2011) Replica exchange with solute scaling: a more efficient version of replica exchange with solute tempering (REST2). J Phys Chem B 115:9431-8
Young, Tom; Hua, Lan; Huang, Xuhui et al. (2010) Dewetting transitions in protein cavities. Proteins 78:1856-69
Wang, Lingle; Friesner, Richard A; Berne, B J (2010) Hydrophobic interactions in model enclosures from small to large length scales: non-additivity in explicit and implicit solvent models. Faraday Discuss 146:247-62; discussion 283-98, 395-
Wang, Lingle; Friesner, Richard A; Berne, B J (2010) Competition of electrostatic and hydrophobic interactions between small hydrophobes and model enclosures. J Phys Chem B 114:7294-301
Perez-Jimenez, Raul; Li, Jingyuan; Kosuri, Pallav et al. (2009) Diversity of chemical mechanisms in thioredoxin catalysis revealed by single-molecule force spectroscopy. Nat Struct Mol Biol 16:890-6

Showing the most recent 10 out of 22 publications