Analysis of Procollagen C-Proteinase Enhancer-1/Glycosaminoglycan Binding Sites and of the Potential Role of Calcium Ions in the Interaction

PCPE2: Expression of multifunctional extracellular glycoprotein associated with diverse cellular functions

Procollagen C-endopeptidase enhancer 2, known as PCPE2 or PCOC2 (gene name, PCOLCE2) is a glycoprotein that resides in the extracellular matrix, and is similar in domain organization to PCPE1/PCPE, PCOC1 (PCOLCE1/PCOLCE). Due to the many similarities between the two related proteins, PCPE2 has been assumed to have biological functions similar to PCPE. PCPE is a well-established enhancer of procollagen processing activating the enzyme, BMP-1. However, reports show that PCPE2 has a strikingly different tissue expression profile compared to PCPE. With that in mind and given the paucity of published studies on PCPE2, this review examines the current literature citing PCPE2 and its association with specific cell types and signaling pathways. Additionally, this review will present a brief history of PCPE2's discovery, highlighting structural and functional similarities and differences compared to PCPE. Considering the widespread use of RNA sequencing techniques to examine associations between cell-specific gene expression and disease states, we will show that PCPE2 is repeatedly found as a differentially regulated gene (DEG) significantly associated with a number of cellular processes, well beyond the scope of procollagen fibril processing.

PCPE-2 (procollagen C-proteinase enhancer-2): the NON-IDENTICAL twin of PCPE-1

PCPE-2 was discovered at the beginning of this century, and was soon identified as a close homolog of PCPE-1 (procollagen C-proteinase enhancer 1). After the demonstration that it could also stimulate the proteolytic maturation of fibrillar procollagens by BMP-1/tolloid-like proteinases (BTPs), PCPE-2 did not attract much attention as it was thought to fulfill the same functions as PCPE-1 which was already well-described. However, the tissue distribution of PCPE-2 shows both common points and significant differences with PCPE-1, suggesting that their activities are not fully overlapping. Also, the recently established connections between PCPE-2 (gene name PCOLCE2) and several important diseases such as atherosclerosis, inflammatory diseases and cancer have highlighted the need for a thorough reappraisal of the in vivo roles of this regulatory protein. In this context, the recent finding that, while retaining the ability to bind fibrillar procollagens and to activate their C-terminal maturation, PCPE-2 can also bind BTPs and inhibit their activity has substantially extended its potential functions. In this review, we describe the current knowledge about PCPE-2 with a focus on collagen fibrillogenesis, lipid metabolism and inflammation, and discuss how we could further advance our understanding of PCPE-2-dependent biological processes.

In silico approaches for better understanding cysteine cathepsin-glycosaminoglycan interactions

Cysteine cathepsins constitute the largest cathepsin family, with 11 proteases in human that are present primarily within acidic endosomal and lysosomal compartments. They are involved in the turnover of intracellular and extracellular proteins. They are synthesized as inactive procathepsins that are converted to mature active forms. Cathepsins play important roles in physiological and pathological processes and, therefore, receive increasing attention as potential therapeutic targets. Their maturation and activity can be regulated by glycosaminoglycans (GAGs), long linear negatively charged polysaccharides composed of recurring dimeric units. In this review, we summarize recent computational progress in the field of (pro)cathepsin-GAG complexes analyses.

Conformational preferences of heparan sulfate to recognize the CXCL8 dimer in aqueous medium: degree of sulfation and hydrogen bonds

The sulfation pattern and epimerization of the long-chain sulfated polysaccharide heparan sulfate (HS) cause structural diversity and regulate various physiological and pathological processes when binding with proteins. In this work, we performed a series of molecular dynamics simulations of three variants of the octadecasaccharide HS with varying sulfation positions in aqueous medium in their free forms and in the presence of the chemokine CXCL8 dimer. The free energy of binding depicts the sulfation at the 6-O position of GlcNAc (HS6S), and both 3-O and 6-O positions of GlcNAc (HS3S6S) of HS variants are more likely to bind with the CXCL8 dimer than the triply sulfated HS2S3S6S, which is sulfated at the 2-O position of GlcUA additionally along with 3-O and 6-O positions of GlcNAc. Binding between HS and CXCL8 was driven by electrostatic and van der Waals interactions predominantly regardless of the sulfation pattern; however, unfavorable entropic contribution suppressed the interaction between HS and CXCL8. The contribution of different amino acid residues to the binding energetics suggested that basic amino acids line up the binding site of CXCL8. This study further acknowledges the role of interfacial water that is structured and bound with HS through hydrogen bonds, exhibiting differential hydrogen bond relaxation dynamics compared to when the HS molecules are free. Moreover, this study identifies that with the increase in sulfation, the HS-water hydrogen bond relaxation occurs faster with the complexation, while the reverse trend is followed in their free forms. Significant structural adaptation of the different sulfated HS molecules, as verified from the free energy landscapes generated from various reaction coordinates, root-mean-square-deviations, end-to-end distances, including ring pucker angles, dihedral flexibility, and the high conformational entropy cost arising from the glycosidic bonds, suggests that the different sulfated variants of HS undergo significant structural transformation to bind with CXCL8. The presence of a CXCL8 dimer imposes the bound forms of HS to adopt non-linear structures with skew-boat conformations. The atomistic details of the study would help in understanding the selectivity and conformational diversity, as well as the role of solvents in the recognition of CXCL8 by different sulfated variants of HS molecules.

Sulfation Effects of Chondroitin Sulfate to Bind a Chemokine in Aqueous Medium: Conformational Heterogeneity and Dynamics from Molecular Simulation

The sulfation patterns and degree of sulfation of chondroitin sulfate (CS), an important class of glycosaminoglycans (GAG), and their interactions with chemokines are accountable for various diseases. To realize the underlying mechanism of such complex biological phenomena at a molecular level and their application in rational drug design, a study on conformations and dynamics of CSs is necessary. To explore this, in this study, we performed a series of atomistic molecular dynamics (MD) simulations with different sulfated variants of octadecasaccharide CS, like CS-C, CS-E, and CS-T, in their free forms and when bound to the protein chemokine CXCL8 dimer in an aqueous medium. The calculated binding free energy of CSs with the CXCL8 dimer is favorable, and the degree of sulfation favors the complexation process further with prominent hydrophobic and hydrogen-bonded interactions. We find that the recognition is associated with the configurational entropy loss of the CS molecules as calculated from the Gaussian mixture approach, which supports that the degree of sulfation regulates the process. Cluster analysis through the k-means algorithm and end-to-end distance measurement revealed that although the free CS molecules adopted linear conformations, the nonlinear conformations during binding with protein were noted. Adaptation of nonlinear forms in the bound forms is noteworthy for the less-sulfated CS-C and CS-E. Apart from favorable 4C1 conformations, the occasional appearance of skew-boat forms from the free-energy map of ring pucker for the GlcUA unit was observed, which remains unaffected by the sulfation. We find that during recognition, the average relaxation time of intra-CS and inter-CS-CXCL8 hydrogen bonds (HBs) is about a magnitude lesser than that of CS-water HBs, most prominent on the involvement of higher sulfated CS-T analogues. The translational motion of surrounded water molecules in CSs exhibited sublinear diffusion, and the degree of sublinearity increases around the heavily sulfated molecules due to the hindrance created by them as well as the presence of the chemokine and exhibited markedly slow heterogeneous diffusion.

Macular Edema in Central Retinal Vein Occlusion Correlates With Aqueous Fibrinogen Alpha Chain

Purpose

The global protein profile of the aqueous humor has been found to correlate with the severity of retinal vascular disease. Studying the aqueous humor in central retinal vein occlusion (CRVO) with proteomic techniques may bring insights to the molecular mechanisms underlying the condition.

Methods

Aqueous humor samples from treatment naïve patients with CRVO complicated by macular edema (n = 28) and age-matched controls (n = 20) were analyzed by label-free quantification liquid chromatography - tandem mass spectrometry. Best corrected visual acuity (BCVA) was measured as logMAR, and the severity of macular edema was evaluated as central retinal thickness (CRT) with optical coherence tomography. Control samples were obtained prior to cataract surgery. Significantly changed proteins were identified by a permutation-based calculation with a false discovery rate of 0.05.

Results

A total of 177 proteins were differentially expressed in CRVO. Regulated proteins were involved in complement activation, innate immune response, blood coagulation, and cell adhesion. Upregulated proteins that correlated with BCVA and CRT included fibrinogen alpha, beta, and gamma chains, fibronectin, Ig lambda-6 chain C region, Ig alpha-1 chain C region, and complement C7. Downregulated proteins that correlated negatively with BCVA, and CRT, included procollagen C-endopeptidase enhancer 1, clusterin, opticin, reelin, fibrillin-1, and cadherin-2. Monocyte differentiation antigen CD14 and lipopolysaccharide-binding protein were increased in CRVO.

Conclusions

Fibrinogen chains, fibronectin, and immunoglobulin components correlated with BCVA and CRT, suggesting a multifactorial response. Protective anti-angiogenic proteins, including procollagen C-endopeptidase enhancer 1, clusterin, and opticin, were downregulated in CRVO and correlated negatively with BCVA and CRT.

Impact of calcium ions on the structural and dynamic properties of heparin oligosaccharides by computational analysis

Heparin (HP) belongs to glycosaminoglycans (GAGs), anionic linear polysaccharides composed of repetitive disaccharide units. They are key players in many biological processes occurring in the extracellular matrix and at the cell surface. GAGs are challenging molecules for computational research due to their high chemical heterogeneity, flexibility, periodicity, pseudosymmetry, predominantly electrostatics-driven nature of interactions with their protein partners and potential multipose binding. The molecular mechanisms underlying GAG interactions mediated by divalent ions, which are important for GAG binding to several proteins, are not well understood. The goal of this study was to characterize the binding of Ca²⁺ to two HP oligosaccharides of different lengths (dp10 and dp18, dp: degree of polymerization) and their impact on HP conformational space and their dynamic behavior with the use of molecular dynamics (MD)-based approaches with two Ca²⁺ parameter sets. MD data suggested that the flexibility of the monosaccharides, the glycosidic linkages and ring puckering were not affected by the presence of Ca²⁺, in contrast to H-bond propensities and the calculated R_g for a fraction of the oligosaccharide populations in both dp10 and dp18. Moreover, the essential differences in the data obtained by using two Ca²⁺ parameter sets were reported.

Pcpe2, a Novel Extracellular Matrix Protein, Regulates Adipocyte SR-BI-Mediated High-Density Lipoprotein Uptake

Objective

To investigate the role of adipocyte Pcpe2 (procollagen C-endopeptidase enhancer 2) in SR-BI (scavenger receptor class BI)-mediated HDL-C (high-density lipoprotein cholesterol) uptake and contributions to adipose lipid storage.

Approach and Results

Pcpe2, a glycoprotein devoid of intrinsic proteolytic activity, is believed to participate in extracellular protein-protein interactions, supporting SR-BI- mediated HDL-C uptake. In published studies, Pcpe2 deficiency increased the development of atherosclerosis by reducing SR-BI-mediated HDL-C catabolism, but the biological impact of this deficiency on adipocyte SR-BI-mediated HDL-C uptake is unknown. Differentiated cells from Ldlr-/-/Pcpe2-/- (Pcpe2-/-) mouse adipose tissue showed elevated SR-BI protein levels, but significantly reduced HDL-C uptake compared to Ldlr-/- (control) adipose tissue. SR-BI-mediated HDL-C uptake was restored by preincubation of cells with exogenous Pcpe2. In diet-fed mice lacking Pcpe2, significant reductions in visceral, subcutaneous, and brown adipose tissue mass were observed, despite elevations in plasma triglyceride and cholesterol concentrations. Significant positive correlations exist between adipose mass and Pcpe2 expression in both mice and humans.

Conclusions

Overall, these findings reveal a novel and unexpected function for Pcpe2 in modulating SR-BI expression and function as it relates to adipose tissue expansion and cholesterol balance in both mice and humans.

Effect of Tetraphenylborate on Physicochemical Properties of Bovine Serum Albumin

The binding interactions of bovine serum albumin (BSA) with tetraphenylborate ions ([B(Ph)4]-) have been investigated by a set of experimental methods (isothermal titration calorimetry, steady-state fluorescence spectroscopy, differential scanning calorimetry and circular dichroism spectroscopy) and molecular dynamics-based computational approaches. Two sets of structurally distinctive binding sites in BSA were found under the experimental conditions (10 mM cacodylate buffer, pH 7, 298.15 K). The obtained results, supported by the competitive interactions experiments of SDS with [B(Ph)4]- for BSA, enabled us to find the potential binding sites in BSA. The first site is located in the subdomain I A of the protein and binds two [B(Ph)4]- ions (logK(ITC)1 = 7.09 ± 0.10; ΔG(ITC)1 = -9.67 ± 0.14 kcal mol-1; ΔH(ITC)1 = -3.14 ± 0.12 kcal mol-1; TΔS(ITC)1 = -6.53 kcal mol-1), whereas the second site is localized in the subdomain III A and binds five ions (logK(ITC)2 = 5.39 ± 0.06; ΔG(ITC)2 = -7.35 ± 0.09 kcal mol-1; ΔH(ITC)2 = 4.00 ± 0.14 kcal mol-1; TΔS(ITC)2 = 11.3 kcal mol-1). The formation of the {[B(Ph)4]-}-BSA complex results in an increase in the thermal stability of the alfa-helical content, correlating with the saturation of the particular BSA binding sites, thus hindering its thermal unfolding.

Investigation of the structure of regulatory proteins interacting with glycosaminoglycans by combining NMR spectroscopy and molecular modeling - the beginning of a wonderful friendship

The interaction of regulatory proteins with extracellular matrix or cell surface-anchored glycosaminoglycans (GAGs) plays important roles in molecular recognition, wound healing, growth, inflammation and many other processes. In spite of their high biological relevance, protein-GAG complexes are significantly underrepresented in structural databases because standard tools for structure determination experience difficulties in studying these complexes. Co-crystallization with subsequent X-ray analysis is hampered by the high flexibility of GAGs. NMR spectroscopy experiences difficulties related to the periodic nature of the GAGs and the sparse proton network between protein and GAG with distances that typically exceed the detection limit of nuclear Overhauser enhancement spectroscopy. In contrast, computer modeling tools have advanced over the last years delivering specific protein-GAG docking approaches successfully complemented with molecular dynamics (MD)-based analysis. Especially the combination of NMR spectroscopy in solution providing sparse structural constraints with molecular docking and MD simulations represents a useful synergy of forces to describe the structure of protein-GAG complexes. Here we review recent methodological progress in this field and bring up examples where the combination of new NMR methods along with cutting-edge modeling has yielded detailed structural information on complexes of highly relevant cytokines with GAGs.

Evaluation of replica exchange with repulsive scaling approach for docking glycosaminoglycans

Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.

Further analyses of APRIL/APRIL-Receptor/Glycosaminoglycan interactions by biochemical assays linked to computational studies

A proliferation-inducing ligand (APRIL) is a member of the tumor necrosis factor superfamily. APRIL is quite unique in this superfamily for at least for two reasons: i) it binds to glycosaminoglycans (GAGs) via its positively charged N-terminus; ii) one of its signaling receptor, the transmembrane activator CAML interactor (TACI) was also reported to bind GAGs. Here, as provided by biochemical evidences with the use of an APRIL deletion mutant linked to computational studies, APRIL-GAG interaction involved other regions than the APRIL N-terminus. Preferential interaction of APRIL with heparin followed by chondroitin sulfate E were confirmed by in silico analysis. Both computational and experimental approaches did not reveal heparan sulfate binding to TACI. Together, computational results corroborated experiments contributing with atomistic details to the knowledge on this biologically relevant trimolecular system. Additionally, a high-throughput rigorous analysis of the free energy calculations data was performed to critically evaluate the applied computational methodologies.

Computational insights into the role of calcium ions in protein–glycosaminoglycan systems

The prediction power of computational methodologies for studying the role of ions in protein–glycosaminoglycan interactions was critically assessed.

Glycosaminoglycan-Protein Interactions: The First Draft of the Glycosaminoglycan Interactome

The six mammalian glycosaminoglycans (GAGs), chondroitin sulfate, dermatan sulfate, heparin, heparan sulfate, hyaluronan, and keratan sulfate, are linear polysaccharides. Except for hyaluronan, they are sulfated to various extent, and covalently attached to proteins to form proteoglycans. GAGs interact with growth factors, morphogens, chemokines, extracellular matrix proteins and their bioactive fragments, receptors, lipoproteins, and pathogens. These interactions mediate their functions, from embryonic development to extracellular matrix assembly and regulation of cell signaling in various physiological and pathological contexts such as angiogenesis, cancer, neurodegenerative diseases, and infections. We give an overview of GAG-protein interactions (i.e., specificity and chemical features of GAG- and protein-binding sequences), and review the available GAG-protein interaction networks. We also provide the first comprehensive draft of the GAG interactome composed of 832 biomolecules (827 proteins and five GAGs) and 932 protein-GAG interactions. This network is a scaffold, which in the future should integrate structures of GAG-protein complexes, quantitative data of the abundance of GAGs in tissues to build tissue-specific interactomes, and GAG interactions with metal ions such as calcium, which plays a major role in the assembly of the extracellular matrix and its interactions with cells. This contextualized interactome will be useful to identify druggable GAG-protein interactions for therapeutic purpose.