SiMPLOD, a Structure-Integrated Database of Collagen Lysyl Hydroxylase (LH/PLOD) Enzyme Variants

Enzyme Databases in the Era of Omics and Artificial Intelligence

Enzyme research is important for the development of various scientific fields such as medicine and biotechnology. Enzyme databases facilitate this research by providing a wide range of information relevant to research planning and data analysis. Over the years, various databases that cover different aspects of enzyme biology (e.g., kinetic parameters, enzyme occurrence, and reaction mechanisms) have been developed. Most of the databases are curated manually, which improves reliability of the information; however, such curation cannot keep pace with the exponential growth in published data. Lack of data standardization is another obstacle for data extraction and analysis. Improving machine readability of databases is especially important in the light of recent advances in deep learning algorithms that require big training datasets. This review provides information regarding the current state of enzyme databases, especially in relation to the ever-increasing amount of generated research data and recent advancements in artificial intelligence algorithms. Furthermore, it describes several enzyme databases, providing the reader with necessary information for their use.

Unleashing the Potential of 1,3-Diketone Analogues as Selective LH2 Inhibitors

Lysyl hydroxylase 2 (LH2) catalyzes the formation of highly stable hydroxylysine aldehyde-derived collagen cross-links (HLCCs), thus promoting lung cancer metastasis through its capacity to modulate specific types of collagen cross-links within the tumor stroma. Using 1 and 2 from our previous high-throughput screening (HTS) as lead probes, we prepared a series of 1,3-diketone analogues, 1-18, and identified 12 and 13 that inhibit LH2 with IC50's of approximately 300 and 500 nM, respectively. Compounds 12 and 13 demonstrate selectivity for LH2 over LH1 and LH3. Quantum mechanics/molecular mechanics (QM/MM) modeling indicates that the selectivity of 12 and 13 may stem from noncovalent interactions like hydrogen bonding between the morpholine/piperazine rings with the LH2-specific Arg661. Treatment of 344SQ WT cells with 13 resulted in a dose-dependent reduction in their migration potential, whereas the compound did not impede the migration of the same cell line with an LH2 knockout (LH2KO).

Identification of Regulatory Molecular "Hot Spots" for LH/PLOD Collagen Glycosyltransferase Activity

Hydroxylysine glycosylations are post-translational modifications (PTMs) essential for the maturation and homeostasis of fibrillar and non-fibrillar collagen molecules. The multifunctional collagen lysyl hydroxylase 3 (LH3/PLOD3) and the collagen galactosyltransferase GLT25D1 are the human enzymes that have been identified as being responsible for the glycosylation of collagen lysines, although a precise description of the contribution of each enzyme to these essential PTMs has not yet been provided in the literature. LH3/PLOD3 is thought to be capable of performing two chemically distinct collagen glycosyltransferase reactions using the same catalytic site: an inverting beta-1,O-galactosylation of hydroxylysines (Gal-T) and a retaining alpha-1,2-glucosylation of galactosyl hydroxylysines (Glc-T). In this work, we have combined indirect luminescence-based assays with direct mass spectrometry-based assays and molecular structure studies to demonstrate that LH3/PLOD3 only has Glc-T activity and that GLT25D1 only has Gal-T activity. Structure-guided mutagenesis confirmed that the Glc-T activity is defined by key residues in the first-shell environment of the glycosyltransferase catalytic site as well as by long-range contributions from residues within the same glycosyltransferase (GT) domain. By solving the molecular structures and characterizing the interactions and solving the molecular structures of human LH3/PLOD3 in complex with different UDP-sugar analogs, we show how these studies could provide insights for LH3/PLOD3 glycosyltransferase inhibitor development. Collectively, our data provide new tools for the direct investigation of collagen hydroxylysine PTMs and a comprehensive overview of the complex network of shapes, charges, and interactions that enable LH3/PLOD3 glycosyltransferase activities, expanding the molecular framework and facilitating an improved understanding and manipulation of glycosyltransferase functions in biomedical applications.

Modeling of Protein Complexes

The recent advances in structural biology, combined with continuously increasing computational capabilities and development of advanced softwares, have drastically simplified the workflow for protein homology modeling. Modeling of individual proteins is nowadays quick and straightforward for a large variety of protein targets, thanks to guided pipelines relying on advanced computational tools and user-friendly interfaces, which have extended and promoted the use of modeling also to scientists not focusing on molecular structures of proteins. Nevertheless, construction of models of multi-protein complexes remains quite challenging for the non-experts, often due to the usage of specific procedures depending on the system under investigation and the need for experimental validation approaches to strengthen the generated output.In this chapter, we provide a brief overview of the approaches enabling generation of multi-protein complex models starting from homology models of individual protein components. Using real-life examples, we include two examples to guide the reader in the generation of homomeric and heteromeric protein models.

Comparative genomic and biochemical analyses identify a collagen galactosylhydroxylysyl glucosyltransferase from Acanthamoeba polyphaga mimivirus

Humans and Acanthamoeba polyphaga mimivirus share numerous homologous genes, including collagens and collagen-modifying enzymes. To explore this homology, we performed a genome-wide comparison between human and mimivirus using DELTA-BLAST (Domain Enhanced Lookup Time Accelerated BLAST) and identified 52 new putative mimiviral proteins that are homologous with human proteins. To gain functional insights into mimiviral proteins, their human protein homologs were organized into Gene Ontology (GO) and REACTOME pathways to build a functional network. Collagen and collagen-modifying enzymes form the largest subnetwork with most nodes. Further analysis of this subnetwork identified a putative collagen glycosyltransferase R699. Protein expression test suggested that R699 is highly expressed in Escherichia coli, unlike the human collagen-modifying enzymes. Enzymatic activity assay and mass spectrometric analyses showed that R699 catalyzes the glucosylation of galactosylhydroxylysine to glucosylgalactosylhydroxylysine on collagen using uridine diphosphate glucose (UDP-glucose) but no other UDP-sugars as a sugar donor, suggesting R699 is a mimiviral collagen galactosylhydroxylysyl glucosyltransferase (GGT). To facilitate further analysis of human and mimiviral homologous proteins, we presented an interactive and searchable genome-wide comparison website for quickly browsing human and Acanthamoeba polyphaga mimivirus homologs, which is available at RRID Resource ID: SCR_022140 or https://guolab.shinyapps.io/app-mimivirus-publication/ .

Overexpressing PLOD Family Genes Predict Poor Prognosis in Pancreatic Cancer

Background

Pancreatic cancer is a common malignant tumor. Multiple studies have shown that procollagen lysyl-hydroxylase (PLOD) family genes were closely related to tumor progression and metastasis in a variety of human cancers. This study aimed to explore the prognosis and biological role of PLOD family genes in pancreatic adenocarcinoma (PAAD).

Methods

GEPIA, GEO, HPA, CCLE, Kaplan-Meier plotter, cBioPortal, LinkedOmics, DAVID6.8, STRING, and TIMER were employed to determine the prognostic values and biological function of PLOD family members in PAAD.

Results

The mRNA and protein expression patterns of PLOD family members were noticeably up-regulated in PAAD compared with normal tissues. PLOD family gene expression was also up-regulated in pancreatic cancer cell lines. PLOD1 was correlated with histological and pathological grades of pancreatic cancer. PLOD2 was related to histological grade. The high expression of PLOD1-2 was correlated with the poor overall survival rate and relapse-free survival rate in patients with PAAD. Additionally, PLODs showed high sensitivity and specificity in distinguishing pancreatic cancer from normal tissues. Through the functional enrichment analysis of PLOD-related genes in PAAD, we found that PLODs were enriched in collagen fiber tissue structure, lysine degradation, and collagen biosynthesis. Pathway analysis confirmed that PLODs regulated the proliferation, migration, and metastasis of pancreatic cancer through the RalGEF-Ral signaling pathway. Furthermore, the level of expression of PLOD1-2 was positively correlated with the activity of tumor-infiltrating immune cells, including CD8+T cells, neutrophils, macrophages, and dendritic cells. The level of expression of PLOD3 was inversely correlated with the level of infiltration of CD8+T cells. PLOD1 and PLOD2 were highly expressed in pancreatic cancer tissues with TP53 and KRAS mutations, respectively. However, the level of expression of PLOD3 in SMAD4 wild-type pancreatic cancer was increased.

Conclusion

The findings showed that individual PLOD genes or PLOD family genes could be potential prognostic biomarkers for PAAD.

PLOD1 acts as a tumor promoter in glioma via activation of the HSF1 signaling pathway

Procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 (PLOD1) is a collagen-related lysyl hydroxylase and its prognostic value in glioma patients was verified. However, its biological function in glioma has yet to be fully investigated. The PLOD1 mRNA status and clinical significance in gliomas were assessed via the GEPIA database. Overexpression or targeted depletion of PLOD1 was carried out in the human glioma cell line U87 and verified by western blotting. CCK8 and colony formation assays were implemented to examine the impact of PLOD1 on the proliferative and colony-forming phenotypes of U87 cells. Luciferase reporter assays and HSF1-specific pharmacologic inhibitors (KRIBB11) were employed to determine the regulatory relationship between PLOD1 and heat shock factor 1 (HSF1). High expression of PLOD1 was observed in tissue samples of glioblastoma multiforme (GBM) and brain lower-grade glioma (LGG). GEPIA overall survival further demonstrated that both GBM and LGG patients with high PLOD1 displayed worse clinical outcomes compared with those with low PLOD1. Overexpression and targeted depletion of PLOD1 enhanced and suppressed U87 cell proliferation and colony formation, respectively. Luciferase reporter assays showed that PLOD1 significantly enhanced the transcriptional activity of HSF1 in HEK293T cells. PLOD1 deficiency in U87 cells inhibited HSF1-induced survivin accumulation, whereas KRIBB11 also blocked the PLOD1-overexpressing induced survivin expression. An inhibitor of HSF1 signaling events abolished the increased clonogenic potential caused by PLOD1 overexpression in U87 cells. High expression of PLOD1 can increase the proliferation and colony formation of U87 cells by activating the HSF1 signaling pathway. This study suggested PLOD1/HSF1 as an effective therapeutic target for gliomas.

New mechanistic insights to PLOD1-mediated human vascular disease

Heritable thoracic aortic disease and familial thoracic aortic aneurysm/dissection are important causes of human morbidity/mortality, most without identifiable genetic cause. In a family with familial thoracic aortic aneurysm/dissection, we identified a missense p. (Ser178Arg) variant in PLOD1 segregating with disease, and evaluated PLOD1 enzymatic activity, collagen characteristics and in human aortic vascular smooth muscle cells, studied the effect on function. Comparison with homologous PLOD3 enzyme indicated that the pathogenic variant may affect the N-terminal glycosyltransferase domain, suggesting unprecedented PLOD1 activity. In vitro assays demonstrated that wild-type PLOD1 is capable of processing UDP-glycan donor substrates, and that the variant affects the folding stability of the glycosyltransferase domain and associated enzymatic functions. The PLOD1 substrate lysine was elevated in the proband, however the enzymatic product hydroxylysine and total collagen content was not different, albeit despite collagen fibril narrowing and preservation of collagen turnover. In VSMCs overexpressing wild-type PLOD1, there was upregulation in procollagen gene expression (secretory function) which was attenuated in the variant, consistent with loss-of-function. In comparison, si-PLOD1 cells demonstrated hypercontractility and upregulation of contractile markers, providing evidence for phenotypic switching. Together, the findings suggest that the PLOD1 product is preserved, however newly identified glucosyltransferase activity of PLOD1 appears to be affected by folding stability of the variant, and is associated with compensatory vascular smooth muscle cells phenotypic switching to support collagen production, albeit with less robust fibril girth. Future studies should focus on the impact of PLOD1 folding/variant stability on the tertiary structure of collagen and ECM interactions.

A Human Pan-Cancer System Analysis of Procollagen-Lysine, 2-Oxoglutarate 5-Dioxygenase 3 (PLOD3)

The overexpression of the enzymes involved in the degradation of procollagen lysine is correlated with various tumor entities. Procollagen-lysine, 2-oxoglutarate 5-dioxygenase 3 (PLOD3) expression was found to be correlated to the progression and migration of cancer cells in gastric, lung and prostate cancer. Here, we analyzed the gene expression, protein expression, and the clinical parameters of survival across 33 cancers based on the Clinical Proteomic Tumor Analysis Consortium (CPTAC), function annotation of the mammalian genome 5 (FANTOM5), Gene Expression Omnibus (GEO), Genotype-Tissue Expression (GTEx), Human Protein Atlas (HPA) and The Cancer Genome Atlas (TCGA) databases. Genetic alteration, immune infiltration and relevant cellular pathways were analyzed in detail. PLOD3 expression negatively correlated with survival periods and the infiltration level of CD8⁺ T cells, but positively correlated to the infiltration of cancer associated fibroblasts in diverse cancers. Immunohistochemistry in colon carcinomas, glioblastomas, and soft tissue sarcomas further confirm PLOD 3 expression in human cancer tissue. Moreover, amplification and mutation accounted for the largest proportion in esophageal adenocarcinoma and uterine corpus endometrial carcinoma, respectively; the copy number alteration of PLOD3 appeared in all cancers from TCGA; and molecular mechanisms further proved the effect of PLOD3 on tumorigenesis. In particular, PLOD3 expression appears to have a tumor immunological effect, and is related to multiple immune cells. Furthermore, it is also associated with tumor mutation burden and microsatellite instability in various tumors. PLOD3 acts as an inducer of various cancers, and it could be a potential biomarker for prognosis and targeted treatment.

A collagen glucosyltransferase drives lung adenocarcinoma progression in mice

Cancer cells are a major source of enzymes that modify collagen to create a stiff, fibrotic tumor stroma. High collagen lysyl hydroxylase 2 (LH2) expression promotes metastasis and is correlated with shorter survival in lung adenocarcinoma (LUAD) and other tumor types. LH2 hydroxylates lysine (Lys) residues on fibrillar collagen's amino- and carboxy-terminal telopeptides to create stable collagen cross-links. Here, we show that electrostatic interactions between the LH domain active site and collagen determine the unique telopeptidyl lysyl hydroxylase (tLH) activity of LH2. However, CRISPR/Cas-9-mediated inactivation of tLH activity does not fully recapitulate the inhibitory effect of LH2 knock out on LUAD growth and metastasis in mice, suggesting that LH2 drives LUAD progression, in part, through a tLH-independent mechanism. Protein homology modeling and biochemical studies identify an LH2 isoform (LH2b) that has previously undetected collagen galactosylhydroxylysyl glucosyltransferase (GGT) activity determined by a loop that enhances UDP-glucose-binding in the GLT active site and is encoded by alternatively spliced exon 13 A. CRISPR/Cas-9-mediated deletion of exon 13 A sharply reduces the growth and metastasis of LH2b-expressing LUADs in mice. These findings identify a previously unrecognized collagen GGT activity that drives LUAD progression.

Collagen hydroxylysine glycosylation: non-conventional substrates for atypical glycosyltransferase enzymes

Collagen is a major constituent of the extracellular matrix (ECM) that confers fundamental mechanical properties to tissues. To allow proper folding in triple-helices and organization in quaternary super-structures, collagen molecules require essential post-translational modifications (PTMs), including hydroxylation of proline and lysine residues, and subsequent attachment of glycan moieties (galactose and glucose) to specific hydroxylysine residues on procollagen alpha chains. The resulting galactosyl-hydroxylysine (Gal-Hyl) and less abundant glucosyl-galactosyl-hydroxylysine (Glc-Gal-Hyl) are amongst the simplest glycosylation patterns found in nature and are essential for collagen and ECM homeostasis. These collagen PTMs depend on the activity of specialized glycosyltransferase enzymes. Although their biochemical reactions have been widely studied, several key biological questions about the possible functions of these essential PTMs are still missing. In addition, the lack of three-dimensional structures of collagen glycosyltransferase enzymes hinders our understanding of the catalytic mechanisms producing this modification, as well as the impact of genetic mutations causing severe connective tissue pathologies. In this mini-review, we summarize the current knowledge on the biochemical features of the enzymes involved in the production of collagen glycosylations and the current state-of-the-art methods for the identification and characterization of this important PTM.

Comprehensive analysis of PLOD family members in low-grade gliomas using bioinformatics methods

Low-grade gliomas (LGGs) is a primary invasive brain tumor that grows slowly but is incurable and eventually develops into high malignant glioma. Novel biomarkers for the tumorigenesis and lifetime of LGG are critically demanded to be investigated. In this study, the expression levels of procollagen-lysine, 2-oxoglutarate 5-dioxygenases (PLODs) were analyzed by ONCOMINE, HPA and GEPIA. The GEPIA online platform was applied to evaluate the interrelation between PLODs and survival index in LGG. Furthermore, functions of PLODs and co-expression genes were inspected by the DAVID. Moreover, we used TIMER, cBioportal, GeneMINIA and NetworkAnalyst analysis to reveal the mechanism of PLODs in LGG. We found that expression levels of each PLOD family members were up-regulated in patients with LGG. Higher expression of PLODs was closely related to shorter disease-free survival (DFS) and overall survival (OS). The findings showed that LGG cases with or without alterations were significantly correlated with the OS and DFS. The mechanism of PLODs in LGG may be involved in response to hypoxia, oxidoreductase activity, Lysine degradation and immune cell infiltration. In general, this research has investigated the values of PLODs in LGG, which could serve as biomarkers for diagnosis, prognosis and potential therapeutic targets of LGG patients.

Prolyl and lysyl hydroxylases in collagen synthesis

Collagens are the most abundant proteins in the extracellular matrix. They provide a framework to build organs and tissues and give structural support to make them resistant to mechanical load and forces. Several intra- and extracellular modifications are needed to make functional collagen molecules, intracellular post-translational modifications of proline and lysine residues having key roles in this. In this article, we provide a review on the enzymes responsible for the proline and lysine modifications, that is collagen prolyl 4-hydroxylases, 3-hydroxylases and lysyl hydroxylases, and discuss their biological functions and involvement in diseases.

ATDC5 cells as a model of cartilage extracellular matrix neosynthesis, maturation and assembly

Fibrillar collagens and proteoglycans (PGs) are quantitatively the major constituents of extracellular matrices (ECM). They carry numerous crucial post-translational modifications (PTMs) that tune the resulting biomechanical properties of the corresponding tissues. The mechanisms determining these PTMs remain largely unknown, notably because available established cell lines do not recapitulate much of the complexity of the machineries involved. ATDC5 cells are a model of chondrogenesis widely used for decades, but it remains described mostly at histological and transcriptional levels. Here, we asked to what extent this model recapitulates the events of ECM synthesis and processing occurring in cartilage. Insulin-stimulated ATDC5 cells exhibit up- or down-regulation of more than one-hundred proteins, including a number of known participants in chondrogenesis and major markers thereof. However, they also lack several ECM components considered of significant, yet more subtle, function in cartilage. Still, they assemble the large PG aggrecan and type II collagen, both carrying most of their in vivo PTMs, into an ECM. Remarkably, collagen crosslinking is fully lysyl oxidase (LOX)-dependent. The ATDC5 model recapitulates critical aspects of the cartilage ECM-processing machinery and should be useful to decipher the mechanisms involved. Proteomics data are available via ProteomeXchange with identifier PXD014121. SIGNIFICANCE: The present work provides the first proteome characterization of the ATDC5 chondrogenesis model, which has been used for decades in the field of cartilage biology. The results demonstrate the up- and down-regulation of more than one hundred proteins. Overall, specific drawbacks of the model are pointed out, that will be important to take into consideration for future studies. However, major cartilage components are massively assembled into an extracellular matrix and carry most of their post-translational modifications occurring in cartilage tissue. Unlike other available established cell lines, the ATDC5 model recapitulates major aspects of cartilage biosynthesis and should be useful in investigating the mechanisms that regulate collagen maturation events.

Fibrosis and cancer: A strained relationship

Tumors are characterized by extracellular matrix (ECM) deposition, remodeling, and cross-linking that drive fibrosis to stiffen the stroma and promote malignancy. The stiffened stroma enhances tumor cell growth, survival and migration and drives a mesenchymal transition. A stiff ECM also induces angiogenesis, hypoxia and compromises anti-tumor immunity. Not surprisingly, tumor aggression and poor patient prognosis correlate with degree of tissue fibrosis and level of stromal stiffness. In this review, we discuss the reciprocal interplay between tumor cells, cancer associated fibroblasts (CAF), immune cells and ECM stiffness in malignant transformation and cancer aggression. We discuss CAF heterogeneity and describe its impact on tumor development and aggression focusing on the role of CAFs in engineering the fibrotic tumor stroma and tuning tumor cell tension and modulating the immune response. To illustrate the role of mechanoreciprocity in tumor evolution we summarize data from breast cancer and pancreatic ductal carcinoma (PDAC) studies, and finish by discussing emerging anti-fibrotic strategies aimed at treating cancer.

Collagen cross-linking mediated by lysyl hydroxylase 2: an enzymatic battlefield to combat fibrosis

The hallmark of fibrosis is an excessive accumulation of collagen, ultimately leading to organ failure. It has become evident that the deposited collagen also exhibits qualitative modifications. A marked modification is the increased cross-linking, leading to a stabilization of the collagen network and limiting fibrosis reversibility. Not only the level of cross-linking is increased, but also the composition of cross-linking is altered: an increase is seen in hydroxyallysine-derived cross-links at the expense of allysine cross-links. This results in irreversible fibrosis, as collagen cross-linked by hydroxyallysine is more difficult to degrade. Hydroxyallysine is derived from a hydroxylysine in the telopeptides of collagen. The expression of lysyl hydroxylase (LH) 2 (LH2), the enzyme responsible for the formation of telopeptidyl hydroxylysine, is universally up-regulated in fibrosis. It is expected that inhibition of this enzyme will lead to reversible fibrosis without interfering with the normal repair process. In this review, we discuss the molecular basis of collagen modifications and cross-linking, with an emphasis on LH2-mediated hydroxyallysine cross-links, and their implications for the pathogenesis and treatment of fibrosis.

Pathogenic variants in PLOD3 result in a Stickler syndrome-like connective tissue disorder with vascular complications

Background

Pathogenic PLOD3 variants cause a connective tissue disorder (CTD) that has been described rarely. We further characterise this CTD and propose a clinical diagnostic label to improve recognition and diagnosis of PLOD3-related disease.

Methods

Reported PLOD3 phenotypes were compared with known CTDs utilising data from three further individuals from a consanguineous family with a homozygous PLOD3 c.809C>T; p.(Pro270Leu) variant. PLOD3 mRNA expression in the developing embryo was analysed for tissue-specific localisation. Mouse microarray expression data were assessed for phylogenetic gene expression similarities across CTDs with overlapping clinical features.

Results

Key clinical features included ocular abnormalities with risk for retinal detachment, sensorineural hearing loss, reduced palmar creases, finger contractures, prominent knees, scoliosis, low bone mineral density, recognisable craniofacial dysmorphisms, developmental delay and risk for vascular dissection. Collated clinical features showed most overlap with Stickler syndrome with variable features of Ehlers-Danlos syndrome (EDS) and epidermolysis bullosa (EB). Human lysyl hydroxylase 3/PLOD3 expression was localised to the developing cochlea, eyes, skin, forelimbs, heart and cartilage, mirroring the clinical phenotype of this disorder.

Conclusion

These data are consistent with pathogenic variants in PLOD3 resulting in a clinically distinct Stickler-like syndrome with vascular complications and variable features of EDS and EB. Early identification of PLOD3 variants would improve monitoring for comorbidities and may avoid serious adverse ocular and vascular outcomes.