Modeling sickle hemoglobin fibers as one chain of coarse-grained particles

Sickle cell disease (SCD) is caused by a single point mutation in the beta-chain hemoglobin gene, resulting in the presence of abnormal hemoglobin S (HbS) in the patients' red blood cells (RBCs). In the deoxygenated state, the defective hemoglobin tetramers polymerize forming stiff fibers which distort the cell and contribute to changes in its biomechanical properties. Because the HbS fibers are essential in the formation of the sickle RBC, their material properties draw significant research interests. Here, a solvent-free coarse-grain molecular dynamics (CGMD) model is introduced to simulate single HbS fibers as a chain of particles. First, we show that the proposed model is able to efficiently simulate the mechanical behavior of single HbS fibers. Then, the zippering process between two HbS fibers is studied and the effect of depletion forces is investigated. Simulation results illustrate that depletion forces play a role comparable to direct fiber-fiber interaction via Van der Waals forces. This proposed model can greatly facilitate studies on HbS polymerization, fiber bundle and gel formation as well as interaction between HbS fiber bundles and the RBC membrane.     

Nitric oxide reduces sickle hemoglobin polymerization: potential role of nitric oxide-induced charge alteration in depolymerization

We previously demonstrated that inhaling nitric oxide (NO) increases the oxygen affinity of sickle red blood cells (RBCs) in patients with sickle cell disease (SCD). Our recent studies found that NO lowered the P₅₀ values of sickle hemoglobin (HbS) hemolysates but did not increase methemoglobin (metHb) levels, supporting the role of NO, but not metHb, in the oxygen affinity of HbS. Here we examine the mechanism by which NO increases HbS oxygen affinity. Because anti-sickling agents increase sickle RBC oxygen affinity, we first determined whether NO exhibits anti-sickling properties. The viscosity of HbS hemolysates, measured by falling ball assays, increased upon deoxygenation; NO treatment reduced the increment. Multiphoton microscopic analyses showed smaller HbS polymers in deoxygenated sickle RBCs and HbS hemolysates exposed to NO. These results suggest that NO inhibits HbS polymer formation and has anti-sickling properties. Furthermore, we found that HbS treated with NO exhibits an isoelectric point similar to that of HbA, suggesting that NO alters the electric charge of HbS. NO–HbS adducts had the same elution time as HbA upon high performance liquid chromatography analysis. This study demonstrates that NO may disrupt HbS polymers by abolishing the excess positive charge of HbS, resulting in increased oxygen affinity.     

Potential of isoquercitrin as antisickling agent: a multi-spectroscopic,thermophoresis and molecular modeling approach

Abstract

Sickle cell disease is an inherited disease caused by point mutation in hemoglobin (β-globin gene). Under oxygen saturation, sickle hemoglobin form polymers, leading to rigid erythrocytes. The transition of the blood vessels is altered and initiated by the adhesion of erythrocytes, neutrophils and endothelial cells. Sickle Hemoglobin (HbS) polymerization is a major cause in red blood cells (RBC), promoting sickling and destruction of RBCs. Isoquercitrin, a medicinal bioactive compound found in various medicinal plants, has multiple health benefits. The present study examines the potential of isoquercitrin as an anti-sickle agent, showing a significant decrease in the rate of polymerization as well as sickling of RBCs. Isoquercitrin-induced graded alteration in absorbance and fluorescence of HbS, confirmed their interaction. A negative value of ΔG° strongly suggests that it is a spontaneous exothermic reaction induced by entropy. Negative ΔH° and positive ΔS° predicted that hydrogen and hydrophobic binding forces interfered with a hydrophobic microenvironment of β6Val leading to polymerization inhibition of HbS. HbS-Isoquercitrin complex exhibits helical structural changes leading to destabilization of the HbS polymer as confirmed by CD spectroscopy. MST and DSC results indicate greater changes in thermophoretic mobility and thermal stability of sickle hemoglobin in the presence of isoquercitrin, respectively. These findings were also supported by molecular simulation studies using DOCK6 and GROMACS. Hence, we can conclude that isoquercitrin interacts with HbS through hydrogen bonding, which leads to polymerization inhibition. Consequently, isoquercitrin could potentially be used as a medication for the treatment of sickle cell disease.

Communicated by Ramaswamy H. Sarma     

Molecular insights into the irreversible mechanical behavior of sickle hemoglobin

Sickle cell disease is caused by the amino acid substitution of glutamic acid to valine, which leads to the polymerization of deoxygenated sickle hemoglobin (HbS) into long strands. These strands are responsible for the sickling of red blood cells (RBCs), making blood hyper-coagulable leading to an increased chance of vaso-occlusive crisis. The conformational changes in sickled RBCs traveling through narrow blood vessels in a highly viscous fluid are critical in understanding; however, there are few studies that investigate the origins of the molecular mechanical behavior of sickled RBCs. In this work, we investigate the molecular mechanical properties of HbS molecules. A mechanical model was used to estimate the directional stiffness of an HbS molecule and the results were compared to adult human hemoglobin (HbA). The comparison shows a significant difference in strength between HbS and HbA, as well as anisotropic behavior of the hemoglobin molecules. The results also indicated that the HbS molecule experienced more irreversible mechanical behavior than HbA under compression. Further, we have characterized the elastic and compressive properties of a double stranded sickle fiber using six HbS molecules, and it shows that the HbS molecules are bound to each other through strong inter-molecular forces.     

The kinetics of nucleation and growth of sickle cell hemoglobin fibers

Polymerization of sickle cell hemoglobin (HbS) in deoxy state is one of the basic events in the pathophysiology of sickle cell anemia. For insight into the polymerization process, we monitor the kinetics of nucleation and growth of the HbS polymer fibers. We define a technique for the determination of the rates J and delay times theta of nucleation and the fiber growth rates R of deoxy-HbS fibers, based on photolysis of CO-HbS by laser illumination. We solve numerically time-dependent equations of heat conductance and CO transport, coupled with respective photo-chemical processes, during kinetics experiments under continuous illumination. After calibration with experimentally determined values, we define a regime of illumination ensuring uniform temperature and deoxy-HbS concentration, and fast (within <1 s) egress to steady conditions. With these procedures, data on the nucleation and growth kinetics have relative errors of <5% and are reproducible within 10% in independent experiments. The nucleation rates and delay times have steep, exponential dependencies on temperature. In contrast, the average fiber growth rates only weakly depend on temperature. The individual growth rates vary by up to 40% under identical conditions. These variations are attributed to instability of the coupled kinetics and diffusion towards the growing end of a fiber. The activation energy for incorporation of HbS molecules into a polymer is E(A)=50 kJ mol(-1), a low value indicating the significance of the hydrophobic contacts in the HbS polymer. More importantly, the contrast between the strong theta(T) and weak R(T) dependencies suggests that the homogenous nucleation of HbS polymers occurs within clusters of a precursor phase. This conclusion may have significant consequences for the understanding of the pathophysiology of sickle cell anemia and should be tested in further work.     

Modulation of the self-assembled structure of biomolecules: coarse grained molecular dynamics simulation

The mechanisms governing the self-assembled structure of biomolecules (single chain and bundle of chains) are studied with an AB copolymer model via the coarse grained molecular dynamics simulations. Non-local hydrophobic interaction is found to play a critical role in the pattern formation of the assembled structure of polymer chains. We show that the polymer structure could be controlled by adjusting the balance between local (short range) and non-local (long range) hydrophobic interaction which are influenced by various factors such as the sequences, chain length, stiffness, confinement, and the topology of polymers. In addition, the competition between the intrachain hydrophobic interaction and interchain hydrophobic interaction determines the structural transition of the chain bundles. This work may provide important insights into the fundamental physics in the structure control and the self-assembly of biomolecules for various practical applications.     

Free Heme and the Polymerization of Sickle Cell Hemoglobin

In search of novel control parameters for the polymerization of sickle cell hemoglobin (HbS), the primary pathogenic event of sickle cell anemia, we explore the role of free heme, which may be excessively released in sickle erythrocytes. We show that the concentration of free heme in HbS solutions typically used in the laboratory is 0.02-0.04 mole heme/mole HbS. We show that dialysis of small molecules out of HbS solutions arrests HbS polymerization. The addition of 100-260 μM of free heme to dialyzed HbS solutions leads to rates of nucleation and polymer fiber growth faster by two orders of magnitude than before dialysis. Toward an understanding of the mechanism of nucleation enhancement by heme, we show that free heme at a concentration of 66 μM increases by two orders of magnitude the volume of the metastable clusters of dense HbS liquid, the locations where HbS polymer nuclei form. These results suggest that spikes of the free heme concentration in the erythrocytes of sickle cell anemia patients may be a significant factor in the complexity of the clinical manifestations of sickle cell anemia. The prevention of free heme accumulation in the erythrocyte cytosol may be a novel avenue to sickle cell therapy.     

Hemoglobin synthesis studies of a family with α-thalassemia trait and sickle cell trait

The ratio of total globin alpha to beta chain synthesis was determined in reticulocytes isolated from the blood of the members of a black family, some of whom had sickle cell trait with low blood HbS concentrations (25-30%). The results support the hypothesis that sickle cell trait individuals with low HbS concentrations also carry a gene for alpha-thalassemia.     

Intermolecular contacts within sickle hemoglobin fibers

By combining X-ray crystallographic co-ordinates of sickle hemoglobin (HbS) molecules with three-dimensional reconstructions of electron micrographs of HbS fibers we have synthesized a model for the structure of the clinically relevant HbS fiber. This model largely accounts for the action of 55 point mutations of HbS whose effect on fiber formation has been studied. In addition, it predicts locations at which additional point mutations are likely to affect fiber formation. The number of intermolecular axial contacts decreases with radius until, at the periphery of the fiber, there are essentially no axial contacts. We suggest that this observation accounts for the limited radial growth of the HbS fiber and that a similar mechanism may be a factor in limiting the size of other helical particles. The methodology for the synthesis of the fiber model is applicable to other systems in which X-ray crystallographic and electron microscopic data are available.     

Measuring forces between protein fibers by microscopy

We propose a general scheme for measuring the attraction between mechanically frustrated semiflexible fibers by measuring their thermal fluctuations and shape. We apply this analysis to a system of sickle hemoglobin (HbS) fibers that laterally attract one another. These fibers appear to "zip" together before reaching mechanical equilibrium due to the existence of cross-links into a dilute fiber network. We are also able to estimate the rigidities of the fibers. These rigidities are found to be consistent with sickle hemoglobin "single" fibers 20 nm in diameter, despite recent experiments indicating that fiber bundling sometimes occurs. Our estimate of the magnitude of the interfiber attraction for HbS fibers is in the range 8 +/- 7 kBT/microm, or 4 +/- 3 k(B)T/microm if the fibers are assumed, a priori to be single fibers (such an assumption is fully consistent with the data). This value is sufficient to bind the fibers, overcoming entropic effects, although extremely chemically weak. Our results are compared to models for the interfiber attraction that include depletion and van der Waals forces. This technique should also facilitate a similar analysis of other filamentous protein assembles in the future, including beta-amyloid, actin, and tubulin.     

Self-assembly of surfactant aqueous solution confined in a Janus amphiphilic nanotube

Abstract

In recent years, Janus amphiphilic nanotubes with complex surfaces have been synthesized. However, the self-assembly behaviour of surfactant solutions confined in a Janus amphiphilic nanotube has not been investigated so far. We performed molecular simulations to investigate the morphologies and phase diagrams of a surfactant confined in Janus amphiphilic nanotube consisting of both hydrophobic and hydrophilic surfaces. We derived a phase diagram of the representative snapshots of equilibrium morphologies. Morphologies that were not observed in the bulk eventuated in confined systems. Moreover, the self-assembled structures were found to be dependent on the spatial confinement. Furthermore, the self-assembled structures confined in hydroneutral and Janus amphiphilic nanotubes were compared. The results suggested that the self-assembled structures confined in the Janus amphiphilic nanotube resembled that confined in the hydroneutral nanotube owing to a strong confinement effect. Further developments in controlling the morphologies and self-assemblies will greatly advance their applications of these materials in nanofluidic devices, or for nanopatterning.     

Mixed Exciton-Charge-Transfer States in Photosystem II: Stark Spectroscopy on Site-Directed Mutants

Sickle erythrocytes exhibit abnormal morphology and membrane mechanics under deoxygenated conditions due to the polymerization of hemoglobin S. We employed dissipative particle dynamics to extend a validated multiscale model of red blood cells (RBCs) to represent different sickle cell morphologies based on a simulated annealing procedure and experimental observations. We quantified cell distortion using asphericity and elliptical shape factors, and the results were consistent with a medical image analysis. We then studied the rheology and dynamics of sickle RBC suspensions under constant shear and in a tube. In shear flow, the transition from shear-thinning to shear-independent flow revealed a profound effect of cell membrane stiffening during deoxygenation, with granular RBC shapes leading to the greatest viscosity. In tube flow, the increase of flow resistance by granular RBCs was also greater than the resistance of blood flow with sickle-shape RBCs. However, no occlusion was observed in a straight tube under any conditions unless an adhesive dynamics model was explicitly incorporated into simulations that partially trapped sickle RBCs, which led to full occlusion in some cases.     

Two-step mechanism of homogeneous nucleation of sickle cell hemoglobin polymers

Sickle cell anemia is a debilitating genetic disease that affects hundreds of thousands of babies born each year worldwide. Its primary pathogenic event is the polymerization of a mutant, sickle cell, hemoglobin (HbS); and this is one of a line of diseases (Alzheimer's, Huntington's, prion, etc.) in which nucleation initiates pathophysiology. We show that the homogeneous nucleation of HbS polymers follows a two-step mechanism with metastable dense liquid clusters serving as precursor to the ordered nuclei of the HbS polymer. The evidence comes from data on the rates of fiber nucleation and growth and nucleation delay times, the interaction of fibers with polarized light, and mesoscopic metastable HbS clusters in solution. The presence of a precursor in the HbS nucleation mechanism potentially allows low-concentration solution components to strongly affect the nucleation kinetics. The variations of these concentrations in patients might account for the high variability of the disease in genetically identical patients. In addition, these components can potentially be utilized for control of HbS polymerization and treatment of the disease.     

Anisotropy in sickle hemoglobin fibers from variations in bending and twist

We have studied the variations of twist and bend in sickle hemoglobin fibers. We find that these variations are consistent with an origin in equilibrium thermal fluctuations, which allows us to estimate the bending and torsional rigidities and effective corresponding material moduli. We measure bending by electron microscopy of frozen hydrated fibers and find that the bending persistence length, a measure of the length of fiber required before it starts to be significantly bent due to thermal fluctuations, is 130microm, somewhat shorter than that previously reported using light microscopy. The torsional persistence length, obtained by re-analysis of previously published experiments, is found to be only 2.5microm. Strikingly this means that the corresponding torsional rigidity of the fibers is only 6x10(-27)Jm, much less than their bending rigidity of 5x10(-25)Jm. For (normal) isotropic materials, one would instead expect these to be similar. Thus, we present the first quantitative evidence of a very significant material anisotropy in sickle hemoglobin fibers, as might arise from the difference between axial and lateral contacts within the fiber. We suggest that the relative softness of the fiber with respect to twist deformation contributes to the metastability of HbS fibers: HbS double strands are twisted in the fiber but not in the equilibrium crystalline state. Our measurements inform a theoretical model of the thermodynamic stability of fibers that takes account of both bending and extension/compression of hemoglobin (double) strands within the fiber.     

The entropically favored osmotic "compression" of sickle cell hemoglobin gels

Contrary to the accurate, hard-sphere depiction of monomeric hemoglobin in solution, sickle cell hemoglobin (HbS) polymerization/gelation requires attention to molecular interactions. From the temperature dependence of the osmotic compressibility of HbS gels, we were able to extract the entropy increase for concentrating HbS in this phase. Normalized per mole of water removed, the entropy increase from gel compression DeltaS(gel) is four times the previously measured DeltaS(trans), for the transition from monomeric HbS solution to HbS gel. The positive entropy change cannot emerge from the assembly of hard spheres but can indicate remodeling of HbS fibers driven by release of ordered water. The fourfold difference in DeltaS(gel) and DeltaS(trans) suggests that the act of initial fiber/gel formation from monomeric solution differs from the process of further polymerization due to tighter packing within the gel phase.     

Fiber Depolymerization: Fracture, Fragments, Vanishing Times, and Stochastics in Sickle Hemoglobin

The well-characterized rates, mechanisms, and stochastics of nucleation-dependent polymerization of deoxyhemoglobin S (HbS) are important in governing whether or not vaso-occlusive sickle cell crises will occur. The less well studied kinetics of depolymerization may also be important, for example in achieving full dissolution of polymers in the lungs, in resolution of crises and/or in minimizing gelation-induced cellular damage. We examine depolymerization by microscopic observations on depolymerizing HbS fibers, by Monte Carlo simulations and by analytical characterization of the mechanisms. We show that fibers fracture. Experimental scatter of rates is consistent with stochastic features of the analytical model and Monte Carlo results. We derive a model for the distribution of vanishing times and also show the distribution of fracture-dependent fiber fragment lengths and its time dependence. We describe differences between depolymerization of single fibers and bundles and propose models for bundle dissolution. Our basic model can be extended to dissolution of gels containing many fibers and is also applicable to other reversible linear polymers that dissolve by random fracture and end-depolymerization. Under the model, conditions in which residual HbS polymers exist and facilitate repolymerization and thus pathology can be defined; whereas for normal polymers requiring cyclic polymerization and depolymerization for function, conditions for rapid cycling due to residual aggregates can be identified.     

O2 dependence of K+ transport in sickle cells: the effect of different cell populations and the substituted benzaldehyde 12C79.

The molecular basis of sickle cell disease (SCD) is well known but the pathophysiology is poorly understood. It remains intractable to therapy. Hyperactivity of several membrane transport systems, including the K+-Cl- cotransporter (termed KCC), cause HbS-containing red cells (termed HbS cells) to dehydrate and sickle, leading to the development of sickle cell crises (SCCs). Contrary to normal red cells (HbA cells), KCC in HbS cells is active at low O2 tensions (PO2s), remaining responsive to low pH or urea. Since these stimuli are usually encountered in hypoxic regions, the abnormal O2 dependence increases the contribution of KCC to dehydration, and hence development of SCCs. These differences with HbA cells may be due to the younger population of cells or to polymerization of HbS. We used 86Rb+ as a K+ congener to investigate the activity of KCC at different PO2s, and density gradient separation to investigate different red cell fractions. We found no correlation of O2 dependence with cell fractions. We also used the substituted benzaldehyde 12C79 to increase the O2 affinity of HbS and found that its effect on HbS O2 saturation and cell sickling correlated with that on both Cl--independent and Cl--dependent K+ transport, implying that, at low PO2s, KCC activity correlated with HbS polymerization. The importance of these results to understanding the pathophysiology of SCD, and for the design of chemotherapeutic agents to ameliorate or prevent SCC, is discussed.     

AKAP-Dependent Modulation of BCAM/Lu Adhesion on Normal and Sickle Cell Disease RBCs Revealed by Force Nanoscopy

Human normal and sickle red blood cells (RBCs) adhere with high affinity to the alpha5 chain of laminin (LAMA5) via the basal cell adhesion molecule/Lutheran (BCAM/Lu) receptor, which is implicated in vasoocclusive episodes in sickle cell disease and activated through the cyclic adenosine monophosphate (cAMP) signaling pathway. However, the effect of the cAMP pathway on the expression of active BCAM/Lu receptors at the single-molecule level is unknown. We established an in vitro technique, based on atomic force microscopy, which enables detection of single BCAM/Lu proteins on the RBC surface and measures the unbinding force between BCAM/Lu and LAMA5. We showed that the expression of active BCAM/Lu receptors is higher in homozygous sickle RBCs (SS-RBCs) than normal RBCs and that it is critically dependent on the cAMP signaling pathway on both normal and SS-RBCs. Of importance, we illustrated that A-kinase anchoring proteins are crucial for BCAM/Lu receptor activation. Furthermore, we found that SS-RBCs from hydroxyurea-treated patients show a lower expression of active BCAM/Lu receptors, a lower unbinding force to LAMA5, and insignificant stimulation by epinephrine as compared to SS-RBCs from untreated patients. To our knowledge, these findings may lead to novel antiadhesive targets for vasoocclusive episodes in sickle cell disease.     

AKAP-Dependent Modulation of BCAM/Lu Adhesion on Normal and Sickle Cell Disease RBCs Revealed by Force Nanoscopy

Human normal and sickle red blood cells (RBCs) adhere with high affinity to the alpha5 chain of laminin (LAMA5) via the basal cell adhesion molecule/Lutheran (BCAM/Lu) receptor, which is implicated in vasoocclusive episodes in sickle cell disease and activated through the cyclic adenosine monophosphate (cAMP) signaling pathway. However, the effect of the cAMP pathway on the expression of active BCAM/Lu receptors at the single-molecule level is unknown. We established an in vitro technique, based on atomic force microscopy, which enables detection of single BCAM/Lu proteins on the RBC surface and measures the unbinding force between BCAM/Lu and LAMA5. We showed that the expression of active BCAM/Lu receptors is higher in homozygous sickle RBCs (SS-RBCs) than normal RBCs and that it is critically dependent on the cAMP signaling pathway on both normal and SS-RBCs. Of importance, we illustrated that A-kinase anchoring proteins are crucial for BCAM/Lu receptor activation. Furthermore, we found that SS-RBCs from hydroxyurea-treated patients show a lower expression of active BCAM/Lu receptors, a lower unbinding force to LAMA5, and insignificant stimulation by epinephrine as compared to SS-RBCs from untreated patients. To our knowledge, these findings may lead to novel antiadhesive targets for vasoocclusive episodes in sickle cell disease.     

Molecular crowding limits the role of fetal hemoglobin in therapy for sickle cell disease

The dominant assumption central to most treatments for sickle cell anemia has been that replacement of sickle hemoglobin (HbS) by fetal hemoglobin (HbF) would have major clinical benefit. Using laser photolysis, we have measured polymerization kinetics including rates of homogeneous and heterogeneous nucleation on mixtures of 20% and 30% HbF with HbS. We find that the present model for polymerization, including molecular crowding, can accurately predict the rates of such mixtures, by using the single assumption that no significant amount of HbF enters the polymer. The effects of replacing HbS by HbF on the rates of polymer formation are found to be significantly lower than previous measurements appeared to indicate because the impact of the replacement is also highly dependent on the total hemoglobin concentration. This is because the molecular crowding of non-polymerizing HbF offsets substantially the effects of decreasing the concentration of HbS concentration, an effect that increases with concentration. Most strikingly, the demonstrated benefit of hydroxyurea therapy in slowing the kinetics of intracellular polymerization cannot be primarily due to enhanced HbF, but must have some other origin, which could itself represent a promising therapeutic approach.