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.     

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.     

The role of pH on instability and aggregation of sickle hemoglobin solutions

Understanding the physical basis of protein aggregation covers strong physical and biomedical interests. Sickle hemoglobin (HbS) is a point-mutant form of normal human adult hemoglobin (HbA). It is responsible for the first identified "molecular disease," as its propensity to aggregation is responsible for sickle cell disease. At moderately higher than physiological pH value, this propensity is inhibited: The rate of aggregate nucleation becomes exceedingly small and solubility after polymerization increases. These order-of-magnitude effects on polymer nucleation rates and concurrent relatively modest changes of solubility after polymerization are here shown to be related to both pH-induced changes of location and shape of the liquid-liquid demixing (LLD) region. This allows establishment of a self-consistent contact between the thermodynamics of the solution as such (i.e., the LLD region), the kinetics of fiber nucleation, the theory of percolation, and the thermodynamics of gelation. The observed pH-induced changes are largely attributable to strong perturbations of hydrophobic hydration configurations and related free energy by electric charges. Similar mechanisms of effective control of aggregate nucleation rates by means of agents such as cosolutes, pH, salts, and additives, shifting the LLD and associated regions of anomalous fluctuations, promise to be relevant to the whole field of protein aggregation pathologies.     

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.     

Effect of T-R conformational change on sickle-cell hemoglobin interactions and aggregation

We compare the role of a conformational switch and that of a point mutation in the thermodynamic stability of a protein solution and in the consequent propensity toward aggregation. We study sickle-cell hemoglobin (HbS), the beta6 Glu-Val point mutant of adult human hemoglobin (HbA), in its R (CO-liganded) conformation, and compare its aggregation properties to those of both HbS and HbA in their T (unliganded) conformation. Static and dynamic light scattering measurements performed for various hemoglobin concentrations showed critical divergences with mean field exponents as temperature was increased. This allowed determining spinodal data points T(S)(c) by extrapolation. These points were fitted to theoretical expressions of the T(S)(c) spinodal line, which delimits the region where the homogeneous solution becomes thermodynamically unstable against demixing in two sets of denser and dilute mesoscopic domains, while remaining still liquid. Fitting provided model-free numerical values of enthalpy and entropy parameters measuring the stability of solutions against demixing, namely, 93.2 kJ/mol and 314 J/ degrees K-mol, respectively. Aggregation was observed also for R-HbS, but in amorphous form and above physiological temperatures close to the spinodal, consistent with the role played in nucleation by anomalous fluctuations governed by the parameter epsilon = (T - T(S))/T(S). Fourier transform infrared (FTIR) and optical spectroscopy showed that aggregation is neither preceded nor followed by denaturation. Transient multiple interprotein contacts occur in the denser liquid domains for R-HbS, T-HbS, and T-HbA. The distinct effects of their specific nature and configurations, and those of desolvation on the demixing and aggregation thermodynamics, and on the aggregate structure are highlighted.     

Time scale of protein aggregation dictated by liquid-liquid demixing

The growing impact of protein aggregation pathologies, together with the current high need for extensive information on protein structures are focusing much interest on the physics underlying the nucleation and growth of protein aggregates and crystals. Sickle Cell Hemoglobin (HbS), a point-mutant form of normal human Hemoglobin (HbA), is the first recognized and best-studied case of pathologically aggregating protein. Here we reanalyze kinetic data on nucleation of deoxy-HbS aggregates by referring them to the (concentration-dependent) temperature T(s) characterizing the occurrence of the phase transition of liquid-liquid demixing (LLD) of the solution. In this way, and by appropriate scaling of kinetic data at different concentrations, so as to normalize their spans, the apparently disparate sets of data are seen to fall on a master curve. Expressing the master curve vs. the parameter epsilon = (T - T(s)) / T(s), familiar from phase transition theory, allows eliciting the role of anomalously large concentration fluctuations associated with the LLD phase transition and also allows decoupling quantitatively the role of such fluctuations from that of microscopic, inter-protein interactions leading to nucleation. Referring to epsilon shows how in a narrow temperature span, that is at T - T(s), nucleation kinetics can undergo orders-of-magnitude changes, unexpected in terms of ordinary chemical kinetics. The same is true for similarly small changes of other parameters (pH, salts, precipitants), capable of altering T(s) and consequently epsilon. This offers the rationale for understanding how apparently minor changes of parameters can dramatically affect protein aggregation and related diseases.     

Vapour-to-liquid nucleation in cone pores

In this work, we perform a systematic study on the vapour-to-liquid nucleation in a cone pore by using comparatively the classical nucleation theory (CNT) and the constrained lattice density functional theory (constrained LDFT). Three different nucleation scenarios are observed depending on the contact angle θ and apex angle α: the spontaneous phase transition, nucleation with a critical nucleus in Cassie state and nucleation with a critical nucleus in Wenzel state. We also found that both the diagram for nucleation mechanisms and the reduced nucleation barriers with respect to the homogeneous nucleation given by the CNT are at least qualitatively consistent with those from the constrained LDFT. For an increasingly small critical nucleus, the difference between nucleation behaviours from two methods becomes significant due to the macroscopic approximations embedded in CNT, which fails to describe the curvature dependence of surface tension, the line tension effect and the space confinement effect.     

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.     

Charge centers and formation of the protein folding core

Electrostatic interactions are important for protein folding. At low resolution, the electrostatic field of the whole molecule can be described in terms of charge center(s). To study electrostatic effects, the centers of positive and negative charge were calculated for 20 small proteins of known structure, for which hydrogen exchange cores had been determined experimentally. Two observations seem to be important. First, in all 20 proteins studied 30-100% of the residues forming hydrogen exchange core(s) were clustered around the charge centers. Moreover, in each protein more than half of the core sequences are located near the centers of charge. Second, the general architecture of all-alpha proteins from the set seems to be stabilized by interactions of residues surrounding the charge centers. In most of the alpha-beta proteins, either or both of the centers are located near a pair of consecutive strands, and this is even more characteristic for alpha/Beta and all-beta structures. Consecutive strands are very probable sites of early folding events. These two points lead to the conclusion that charge centers, defined solely from the structure of the folded protein may indicate the location of a protein's hydrogen exchange/folding core. In addition, almost all the proteins contain well-conserved continuous hydrophobic sequences of three or more residues located in the vicinity of the charge centers. These hydrophobic sequences may be primary nucleation sites for protein folding. The results suggest the following scheme for the order of events in folding: local hydrophobic nucleation, electrostatic collapse of the core, global hydrophobic collapse, and slow annealing to the native state. This analysis emphasizes the importance of treating electrostatics during protein-folding simulations.     

A model for the sickle hemoglobin fiber using both mutation sites

The standard molecular model of the fiber of the sickle hemoglobin (HbS: beta6 Glu-->Val) has been revised to allow both beta6 mutation sites to participate in intermolecular contacts, rather than only one beta6 site as previously thought, for four molecules per 14-molecule fiber cross section. This structure accurately predicts the copolymerization of hybridized mixtures of HbS with HbA or HbC (beta6 Glu-->Lys), which could not be reconciled with prior models in which only half the beta6 sites were required for assembly. This model suggests new contacts within the fiber and raises the question of whether these cross-linked double strands could possess added stability important in such processes as nucleation.     

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.     

Polymerization of deoxy-sickle cell hemoglobin in high-phosphate buffer

Deoxy-sicklecell hemoglobin (HbS) polymerizes in 0.05 M phosphate buffer to form long helical fibers. The reaction typically occurs when the concentration of HbS is about 165 mg/ml. Polymerization produces a variety of polymorphic forms. The structure of the fibers can be probed by using site-directed mutants to examine the effect of altering the residues involved in intermolecular interactions. Polymerization can also be induced in the presence of 1.5 M phosphate buffer. Under these conditions polymerization occurs at much lower concentrations (ca. 5 mg/ml), which is advantageous when site-directed mutants are being used because only small quantities of the mutants are available. We have characterized the structure of HbS polymers formed in 1.5 M phosphate to determine how their structures are related to the polymers formed under more physiological conditions. Under both sets of conditions fibers are the first species to form. At pHs between 6.7 and 7.3 fibers initially form bundles and then crystals. At lower pHs fibers form macrofibers and then crystals. Fourier transforms of micrographs of the polymers formed in 1.5 M phosphate display the 32- and 64-A(-1) periodicity characteristic of fibers formed in 0.05 M phosphate buffer. The 64-A(-1) layer line is less prominent in Fourier transforms of negatively stained fibers formed in 1.5 M phosphate possibly because salt interferes with staining of the fibers. However, micrographs and Fourier transforms of frozen hydrated fibers formed in high and low phosphate display the same periodicities. Under both sets of reaction conditions HbS polymers form crystals with the same unit cell parameters as Wishner-Love crystals (a = 64 A, b = 185 A, c = 53 A). Some of the polymerization intermediates were examined in the frozen-hydrated state in order to determine whether their structures were significantly perturbed by negative staining. We have also carried out reconstructions of the frozen-hydrated fibers in high and low phosphate to compare their molecular coordinates. The helical projection of the reconstructions in low phosphate shows the expected 14-strand structure. In high phosphate the 14-strand fibers are also formed and their molecular coordinates are the same (within experimental error) as those of fibers formed in 0.05 M phosphate. In addition, the reconstructions of high-phosphate fibers reveal a new minor variant of fiber containing 10 strands. The polymerization products in 1.5 M phosphate buffer were generally indistinguishable from those formed in 0.05 M phosphate buffer. Micrographs of frozen hydrated specimens have facilitated the interpretation of previously published micrographs using negative staining.     

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.     

Metastable mesoscopic clusters in solutions of sickle-cell hemoglobin

Sickle cell hemoglobin (HbS) is a mutant, whose polymerization while in deoxy state in the venous circulation underlies the debilitating sickle cell anemia. It has been suggested that the nucleation of the HbS polymers occurs within clusters of dense liquid, existing in HbS solutions. We use dynamic light scattering with solutions of deoxy-HbS, and, for comparison, of oxy-HbS and oxy-normal adult hemoglobin, HbA. We show that solutions of all three Hb variants contain clusters of dense liquid, several hundred nanometers in size, which are metastable with respect to the Hb solutions. The clusters form within a few seconds after solution preparation and their sizes and numbers remain relatively steady for up to 3 h. The lower bound of the cluster lifetime is 15 ms. The clusters exist in broad temperature and Hb concentration ranges, and occupy 10(-5)-10(-2) of the solution volume. The results on the cluster properties can serve as test data for a potential future microscopic theory of cluster stability and kinetics. More importantly, if the clusters are a part of the nucleation mechanism of HbS polymers, the rate of HbS polymerization can be controlled by varying the cluster properties.     

物理环境影响蛋白质晶体形核的研究进展

物理环境是影响蛋白质晶体形核的重要因素。文中回顾了各种物理环境如光、电场、超声波、磁场、微重力、温度、机械振动、异相形核界面对蛋白质晶体形核的影响,并对各物理环境下蛋白质晶体形核的可能机制进行探讨,展望了利用物理环境影响蛋白质晶体形核的研究前景。     

Heterogeneous nucleation in sickle hemoglobin: experimental validation of a structural mechanism

Sickle hemoglobin polymerizes by two types of nucleation: homogeneous nucleation of aggregates in solution, and heterogeneous nucleation on preexisting polymers. It has been proposed that the same contact that is made in the interior of the polymer between the mutant site beta6 and its receptor pocket on an adjacent molecule is the primary contact site for the heterogeneous nucleus. We have constructed cross-linked hybrid molecules in which one beta-subunit is from HbA with Glu at beta6, and the other is from HbS with a Val at beta6. We measured solubility (using sedimentation) and polymerization kinetics (using laser photolysis) on cross-linked hybrids, and cross-linked HbS as controls. We find approximately 4000 times less heterogeneous nucleation in the cross-linked AS molecules than in cross-linked HbS, in strong confirmation of the proposal. In addition, changes in stability of the nucleus support a further proposal that more than one beta6 contact is involved in the homogeneous nucleus.     

Microtubule nucleation: The waltz between γ-tubulin ring complex and associated proteins

Microtubules are essential cytoskeletal elements assembled from αβ-tubulin dimers. In high eukaryotes, microtubule nucleation, the de novo assembly of a microtubule from its minus end, is initiated by the γ-tubulin ring complex (γ-TuRC). Despite many years of research, the structural and mechanistic principles of the microtubule nucleation machinery remained poorly understood. Only recently, cryoelectron microscopy studies uncovered the molecular organization and potential activation mechanisms of γ-TuRC. In vitro assays further deciphered the spatial and temporal cooperation between γ-TuRC and additional factors, for example, the augmin complex, the phase separation protein TPX2, and the microtubule polymerase XMAP215. These breakthroughs deepen our understanding of microtubule nucleation mechanisms and will link the assembly of individual microtubules to the organization of cellular microtubule networks.