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; Stirnemann, Guillaume; Berne, B J (2013) When does trimethylamine N-oxide fold a polymer chain and urea unfold it? J Phys Chem B 117:8723-32
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
Stirnemann, Guillaume; Giganti, David; Fernandez, Julio M et al. (2013) Elasticity, structure, and relaxation of extended proteins under force. Proc Natl Acad Sci U S A 110:3847-52
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
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
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
Abel, Robert; Wang, Lingle; Friesner, Richard A et al. (2010) A displaced-solvent functional analysis of model hydrophobic enclosures. J Chem Theory Comput 6:2924-2934
Li, Jingyuan; Fernandez, Julio M; Berne, B J (2010) Water's role in the force-induced unfolding of ubiquitin. Proc Natl Acad Sci U S A 107:19284-9

Showing the most recent 10 out of 28 publications