Matching Items (2)
150032-Thumbnail Image.png
Description
Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into the free energy of the system, the dynamical stability of

Molecular dynamics (MD) simulations provide a particularly useful approach to understanding conformational change in biomolecular systems. MD simulations provide an atomistic, physics-based description of the motions accessible to biomolecular systems on the pico- to micro-second timescale, yielding important insight into the free energy of the system, the dynamical stability of contacts and the role of correlated motions in directing the motions of the system. In this thesis, I use molecular dynamics simulations to provide molecular mechanisms that rationalize structural, thermodynamic, and mutation data on the interactions between the lac repressor headpiece and its O1 operator DNA as well as the ERK2 protein kinase. I performed molecular dynamics simulations of the lac repressor headpiece - O1 operator complex at the natural angle as well as at under- and overbent angles to assess the factors that determine the natural DNA bending angle. I find both energetic and entropic factors contribute to recognition of the natural angle. At the natural angle the energy of the system is minimized by optimization of protein-DNA contacts and the entropy of the system is maximized by release of water from the protein-DNA interface and decorrelation of protein motions. To identify the mechanism by which mutations lead to auto-activation of ERK2, I performed a series of molecular dynamics simulations of ERK1/2 in various stages of activation as well as the constitutively active Q103A, I84A, L73P and R65S ERK2 mutants. My simulations indicate the importance of domain closure for auto-activation and activity regulation. My results enable me to predict two loss-of-function mutants of ERK2, G83A and Q64C, that have been confirmed in experiments by collaborators. One of the powerful capabilities of MD simulations in biochemistry is the ability to find low free energy pathways that connect and explain disparate structural data on biomolecular systems. An extention of the targeted molecular dynamics technique using constraints on internal coordinates will be presented and evaluated. The method gives good results for the alanine dipeptide, but breaks down when applied to study conformational changes in GroEL and adenylate kinase.
ContributorsBarr, Daniel Alan (Author) / van der Vaart, Arjan (Thesis advisor) / Matyushov, Dmitry (Committee member) / Wolf, George (Committee member) / Shumway, John (Committee member) / Arizona State University (Publisher)
Created2011
155215-Thumbnail Image.png
Description
Proteins are essential for most biological processes that constitute life. The function of a protein is encoded within its 3D folded structure, which is determined by its sequence of amino acids. A variation of a single nucleotide in the DNA during transcription (nSNV) can alter the amino acid sequence (i.e.,

Proteins are essential for most biological processes that constitute life. The function of a protein is encoded within its 3D folded structure, which is determined by its sequence of amino acids. A variation of a single nucleotide in the DNA during transcription (nSNV) can alter the amino acid sequence (i.e., a mutation in the protein sequence), which can adversely impact protein function and sometimes cause disease. These mutations are the most prevalent form of variations in humans, and each individual genome harbors tens of thousands of nSNVs that can be benign (neutral) or lead to disease. The primary way to assess the impact of nSNVs on function is through evolutionary approaches based on positional amino acid conservation. These approaches are largely inadequate in the regime where positions evolve at a fast rate. We developed a method called dynamic flexibility index (DFI) that measures site-specific conformational dynamics of a protein, which is paramount in exploring mechanisms of the impact of nSNVs on function. In this thesis, we demonstrate that DFI can distinguish the disease-associated and neutral nSNVs, particularly for fast evolving positions where evolutionary approaches lack predictive power. We also describe an additional dynamics-based metric, dynamic coupling index (DCI), which measures the dynamic allosteric residue coupling of distal sites on the protein with the functionally critical (i.e., active) sites. Through DCI, we analyzed 200 disease mutations of a specific enzyme called GCase, and a proteome-wide analysis of 75 human enzymes containing 323 neutral and 362 disease mutations. In both cases we observed that sites with high dynamic allosteric residue coupling with the functional sites (i.e., DARC spots) have an increased susceptibility to harboring disease nSNVs. Overall, our comprehensive proteome-wide analysis suggests that incorporating these novel position-specific conformational dynamics based metrics into genomics can complement current approaches to increase the accuracy of diagnosing disease nSNVs. Furthermore, they provide mechanistic insights about disease development. Lastly, we introduce a new, purely sequence-based model that can estimate the dynamics profile of a protein by only utilizing coevolution information, eliminating the requirement of the 3D structure for determining dynamics.
ContributorsButler, Brandon Mac (Author) / Ozkan, S. Banu (Thesis advisor) / Vaiana, Sara (Committee member) / Ghirlanda, Giovanna (Committee member) / Ros, Robert (Committee member) / Arizona State University (Publisher)
Created2016