Jeff Wereszczynski

Molecular recognition plays a central role in biochemical processes. Although well studied, understanding the mechanisms of recognition is inherently difficult due to the range of potential interactions, the molecular rearrangement associated with binding, and the time and length scales involved. Computational methods have the potential for not only complementing experiments that have been performed, but also in guiding future ones through their predictive abilities. In this review, we discuss how molecular dynamics (MD) simulations may be used in advancing our understanding of the thermodynamics that drive biomolecular recognition. We begin with a brief review of the statistical mechanics that form a basis for these methods. This is followed by a description of some of the most commonly used methods: thermodynamic pathways employing alchemical transformations and potential of mean force calculations, along with end-point calculations for free energy differences, and harmonic and quasi-harmonic analysis for entropic calculations. Finally, a few of the fundamental findings that have resulted from these methods are discussed, such as the role of configurational entropy and solvent in intermolecular interactions, along with selected results of the model system T4 lysozyme to illustrate potential and current limitations of these methods.

Histone tails harbor a plethora of post-translational modifications that direct the function of chromatin regulators, which recognize them through effector domains. Effector domain/histone interactions have been broadly studied, but largely using peptide fragments of histone tails. Here, we extend these studies into the nucleosome context and find that the conformation adopted by the histone H3 tails is inhibitory to BPTF PHD finger binding. Using NMR spectroscopy and MD simulations, we show that the H3 tails interact robustly but dynamically with nucleosomal DNA, substantially reducing PHD finger association. Altering the electrostatics of the H3 tail via modification or mutation increases accessibility to the PHD finger, indicating that PTM crosstalk can regulate effector domain binding by altering nucleosome conformation. Together, our results demonstrate that the nucleosome context has a dramatic impact on signaling events at the histone tails, and highlights the importance of studying histone binding in the context of the nucleosome.

Hexasomes and tetrasomes are intermediates in nucleosome assembly and disassembly. Their formation is promoted by histone chaperones, ATP-dependent remodelers, and RNA polymerase II. In addition, hexasomes are maintained in transcribed genes and could be an important regulatory factor. While nucleosome composition has been shown to affect the structure and accessibility of DNA, its influence on histone tails is largely unknown. Here, we investigate the conformational dynamics of the H3 tail in the hexasome and tetrasome. Using a combination of NMR spectroscopy, MD simulations, and trypsin proteolysis, we find that the conformational ensemble of the H3 tail is regulated by nucleosome composition. As has been found for the nucleosome, the H3 tails bind robustly to DNA within the hexasome and tetrasome, but upon loss of the H2A/H2B dimer, we determined that the adjacent H3 tail has an altered conformational ensemble, increase in dynamics, and increase in accessibility. Similar to observations of DNA dynamics, this is seen to be asymmetric in the hexasome. Our results indicate that nucleosome composition has the potential to regulate chromatin signaling and ultimately help shape the chromatin landscape.