Homology modeling by the ICM method

High-level nitrofurantoin resistance in a clinical isolate of Klebsiella pneumoniae: a comparative genomics and metabolomics analysis

Nitrofurantoin is a commonly used chemotherapeutic agent in the treatment of uncomplicated urinary tract infections caused by the problematic multidrug resistant Gram-negative pathogen Klebsiella pneumoniae. The present study aims to elucidate the mechanism of nitrofurantoin action and high-level resistance in K. pneumoniae using whole-genome sequencing (WGS), qPCR analysis, mutation structural modeling and untargeted metabolomic analysis. WGS profiling of evolved highly resistant mutants (nitrofurantoin minimum inhibitory concentrations > 256 mg/L) revealed modified expression of several genes related to membrane transport (porin ompK36 and efflux pump regulator oqxR) and nitroreductase activity (ribC and nfsB, involved in nitrofurantoin reduction). Untargeted metabolomics analysis of total metabolites extracted at 1 and 4 h post-nitrofurantoin treatment revealed that exposure to the drug caused a delayed effect on the metabolome which was most pronounced after 4 h. Pathway enrichment analysis illustrated that several complex interrelated metabolic pathways related to nitrofurantoin bacterial killing (aminoacyl-tRNA biosynthesis, purine metabolism, central carbohydrate metabolism, and pantothenate and CoA biosynthesis) and the development of nitrofurantoin resistance (riboflavin metabolism) were significantly perturbed. This study highlights for the first time the key role of efflux pump regulator oqxR in nitrofurantoin resistance and reveals global metabolome perturbations in response to nitrofurantoin, in K. pneumoniae.IMPORTANCEA quest for novel antibiotics and revitalizing older ones (such as nitrofurantoin) for treatment of difficult-to-treat Gram-negative bacterial infections has become increasingly popular. The precise antibacterial activity of nitrofurantoin is still not fully understood. Furthermore, although the prevalence of nitrofurantoin resistance remains low currently, the drug's fast-growing consumption worldwide highlights the need to comprehend the emerging resistance mechanisms. Here, we used multidisciplinary techniques to discern the exact mechanism of nitrofurantoin action and high-level resistance in Klebsiella pneumoniae, a common cause of urinary tract infections for which nitrofurantoin is the recommended treatment. We found that the expression of multiple genes related to membrane transport (including active efflux and passive diffusion of drug molecules) and nitroreductase activity was modified in nitrofurantoin-resistant strains, including oqxR, the transcriptional regulator of the oqxAB efflux pump. Furthermore, complex interconnected metabolic pathways that potentially govern the nitrofurantoin-killing mechanisms (e.g., aminoacyl-tRNA biosynthesis) and nitrofurantoin resistance (riboflavin metabolism) were significantly inhibited following nitrofurantoin treatment. Our study could help inform the improvement of nitrofuran derivatives, the development of new pharmacophores, or drug combinations to support the resurgence of nitrofurantoin in the management of multidrug resistant K. pneumouniae infection.

Software and Databases for Protein-Protein Docking

Protein-protein interactions (PPIs) provide valuable insights for understanding the principles of biological systems and for elucidating causes of incurable diseases. One of the techniques used for computational prediction of PPIs is protein-protein docking calculations, and a variety of software has been developed. This chapter is a summary of software and databases used for protein-protein docking.

Crystal Structures of Bacterial Pectin Methylesterases Pme8A and PmeC2 from Rumen Butyrivibrio

Pectin is a complex polysaccharide that forms a substantial proportion of the plant's middle lamella of forage ingested by grazing ruminants. Methanol in the rumen is derived mainly from methoxy groups released from pectin by the action of pectin methylesterase (PME) and is subsequently used by rumen methylotrophic methanogens that reduce methanol to produce methane (CH4). Members of the genus Butyrivibrio are key pectin-degrading rumen bacteria that contribute to methanol formation and have important roles in fibre breakdown, protein digestion, and the biohydrogenation of fatty acids. Therefore, methanol release from pectin degradation in the rumen is a potential target for CH4 mitigation technologies. Here, we present the crystal structures of PMEs belonging to the carbohydrate esterase family 8 (CE8) from Butyrivibrio proteoclasticus and Butyrivibrio fibrisolvens, determined to a resolution of 2.30 Å. These enzymes, like other PMEs, are right-handed β-helical proteins with a well-defined catalytic site and reaction mechanisms previously defined in insect, plant, and other bacterial pectin methylesterases. Potential substrate binding domains are also defined for the enzymes.

CL-705G: a novel chemical Kir6.2-specific KATP channel opener

Background: KATP channels have diverse roles, including regulation of insulin secretion and blood flow, and protection against biological stress responses and are excellent therapeutic targets. Different subclasses of KATP channels exist in various tissue types due to the unique assemblies of specific pore-forming (Kir6.x) and accessory (SURx) subunits. The majority of pharmacological openers and blockers act by binding to SURx and are poorly selective against the various KATP channel subclasses. Methods and Results: We used 3D models of the Kir6.2/SUR homotetramers based on existing cryo-EM structures of channels in both the open and closed states to identify a potential agonist binding pocket in a functionally critical area of the channel. Computational docking screens of this pocket with the Chembridge Core chemical library of 492,000 drug-like compounds yielded 15 top-ranked "hits", which were tested for activity against KATP channels using patch clamping and thallium (Tl+) flux assays with a Kir6.2/SUR2A HEK-293 stable cell line. Several of the compounds increased Tl+ fluxes. One of them (CL-705G) opened Kir6.2/SUR2A channels with a similar potency as pinacidil (EC50 of 9 µM and 11 μM, respectively). Remarkably, compound CL-705G had no or minimal effects on other Kir channels, including Kir6.1/SUR2B, Kir2.1, or Kir3.1/Kir3.4 channels, or Na+ currents of TE671 medulloblastoma cells. CL-705G activated Kir6.2Δ36 in the presence of SUR2A, but not when expressed by itself. CL-705G activated Kir6.2/SUR2A channels even after PIP2 depletion. The compound has cardioprotective effects in a cellular model of pharmacological preconditioning. It also partially rescued activity of the gating-defective Kir6.2-R301C mutant that is associated with congenital hyperinsulinism. Conclusion: CL-705G is a new Kir6.2 opener with little cross-reactivity with other channels tested, including the structurally similar Kir6.1. This, to our knowledge, is the first Kir-specific channel opener.

AlphaFold2 and its applications in the fields of biology and medicine

AlphaFold2 (AF2) is an artificial intelligence (AI) system developed by DeepMind that can predict three-dimensional (3D) structures of proteins from amino acid sequences with atomic-level accuracy. Protein structure prediction is one of the most challenging problems in computational biology and chemistry, and has puzzled scientists for 50 years. The advent of AF2 presents an unprecedented progress in protein structure prediction and has attracted much attention. Subsequent release of structures of more than 200 million proteins predicted by AF2 further aroused great enthusiasm in the science community, especially in the fields of biology and medicine. AF2 is thought to have a significant impact on structural biology and research areas that need protein structure information, such as drug discovery, protein design, prediction of protein function, et al. Though the time is not long since AF2 was developed, there are already quite a few application studies of AF2 in the fields of biology and medicine, with many of them having preliminarily proved the potential of AF2. To better understand AF2 and promote its applications, we will in this article summarize the principle and system architecture of AF2 as well as the recipe of its success, and particularly focus on reviewing its applications in the fields of biology and medicine. Limitations of current AF2 prediction will also be discussed.

Structure-activity relationship studies of benzothiazole-phenyl analogs as multi-target ligands to alleviate pain without affecting normal behavior

Soluble epoxide hydrolase (sEH) and fatty acid amide hydrolase (FAAH) are potential targets for several diseases. Previous studies have reported that concomitant selective inhibition of sEH and FAAH produced antinociception effects in an animal model of pain. However, the co-administration of a selective sEH inhibitor and a selective FAAH inhibitor might produce serious side effects due to drug-drug interactions that could complicate drug development in the long term. Thus, discovering dual sEH/FAAH inhibitors, single small molecules that can simultaneously inhibit both sEH and FAAH, would be a significant accomplishment in the medicinal chemistry field. Herein, we report the synthesis and biological evaluation of benzothiazole-phenyl-based analogs as potential dual sEH/FAAH inhibitors. This work represents a follow-up structure-activity relationship (SAR) and metabolic-stability studies of our best dual sEH/FAAH inhibitor identified previously, as well as in vivo evaluation of its effects on voluntary locomotor behavior in rats. Our SAR study indicates that trifluoromethyl groups on the aromatic rings are well tolerated by the targeted enzymes when placed at the ortho and para positions; however, they, surprisingly, did not improve metabolic stability in liver microsomes. Our behavioral studies indicate that doses of dual sEH/FAAH inhibitors that alleviate pain do not depress voluntary behavior in naïve rats, which is a common side effect of currently available analgesic drugs (e.g., opioids). Thus, dual sEH/FAAH inhibitors may be a safe and effective approach to treat pain.

Allosteric regulation of CAD modulates de novo pyrimidine synthesis during the cell cycle

Metabolism is a fundamental cellular process that is coordinated with cell cycle progression. Despite this association, a mechanistic understanding of cell cycle phase-dependent metabolic pathway regulation remains elusive. Here we report the mechanism by which human de novo pyrimidine biosynthesis is allosterically regulated during the cell cycle. Combining traditional synchronization methods and metabolomics, we characterize metabolites by their accumulation pattern during cell cycle phases and identify cell cycle phase-dependent regulation of carbamoyl-phosphate synthetase 2, aspartate transcarbamylase and dihydroorotase (CAD), the first, rate-limiting enzyme in de novo pyrimidine biosynthesis. Through systematic mutational scanning and structural modelling, we find allostery as a major regulatory mechanism that controls the activity change of CAD during the cell cycle. Specifically, we report evidence of two Animalia-specific loops in the CAD allosteric domain that involve sensing and binding of uridine 5'-triphosphate, a CAD allosteric inhibitor. Based on homology with a mitochondrial carbamoyl-phosphate synthetase homologue, we identify a critical role for a signal transmission loop in regulating the formation of a substrate channel, thereby controlling CAD activity.

Multi-Omic Profiling, Structural Characterization, and Potent Inhibitor Screening of Evasion-Related Proteins of a Parasitic Nematode, Haemonchus contortus, Surviving Vaccine Treatment

The emergence of drug-resistant parasitic nematodes in both humans and livestock calls for development of alternative and cost-effective control strategies. Barbervax^® is the only registered vaccine for the economically important ruminant strongylid Haemonchus contortus. In this study, we compared the microbiome, genome-wide diversity, and transcriptome of H. contortus adult male populations that survived vaccination with an experimental vaccine after inoculation in sheep. Our genome-wide SNP analysis revealed 16 putative candidate vaccine evasion genes. However, we did not identify any evidence for changes in microbial community profiling based on the 16S rRNA gene sequencing results of the vaccine-surviving parasite populations. A total of fifty-eight genes were identified as significantly differentially expressed, with six genes being long non-coding (lnc) RNAs and none being putative candidate SNP-associated genes. The genes that highly upregulated in surviving parasites from vaccinated animals were associated with GO terms belonging to predominantly molecular functions and a few biological processes that may have facilitated evasion or potentially lessened the effect of the vaccine. These included five targets: astacin (ASTL), carbonate dehydratase (CA2), phospholipase A2 (PLA2), glutamine synthetase (GLUL), and fatty acid-binding protein (FABP3). Our tertiary structure predictions and modelling analyses were used to perform in silico searches of all published and commercially available inhibitor molecules or substrate analogs with potential broad-spectrum efficacy against nematodes of human and veterinary importance.

Homology Modeling of Transporter Proteins

Membrane transporter proteins are divided into channels/pores and carriers and constitute protein families of physiological and pharmacological importance. Several presently used therapeutic compounds elucidate their effects by targeting membrane transporter proteins, including anti-arrhythmic, anesthetic, antidepressant, anxiolytic and diuretic drugs. The lack of three-dimensional structures of human transporters hampers experimental studies and drug discovery. In this chapter, the use of homology modeling for generating structural models of membrane transporter proteins is reviewed. The increasing number of atomic resolution structures available as templates, together with improvements in methods and algorithms for sequence alignments, secondary structure predictions, and model generation, in addition to the increase in computational power have increased the applicability of homology modeling for generating structural models of transporter proteins. Different pitfalls and hints for template selection, multiple-sequence alignments, generation and optimization, validation of the models, and the use of transporter homology models for structure-based virtual ligand screening are discussed.

Control of cell surface expression of GABAA receptors by a conserved region at the end of the N-terminal extracellular domain of receptor subunits

Type A γ-aminobutyric acid receptors (GABA_ARs) represent a family of pentameric GABA-gated Cl^-/HCO₃^- ion channels which mediate inhibitory transmission in the central nervous system. Cell surface expression of GABA_ARs, a prerequisite for their function, is dependent on the appropriate assembly of the receptor subunits and their transient interactions with molecular chaperones within the endoplasmic reticulum (ER) and Golgi apparatus. Here, we describe a highly conserved amino acid sequence within the extracellular N-terminal domain of the receptor subunits adjoining the first transmembrane domain as a region important for GABA_AR processing within the ER. Modifications of this region in the α1, β3, and γ2 subunits using insertion or site-directed mutagenesis impaired GABA_AR trafficking to the cell surface in heterologous cell systems although they had no effect on the subunit assembly. We found that mutated receptors accumulated in the ER where they were shown to associate with chaperones calnexin, BiP, and Grp94. However, their surface expression was increased when ER-associated degradation or proteosome function was inhibited, while modulation of ER calcium stores had little effect. When compared to the wt, mutated receptors showed decreased interaction with calnexin, similar binding to BiP, and increased association with Grp94. Structural modeling of calnexin interaction with the wt or mutated GABA_AR revealed that disruption in structure caused by mutations in the conserved region adjoining the first transmembrane domain may impair calnexin binding. Thus, this previously uncharacterized region plays an important role in intracellular processing of GABA_ARs at least in part by stabilizing their interaction with calnexin.

Characterization of functionally deficient SIM2 variants found in patients with neurological phenotypes

Single-Minded 2 (SIM2) is a neuron enriched basic Helix-Loop-Helix/PER-ARNT-SIM (bHLH/PAS) transcription factor essential for mammalian survival. SIM2 is located within the Down Syndrome Critical Region (DSCR) of chromosome 21, and manipulation in mouse models suggests Sim2 may play a role in brain development and function. During screening of a clinical exome sequencing database, nine SIM2 non-synonymous mutations were found which were subsequently investigated for impaired function using cell-based reporter gene assays. A number of these human variants attenuated abilities to activate transcription and were further characterized to determine the mechanisms underpinning their deficiencies. These included impaired partner protein dimerization, reduced DNA binding and reduced expression and nuclear localization. This study highlighted several SIM2 variants found in patients with disabilities and validated a candidate set as potentially contributing to pathology.

Identification and Functional Characterization of Peptides With Antimicrobial Activity From the Syphilis Spirochete, Treponema pallidum

The etiological agent of syphilis, Treponema pallidum ssp. pallidum, is a highly invasive “stealth” pathogen that can evade the host immune response and persist within the host for decades. This obligate human pathogen is adept at establishing infection and surviving at sites within the host that have a multitude of competing microbes, sometimes including pathogens. One survival strategy employed by bacteria found at polymicrobial sites is elimination of competing microorganisms by production of antimicrobial peptides (AMPs). Antimicrobial peptides are low molecular weight proteins (miniproteins) that function directly via inhibition and killing of microbes and/or indirectly via modulation of the host immune response, which can facilitate immune evasion. In the current study, we used bioinformatics to show that approximately 7% of the T. pallidum proteome is comprised of miniproteins of 150 amino acids or less with unknown functions. To investigate the possibility that AMP production is an unrecognized defense strategy used by T. pallidum during infection, we developed a bioinformatics pipeline to analyze the complement of T. pallidum miniproteins of unknown function for the identification of potential AMPs. This analysis identified 45 T. pallidum AMP candidates; of these, Tp0451a and Tp0749 were subjected to further bioinformatic analyses to identify AMP critical core regions (AMPCCRs). Four potential AMPCCRs from the two predicted AMPs were identified and peptides corresponding to these AMPCCRs were experimentally confirmed to exhibit bacteriostatic and bactericidal activity against a panel of biologically relevant Gram-positive and Gram-negative bacteria. Immunomodulation assays performed under inflammatory conditions demonstrated that one of the AMPCCRs was also capable of differentially regulating expression of two pro-inflammatory chemokines [monocyte chemoattractant protein-1 (MCP-1) and interleukin-8 (IL-8)]. These findings demonstrate proof-of-concept for our developed AMP identification pipeline and are consistent with the novel concept that T. pallidum expresses AMPs to defend against competing microbes and modulate the host immune response.

Computational Prediction of Phosphoinositide Binding to Hyperpolarization-Activated Cyclic-Nucleotide Gated Channels

Protein-lipid interactions are key regulators of ion channel function. Numerous ion channels, including hyperpolarization-activated cyclic-nucleotide gated (HCN) channels have been shown to be regulated by phosphoinositides (PIPs), with important implications in cardiac and neuronal function. Specifically, PIPs have been shown to enhance HCN activation. Using computational approaches, we aim to identify potential binding sites for HCN1-PIP interactions. Computational docking and coarse-grained simulations indicate that PIP binding to HCN1 channels is not well coordinated, but rather occurs over a broad surface of charged residues primarily in the HCN-domain, S2 and S3 helices that can be loosely organized in 2 or 3 overlapping clusters. Thus, PIP-HCN1 interactions are more resembling of electrostatic interactions that occur in myristoylated alanine-rich C kinase substrate (MARCKS) proteins, than the specifically coordinated interactions that occur in pleckstrin homology domains (PH domains) or ion channels such as inward rectifier potassium (Kir) channels. Our results also indicate that phosphatidylinositol (PI) interactions with HCN1 are even lower affinity, explaining why unphosphorylated PI have no effect on HCN1 activation unlike phosphorylated PIPs.

Biophysical Compatibility of a Heterotrimeric Tyrosinase-TYRP1-TYRP2 Metalloenzyme Complex

Tyrosinase (TYR) is a copper-containing monooxygenase central to the function of melanocytes. Alterations in its expression or activity contribute to variations in skin, hair and eye color, and underlie a variety of pathogenic pigmentary phenotypes, including several forms of oculocutaneous albinism (OCA). Many of these phenotypes are linked to individual missense mutations causing single nucleotide variants and polymorphisms (SNVs) in TYR. We previously showed that two TYR homologues, TYRP1 and TYRP2, modulate TYR activity and stabilize the TYR protein. Accordingly, to investigate whether TYR, TYRP1, and TYRP2 are biophysically compatible with various heterocomplexes, we computationally docked a high-quality 3D model of TYR to the crystal structure of TYRP1 and to a high-quality 3D model of TYRP2. Remarkably, the resulting TYR-TYRP1 heterodimer was complementary in structure and energy with the TYR-TYRP2 heterodimer, with TYRP1 and TYRP2 docking to different adjacent surfaces on TYR that apposed a third realistic protein interface between TYRP1-TYRP2. Hence, the 3D models are compatible with a heterotrimeric TYR-TYRP1-TYRP2 complex. In addition, this heterotrimeric TYR-TYRP1-TYRP2 positioned the C-terminus of each folded enzymatic domain in an ideal position to allow their C-terminal transmembrane helices to form a putative membrane embedded three-helix bundle. Finally, pathogenic TYR mutations causing OCA1A, which also destabilize TYR biochemically, cluster on an unoccupied protein interface at the periphery of the heterotrimeric complex, suggesting that this may be a docking site for OCA2, an anion channel. Pathogenic OCA2 mutations result in similar phenotypes to those produced by OCA1A TYR mutations. While this complex may be difficult to detect in vitro, due to the complex environment of the vertebrate cellular membranous system, our results support the existence of a heterotrimeric complex in melanogenesis.

Anti-V2 antibodies virus vulnerability revealed by envelope V1 deletion in HIV vaccine candidates

The efficacy of ALVAC-based HIV and SIV vaccines in humans and macaques correlates with antibodies to envelope variable region 2 (V2). We show here that vaccine-induced antibodies to SIV variable region 1 (V1) inhibit anti-V2 antibody-mediated cytotoxicity and reverse their ability to block V2 peptide interaction with the α4β7 integrin. SIV vaccines engineered to delete V1 and favor an α helix, rather than a β sheet V2 conformation, induced V2-specific ADCC correlating with decreased risk of SIV acquisition. Removal of V1 from the HIV-1 clade A/E A244 envelope resulted in decreased binding to antibodies recognizing V2 in the β sheet conformation. Thus, deletion of V1 in HIV envelope immunogens may improve antibody responses to V2 virus vulnerability sites and increase the efficacy of HIV vaccine candidates.

A receptor-like protein mediates plant immune responses to herbivore-associated molecular patterns

Herbivory is fundamental to the regulation of both global food webs and the extent of agricultural crop losses. Induced plant responses to herbivores promote resistance and often involve the perception of specific herbivore-associated molecular patterns (HAMPs); however, precisely defined receptors and elicitors associated with herbivore recognition remain elusive. Here, we show that a receptor confers signaling and defense outputs in response to a defined HAMP common in caterpillar oral secretions (OS). Staple food crops, including cowpea (Vigna unguiculata) and common bean (Phaseolus vulgaris), specifically respond to OS via recognition of proteolytic fragments of chloroplastic ATP synthase, termed inceptins. Using forward-genetic mapping of inceptin-induced plant responses, we identified a corresponding leucine-rich repeat receptor, termed INR, specific to select legume species and sufficient to confer inceptin-induced responses and enhanced defense against armyworms (Spodoptera exigua) in tobacco. Our results support the role of plant immune receptors in the perception of chewing herbivores and defense.

Differential activities of maize plant elicitor peptides as mediators of immune signaling and herbivore resistance

Plant elicitor peptides (Peps) are conserved regulators of defense responses and models for the study of damage-associated molecular pattern-induced immunity. Although present as multigene families in most species, the functional relevance of these multigene families remains largely undefined. While Arabidopsis Peps appear largely redundant in function, previous work examining Pep-induced responses in maize (Zm) implied specificity of function. To better define the function of individual ZmPeps and their cognate receptors (ZmPEPRs), activities were examined by assessing changes in defense-associated phytohormones, specialized metabolites and global gene expression patterns, in combination with heterologous expression assays and analyses of CRISPR/Cas9-generated knockout plants. Beyond simply delineating individual ZmPep and ZmPEPR activities, these experiments led to a number of new insights into Pep signaling mechanisms. ZmPROPEP and other poaceous precursors were found to contain multiple active Peps, a phenomenon not previously observed for this family. In all, seven new ZmPeps were identified and the peptides were found to have specific activities defined by the relative magnitude of their response output rather than by uniqueness. A striking correlation was observed between individual ZmPep-elicited changes in levels of jasmonic acid and ethylene and the magnitude of induced defense responses, indicating that ZmPeps may collectively regulate immune output through rheostat-like tuning of phytohormone levels. Peptide structure-function studies and ligand-receptor modeling revealed structural features critical to the function of ZmPeps and led to the identification of ZmPep5a as a potential antagonist peptide able to competitively inhibit the activity of other ZmPeps, a regulatory mechanism not previously observed for this family.

Riproximin Exhibits Diversity in Sugar Binding, and Modulates some Metastasis-Related Proteins with Lectin like Properties in Pancreatic Ductal Adenocarcinoma

Riproximin (Rpx) is a type II ribosome-inactivating protein with specific anti-proliferative activity. It was purified from Ximenia americana by affinity chromatography using a resin coupled with lactosyl residues. The same technique facilitated isolation of proteins with lectin-like properties from human Suit2-007 and rat ASML pancreatic cancer cells, which were termed lactosyl-sepharose binding proteins (LSBPs). The role of these proteins in cancer progression was investigated at mRNA level using chip array data of Suit2-007 and ASML cells re-isolated from nude rats. These data compared significant mRNA expression changes when relating primary (pancreas) and metastatic (liver) sites following orthotopic and intraportal implantation of Pancreatic Ductal Adenocarcinoma (PDAC) cells, respectively. The affinity of Rpx to 13 simple sugar structures was modeled by docking experiments, the ranking of which was principally confirmed by NMR-spectroscopy. In addition, Rpx and LSBPs were evaluated for anti-proliferative activity and their cellular uptake was assessed by fluorescence microscopy. From 13 monosaccharides evaluated, open-chain rhamnose, β-d-galactose, and α-l-galactopyranose showed the highest affinities for site 1 of Rpx’s B-chain. NMR evaluation yielded a similar ranking, as galactose was among the best binders. Both, Rpx and LSBPs reduced cell proliferation in vitro, but their anti-proliferative effects were decreased by 15–20% in the presence of galactose. The program “Ingenuity Pathway Analysis” identified 2,415 genes showing significantly modulated mRNA expression following exposure of Suit2-007 cells to Rpx in vitro. These genes were then matched to those 1,639 genes, which were significantly modulated in the rat model when comparing primary and metastatic growth of Suit2-007 cells. In this overlap analysis, LSBP genes were considered separately. The potential suitability of Rpx for treating metastatic Suit2-007 PDAC cells was reflected by those genes, which were modulated by Rpx in a way opposite to that observed in cancer progression. Remarkably, these were 14% of all genes modulated during cancer progression, but 71% of the respective LSBP gene subgroup. Based on these findings, we predict that Rpx has the potential to treat PDAC metastasis by modulating genes involved in metastatic progression, especially by targeting LSBPs.

Dynamin-2 R465W mutation induces long range perturbation in highly ordered oligomeric structures

High order oligomers are crucial for normal cell physiology, and protein function perturbed by missense mutations underlies several autosomal dominant diseases. Dynamin-2 is one of such protein forming helical oligomers that catalyze membrane fission. Mutations in this protein, where R465W is the most frequent, cause dominant centronuclear myopathy, but the molecular mechanisms underpinning the functional modifications remain to be investigated. To unveil the structural impact of this mutation in dynamin-2, we used full-atom molecular dynamics simulations and coarse-grained models and built dimers and helices of wild-type (WT) monomers, mutant monomers, or both WT and mutant monomers combined. Our results show that the mutation R465W causes changes in the interactions with neighbor amino acids that propagate through the oligomer. These new interactions perturb the contact between monomers and favor an extended conformation of the bundle signaling element (BSE), a dynamin region that transmits the conformational changes from the GTPase domain to the rest of the protein. This extended configuration of the BSE that is only relevant in the helices illustrates how a small change in the microenvironment surrounding a single residue can propagate through the oligomer structures of dynamin explaining how dominance emerges in large protein complexes.

Proline Hydroxylation Primes Protein Kinases for Autophosphorylation and Activation

Activation of dual-specificity tyrosine-phosphorylation-regulated kinases 1A and 1B (DYRK1A and DYRK1B) requires prolyl hydroxylation by PHD1 prolyl hydroxylase. Prolyl hydroxylation of DYRK1 initiates a cascade of events leading to the release of molecular constraints on von Hippel-Lindau (VHL) ubiquitin ligase tumor suppressor function. However, the proline residue of DYRK1 targeted by hydroxylation and the role of prolyl hydroxylation in tyrosine autophosphorylation of DYRK1 are unknown. We found that a highly conserved proline in the CMGC insert of the DYRK1 kinase domain is hydroxylated by PHD1, and this event precedes tyrosine autophosphorylation. Mutation of the hydroxylation acceptor proline precludes tyrosine autophosphorylation and folding of DYRK1, resulting in a kinase unable to preserve VHL function and lacking glioma suppression activity. The consensus proline sequence is shared by most CMGC kinases, and prolyl hydroxylation is essential for catalytic activation. Thus, formation of prolyl-hydroxylated intermediates is a novel mechanism of kinase maturation and likely a general mechanism of regulation of CMGC kinases in eukaryotes.

Theoretical investigation of AhR binding property with relevant structural requirements for AhR-mediated toxicity of polybrominated diphenyl ethers

Polybrominated diphenyl ethers (PBDEs) are more frequently suspected with the induction of toxicity via signal transduction pathway of cytosolic aryl hydrocarbon receptor (AhR), the initial binding to which is assumed to be an essential prerequisite during the ligand-dependent activation. However, the AhR binding property and associated toxicity of PBDEs is yet to be clearly known for lacking insights into the structural requirements at molecular level. To understand the AhR binding property of PBDEs, the ligand binding domain (LBD) of AhR was simulatively developed on homologous protein after basic validation of geometrical rationality and the binding interaction profile was visually described using molecular docking approach. For AhR binding, the offset or edge-on π-π stackings with aromatic motifs including Phe289, Phe345 and His285 were shown to be structurally required whereas the electrostatic attraction validated for AhR binding to dioxins might be less effective for 2,2',3,4,4'-pentabromodiphenyl ether (BDE-85). Besides the demands of less steric hindrance from alanines and weak formulation of hydrogen bonds, the dispersion force through large contact and polarization of S-π electrons seemed to be impactful when BDE-85 were closer to Cys327, Met334 or Met342. With theoretical computation of AhR binding energies, the more significant correlativity with bioassays was derived especially for the lowly/moderately brominated congeners, and could be used to predict the AhR binding affinity on certain degree. The informative results would thus not only help well understand the molecular basis of AhR-mediated toxicity but give an approach for accelerative evaluation of AhR binding and toxicity of PBDEs.

Palmitoylation of the KATP channel Kir6.2 subunit promotes channel opening by regulating PIP2 sensitivity

A physiological role for long-chain acyl-CoA esters to activate ATP-sensitive K⁺ (K_ATP) channels is well established. Circulating palmitate is transported into cells and converted to palmitoyl-CoA, which is a substrate for palmitoylation. We found that palmitoyl-CoA, but not palmitic acid, activated the channel when applied acutely. We have altered the palmitoylation state by preincubating cells with micromolar concentrations of palmitic acid or by inhibiting protein thioesterases. With acyl-biotin exchange assays we found that Kir6.2, but not sulfonylurea receptor (SUR)1 or SUR2, was palmitoylated. These interventions increased the K_ATP channel mean patch current, increased the open time, and decreased the apparent sensitivity to ATP without affecting surface expression. Similar data were obtained in transfected cells, rat insulin-secreting INS-1 cells, and isolated cardiac myocytes. Kir6.2ΔC36, expressed without SUR, was also positively regulated by palmitoylation. Mutagenesis of Kir6.2 Cys¹⁶⁶ prevented these effects. Clinical variants in KCNJ11 that affect Cys¹⁶⁶ had a similar gain-of-function phenotype, but was more pronounced. Molecular modeling studies suggested that palmitoyl-C166 and selected large hydrophobic mutations make direct hydrophobic contact with Kir6.2-bound PIP₂ Patch-clamp studies confirmed that palmitoylation of Kir6.2 at Cys¹⁶⁶ enhanced the PIP₂ sensitivity of the channel. Physiological relevance is suggested since palmitoylation blunted the regulation of K_ATP channels by α1-adrenoreceptor stimulation. The Cys¹⁶⁶ residue is conserved in some other Kir family members (Kir6.1 and Kir3, but not Kir2), which are also subject to regulated palmitoylation, suggesting a general mechanism to control the open state of certain Kir channels.

Discovery of a multipotent chaperone, 1-(2,6-Difluorobenzylamino)-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol with the inhibitory effects on the proliferation of prion, cancer as well as influenza virus

We previously discovered three carbazole derivatives, GJP14 (1-piperidinylmethyl-2-(1-oxo-6-methyl-1,2,3,4-tetrahydrocarbazol-9-yl)-ethan-1-ol) with anti-prion activity, GJC29 (benzylamino-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol) with anti-cancer activity, and THC19 (1-piperidinylmethyl-2-(1,2,3,4-tetrahydrocarnazol-9-yl)-ethan-1-ol) with anti-influenza virus activity. During optimization of GJP14 for the anti-prion activity, we discovered a compound, 1-(2,6-difluorobenzylamino)-3-(1,2,3,4-tetrahydrocarbazol-9-yl)-propan-2-ol, termed 5Y, had the most strong anti-prion activity among a series of newly synthesized derivatives. Intriguingly, we noticed that 5Y had also the most strong anti-colon cancer as well as the anti-influenza virus activities among derivatives. No significant toxicity of 5Y was observed. These results demonstrate that 5Y is a multipotent lead compound with unusually wide spectrum, and may be applicable to therapeutics targeting multiple diseases.Abbreviations: MoPrP: mouse prion protein of amino acid residues of 23-231; PrP^C: cellular form of prion protein; PrP^Sc: scrapie form of prion protein.

An Improved Method for Modeling Voltage-Gated Ion Channels at Atomic Accuracy Applied to Human Cav Channels

Voltage-gated ion channels (VGICs) are associated with hundreds of human diseases. To date, 3D structural models of human VGICs have not been reported. We developed a 3D structural integrity metric to rank the accuracy of all VGIC structures deposited in the PDB. The metric revealed inaccuracies in structural models built from recent single-particle, non-crystalline cryo-electron microscopy maps and enabled the building of highly accurate homology models of human Ca_v channel α₁ subunits at atomic resolution. Human Ca_v Mendelian mutations mostly located to segments involved in the mechanism of voltage sensing and gating within the 3D structure, with multiple mutations targeting equivalent 3D structural locations despite eliciting distinct clinical phenotypes. The models also revealed that the architecture of the ion selectivity filter is highly conserved from bacteria to humans and between sodium and calcium VGICs.

Characterization of drug binding within the HCN1 channel pore

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels mediate rhythmic electrical activity of cardiac pacemaker cells, and in neurons play important roles in setting resting membrane potentials, dendritic integration, neuronal pacemaking, and establishing action potential threshold. Block of HCN channels slows the heart rate and is currently used to treat angina. However, HCN block also provides a promising approach to the treatment of neuronal disorders including epilepsy and neuropathic pain. While several molecules that block HCN channels have been identified, including clonidine and its derivative alinidine, lidocaine, mepivacaine, bupivacaine, ZD7288, ivabradine, zatebradine, and cilobradine, their low affinity and lack of specificity prevents wide-spread use. Different studies suggest that the binding sites of these inhibitors are located in the inner vestibule of HCN channels, but the molecular details of their binding remain unknown. We used computational docking experiments to assess the binding sites and mode of binding of these inhibitors against the recently solved atomic structure of human HCN1 channels, and a homology model of the open pore derived from a closely related CNG channel. We identify a possible hydrophobic groove in the pore cavity that plays an important role in conformationally restricting the location and orientation of drugs bound to the inner vestibule. Our results also help explain the molecular basis of the low-affinity binding of these inhibitors, paving the way for the development of higher affinity molecules.

Homology modeling in drug discovery: Overview, current applications, and future perspectives

Homology modeling is one of the computational structure prediction methods that are used to determine protein 3D structure from its amino acid sequence. It is considered to be the most accurate of the computational structure prediction methods. It consists of multiple steps that are straightforward and easy to apply. There are many tools and servers that are used for homology modeling. There is no single modeling program or server which is superior in every aspect to others. Since the functionality of the model depends on the quality of the generated protein 3D structure, maximizing the quality of homology modeling is crucial. Homology modeling has many applications in the drug discovery process. Since drugs interact with receptors that consist mainly of proteins, protein 3D structure determination, and thus homology modeling is important in drug discovery. Accordingly, there has been the clarification of protein interactions using 3D structures of proteins that are built with homology modeling. This contributes to the identification of novel drug candidates. Homology modeling plays an important role in making drug discovery faster, easier, cheaper, and more practical. As new modeling methods and combinations are introduced, the scope of its applications widens.

Triazolo[4,5-d]pyrimidines as Validated General Control Nonderepressible 2 (GCN2) Protein Kinase Inhibitors Reduce Growth of Leukemia Cells

Cellular stress signals activate adaptive signaling pathways of the mammalian integrated stress response (ISR), of which the unfolded protein response (UPR) is a subset. These pathways converge at the phosporylation of eIF2α. Drug-like, potent and selective chemical inhibitors (valid chemical probes) targeting major ISR kinases have been previously identified, with the exception of GCN2. We synthesized and evaluated a series of GCN2 inhibitors based on a triazolo[4,5-d]pyrimidine scaffold. Several compounds potently inhibited GCN2 in vitro and displayed good selectivity over the related kinases PERK, HRI, and IRE1. The compounds inhibited phosporylation of eIF2α in HEK293T cells with an IC₅₀ < 150 nM, validating them as chemical probes for cellular studies. These probes were screened against the National Cancer Institute NCI-60 human cancer cell line panel. Uniform growth inhibition was observed in the leukemia group of cell lines. Growth inhibition in the most sensitive cell lines coincided with high GCN2 mRNA expression levels. Oncomine analysis revealed high GCN2 expression accompanied by lower asparagine synthetase (ASNS) expression in patient-derived acute lymphoblastic leukemias with B-Cell origins (B-ALL) as well. Notably, asparaginase, which depletes amino acids and triggers GCN2 activity, is a licensed, first-line B-ALL treatment. Thus, we hypothesize that leukemias exhibiting high GCN2 expression and low ASNS expression may be susceptible to pharmacologic GCN2 inhibition.

Protein-RNA Docking Using ICM

Protein-RNA interactions play an important role in many biological processes. Computational methods such as docking have been developed to complement existing biophysical and structural biology techniques. Computational prediction of protein-RNA complex structures includes two steps: generating candidate structures from the individual protein and RNA parts and scoring the generated poses to pick out the correct one. In this work, we considered three recently developed data sets of protein-RNA complexes to evaluate and improve the performance of the FFT-based rigid-body docking algorithm implemented in the ICM package. An electrostatic term describing interactions between negatively charged phosphate groups and positively charged protein residues was added to the energy function used during the docking step to take into account the greater role that electrostatic interactions play in protein-RNA complexes. Next, the docking results were used to optimize a scoring function including van der Waals, electrostatic, and solvation terms. This optimization yielded a much smaller weight for the solvation term indicating that solvation energy may be less important for the scoring of protein-RNA structures. Rescoring of the generated poses with the new scoring function led to much higher success rates, while pose clustering by contact fingerprints produced further improvements, achieving a success rate of 0.66 for the top 100 structures.

Novel piperidine-derived amide sEH inhibitors as mediators of lipid metabolism with improved stability

We have previously identified and reported several potent piperidine-derived amide inhibitors of the human soluble epoxide hydrolase (sEH) enzyme. The inhibition of this enzyme leads to elevated levels of epoxyeicosatrienoic acids (EETs), which are known to possess anti-inflammatory, vasodilatory, and anti-fibrotic effects. Herein, we report the synthesis of 9 analogs of the lead sEH inhibitor and the follow-up structure-activity relationship and liver microsome stability studies. Our findings show that isosteric modifications that lead to significant alterations in the steric and electronic properties at a specific position in the molecule can reduce the efficacy by up to 75-fold. On the other hand, substituting hydrogen with deuterium produces a notable increase (∼30%) in the molecules' half-lives in both rat and human microsomes, while maintaining sEH inhibition potency. These data highlight the utility of isosteric replacement for improving bioavailability, and the newly-synthesized inhibitor structures may thus, serve as a starting point for preclinical development. Our docking study reveals that in the catalytic pocket of sEH, these analogs are in proximity of the key amino acids involved in hydrolysis of EETs.

Convergence of Wnt, growth factor, and heterotrimeric G protein signals on the guanine nucleotide exchange factor Daple

Cellular proliferation, differentiation, and morphogenesis are shaped by multiple signaling cascades, and their dysregulation plays an integral role in cancer progression. Three cascades that contribute to oncogenic potential are those mediated by Wnt proteins and the receptor Frizzled (FZD), growth factor receptor tyrosine kinases (RTKs), and heterotrimeric G proteins and associated GPCRs. Daple is a guanine nucleotide exchange factor (GEF) for the G protein G_αi Daple also binds to FZD and the Wnt/FZD mediator Dishevelled (Dvl), and it enhances β-catenin-independent Wnt signaling in response to Wnt5a-FZD7 signaling. We identified Daple as a substrate of multiple RTKs and non-RTKs and, hence, as a point of convergence for the three cascades. We found that phosphorylation near the Dvl-binding motif in Daple by both RTKs and non-RTKs caused Daple/Dvl complex dissociation and augmented the ability of Daple to bind to and activate G_αi, which potentiated β-catenin-independent Wnt signals and stimulated epithelial-mesenchymal transition (EMT) similarly to Wnt5a/FZD7 signaling. Although Daple acts as a tumor suppressor in the healthy colon, the concurrent increased abundance of Daple and epidermal growth factor receptor (EGFR) in colorectal tumors was associated with poor patient prognosis. Thus, the Daple-dependent activation of G_αi and the Daple-dependent enhancement of β-catenin-independent Wnt signals are not only stimulated by Wnt5a/FZD7 to suppress tumorigenesis but also hijacked by growth factor-activated RTKs to enhance tumor progression. These findings identify a cross-talk paradigm among growth factor RTKs, heterotrimeric G proteins, and the Wnt/FZD pathway in cancer.

Targeting of AMP-activated protein kinase: prospects for computer-aided drug design

INTRODUCTION

Dysregulation of energy homeostasis has been implicated in a number of human chronic diseases including diabetes, obesity, cancer, and inflammation. Given the functional attributes as a central regulator of energy homeostasis, AMP-activated protein kinase (AMPK) is emerging as a therapeutic target for these diseases, and lines of evidence have highlighted the need for rational and robust screening systems for identifying specific AMPK modulators with a therapeutic potential for preventing and/or curing these diseases. Areas covered: Here, the authors review the recent advances in the understanding of three-dimensional structures of AMPK in relationship with the regulatory mechanisms, potentials of AMPK as a therapeutic target in human chronic diseases, and prospects of computer-based drug design for AMPK. Expert opinion: Accumulating information of AMPK structure has provided us with deep insight into the molecular basis underlying the regulatory mechanisms, and further discloses several structural domains, which can be served for a target site for computer-based drug design. Molecular docking and simulations provides useful information about the binding sites between potent drugs and AMPK as well as a rational screening format to discover isoform-specific AMPK modulators. For these reasons, the authors suggest that computer-aided virtual screening methods hold promise as a rational approach for discovering more specific AMPK modulators.

AMP-activated protein kinase fortifies epithelial tight junctions during energetic stress via its effector GIV/Girdin

Loss of epithelial polarity impacts organ development and function; it is also oncogenic. AMPK, a key sensor of metabolic stress stabilizes cell-cell junctions and maintains epithelial polarity; its activation by Metformin protects the epithelial barrier against stress and suppresses tumorigenesis. How AMPK protects the epithelium remains unknown. Here, we identify GIV/Girdin as a novel effector of AMPK, whose phosphorylation at a single site is both necessary and sufficient for strengthening mammalian epithelial tight junctions and preserving cell polarity and barrier function in the face of energetic stress. Expression of an oncogenic mutant of GIV (cataloged in TCGA) that cannot be phosphorylated by AMPK increased anchorage-independent growth of tumor cells and helped these cells to evade the tumor-suppressive action of Metformin. This work defines a fundamental homeostatic mechanism by which the AMPK-GIV axis reinforces cell junctions against stress-induced collapse and also provides mechanistic insight into the tumor-suppressive action of Metformin.

The antiandrogen flutamide is a novel aryl hydrocarbon receptor ligand that disrupts bile acid homeostasis in mice through induction of Abcc4

Flutamide (FLU), an oral, nonsteroidal antiandrogen drug used in the treatment of prostate cancer, is associated with idiosyncratic hepatotoxicity that sometimes causes severe liver damage, including cholestasis, jaundice, and liver necrosis. To understand the mechanism of toxicity, a combination of aryl hydrocarbon receptor (Ahr)-deficient (Ahr^-/-) mice, primary hepatocytes, luciferase reporter gene assays, in silico ligand docking and ultra-performance chromatography-quadrupole time-of-flight mass spectrometry-based metabolomics was used. A significant increase of liver weights, and liver and serum bile acid levels was observed after FLU treatment, indicating hepatomegaly and disrupted bile acid homeostasis. Expression of the AhR gene battery was markedly increased in livers of wild-type mice Ahr^+/+ treated with FLU, while no change was noted in Ahr^-/- mice. Messenger RNAs encoded by AhR target genes were induced in primary mouse hepatocytes cultured with FLU, which confirmed the in vivo results. Ligand-docking analysis further predicted that FLU is an AhR agonist ligand which was confirmed by luciferase reporter gene assays. Multivariate data analysis showed that bile acids were responsible for the separation of vehicle- and FLU-treated Ahr^+/+ mice, while there was no separation in Ahr^-/- mice. Expression of mRNA encoding the bile acid transporter ABCC4 was increased and farnesoid X receptor signaling was inhibited in the livers of Ahr^+/+ mice, but not in Ahr^-/- mice treated with FLU, in agreement with the observed downstream metabolic alterations. These findings provide new insights into the mechanism of liver injury caused by FLU treatment involving activation of AhR and the alterations of bile acid homeostasis, which could guide clinical application.

Visualizing the phage T4 activated transcription complex of DNA and E. coli RNA polymerase

The ability of RNA polymerase (RNAP) to select the right promoter sequence at the right time is fundamental to the control of gene expression in all organisms. However, there is only one crystallized structure of a complete activator/RNAP/DNA complex. In a process called σ appropriation, bacteriophage T4 activates a class of phage promoters using an activator (MotA) and a co-activator (AsiA), which function through interactions with the σ⁷⁰ subunit of RNAP. We have developed a holistic, structure-based model for σ appropriation using multiple experimentally determined 3D structures (Escherichia coli RNAP, the Thermus aquaticus RNAP/DNA complex, AsiA /σ⁷⁰ Region 4, the N-terminal domain of MotA [MotA^NTD], and the C-terminal domain of MotA [MotA^CTD]), molecular modeling, and extensive biochemical observations indicating the position of the proteins relative to each other and to the DNA. Our results visualize how AsiA/MotA redirects σ, and therefore RNAP activity, to T4 promoter DNA, and demonstrate at a molecular level how the tactful interaction of transcriptional factors with even small segments of RNAP can alter promoter specificity. Furthermore, our model provides a rational basis for understanding how a mutation within the β subunit of RNAP (G1249D), which is far removed from AsiA or MotA, impairs σ appropriation.

Computational studies of human class V alcohol dehydrogenase - the odd sibling

Background

All known attempts to isolate and characterize mammalian class V alcohol dehydrogenase (class V ADH), a member of the large ADH protein family, at the protein level have failed. This indicates that the class V ADH protein is not stable in a non-cellular environment, which is in contrast to all other human ADH enzymes. In this report we present evidence, supported with results from computational analyses performed in combination with earlier in vitro studies, why this ADH behaves in an atypical way.

Results

Using a combination of structural calculations and sequence analyses, we were able to identify local structural differences between human class V ADH and other human ADHs, including an elongated β-strands and a labile α-helix at the subunit interface region of each chain that probably disturb it. Several amino acid residues are strictly conserved in class I–IV, but altered in class V ADH. This includes a for class V ADH unique and conserved Lys51, a position directly involved in the catalytic mechanism in other ADHs, and nine other class V ADH-specific residues.

Conclusions

In this study we show that there are pronounced structural changes in class V ADH as compared to other ADH enzymes. Furthermore, there is an evolutionary pressure among the mammalian class V ADHs, which for most proteins indicate that they fulfill a physiological function. We assume that class V ADH is expressed, but unable to form active dimers in a non-cellular environment, and is an atypical mammalian ADH. This is compatible with previous experimental characterization and present structural modelling. It can be considered the odd sibling of the ADH protein family and so far seems to be a pseudoenzyme with another hitherto unknown physiological function.

Electronic supplementary material

The online version of this article (doi:10.1186/s12858-016-0072-y) contains supplementary material, which is available to authorized users.

An antiarrhythmic agent as a promising lead compound for targeting the hEAG1 ion channel in cancer therapy: insights from molecular dynamics simulations

Experimental evidence suggests that hERG and hEAG potassium channels may serve as important cancer therapy targets because either of the channel blockade or inactivation by different methods leads to inhibition of cancer cells growth and proliferation. However, there is no known hEAG specific blocker, and hERG blockade leads to adverse cardiac side effects, although it is currently used in treating certain types of arrhythmias. There have been some attempts to explain the channels blockade by clofilium, an antiarrhythmic agent, and the results lead to different possible binding modes. This study investigates for the first time the potential of using clofilium as a lead compound for finding a novel cancer therapy agent which may target ion channels. The implied findings from a comparative assessment of literature studies were verified using molecular dynamics simulations. The results indicate a particular structural difference between the two channels that could provide a novel and realistic way of using clofilium analogs which may target the hEAG1 ion channel in cancer therapy.

PR Domain-containing Protein 7 (PRDM7) Is a Histone 3 Lysine 4 Trimethyltransferase

PR domain-containing protein 7 (PRDM7) is a primate-specific histone methyltransferase that is the result of a recent gene duplication of PRDM9. The two proteins are highly homologous, especially in the catalytic PR/SET domain, where they differ by only three amino acid residues. Here we report that PRDM7 is an efficient methyltransferase that selectively catalyzes the trimethylation of H3 lysine 4 (H3K4) both in vitro and in cells. Through selective mutagenesis we have dissected the functional roles of each of the three divergent residues between the PR domains of PRDM7 and PRDM9. These studies indicate that after a single serine to tyrosine mutation at residue 357 (S357Y), PRDM7 regains the substrate specificities and catalytic activities similar to its evolutionary predecessor, including the ability to efficiently methylate H3K36.

Functional Complementation Assay for 47 MUTYH Variants in a MutY-Disrupted Escherichia coli Strain

MUTYH-associated polyposis (MAP) is an adenomatous polyposis transmitted in an autosomal-recessive pattern, involving biallelic inactivation of the MUTYH gene. Loss of a functional MUTYH protein will result in the accumulation of G:T mismatched DNA caused by oxidative damage. Although p.Y179C and p.G396D are the two most prevalent MUTYH variants, more than 200 missense variants have been detected. It is difficult to determine whether these variants are disease-causing mutations or single-nucleotide polymorphisms. To understand the functional consequences of these variants, we generated 47 MUTYH gene variants via site-directed mutagenesis, expressed the encoded proteins in MutY-disrupted Escherichia coli, and assessed their abilities to complement the functional deficiency in the E. coli by monitoring spontaneous mutation rates. Although the majority of variants exhibited intermediate complementation relative to the wild type, some variants severely interfered with this complementation. However, some variants retained functioning similar to the wild type. In silico predictions of functional effects demonstrated a good correlation. Structural prediction of MUTYH based on the MutY protein structure allowed us to interpret effects on the protein stability or catalytic activity. These data will be useful for evaluating the functional consequences of missense MUTYH variants detected in patients with suspected MAP.

Mutations in Known Monogenic High Bone Mass Loci Only Explain a Small Proportion of High Bone Mass Cases

High bone mass (HBM) can be an incidental clinical finding; however, monogenic HBM disorders (eg, LRP5 or SOST mutations) are rare. We aimed to determine to what extent HBM is explained by mutations in known HBM genes. A total of 258 unrelated HBM cases were identified from a review of 335,115 DXA scans from 13 UK centers. Cases were assessed clinically and underwent sequencing of known anabolic HBM loci: LRP5 (exons 2, 3, 4), LRP4 (exons 25, 26), SOST (exons 1, 2, and the van Buchem's disease [VBD] 52-kb intronic deletion 3'). Family members were assessed for HBM segregation with identified variants. Three-dimensional protein models were constructed for identified variants. Two novel missense LRP5 HBM mutations ([c.518C>T; p.Thr173Met], [c.796C>T; p.Arg266Cys]) were identified, plus three previously reported missense LRP5 mutations ([c.593A>G; p.Asn198Ser], [c.724G>A; p.Ala242Thr], [c.266A>G; p.Gln89Arg]), associated with HBM in 11 adults from seven families. Individuals with LRP5 HBM (∼prevalence 5/100,000) displayed a variable phenotype of skeletal dysplasia with increased trabecular BMD and cortical thickness on HRpQCT, and gynoid fat mass accumulation on DXA, compared with both non-LRP5 HBM and controls. One mostly asymptomatic woman carried a novel heterozygous nonsense SOST mutation (c.530C>A; p.Ser177X) predicted to prematurely truncate sclerostin. Protein modeling suggests the severity of the LRP5-HBM phenotype corresponds to the degree of protein disruption and the consequent effect on SOST-LRP5 binding. We predict p.Asn198Ser and p.Ala242Thr directly disrupt SOST binding; both correspond to severe HBM phenotypes (BMD Z-scores +3.1 to +12.2, inability to float). Less disruptive structural alterations predicted from p.Arg266Cys, p.Thr173Met, and p.Gln89Arg were associated with less severe phenotypes (Z-scores +2.4 to +6.2, ability to float). In conclusion, although mutations in known HBM loci may be asymptomatic, they only account for a very small proportion (∼3%) of HBM individuals, suggesting the great majority are explained by either unknown monogenic causes or polygenic inheritance.

P16INK4a Upregulation Mediated by SIX6 Defines Retinal Ganglion Cell Pathogenesis in Glaucoma

Glaucoma, a blinding neurodegenerative disease, whose risk factors include elevated intraocular pressure (IOP), age, and genetics, is characterized by accelerated and progressive retinal ganglion cell (RGC) death. Despite decades of research, the mechanism of RGC death in glaucoma is still unknown. Here, we demonstrate that the genetic effect of the SIX6 risk variant (rs33912345, His141Asn) is enhanced by another major POAG risk gene, p16INK4a (cyclin-dependent kinase inhibitor 2A, isoform INK4a). We further show that the upregulation of homozygous SIX6 risk alleles (CC) leads to an increase in p16INK4a expression, with subsequent cellular senescence, as evidenced in a mouse model of elevated IOP and in human POAG eyes. Our data indicate that SIX6 and/or IOP promotes POAG by directly increasing p16INK4a expression, leading to RGC senescence in adult human retinas. Our study provides important insights linking genetic susceptibility to the underlying mechanism of RGC death and provides a unified theory of glaucoma pathogenesis.

Lanosterol reverses protein aggregation in cataracts

The human lens is comprised largely of crystallin proteins assembled into a highly ordered, interactive macro-structure essential for lens transparency and refractive index. Any disruption of intra- or inter-protein interactions will alter this delicate structure, exposing hydrophobic surfaces, with consequent protein aggregation and cataract formation. Cataracts are the most common cause of blindness worldwide, affecting tens of millions of people, and currently the only treatment is surgical removal of cataractous lenses. The precise mechanisms by which lens proteins both prevent aggregation and maintain lens transparency are largely unknown. Lanosterol is an amphipathic molecule enriched in the lens. It is synthesized by lanosterol synthase (LSS) in a key cyclization reaction of a cholesterol synthesis pathway. Here we identify two distinct homozygous LSS missense mutations (W581R and G588S) in two families with extensive congenital cataracts. Both of these mutations affect highly conserved amino acid residues and impair key catalytic functions of LSS. Engineered expression of wild-type, but not mutant, LSS prevents intracellular protein aggregation of various cataract-causing mutant crystallins. Treatment by lanosterol, but not cholesterol, significantly decreased preformed protein aggregates both in vitro and in cell-transfection experiments. We further show that lanosterol treatment could reduce cataract severity and increase transparency in dissected rabbit cataractous lenses in vitro and cataract severity in vivo in dogs. Our study identifies lanosterol as a key molecule in the prevention of lens protein aggregation and points to a novel strategy for cataract prevention and treatment.

All-Atom Internal Coordinate Mechanics (ICM) Force Field for Hexopyranoses and Glycoproteins

We present an extension of the all-atom internal-coordinate force field, ICMFF, that allows for simulation of heterogeneous systems including hexopyranose saccharides and glycan chains in addition to proteins. A library of standard glycan geometries containing α- and β-anomers of the most common hexapyranoses, i.e., d-galactose, d-glucose, d-mannose, d-xylose, l-fucose, N-acetylglucosamine, N-acetylgalactosamine, sialic, and glucuronic acids, is created based on the analysis of the saccharide structures reported in the Cambridge Structural Database. The new force field parameters include molecular electrostatic potential-derived partial atomic charges and the torsional parameters derived from quantum mechanical data for a collection of minimal molecular fragments and related molecules. The ϕ/ψ torsional parameters for different types of glycosidic linkages are developed using model compounds containing the key atoms in the full carbohydrates, i.e., glycosidic-linked tetrahydropyran-cyclohexane dimers. Target data for parameter optimization include two-dimensional energy surfaces corresponding to the ϕ/ψ glycosidic dihedral angles in the disaccharide analogues, as determined by quantum mechanical MP2/6-31G** single-point energies on HF/6-31G** optimized structures. To achieve better agreement with the observed geometries of glycosidic linkages, the bond angles at the O-linkage atoms are added to the internal variable set and the corresponding bond bending energy term is parametrized using quantum mechanical data. The resulting force field is validated on glycan chains of 1-12 residues from a set of high-resolution X-ray glycoprotein structures based on heavy atom root-mean-square deviations of the lowest-energy glycan conformations generated by the biased probability Monte Carlo (BPMC) molecular mechanics simulations from the native structures. The appropriate BPMC distributions for monosaccharide-monosaccharide and protein-glycan linkages are derived from the extensive analysis of conformational properties of glycoprotein structures reported in the Protein Data Bank. Use of the BPMC search leads to significant improvements in sampling efficiency for glycan simulations. Moreover, good agreement with the X-ray glycoprotein structures is achieved for all glycan chain lengths. Thus, average/median RMSDs are 0.81/0.68 Å for one-residue glycans and 1.32/1.47 Å for three-residue glycans. RMSD from the native structure for the lowest-energy conformation of the 12-residue glycan chain (PDB ID 3og2) is 1.53 Å. Additionally, results obtained for free short oligosaccharides using the new force field are in line with the available experimental data, i.e., the most populated conformations in solution are predicted to be the lowest energy ones. The newly developed parameters allow for the accurate modeling of linear and branched hexopyranose glycosides in heterogeneous systems.

Modulation of Macrophage Gene Expression via Liver X Receptor α Serine 198 Phosphorylation

In mouse models of atherosclerosis, normalization of hyperlipidemia promotes macrophage emigration and regression of atherosclerotic plaques in part by liver X receptor (LXR)-mediated induction of the chemokine receptor CCR7. Here we report that LXRα serine 198 (S198) phosphorylation modulates CCR7 expression. Low levels of S198 phosphorylation are observed in plaque macrophages in the regression environment where high levels of CCR7 expression are observed. Consistent with these findings, CCR7 gene expression in human and mouse macrophages cell lines is induced when LXRα at S198 is nonphosphorylated. In bone marrow-derived macrophages (BMDMs), we also observed induction of CCR7 by ligands that promote nonphosphorylated LXRα S198, and this was lost in LXR-deficient BMDMs. LXRα occupancy at the CCR7 promoter is enhanced and histone modifications associated with gene repression are reduced in RAW264.7 cells expressing nonphosphorylated LXRα (RAW-LXRα S198A) compared to RAW264.7 cells expressing wild-type (WT) phosphorylated LXRα (RAW-LXRα WT). Expression profiling of ligand-treated RAW-LXRα S198A cells compared to RAW-LXRα WT cells revealed induction of cell migratory and anti-inflammatory genes and repression of proinflammatory genes. Modeling of LXRα S198 in the nonphosphorylated and phosphorylated states identified phosphorylation-dependent conformational changes in the hinge region commensurate with the presence of sites for protein interaction. Therefore, gene transcription is regulated by LXRα S198 phosphorylation, including that of antiatherogenic genes such as CCR7.

Exploring a 2-naphthoic acid template for the structure-based design of P2Y14 receptor antagonist molecular probes

The P2Y₁₄ receptor (P2Y₁₄R), one of eight P2Y G protein-coupled receptors (GPCR), is involved in inflammatory, endocrine, and hypoxic processes and is an attractive pharmaceutical target. The goal of this research is to develop high-affinity P2Y₁₄R fluorescent probes based on the potent and highly selective antagonist 4-(4-(piperidin-4-yl)-phenyl)-7-(4-(trifluoromethyl)-phenyl)-2-naphthoic acid (6, PPTN). A model of hP2Y₁₄R based on recent hP2Y₁₂R X-ray structures together with simulated antagonist docking suggested that the piperidine ring is suitable for fluorophore conjugation while preserving affinity. Chain-elongated alkynyl or amino derivatives of 6 for click or amide coupling were synthesized, and their antagonist activities were measured in hP2Y₁₄R-expressing CHO cells. Moreover, a new Alexa Fluor 488 (AF488) containing derivative 30 (MRS4174, K_i = 80 pM) exhibited exceptionally high affinity, as compared to 13 nM for the alkyne precursor 22. A flow cytometry assay employing 30 as a fluorescent probe was used to quantify specific binding to P2Y₁₄R. Known P2Y receptor ligands inhibited binding of 30 with properties consistent with their previously established receptor selectivities and affinities. These results illustrate that potency in this series of 2-naphthoic acid derivatives can be preserved by chain functionalization, leading to highly potent fluorescent molecular probes for P2Y₁₄R. Such conjugates will be useful tools in expanding the SAR of this receptor, which still lacks chemical diversity in its collective ligands. This approach demonstrates the predictive power of GPCR homology modeling and the relevance of newly determined X-ray structures to GPCR medicinal chemistry.

Studies of synthetic chalcone derivatives as potential inhibitors of secretory phospholipase A2, cyclooxygenases, lipoxygenase and pro-inflammatory cytokines

Arachidonic acid metabolism leads to the generation of key lipid mediators which play a fundamental role during inflammation. The inhibition of enzymes involved in arachidonic acid metabolism has been considered as a synergistic anti-inflammatory effect with enhanced spectrum of activity. A series of 1,3-diphenyl-2-propen-1-one derivatives were investigated for anti-inflammatory related activities involving inhibition of secretory phospholipase A2, cyclooxygenases, soybean lipoxygenase, and lipopolysaccharides-induced secretion of interleukin-6 and tumor necrosis factor-alpha in mouse RAW264.7 macrophages. The results from the above mentioned assays exhibited that the synthesized compounds were effective inhibitors of pro-inflammatory enzymes and cytokines. The results also revealed that the chalcone derivatives with 4-methlyamino ethanol substitution seem to be significant for inhibition of enzymes and cytokines. Molecular docking experiments were carried out to elucidate the molecular aspects of the observed inhibitory activities of the investigated compounds. Present findings increase the possibility that these chalcone derivatives might serve as a beneficial starting point for the design and development of improved anti-inflammatory agents.

Structural basis for activation of trimeric Gi proteins by multiple growth factor receptors via GIV/Girdin

GIV, a guanidine exchange factor for trimeric Gi, contains a unique domain that functions like a SH2 domain. GIV's SH2-like domain binds autophosphorylated RTKs. Binding of GIV's SH2 to RTKs enables the receptors to activate trimeric Gi. Inhibition of GIV:RTK interaction abolishes GIV-dependent Akt enhancement downstream of RTKs.

A long-standing issue in the field of signal transduction is to understand the cross-talk between receptor tyrosine kinases (RTKs) and heterotrimeric G proteins, two major and distinct signaling hubs that control eukaryotic cell behavior. Although stimulation of many RTKs leads to activation of trimeric G proteins, the molecular mechanisms behind this phenomenon remain elusive. We discovered a unifying mechanism that allows GIV/Girdin, a bona fide metastasis-related protein and a guanine-nucleotide exchange factor (GEF) for Gαi, to serve as a direct platform for multiple RTKs to activate Gαi proteins. Using a combination of homology modeling, protein–protein interaction, and kinase assays, we demonstrate that a stretch of ∼110 amino acids within GIV C-terminus displays structural plasticity that allows folding into a SH2-like domain in the presence of phosphotyrosine ligands. Using protein–protein interaction assays, we demonstrated that both SH2 and GEF domains of GIV are required for the formation of a ligand-activated ternary complex between GIV, Gαi, and growth factor receptors and for activation of Gαi after growth factor stimulation. Expression of a SH2-deficient GIV mutant (Arg 1745→Leu) that cannot bind RTKs impaired all previously demonstrated functions of GIV—Akt enhancement, actin remodeling, and cell migration. The mechanistic and structural insights gained here shed light on the long-standing questions surrounding RTK/G protein cross-talk, set a novel paradigm, and characterize a unique pharmacological target for uncoupling GIV-dependent signaling downstream of multiple oncogenic RTKs.

Characterization of human variants in obesity-related SIM1 protein identifies a hot-spot for dimerization with the partner protein ARNT2

Several non-synonymous variants found in obese patients disrupt function of the transcription factor SIM1 (single-minded 1) by impairing binding to an essential partner protein. The clustering of these variants reveals a mutational hot-spot critical for function of SIM1 and the related protein SIM2.

The second extracellular loop of the adenosine A1 receptor mediates activity of allosteric enhancers

Allosteric enhancers of the adenosine A1 receptor amplify signaling by orthosteric agonists. Allosteric enhancers are appealing drug candidates because their activity requires that the orthosteric site be occupied by an agonist, thereby conferring specificity to stressed or injured tissues that produce adenosine. To explore the mechanism of allosteric enhancer activity, we examined their action on several A1 receptor constructs, including (1) species variants, (2) species chimeras, (3) alanine scanning mutants, and (4) site-specific mutants. These findings were combined with homology modeling of the A1 receptor and in silico screening of an allosteric enhancer library. The binding modes of known docked allosteric enhancers correlated with the known structure-activity relationship, suggesting that these allosteric enhancers bind to a pocket formed by the second extracellular loop, flanked by residues S150 and M162. We propose a model in which this vestibule controls the entry and efflux of agonists from the orthosteric site and agonist binding elicits a conformational change that enables allosteric enhancer binding. This model provides a mechanism for the observations that allosteric enhancers slow the dissociation of orthosteric agonists but not antagonists.

A yeast screening method to decipher the interaction between the adenosine A2B receptor and the C-terminus of different G protein α-subunits

The expression of human G protein-coupled receptors (GPCRs) in Saccharomyces cerevisiae containing chimeric yeast/mammalian Gα subunits provides a useful tool for the study of GPCR activation. In this study, we used a one-GPCR-one-G protein yeast screening method in combination with molecular modeling and mutagenesis studies to decipher the interaction between GPCRs and the C-terminus of different α-subunits of G proteins. We chose the human adenosine A2B receptor (hA2BR) as a paradigm, a typical class A GPCR that shows promiscuous behavior in G protein coupling in this yeast system. The wild-type hA2BR and five mutant receptors were expressed in 8 yeast strains with different humanized G proteins, covering the four major classes: Gαi, Gαs, Gαq, and Gα12. Our experiments showed that a tyrosine residue (Y) at the C-terminus of the Gα subunit plays an important role in controlling the activation of GPCRs. Receptor residues R103(3.50) and I107(3.54) are vital too in G protein-coupling and the activation of the hA2BR, whereas L213(IL3) is more important in G protein inactivation. Substitution of S235(6.36) to alanine provided the most divergent G protein-coupling profile. Finally, L236(6.37) substitution decreased receptor activation in all G protein pathways, although to a different extent. In conclusion, our findings shed light on the selectivity of receptor/G protein coupling, which may help in further understanding GPCR signaling.