Coarse‐grained and all‐atom modeling of structural states and transitions in hemoglobin

Hemoglobin (Hb), an oxygen‐binding protein composed of four subunits (α1, α2, β1, and β2), is a well‐known example of allosteric proteins that are capable of cooperative ligand binding. Despite decades of studies, the structural basis of its cooperativity remains controversial. In this study, we have integrated coarse‐grained (CG) modeling, all‐atom simulation, and structural data from X‐ray crystallography and wide‐angle X‐ray scattering (WAXS), aiming to probe dynamic properties of the two structural states of Hb (T and R state) and the transitions between them. First, by analyzing the WAXS data of unliganded and liganded Hb, we have found that the structural ensemble of T or R state is dominated by one crystal structure of Hb with small contributions from other crystal structures of Hb. Second, we have used normal mode analysis to identify two distinct quaternary rotations between the α1β1 and α2β2 dimer, which drive the transitions between T and R state. We have also identified the hot‐spot residues whose mutations are predicted to greatly change these quaternary motions. Third, we have generated a CG transition pathway between T and R state, which predicts a clear order of quaternary and tertiary changes involving α and β subunits in Hb. Fourth, we have used the accelerated molecular dynamics to perform an all‐atom simulation starting from the T state of Hb, and we have observed a transition toward the R state of Hb. Further analysis of crystal structural data and the all‐atom simulation trajectory has corroborated the order of quaternary and tertiary changes predicted by CG modeling. Proteins 2013. © 2012 Wiley Periodicals, Inc.     

A genetic adaptive pattern＊low hemoglobin concentration in the himalayan highlanders

Mean hemoglobin （Hb） concentration of about 3 500 subjects derived from 17 studies of Himalayan highlanders （Tibetans, Sherpas, and Ladakhis） was compared with lowlanders （Chinese Han, Indian Tamils） lived in the Himalayas, and European climbers during Everest expeditions as well as Andean natives. The results found that Hb concentration in Himalayan highlanders was systemically lower than those reported for Andean natives and lowland immigrants. These comparative data demonstrated that a healthy native population may successfully reside at high altitude without a significant elevation in Hb, and the lower Hb levels of Himalayan highlanders than those of migrated lowlanders and Andean natives are an example of favourable adaptation over the generations. In addition, excessive polycythemia has frequently been used as a marker of chronic mountain sickness （CMS）. Altitude populations who have a higher Hb concentration also have a higher incidence of CMS. The low Hb in Himalayans suggested as showing adaptation over many generations in Tibetan stock. Recent work in Tibet, suggested that Tibetans there may have adapted to high altitude as a result of evolutionary pressure selecting for genes which give an advantage at altitude. All of the population genomic and statistical analysis indicated that EPAS1 and EGLN1 are mostly likely responsible for high altitude adaptation and closely related to low Hb concentration in Tibetans. These data supported the hypothesis that Himalayan highlanders have evolved a genetically different erythropoietic response to chronic hypoxia by virtue of their much longer exposure to high altitude.     

Monoamine-dependent Production of Reactive Oxygen Species Catalyzed by Pseudoperoxidase Activity of Human Hemoglobin

Hemoglobin (Hb) solution-based blood substitutes are being developed as oxygen-carrying agents for the prevention of ischemic tissue damage and low blood volume-shock. However, the cell-free Hb molecule has intrinsic toxicity to the tissue since harmful reactive oxygen species (ROS) are readily produced during autoxidation of Hb from the ferrous state to the ferric state, and the cell-free Hb also causes distortion in the oxidant/antioxidant balance in the tissues. There may be further hindering dangers in the use of free Hb as a blood substitute. It has been reported that Hb has peroxidase-like activity oxidizing peroxidase substrates such as aromatic amines. Here we observed the Hb-catalyzed ROS production coupled to oxidation of a neurotransmitter precursor, β-phenylethylamine (PEA). Addition of PEA to Hb solution resulted in generation of superoxide anion (O₂^??). We also observed that PEA increases the Hb-catalyzed monovalent oxidation of ascorbate to ascorbate free radicals (Asc^?). The O₂^?? generation and Asc^? formation were detected by O₂^??-specific chemiluminescence of the Cypridina lucigenin analog and electron spin resonance spectroscopy, respectively. PEA-dependent O₂^?? production and monovalent oxidation of ascorbate in the Hb solution occurred without addition of H₂O₂, but a trace of H₂O₂ added to the system greatly increased the production of both O₂^?? and Asc^?. Addition of GSH completely inhibited the PEA-dependent production of O₂^?? and Asc^? in Hb solution. We propose that the O₂^?? generation and Asc^? formation in the Hb solution are due to the pseudoperoxidase activity-dependent oxidation of PEA and resultant ROS may damage tissues rich in monoamines, if the Hb-based blood substitutes were circulated without addition of ROS scavengers such as thiols.     

Construction of novel chimeric proteins through the truncation of SEC2 and Sak from Staphylococcus aureus

It is an usual clinical phenomenon that cancer patients are prone to thrombosis. Until now, there have been no efficient methods or appropriate drugs to prevent and cure tumor thrombus. Therefore, the construction of a bifunctional chimeric protein for the treatment of cancer, complicated with thrombosis, is of great significance. Utilizing the superantigenic activity of staphylococcal enterotoxin C2 (SEC2) and the thrombolytic activity of staphylokinase (Sak), Sak-linker-SEC2 and SEC2-linker-Sak were constructed which had good anti-tumor and thrombolytic activities at the same time. Due to the intrinsic emetic activity of SEC2 and high molecular weight (MW) of chimeric proteins (44?kDa), their clinical applications will be restricted. In this study, novel chimeric proteins including ΔSEC2–ΔSak and ΔSak–ΔSEC2 were constructed through the truncation of SEC2 and Sak without 9-Ala linker and His-tag. Compared with the former, both the truncated proteins preserved nearly the same anti-tumor and thrombolytic activities. In addition, their MWs were only 29?kDa and their immunoreactivities were slightly lower than that of Sak-linker-SEC2 and SEC2-linker-Sak, respectively. Therefore, the novel chimeric proteins possessed merits and characteristics, such as low MS, low immunogenicity, and difunctionality which the former had not. It will be of great interest if the above-mentioned proteins can be used to cure Trousseau syndrome in clinic.     

Midpoint attractors and species richness: Modelling the interaction between environmental drivers and geometric constraints

We introduce a novel framework for conceptualising, quantifying and unifying discordant patterns of species richness along geographical gradients. While not itself explicitly mechanistic, this approach offers a path towards understanding mechanisms. In this study, we focused on the diverse patterns of species richness on mountainsides. We conjectured that elevational range midpoints of species may be drawn towards a single midpoint attractor – a unimodal gradient of environmental favourability. The midpoint attractor interacts with geometric constraints imposed by sea level and the mountaintop to produce taxon‐specific patterns of species richness. We developed a Bayesian simulation model to estimate the location and strength of the midpoint attractor from species occurrence data sampled along mountainsides. We also constructed midpoint predictor models to test whether environmental variables could directly account for the observed patterns of species range midpoints. We challenged these models with 16 elevational data sets, comprising 4500 species of insects, vertebrates and plants. The midpoint predictor models generally failed to predict the pattern of species midpoints. In contrast, the midpoint attractor model closely reproduced empirical spatial patterns of species richness and range midpoints. Gradients of environmental favourability, subject to geometric constraints, may parsimoniously account for elevational and other patterns of species richness.     

GbEXPATR,a species‐specific expansin,enhances cotton fibre elongation through cell wall restructuring

Cotton provides us the most important natural fibre. High fibre quality is the major goal of cotton breeding, and introducing genes conferring longer, finer and stronger fibre from Gossypium barbadense to Gossypium hirsutum is an important breeding strategy. We previously analysed the G. barbadense fibre development mechanism by gene expression profiling and found two homoeologous fibre‐specific α‐expansins from G. barbadense, GbEXPA2 and GbEXPATR. GbEXPA2 (from the D_T genome) is a classical α‐expansin, while its homoeolog, GbEXPATR (A_T genome), encodes a truncated protein lacking the normal C‐terminal polysaccharide‐binding domain of other α‐expansins and is specifically expressed in G. barbadense. Silencing EXPA in G. hirsutum induced shorter fibres with thicker cell walls. GbEXPA2 overexpression in G. hirsutum had no effect on mature fibre length, but produced fibres with a slightly thicker wall and increased crystalline cellulose content. Interestingly, GbEXPATR overexpression resulted in longer, finer and stronger fibres coupled with significantly thinner cell walls. The longer and thinner fibre was associated with lower expression of a number of secondary wall‐associated genes, especially chitinase‐like genes, and walls with lower cellulose levels but higher noncellulosic polysaccharides which advocated that a delay in the transition to secondary wall synthesis might be responsible for better fibre. In conclusion, we propose that α‐expansins play a critical role in fibre development by loosening the cell wall; furthermore, a truncated form, GbEXPATR, has a more dramatic effect through reorganizing secondary wall synthesis and metabolism and should be a candidate gene for developing G. hirsutum cultivars with superior fibre quality.     

Heme Orientation of Cavity Mutant Hemoglobins (His F8 → Gly) in Either α or β Subunits: Circular Dichroism, 1H NMR,and Resonance Raman Studies

Native human adult hemoglobin (Hb A) has mostly normal orientation of heme, whereas recombinant Hb A (rHb A) expressed in E. coli contains both normal and reversed orientations of heme. Hb A with the normal heme exhibits positive circular dichroism (CD) bands at both the Soret and 260‐nm regions, while rHb A with the reversed heme shows a negative Soret and decreased 260‐nm CD bands. In order to examine involvement of the proximal histidine (His F8) of either α or β subunits in determining the heme orientation, we prepared two cavity mutant Hbs, rHb(αH87G) and rHb(βH92G), with substitution of glycine for His F8 in the presence of imidazole. CD spectra of both cavity mutant Hbs did not show a negative Soret band, but instead exhibited positive bands with strong intensity at the both Soret and 260‐nm regions, suggesting that the reversed heme scarcely exists in the cavity mutant Hbs. We confirmed by ¹H NMR and resonance Raman (RR) spectroscopies that the cavity mutant Hbs have mainly the normal heme orientation in both the mutated and native subunits. These results indicate that the heme Fe‐His F8 linkage in both α and β subunits influences the heme orientation, and that the heme orientation of one type of subunit is related to the heme orientation of the complementary subunits to be the same. The present study showed that CD and RR spectroscopies also provided powerful tools for the examination of the heme rotational disorder of Hb A, in addition to the usual ¹H NMR technique. Chirality 28:585–592, 2016. © 2016 Wiley Periodicals, Inc.     

Truncated O‐glycans promote epithelial‐to‐mesenchymal transition and stemness properties of pancreatic cancer cells

Aberrant expression of Sialyl‐Tn (STn) antigen correlates with poor prognosis and reduced patient survival. We demonstrated that expression of Tn and STn in pancreatic ductal adenocarcinoma (PDAC) is due to hypermethylation of Co re 1 s ynthase specific m olecular c haperone (COSMC) and enhanced the malignant properties of PDAC cells with an unknown mechanism. To explore the mechanism, we have genetically deleted COSMC in PDAC cells to express truncated O‐glycans (SimpleCells, SC) which enhanced cell migration and invasion. Since epithelial‐to‐mesenchymal transition (EMT) play a vital role in metastasis, we have analysed the induction of EMT in SC cells. Expressions of the mesenchymal markers were significantly high in SC cells as compared to WT cells. Equally, we found reduced expressions of the epithelial markers in SC cells. Re‐expression of COSMC in SC cells reversed the induction of EMT. In addition to this, we also observed an increased cancer stem cell population in SC cells. Furthermore, orthotopic implantation of T3M4 SC cells into athymic nude mice resulted in significantly larger tumours and reduced animal survival. Altogether, these results suggest that aberrant expression of truncated O‐glycans in PDAC cells enhances the tumour aggressiveness through the induction of EMT and stemness properties.     

The aspartimide problem persists: Fluorenylmethyloxycarbonyl‐solid‐phase peptide synthesis (Fmoc‐SPPS) chain termination due to formation of N‐terminal piperazine‐2,5‐diones

Aspartimide (Asi) formation is a notorious side reaction in peptide synthesis that is well characterized and described in literature. In this context, we observed significant amounts of chain termination in Fmoc‐SPPS while synthesizing the N‐terminal Xaa‐Asp‐Yaa motif. This termination was caused by the formation of piperazine‐2,5‐diones. We investigated this side reaction using a linear model peptide and independently synthesizing its piperazine‐2,5‐dione derivative. Nuclear magnetic resonance (NMR) data of the side product present in the crude linear peptide proves that exclusively the six‐membered ring is formed whereas the theoretically conceivable seven‐membered 1,4‐diazepine‐2,5‐dione is not found. We propose a mechanism where nucleophilic attack of the N‐terminal amino function takes place at the α‐carbon of the carbonyl group of the corresponding Asi intermediate. In addition, we systematically investigated the impact of (a) different adjacent amino acid residues, (b) backbone protection, and (c) side chain protection of flanking amino acids. The side reaction is directly related to the Asi intermediate. Hence, hindering or avoiding Asi formation reduces or completely suppresses this side reaction.