Structural Evidence of a Fourth Gla Residue in Fish Osteocalcin: Biological Implications,

The Landscape of Osteocalcin Proteoforms Reveals Distinct Structural and Functional Roles of Its Carboxylation Sites

Human osteocalcin (OC) undergoes reversible, vitamin K-dependent γ-carboxylation at three glutamic acid residues, modulating its release from bones and its hormonal roles. A complete understanding of OC roles and structure-activity relationships is still lacking, as only uncarboxylated and few differently carboxylated variants have been considered so far. To fill this lack of knowledge, a comprehensive experimental and computational investigation of the structural properties and calcium-binding activity of all the OC variants is reported here. Such a comparative study indicates that the carboxylation sites are not equivalent and differently affect the OC structure and interaction with calcium, properties that are relevant for the modulation of OC functions. This study also discloses cooperative effects and provides structural and mechanistic interpretation. The disclosed peculiar features of each carboxylated proteoform strongly suggest that considering all eight possible OC variants in future studies may help rationalize some of the conflicting hypotheses observed in the literature.

Recent advances in the crosstalk between adipose, muscle and bone tissues in fish

Control of tissue metabolism and growth involves interactions between organs, tissues, and cell types, mediated by cytokines or direct communication through cellular exchanges. Indeed, over the past decades, many peptides produced by adipose tissue, skeletal muscle and bone named adipokines, myokines and osteokines respectively, have been identified in mammals playing key roles in organ/tissue development and function. Some of them are released into the circulation acting as classical hormones, but they can also act locally showing autocrine/paracrine effects. In recent years, some of these cytokines have been identified in fish models of biomedical or agronomic interest. In this review, we will present their state of the art focusing on local actions and inter-tissue effects. Adipokines reported in fish adipocytes include adiponectin and leptin among others. We will focus on their structure characteristics, gene expression, receptors, and effects, in the adipose tissue itself, mainly regulating cell differentiation and metabolism, but in muscle and bone as target tissues too. Moreover, lipid metabolites, named lipokines, can also act as signaling molecules regulating metabolic homeostasis. Regarding myokines, the best documented in fish are myostatin and the insulin-like growth factors. This review summarizes their characteristics at a molecular level, and describes both, autocrine effects and interactions with adipose tissue and bone. Nonetheless, our understanding of the functions and mechanisms of action of many of these cytokines is still largely incomplete in fish, especially concerning osteokines (i.e., osteocalcin), whose potential cross talking roles remain to be elucidated. Furthermore, by using selective breeding or genetic tools, the formation of a specific tissue can be altered, highlighting the consequences on other tissues, and allowing the identification of communication signals. The specific effects of identified cytokines validated through in vitro models or in vivo trials will be described. Moreover, future scientific fronts (i.e., exosomes) and tools (i.e., co-cultures, organoids) for a better understanding of inter-organ crosstalk in fish will also be presented. As a final consideration, further identification of molecules involved in inter-tissue communication will open new avenues of knowledge in the control of fish homeostasis, as well as possible strategies to be applied in aquaculture or biomedicine.

Atomic Details of Biomineralization Proteins Inspiring Protein Design and Reengineering for Functional Biominerals

Biominerals are extraordinary materials that provide organisms with a variety of functions to support life. The synthesis of biominerals and organization at the macroscopic level is a consequence of the interactions of these materials with proteins. The association of biominerals and proteins is very ancient and has sparked a wealth of research across biological, medical and material sciences. Calcium carbonate, hydroxyapatite, and silica represent widespread natural biominerals. The atomic details of the interface between macromolecules and these biominerals is very intriguing from a chemical perspective, considering the association of chemical entities that are structurally different. With this review I provide an overview of the available structural studies of biomineralization proteins, explored from the Protein Data Bank (wwPDB) archive and scientific literature, and of how these studies are inspiring the design and engineering of proteins able to synthesize novel biominerals. The progression of this review from classical template proteins to silica polymerization seeks to benefit researchers involved in various interdisciplinary aspects of a biomineralization project, who need background information and a quick update on advances in the field. Lessons learned from structural studies are exemplary and will guide new projects for the imaging of new hybrid biomineral/protein superstructures at the atomic level.

The carboxylation status of osteocalcin has important consequences for its structure and dynamics

BACKGROUND

The carboxylation status of Osteocalcin (Ocn) not only influences formation and structure in bones but also has important endocrine functions affecting energy metabolism and expenditure. In this study, the role of γ-carboxylation of the glutamate residues in the structure-dynamics-function relationship in Ocn is investigated.

METHODS

Three forms of Ocn, differentially carboxylated at the Glu-17, 21 and 24 residues, along with a mutated form of Ocn carrying Glu/Ala mutations, are modeled and simulated using molecular dynamics (MD) simulation in the presence of calcium ions.

RESULTS

Characterization of the global conformational dynamics of Ocn, described in terms of the orientational variations within its 3-helical domain, highlights large structural variations in the non-carboxylated osteocalcin (nOcn). The bi-carboxylated Ocn (bOcn) and tri-carboxylated (tOcn) species, in contrast, display relatively rigid tertiary structures, with the dynamics of most regions strongly correlated. Radial distribution functions calculated for both bOcn and tOcn show long-range ordering of the calcium ion distribution around the carboxylated glutamate (γGlu) residues, likely playing an important role in promoting stability of these Ocns. Additionally, the same calcium ions are observed to coordinate with neighboring γGlu, better shielding their negative charges and in turn stabilizing these systems more than do the singly coordinating calcium ions observed in the case of nOcn. bOcn is also found to exhibit a more helical C-terminal structure, that has been shown to activate its cellular receptor GPRC6A, highlighting the allosteric role of Ocn carboxylation in modulating the stability and binding potential of the active C-terminal.

CONCLUSIONS

The carboxylation status of Ocn as well and its calcium coordination appear to have a direct influence on Ocn structure and dynamics, possibly leading to the known differences in Ocn biological function.

GENERAL SIGNIFICANCE

Modification of Ocn sequence or its carboxylation state may provide the blueprint for developing high-affinity peptides targeting its cellular receptor GPRC6A, with therapeutic potential for treatment of metabolic disorders.

Computational Nanomechanics of Noncollagenous Interfibrillar Interface in Bone

The noncollagenous interfibrillar interface in bone provides the critical function of transferring loads among collagen fibrils and their bundles, with adhesive mechanisms at this site thus significantly contributing to the mechanical properties of bone. Motivated by the experimental observations and hypotheses, a computational study is presented to elucidate the critical roles of two major proteins at the nanoscale interfibrillar interface, that is, osteopontin (OPN) and osteocalcin (OC) in bone. This study reveals the extremely high interfacial toughness of the OPN/OC composite. The previously proposed hypothesis of sacrificial bonds in the extracellular organic matrix is tested, and the remarkable mechanical properties of the nanoscale bone interface are attributed to the collaborative interactions between the OPN and OC proteins.

Deformation Mechanisms of "Two-Part" Natural Adhesive in Bone Interfibrillar Nano-Interfaces

Noncollagenous proteins at nanoscale interfaces in bone are less than 2-3% of bone content by weight, while they contribute more than 30% to fracture toughness. Major gaps in quantitative understanding of noncollagenous proteins' role in the interfibrillar interfaces, largely because of the limitation of probing their nanoscale dimension, have resulted in ongoing controversies and several outstanding hypotheses on their role and function, arguably going back to centuries ago to the original work from Galileo. Our results from the first detailed computational model of the nano-interface in the bone reveal "synergistic" deformation mechanism of a "double-part" natural glue, that is, noncollagenous osteopontin and osteocalcin at the interfibrillar interface. Specifically, through strong anchoring and formation of dynamic binding sites on mineral nanoplatelets, the nano-interface can sustain a large nonlinear deformation with ductility approaching 5000%. This large deformation results in an outstanding specific energy to failure exceeding ∼350 J/g, which is larger than the most known tough materials (such as Kevlar, spider silk, and so forth.).

Evidence for Osteocalcin Binding and Activation of GPRC6A in β-Cells

The possibility that G protein-coupled receptor family C member A (GPRC6A) is the osteocalcin (Ocn)-sensing G protein-coupled receptor that directly regulates pancreatic β-cell functions is controversial. In the current study, we found that Ocn and an Ocn-derived C-terminal hexapeptide directly activate GPRC6A-dependent ERK signaling in vitro. Computational models probe the structural basis of Ocn binding to GPRC6A and predict that the C-terminal hexapeptide docks to the extracellular side of the transmembrane domain of GPRC6A. Consistent with the modeling, mutations in the computationally identified binding pocket of GPRC6A reduced Ocn and C-terminal hexapeptide activation of this receptor. In addition, selective deletion of Gprc6a in β-cells (Gprc6a(β)(-cell-cko)) by crossing Gprc6a(flox/flox) mice with Ins2-Cre mice resulted in reduced pancreatic weight, islet number, insulin protein content, and insulin message expression. Both islet size and β-cell proliferation were reduced in Gprc6a(β)(-cell-cko) compared with control mice. Gprc6a(β)(-cell-cko) exhibited abnormal glucose tolerance, but normal insulin sensitivity. Islets isolated from Gprc6a(β)(-cell-cko) mice showed reduced insulin simulation index in response to Ocn. These data establish the structural basis for Ocn direct activation of GPRC6A and confirm a role for GPRC6A in regulating β-cell proliferation and insulin secretion.

Matrix Gla protein and osteocalcin: from gene duplication to neofunctionalization

Osteocalcin (OC or bone Gla protein, BGP) and matrix Gla protein (MGP) are two members of the growing family of vitamin K-dependent (VKD) proteins. They were the first VKD proteins found not to be involved in coagulation and synthesized outside the liver. Both proteins were isolated from bone although it is now known that only OC is synthesized by bone cells under normal physiological conditions, but since both proteins can bind calcium and hydroxyapatite, they can also accumulate in bone. Both OC and MGP share similar structural features, both in terms of protein domains and gene organization. OC gene is likely to have appeared from MGP through a tandem gene duplication that occurred concomitantly with the appearance of the bony vertebrates. Despite their relatively close relationship and the fact that both can bind calcium and affect mineralization, their functions are not redundant and they also have other unrelated functions. Interestingly, these two proteins appear to have followed quite different evolutionary strategies in order to acquire novel functionalities, with OC following a gene duplication strategy while MGP variability was obtained mostly by the use of multiple promoters and alternative splicing, leading to proteins with additional functional characteristics and alternative gene regulatory pathways.

A novel anabolic agent: a simvastatin analogue without HMG-CoA reductase inhibitory activity

For the first time, structural information regarding the role of simvastatin in bone anabolism is described, and a bone-specific statin is introduced. Polyaspartate-conjugated simvastatin was synthesized by solid-phase synthesis with the assistance of microwave irradiation. It displays significant bone targeting and bone formation with less toxicity than simvastatin.

Teleost fish osteocalcin 1 and 2 share the ability to bind the calcium mineral phase

The occurrence of a second osteocalcin (OC2) has been reported in teleost fish, where it coexists with OC1 in some species. While it has been proposed that OC2 gene originated from OC1 through the fish whole-genome duplication event, little information is available on its molecular function and physiological role. The present study brings biological data supporting the presence of OC2 in the mineral phase of teleost fish bone and its association with the mineral phase together with OC1. The occurrence of OC2 forms with different levels of phosphorylation or γ-carboxylation, and with amino acid substitutions was observed. Comparative analysis of mature peptide sequences revealed the high conservation existing between OC1 and OC2, in particular within the core γ-carboxyglutamic acid domain, and suggests that both protein forms may have the same function, i.e., binding of calcium ions or hydroxyapatite crystals.

Development of an enzyme-linked immunosorbent assay for detection of chicken osteocalcin and its use in evaluation of perch effects on bone remodeling in caged White Leghorns

Osteocalcin (OC) is a sensitive biochemical marker for evaluating bone turnover in mammals. The role of avian OC is less clear because of the need for a chicken assay. Our objectives were to develop an assay using indirect competitive ELISA for detecting chicken serum OC and use the assay to examine the effects of perches on bone remodeling in caged hens. Anti-chicken OC polyclonal antibody was produced by immunization of rabbits with a recombinant OC from Escherichia coli. Chicken OC extracted from bone was used as a coated protein, and purified chicken OC was used for calibration. The limit of detection of the developed OC ELISA was 0.13 ng/mL. The intra- and interassay CV were <7 and <12%, respectively. The sensitivity of the developed OC ELISA was compared with a commercial Rat-Mid OC ELISA in laying hens housed in conventional cages with or without perches. Serum samples were collected from 71-wk-old White Leghorn hens subjected to 4 treatments. Treatment 1 was control chickens that never had access to perches during their life cycle. Treatment 2 chickens had perches during the pullet phase (0 to 16.9 wk of age), whereas treatment 3 chickens had perches only during the egg-laying phase of the life cycle (17 to 71 wk of age). Treatment 4 chickens always had access to perches (0 to 71 wk of age). Correlation between the 2 assays was 0.62 (P < 0.0001). Levels of serum OC using the developed chicken ELISA were higher than that detected using the Rat-Mid ELISA (P < 0.0001). Results from the chicken ELISA assay showed that hens with perch access had higher concentrations of serum OC than hens without perches during egg laying (P = 0.04). Pullet access to perches did not affect serum OC levels in 71-wk-old hens (P = 0.15). In conclusion, a chicken OC ELISA has been validated that is sensitive and accurate with adequate discriminatory power for measuring bone remodeling in chickens.

Energy metabolism and the skeleton: Reciprocal interplay

The relation between bone remodelling and energy expenditure is an intriguing, and yet unexplained, challenge of the past ten years. In fact, it was only in the last few years that the skeleton was found to function, not only in its obvious roles of body support and protection, but also as an important part of the endocrine system. In particular, bone produces different hormones, like osteocalcin (OC), which influences energy expenditure in humans. The undercarboxylated form of OC has a reduced affinity for hydroxyapatite; hence it enters the systemic circulation more easily and exerts its metabolic functions for the proliferation of pancreatic β-cells, insulin secretion, sensitivity, and glucose tolerance. Leptin, a hormone synthesized by adipocytes, also has an effect on both bone remodelling and energy expenditure; in fact it inhibits appetite through hypothalamic influence and, in bone, stimulates osteoblastic differentiation and inhibits apoptosis. Leptin and serotonin exert opposite influences on bone mass accrual, but several features suggest that they might operate in the same pathway through a sympathetic tone. Serotonin, in fact, acts via two opposite pathways in controlling bone remodelling: central and peripheral. Serotonin product by the gastrointestinal tract (95%) augments bone formation by osteoblast, whereas brain-derived serotonin influences low bone mineral density and its decrease leads to an increase in bone resorption parameters. Finally, amylin (AMY) acts as a hormone that alters physiological responses related to feeding, and plays a role as a growth factor in bone. In vitro AMY stimulates the proliferation of osteoblasts, and osteoclast differentiation. Here we summarize the evidence that links energy expenditure and bone remodelling, with particular regard to humans.

De novo design of peptide-calcite biomineralization systems

Many organisms produce complex, hierarchically structured, inorganic materials via protein-influenced crystal growth--a process known as biomineralization. Understanding this process would shed light on hard-tissue formation and guide efforts to develop biomaterials. We created and tested a computational method to design protein-biomineralization systems. The algorithm folds a protein from a fully extended structure and simultaneously optimizes the fold, orientation, and sequence of the protein adsorbed to a crystal surface. We used the algorithm to design peptides (16 residues) to modify calcite (CaCO(3)) crystallization. We chemically synthesized six peptides that were predicted to bind different states of a calcite growth plane. All six peptides dramatically affected calcite crystal growth (as observed by scanning electron microscopy), and the effects were dependent on the targeted state of the {001} growth plane. Additionally, we synthesized and assayed scrambled variants of all six designed peptides to distinguish cases where sequence composition determines the interactions versus cases where sequence order (and presumably structure) plays a role. Scrambled variants of negatively charged peptides also had dramatic effects on calcite crystallization; in contrast, scrambled variants of positively charged peptides had a variable effect on crystallization, ranging from dramatic to mild. Special emphasis is often placed on acidic protein residues in calcified tissue mineralization; the work presented here suggests an important role for basic residues as well. In particular, this work implicates a potential role for basic residues in sequence-order specificity for peptide-mineral interactions.

Bone and glucose metabolism: a two-way street

Evidence from rodent models indicates that undercarboxylated osteocalcin (ucOC), a product of osteoblasts, is a hormone affecting insulin production by the pancreas and insulin sensitivity in peripheral tissues, at least in part through enhanced secretion of adiponectin from adipocytes. Clinical research to test whether this relationship is found in humans is just beginning to emerge. Cross-sectional studies confirm associations between total osteocalcin (OC), ucOC and glucose metabolism but cannot distinguish causality. To date, longitudinal studies have not provided a consistent picture of the effects of ucOC or OC on fasting glucose and insulin sensitivity. Further exploration into the physiological and mechanistic effects of ucOC and OC, in rodent models and clinical studies, is necessary to determine to what extent the skeleton regulates energy metabolism in humans.

Solution- and adsorbed-state structural ensembles predicted for the statherin-hydroxyapatite system

We have developed a multiscale structure prediction technique to study solution- and adsorbed-state ensembles of biomineralization proteins. The algorithm employs a Metropolis Monte Carlo-plus-minimization strategy that varies all torsional and rigid-body protein degrees of freedom. We applied the technique to fold statherin, starting from a fully extended peptide chain in solution, in the presence of hydroxyapatite (HAp) (001), (010), and (100) monoclinic crystals. Blind (unbiased) predictions capture experimentally observed macroscopic and high-resolution structural features and show minimal statherin structural change upon adsorption. The dominant structural difference between solution and adsorbed states is an experimentally observed folding event in statherin's helical binding domain. Whereas predicted statherin conformers vary slightly at three different HAp crystal faces, geometric and chemical similarities of the surfaces allow structurally promiscuous binding. Finally, we compare blind predictions with those obtained from simulation biased to satisfy all previously published solid-state NMR (ssNMR) distance and angle measurements (acquired from HAp-adsorbed statherin). Atomic clashes in these structures suggest a plausible, alternative interpretation of some ssNMR measurements as intermolecular rather than intramolecular. This work demonstrates that a combination of ssNMR and structure prediction could effectively determine high-resolution protein structures at biomineral interfaces.

The (1)H NMR structure of bovine Pb(2+)-osteocalcin and implications for lead toxicity

Structural information on the effect of Pb(2+) on proteins under physiologically relevant conditions is largely unknown. We have previously shown that low levels of lead increased the amount of osteocalcin bound to hydroxyapatite (BBA 1535:153). This suggested that lead induced a more compact structure in the protein. We have determined the 3D structure of Pb(2+)-osteocalcin (49 amino acids), a bone protein from a target tissue, using (1)H 2D NMR techniques. Lead, at a stoichiometry of only 1:1, induced a similar fold in the protein as that induced by Ca(2+) at a stoichiometry of 3:1. The structure consisted of an unstructured N-terminus and an ordered C-terminal consisting of a hydrophobic core (residues 16-49). The genetic algorithm-molecular dynamics simulation predicted the lead ion was coordinated by the Gla 24 and Gla 21 residues. It is proposed that mineral binding occurs via uncoordinated Gla oxygen ions binding to calcium in hydroxyapatite. A comparison of Pb(2+)- and Ca(2+)-osteocalcin suggests Pb(2+), at a lower stoichiometry, may induce similar conformational changes in proteins and subsequent molecular processes normally controlled by calcium alone. This may contribute to a molecular mechanism of lead toxicity for calcium binding proteins. Lead exposure may alter the amount of mineral bound osteocalcin and contribute to abnormal bone remodeling.

Identification of an osteocalcin isoform in fish with a large acidic prodomain

Osteocalcin is a small, secreted bone protein whose gene consists of four exons. In the course of analyzing the structure of fish osteocalcin genes, we recently found that the spotted green pufferfish has two possible exon 2 structures, one of 15 bp and the other of 324 bp. Subsequent analysis of the pufferfish cDNA showed that only the transcript with a large exon 2 exists. Exon 2 codes for the osteocalcin propeptide, and exon 2 of pufferfish osteocalcin is approximately 3.4-fold larger than exon 2 previously found in other vertebrate species. We have termed this new pufferfish osteocalcin isoform OC2. Additional studies showed that the OC2 isoform is restricted to a unique fish taxonomic group, the Osteichthyes; OC2 is the only osteocalcin isoform found so far in six Osteichthyes species, whereas both OC1 and OC2 isoforms coexist in zebrafish and rainbow trout. The larger size of the OC2 propeptide is due to an acidic region that is likely to be highly phosphorylated and has no counterpart in the OC1 propeptide. We propose 1) that OC1 and OC2 are encoded by distinct genes that originated from a duplication event that probably occurred in the teleost fish lineage soon after divergence from tetrapods and 2) that the novel OC2 propeptide could be, if secreted, a phosphoprotein that participates in the regulation of biomineralization through its large acidic and phosphorylated propeptide.

Evolution of matrix and bone gamma-carboxyglutamic acid proteins in vertebrates

The evolution of calcified tissues is a defining feature in vertebrate evolution. Investigating the evolution of proteins involved in tissue calcification should help elucidate how calcified tissues have evolved. The purpose of this study was to collect and compare sequences of matrix and bone gamma-carboxyglutamic acid proteins (MGP and BGP, respectively) to identify common features and determine the evolutionary relationship between MGP and BGP. Thirteen cDNAs and genes were cloned using standard methods or reconstructed through the use of comparative genomics and data mining. These sequences were compared with available annotated sequences (a total of 48 complete or nearly complete sequences, 28 BGPs and 20 MGPs) have been identified across 32 different species (representing most classes of vertebrates), and evolutionarily conserved features in both MGP and BGP were analyzed using bioinformatic tools and the Tree-Puzzle software. We propose that: 1) MGP and BGP genes originated from two genome duplications that occurred around 500 and 400 million years ago before jawless and jawed fish evolved, respectively; 2) MGP appeared first concomitantly with the emergence of cartilaginous structures, and BGP appeared thereafter along with bony structures; and 3) BGP derives from MGP. We also propose a highly specific pattern definition for the Gla domain of BGP and MGP.