Bridging the gap between chemistry, physiology, and evolution: Quantifying the functionality of sperm whale myoglobin mutants

Controlled oxygen delivery to power tissue regeneration

Oxygen plays a crucial role in human embryogenesis, homeostasis, and tissue regeneration. Emerging engineered regenerative solutions call for novel oxygen delivery systems. To become a reality, these systems must consider physiological processes, oxygen release mechanisms and the target application. In this review, we explore the biological relevance of oxygen at both a cellular and tissue level, and the importance of its controlled delivery via engineered biomaterials and devices. Recent advances and upcoming trends in the field are also discussed with a focus on tissue-engineered constructs that could meet metabolic demands to facilitate regeneration.

Hydrogel oxygen reservoirs increase functional integration of neural stem cell grafts by meeting metabolic demands

Injectable biomimetic hydrogels have great potential for use in regenerative medicine as cellular delivery vectors. However, they can suffer from issues relating to hypoxia, including poor cell survival, differentiation, and functional integration owing to the lack of an established vascular network. Here we engineer a hybrid myoglobin:peptide hydrogel that can concomitantly deliver stem cells and oxygen to the brain to support engraftment until vascularisation can occur naturally. We show that this hybrid hydrogel can modulate cell fate specification within progenitor cell grafts, resulting in a significant increase in neuronal differentiation. We find that the addition of myoglobin to the hydrogel results in more extensive innervation within the host tissue from the grafted cells, which is essential for neuronal replacement strategies to ensure functional synaptic connectivity. This approach could result in greater functional integration of stem cell-derived grafts for the treatment of neural injuries and diseases affecting the central and peripheral nervous systems.

Common and unique strategies of myoglobin evolution for deep-sea adaptation of diving mammals

Myoglobin (Mb) is highly concentrated in the myocytes of diving mammals such as whales and seals, in comparison with land animals, and its molecular evolution has played a crucial role in their deep-sea adaptation. We previously resurrected ancestral whale Mbs and demonstrated the evolutional strategies for higher solubility under macromolecular crowding conditions. Pinnipeds, such as seals and sea lions, are also expert diving mammals with Mb-rich muscles. In the present study, we resurrected ancestral pinniped Mbs and investigated their biochemical and structural properties. Comparisons between pinniped and whale Mbs revealed the common and distinctive strategies for the deep-sea adaptation. The overall evolution processes, gaining precipitant tolerance and improving thermodynamic stability, were commonly observed. However, the strategies for improving the folding stability differed, and the pinniped Mbs exploited the shielding of hydrophobic surfaces more effectively than the whale Mbs.

Adaptations to environmental change: Globin superfamily evolution in Antarctic fishes

The ancient origins and functional versatility of globins make them ideal subjects for studying physiological adaptation to environmental change. Our goals in this review are to describe the evolution of the vertebrate globin gene superfamily and to explore the structure/function relationships of hemoglobin, myoglobin, neuroglobin and cytoglobin in teleost fishes. We focus on the globins of Antarctic notothenioids, emphasizing their adaptive features as inferred from comparisons with human proteins. We dedicate this review to Guido di Prisco, our co-author, colleague, friend, and husband of C.V. Ever thoughtful, creative, and enthusiastic, Guido spearheaded study of the structure, function, and evolution of the hemoglobins of polar fishes - this review is testimony to his wide-ranging contributions. Throughout his career, Guido inspired younger scientists to embrace polar biological research, and he challenged researchers of all ages to explore evolutionary adaptation in the context of global climate change. Beyond his scientific contributions, we will miss his warmth, his culture, and his great intellect. Guido has left an outstanding legacy, one that will continue to inspire us and our research.

Myoglobin oxygen affinity in aquatic and terrestrial birds and mammals

Myoglobin (Mb) is an oxygen binding protein found in vertebrate skeletal muscle, where it facilitates intracellular transport and storage of oxygen. This protein has evolved to suit unique physiological needs in the muscle of diving vertebrates that express Mb at much greater concentrations than their terrestrial counterparts. In this study, we characterized Mb oxygen affinity (P50) from 25 species of aquatic and terrestrial birds and mammals. Among diving species, we tested for correlations between Mb P50 and routine dive duration. Across all species examined, Mb P50 ranged from 2.40 to 4.85 mmHg. The mean P50 of Mb from terrestrial ungulates was 3.72±0.15 mmHg (range 3.70-3.74 mmHg). The P50 of cetaceans was similar to terrestrial ungulates ranging from 3.54 to 3.82 mmHg, with the exception of the melon-headed whale, which had a significantly higher P50 of 4.85 mmHg. Among pinnipeds, the P50 ranged from 3.23 to 3.81 mmHg and showed a trend for higher oxygen affinity in species with longer dive durations. Among diving birds, the P50 ranged from 2.40 to 3.36 mmHg and also showed a trend of higher affinities in species with longer dive durations. In pinnipeds and birds, low Mb P50 was associated with species whose muscles are metabolically active under hypoxic conditions associated with aerobic dives. Given the broad range of potential globin oxygen affinities, Mb P50 from diverse vertebrate species appears constrained within a relatively narrow range. High Mb oxygen affinity within this range may be adaptive for some vertebrates that make prolonged dives.

The influence of selection for protein stability on dN/dS estimations

Understanding the relative contributions of various evolutionary processes-purifying selection, neutral drift, and adaptation-is fundamental to evolutionary biology. A common metric to distinguish these processes is the ratio of nonsynonymous to synonymous substitutions (i.e., dN/dS) interpreted from the neutral theory as a null model. However, from biophysical considerations, mutations have non-negligible effects on the biophysical properties of proteins such as folding stability. In this work, we investigated how stability affects the rate of protein evolution in phylogenetic trees by using simulations that combine explicit protein sequences with associated stability changes. We first simulated myoglobin evolution in phylogenetic trees with a biophysically realistic approach that accounts for 3D structural information and estimates of changes in stability upon mutation. We then compared evolutionary rates inferred directly from simulation to those estimated using maximum-likelihood (ML) methods. We found that the dN/dS estimated by ML methods (ωML) is highly predictive of the per gene dN/dS inferred from the simulated phylogenetic trees. This agreement is strong in the regime of high stability where protein evolution is neutral. At low folding stabilities and under mutation-selection balance, we observe deviations from neutrality (per gene dN/dS > 1 and dN/dS < 1). We showed that although per gene dN/dS is robust to these deviations, ML tests for positive selection detect statistically significant per site dN/dS > 1. Altogether, we show how protein biophysics affects the dN/dS estimations and its subsequent interpretation. These results are important for improving the current approaches for detecting positive selection.

Computing stability effects of mutations in human superoxide dismutase 1

Protein stability is affected in several diseases and is of substantial interest in efforts to correlate genotypes to phenotypes. Superoxide dismutase 1 (SOD1) is a suitable test case for such correlations due to its abundance, stability, available crystal structures and thermochemical data, and physiological importance. In this work, stability changes of SOD1 mutations were computed with five methods, CUPSAT, I-Mutant2.0, I-Mutant3.0, PoPMuSiC, and SDM, with emphasis on structural sensitivity as a potential issue in structure-based protein calculation. The large correlation between experimental literature data of SOD1 dimers and monomers (r = 0.82) suggests that mutations in separate protein monomers are mostly additive. PoPMuSiC was most accurate (typical MAE ~ 1 kcal/mol, r ~ 0.5). The relative performance of the methods was not very structure-dependent, and the more accurate methods also displayed less structural sensitivity, with the standard deviation from different high-resolution structures down to ~0.2 kcal/mol. Structures of variable resolution and number of protein copies locally affected specific sites, emphasizing the use of state-relevant crystal structures when such sites are of interest, but had little impact on overall batch estimates. Protein-interaction effects (as a mimic of crystal packing) were small for the more accurate methods. Thus, batch computations, relevant to, e.g., comparisons of disease/nondisease mutant sets or different clades in phylogenetic trees, are much more significant than single mutant calculations and may be the only meaningful way to computationally bridge the genotype-phenotype gap of proteomics. Finally, mutations involving glycine were most difficult to model, of relevance to future method improvement. This could be due to structure changes (glycine has a low structural propensity) or water colocalization with glycine.

Unfolding simulations of holomyoglobin from four mammals: identification of intermediates and β-sheet formation from partially unfolded states

Myoglobin (Mb) is a centrally important, widely studied mammalian protein. While much work has investigated multi-step unfolding of apoMb using acid or denaturant, holomyoglobin unfolding is poorly understood despite its biological relevance. We present here the first systematic unfolding simulations of holoMb and the first comparative study of unfolding of protein orthologs from different species (sperm whale, pig, horse, and harbor seal). We also provide new interpretations of experimental mean molecular ellipticities of myoglobin intermediates, notably correcting for random coil and number of helices in intermediates. The simulated holoproteins at 310 K displayed structures and dynamics in agreement with crystal structures (R g ~1.48-1.51 nm, helicity ~75%). At 400 K, heme was not lost, but some helix loss was observed in pig and horse, suggesting that these helices are less stable in terrestrial species. At 500 K, heme was lost within 1.0-3.7 ns. All four proteins displayed exponentially decaying helix structure within 20 ns. The C- and F-helices were lost quickly in all cases. Heme delayed helix loss, and sperm whale myoglobin exhibited highest retention of heme and D/E helices. Persistence of conformation (RMSD), secondary structure, and ellipticity between 2-11 ns was interpreted as intermediates of holoMb unfolding in all four species. The intermediates resemble those of apoMb notably in A and H helices, but differ substantially in the D-, E- and F-helices, which interact with heme. The identified mechanisms cast light on the role of metal/cofactor in poorly understood holoMb unfolding. We also observed β-sheet formation of several myoglobins at 500 K as seen experimentally, occurring after disruption of helices to a partially unfolded, globally disordered state; heme reduced this tendency and sperm-whale did not display any sheet propensity during the simulations.

Positively selected sites in cetacean myoglobins contribute to protein stability

Since divergence ∼50 Ma ago from their terrestrial ancestors, cetaceans underwent a series of adaptations such as a ∼10-20 fold increase in myoglobin (Mb) concentration in skeletal muscle, critical for increasing oxygen storage capacity and prolonging dive time. Whereas the O2-binding affinity of Mbs is not significantly different among mammals (with typical oxygenation constants of ∼0.8-1.2 µM(-1)), folding stabilities of cetacean Mbs are ∼2-4 kcal/mol higher than for terrestrial Mbs. Using ancestral sequence reconstruction, maximum likelihood and bayesian tests to describe the evolution of cetacean Mbs, and experimentally calibrated computation of stability effects of mutations, we observe accelerated evolution in cetaceans and identify seven positively selected sites in Mb. Overall, these sites contribute to Mb stabilization with a conditional probability of 0.8. We observe a correlation between Mb folding stability and protein abundance, suggesting that a selection pressure for stability acts proportionally to higher expression. We also identify a major divergence event leading to the common ancestor of whales, during which major stabilization occurred. Most of the positively selected sites that occur later act against other destabilizing mutations to maintain stability across the clade, except for the shallow divers, where late stability relaxation occurs, probably due to the shorter aerobic dive limits of these species. The three main positively selected sites 66, 5, and 35 undergo changes that favor hydrophobic folding, structural integrity, and intra-helical hydrogen bonds.

Functional properties of myoglobins from five whale species with different diving capacities

Whales show an exceptionally wide range of diving capabilities and many express high amounts of the O(2) carrier protein myoglobin (Mb) in their muscle tissues, which increases their aerobic diving capacity. Although previous studies have mainly focused on the muscle Mb concentration and O(2) carrying capacity as markers of diving behavior in whales, it still remains unexplored whether whale Mbs differ in their O(2) affinities and nitrite reductase and peroxidase enzymatic activities, all functions that could contribute to differences in diving capacities. In this study, we have measured the functional properties of purified Mbs from five toothed whales and two baleen whales and have examined their correlation with average dive duration. Results showed that some variation in functional properties exists among whale Mbs, with toothed whale Mbs having higher O(2) affinities and nitrite reductase activities (similar to those of horse Mb) compared with baleen whale Mbs. However, these differences did not correlate with average dive duration. Instead, a significant correlation was found between whale Mb concentration and average duration and depth of dives, and between O(2) affinity and nitrite reductase activity when including horse Mb. Despite the fact that the functional properties showed little species-specific differences in vitro, they may still contribute to enhancing diving capacity as a result of the increased muscle Mb concentration found in extreme divers. In conclusion, Mb concentration rather than specific functional reactivities may support whale diving performance.