Recent progress in the reprogramming of nonribosomal peptide synthetases.

Nonribosomal peptide synthetases (NRPSs) biosynthesize nonribosomal peptide (NRP) natural products, which belong to the most promising resources for drug discovery and development because of their wide range of therapeutic applications. The results of genetic, biochemical, and bioinformatics analyses have enhanced our understanding of the mechanisms of the NRPS machinery. A major goal in NRP biosynthesis is to reprogram the NRPS machinery to enable the biosynthetic production of designed peptides. Reprogramming strategies for the NRPS machinery have progressed considerably in recent years, thereby increasing the yields and generating modified peptides. Here, the recent progress in NRPS reprogramming and its application in peptide synthesis are described.

System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study.

The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.

Fragment Molecular Orbital Based Affinity Prediction toward Pyruvate Dehydrogenase Kinases: Insights into the Charge Transfer in Hydrogen Bond Networks.

The fragment molecular orbital (FMO) method is a fast quantum-mechanics method that divides systems into pieces of fragments and performs ab initio calculations. The method has been expected to improve the accuracy of describing protein-ligand interactions by incorporating electronic effects. In this article, FMO calculation with solvation methods were applied to the affinity prediction at the ATP-binding site of PDHK4. As the ionized aspartic acid lies at the center and is involved in the complex hydrogen bond networks, this system has turned out to be a difficult target to describe by traditional molecular-mechanics method. In the FMO calculation with the polarizable continuum model (PCM) solvation method, a considerable amount of charge (-0.27e) was transferred from the ionized aspartate to the surrounding residues. We found that using FMO with the PCM solvation method was important to increase the correlation, and by incorporating the ligand deformation energy, the correlation was improved to R = 0.81 for whole twelve compounds and R = 0.91 without one outlier compound.

Design, Synthesis and Structure-Activity Relationship Studies of Protein Kinase CK2 Inhibitors Containing a Purine Scaffold.

Protein kinase CK2 (CK2) is involved in the suppression of gene expression, protein synthesis, cell proliferation, and apoptosis, thus making it a target protein for the development of therapeutics toward cancer, nephritis, and coronavirus disease 2019. Using the solvent dipole ordering-based method for virtual screening, we identified and designed new candidate CK2α inhibitors containing purine scaffolds. Virtual docking experiments supported by experimental structure-activity relationship studies identified the importance of the 4-carboxyphenyl group at the 2-position, a carboxamide group at the 6-position, and an electron-rich phenyl group at the 9-position of the purine scaffold. Docking studies based on the crystal structures of CK2α and inhibitor (PDBID: 5B0X) successfully predicted the binding mode of 4-(6-carbamoyl-8-oxo-9-phenyl-8,9-dihydro-7H-purin-2-yl) benzoic acid (11), and the results were used to design stronger small molecule targets for CK2α inhibition. Interaction energy analysis suggested that 11 bound around the hinge region without the water molecule (W1) near Trp176 and Glu81 that is frequently reported in crystal structures of CK2α inhibitor complexes. X-ray crystallographic data for 11 bound to CK2α was in very good agreement with the docking experiments, and consistent with activity. From the structure-activity relationship (SAR) studies presented here, 4-(6-Carbamoyl-9-(4-(dimethylamino)phenyl)-8-oxo-8,9-dihydro-7H-purin-2-yl) benzoic acid (12) was identified as an improved active purine-based CK2α inhibitor with an IC50 of 4.3 µM. These active compounds with an unusual binding mode are expected to inspire new CK2α inhibitors and the development of therapeutics targeting CK2 inhibition.

Bivalent binding mode of an amino-pyrazole inhibitor indicates the potentials for CK2α1-selective inhibitors.

Casein kinase 2 (CK2) is a vital protein kinase that consists of two catalytic subunits (CK2α1 and/or CK2α2) and two regulatory subunits (CK2ß). CK2α1 is a drug target for nephritis and cancers, while CK2α2 is a serious off-target because its inhibition causes testicular toxicity. High similarity between the isozymes CK2α1 and CK2α2 make it difficult to design CK2α1-specific inhibitors. Herein, the crystal structures of CK2α1 and CK2α2 complexed with a 3-amino-pyrazole inhibitor revealed the remarkable differences in the protein-inhibitor interaction modes. This inhibitor bound to the ATP binding sites of both isozymes in apparently distinct orientations. In addition, another molecule of this inhibitor bound to CK2α1, but not to CK2α2, at the CK2ß protein-protein interface. Binding energy calculations and biochemical experiments suggested that this inhibitor possesses the conventional ATP-competitive characteristics with moderate allosteric function in a molecular glue mechanism. These results will assist the potential design of potent and selective CK2α1 inhibitors.

Elongation of the side chain by linear alkyl groups increases the potency of salacinol, a potent α-glucosidase inhibitor from the Ayurvedic traditional medicine "Salacia," against human intestinal maltase.

Four chain-extended analogs (12a-12d) and two related de-O-sulfonated analogs (13a and 13c) by introducing alkyl groups (a: R = C3H7, b R = C6H13, c: R = C8H17, d: R = C10H21) to the side chains of salacinol (1), a natural α-glucosidase inhibitor from Ayurvedic traditional medicine "Salacia", were synthesized. The α-glucosidase inhibitory activities of all the synthesized analogs were evaluated in vitro. Against human intestinal maltase, the inhibitory activities of 12a and 13a with seven-carbon side chain were equal to that of 1. In contrast, analogs (12b-12d, and 13c) exhibited higher level of inhibitory activity against the same enzyme than 1 and had equal or higher potency than those of the clinically used anti-diabetics, voglibose, acarbose, and miglitol. Thus, elongation of the side chains of 1 was effective for specifically increasing the inhibitory activity against human intestinal maltase.

Fragment-based lead discovery to identify novel inhibitors that target the ATP binding site of pyruvate dehydrogenase kinases.

A fragment-based lead discovery approach was applied to Pyruvate Dehydrogenase Kinases (PDHKs) to discover inhibitors against the ATP binding site with novel chemotypes. X-ray fragment screening toward PDHK4 provided a fragment hit 1 with a characteristic interaction in a deep pocket of the ATP binding site. While known inhibitors utilize several water molecules in a deep pocket to form water-mediated hydrogen bond interactions, the fragment hit binds deeper in the pocket with a hydrophobic group. Displacement of a remaining water molecule in the pocket led to the identification of lead compound 7 with a notable improvement in inhibition potency. This lead compound possessed high ligand efficiency (LE) and showed decent selectivity profile. Two additional lead compounds 10 and 13 with new scaffolds with tricyclic and bicyclic cores were generated by merging structural information of another fragment hit 2. The characteristic interaction of these novel inhibitors in a deep pocket provides new structural insights about PDHKs ATP binding site and opens a novel direction for the development of PDHKs inhibitors.

Activity, Binding, and Modeling Studies of a Reprogrammed Aryl Acid Adenylation Domain with an Enlarged Substrate Binding Pocket.

The gatekeeping adenylation (A) domain of the non-ribosomal peptide synthetase (NRPS) selectively incorporates specific proteinogenic/non-proteinogenic amino acid into a growing peptide chain. The EntE of the enterobactin NRPS is a discrete aryl acid A-domain with 2,3-dihydroxybenzoic acid (DHB) substrate specificity. Reprogrammed EntE N235G variant possesses an enlarged substrate recognition site, and is capable of accepting non-native aryl acids. Biochemical characterization of this unique substrate recognition site should provide a better understanding of activi-site microenvironments. Here, we synthesized a non-hydrolysable adenylate analogue with 2-aminobenzoic acid (2-ABA), 3-aminobenzoic acid (3-ABA), and 4-aminobenzoic acid (4-ABA) and used them to calculate the apparent inhibition constants (Kiapp.). Dose-response experiments using 3-ABA-sulfamoyladenosine (AMS) provided Kiapp. values of 596 nM for wild-type EntE and 2.4 nM for the N235G variants. These results suggest that 3-amino group of benzoic acid plays an important role in substrate recognition by the N235G variant. These findings would help designing aryl acid substrates with substituents at the 2- and 3-positions.

Precise Probing of Residue Roles by NRPS Code Swapping: Mutation, Enzymatic Characterization, Modeling, and Substrate Promiscuity of Aryl Acid Adenylation Domains.

Aryl acids are most commonly found in iron-scavenging siderophores but are not limited to them. The nonribosomal peptide synthetase (NRPS) codes of aryl acids remain poorly elucidated relative to those of amino acids. Here, we defined more precisely the role of active-site residues in aryl acid adenylation domains (A-domains) by gradually grafting the NRPS codes used for salicylic acid (Sal) into an archetypal aryl acid A-domain, EntE [specific for the substrate 2,3-dihydroxybenzoic acid (DHB)]. Enzyme kinetics and modeling studies of these EntE variants demonstrated that the NRPS code residues at positions 236, 240, and 339 collectively regulate the substrate specificity toward DHB and Sal. Furthermore, the EntE variants exhibited the ability to activate the non-native aryl acids 3-hydroxybenzoic acid, 3-aminobenzoic acid, 3-fluorobenzoic acid, and 3-chlorobenzoic acid. These studies enhance our knowledge of the NRPS codes of aryl acids and could be exploited to reprogram aryl acid A-domains for non-native aryl acids.

A new chemosensor for cyanide in blood based on the Pd complex of 2-(5-bromo-2-pyridylazo)-5-[N-n-propyl-N-(3-sulfopropyl)amino]phenol.

A new indirect chemosensor for the detection of cyanide in blood is developed. 2-(5-Bromo-2-pyridylazo)-5-[N-n-propyl-N-(3-sulfopropyl)amino]phenol, a yellow dye, forms a blue-coloured complex with palladium ions. The yellow colour of this complex is regained upon reaction with cyanide ions. The complex shows high selectivity for the detection of cyanide over 16 other anions. The system was applied to two different methods for the detection of cyanide in human whole blood. As a quantitative absorbance method, blood samples were mixed with acid, and the resulting vaporised hydrogen cyanide was absorbed in an alkaline solution containing the complex in a Conway cell. The resulting absorbance response of the solution at 450 nm is linear over the range 4-40 µM (R2 = 1.000), and the limit of detection is 0.6 µM. Furthermore, the complex-soaked paper is applicable as a test strip for cyanide detection. When a test strip is used with 0.5 mL of blood, the limit of detection is 15 µM. The detection limits of these two methods are below the toxic blood cyanide concentration (19 µM). Therefore, both methods allow the quantification and screening of cyanide in blood samples. Furthermore, the test strip is low cost and enables on-site analysis.

Structural insights for producing CK2α1-specific inhibitors.

Casein kinase 2 catalytic subunit (CK2α) is classified into two subtypes CK2α1 and CK2α2. CK2α1 is a drug discovery target, whereas CK2α2 is an off-target of CK2α1 inhibitors. High amino acid sequence homology between these subtypes hampers efforts to produce ATP competitive inhibitors that are highly selective to CK2α1. Hematein was identified previously as a non-ATP-competitive inhibitor for CK2α1, whereas this compound acts as an ATP competitive CK2α2 inhibitor. Crystal structures of CK2α1 and CK2α2 in complex with hematein revealed distinct binding features that provide structural insights for producing CK2α1-selective inhibitors.

An Engineered Aryl Acid Adenylation Domain with an Enlarged Substrate Binding Pocket.

Adenylation (A) domains act as the gatekeepers of non-ribosomal peptide synthetases (NRPSs), ensuring the activation and thioesterification of the correct amino acid/aryl acid building blocks. Aryl acid building blocks are most commonly observed in iron-chelating siderophores, but are not limited to them. Very little is known about the reprogramming of aryl acid A-domains. We show that a single asparagine-to-glycine mutation in an aryl acid A-domain leads to an enzyme that tolerates a wide range of non-native aryl acids. The engineered catalyst is capable of activating non-native aryl acids functionalized with nitro, cyano, bromo, and iodo groups, even though no enzymatic activity of wild-type enzyme was observed toward these substrates. Co-crystal structures with non-hydrolysable aryl-AMP analogues revealed the origins of this expansion of substrate promiscuity, highlighting an enlargement of the substrate binding pocket of the enzyme. Our findings may be exploited to produce diversified aryl acid containing natural products and serve as a template for further directed evolution in combinatorial biosynthesis.

Structure-Activity Relationships and Docking Studies of Hydroxychavicol and Its Analogs as Xanthine Oxidase Inhibitors.

Hydroxychavicol (HC), which is obtained from the leaves of Piper betle LINN. (Piperaceae), inhibits xanthine oxidase (XO) with an IC50 value of 16.7 µM, making it more potent than the clinically used allopurinol (IC50=30.7 µM). Herein, a structure-activity relationship analysis of the polar part analogs of HC was conducted and an inhibitor was discovered with a potency 13 times that of HC. Kinetic studies have revealed that HC and its active analog inhibit XO in an uncompetitive manner. The binding structure prediction of these inhibitor molecules to the XO complex with xanthine suggested that both compounds (HC and its analog) could simultaneously form hydrogen bonds with xanthine and XO.

Synthesis and Functional Characterization of Novel Sialyl LewisX Mimic-Decorated Liposomes for E-selectin-Mediated Targeting to Inflamed Endothelial Cells.

Sialyl LewisX (sLeX) is a natural ligand of E-selectin that is overexpressed by inflamed and tumor endothelium. Although sLeX is a potential ligand for drug targeting, synthesis of the tetrasaccharide is complicated with many reaction steps. In this study, structurally simplified novel sLeX analogues were designed and linked with 1,2-distearoyl-sn-glycero-3-phosphoethanolamine-polyethylene glycol-2000 (DSPE-PEG) for E-selectin-mediated liposomal delivery. The sLeX structural simplification strategies include (1) replacement of the Gal-GlcNAc disaccharide unit with lactose to reduce many initial steps and (2) substitution of neuraminic acid with a negatively charged group, i.e., 3'-sulfo, 3'-carboxymethyl (3'-CM), or 3'-(1-carboxy)ethyl (3'-CE). While all the liposomes developed were similar in particle size and charge, the 3'-CE sLeX mimic liposome demonstrated the highest uptake in inflammatory cytokine-treated human umbilical vein endothelial cells (HUVECs), being even more potent than native sLeX-decorated liposomes. Inhibition studies using antiselectin antibodies revealed that their uptake was mediated primarily by overexpressed E-selectin on inflamed HUVECs. Molecular dynamics simulations were performed to gain mechanistic insight into the E-selectin binding differences among native and mimic sLeX. The terminally branched methyl group of the 3'-CE sLeX mimic oriented and faced the bulk hydrophilic solution during E-selectin binding. Since this state is entropically unfavorable, the 3'-CE sLeX mimic molecule might be pushed toward the binding pocket of E-selectin by a hydrophobic effect, leading to a higher probability of hydrogen-bond formation than native sLeX and the 3'-CM sLeX mimic. This corresponded with the fact that the 3'-CE sLeX mimic liposome exhibited much greater uptake than the 3'-CM sLeX mimic liposome.

Identification of selective inhibitors of sphingosine kinases 1 and 2 through a structure-activity relationship study of 4-epi-jaspine B.

We recently reported that 4-epi-jaspine B exhibits potent inhibitory activity towards sphingosine kinases (SphKs). In this study, we investigated the effects of modifying the 2-alkyl group, as well as the functional groups on the THF ring of 4-epi-jaspine B using a diversity-oriented synthesis approach based on a late-stage cross metathesis reaction. The introduction of a p-phenylene tether to the alkyl group was favored in most cases, whereas the replacement of a carbon atom with an oxygen atom led to a decrease in the inhibitory activity. Furthermore, the introduction of a bulky alkyl group at the terminus led to a slight increase in the inhibitory activity of this series towards SphKs compared with 4-epi-jaspine B (the Q values of compound 13 for SphK1 and SphK2 were 0.2 and 0.4, respectively). Based on this study, we identified two isoform selective inhibitors, including the m-phenylene derivative 4 [IC50 (SphK1) ≥30µM; IC50 (SphK2)=2.2µM] and the methyl ether derivative 22 [IC50 (SphK1)=4.0µM; IC50 (SphK2) ≥30µM].

Structure-activity relationship study of 4-(thiazol-5-yl)benzoic acid derivatives as potent protein kinase CK2 inhibitors.

Two classes of modified analogs of 4-(thiazol-5-yl)benzoic acid-type CK2 inhibitors were designed. The azabenzene analogs, pyridine- and pyridazine-carboxylic acid derivatives, showed potent protein kinase CK2 inhibitory activities [IC50 (CK2α)=0.014-0.017µM; IC50 (CK2α')=0.0046-0.010µM]. Introduction of a 2-halo- or 2-methoxy-benzyloxy group at the 3-position of the benzoic acid moiety maintained the potent CK2 inhibitory activities [IC50 (CK2α)=0.014-0.016µM; IC50 (CK2α')=0.0088-0.014µM] and led to antiproliferative activities [CC50 (A549)=1.5-3.3µM] three to six times higher than those of the parent compound.

Hydrophobic substituents increase the potency of salacinol, a potent α-glucosidase inhibitor from Ayurvedic traditional medicine 'Salacia'.

Using an in silico method, seven analogs bearing hydrophobic substituents (8a: Me, 8b: Et, 8c: n-Pent, 8d: n-Hept, 8e: n-Tridec, 8f: isoBu and 8g: neoPent) at the 3'-O-position in salacinol (1), a highly potent natural α-glucosidase inhibitor from Ayurvedic traditional medicine 'Salacia', were designed and synthesized. In order to verify the computational SAR assessments, their α-glucosidase inhibitory activities were evaluated in vitro. All analogs (8a-8g) exhibited an equal or considerably higher level of inhibitory activity against rat small intestinal α-glucosidases compared with the original sulfonate (1), and were as potent as or higher in potency than the clinically used anti-diabetics, voglibose, acarbose or miglitol. Their activities against human maltase exhibited good relationships to the results obtained with enzymes of rat origin. Among the designed compounds, the one with a 3'-O-neopentyl moiety (8g) was most potent, with an approximately ten fold increase in activity against human maltase compared to 1.

In vitro activity of five quinolones and analysis of the quinolone resistance-determining regions of gyrA, gyrB, parC, and parE in Ureaplasma parvum and Ureaplasma urealyticum clinical isolates from perinatal patients in Japan.

Ureaplasma spp. cause several disorders, such as nongonococcal urethritis, miscarriage, and preterm delivery with lung infections in neonates, characterized by pathological chorioamnionitis in the placenta. Although reports on antibiotic resistance in Ureaplasma are on the rise, reports on quinolone-resistant Ureaplasma infections in Japan are limited. The purpose of this study was to determine susceptibilities to five quinolones of Ureaplasma urealyticum and Ureaplasma parvum isolated from perinatal samples in Japan and to characterize the quinolone resistance-determining regions in the gyrA, gyrB, parC, and parE genes. Out of 28 clinical Ureaplasma strains, we isolated 9 with high MICs of quinolones and found a single parC gene mutation, resulting in the change S83L. Among 158 samples, the ParC S83L mutation was found in 37 samples (23.4%), including 1 sample harboring a ParC S83L-GyrB P462S double mutant. Novel mutations of ureaplasmal ParC (S83W and S84P) were independently found in one of the samples. Homology modeling of the ParC S83W mutant suggested steric hindrance of the quinolone-binding pocket (QBP), and de novo prediction of peptide structures revealed that the ParC S84P may break/kink the formation of the α4 helix in the QBP. Further investigations are required to unravel the extent and mechanism of antibiotic resistance of Ureaplasma spp. in Japan.

The use of many-body expansions and geometry optimizations in fragment-based methods.

Conspectus Chemists routinely work with complex molecular systems: solutions, biochemical molecules, and amorphous and composite materials provide some typical examples. The questions one often asks are what are the driving forces for a chemical phenomenon? How reasonable are our views of chemical systems in terms of subunits, such as functional groups and individual molecules? How can one quantify the difference in physicochemical properties of functional units found in a different chemical environment? Are various effects on functional units in molecular systems additive? Can they be represented by pairwise potentials? Are there effects that cannot be represented in a simple picture of pairwise interactions? How can we obtain quantitative values for these effects? Many of these questions can be formulated in the language of many-body effects. They quantify the properties of subunits (fragments), referred to as one-body properties, pairwise interactions (two-body properties), couplings of two-body interactions described by three-body properties, and so on. By introducing the notion of fragments in the framework of quantum chemistry, one obtains two immense benefits: (a) chemists can finally relate to quantum chemistry, which now speaks their language, by discussing chemically interesting subunits and their interactions and (b) calculations become much faster due to a reduced computational scaling. For instance, the somewhat academic sounding question of the importance of three-body effects in water clusters is actually another way of asking how two hydrogen bonds affect each other, when they involve three water molecules. One aspect of this is the many-body charge transfer (CT), because the charge transfers in the two hydrogen bonds are coupled to each other (not independent). In this work, we provide a generalized view on the use of many-body expansions in fragment-based methods, focusing on the general aspects of the property expansion and a contraction of a many-body expansion in a formally two-body series, as exemplified in the development of the fragment molecular orbital (FMO) method. Fragment-based methods have been very successful in delivering the properties of fragments, as well as the fragment interactions, providing insights into complex chemical processes in large molecular systems. We briefly review geometry optimizations performed with fragment-based methods and present an efficient geometry optimization method based on the combination of FMO with molecular mechanics (MM), applied to the complex of a subunit of protein kinase 2 (CK2) with a ligand. FMO results are discussed in comparison with experimental and MM-optimized structures.

Possible Contribution of a Diverticulum to the Development and Rupture of Colonic Lymphangioma.

A 55-year-old Japanese man with a history of diverticulitis underwent colonoscopy for careful evaluation of progressive anemia. A 5-mm depressed lesion oozing spontaneously was observed at the hepatic flexure. On suspicion of depressed-type of cancer, right-sided hemicolectomy was performed. Histopathological examination indicated a collapsed lymphangioma exactly over a diverticulum, which had previously been complicated diverticulitis. The colonic mucosa and lymphangioma prolapsed beyond the subserosal layer via the muscularis propria defect, resulting in a depressed lesion and mucosal laceration with hemorrhage. This case suggests the contribution of a colonic diverticulum to the development and rupture of lymphangioma, which needed to be distinguished from depressed-type colon cancer.