Skip to main content

First Principles Multiscale Analysis of Biochemical Processes: Signal Transduction and Spectroscopic Analysis of Local Structure-Folding Relations in Membrane Proteins


EMSL Project ID
20904

Abstract

We propose computational work in two application areas: first principles simulation of posttranslation modification in signal transduction in eukaryote cells and the use of spectroscopic methods to probe the mechanisms of protein folding in membrane proteins. Substantial effort will also be devoted to the necessary development of the multiscale high performance computational tools to address problems encountered in the course of our project. The proposed computational work will be carried out in close collaboration with experimental groups (Adams, Cole, and Taylor, signal transduction program; Kim, protein folding work).
The most important posttranslational modification affecting cellular function is the phosphorylation of residues in target proteins carried out by the kinase enzymes. The structure-function relationship among these enzymes has been an active topic of research for many years. The structure of their active site is highly conserved implying that the mechanism for phosphorylation is similar across the family. Despite intense experimental work the nature of the reaction process and the role of the conserved elements in the active site are not well understood. These issues will be addressed in this projects using recently developed QM/MM module(1) in the NWChem computational chemistry package to simulate the phosphorylation reactions. The atomic level accuracy provided by QM/MM methods can provide much needed interpretation for the considerable amount of structural, kinetic, chemical modification and mutagenesis data accumulated for these systems. The opposing process of dephosphorylation is carried out by phosphatases. Analysis of biochemical pathways with temporal and spatial resolution has become a major objective of "systems biology". Ab-initio methods are needed to supplement the sparse experimental data on rate constants to enable kinetic simulations that involve the actual molecular species of signaling networks. As part of this program we will investigate the modulation of cellular signaling through the interaction of phosphatases with reactive oxygen species (ROS); a system well-suited to this approach because (1) the fundamental reaction, oxidation of a catalytic cysteine residue, can be examined by QM/MM methods and (2) kinetic signaling data are available to test simulations results based on calculated rate parameters.
Another part of our project involves a combined experimental and computational investigation of membrane protein folding via fluorescence spectroscopy. Membrane proteins constitute around 30% of all cell proteins and are major targets of biological research. They facilitate the transfer of materials and information between cells and their environment and are essential to the integrity of biological interfaces. Despite their significance, the fundamental assembly mechanism of these ubiquitous biomolecules remains a mystery. Our program will focus on one of the best characterized membrane protein folding reactions, that of Outer Membrane Protein A (OmpA). This bacterial integral protein is the most abundant protein of the outer membrane of E. Coli with proposed functions that include acting as a nonspecific pore and receptor, and providing structural integrity to the biological membrane. In addition, OmpA belongs to a large family of pore-forming -barrel membrane proteins and hence, elucidation of the assembly mechanism of OmpA may shed light on the important porin proteins implicated in transport of toxic substances during bioremediation. OmpA is an ideal model system; it is one of the few membrane proteins that reversibly folds from the fully denatured state. High resolution x-ray and NMR structures are available identifying the residues and their environments in the protein. Because the protein folding can be reversibly controlled and the fluorescence is sensitive to the local folding environment we will obtain direct spectroscopic information about the local structure in various folded environments. Our computational approach is unique in providing the highly accurate finite temperature excited state properties of chromophores in large proteins that are required for the analysis of the fluorescence data. This new theoretical tool is based on the integration of coupled cluster methods with molecular dynamics simulations in the NWChem program.
In order make progress in these research areas our development of new multiscale computational tools for application to problems requiring high performance computing will have to continue. We believe that are entering a new era of discovery in which the information from necessarily high level experimental observations can be greatly expanded by appropriate multiscale theoretical interpretation beginning at the atomic first principles level. The applications proposed will also drive the development of new methods that will be included in NWChem software package. The availability of such tools benefits other areas important to the mission of the Department of Energy including biogeochemistry, subsurface science, waste isolation, interfacial chemistry and catalysis.

Project Details

Project type
Capability Research
Start Date
2006-09-30
End Date
2009-09-30
Status
Closed

Team

Principal Investigator

John Weare
Institution
University of California, San Diego

Team Members

Ryan Daly
Institution
University of California, San Diego

Philip Cheng
Institution
University of California, San Diego

Mitchell Luskin
Institution
University of Minnesota

Emilie Cauet
Institution
Université Libre de Bruxelles

Sergio Wong
Institution
University of California, San Diego

Wei Wang
Institution
University of California, San Diego

Peter Langhoff
Institution
University of California, San Diego

Karol Kowalski
Institution
Pacific Northwest National Laboratory

Judy Kim
Institution
University of California, San Diego

Weinan E
Institution
Princeton University

Phillip Cole
Institution
Johns Hopkins University

Scott Baden
Institution
University of California, San Diego

Julie Mitchell
Institution
University of Wisconsin, Madison

Brigitta Elsasser
Institution
University of Salzburg

Susan Taylor
Institution
University of California, San Diego

Joseph Adams
Institution
University of California, San Diego

T. Straatsma
Institution
Oak Ridge National Laboratory

Eric Bylaska
Institution
Pacific Northwest National Laboratory

J Mccammon
Institution
University of California, San Diego

John Miller
Institution
Washington State University Tri-Cities

Marat Valiev
Institution
Environmental Molecular Sciences Laboratory

Related Publications

Bylaska EJ, M Holst, and JH Weare. 2009. "Adaptive Finite Element Method for Solving the Exact Kohn-Sham Equation of Density Functional Theory." Journal of Chemical Theory and Computation 5(4):937-948. doi:10.1021/ct800350j
Cauet EL, M Valiev, and JH Weare. 2010. "Vertical Ionization Potentials of Nucleobases in a Fully Solvated DNA Environment." Journal of Physical Chemistry B 114(17):5886-5894.
Elsasser BM, and G Fels. 2010. "Atomistic Details of the Associative Phosphodiester Cleavage in Human Ribonuclease H." Physical Chemistry Chemical Physics. PCCP 12(36):11081-11088. doi:10.1039/c001097a
Elsasser BM, M Valiev, and JH Weare. 2009. "A Dianionic Phosphorane Intermediate and Transition States in an Associative AN+DN Mechanism for the RibonucleaseA Hydrolysis Reaction." Journal of the American Chemical Society 131(11):3869-3871. doi:10.1021/ja807940y
Valiev M, J Yang, J Adams, SS Taylor, and JH Weare. 2007. "Phosphorylation Reaction in cAPK Protein Kinase - Free Energy Quantum Mechanic/Molecular Mechanics Simulations." Journal of Physical Chemistry B 111(47):13455-13464. doi:10.1021/jp074853q