De Novo Discovery of Pseudo-Natural Prenylated Macrocyclic Peptide Ligands.

Prenylation of peptides is widely observed in the secondary metabolites of diverse organisms, granting peptides unique chemical properties distinct from proteinogenic amino acids. Discovery of prenylated peptide agents has largely relied on isolation or genome mining of naturally occurring molecules. To devise a platform technology for de novo discovery of artificial prenylated peptides targeting a protein of choice, here we have integrated the thioether-macrocyclic peptide (teMP) library construction/selection technology, so-called RaPID (Random nonstandard Peptides Integrated Discovery) system, with a Trp-C3-prenyltransferase KgpF involved in the biosynthesis of a prenylated natural product. This unique enzyme exhibited remarkably broad substrate tolerance, capable of modifying various Trp-containing teMPs to install a prenylated residue with tricyclic constrained structure. We constructed a vast library of prenylated teMPs and subjected it to in vitro selection against a phosphoglycerate mutase. This selection platform has led to the identification of a pseudo-natural prenylated teMP inhibiting the target enzyme with an IC50 of 30 nM. Importantly, the prenylation was essential for the inhibitory activity, enhanced serum stability, and cellular uptake of the peptide, highlighting the benefits of peptide prenylation. This work showcases the de novo discovery platform for pseudo-natural prenylated peptides, which is readily applicable to other drug targets.

Genetic and functional diversity of ß-N-acetylgalactosamine-targeting glycosidases expanded by deep-sea metagenome analysis.

ß-N-Acetylgalactosamine-containing glycans play essential roles in several biological processes, including cell adhesion, signal transduction, and immune responses. ß-N-Acetylgalactosaminidases hydrolyze ß-N-acetylgalactosamine linkages of various glycoconjugates. However, their biological significance remains ambiguous, primarily because only one type of enzyme, exo-ß-N-acetylgalactosaminidases that specifically act on ß-N-acetylgalactosamine residues, has been documented to date. In this study, we identify four groups distributed among all three domains of life and characterize eight ß-N-acetylgalactosaminidases and ß-N-acetylhexosaminidase through sequence-based screening of deep-sea metagenomes and subsequent searching of public protein databases. Despite low sequence similarity, the crystal structures of these enzymes demonstrate that all enzymes share a prototype structure and have diversified their substrate specificities (oligosaccharide-releasing, oligosaccharide/monosaccharide-releasing, and monosaccharide-releasing) through the accumulation of mutations and insertional amino acid sequences. The diverse ß-N-acetylgalactosaminidases reported in this study could facilitate the comprehension of their structures and functions and present evolutionary pathways for expanding their substrate specificity.

Structure of the human Bre1 complex bound to the nucleosome.

Histone H2B monoubiquitination (at Lys120 in humans) regulates transcription elongation and DNA repair. In humans, H2B monoubiquitination is catalyzed by the heterodimeric Bre1 complex composed of Bre1A/RNF20 and Bre1B/RNF40. The Bre1 proteins generally function as tumor suppressors, while in certain cancers, they facilitate cancer cell proliferation. To obtain structural insights of H2BK120 ubiquitination and its regulation, we report the cryo-electron microscopy structure of the human Bre1 complex bound to the nucleosome. The two RING domains of Bre1A and Bre1B recognize the acidic patch and the nucleosomal DNA phosphates around SHL 6.0-6.5, which are ideally located to recruit the E2 enzyme and ubiquitin for H2BK120-specific ubiquitination. Mutational experiments suggest that the two RING domains bind in two orientations and that ubiquitination occurs when Bre1A binds to the acidic patch. Our results provide insights into the H2BK120-specific ubiquitination by the Bre1 proteins and suggest that H2B monoubiquitination can be regulated by nuclesomal DNA flexibility.

A Compact Reprogrammed Genetic Code for De Novo Discovery of Proteolytically Stable Thiopeptides.

Thiopeptides make up a group of structurally complex peptidic natural products holding promise in bioengineering applications. The previously established thiopeptide/mRNA display platform enables de novo discovery of natural product-like thiopeptides with designed bioactivities. However, in contrast to natural thiopeptides, the discovered structures are composed predominantly of proteinogenic amino acids, which results in low metabolic stability in many cases. Here, we redevelop the platform and demonstrate that the utilization of compact reprogrammed genetic codes in mRNA display libraries can lead to the discovery of thiopeptides predominantly composed of nonproteinogenic structural elements. We demonstrate the feasibility of our designs by conducting affinity selections against Traf2- and NCK-interacting kinase (TNIK). The experiment identified a series of thiopeptides with high affinity to the target protein (the best KD = 2.1 nM) and kinase inhibitory activity (the best IC50 = 0.15 µM). The discovered compounds, which bore as many as 15 nonproteinogenic amino acids in an 18-residue macrocycle, demonstrated high metabolic stability in human serum with a half-life of up to 99 h. An X-ray cocrystal structure of TNIK in complex with a discovered thiopeptide revealed how nonproteinogenic building blocks facilitate the target engagement and orchestrate the folding of the thiopeptide into a noncanonical conformation. Altogether, the established platform takes a step toward the discovery of thiopeptides with high metabolic stability for early drug discovery applications.

Switching Prenyl Donor Specificities of Cyanobactin Prenyltransferases.

Prenyltransferases in cyanobactin biosynthesis are of growing interest as peptide alkylation biocatalysts, but their prenylation modes characterized so far have been limited to dimethylallylation (C5) or geranylation (C10). Here we engaged in structure-guided engineering of the prenyl-binding pocket of a His-C2-geranyltransferase LimF to modulate its prenylation mode. Contraction of the pocket by a single mutation led to a His-C2-dimethylallyltransferase. More importantly, pocket expansion by a double mutation successfully repurposed LimF for farnesylation (C15), which is an unprecedented mode in this family. Furthermore, the obtained knowledge of the essential residues to construct the farnesyl-binding pocket has allowed for rational design of a Tyr-O-farnesyltransferase by a triple mutation of a Tyr-O-dimethylallyltransferase PagF. These results provide an approach to manipulate the prenyl specificity of cyanobactin prenyltransferases, broadening the chemical space covered by this class of enzymes and expanding the toolbox of peptide alkylation biocatalysts.

Biallelic structural variations within FGF12 detected by long-read sequencing in epilepsy.

We discovered biallelic intragenic structural variations (SVs) in FGF12 by applying long-read whole genome sequencing to an exome-negative patient with developmental and epileptic encephalopathy (DEE). We also found another DEE patient carrying a biallelic (homozygous) single-nucleotide variant (SNV) in FGF12 that was detected by exome sequencing. FGF12 heterozygous recurrent missense variants with gain-of-function or heterozygous entire duplication of FGF12 are known causes of epilepsy, but biallelic SNVs/SVs have never been described. FGF12 encodes intracellular proteins interacting with the C-terminal domain of the alpha subunit of voltage-gated sodium channels 1.2, 1.5, and 1.6, promoting excitability by delaying fast inactivation of the channels. To validate the molecular pathomechanisms of these biallelic FGF12 SVs/SNV, highly sensitive gene expression analyses using lymphoblastoid cells from the patient with biallelic SVs, structural considerations, and Drosophila in vivo functional analysis of the SNV were performed, confirming loss-of-function. Our study highlights the importance of small SVs in Mendelian disorders, which may be overlooked by exome sequencing but can be detected efficiently by long-read whole genome sequencing, providing new insights into the pathomechanisms of human diseases.

Brain monoamine vesicular transport disease caused by homozygous SLC18A2 variants: A study in 42 affected individuals.

PURPOSE: Brain monoamine vesicular transport disease is an infantile-onset movement disorder that mimics cerebral palsy. In 2013, the homozygous SLC18A2 variant, p.Pro387Leu, was first reported as a cause of this rare disorder, and dopamine agonists were efficient for treating affected individuals from a single large family. To date, only 6 variants have been reported. In this study, we evaluated genotype-phenotype correlations in individuals with biallelic SLC18A2 variants. METHODS: A total of 42 affected individuals with homozygous SLC18A2 variant alleles were identified. We evaluated genotype-phenotype correlations and the missense variants in the affected individuals based on the structural modeling of rat VMAT2 encoded by Slc18a2, with cytoplasm- and lumen-facing conformations. A Caenorhabditis elegans model was created for functional studies. RESULTS: A total of 19 homozygous SLC18A2 variants, including 3 recurrent variants, were identified using exome sequencing. The affected individuals typically showed global developmental delay, hypotonia, dystonia, oculogyric crisis, and autonomic nervous system involvement (temperature dysregulation/sweating, hypersalivation, and gastrointestinal dysmotility). Among the 58 affected individuals described to date, 16 (28%) died before the age of 13 years. Of the 17 patients with p.Pro237His, 9 died, whereas all 14 patients with p.Pro387Leu survived. Although a dopamine agonist mildly improved the disease symptoms in 18 of 21 patients (86%), some affected individuals with p.Ile43Phe and p.Pro387Leu showed milder phenotypes and presented prolonged survival even without treatment. The C. elegans model showed behavioral abnormalities. CONCLUSION: These data expand the phenotypic and genotypic spectra of SLC18A2-related disorders.

Structural basis of transcription regulation by CNC family transcription factor, Nrf2.

Several basic leucine zipper (bZIP) transcription factors have accessory motifs in their DNA-binding domains, such as the CNC motif of CNC family or the EHR motif of small Maf (sMaf) proteins. CNC family proteins heterodimerize with sMaf proteins to recognize CNC-sMaf binding DNA elements (CsMBEs) in competition with sMaf homodimers, but the functional role of the CNC motif remains elusive. In this study, we report the crystal structures of Nrf2/NFE2L2, a CNC family protein regulating anti-stress transcriptional responses, in a complex with MafG and CsMBE. The CNC motif restricts the conformations of crucial Arg residues in the basic region, which form extensive contact with the DNA backbone phosphates. Accordingly, the Nrf2-MafG heterodimer has approximately a 200-fold stronger affinity for CsMBE than canonical bZIP proteins, such as AP-1 proteins. The high DNA affinity of the CNC-sMaf heterodimer may allow it to compete with the sMaf homodimer on target genes without being perturbed by other low-affinity bZIP proteins with similar sequence specificity.

De Novo Discovery of Thiopeptide Pseudo-natural Products Acting as Potent and Selective TNIK Kinase Inhibitors.

Bioengineering of ribosomally synthesized and post-translationally modified peptides (RiPPs) is an emerging approach to explore the diversity of pseudo-natural product structures for drug discovery purposes. However, despite the initial advances in this area, bioactivity reprogramming of multienzyme RiPP biosynthetic pathways remains a major challenge. Here, we report a platform for de novo discovery of functional thiopeptides based on reengineered biosynthesis of lactazole A, a RiPP natural product assembled by five biosynthetic enzymes. The platform combines in vitro biosynthesis of lactazole-like thiopeptides and mRNA display to prepare and screen large (≥1012) combinatorial libraries of pseudo-natural products. We demonstrate the utility of the developed protocols in an affinity selection against Traf2- and NCK-interacting kinase (TNIK), a protein involved in several cancers, which yielded a plethora of candidate thiopeptides. Of the 11 synthesized compounds, 9 had high affinities for the target kinase (best KD = 1.2 nM) and 10 inhibited its enzymatic activity (best Ki = 3 nM). X-ray structural analysis of the TNIK/thiopeptide interaction revealed the unique mode of substrate-competitive inhibition exhibited by two of the discovered compounds. The thiopeptides internalized to the cytosol of HEK293H cells as efficiently as the known cell-penetrating peptide Tat (4-6 µM). Accordingly, the most potent compound, TP15, inhibited TNIK in HCT116 cells. Altogether, our platform enables the exploration of pseudo-natural thiopeptides with favorable pharmacological properties in drug discovery applications.

Structural basis of the regulation of the normal and oncogenic methylation of nucleosomal histone H3 Lys36 by NSD2.

Dimethylated histone H3 Lys36 (H3K36me2) regulates gene expression, and aberrant H3K36me2 upregulation, resulting from either the overexpression or point mutation of the dimethyltransferase NSD2, is found in various cancers. Here we report the cryo-electron microscopy structure of NSD2 bound to the nucleosome. Nucleosomal DNA is partially unwrapped, facilitating NSD2 access to H3K36. NSD2 interacts with DNA and H2A along with H3. The NSD2 autoinhibitory loop changes its conformation upon nucleosome binding to accommodate H3 in its substrate-binding cleft. Kinetic analysis revealed that two oncogenic mutations, E1099K and T1150A, increase NSD2 catalytic turnover. Molecular dynamics simulations suggested that in both mutants, the autoinhibitory loop adopts an open state that can accommodate H3 more often than the wild-type. We propose that E1099K and T1150A destabilize the interactions that keep the autoinhibitory loop closed, thereby enhancing catalytic turnover. Our analyses guide the development of specific inhibitors of NSD2.

An Ultrapotent and Selective Cyclic Peptide Inhibitor of Human ß-Factor XIIa in a Cyclotide Scaffold.

Cyclotides are plant-derived peptides with complex structures shaped by their head-to-tail cyclic backbone and cystine knot core. These structural features underpin the native bioactivities of cyclotides, as well as their beneficial properties as pharmaceutical leads, including high proteolytic stability and cell permeability. However, their inherent structural complexity presents a challenge for cyclotide engineering, particularly for accessing libraries of sufficient chemical diversity to design potent and selective cyclotide variants. Here, we report a strategy using mRNA display enabling us to select potent cyclotide-based FXIIa inhibitors from a library comprising more than 1012 members based on the cyclotide scaffold of Momordica cochinchinensis trypsin inhibitor-II (MCoTI-II). The most potent and selective inhibitor, cMCoFx1, has a pM inhibitory constant toward FXIIa with greater than three orders of magnitude selectivity over related serine proteases, realizing specific inhibition of the intrinsic coagulation pathway. The cocrystal structure of cMCoFx1 and FXIIa revealed interactions at several positions across the contact interface that conveyed high affinity binding, highlighting that such cyclotides are attractive cystine knot scaffolds for therapeutic development.

Pathogenic variants in the survival of motor neurons complex gene GEMIN5 cause cerebellar atrophy.

Cerebellar ataxia is a genetically heterogeneous disorder. GEMIN5 encoding an RNA-binding protein of the survival of motor neuron complex, is essential for small nuclear ribonucleoprotein biogenesis, and it was recently reported that biallelic loss-of-function variants cause neurodevelopmental delay, hypotonia, and cerebellar ataxia. Here, whole-exome analysis revealed compound heterozygous GEMIN5 variants in two individuals from our cohort of 162 patients with cerebellar atrophy/hypoplasia. Three novel truncating variants and one previously reported missense variant were identified: c.2196dupA, p.(Arg733Thrfs*6) and c.1831G > A, p.(Val611Met) in individual 1, and c.3913delG, p.(Ala1305Leufs*14) and c.4496dupA, p.(Tyr1499*) in individual 2. Western blotting analysis using lymphoblastoid cell lines derived from both affected individuals showed significantly reduced levels of GEMIN5 protein. Zebrafish model for null variants p.(Arg733Thrfs*6) and p.(Ala1305Leufs*14) exhibited complete lethality at 2 weeks and recapitulated a distinct dysplastic phenotype. The phenotypes of affected individuals and the zebrafish mutant models strongly suggest that biallelic loss-of-function variants in GEMIN5 cause cerebellar atrophy/hypoplasia.

De novo ARF3 variants cause neurodevelopmental disorder with brain abnormality.

An optimal Golgi transport system is important for mammalian cells. The adenosine diphosphate (ADP) ribosylation factors (ARF) are key proteins for regulating cargo sorting at the Golgi network. In this family, ARF3 mainly works at the trans-Golgi network (TGN), and no ARF3-related phenotypes have yet been described in humans. We here report the clinical and genetic evaluations of two unrelated children with de novo pathogenic variants in the ARF3 gene: c.200A > T (p.Asp67Val) and c.296G > T (p.Arg99Leu). Although the affected individuals presented commonly with developmental delay, epilepsy and brain abnormalities, there were differences in severity, clinical course and brain lesions. In vitro subcellular localization assays revealed that the p.Arg99Leu mutant localized to Golgi apparatus, similar to the wild-type, whereas the p.Asp67Val mutant tended to show a disperse cytosolic pattern together with abnormally dispersed Golgi localization, similar to that observed in a known dominant negative variant (p.Thr31Asn). Pull-down assays revealed that the p.Asp67Val had a loss-of-function effect and the p.Arg99Leu variant had increased binding of the adaptor protein, Golgi-localized, γ-adaptin ear-containing, ARF-binding protein 1 (GGA1), supporting the gain of function. Furthermore, in vivo studies revealed that p.Asp67Val transfection led to lethality in flies. In contrast, flies expressing p.Arg99Leu had abnormal rough eye, as observed in the gain-of-function variant p.Gln71Leu. These data indicate that two ARF3 variants, the possibly loss-of-function p.Asp67Val and the gain-of-function p.Arg99Leu, both impair the Golgi transport system. Therefore, it may not be unreasonable that they showed different clinical features like diffuse brain atrophy (p.Asp67Val) and cerebellar hypoplasia (p.Arg99Leu).

Systematic analysis of exonic germline and postzygotic de novo mutations in bipolar disorder.

Bipolar disorder is a severe mental illness characterized by recurrent manic and depressive episodes. To better understand its genetic architecture, we analyze ultra-rare de novo mutations in 354 trios with bipolar disorder. For germline de novo mutations, we find significant enrichment of loss-of-function mutations in constrained genes (corrected-P = 0.0410) and deleterious mutations in presynaptic active zone genes (FDR = 0.0415). An analysis integrating single-cell RNA-sequencing data identifies a subset of excitatory neurons preferentially expressing the genes hit by deleterious mutations, which are also characterized by high expression of developmental disorder genes. In the analysis of postzygotic mutations, we observe significant enrichment of deleterious ones in developmental disorder genes (P = 0.00135), including the SRCAP gene mutated in two unrelated probands. These data collectively indicate the contributions of both germline and postzygotic mutations to the risk of bipolar disorder, supporting the hypothesis that postzygotic mutations of developmental disorder genes may contribute to bipolar disorder.

OTUD5 Variants Associated With X-Linked Intellectual Disability and Congenital Malformation.

BACKGROUND: X-linked intellectual disability (XLID), which occurs predominantly in males, is a relatively common and genetically heterogeneous disorder in which over 100 mutated genes have been reported. The OTUD5 gene at Xp11.23 encodes ovarian tumor deubiquitinase 5 protein, which is a deubiquitinating enzyme member of the ovarian tumor family. LINKage-specific-deubiquitylation-deficiency-induced embryonic defects (LINKED) syndrome, arising from pathogenic OTUD5 variants, was recently reported as a new XLID with additional congenital anomalies. METHODS: We investigated three affected males (49- and 47-year-old brothers [Individuals 1 and 2] and a 2-year-old boy [Individual 3]) from two families who showed developmental delay. Their common clinical features included developmental delay, hypotonia, short stature, and distinctive facial features, such as telecanthus and a depressed nasal bridge. Individuals 1 and 2 showed epilepsy and brain magnetic resonance imaging showed a thin corpus callosum and mild ventriculomegaly. Individual 3 showed congenital malformations, including tetralogy of Fallot, hypospadias, and bilateral cryptorchidism. To identify the genetic cause of these features, we performed whole-exome sequencing. RESULTS: A hemizygous OTUD5 missense variant, c.878A>T, p.Asn293Ile [NM_017602.4], was identified in one family with Individuals 1 and 2, and another missense variant, c.1210 C>T, p.Arg404Trp, in the other family with Individual 3, respectively. The former variant has not been registered in public databases and was predicted to be pathogenic by multiple in silico prediction tools. The latter variant p.Arg404Trp was previously reported as a pathogenic OTUD5 variant, and Individual 3 showed a typical LINKED syndrome phenotype. However, Individuals 1 and 2, with the novel variant (p.Asn293Ile), showed no cardiac or genitourinary malformations. CONCLUSIONS: Unlike previous reports of LINKED syndrome, which described early lethality with congenital cardiac anomalies, our three cases are still alive. Notably, the adult brothers with the novel missense OTUD5 variant have lived into their forties. This may be indicative of a milder phenotype as a possible genotype-phenotype correlation. These findings imply a possible long-term prognosis for individuals with this new XLID syndrome, and a wider phenotypic variation than initially thought.

Ribosomal synthesis and de novo discovery of bioactive foldamer peptides containing cyclic ß-amino acids.

Peptides that contain ß-amino acids display stable secondary structures, such as helices and sheets, and are often referred to as foldamers. Cyclic ß2,3-amino acids (cßAAs), such as 2-aminocyclohexanecarboxylic acid (2-ACHC), are strong helix/turn inducers due to their restricted conformations. Here we report the ribosomal synthesis of foldamer peptides that contain multiple, up to ten, consecutive cßAAs via genetic code reprogramming. We also report the de novo discovery of macrocyclic cßAA-containing peptides capable of binding to a protein target. As a demonstration, potent binders with low-to-subnanomolar KD values were identified for human factor XIIa (hFXIIa) and interferon-gamma receptor 1, from a library of their 1012 members. One of the anti-hFXIIa macrocyclic peptides that exhibited a high inhibitory activity and serum stability was co-crystallized with hFXIIa. The X-ray structure revealed that it adopts an antiparallel ß-sheet structure induced by a (1S,2S)-2-ACHC residue via the formation of two γ-turns. This work demonstrates the potential of this platform to explore the previously inaccessible sequence space of cßAA-containing peptides.

Development and Structural Evaluation of N-Alkylated trans-2-Phenylcyclopropylamine-Based LSD1 Inhibitors.

Lysine-specific demethylase 1 (LSD1) is a flavin adenine dinucleotide (FAD)-dependent enzyme that catalyzes the demethylation of histone H3 and regulates gene expression. Because it is implicated in the regulation of diseases such as acute myeloid leukemia, potent LSD1-specific inhibitors have been pursued. Trans-2-phenylcyclopropylamine (2-PCPA)-based inhibitors featuring substitutions on the amino group have emerged, with sub-micromolar affinities toward LSD1 and high selectivities over monoamine oxidases (MAOs). We synthesized two N-alkylated 2-PCPA-based LSD1 inhibitors, S2116 and S2157, based on the previously developed S2101. S2116 and S2157 exhibited enhanced potency for LSD1 by 2.0- to 2.6-fold, as compared with S2101. In addition, they exhibited improved selectivity over MAOs. Structural analyses of LSD1 co-crystallized with S2101, S2116, S2157, or another N-alkylated inhibitor (FCPA-MPE) confirmed that the N-substituents enhance the potency of a 2-PCPA-based inhibitor of LSD1, without constituting the adduct formed with FAD.

A novel PAK1 variant causative of neurodevelopmental disorder with postnatal macrocephaly.

p21-activated kinases (PAKs) are protein serine/threonine kinases stimulated by Rho-family p21 GTPases such as CDC42 and RAC. PAKs have been implicated in several human disorders, with pathogenic variants in PAK3 associated with intellectual disability and several PAK members, especially PAK1 and PAK4, overexpressed in human cancer. Recently, de novo PAK1 variants were reported to be causative of neurodevelopmental disorder (ND) with secondary macrocephaly in three patients. We herein report a fourth patient with ND, epilepsy, and macrocephaly caused by a de novo PAK1 missense variant. Two previously reported missense PAK1 variants functioned as activating alleles by reducing PAK1 homodimerization. To examine the pathogenicity of the identified novel p.Ser110Thr variant, we carried out in silico structural analysis. Our findings suggest that this variant also prevents PAK1 homodimerization, leading to constitutive PAK1 activation.

The recurrent postzygotic pathogenic variant p.Glu47Lys in RHOA causes a novel recognizable neuroectodermal phenotype.

RHOA is a member of the Rho family of GTPases that are involved in fundamental cellular processes including cell adhesion, migration, and proliferation. RHOA can stimulate the formation of stress fibers and focal adhesions and is a key regulator of actomyosin dynamics in various tissues. In a Genematcher-facilitated collaboration, we were able to identify four unrelated individuals with a specific phenotype characterized by hypopigmented areas of the skin, dental anomalies, body asymmetry, and limb length discrepancy due to hemihypotrophy of one half of the body, as well as brain magnetic resonance imaging (MRI) anomalies. Using whole-exome and ultra-deep amplicon sequencing and comparing genomic data of affected and unaffected areas of the skin, we discovered that all four individuals carried the identical RHOA missense variant, c.139G>A; p.Glu47Lys, in a postzygotic state. Molecular modeling and in silico analysis of the affected p.Glu47Lys residue in RHOA indicated that this exchange is predicted to specifically alter the interaction of RHOA with its downstream effectors containing a PKN-type binding domain and thereby disrupts its ability to activate signaling. Our findings indicate that the recurrent postzygotic RHOA missense variant p.Glu47Lys causes a specific mosaic disorder in humans.