(Left to Right) Po-Chao Wen, Nandan Haloi, and Emad Tajkhorshid of the University of Illinois at Urbana-Champaign
Supercomputer simulations run on the Stampede2 system of the Texas Advanced Computing Center (TACC) have helped scientists studying how cells work understand the way in which protein molecules in a sense dance with a partner to respond to signals and regulate each other’s activities.
Crucial to giving cells energy for life is the migration of a compound called adenosine triphosphate (ATP) out of the cell’s powerhouse, the mitochondria. And critical for this flow out to the power-hungry parts of the cell is the interaction between a protein enzyme called hexokinase-II (HKII) and proteins in the voltage-dependent anion channel (VDAC) found on the outer membrane of the mitochondria.
HPC simulations have revealed for the first time how VDAC binds to HKII. The work was supported by allocations awarded by the Extreme Science and Engineering Discovery Environment (XSEDE) on Stampede2. XSEDE is funded by the National Science Foundation.
This basic research in how proteins interact out of the cell’s powerhouses, the mitochondria, will help researchers understand the molecular basis of diseases such as cancer.
“We had strong evidence that they bind, but we didn’t know how they bind to each other,” said Emad Tajkhorshid, the J. Woodland Hastings Endowed Chair in Biochemistry at the University of Illinois at Urbana-Champaign. “That was the million-dollar question.”
Tajkhorshid co-authored a study published in Nature Communications Biology, June 2021. The study found that when the enzyme and the channel proteins bind to each other, the conduction of the channel changes and partially blocks the flow of ATP. Simulations on the Stampede2 system of TACC revealed this binding.
What’s more, the XSEDE-allocated Ranch system at TACC holds the off-site permanent file storage for the study’s data.
HPC simulations reveal how the cell’s mitochondrial voltage-dependent anion channel (VDAC) binds to the enzyme hexokinase-II (HKII). Artistic depiction of membrane binding of the cytosolic enzyme Hexokinase (light blue), followed by its complex formation with the integral membrane protein VDAC (dark blue), at the surface of the outer membrane of mitochondria. ATP (red) are phosphorylated by HKII. Credit: Haloi, N., Wen, PC., Cheng, Q. et al.
“If it wasn’t for XSEDE, we wouldn’t be studying many of these complex projects and biological systems because you simply can’t afford running the simulation. They usually require long simulations, and we need multiple copies of these simulations to be scientifically convincing. Without XSEDE it is impossible. We would have to go back to studying smaller systems,” Tajkhorshid said.
This work has implications for a deeper understanding of not only healthy cells, but cancer cells.
Basically, a cell needs ATP to metabolize glucose. It uses the ‘P’ to convert glucose to glucose phosphate, giving it a ‘handle’ that the cell can use. Hexokinase-II makes the conversion happen, binding at the mitochondrial channel to gobble up the ATP and phosphorylate it.
“We showed how the phosphorylation affects this process of binding between the two proteins. That was also verified experimentally,” Tajkhorshid said.
The VDAC channel is critical for efficient delivery of ATP directly to hexokinase. “It can work like a double-edged sword. For a healthy cell it’s good. For a cancer cell, it also helps the cell to promote and to proliferate,” he said.
Tajkhorshid’s team developed the most detailed and sophisticated model yet of the complex formed by the binding of HKII and VDAC, combining highest-resolution all-atom molecular dynamics simulations with coarser Brownian dynamics techniques. The system size of the VDAC-HKII complex was about 700,000 atoms, including the membrane. It’s about one-fifth of the diameter of the COVID-19 virus.
“What stands out in our approach is that we actually considered the cellular background of this interaction,” said Po-Chao Wen, a post-doctoral research associate at the NIH Center for Macromolecular Modeling and Bioinformatics, University of Illinois at Urbana-Champaign.
Wen explained that their simulation design started with the hypothesis that the VDAC protein in the outer membrane might interact at all with HKII, which is localized at a different part of the cell called the cytosol. They speculated that HKII should bind first to the membrane, drifting on it until it reaches a VDAC protein.
VDAC sitting on the membrane has already been well modeled, and the researchers built on this knowledge to break down the modelling of the HKII-VDAC complex into three parts, which initially focused on HKII.
To study how HKII binds to the mitochondrial outer membrane, they used all-atom molecular dynamics and a tool developed by their center called the Highly Mobile Membrane Model (HMMM), which deals with the membrane interaction.
Stampede2 (left) and Ranch (right). Credit: TACC
They next used Brownian dynamics to study how HKII drifts on the membrane to meet VDAC, creating many encounter/collision events between a sitting VDAC and a drifting HKII on a planar membrane.
“Then we used all-atom molecular dynamics to get a more refined model and specific size of the interaction to look for this particular protein-protein interaction,” Wen added. This helped them find the most stable complex of the two proteins formed.
“It seemed almost impossible when we started this process, because of the long timescales of milliseconds to seconds of the all-atom simulations,” said study co-author Nandan Haloi, a PhD student also at the center.
“These are really expensive calculations, which would require millions of dollars to set up independently. And you need to run on parallel supercomputers using our NAMD code, otherwise we could not reach the time scales that we needed,” Tajkhorshid said.
“We’re extremely happy with TACC and their support for not only just this work, but most of our projects and also our software development, tuning the software and making it faster. TACC has been wonderful in supporting us,” Tajkhorshid said.
TACC scientists work with the NIH Center for Macromolecular Modeling and Bioinformatics to constantly optimize the NAMD software, currently used by thousands of researchers.
The next steps in the research include more ambitious systems such as the fusion of two cells, important in understanding how neurons in the brain fire signals to each other; and how a virus such as the novel coronavirus fuses to the host cell.
Tajkhorshid’s group was awarded a Leadership Resource Allocation on the NSF-funded flagship supercomputer Frontera at TACC to investigate some of these ambitious projects.
Said Tajkhorshid: “We like to look at our work as a computational microscope that allows one to look at molecular systems and processes, how molecules come together, how they move, and how they change their structure to accomplish a particular function that people have been indirectly measuring experimentally. Supercomputers are essential in providing this level of detail, which we can use to understanding the molecular basis of diseases, drug discovery, and more.”
The study, “Structural basis of complex formation between mitochondrial anion channel VDAC1 and Hexokinase-II,” was published June 2021 in the journal Nature Communications Biology. The study authors are Nandan Haloi, Po-Chao Wen, and Emad Tajkhorshid of the University of Illinois at Urbana-Champaign; Qunli Cheng, Meiying Yang, Gayathri Natarajan, Amadou K. S. Camara, and Wai-Meng Kwok of the Medical College of Wisconsin.
This research was supported by the National Institutes of Health grants R01-HL131673 and P41- GM104601. Simulations in this study were performed using allocations at National Science Foundation Supercomputing Centers (XSEDE grant number MCA06N060) and the Blue Waters supercomputer of National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign.
source: Jorge Salazar, TACC
