Dataset Information

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation.

ABSTRACT: We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(beta)146, and they sum up to a total length of 5.6 micros. We observe spontaneous and reproducible T-->R quaternary transitions of the Hb tetramer and tertiary transitions of the alpha and beta subunits, as detected from principal component projections, from an RMSD measure, and from rigid body rotation analysis. The simulations reveal a marked asymmetry between the alpha and beta subunits. Using the mutual information as correlation measure, we find that the beta subunits are substantially more strongly linked to the quaternary transition than the alpha subunits. In addition, the tertiary populations of the alpha and beta subunits differ substantially, with the beta subunits showing a tendency towards R, and the alpha subunits showing a tendency towards T. Based on the simulation results, we present a transition pathway for coupled quaternary and tertiary transitions between the R and T conformations of Hb.

SUBMITTER: Hub JS

PROVIDER: S-EPMC2865513 | biostudies-literature | 2010 May

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation.

Hub Jochen S JS Kubitzki Marcus B MB de Groot Bert L BL

PLoS computational biology 20100506 5

We present molecular dynamics simulations of unliganded human hemoglobin (Hb) A under physiological conditions, starting from the R, R2, and T state. The simulations were carried out with protonated and deprotonated HC3 histidines His(beta)146, and they sum up to a total length of 5.6 micros. We observe spontaneous and reproducible T-->R quaternary transitions of the Hb tetramer and tertiary transitions of the alpha and beta subunits, as detected from principal component projections, from an RMS ...[more]

PMID: 20463873

Similar Datasets

Project description:Research conducted both prior to and after the emergence of the COVID-19 pandemic reveals a notable rise in human exposure to cleaning products, hand sanitizers, and personal care items. Moreover, there has been a corresponding increase in the environmental release of these chemicals. Cleaning and disinfecting products often contain quaternary ammonium compounds (QACs) with alkyl chains as long as 8-12 carbon atoms. The attachment of quaternary ammonium surfactants to the membrane resulted in the deformation of the bilayer and membrane disruption. Before interactions with cell membranes, these surfactant molecules may form different aggregates depending on their architecture. Interaction of surfactant monomers or clusters with the cell membrane changes the physiochemical properties of the biomembranes. To investigate this interaction and its influence on membrane properties, we conducted molecular dynamics simulations of cationic quaternary ammonium surfactants interacting with dipalmitoylphosphatidylcholine (DPPC) membranes. Our simulations revealed significant interactions between the surfactants and the phospholipids, leading to substantial alterations in the structure of the bilayer. The results are compared with the simulated anionic (SDS) and nonionic surfactants/bilayer systems. Various aspects were considered, including the aggregation process, migration behavior, and eventual equilibrium of these molecules at the interface between the membrane and water. This analysis used various techniques such as density profiles, distribution functions, cluster analysis, order parameters, hydrogen bonding (H-bonding), and mean-square displacements. The results indicate that while surfactants with shorter alkyl tails (N-(2-hydroxyethyl)-N,N-dimethyloctan-1-aminium chloride (HEDMOAC)) make strong hydrogen bonds with the phosphate group and ester oxygen of the phosphatidylcholine bilayer and enter toward the bilayer in the monomer form, surfactants of longer alkyl tails aggregated on the membrane head-water interface and interact minimally with the head groups of the DPPC bilayer. For DDEDMEAC, a quaternary ammonium surfactant with a hydrophobic alkyl chain consisting of two decanoate groups, alteration of the structural and dynamical properties of the bilayer is expected to be governed by two different factors. First, the structural order of DPPC increases as surfactant aggregates interact with the membrane head group. Second, the decrease in the order of the bilayer occurs due to the insertion of surfactant monomers within the hydrophobic region of the bilayer. Strong interactions between constituents of tetraoctylammonium bromide (TOABr) and lipid head groups lead to a reduction in interlipid interactions and order, which further results in increased porosity of cellular membranes. Understanding the extent of these interactions plays a pivotal role in the toxicological assessment of these surfactants.

Project description:Hemoglobin is one of the most widely studied proteins genetically, biochemically, and structurally. It is an oxygen carrying tetrameric protein that imparts the characteristic red color to blood. Each chain of hemoglobin harbors a heme group embedded in a hydrophobic pocket. Several studies have investigated structural variations present in mammalian hemoglobin and their functional implications. However, camel hemoglobin has not been thoroughly explored, especially from a structural perspective. Importantly, very little is known about how the heme group interacts with hemoglobin under varying conditions of osmolarity and temperature. Several experimental studies have indicated that the tense (T) state is more stable than the relaxed (R) state of hemoglobin under normal physiological conditions. Despite the fact that R state is less stable than the T state, no extensive structural dynamics studies have been performed to investigate global quaternary transitions of R state hemoglobin under normal physiological conditions. To evaluate this, several 500 ns all-atom molecular dynamics simulations were performed to get a deeper understanding of how camel hemoglobin behaves under stress, which it is normally exposed to, when compared to human hemoglobin. Notably, camel hemoglobin was more stable under physiological stress when compared to human hemoglobin. Additionally, when compared to camel hemoglobin, cofactor-binding regions of hemoglobin also exhibited more fluctuations in human hemoglobin under the conditions studied. Several differences were observed between the residues of camel and human hemoglobin that interacted with heme. Importantly, distal residues His58 of α hemoglobin and His63 of β hemoglobin formed more sustained interactions, especially at higher temperatures, in camel hemoglobin. These residues are important for oxygen binding to hemoglobin. Thus, this work provides insights into how camel and human hemoglobin differ in their interactions under stress.

Project description:Extra unmatched nucleotides (single base bulges) are common structural motifs in folded RNA molecules and can participate in RNA-ligand binding and RNA tertiary structure formation. Often these processes are associated with conformational transitions in the bulge region such as flipping out of the bulge base from an intrahelical stacked toward a looped out state. Knowledge of the flexibility of bulge structures and energetics of conformational transitions is an important prerequisite to better understand the function of this RNA motif. Molecular dynamics simulations were performed on single uridine and adenosine bulge nucleotides at the center of eight basepair RNA molecules and indicated larger flexibility of the bulge bases compared to basepaired regions. The umbrella sampling method was applied to study the bulge base looping out process and accompanying conformational and free energy changes. Looping out toward the major groove resulted in partial disruption of adjacent basepairs and was found to be less favorable compared to looping out toward the minor groove. For both uridine and adenosine bulges, a positive free energy change for full looping out was obtained which was approximately 1.5 kcal mol-1 higher in the case of the adenosine compared to the uridine bulge system. The simulations also indicated stable partially looped out states with the bulge bases located in the RNA minor groove and forming base triples with 5'-neighboring basepairs. In the case of the uridine bulge this state was more stable than the intrahelical stacked bulge structure. Induced looping out toward the minor groove involved crossing of an energy barrier of approximately 3.5 kcal mol-1 before reaching the base triple state. A continuum solvent analysis of intermediate bulge states indicated that electrostatic interactions stabilize looped out and base triple states, whereas van der Waals interactions and nonpolar contributions favor the stacked bulge conformation.

Dataset Information

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation.

Publications

Spontaneous quaternary and tertiary T-R transitions of human hemoglobin in molecular dynamics simulation.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets