Related Experiment Video
Updated: Feb 7, 2026

Analyzing Melts and Fluids from Ab Initio Molecular Dynamics Simulations with the UMD Package
Published on: September 17, 2021
Valid molecular dynamics simulations of human hemoglobin require a surprisingly large box size
Krystel El Hage1, Florent Hédin1, Prashant K Gupta1
1Department of Chemistry, University of Basel, Basel, Switzerland.
Abstract:
Recent molecular dynamics (MD) simulations of human hemoglobin (Hb) give results in disagreement with experiment. Although it is known that the unliganded (T[Formula: see text]) and liganded (R[Formula: see text]) tetramers are stable in solution, the published MD simulations of T[Formula: see text] undergo a rapid quaternary transition to an R-like structure. We show that T[Formula: see text] is stable only when the periodic solvent box contains ten times more water molecules than the standard size for such simulations. The results suggest that such a large box is required for the hydrophobic effect, which stabilizes the T[Formula: see text] tetramer, to be manifested. Even in the largest box, T[Formula: see text] is not stable unless His146 is protonated, providing an atomistic validation of the Perutz model. The possibility that extra large boxes are required to obtain meaningful results will have to be considered in evaluating existing and future simulations of a wide range of systems.
Related Concept Videos
Requirements for Human Life
Oxygen
Atmospheric air is only about 20 percent oxygen, but that oxygen is a key component of the chemical reactions that keep the body alive, including the reactions that produce ATP. Brain cells are susceptible to a lack of oxygen because they require a...
Hemoglobin
When all four heme groups are bound to oxygen, the resulting molecule is called oxyhemoglobin. As a result, arterial blood...
Basic Postulates of Kinetic Molecular Theory: Particle Size, Energy, and Collision
Molecular Models
Reliability and Validity
Dynamic Equilibrium

