Identification and characterization of a novel heparan sulfate-binding domain in Activin A longest variants and implications for function


Authors: Evan Yang aff001;  Christina Mundy aff001;  Eric F. Rappaport aff002;  Maurizio Pacifici aff001;  Paul C. Billings aff001
Authors place of work: Translational Research Program in Pediatric Orthopaedics, Division of Orthopaedic Surgery, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America aff001;  Molecular Genetics Core, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America aff002
Published in the journal: PLoS ONE 14(9)
Category: Research Article
doi: 10.1371/journal.pone.0222784

Summary

Activins regulate numerous processes including inflammation and are synthesized as precursors consisting of a long N-terminal pro-region and a mature protein. Genomic human databases currently list three activin A (Act A) variants termed X1, X2 and X3. The X3 variant is the shortest, lacks N-terminal segments present in X1 and X2, and has been the focus of most past literature. Here, we asked whether these variants are expressed by human cells and tissues and what structural features are contained within their pro-regions. Human monocytic-like cells THP1 and U937 expressed X1 and X2 variants after exposure to phorbol ester or granulocyte-macrophage colony-stimulating factor, while X2 transcripts were present in placenta. Expression vectors encoding full length X2 or X3 variants resulted in production and secretion of biologically active Act A from cultured cells. Previous studies reported a putative HS-binding domain (HBD) in the X3 pro-region. Here, we identified a novel HBD with consensus HS-binding motifs near the N-terminal end of X1 and X2 pro-regions. Peptides encompassing this new domain interacted with substrate-bound HS with nanomolar affinity, while peptides from putative X3 HBD did not. In good agreement, full length X2 pro-region interacted with heparin-agarose, while the X3 pro-region did not. In sum, our study reveals that Act A variants are expressed by inflammatory cells and placenta and yield biological activity. The high affinity HBD in X1 and X2 pro-region and its absence in X3 could greatly influence overall Act A distribution, availability and activity in physiological and pathological circumstances.

Keywords:

Biology and life sciences – Molecular biology – Molecular biology techniques – Artificial gene amplification and extension – Polymerase chain reaction – Macromolecular structure analysis – Protein structure – Protein structure prediction – Biochemistry – Proteins – Protein structure – Protein structure prediction – Extracellular matrix proteins – Protein interactions – Protein domains – Research and analysis methods – Molecular biology techniques – Artificial gene amplification and extension – Polymerase chain reaction – Database and informatics methods – Bioinformatics – Sequence analysis – Sequence motif analysis – Medicine and health sciences – Pharmacology – Drugs – Heparin

Introduction

Activins are members of the transforming growth factor-β (TGF-β) superfamily of signaling and growth factor proteins that include bone morphogenetic proteins (BMPs) and growth and differentiation factors (GDFs) [1, 2]. Activin A (Act A) was originally purified from gonadal fluids and found to stimulate release of pituitary follicle stimulating hormone [3, 4]. Since then, it has become apparent that Act A regulates many physiologic processes and events including cell proliferation and differentiation, wound repair and immune responses [58]. Changes in Act A availability or function contribute to various pathologies [9].

As other TGF-β superfamily members, Act A is initially synthesized as a pre-pro-protein that encompasses a large N-terminal pro-region and a smaller mature protein portion which is subsequently dimerized by a C-terminal disulfide bond [10, 11]. While transitioning through the secretory pathway, the precursors undergo processing by furin-like proteases that cleave them at a prescribed and obligatory cluster of basic amino acids, but do not disrupt the association of pro-region and mature protein [8, 12, 13]. The processed pro-region/mature protein dimeric complexes are finally secreted, and evidence indicates that they can exert autocrine or paracrine roles in multiple tissue and organ functions [8, 9]. Previous studies have examined activin interactions with protein antagonists including follistatin (Fs) [14] and the ability of activins to form diffusion gradients [15], presumably via interactions with extracellular matrix components [16].

The NCBI database lists three human Act A variants, termed X1, X2, and X3. These variants were initially identified and sequenced by the Mammalian Gene Collection (MGC) Program Team as expressed sequence tags (ESTs) [17], but their differences in expression or function have yet to be examined empirically in human cells and tissues. The three variants are predicted to contain the same mature protein, but differ in the length of their pro-region as a result of alternative splicing and translation start site. Much of the past and recent literature has focused on the shortest X3 variant that lacks the longer N-terminal segments characterizing the X1 and X2 variants [8, 11]. In vitro studies have indicated that bioactivity of the entire X3 Act A protein complex is similar to that of the mature protein alone, suggesting that the X3 pro-region may have roles other than influencing ligand availability or receptor-ligand interactions [18]. A previous structural analysis of the X3 variant reported the presence of a putative domain in the middle of the pro-region seemingly able to interact with heparan sulfate (HS) and heparin [19]. To date, similar analyses have not been reported for the X1 and X2 variants.

In other TFG-β superfamily members and BMPs in particular, there is ample evidence that their distribution, availability, bioactivity and turn over are tightly regulated by interactions with HS-rich proteoglycans on the cell surface and extracellular matrix [20, 21]. In pioneering studies, Cardin and Weintraub were among the first investigators to identify and characterize protein domain(s) responsible for interaction with HS [22]. Comparative analyses on vitronectin, platelet factor-4, apoliprotein A and apoliprotein B led them to identify two HS-binding amino acid sequences -XBBXBX and XBBBXXBX-, thus named CW motifs, where B represents basic residues (usually Arg or Lys) interacting with reciprocally-charged sulfated carbohydrate residues and X represents non-charged residues. Subsequent studies by multiple investigators have verified these original findings and identified similar motifs or variations in many other proteins [23, 24]. Notably, the HS chains themselves are endowed with remarkable variability in their patterns of: O-sulfation at C2, C3 and/or C6 positions and N-deacetylation and N-sulfation in glucosamine residues, and C5 epimerization of D-glucuronic acid to L-iduronic acid residues [25, 26]. These overall modifications create short segments with specific sulfation patterns and sugar structure that exhibit binding avidity toward different CW domains on various proteins [27].

One example of the fundamental importance of HS-protein interactions was provided by studies showing that mutation of basic residues within the HS-binding domain (HBD) of BMP2 and BMP4 reduced their ability to interact with HS, but greatly increased their biological activity and potency as revealed by in vitro and in vivo assays [28, 29]. These studies showed that at least for certain BMPs, interactions with HS normally serve to restrain, delimit and define their biological activities [30]. Further, we recently showed that while the HBD in BMP2 and BMP4 resides near the N-terminus of the mature protein, that domain resides at the C-terminus in BMP5, BMP6 and BMP7 [31]. Structural analyses suggested that the HBDs in these two BMP subgroups have distinct structural characteristics and spatial orientation, indicating that HBD location and structure are protein-specific [32]. In good agreement, we also found that the congenital HS deficiency characteristic of the pediatric disorder Hereditary Multiple Exostoses (HME) triggers excessive and ectopic canonical BMP signaling in growth plate perichondrium and leads to exostosis (osteochondroma) formation [33, 34], demonstrating that aberrations in normal HS-protein interactions can have severe pathological consequences.

In the present study, we have carried out a detailed examination of Act A variants and in particular, asked whether the X1 and X2 variants are expressed by human cells and tissues and what possible additional structural features may be provided by their longer N-terminal pro-regions. Our data reveal for the first time that the X1 and X2 variants are in fact expressed by human cells and placenta and contain a previously unidentified high affinity HBD toward the N-terminal end of their pro-region that is absent in the X3 pro-region. The diverse structural features of the variants could influence overall Act A availability, range and activity in physiological and pathological circumstances.

Results

Human inhibin genes and expression variants

Previous studies have provided general information and insights into the genomic organization of the human inhibinβ genes, designated INHβA, INHβB, INHβC and INHβE, that encode Act A, activin B (Act B), activin C (Act C) and activin E (Act E), respectively [3537]. To gain further structural details, we relied on the UCSC Genome Browser (genome.ucsc.edu) and NCBI sequence and gene structure databases. Computer-driven algorithms indicated that human INHβA resides on the short arm of chromosome 7 and is characterized by 7 coding exons spanning an overall genomic segment of about 25Kb (Fig 1, Tables 13 and S1 Fig).

<h2>Human <i>INHβA</i> genomic organization.</h2>
Fig. 1.

Human INHβA genomic organization.


Human INHβA is located on Chr 7 and includes seven coding exons. Arrows indicate AUG translation start site in each variant. The X1 variant is encoded by exons 2, 3, 6, 7, the X2 variant is encoded by exons 1, 6 and 7, and the X3 variant is encoded only by exons 6 and 7 but utilizes an alternative AUG start site within exon 6. This information was obtained from the UCSC Genome Browser (genome.ucsc.edu) and NCBI GenBank. Lines connecting each exon denote introns and are not drawn to scale. Also see Tables 13.

Tab. 1.

Chromosomal location and genomic organization of human activin genes.

<h2>Chromosomal location and genomic organization of human activin genes.</h2>
Tab. 2.

Exon locations within human activin genes.

<h2>Exon locations within human activin genes.</h2>
Tab. 3.

Human activin isoforms and variants.

<h2>Human activin isoforms and variants.</h2>

The NCBI database lists three variants for human INHβA designated X1, X2 and X3 and alignment of the nucleotide and amino acid sequences indicate that the variants are the product of alternative splicing and translation start sites (Fig 1). The X1 variant is encoded by exons 2, 3, 6, and 7 and the X2 variant is encoded by exons 1, 6, and 7. The X3 variant includes exons 6 and 7 only, but results from the use of an alternative translation start site within exon 6 (Fig 1, arrow in X3). The mature protein is encoded by approximately the 3’ half of exon 7 and is common to all Act A variants, which differ in the length of their pro-regions (Fig 1 and S2 Fig). It remains to be clarified whether exons 4 and 5 are incorporated into other as yet undiscovered Act A variants, if other potential exons reside upstream of coding exon 1, and what are the promoter/enhancer elements regulating INHβA gene expression [1].

It should be noted that the above transcript/protein predictions of variants are based on computational genomic analyses that were conducted initially by the MGC team using ESTs to explore the transcriptome and then subsequently confirmed by whole genome sequencing and bioinformatics [17]. Much of the past and recent extensive literature on Act A has utilized recombinant forms of the X3 variant [14, 18], but currently there is little information about the expression and structural characteristics of the X1 and X2 variants. Thus, we first asked whether these variants are expressed by myeloid- lineage cells that are amongst the cell types known to produce Act A [38, 39]. For these studies, we utilized the human monocytic-like cell lines THP1 and U937 that have the capacity to differentiate into macrophages following appropriate stimulation [40]. For initial experiments, cells were treated with phorbol 12-Myristate 13-Acetate (PMA) for 4 hrs to induce differentiation, total cellular RNAs were isolated and processed for cDNA synthesis and semi-quantitative PCR. We used unique forward primers anchored at the ATG start sites of X1 or X2 variants and combined with a common reverse primer that includes the TAG stop codon in exon 7 (Table 4). We found that PMA stimulated expression of transcripts encoding the X2 variant in U937 cells, whose identity was confirmed by Sanger sequencing of the purified PCR product (Fig 2A, lane 2, S5 Fig); however, expression of the X1 variant was undetectable (Fig 2A, lane 8). Because the X3 variant shares its entire composition with the X1 and X2 variants and does not have a unique 5’ end (Fig 1), we did not test its expression using this approach. Granulocyte macrophage-colony stimulating factor (GM-CSF) induces strong expression of Act A in purified human monocytes [41]. Thus, U937 and THP1 cells were untreated or treated with GM-CSF (1 ng/ml; 8 hr) and the resulting RNA samples were processed for PCR analysis as above. Interestingly, GM-CSF upregulated the X2 variant in both cell lines (Fig 2A, Lanes 3, 6), but only stimulated expression of the X1 variant in U937 cells (Fig 2A, lane 9). To make sure both cell lines have the cell surface machinery to respond to GM-CSF, we subjected them to FACS analysis using antibodies to its receptor CD116 or pre-immune IgGs as control. Both cell lines exhibited readily detectable levels of CD116 (Fig 2B and 2C).

Tab. 4.

PCR Primers used to amplify human activin transcripts.

<h2>PCR Primers used to amplify human activin transcripts.</h2>
<h2>Human monocytic-like cells express Act A X1 and/or X2 variants and CD116.</h2>
Fig. 2.

Human monocytic-like cells express Act A X1 and/or X2 variants and CD116.


A, Cells were treated with PMA or GM-CSF. Following treatment, RNAs were processed for PCR analysis of Act A variant expression. Amplified bands were purified from the gels and subjected to Sanger sequencing to confirm transcript identity. Treatment of U937 cells with PMA induced expression of Act A-X2 (lane 2), while treatment with GM-CSF induced expression of both Act A-X2 and -X1 (lanes 3, 9). Treatment of THP-1 cells with GM-CSF induced expression of Act A-X2 (lane 6). Lane 10, Fisher exACTGene DNA ladder. B and C, Scatter plot and histogram depicting CD116 detection. U937 and THP1 cells were incubated with anti-CD116-PE antibody or IgG-PE control. The cells were washed, and antibody binding was determined by FACS. Note that both cell lines contain cell surface CD116. Statistical Analysis (Student’s t test) for B: IgG vs CD116: ****: P< = 0.0001, *: P< = 0.05.

To determine which Act A variant(s) are expressed under normal physiological conditions, we prepared cDNA from human placental RNA and processed it for PCR using different primer sets. A first set amplified an amplicon covering exons 6 and 7, present in all three Act A variants (Fig 3A and 3C, black arrows). Next, we used primer sets containing unique forward primers specific for the 5’ ends of X1 or X2 variants (Fig 3B, white arrows). This yielded a product encoding the 5’ region of the X2 variant, but no amplification product homologous to the X1 5’ region was obtained (not shown). Thus, the Act A X2 variant is expressed in human placenta.

<h2>Normal human placenta tissues contain transcripts encoding Act A X2 variant.</h2>
Fig. 3.

Normal human placenta tissues contain transcripts encoding Act A X2 variant.


cDNA was prepared from human placenta RNA and used as template to amplify Act A transcripts by PCR. The PCR products were size fractionated on 1% agarose gels, and bands were purified and subjected to Sanger sequencing to confirm transcript identity. A, amplification of exons 6 and 7 using primers indicated by black arrows, yielded a predicted product of 1,443 bp. B, amplification of 5’ end of Act A-X2, using primers indicated by clear arrows, yielded a predicted amplicon of 548 bp. C, schematic representation of X2 variant and primer sets (White and Black arrows) used to amplify products shown in A and B. Size markers; Fisher exACTGene DNA ladder is on the left, lane 10. Also see Fig 1.

Together, the data provide new evidence that Act A X1 and X2 variants are expressed by differentiating human monocytic cell lines and human tissue. Because the X1 and/or X2 variants are expressed in response to specific chemical or biological cues, it is possible that they may be expressed and utilized in vivo during specific stages of development and under certain pathologic conditions.

AD293 and COS1 cells secrete biologically active Act A following transfection with X2 and X3 expression vectors

Next, we explored the biological activity of Act A produced and secreted by cultured cells. The entire open reading frame encoding each variant was cloned into the mammalian expression vector pcDNA3.1+, which utilizes the CMV promotor to drive protein expression. Following verification of constructs by DNA sequencing, each expression construct along with a construct expressing green fluorescent protein (GFP) was transfected into AD293 cells using FuGENE 6 Transfection Reagent (Promega), resulting in a transfection efficiency of > 90% (Fig 4B). Three days after transfection, conditioned medium was harvested and the remaining cells were scraped from the dishes, RNA was isolated, treated with DNAse and subsequently used for cDNA synthesis. Cells transfected with each Act A variant expressed the appropriate mRNA (Fig 4C, lanes 3–5). Removal of plasmid DNA was confirmed using PCR primers specific for beta-lactamase, which encodes the AmpR gene on pcDNA3.1 expression plasmid (Fig 4C, right). This amplicon was readily amplified from plasmid DNA (Fig 4C, lane 6), but was absent in cDNAs prepared from cells (Fig 4C, lanes 7–10).

<h2>Expression and secretion of Act A variants.</h2>
Fig. 4.

Expression and secretion of Act A variants.


Expression constructs for full length X1, X2 and X3 variants were cloned into pcDNA 3.1 and transfected into AD293 and COS1 cells. A, diagram of Act A expression cassette with His tag on the 3’ end. B, GFP detection indicates that transfection efficiency was > 90%. C, RT-PCR analysis (left panel) demonstrating that X1, X2 and X3 transcripts were efficiently expressed following transfection (lanes 3–5), producing an amplicon identical to that of the expression plasmid (lane 1). Untransfected cells are shown in lane 2. RT-PCR analysis (right panel) indicating that plasmid DNA was efficiently removed from our RNA preps (lanes 8–10) but was amplified from plasmid DNA itself (lane 6). D, immunoblots of conditioned medium samples separated under reducing conditions and probed with anti-His tag antibody. Note that samples from AD293 and COS1 cells transfected with X2 or X3 expression constructs exhibit an immunoreactive band ~13 kDa that has the predicted size of mature Act A monomer (lanes 2, 3, 5, 6, 10 and 11). The band was undetectable in conditioned medium from cells transfected with the X1 vector (lanes 1, 4 and 9) or un-transfected cells (lanes 7 and 8). Lanes 4–6 are a longer exposure of lanes 1–3. E, THP1 cells were incubated for 2 hrs with conditioned media from AD293 cells that were untransfected (U) or transfected with X2 or X3 constructs. The cells were processed for q-PCR analysis of TIEG expression, an early response gene induced by Act A and TGF-β. Note that TIEG expression is up-regulated by treatment with recombinant human Act A (rhActA; lane 2, positive control) as well as conditioned medium from X2 or X3 expressing cells (lanes 4 and 6). This induction was blocked by pre-incubation with a neutralizing antibody to Act A (lanes 3, 5 and 7). Conditioned medium from untransfected cells elicited no TIEG expression (U, lane 1). Statistical Analysis (Student’s t test): 1 vs 2, 1 vs 4 and 1 vs 6; *: P< = 0.05, ***: P< = 0.001.

To assess protein expression, the conditioned media sample from each transfected cell population was incubated with Talon resin, and the bound proteins were eluted with 0.5M Imidazole and analyzed on immunoblots. AD293 Cells transfected with X2 or X3 expression vectors both secreted a protein of about 13 kDa, consistent with the mass of mature Act A protein under reducing conditions (Fig 4D, lanes 2, 3, 5 and 6; note that lanes 4–6 are over-exposed images of lanes 1–3). This protein band was more prominent in X2 versus X3 transfected cells, but was not detectable in cells transfected with X1 construct (Fig 4D, lanes 1 and 4). To verify these results, Cos1 cells were transfected with each Act A variant expression construct. Again, immunoblot analysis of conditioned media samples showed secretion of mature Act A ligand only by cells transfected with X2 or X3 variant constructs (Fig 4D, lanes 10 and 11). To test whether the secreted proteins were biologically active, THP1 cells were treated with commercially available recombinant Act A (rhAct A) or conditioned medium samples obtained from cells transfected with the Act A expression constructs. After 2 hr treatment, RNA was extracted, cDNA was prepared and expression of TIEG -an early response gene upregulated by Act A or TGF-β was assessed [42]. As expected, TIEG was rapidly upregulated when cells were treated with rhAct A used as positive control (Fig 4E, lane 2). A similar strong up-regulation was observed when the cells were exposed to the conditioned medium from cells transfected with X2 and X3 expression constructs (Fig 4E, lanes 4 and 6). Medium from un-transfected cells produced no response (Fig 4E, lane 1). TIEG induction was blocked when the cells were co-treated with a monoclonal Act A neutralizing antibody (Fig 4E, lanes 3, 5 and 7). Thus, AD293 and Cos1 cells transfected with X2 and X3 variant expression vectors secrete biologically active Act A.

Act A X1 and X2 variants contain a previously unsuspected HS-binding domain (HBD)

Many TFG-β superfamily members are known to establish high affinity interactions with HS-rich proteoglycans on the cell surface and extracellular matrix that regulate protein distribution, availability, bioactivity and turn over [20, 21]. A study published in 2010 analyzed the HS-binding properties of Act A X3 variant and reported the identification of a putative HBD located toward the middle of the pro-region and encompassing amino acids 259 to 272 (Fig 5A) [19]. This segment corresponds to amino acids 381–394 and 316–329 in the X1 and X2 variants, respectively (Fig 5A) and contains several clusters of lysine residues (Fig 5A–5D, highlighted in green). Notably, this domain also contains four glutamic acid residues, which are predicted to interfere with interactions with the anionic sulfate groups in HS. Thus, we closely reexamined the entire amino acid sequence of the X3 variant and compared it to the X1 and X2 variants. We noted the presence of another basic residue-rich domain located near the N-terminal end of the X1 and X2 pro-regions that closely resembled a CW motif (Fig 5A, highlighted in light blue). The domain spans amino acids 78–102 in X1 and amino acids 13–37 in X2, but is absent in X3 (Fig 5A–5D).

<h2>Schematic representation of human Act A variants.</h2>
Fig. 5.

Schematic representation of human Act A variants.


A, the three variants, X1, X2, and X3 listed in the NCBI protein database are depicted as a linear segment proportional to their protein sizes that are 548, 483 and 426 amino acids long, respectively. The box colored in aqua depicts the location of the novel high affinity HBD identified in the X1 and X2 variants. The green box depicts the location of the previously-reported region in the X3 variant (presumed to bind to HS) that is also present in the X1 and X2 variants. B-D, the entire amino acid sequences of the X1, X2 and X3 variants are shown in which the highlights in aqua and green depict the protein segments as in A, while the highlight in yellow depicts the mature protein portion shared by all variants.

To determine whether the newly identified domain is in fact endowed with HS binding capacity, we prepared a synthetic peptide encompassing this region, which was biotinylated on the N-terminus. The domain is identical in X1 and X2 variants and denoted as X213-37 (Fig 6A). In addition, we prepared a biotinylated peptide encompassing the previously described putative HS-binding region in X3 pro-region [19], which is designated X3259-272 (Fig 6B). The peptides were tested in solid-phase binding assays, employing 96 well plates pre-coated with HS, following standard protocols established in our previous studies [31]. We found that the X213-37 peptide interacted with HS with high affinity, exhibiting a Kd of 4 ± 2 nM (n = 3); binding was blocked by pre-incubation of the peptide with excess heparin serving as competitor (Fig 6C). As we anticipated given its amino acid composition, the X3259-272 failed to interact with HS at any concentration tested (Fig 6C). Changing incubation conditions or lengthening assay time did not change outcomes.

<h2>Peptides corresponding to the novel HBD bind HS with high affinity.</h2>
Fig. 6.

Peptides corresponding to the novel HBD bind HS with high affinity.


96-well plates were coated with HS and incubated with tetrameric peptide complexes (HBC) for 2 hrs. Plates were washed, and bound peptides were assessed by addition of OPD substrate. A and B, schematic diagram depict the amino acid sequences of Act A- X213-37 and Act A- X3259-272 peptides and their location within the full-length proteins, respectively. C, solid phase binding assays of Act A peptides to HS. Note that Act A- X213-37 binds HS with high affinity and its binding is blocked by addition of soluble heparin. In contrast, Act A-X3259-272 displays minimal binding. Inset, double reciprocal plot that yields a Kd of 4 nM (n = 3) for Act A- X213-37. NA-HRP backbone displays minimal binding to HS.

Solid-phase assays are widely used to test macromolecular interactions, but may not fully reflect biological conditions. Thus, we tested the peptides for interaction with the cell surface. U937 cells were lightly fixed with formalin to preserve the cell surface and then incubated with biotinylated peptides complexed with NA-488 [31]. After careful rinsing, cells were processed for FACS analysis. The X213-37 peptide readily and efficiently interacted with the cell surface, and binding was blocked by pre-treatment with heparin (Fig 7A and 7B). However, X3259-272 peptide failed to interact with the cells (Fig 7A and 7C), consistent with results from the solid phase binding assay (Fig 6C). To verify assay reliability, we concurrently tested a peptide that encompasses the HBD of BMP4, characterized in a previous study [31]. As expected, the BMP4 peptide effectively interacted with the cell surface and was competed out by excess heparin (Fig 7A). No fluorescence signal was generated by untreated control cells or by cells incubated with NA-488 in the absence of peptides.

<h2>HBD peptides bind to the cell surface with avidity.</h2>
Fig. 7.

HBD peptides bind to the cell surface with avidity.


A, U937 cells were incubated with peptide complexes consisting of biotinylated peptides pre-incubated with NA-488. Following incubation, the cells were washed, and bound peptide was assessed by FACS analysis. Note that Act A- X213-37 readily binds to the cell surface, and its interaction is blocked by soluble heparin competitor. In contrast, Act A-X3259-272 peptide fails to bind. Binding specificity was verified using a previously characterized peptide encompassing the high affinity HBD from BMP4. B, representative FACS plots show binding of Act A- X213-37 in the absence (blue line) or presence (red line) of heparin competitor. C, representative FACS plots showing binding of Act A- X213-37 (blue line), lack of binding by Act A-X3259-272 (red line) and background signal generated by U937 cells alone (green line). Statistical Analysis (Student’s t test): NA488 vs ActA X213-37, NA488 vs BMP4: ****: P< = 0.0001, **: P< = 0.01.

To further verify the results, we prepared His-tagged expression constructs in pcDNA3.1 vector encompassing the entire pro-region of X2 and X3 variants and transfected each into AD293 cells. Each construct was efficiently transcribed into mRNA (Fig 8B). Conditioned media samples were incubated with Talon resin or heparin agarose, bound protein was eluted and analyzed on immunoblots with His-tag antibodies. Both pro-regions were produced and secreted by the cells (Fig 8C, lanes 2 and 3). In good agreement with the data above, the secreted X2 pro-region interacted with heparin agarose, while the X3 pro-region did not (Fig 8C, lanes 5 and 6).

<h2>Expression and binding characteristics of Act A X2 and X3 pro-regions.</h2>
Fig. 8.

Expression and binding characteristics of Act A X2 and X3 pro-regions.


A, diagram of Act A pro-region expression cassette with HIS tag on 3’ end. Cassette was cloned into pcDNA 3.1. Expression constructs encoding the entire pro-region of Act A X1, X2 or X3 variants were transfected into AD293. Cells were processed for RT-PCR analysis four days later. B, RT-PCR product analysis shows that each variant was expressed efficiently (lanes 3–5) and the amplification product matched that produced by expression vector itself (lane 1). RNA from un-transfected cells (U, lane 2) produced a faint band likely representing endogenous Act A transcripts. C, conditioned media samples from each cell population above were incubated with TALON resin or heparin agarose (30 min, 4C). After rinsing, bound proteins were eluted and analyzed by immunoblot with anti-His tag antibody. Note the presence of a thick protein band of ~ 40 kDa (exhibiting the predicted size of the pro-region) in samples from both X2 or X3 pro-region-expressing cells (lanes 2 and 3) but not in samples from un-transfected cells (U, lane 1). Note also that only the X2 encoded pro-region binds heparin agarose (lane 5), while X3 encoded pro-region does not (lane 6).

3D peptide and protein conformation analyses

The above data predict that the X213-37 HBD should possess a 3D configuration and amino acid distribution able to create a surface(s) amenable to interactions with HS. However, these features and properties should be absent or disrupted in the segment encompassing the X3259-272 segment. To test these predictions, we utilized the I-TASSER server [43] as in previous structural studies on BMP family members [20, 31]. We found that the X213-37 domain displayed a helical structure Fig 9A). When examined using helical wheel diagrams, it displayed basic residues (i.e., Lys) mostly on one side, thus generating a cationic surface on one face of the helix, conducive to HS binding (Fig 9A, 9C and 9E). The X3259-272 segment also displayed a helical configuration, but structural views and helical wheel projections showed that its residues did not elicit a cationic surface suitable for HS interaction (Fig 9B, 9D and 9F). Again, the presence of negative charges from the 4 glutamic acid residues within this sequence are predicted to interfere with HS binding (Fig 9F, highlighted in red).

<h2>Peptides display a helical configuration.</h2>
Fig. 9.

Peptides display a helical configuration.


Three dimensional analyses and structure predictions were carried out using the I-TASSER server. A-B, side views of peptides oriented with the N-terminus labeled (N) on the left indicate that both the Act A-X213-37 and Act A-X3259-272 peptides are predicted to possess a helical configuration. C-D, side-view of models depicted in A-B are rotated 90° and provide a frontal view down the axis of the alpha helix, beginning with the N-terminus. Acidic and basic residues are highlighted as follows: Glu; red; Lys, purple; and Arg, cyan. Overall amino acid sequences are noted below, with putative CW motifs underlined. Note in C that the distribution of Lys and Arg residues in Act A-X213-37 elicits uniformly basic surfaces amenable to HS binding. In contrast, the 4 Glu residues in Act A-X3259-272 would disrupt such HS interactions (D). E-F, helical wheel diagrams illustrate the stereological position of amino acids in sequences underlined in C and D, respectively. It is evident that the glutamic acid residues (in red) prevent formation of a uniform basic surface around the helix.

Next, we used the I-TASSER server to visualize the relative location, orientation and configuration of the X213-37 and X3259-272 domains within the full-length X2 variant. The generated structures revealed that the domains occupy locations on opposite sides of the protein (Fig 10 and S1 and S2 Videos). As also seen in previous crystallography studies of the full-length X3 variant [18], the pro-region is physically separated from the mature protein and connected to it via the latency lasso (Fig 10B). To more clearly determine whether the X213-37 domain created an exposed HS-binding surface, we analyzed the structures at different rotational angles. Orthogonal views revealed that the domain created a cationic surface facing away from the protein moiety and likely accessible for interactions with HS (see S6 and S7 V\. In contrast, the X3259-272 segment was also positioned on the outer perimeter of the protein, but its acidic amino acid composition prevented it from generating a suitable cationic surface (Figs 9 and 10).

<h2>Predicted structure of Act A variants and location of putative HBDs.</h2>
Fig. 10.

Predicted structure of Act A variants and location of putative HBDs.


Analyses were carried out with the aid of the I-TASSER server for protein structure and function prediction and visualized using Chimera. Models are based on the crystal structure of pro-TGF-β1 (PDB: 3RJR). A, ribbon structure of full length Act A X2 variant monomer in which the pro-region is in orange and the mature protein is in green. Note the location of the high affinity X213-37 HBD (bracketed) and that of the previously reported X3258-272 segment (bracketed). Note that the former appears to be exposed to the surface of the monomer that would favor interactions with HS. B, same ribbon structure, but viewed after a 90° rotation to provide further stereological insights. Also see S1 and S2 Videos in Supporting Information.

Discussion

Our data reveal that the X1 and X2 variants of Act A are expressed by human inflammatory cells in response to PMA and GM-CSF and in human placenta. The variants differ from each other at the N-terminal end of their pro-region, but share the bulk of the pro-region and the entire mature protein. Our data also show for the first time that the X1 and X2 variants possess a domain in the pro-region that exhibits typical CW characteristics and binds to HS and the cell surface readily and with high affinity. As suggested by 3D stereological analysis, the basic residues in the domain are organized to produce an anionic surface that is adept at HS binding and oriented toward the outer perimeter of the pro-region. Because Act A is assembled and secreted as a dimer held together by an intermolecular disulfide bond [18], each dimer would possess two such HBDs oriented away from each other and potentially functioning autonomously. This is reminiscent of the situation in certain BMP family members such as BMP2 and BMP4 in which the two HBDs in each dimer are distinct in terms of location, action and interactions and act in a bivalent manner in binding assays [29]. However, this appears not be the case for BMP5, BMP6 and BMP7 where the two HBDs come together to form a single central cationic surface within each dimer, likely as a result of their location at the C-terminus of each of the two closely juxtaposed mature proteins [31]. These observations and considerations highlight the fact that the HBDs in different proteins differ in location, organization, specific amino acid sequences and number of basic clusters, though they all share fundamental traits of CW motifs. The significant structural variability amongst HBDs in different proteins are likely to have physiological significance in terms of protein availability, distribution and orientation within the extracellular matrix or cell surface as well as with respect to protein interactions with specific HS-proteoglycans and even endogenous protein competitors and antagonists [20]. This conclusion is also sustained by the fact that the domains in many such proteins are evolutionarily conserved [20], but details about their functional significance still need to be fully clarified and understood.

In the standard binding assays used here, peptide X3259-272 spanning the previously reported HS-binding region of the X3 variant [19] failed to interact with substrate-bound HS or the cell surface (Figs 6 and 7). Moreover, the entire pro-region of X3 failed to interact with heparin, while the pro-region of the X2 variant did interact (Fig 8). These results are not surprising given that the clusters of basic residue clusters within X3259-272 contain 4 intervening glutamic acid residues, which would be predicted to prevent the formation of a uniform cationic HS-binding surface, as our 3D data indicate. How then can our data be reconciled with those of Li et al. [19]? A major difference is that in our solid phase and FACs assays (Figs 6 and 7) we made use of peptides corresponding to specific protein domains, while Li et al. made use of full length recombinant X3 variant and perlecan in their solid phase and pulldown experiments. It is possible that full length proteins are able to establish a number of diverse macromolecular interactions with each other in addition to those involving specific segments or motifs interacting with HS chains [44]. HS-rich glypican family members have been shown to interact with Wnts, hedgehogs and other signaling proteins in complex manners, also involving their core protein [45, 46]. In addition, differences in experimental conditions including salt concentration, stringency, washing conditions, buffers, detergents and other parameters could influence outcomes in different studies. We have tested HS derived from different sources including kidney, basement membranes of Engelbreth-Holm-Swarm mouse sarcoma and obtained from different vendors, but the results were the same. Nonetheless, the inhibition of X3 variant/perlecan interactions by heparitinase pre-treatment or mutation of basic residue clusters observed by Li et al. [19] does suggest that those clusters could have some binding roles, despite the presence of acidic residues. However, their affinity and binding strengths for HS could be far lower than those of the new HBD, X213-37, described here. It is worth recalling that the three major follistatin variants differ in their ability to bind HS, with the shorter variant Fs288 exhibiting the highest affinity and the longest variant Fs315 exhibiting the lowest, likely due to the presence of acidic residues [47]. Hence, future studies will be needed to compare side-by-side the binding affinities and kinetics of the X1, X2 and X3 variants and their respective mutants to HS-rich proteoglycans, define the binding hierarchy of different domains, and decipher their roles and underlying mechanisms.

The biosynthesis and elongation of HS chains at serine residues on the proteoglycan core proteins are complex processes involving over 25 enzymes. As indicated above, these step-wise biosynthetic steps results in the creation of segments with specific sulfation and sugar structure patterns, also referred to as ‘sulfation codes’ [26]. Following secretion, the chains can be further modified by tissue-specific sulfatases that selectively remove sulfate groups at specific positions. The resulting diversity and complexity of sulfated segments along the HS chains are thought to serve as specific sites of protein interactions [20, 21]. Complementarily, the diverse HBDs identified in different proteins likely add specificity and enable each protein to recognize and bind distinct segments along HS chains with specific sulfation codes. Some of the most characterized protein-HS interactions are those by members of the fibroblast growth factor (FGF) family of signaling molecules [48, 49]. FGFs contain a nine amino acid domain or ‘glycine box’ that mediates the ability of full-length proteins to selectively bind octasaccharides with different sulfation patterns [48]. Our group recently demonstrated that the HBD present in BMP2 and BMP4 competitively inhibits binding of full-length BMP2 protein to HS, while the HBD in BMP5 does not [31]. While we do not know the preferred sulfation code recognized by ActAX213-37 HBD, it may enable the pro-region/mature protein complex for X1 and X2 variants to specifically interact with HS-rich proteoglycans on the cell surface and extracellular matrix, thereby regulating signaling, turn over and function.

As pointed out above, past and recent studies on Act A biology and function have made use of the X3 variant or the mature protein but to our knowledge, the X1 or X2 variants have not been tested. For example, Harrington et al. used recombinant Act A mature protein to clarify the structural basis by which Fs inhibits Act A signaling [14]. They found that Fs dimers wrap themselves around the back of activin dimer wings reorienting their own EGF- and Kazal-like subdomains. Wang et al. carried out a detailed analysis of mechanisms involved in initiation of Act A biological activity, using the X3 variant [18]. They found that intracellular furin cleavage is required to turn on the biological activity of the newly-synthesized X3 variant; as already pointed out above, the signaling potency of the furin-cleaved X3 complex, containing the pro-region and active ligand, was nerly identical to that of mature protein. Notably, this group found that removal of residues 259 to 272, encompassing the X3259-272 domain described by Li et al. [19], had no effect on biological activity, pointing to a minor if any role of that domain as our data also suggest. Quaternary crystal structure analysis indicated that the two X3 chains form a cross-armed, domain-swapped dimer in which the two central disulfide-linked mature domains wrap around each other and create a connecting bridge to the reminder of the outward oriented globular pro-regions, a configuration different from that in other superfamily members such as TGF-β1. Lastly, Johnson et al. compared the biological activity and half-life of recombinant (cleaved) X3 variant versus mature protein after intravenous injection in male rats [50]. They found that both proteins elicited similar biological activity (monitored by follicle stimulating hormone levels), but the circulating half-life of the entire X3 variant was approximately twice as long as the mature protein (t1/2 fast = 12.5 min vs 5.5 min, respectively). Undoubtedly, the above studies and related efforts [5] have provided numerous important and novel insights into Act A biology. However, in view of the data described here, one has to wonder whether testing the X3 variant and mature protein only can actually uncover and clarify the overall physiologic spectrum of Act A activities and underlying mechanisms. Given that the X1 and X2 variants have a more extended N-terminal pro-region and possess a high affinity HBD that is absent in the X3 variant, the X1 and X2 variants could behave and act differently than the X3 variant or mature protein. For instance, the half-life of systemically-injected X1 and X2 variants is likely to be significantly different compared with that of X3 variant or mature protein, given their expected ability to interact with multiple constituents in both blood and tissues through their high affinity HBD. Interactions of X1 and X2 variants with follistatin could also differ from those of X3 due to their larger structure and presence of the high affinity HBD (X213-37) which could affect the binding of Fs288 or Fs315 variants. Within developing embryonic tissues, Act A is thought to form signaling gradients, but this important process has mainly been studied using exogenous mature protein [51]. It will be interesting to reexamine these issues using the X1 and X2 variants and determine endogenous Act A variant expression, diffusion and gradient formation. In sum, it is possible that previous observations and conclusions derived from studies utilizing mature Act A or the X3 variant may need to be verified and extended by also focusing on the X1 and X2 variants.

In the musculoskeletal research field, interest in Act A biology has been re-invigorated by recent evidence of its pathogenic roles in the pediatric congenital disorder Fibrodysplasia Ossificans Progressiva (FOP) [52]. FOP is characterized by episodic formation of extraskeletal cartilage and bone tissue–a process collectively termed heterotopic ossification (HO)—at multiple sites, eventually involving much of the body [53]. HO is preceded by local inflammation and over time, the accumulated ectopic tissue causes multiple health problems and impediments—including skeletal deformities, chronic pain, growth impairment and joint ankylosis—that increasingly interfere with basic daily functions and can lead to breathing difficulties and premature death. FOP patients bear a missense activating mutation in ACVR1 that encodes the BMP receptor ALK2 [54]. It was previously thought that activity of mutant ALK2 would be boosted by BMP ligands present at HO initiation sites; given that ALK2-driven BMP signaling mediated by canonical phosphorylated SMAD1/5/8 is chondrogenic, its enhanced activity would lead to local chondrogenesis followed by osteogenesis and ultimately, HO development [55]. However, recent studies have indicated that unlike wild type ALK2, mutant ALK2 can also respond to Act A and elicits pSMAD1/5/8 signaling, thus directly linking inflammation to ectopic chondrogenesis and HO formation [52]. These interesting and important studies have left many questions unanswered. It is assumed that inflammatory cells invariably recruited to the prospective site of HO formation are the source of Act A, but to date, this has not been demonstrated. Further, it is unclear which Act A variant(s) are expressed and how the mature protein ligand becomes available for interaction with mutant ALK2 in skeletogenic cells. Likewise, Act A normally signals via pSMAD2/3 [2], and it remains unknown how its interactions with mutant ALK2 would elicit a switch to pSMAD1/5/8 signaling and what combination(s) of type I and type II receptors are involved. Given our demonstration of the high affinity HBD in the pro-region of X1 and X2 Act A variants, it will be interesting to consider the possibility that synthetic mimetic HS oligomers [56] could be used to sequester Act A precursors, rendering the mature protein less available for action at sites of HO formation.

Experimental procedures

Reagents

Human Placental (AM7950), Liver (AM7960) and Lung (AM7968) RNA, NeutrAvidin (NA), NA-HRP and NA-DyLight 488 (NA-488) and Fisher exACTGene DNA Ladder (100–10,000 bp) were obtained from Thermo-Fisher. Heparin was obtained from Sagent Pharmaceuticals and HS was obtained from Millipore-Sigma (www.sigmaaldrich.com). Phorbol 12-myristate 13-acetate (PMA; aka 12-O-Tetradecanoylphorbol-13-Acetate) was obtained from Cell Signaling (#4174). All DNA Oligonucleotide primers were obtained from Integrated DNA Technologies (www.idtdna.com) and are listed in Table 4.

RNA, PCR & DNA sequencing

Total RNA was isolated using TRIzol Reagent (Thermo-Fisher) following the manufacturers protocol. Five μg of glycogen (Thermo-Fisher AM9510) was added to the aqueous phase prior to the addition of isopropanol to facilitate RNA precipitation. For some experiments, RNA quantification was determined using a Nanodrop spectrophotometer. cDNA was prepared from 2 μg of purified RNA, using oligo dT as primer, supplied with the Verso cDNA synthesis kit (Thermo-Fisher # AB1453A) or SuperScript IV Reverse Transcriptase (Thermo-Fisher #18090010) following the manufactures protocol. All PCR reactions were performed using Act A specific primers and GoTaq Green Master Mix (Promega) following the manufacturers protocols. PCR reactions were optimized by gradient PCR, using an Applied Biosystems Veriti Thermal Cycler. PCR products were size fractionated on 1% agarose gels, the bands were isolated, DNA was extracted from gel slices using Micro Bio-Spin columns (Bio Rad# 7326204) and sequenced. All DNA sequencing was performed by the NAPCore Facility at CHOP.

Design of Act A protein expression constructs

The entire open reading frame for Act A variants X1, X2 and X3 were cloned into the EcoRI/XbaI site of the mammalian expression vector pcDNA3.1 (+). A poly-His tag was added to the 3’ end (i.e., C-terminus) to facilitate purification via immobilized metal affinity chromatography (IMAC). The pro-region from each variant was cloned in an analogous manner, except that the His tag was placed before the Furin cleavage site (QSEDHPHHHHHHHH) (see Table 5). Each expression construct was analyzed by Sanger DNA sequencing for verification. For transfection, 2.5 μg of plasmid DNA plus 7.5 μl FuGENE 6 Transfection Reagent (Promega) was added to 0.25 ml Opti-MEM medium and incubated for 15 minutes at room temperature. The transfection mix was subsequently added to cells growing in medium containing 10% fetal calf serum and 24 hr later, the medium was changed to Serum-free media containing Insulin-Transferrin-Sodium Selenite supplement (Roche 11074547001) and G418 (100 μg/ml). Cells were grown for 3 days, the medium was collected, clarified by centrifugation and incubated with TALON Metal Affinity Resin (Takara Bio USA). The resin was washed with PBS, transferred to a Micro Bio-Spin 6 Column (Bio-Rad); bound protein was eluted with 0.5M Imidazole. Purified proteins were analyzed on SDS-polyacrylamide gels/Western blots under reducing conditions and probed with anti-His Tag Monoclonal Antibody (HIS.H8). Western blots were developed using HRP- or IRDye- tagged Secondary Antibodies (Li-Cor).

Tab. 5.

Human activin A full-length and pro-regions.

<h2>Human activin A full-length and pro-regions.</h2>

mRNA expression of Act A constructs

To evaluate expression of Act A expression constructs, transfected cells were washed and pelleted. RNA was isolated by combining Trizol and Direct-zol RNA MiniPrep Kit (Zymo Research #2050) following the manufacturers protocols. During isolation, RNA is treated with DNAase, to remove any contaminating plasmid or genomic DNA. RNA concentration and purity was assessed using a NanoDrop 2000 spectrophotometer. Act A mRNA expression was evaluated using forward and reverse primers by RT-PCR.

Biological activity of secreted Act A

To assess biological activity of the expressed Act A variants, THP1 cells were incubated with 3 ml of conditioned medium obtained from cells transfected with each Act A variant for 2 hr, in the absence or presences of an ActA neutralizing antibody (BioLegend # 94746) the cells were pelleted and total RNA was isolated using TRIzol Reagent as described above. Induction of Kruppel-like factor 10 (KLF10, TIEG; accession NM_005655) mRNA, an early response gene induced by Act A and TGF-β was assessed by qPCR.

Peptides

Peptides encompassing the heparin binding domains, contained the putative Cardin-Weintraub (CW) motif and a minimum of 2 flanking amino acids on the N- and C- terminal sides of the domain. Three Glycine residues were added to the N-terminus of each peptide to provide a flexible linker between an N-terminal biotin tag and the peptide itself (See Table 6). Peptides were purified by reverse phase HPLC and peptide mass was confirmed by MALDI-TOF mass spectrometry. All peptides were synthesized and purified by Peptide 2 (Chantilly, VA). Where indicated, biotinylated peptides were oligomerized into tetramers by incubation with NAHRP (each avidin molecule binding 4 biotins ([57]) as described [31].

Tab. 6.

Activin peptides.

<h2>Activin peptides.</h2>

Cells

U937 (human pro-monocytic) and THP-1 (human monocytic) cell lines were grown in RPMI containing 10% fetal bovine serum in a humidified atmosphere of 5% CO2. Cells were treated with PMA (100 ng/ml; 160 nM) for 4 hr (23) or GM-CSF (1 ng/ml) for 8 hr in RPMI containing 1% fetal bovine serum.

FACS analysis

Expression of CD116 (α subunit of GM-CSF receptor). U937 and THP1 cells were washed with PBS and fixed with 2% buffered formalin on ice for 20 min. The cells were washed with PBS and blocked by incubation in PBS, 1% BSA for 20 min. Cells were incubated with anti-Human CD116-PE (BD 551373) or IgG-ĸ-PE isotype control (BD 554680) on ice for 2 hr, washed with PBS and analyzed on a BD Accuri Flow Cytometer (FL2 channel), located in the Flow Cytometry Core Laboratory at CHOP.

Peptide binding to intact cells

U937 cells were washed with PBS and fixed with 2% buffered formalin for 20 min on ice. The cells were washed with PBS and blocked by incubation in PBS, 1% BSA for 20 min on ice. Approximately 106 cells (100 μl) were incubated with peptide tetramers consisting of biotinylated peptide complexed with NA-488. Following incubation for 2 hr on ice, the cells were washed and analyzed on a BD Accuri Flow Cytometer (FL1 channel), located in the Flow Cytometry Core Laboratory at CHOP.

In Silico modeling

Amino acid sequence alignments were performed using the T-coffee alignment tool (tcoffee.crg.cat). Helical wheel projections were constructed using using NetWheels: Peptides Helical Wheel and Net projections maker (lbqp.unb.br/NetWheels). Secondary structure predictions were carried out using the I-TASSER server for protein structure and function prediction [43] and the resulting structures were visualized using Chimera (www.cgl.ucsf.edu/chimera). The structures of Act A-X1 and Act A-X2 are based on the crystal structure of pro-TGF-β1 (PDB: 3RJR).

Solid phase binding assays

Nunc MaxiSorp flat bottom 96 well plates were coated with HS in 50 mM carbonate buffer (pH 9.4) overnight at 4°C as described [31]. All binding assays were carried out in 1 x PBS, 0.1% tween 20 (PBST) containing 1% bovine serum albumin. Plates were incubated with peptides for a minimum of 2 hr at room temperature with gentle shaking. At the termination of the assay, the plates were washed 3 times with PBST and a final rinse with PBS. Plates were developed by addition of HRP substrate O-phenylenediamine dihydrochloride (OPD) in Citrate-Phosphate buffer (25 mM citric acid, 50 mM sodium phosphate, pH 5) and read at 450 nm in a Bio-Tek Synergy HT plate reader. Very low levels of background binding were exhibited by NA-HRP alone or when complexed with biotinylated-BSA.

Statistical analysis

Statistical evaluation of data was performed using Student’s t-test employing GraphPad Prism software (www.graphpad.com).

Supporting information

S1 Fig [docx]
Genomic organization of human Activin B, C & E.

S2 Fig [sp]
Protein Sequences of full-length human Activins.

S3 Fig [arg]
Border region between prodomain (light gray) and active ligand (dark gray) of Activins.

S4 Fig [2015]
Alignment and structure of human Activin ligands.

S5 Fig [docx]
Sanger sequencing results.

S1 Video [mp4]
ActA X1 mp4.

S2 Video [mp4]
ActA X2 mp4.

S1 Raw Images [pdf]


Zdroje

1. Makanji Y, Zhu J, Mishra R, Holmquist C, Wong WPS, Schartz NB, et al. Inhibin at 90: from discovery to clinical application, a historical review. Endocrine Rev. 2014;35:747–94.

2. Salazar VS, Gamer LW, Rosen V. BMP signaling in skeletal development, disease and repair. Nat Rev Endocrinology. 2016;12:203–21. doi: 10.1038/nrendo.2016.12 26893264

3. Schwartz NB, Channing CP. Evidence for ovarian "inhibin": suppression of the secondary rise in serum stimulating hormone levels in proestrous rats by injection of porcine follicular fluid. Proc Natl Acad Sci USA. 1977;74:5721–4. doi: 10.1073/pnas.74.12.5721 271996

4. Vale W, Rivier J, Vaughan J, McClintock R, Corrigan A, Woo W, Karr D, et al. Purification and characterization of an FSH releasing protein from porcine ovarian follicular fluid. Nature. 1986;321:776–9. doi: 10.1038/321776a0 3012369

5. de Kretser DM, O’Hehir RE, Hardy CL, Hedger MP. The roles of activin A and its binding protein, follistatin, in inflammation and tissue repair. Mol Cell Endocrinol. 2012;359:101–6. doi: 10.1016/j.mce.2011.10.009 22037168

6. Jones KL, de Kretser DM, Patella S, Phillips DJ. Activin A and follistatin in systemic inflammation. Mol Cell Endocrinol. 2004;225:119–25. doi: 10.1016/j.mce.2004.07.010 15451576

7. Murry CE, Keller G. Differentiation of embryonic stem cells to clinically relevant populations: lessons from embryonic development. Cell. 2008;132:661–80. doi: 10.1016/j.cell.2008.02.008 18295582

8. Xia Y, Schneyer AL. The biology of activin: recent advances in structure, regulation and function. J Endocrinol. 2009;202:1–12. doi: 10.1677/JOE-08-0549 19273500

9. Antsiferova M, Werner S. The bright and dark sides of activin in wound healing and cancer. J Cell Sci. 2012;125:3929–37. doi: 10.1242/jcs.094789 22991378

10. Gray AM, Mason AJ. Requirement for activin A and transforming growth factor-beta 1 pro-regions in homodimer assembly. Science. 1990;247:1328–30. doi: 10.1126/science.2315700 2315700

11. Walton KL, Makanji Y, Harrison CA. New insights into the mechanisms of activin action and inhibition. Mol Cell Endocrinol. 2012;359:2–12. doi: 10.1016/j.mce.2011.06.030 21763751

12. Gentry LE, Lioubin MN, Purchio AF, Marquardt H. Molecular events in the processing of recombinant type I pre-pro-transforming growth factor beta to the mature polypeptide. Mol Cell Biol. 1988;8:4162–8. doi: 10.1128/mcb.8.10.4162 3185545

13. Mason AJ, Farnworth PG, Sullivan J. Characterization and determination of the biological activities of nonclevable high molecular wieght forms of inhibin A and activin A. Mol Endocrinol. 1996;10:1055–65. doi: 10.1210/mend.10.9.8885240 8885240

14. Harrington AE, Morris-Triggs SA, Ruotolo BT, Robinson CV, Ohnuma S, Hyvonen M. Structural basis for the inhibtion of activin signaling by follistatin. EMBO J. 2006;25:1035–45. doi: 10.1038/sj.emboj.7601000 16482217

15. McDowell N, Zorn AM, Crease DJ, Gurdon JB. Activin has direct long-range signaling activity and can form a concentration gradient by diffusion. Curr Biol. 1997;7:671–81. doi: 10.1016/s0960-9822(06)00294-6 9285724

16. Iozzo RV, Schaefer L. Proteoglycan form and function: a comprehensive nomenclature for proteoglycans. Matrix Biol. 2015;42:11–55. doi: 10.1016/j.matbio.2015.02.003 25701227

17. Nagaraj SH, Gasser RB, Ranganathan S. A hitchhiker’s guide to expressed sequence tag (EST) analysis. Brief Bioinform. 2007;8:6–21. doi: 10.1093/bib/bbl015 16772268

18. Wang X, Fischer G, Hyvonen M. Structure and activation of pro-activin A. Nat Communications. 2016;7:12052.

19. Li S, Shimono C, Norioka N, Nakano I, Okubo T, Yagi Y, et al. Activin A binds to perlecan through its pro-region that has heparin/heparan sulfate binding activity. J Biol Chem. 2010;285:36645–55. doi: 10.1074/jbc.M110.177865 20843788

20. Billings PC, Pacifici M. Interactions of signaling proteins, growth factors and other proteins with heparan sulfate: mechanisms and mysteries. Connect Tissue Res. 2015;56:272–80. doi: 10.3109/03008207.2015.1045066 26076122

21. Xu D, Esko JD. Demystifying heparan sulfate-protein interactions. Annu Rev Biochem. 2014;83:129–57. doi: 10.1146/annurev-biochem-060713-035314 24606135

22. Cardin AD, Weintraub HJ. Molecular modeling of protein-glycosaminoglycan interactions. Arterioschler Thromb Vasc Biol. 1989;9:21–32.

23. Fromm JR, Hileman RE, Caldwell EE, Weiler JM, Linhardt RJ. Pattern and spacing of basic amino acids in heparin binding sites. Arch Biochem Biophys. 1997;343:92–100. doi: 10.1006/abbi.1997.0147 9210650

24. Verrecchio A, Germann MW, Schick BP, Kung B, Twardowski T, San Antonio JD. Design of peptides with high affinities for heparin and endothelial cell proteoglycans. J Biol Chem. 2000;275:7701–7. doi: 10.1074/jbc.275.11.7701 10713081

25. Bulow HE, Hobert O. The molecular diversity of glycosaminoglycans shapes animal development. Annu Rev Cell Dev Biol. 2006;22:375–407. doi: 10.1146/annurev.cellbio.22.010605.093433 16805665

26. Sarrazin S, Lamanna WC, Esko JD. Heparan sulfate proteoglycans. Cold Spring Harb Perspec Biol. 2011;3:a004952.

27. Schultz V, Suflita M, Liu X, Zhang X, Yu Y, Li L, et al. Heparan sulfate domains required for fibroblast growth factor 1 and 2 signaling through fibroblast growth factor receptor 1c. J Biol Chem. 2017;292:2495–509. doi: 10.1074/jbc.M116.761585 28031461

28. Ohkawara B, Iemura S, ten Dijke P, Ueno N. Action range of BMP is defined by its N-terminal basic amino acid core. Curr Biol. 2002;12:205–9.

29. Ruppert R, Hoffmann E, Sebald W. Human bone morphogenetic protein 2 contains a heparin-binding site which modifies its biological activity. Eur J Biochem. 1996;237:295–302. doi: 10.1111/j.1432-1033.1996.0295n.x 8620887

30. Mundy C, Yang E, Takano H, Billings PC, Pacifici M. Heparan sulfate antagonism alters bone morphogenetic protein signaling and receptor dynamics, suggetsing a mechanism in Hereditary Multiple Exostoses. J Biol Chem. 2018; doi: 10.1074/jbc.RA117.000264 29622677

31. Billings PC, Yang E, Mundy C, Pacifici M. Domains with highest heparan sulfate-binding affinity reside at opposite ends in BMP2/4 versus BMP5/6/7: Implications for function. J Biol Chem. 2018;293:14371–83. doi: 10.1074/jbc.RA118.003191 30082319

32. Gandhi NS, Mancera RL. Prediction of heparin binding sites in bone morphogenetic proteins (BMPs). Biochim Biophys Acta. 2012;1824:1374–81. doi: 10.1016/j.bbapap.2012.07.002 22824487

33. Huegel J, Mundy C, Sgariglia F, Nygren P, Billings PC, Yamaguchi Y, et al. Perichondrium phenotype and border function are regulated by Ext1 and heparan sulfate in developing long bones: A mechanism likely deranged in Hereditary Multiple Exostoses. Dev Biol. 2013;377:100–12. doi: 10.1016/j.ydbio.2013.02.008 23458899

34. Huegel J, Sgariglia F, Enomoto-Iwamoto M, Koyama E, Dormans JP, Pacifici M. Heparan sulfate in skeletal development, growth, and pathology: the case of Hereditary Multiple Exostoses. Dev Dyn. 2013;242:1021–32. doi: 10.1002/dvdy.24010 23821404

35. Barton DE, Yang-Feng TL, Mason AJ, Seeburg PH, Francke U. Mapping of genes for inhibin subunits alpha, beta A, and beta B on human and mouse chromosomes and studies in jsd mice. Genomics. 1989;5:91–9. doi: 10.1016/0888-7543(89)90091-8 2767687

36. Mayo KE. Inhibin and activin: molecular aspects of regulation and function. Trends Endocrinol Metab. 1994;5:407–15. 18407237

37. Stewart AG, Milborrow HM, Ring JM, Crowther CE, Forage RG. Human inhibin genes. Genomic chracterisation and sequencing. FEBS Lett. 1986;206:329–34. doi: 10.1016/0014-5793(86)81006-7 3758355

38. Abe M, Shintani Y, Eto Y, Harada K, Kosaka M, Matsumoto T. Potent induction of activin A secretion from monocytes and bone marrow stromal fibroblasts by cogante interaction with activated T cells. J Leukoc Biol. 2002;72:347–52. 12149426

39. Sierra-Filardi E, Puig-Kroger A, Blanco FJ, Nieto C, Bragado R, Palomero MI, et al. Activin A skews macrophage polarization by promoting a proinflammatory phenotype and inhibiting the acquisition of anti-inflammatory macrophage markers. Blood. 2011;117:5092–101. doi: 10.1182/blood-2010-09-306993 21389328

40. Farokhzad OC, Teodoridis JM, Park HL, Arnaout MA, Shelley CS. CD43 gene expression is mediated by a nuclear factor which binds pyrimidine-rich single-stranded DNA. Nucleic Acids Res. 2000;28:2256–67. doi: 10.1093/nar/28.11.2256 10871347

41. Yu J, Shao L-E, Frigon NL, Lofgren J, Schwall R. Induced expression of a new cytokine, activin A, in human monocytes: inhibition by glucocorticoids and retinoic acid. Immunology. 1996;88:368–74. doi: 10.1046/j.1365-2567.1996.d01-675.x 8774352

42. Hefferan T, Subramaniam M, Hhosla S, Riggs BL, Spelsberg TC. Cytokine-specific induction of the TGF-beta inducible early gene (TIEG): regulation by specific members of the TGF-beta family. J Cell Biochem. 2000;78:380–90. 10861837

43. Yang J, Yan R, Roy A, Xu D, Poisson J, Zhang Y. The I-TASSER suite: protein structure and function prediction. Nat Methods. 2015;12:7–8. doi: 10.1038/nmeth.3213 25549265

44. Whalen DM, Malinauskas T, Gilbert RJC, Siebold C. Structural insights into proteoglycan-shaped hedgehog signaling. Proc Natl Acad Sci USA. 2013;110:16420–5. doi: 10.1073/pnas.1310097110 24062467

45. Capurro MI, Xu P, Shi W, Li F, Jia A, Filmus J. Glypican-3 inhibits hedgehog signaling during development by competing with patched for hedgehog binding. Dev Cell. 2008;14:700–11. doi: 10.1016/j.devcel.2008.03.006 18477453

46. Kirkpatrick CA, Knox SM, Staatz WD, Fox B, Lercher DM, Selleck SB. The function of a Drosophila glypican does not depend entirely on heparan sulfate modification. Dev Biol. 2006;300:570–82. doi: 10.1016/j.ydbio.2006.09.011 17055473

47. Sugino K, Kurosawa N, Nakamura T, Takio K, Shimasaki S, Ling N, et al. Molecular heterogeneity of follistatin, an activin-binding protein. Higher affinity of the carboxyl-terminal truncated forms for heparan sulphate proteoglyacns on the avorian granulosa cells. J Biol Chem. 1993;268:15579–87. 8340384

48. Ashikari-Hada S, Habuchi H, Kariya Y, Itoh N, Reddi AH, Kimata K. Characterization of growth factor-binding structures in heparin/heparan sulfate using an octasaccharide library. J Biol Chem. 2004;279:12346–54. doi: 10.1074/jbc.M313523200 14707131

49. Ye S, Luo Y, Lu W, Jones RB, Linhardt RJ, Capila I, et al. Structural basis for interaction of FGF-1, FGF-2 and FGF-7 with different heparan sulfate motifs. Biochemistry. 2001;40:14429–39. doi: 10.1021/bi011000u 11724555

50. Johnson KE, Makanji Y, Temple-Smith P, Kelly EK, Barton PA, Al-Musawi SL, et al. Biological activity and in vivo half-life of pro-activin A in male rats. Mol Cell Endocrinol. 2016;422:84–92. doi: 10.1016/j.mce.2015.12.007 26687063

51. Hagemann AI, Xu X, Nentwich O, Hyvonen M, Smith JC. Rab5-mediated endocytosis of activin is not required for gene activation or long-range signaling in Xenopous. Development. 2009;136:2803–13. doi: 10.1242/dev.034124 19605501

52. Hatsell SJ, Idone V, Alessi Wolken DM, Huang L, Kim HJ, Wang LC, et al. ACVR1R206H receptor mutation causes fibrodysplasia ossificans progressiva by imparting responsiveness to activin A. Science Trans Med. 2015;7:303ra137.

53. Shore E, Kaplan FS. Insights from a rare genetic disorder of extra-skeletal bone formation fibrodysplasia ossificans progressiva (FOP). Bone. 2008;43:427–33. doi: 10.1016/j.bone.2008.05.013 18590993

54. Shore E, Xu M, Feldman GJ, Fenstermacher DA, Consortium TFIR, Brown MA, et al. A recurrent mutation in the BMP type I receptor ACVR1 causes inherited and sporadic fibrodysplasia ossificans progressiva. Nature Genet. 2006;38:525–7. doi: 10.1038/ng1783 16642017

55. Pacifici M, Shore EM. Common mutations in ALK2/ACVR1, a multi-facet receptor, have roles in distinct musculoskeletal and neural orphan disorders. Cytokine & Growth Factor Rev. 2016;27:93–104.

56. Khachigian LM, Parish CR. Phosphomannopentoaose sulfate (PI-88): heparan sulfate mimetic with clinical potential in multiple vascular pathologies. Cardiovasc Drug Rev. 2004;22:1–6. 14978514

57. Dundas CM, Demonte D, Park S. Streptavidin-biotin technology: improvements and innovations in chemical and biological applications. Appl Microbiol Biotechnol. 2013;97:9343–53. doi: 10.1007/s00253-013-5232-z 24057405


Článek vyšel v časopise

PLOS One


2019 Číslo 9