Researchers use supercomputer to track pathways in myoglobin

June 30, 2008

Blacksburg, Va. - Some 50 years ago, after decades of effort, John Kendrew determined the structure of the small globular protein, myoglobin, which is responsible for oxygen storage in cells. For this discovery, he shared the Nobel Prize in chemistry with Max Perutz, who did similar work on hemoglobin. But a mystery remained: Exactly what paths does oxygen follow as it moves in and out of myoglobin?

An interdisciplinary team led by researchers at Virginia Tech has provided a computational solution to the decades-old puzzle.

Oxygen binds to an iron atom deep within myoglobin. How does the oxygen reach the iron through the solid protein wall? Thirty years of experimental and theoretical work have indicated that, as the protein fluctuates, transient "tunnels" may open up to form a network of pathways. But the exact nature and location of these dynamic pathways at atomic resolution have remained a mystery - until now, as reported in the online Early Edition of the Proceedings of the National Academy of Sciences (PNAS) the week of June 30, 2008.

"Every decade has seen breakthroughs in solving the puzzles this protein presented. Yet, exactly how small ligands such as oxygen or carbon monoxide (CO) migrate to the deeply buried binding site and out again was never completely resolved, in my opinion," said Alexey Onufriev of Blacksburg, Va., assistant professor of computer science and of physics at Virginia Tech.

He began his myoglobin research in 1997 as a postdoctoral researcher in chemistry faculty member Michael Prisant's lab at Duke University. According to Onufriev, Prisant suggested a computational geometry volume analysis. "That effort showed us what to look for," said Onufriev.

Ten years later, the Virginia Tech researchers proposed a multi-level computational solution to provide a detailed, complete, atom-level picture of ligand migration in myoglobin. "It ties together much previous experimental and theoretical understanding," said Onufriev.

The study is likely to be the first in which molecular dynamics succeeded to directly simulate the ligand migration all the way in and out of the protein at room temperature at atomic resolution.

The simulations provided insights into several long-standing questions about myoglobin. They revealed that ligand migration occurs along two pathways with many branches and that ligands move into myoglobin using the same pathways that they use to move out of the protein. The simulations showed it is mainly local and short-lived fluctuations that are responsible for the opening and closing of these dynamic pathways. The simulations clearly show the pathways occur mostly between and around alpha helices - the protein scaffolds. Like people looking for an exit from an office complex, they cannot go through the walls or columns. Ligands pause at transient docking sites for 10s of nanoseconds. "Some of these sites were previously known but we found more," Onufriev said. "They are where pathways fork."

Another question was how the ligands moved. Do they force their way through the protein like a child in a play land ball pit or do they just go where the doors open? "Ligands follow the openings as they occur in the protein," Onufriev said.


The work has also advanced methodology for the study of proteins. Myoglobin is the model protein - the guinea pig - of structural biology. "It has been thoroughly studied and has been key in elucidating many basic principles of the structure-to-function paradigm," said Onufriev.

The group used tried-and-true molecular dynamic (MD) simulations complemented by a relatively simple, but robust computational geometry analysis of structural fluctuations available for ligand migration in myoglobin. "Our solution is on the edge of what is currently possible computation-wise, so Virginia Tech's SystemX supercomptuer was very helpful," Onufriev said.

The researchers used CO because it is easier to work with than oxygen and because most past experiments had also used CO. "The consensus is that small non-polar ligands such as oxygen, CO and nitric oxide follow basically the same pathways as oxygen ligands," Onufriev said.

One type of computational experiment initially puts ligands where they are bound, cuts the bond, and observes some of them wander out of the myoglobin protein. "Alternatively, we can place the ligand in water outside of the protein and occasionally observe it enter the protein and make its way all the way to the binding site," said Onufriev. "There is enough computational capacity in modern supercomputers to simulate enough of these events at room temperature to make statistically justified conclusions," he said.

"Another type of computation uses simple but highly-effective techniques from computational geometry and graph theory to look for space available for ligand migration based on free volume in the protein. For this one, we used snapshots from an MD simulation to find cavities in the myoglobin protein big enough to hold the ligand. Then we 'connected the dots' by computing the union of the cavities across all the snapshots and we clearly saw the pathways. These looked very similar to the ones we saw directly by following individual ligand trajectories," Onufriev said. "Over all, both methods complemented each other to provide a clear, simple picture of what was going on."

The methodology developed in this work can be used to study other systems. "There are many biomedically relevant systems that are much less understood, particularly as it concerns pathways, which are critical to protein function," said Onufriev.
Article and authors

The article, "Atomic level computational identification of ligand migration pathways between solvent and binding site in myoglobin," was co-authored by Jory Z. Ruscio, Deept Kumar, Maulik Shukla., Michael G. Prisant, T. M. Murali, and Alexey Onufriev.

J. Ruscio, of Leesburg, Va., a Ph.D. student in the Genetics, Bioinformatics, and Computational Biology graduate program, and computer science students Kumar, of Lucknow, India, and Shukla, of Ahmedabad, India, became interested in the myoglobin puzzle as students in one of Onufriev's classes, Ruscio said.

"This work would have passed out of view if Alexey (Onufriev) had not been stubborn, tenacious, and dedicated," said Prisant.

Ruscio and Onufriev emphasize that the success of this work was due to the interdisciplinary nature of the group. "This project would not have come together as well if we were missing any part of the team," Ruscio said. "Each of us brought our expertise to use computational methods, both pure computer science and biophysical, to study the migration pathways on a resolution that is not quite attainable to experimentalists. And we were very lucky to have a resource such as Virginia Tech's own supercomputer System X," she said.

The myoglobin project became part of Ruscio's Ph.D. research on Applications of Molecular Modeling to Protein Structure and Function. She is now a postdoctoral researcher at the University of California, Berkeley.

Kumar, a Ph.D. student, and Shukla, a master's degree student, have also since graduated. Kumar is an analyst for a startup operating out of San Francisco. Shukla is a software engineer with the Virginia Bioinformatics Institute at Virginia Tech.

Prisant, of Durham, N.C., is a visiting professor with the biochemistry department at Duke University Medical Center. Murali, originally from Chennai, India., is an assistant professor of computer science at Virginia Tech. Onufriev is originally from Moscow, Russia.

Learn more about Alexey Onufriev's work at Learn more about the 1962 Nobel Prize in Chemistry here:

Virginia Tech

Related Protein Articles from Brightsurf:

The protein dress of a neuron
New method marks proteins and reveals the receptors in which neurons are dressed

Memory protein
When UC Santa Barbara materials scientist Omar Saleh and graduate student Ian Morgan sought to understand the mechanical behaviors of disordered proteins in the lab, they expected that after being stretched, one particular model protein would snap back instantaneously, like a rubber band.

Diets high in protein, particularly plant protein, linked to lower risk of death
Diets high in protein, particularly plant protein, are associated with a lower risk of death from any cause, finds an analysis of the latest evidence published by The BMJ today.

A new understanding of protein movement
A team of UD engineers has uncovered the role of surface diffusion in protein transport, which could aid biopharmaceutical processing.

A new biotinylation enzyme for analyzing protein-protein interactions
Proteins play roles by interacting with various other proteins. Therefore, interaction analysis is an indispensable technique for studying the function of proteins.

Substituting the next-best protein
Children born with Duchenne muscular dystrophy have a mutation in the X-chromosome gene that would normally code for dystrophin, a protein that provides structural integrity to skeletal muscles.

A direct protein-to-protein binding couples cell survival to cell proliferation
The regulators of apoptosis watch over cell replication and the decision to enter the cell cycle.

A protein that controls inflammation
A study by the research team of Prof. Geert van Loo (VIB-UGent Center for Inflammation Research) has unraveled a critical molecular mechanism behind autoimmune and inflammatory diseases such as rheumatoid arthritis, Crohn's disease, and psoriasis.

Resurrecting ancient protein partners reveals origin of protein regulation
After reconstructing the ancient forms of two cellular proteins, scientists discovered the earliest known instance of a complex form of protein regulation.

Sensing protein wellbeing
The folding state of the proteins in live cells often reflect the cell's general health.

Read More: Protein News and Protein Current Events is a participant in the Amazon Services LLC Associates Program, an affiliate advertising program designed to provide a means for sites to earn advertising fees by advertising and linking to