Novel Fragment Inhibitors of PYCR1 from Docking-Guided X-ray Crystallography.

The proline biosynthetic enzyme Δ1-pyrroline-5-carboxylate (P5C) reductase 1 (PYCR1) is one of the most consistently upregulated enzymes across multiple cancer types and central to the metabolic rewiring of cancer cells. Herein, we describe a fragment-based, structure-first approach to the discovery of PYCR1 inhibitors. Thirty-seven fragment-like carboxylic acids in the molecular weight range of 143-289 Da were selected from docking and then screened using X-ray crystallography as the primary assay. Strong electron density was observed for eight compounds, corresponding to a crystallographic hit rate of 22%. The fragments are novel compared to existing proline analog inhibitors in that they block both the P5C substrate pocket and the NAD(P)H binding site. Four hits showed inhibition of PYCR1 in kinetic assays, and one has lower apparent IC50 than the current best proline analog inhibitor. These results show proof-of-concept for our inhibitor discovery approach and provide a basis for fragment-to-lead optimization.

Expression and kinetic characterization of PYCR3.

PYCRs are proline biosynthetic enzymes that catalyze the NAD(P)H-dependent reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline in humans. PYCRs - especially PYCR1 - are upregulated in many types of cancers and have been implicated in the altered metabolism of cancer cells. Of the three isoforms of PYCR, PYCR3 remains the least studied due in part to the lack of a robust recombinant expression. Herein, we describe a procedure for the expression of soluble SUMO-PYCR3 in Escherichia coli, purification of the fusion protein, and removal of the SUMO tag. PYCR3 is active with either NADPH or NADH as the coenzyme. Bi-substrate kinetic measurements obtained by varying the concentrations of both L-P5C and NADH, along with product inhibition data for l-proline, suggest a random ordered bi bi mechanism. A panel of 19 proline analogs was screened for inhibition, and the kinetics of competitive inhibition (with L-P5C) were measured for five of the compounds screened, including N-formyl-l-proline, a validated inhibitor of PYCR1. N-formyl-l-proline was found to be ten times more selective for PYCR1 over PYCR3. The SUMO-PYCR3 expression system should be useful for testing the isoform specificity of PYCR1 inhibitors.

Disease variants of human Δ1-pyrroline-5-carboxylate reductase 2 (PYCR2).

Pyrroline-5-carboxylate reductase (PYCR in humans) catalyzes the final step of l-proline biosynthesis by catalyzing the reduction of L-Δ1-pyrroline-5-carboxylate (L-P5C) to l-proline using NAD(P)H as the hydride donor. In humans, three isoforms PYCR1, PYCR2, and PYCR3 are known. Recent genome-wide association and clinical studies have revealed that homozygous mutations in human PYCR2 lead to postnatal microcephaly and hypomyelination, including hypomyelinating leukodystrophy type 10. To uncover biochemical and structural insights into human PYCR2, we characterized the steady-state kinetics of the wild-type enzyme along with two protein variants, Arg119Cys and Arg251Cys, that were previously identified in patients with microcephaly and hypomyelination. Kinetic measurements with PYCR2 suggest a sequential binding mechanism with L-P5C binding before NAD(P)H and NAD(P)+ releasing before L-Pro. Both disease-related variants are catalytically impaired. Depending on whether NADPH or NADH was used, the catalytic efficiency of the R119C protein variant was 40 or 366 times lower than that of the wild-type enzyme, while the catalytic efficiency of the R251C protein variant was 7 or 26 times lower than that of the wild-type enzyme. In addition, thermostability and circular dichroism measurements suggest that the R251C protein variant has a pronounced folding defect. These results are consistent with the involvement of Arg119Cys and Arg251Cys in disease pathology.

Structure of Thermococcus litoralis Δ1-pyrroline-2-carboxylate reductase in complex with NADH and L-proline.

L-Hydroxyproline (L-Hyp) is a nonstandard amino acid that is present in certain proteins, in some antibiotics and in the cell-wall components of plants. L-Hyp is the product of the post-translational modification of protein prolines by prolyl hydroxylase enzymes, and the isomers trans-3-hydroxy-L-proline (T3LHyp) and trans-4-hydroxy-L-proline (T4LHyp) are major components of mammalian collagen. T4LHyp follows two distinct degradation pathways in bacteria and mammals, while T3LHyp is metabolized by a two-step metabolic pathway that is conserved in bacteria and mammals, which involves a T3LHyp dehydratase and a Δ1-pyrroline-2-carboxylate (Pyr2C) reductase. In order to shed light on the structure and catalysis of the enzyme involved in the second step of the T3LHyp degradation pathway, the crystal structure of Pyr2C reductase from the archaeon Thermococcus litoralis DSM 5473 complexed with NADH and L-proline is presented. The model allows the mapping of the residues involved in cofactor and product binding and represents a valid model for rationalizing the catalysis of Pyr2C reductases.

Resolving the cofactor-binding site in the proline biosynthetic enzyme human pyrroline-5-carboxylate reductase 1.

Pyrroline-5-carboxylate reductase (PYCR) is the final enzyme in proline biosynthesis, catalyzing the NAD(P)H-dependent reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline. Mutations in the PYCR1 gene alter mitochondrial function and cause the connective tissue disorder cutis laxa. Furthermore, PYCR1 is overexpressed in multiple cancers, and the PYCR1 knock-out suppresses tumorigenic growth, suggesting that PYCR1 is a potential cancer target. However, inhibitor development has been stymied by limited mechanistic details for the enzyme, particularly in light of a previous crystallographic study that placed the cofactor-binding site in the C-terminal domain rather than the anticipated Rossmann fold of the N-terminal domain. To fill this gap, we report crystallographic, sedimentation-velocity, and kinetics data for human PYCR1. Structures of binary complexes of PYCR1 with NADPH or proline determined at 1.9 Å resolution provide insight into cofactor and substrate recognition. We see NADPH bound to the Rossmann fold, over 25 Å from the previously proposed site. The 1.85 Å resolution structure of a ternary complex containing NADPH and a P5C/proline analog provides a model of the Michaelis complex formed during hydride transfer. Sedimentation velocity shows that PYCR1 forms a concentration-dependent decamer in solution, consistent with the pentamer-of-dimers assembly seen crystallographically. Kinetic and mutational analysis confirmed several features seen in the crystal structure, including the importance of a hydrogen bond between Thr-238 and the substrate as well as limited cofactor discrimination.

Effect of R119G Mutation on Human P5CR1 Dynamic Property and Enzymatic Activity.

Pyrroline-5-carboxylate reductase (P5CR1) is a universal housekeeping enzyme that catalyzes the reduction of Δ1-pyrroline-5-carboxylate (P5C) to proline with concomitant oxidation of NAD(P)H to NAD(P)+. The enzymatic cycle between P5C and proline is important for function in amino acid metabolism, apoptosis, and intracellular redox potential balance in mitochondria. Autosomal recessive cutis laxa (ARCL) results from a mutation in P5CR1 encoded by PYCR1. Specifically, the R119G mutation is reported to be linked to ARCL although it has not yet been characterized. We synthesized R119G P5CR1 and compared it to WT P5CR1. Foldx prediction of WT and R119G mutant P5CR1 protein stability suggests that the R119G mutation could significantly reduce protein stability. We also performed enzymatic activity assays to determine how the mutation impacts P5CR1 enzymatic function. The results of these experiments show that mutagenesis of R119 to G decreases P5CR1 catalytic efficiency for 3,4-dehydro-L-proline relative to WT. Mutagenesis and kinetic studies reveal that the activity of the mutant decreases as temperature increases from 5°C to 37°C, with almost no activity at 37°C, indicating that this mutation impairs P5CR1 function in vivo. Conversely, WT P5CR1 retains its activity after incubation at 37°C and has essentially no remaining activity at 75°C. Taken together, our experimental results indicate the R119G mutation could be an involving pathomechanism for ARCL.

Effect of the R119G mutation on human P5CR structure and its interactions with NAD: Insights derived from molecular dynamics simulation and free energy analysis.

Pyrroline-5-carboxylate reductase (P5CR), an enzyme with conserved housekeeping roles, is involved in the etiology of cutis laxa. While previous work has shown that the R119G point mutation in the P5CR protein is involved, the structural mechanism behind the pathology remains to be elucidated. In order to probe the role of the R119G mutation in cutis laxa, we performed molecular dynamics (MD) simulations, essential dynamics (ED) analysis, and Molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) binding free energy calculations on wild type (WT) and mutant P5CR-NAD complex. These MD simulations and ED analyses suggest that the R119G mutation decreases the flexibility of P5CR, specifically in the substrate binding pocket, which could decrease the kinetics of the cofactor entrance and egress. Furthermore, the MM-PBSA calculations suggest the R119G mutant has a lower cofactor binding affinity for NAD than WT. Our study provides insight into the possible role of the R119G mutation during interactions between P5CR and NAD, thus bettering our understanding of how the mutation promotes cutis laxa.

In silico screening, molecular docking, and molecular dynamics studies of SNP-derived human P5CR mutants.

Pyrroline-5-carboxylate reductase (P5CR) encoded by PYCR1 gene is a housekeeping enzyme that catalyzes the reduction of P5C to proline using NAD(P)H as the cofactor. In this study, we used in silico approaches to examine the role of nonsynonymous single-nucleotide polymorphisms in the PYCR1 gene and their putative functions in the pathogenesis of Cutis Laxa. Among the 348 identified SNPs, 15 were predicted to be potentially damaging by both SIFT and PolyPhen tools; of them two SNP-derived mutations, R119G and G206W, have been previously reported to correlate with Cutis Laxa. These two mutations were therefore selected to be mapped to the wild-type (WT) P5CR structure for further structural and functional analyses. The results of comparative computational analyses using I-Mutant and Autodock reveal reductions in both stability and cofactor binding affinity of these two mutants. Comparative molecular dynamics (MD) simulations were performed to evaluate the changes in dynamic properties of P5CR upon mutations. The results reveal that the two mutations enhance the rigidity of P5CR structure, especially that of cofactor binding site, which could result in decreased kinetics of cofactor entrance and egress. Comparison between the structural properties of the WT and mutants during MD simulations shows that the enhanced rigidity of mutants results most likely from the increased number of inter-atomic interactions and the decreased number of dynamic hydrogen bonds. Our study provides novel insight into the deleterious effects of the R119G and G206W mutations on P5CR, and sheds light on the mechanisms by which these mutations mediate Cutis Laxa.

Kinetic and structural characterization of tunnel-perturbing mutants in Bradyrhizobium japonicum proline utilization A.

Proline utilization A from Bradyrhizobium japonicum (BjPutA) is a bifunctional flavoenzyme that catalyzes the oxidation of proline to glutamate using fused proline dehydrogenase (PRODH) and Δ(1)-pyrroline-5-carboxylate dehydrogenase (P5CDH) domains. Recent crystal structures and kinetic data suggest an intramolecular channel connects the two active sites, promoting substrate channeling of the intermediate Δ(1)-pyrroline-5-carboxylate/glutamate-γ-semialdehyde (P5C/GSA). In this work, the structure of the channel was explored by inserting large side chain residues at four positions along the channel in BjPutA. Kinetic analysis of the different mutants revealed replacement of D779 with Tyr (D779Y) or Trp (D779W) significantly decreased the overall rate of the PRODH-P5CDH channeling reaction. X-ray crystal structures of D779Y and D779W revealed that the large side chains caused a constriction in the central section of the tunnel, thus likely impeding the travel of P5C/GSA in the channel. The D779Y and D779W mutants have PRODH activity similar to that of wild-type BjPutA but exhibit significantly lower P5CDH activity, suggesting that exogenous P5C/GSA enters the channel upstream of Asp779. Replacement of nearby Asp778 with Tyr (D778Y) did not impact BjPutA channeling activity. Consistent with the kinetic results, the X-ray crystal structure of D778Y shows that the main channel pathway is not impacted; however, an off-cavity pathway is closed off from the channel. These findings provide evidence that the off-cavity pathway is not essential for substrate channeling in BjPutA.

Δ1-Pyrroline-5-carboxylate reductase from Arabidopsis thaliana: stimulation or inhibition by chloride ions and feedback regulation by proline depend on whether NADPH or NADH acts as co-substrate.

Δ(1)-pyrroline-5-carboxylate (P5C) reductase (P5CR) catalyses the final step of proline synthesis in plants. In Arabidopsis thaliana, protein levels are correlated neither to the corresponding mRNA copy numbers, nor to intracellular proline concentrations. The occurrence of post-translational regulatory mechanisms has therefore been hypothesized, but never assessed. The purification of A. thaliana P5CR was achieved through either a six-step protocol from cultured cells, or heterologous expression of AtP5CR in Escherichia coli. The protein was characterized with respect to structural, kinetic, and biochemical properties. P5CR was able to use either NADPH or NADH as the electron donor, with contrasting affinities and maximum reaction rates. The presence of equimolar concentrations of NADP(+) completely suppressed the NADH-dependent activity, whereas the NADPH-dependent reaction was mildly affected. Proline inhibited only the NADH-dependent reaction. At physiological values, increasing concentrations of salt progressively inhibited the NADH-dependent activity, but were stimulatory of the NADPH-dependent reaction. The biochemical properties of A. thaliana P5CR suggest a complex regulation of enzyme activity by the redox status of the pyridine nucleotide pools, and the concentrations of proline and chloride in the cytosol. Data support a to date underestimated role of P5CR in controlling stress-induced proline accumulation.

Genotype-phenotype spectrum of PYCR1-related autosomal recessive cutis laxa.

Autosomal recessive cutis laxa type 2B (ARCL2B; OMIM # 612940) is a segmental progeroid disorder caused by mutations in PYCR1 encoding pyrroline-5-carboxylate reductase 1, which is part of the conserved proline de novo synthesis pathway. Here we describe 33 patients with PYCR1-related ARCL from 27 families with initial diagnoses varying between wrinkly skin syndrome, gerodermia osteodysplastica, De Barsy syndrome or more severe progeria syndromes. Given the difficult differential diagnosis of ARCL syndromes we performed a systematic comparison of clinical features of PYCR1-related ARCL. Intrauterine growth retardation, a characteristic triangular facial gestalt, psychomotor retardation, and hypotonia were the most relevant distinctive hallmarks of ARCL due to proline de novo synthesis defects. Corneal clouding or cataracts, athetoid movements, and finger contractures were rather rare features, but had a high predictive value. In our cohort we identified 20 different PYCR1 mutations of which seven were novel. Most of the mutations accumulated in exons 4 to 6. Missense alterations of highly conserved residues were most frequent followed by splice site changes and a single nonsense mutation. Analysis of genotype-phenotype correlation revealed that patients with mutations in the first two exons had lower average clinical scores and absent or only mild intellectual disability. Structural analyses predicted interference with PYCR1 multimerization for a subset of missense mutations. These findings have implications for the clinics as well as the pathomechanism of PYCR1-related ARCL.

Coenzyme preference of Streptococcus pyogenes δ1-pyrroline-5-carboxylate reductase: evidence supporting NADPH as the physiological electron donor.

The streptococcal enzyme that catalyzes the last step in proline biosynthesis was heterologously expressed and the recombinant protein was purified to electrophoretic homogeneity and characterized thoroughly. As for δ1-pyrroline-5-carboxylate reductases from other sources, it was able to use either NADH or NADPH as the electron donor in vitro. However, with NADH the activity was markedly inhibited by physiological levels of NADP+. Results also strengthen the possibility that an unusual ordered substrate binding occurs, in which the dinucleotide binds last.

Δ1-pyrroline-5-carboxylate reductase as a new target for therapeutics: inhibition of the enzyme from Streptococcus pyogenes and effects in vivo.

Compounds able to interfere with amino acid biosynthesis have the potential to inhibit cell growth. In both prokaryotic and eukaryotic microorganisms, unless an ornithine cyclodeaminase is present, the activity of δ1-pyrroline-5-carboxylate (P5C) reductase is mandatory to proline production, and the enzyme inhibition should result in amino acid starvation, blocking in turn protein synthesis. The ability of some substituted derivatives of aminomethylenebisphosphonic acid and its analogues to interfere with the activity of the enzyme from the human pathogen Streptococcus pyogenes was investigated. Several compounds were able to suppress activity in the micromolar range of concentrations, with a mechanism of uncompetitive type with respect to the substrate P5C and non-competitive with respect to the electron donor NAD(P)H. The actual occurrence of enzyme inhibition in vivo was supported by the effects of the most active derivatives upon bacterial growth and free amino acid content.

Mutation in pyrroline-5-carboxylate reductase 1 gene in families with cutis laxa type 2.

Autosomal-recessive cutis laxa type 2 (ARCL2) is a multisystem disorder characterized by the appearance of premature aging, wrinkled and lax skin, joint laxity, and a general developmental delay. Cutis laxa includes a family of clinically overlapping conditions with confusing nomenclature, generally requiring molecular analyses for definitive diagnosis. Six genes are currently known to mutate to yield one of these related conditions. We ascertained a cohort of typical ARCL2 patients from a subpopulation isolate within eastern Canada. Homozygosity mapping with high-density SNP genotyping excluded all six known genes, and instead identified a single homozygous region near the telomere of chromosome 17, shared identically by state by all genotyped affected individuals from the families. A putative pathogenic variant was identified by direct DNA sequencing of genes within the region. The single nucleotide change leads to a missense mutation adjacent to a splice junction in the gene encoding pyrroline-5-carboxylate reductase 1 (PYCR1). Bioinformatic analysis predicted a pathogenic effect of the variant on splice donor site function. Skipping of the associated exon was confirmed in RNA from blood lymphocytes of affected homozygotes and heterozygous mutation carriers. Exon skipping leads to deletion of the reductase functional domain-coding region and an obligatory downstream frameshift. PYCR1 plays a critical role in proline biosynthesis. Pathogenicity of the genetic variant in PYCR1 is likely, given that a similar clinical phenotype has been documented for mutation carriers of another proline biosynthetic enzyme, pyrroline-5-carboxylate synthase. Our results support a significant role for proline in normal development.

Purification, characterization and crystallization of pyrroline-5-carboxylate reductase from the hyperthermophilic archeon Sulfolobus Solfataricus.

The gene SSO0495 (proC), which encodes pyrroline-5-carboxylate reductase (P5CR) from the thermoacidophilic archeon Sulfolobus solfataricus P2 (Ss-P5CR), was cloned and expressed. The purified recombinant enzyme catalyzes the thioproline dehydrogenase with concomitant oxidation of NAD(P)H to NAD(P)+. This archeal enzyme has an optimal alkaline pH in this reversible reaction and is thermostable with a half-life of approximately 30 min at 80 degrees C. At pH 9.0, the reverse activation rate is nearly 3-fold higher than at pH 7.0. The homopolymer was characterized by cross-linking and size exclusion gel filtration chromatography. Ss-P5CR was crystallized by the hanging-drop vapor-diffusion method at 37 degrees C. Diffraction data were obtained to a resolution of 3.5A and were suitable for X-ray structure determination.

Crystal structure of human pyrroline-5-carboxylate reductase.

Pyrroline-5-carboxylate reductase (P5CR) is a universal housekeeping enzyme that catalyzes the reduction of Delta(1)-pyrroline-5-carboxylate (P5C) to proline using NAD(P)H as the cofactor. The enzymatic cycle between P5C and proline is very important for the regulation of amino acid metabolism, intracellular redox potential, and apoptosis. Here, we present the 2.8 Angstroms resolution structure of the P5CR apo enzyme, its 3.1 Angstroms resolution ternary complex with NAD(P)H and substrate-analog. The refined structures demonstrate a decameric architecture with five homodimer subunits and ten catalytic sites arranged around a peripheral circular groove. Mutagenesis and kinetic studies reveal the pivotal roles of the dinucleotide-binding Rossmann motif and residue Glu221 in the human enzyme. Human P5CR is thermostable and the crystals were grown at 37 degrees C. The enzyme is implicated in oxidation of the anti-tumor drug thioproline.

Purification, characterization, and crystallization of human pyrroline-5-carboxylate reductase.

Pyrroline-5-carboxylate reductase (P5CR) catalyzes the reduction of Delta1-pyrroline-5-carboxylate (P5C) to proline with concomitant oxidation of NAD(P)H to NAD(P)(+). The enzymatic cycle between P5C and proline is very important in many physiological and pathological processes. Human P5CR was over-expressed in Escherichia coli and purified to homogeneity by chromatography. Enzymatic assays of the wild-type protein were carried out using 3,4-dehydro-L-proline as substrate and NAD(+) as cofactor. The homopolymer was characterized by cross-linking and size exclusion gel filtration chromatography. Human P5CR was crystallized by the hanging-drop vapor-diffusion method at 37 degrees C. Diffraction data were obtained to a resolution of 2.8A and were suitable for high resolution X-ray structure determination.

Purification and characterization of a functionally active Mycobacterium tuberculosis pyrroline-5-carboxylate reductase.

Pyrroline-5-carboxylate reductase (P5CR) plays an important role in the survival of Mycobacterium tuberculosis and is related to virulence of this pathogen. RT-PCR analysis indicated that proC, encoding P5CR, was expressed at the transcriptional level cultured in vitro. The His-rMtP5CR with an N-terminal His-tag (His-rMtP5CR) was firstly purified in Escherichia coli and rMtP5CR was obtained by removal of the N-terminal fusion partner using enterokinase. His-rMtP5CR had considerable beta-pleated sheet analyzed by circular dichroism spectroscopy. The effect of pH, temperature, cations, denaturants, and detergents on the purified enzyme activity and stability was characterized. The N-terminal fusion partner was found to have very little effect on the biochemical properties of P5CR.

Crystal structures of delta1-pyrroline-5-carboxylate reductase from human pathogens Neisseria meningitides and Streptococcus pyogenes.

L-proline is an amino acid that plays an important role in proteins uniquely contributing to protein folding, structure, and stability, and this amino acid serves as a sequence-recognition motif. Proline biosynthesis can occur via two pathways, one from glutamate and the other from arginine. In both pathways, the last step of biosynthesis, the conversion of delta1-pyrroline-5-carboxylate (P5C) to L-proline, is catalyzed by delta1-pyrroline-5-carboxylate reductase (P5CR) using NAD(P)H as a cofactor. We have determined the first crystal structure of P5CR from two human pathogens, Neisseria meningitides and Streptococcus pyogenes, at 2.0 angstroms and 2.15 angstroms resolution, respectively. The catalytic unit of P5CR is a dimer composed of two domains, but the biological unit seems to be species-specific. The N-terminal domain of P5CR is an alpha/beta/alpha sandwich, a Rossmann fold. The C-terminal dimerization domain is rich in alpha-helices and shows domain swapping. Comparison of the native structure of P5CR to structures complexed with L-proline and NADP+ in two quite different primary sequence backgrounds provides unique information about key functional features: the active site and the catalytic mechanism. The inhibitory L-proline has been observed in the crystal structure.

Crystal structures of Delta1-piperideine-2-carboxylate/Delta1-pyrroline-2-carboxylate reductase belonging to a new family of NAD(P)H-dependent oxidoreductases: conformational change, substrate recognition, and stereochemistry of the reaction.

Delta(1)-Piperideine-2-carboxylate/Delta(1)-pyrroline-2-carboxylate reductase from Pseudomonas syringae pv. tomato belongs to a novel sub-class in a large family of NAD(P)H-dependent oxidoreductases distinct from the conventional MDH/LDH superfamily characterized by the Rossmann fold. We have determined the structures of the following three forms of the enzyme: the unliganded form, the complex with NADPH, and the complex with NADPH and pyrrole-2-carboxylate at 1.55-, 1.8-, and 1.7-A resolutions, respectively. The enzyme exists as a dimer, and the subunit consists of three domains; domain I, domain II (NADPH binding domain), and domain III. The core of the NADPH binding domain consists of a seven-stranded predominantly antiparallel beta-sheet fold (which we named SESAS) that is characteristic of the new oxidoreductase family. The enzyme preference for NADPH over NADH is explained by the cofactor binding site architecture. A comparison of the overall structures revealed that the mobile domains I and III change their conformations to produce the catalytic form. This conformational change plays important roles in substrate recognition and the catalytic process. The active site structure of the catalytic form made it possible to identify the catalytic Asp:Ser:His triad and investigate the catalytic mechanism from a stereochemical point of view.