Comparative Protein Structure Modeling Using MODELLER

A divergent protein kinase A regulatory subunit essential for morphogenesis of the human pathogen Leishmania

Parasitic protozoa of the genus Leishmania cycle between the phagolysosome of mammalian macrophages, where they reside as rounded intracellular amastigotes, and the midgut of female sand flies, which they colonize as elongated extracellular promastigotes. Previous studies indicated that protein kinase A (PKA) plays an important role in the initial steps of promastigote differentiation into amastigotes. Here, we describe a novel regulatory subunit of PKA (which we have named PKAR3) that is unique to Leishmania and most (but not all) other Kinetoplastidae. PKAR3 is localized to subpellicular microtubules (SPMT) in the cell cortex, where it recruits a specific catalytic subunit (PKAC3). Promastigotes of pkar3 or pkac3 null mutants lose their elongated shape and become rounded but remain flagellated. Truncation of an N-terminal formin homology (FH)-like domain of PKAR3 results in its detachment from the SPMT, also leading to rounded promastigotes. Thus, the tethering of PKAC3 via PKAR3 at the cell cortex is essential for maintenance of the elongated shape of promastigotes. This role of PKAR3 is reminiscent of PKARIβ and PKARIIβ binding to microtubules of mammalian neurons, which is essential for the elongation of dendrites and axons, respectively. Interestingly, PKAR3 binds nucleoside analogs, but not cAMP, with a high affinity similar to the PKAR1 isoform of Trypanosoma. We propose that these early-diverged protists have re-purposed PKA for a novel signaling pathway that spatiotemporally controls microtubule remodeling and cell shape.

4″-Alkyl EGCG Derivatives Induce Cytoprotective Autophagy Response by Inhibiting EGFR in Glioblastoma Cells

Epidermal growth factor receptor (EGFR)-targeted therapy has been proven vital in the last two decades for the treatment of multiple cancer types, including nonsmall cell lung cancer, glioblastoma, breast cancer and head and neck squamous cell carcinoma. Unfortunately, the majority of approved EGFR inhibitors fall into the drug resistance category because of continuous mutations and acquired resistance. Recently, autophagy has surfaced as one of the emerging underlying mechanisms behind resistance to EGFR-tyrosine kinase inhibitors (TKIs). Previously, we developed a series of 4″-alkyl EGCG (4″-Cn EGCG, n = 6, 8, 10, 12, 14, 16, and 18) derivatives with enhanced anticancer effects and stability. Therefore, the current study hypothesized that 4″-alkyl EGCG might induce cytoprotective autophagy upon EGFR inhibition, and inhibition of autophagy may lead to improved cytotoxicity. In this study, we have observed growth inhibition and caspase-3-dependent apoptosis in 4″-alkyl EGCG derivative-treated glioblastoma cells (U87-MG). We also confirmed that 4″-alkyl EGCG could inhibit EGFR in the cells, as well as mutant L858R/T790M EGFR, through an in vitro kinase assay. Furthermore, we have found that EGFR inhibition with 4″-alkyl EGCG induces cytoprotective autophagic responses, accompanied by the blockage of the AKT/mTOR signaling pathway. In addition, cytotoxicity caused by 4″-C10 EGCG, 4″-C12 EGCG, and 4″-C14 EGCG was significantly increased after the inhibition of autophagy by the pharmacological inhibitor chloroquine. These findings enhance our understanding of the autophagic response toward EGFR inhibitors in glioblastoma cells and suggest a potent combinatorial strategy to increase the therapeutic effectiveness of EGFR-TKIs.

Disease-associated nonsense and frame-shift variants resulting in the truncation of the GluN2A or GluN2B C-terminal domain decrease NMDAR surface expression and reduce potentiating effects of neurosteroids

N-methyl-D-aspartate receptors (NMDARs) play a critical role in normal brain function, and variants in genes encoding NMDAR subunits have been described in individuals with various neuropsychiatric disorders. We have used whole-cell patch-clamp electrophysiology, fluorescence microscopy and in-silico modeling to explore the functional consequences of disease-associated nonsense and frame-shift variants resulting in the truncation of GluN2A or GluN2B C-terminal domain (CTD). This study characterizes variant NMDARs and shows their reduced surface expression and synaptic localization, altered agonist affinity, increased desensitization, and reduced probability of channel opening. We also show that naturally occurring and synthetic steroids pregnenolone sulfate and epipregnanolone butanoic acid, respectively, enhance NMDAR function in a way that is dependent on the length of the truncated CTD and, further, is steroid-specific, GluN2A/B subunit-specific, and GluN1 splice variant-specific. Adding to the previously described effects of disease-associated NMDAR variants on the receptor biogenesis and function, our results improve the understanding of the molecular consequences of NMDAR CTD truncations and provide an opportunity for the development of new therapeutic neurosteroid-based ligands.

Identification of allosteric inhibitor against AKT1 through structure-based virtual screening

AKT (serine/threonine protein kinase) is a potential therapeutic target for many types of cancer as it plays a vital role in cancer progression. Many AKT inhibitors are already in practice under single and combinatorial therapy. However, most of these inhibitors are orthosteric / pan-AKT that are non-selective and non-specific to AKT kinase and their isoforms. Hence, researchers are searching for novel allosteric inhibitors that bind in the interface between pH and kinase domain. In this study, we performed structure-based virtual screening from the afroDB (a diverse natural compounds library) to find the potential inhibitor targeting the AKT1. These compounds were filtered through Lipinski, ADMET properties, combined with a molecular docking approach to obtain the 8 best compounds. Then we performed molecular dynamics simulation for apoprotein, AKT1 with 8 complexes, and AKT1 with the positive control (Miransertib). Molecular docking and simulation analysis revealed that Bianthracene III (hit 1), 10-acetonyl Knipholonecyclooxanthrone (hit 2), Abyssinoflavanone VII (hit 5) and 8-c-p-hydroxybenzyldiosmetin (hit 6) had a better binding affinity, stability, and compactness than the reference compound. Notably, hit 1, hit 2 and hit 5 had molecular features required for allosteric inhibition.

Inhibitors of 3-Hydroxy-3-methylglutaryl Coenzyme A Reductase Decrease the Growth, Ergosterol Synthesis and Generation of petite Mutants in Candida glabrata and Candida albicans

Candida glabrata and Candida albicans, the most frequently isolated candidiasis species in the world, have developed mechanisms of resistance to treatment with azoles. Among the clinically used antifungal drugs are statins and other compounds that inhibit 3-hydroxy-3-methylglutaryl coenzyme A reductase (HMGR), resulting in decreased growth and ergosterol levels in yeasts. Ergosterol is a key element for the formation of the yeast cell membrane. However, statins often cause DNA damage to yeast cells, facilitating mutation and drug resistance. The aim of the current contribution was to synthesize seven series of compounds as inhibitors of the HMGR enzyme of Candida ssp., and to evaluate their effect on cellular growth, ergosterol synthesis and generation of petite mutants of C. glabrata and C. albicans. Compared to the reference drugs (fluconazole and simvastatin), some HMGR inhibitors caused lower growth and ergosterol synthesis in the yeast species and generated fewer petite mutants. Moreover, heterologous expression was achieved in Pichia pastoris, and compounds 1a, 1b, 6g and 7a inhibited the activity of recombinant CgHMGR and showed better binding energy values than for α-asarone and simvastatin. Thus, we believe these are good candidates for future antifungal drug development.

Molecular dynamics simulation study on Bacillus subtilis EngA: the presence of Mg2+ at the active-sites promotes the functionally important conformation

EngA, a GTPase contains two GTP binding domains [GD1, GD2], and the C-terminal KH domain shown to be involved in the later stages of ribosome maturation. Association of EngA to the ribosomal subunit in the intermediate stage of maturation is essential for complete ribosome maturation. However, this association was shown to be dependent on the nucleotide bound combinations. This nucleotide dependent association tendency is attributed to the conformational changes that occur among different nucleotide bound combinations. Therefore, to explore the conformational changes, all-atom molecular dynamics simulations for Bacillus subtilis EngA in different nucleotide bound combinations along with the presence or absence of Mg2+ in the active-sites were carried out. The presence of Mg2+ along with the bound nucleotide at the GD2 active-site dictates the GD2-Sw-II mobility, but the GD1-Sw-II mobility has not shown any nucleotide or Mg2+ dependent movement. However, the GD1-Sw-II secondary conformations are shown to be influenced by the GD2 nucleotide bound state. This allosteric connection between the GD2 active-site and the GD1-Sw-II is also observed through the dynamic network analysis. Further, the exploration of the GD1-KH interface interactions exhibited a more attractive tendency when GD1 is bound to GTP-Mg2+. In addition, the presence of Mg2+ stabilizes active-site water and also increases the distances between the α- and γ- phosphates of the bound GTP. Curiously, three water molecules in the GD1 active-site and only one water molecule in the GD2 active-site are stabilized. This indicates that the probability of GTP hydrolysis is more in GD1 compared to GD2.Communicated by Ramaswamy H. Sarma.

The solution structure of the unbound IgG Fc receptor CD64 resembles its crystal structure: Implications for function

FcγRI (CD64) is the only high-affinity Fcγ receptor found on monocytes, macrophages, eosinophils, neutrophils and dendritic cells. It binds immunoglobulin G (IgG) antibody-antigen complexes at its Fc region to trigger key immune responses. CD64 contains three immunoglobulin-fold extracellular domains (D1, D2 and D3) and a membrane-spanning region. Despite the importance of CD64, no solution structure for this is known to date. To investigate this, we used analytical ultracentrifugation, small-angle X-ray scattering, and atomistic modelling. Analytical ultracentrifugation revealed that CD64 was monomeric with a sedimentation coefficient s020,w of 2.53 S, together with some dimer. Small-angle X-ray scattering showed that its radius of gyration RG was 3.3-3.4 nm and increased at higher concentrations to indicate low dimerization. Monte Carlo modelling implemented in the SASSIE-web package generated 279,162 physically-realistic trial CD64 structures. From these, the scattering best-fit models at the lowest measured concentrations that minimised dimers revealed that the D1, D2 and D3 domains were structurally similar to those seen in three CD64 crystal structures, but showed previously unreported flexibility between D1, D2 and D3. Despite the limitations of the scattering data, the superimposition of the CD64 solution structures onto crystal structures of the IgG Fc-CD64 complex showed that the CD64 domains do not sterically clash with the IgG Fc region, i.e. the solution structure of CD64 was sufficiently compact to allow IgG to bind to its high-affinity Fcγ receptor. This improved understanding may result in novel approaches to inhibit CD64 function, and opens the way for the solution study of the full-length CD64-IgG complex.

Crystal structure and activity of a de novo enzyme, ferric enterobactin esterase Syn-F4

Producing novel enzymes that are catalytically active in vitro and biologically functional in vivo is a key goal of synthetic biology. Previously, we reported Syn-F4, the first de novo protein that meets both criteria. Syn-F4 hydrolyzed the siderophore ferric enterobactin, and expression of Syn-F4 allowed an inviable strain of Escherichia coli (Δfes) to grow in iron-limited medium. Here, we describe the crystal structure of Syn-F4. Syn-F4 forms a dimeric 4-helix bundle. Each monomer comprises two long α-helices, and the loops of the Syn-F4 dimer are on the same end of the bundle (syn topology). Interestingly, there is a penetrated hole in the central region of the Syn-F4 structure. Extensive mutagenesis experiments in a previous study showed that five residues (Glu26, His74, Arg77, Lys78, and Arg85) were essential for enzymatic activity in vivo. All these residues are located around the hole in the central region of the Syn-F4 structure, suggesting a putative active site with a catalytic dyad (Glu26-His74). The complete inactivity of purified proteins with mutations at the five residues supports the putative active site and reaction mechanism. Molecular dynamics and docking simulations of the ferric enterobactin siderophore binding to the Syn-F4 structure demonstrate the dynamic property of the putative active site. The structure and active site of Syn-F4 are completely different from native enterobactin esterase enzymes, thereby demonstrating that proteins designed de novo can provide life-sustaining catalytic activities using structures and mechanisms dramatically different from those that arose in nature.

FA Sliding as the Mechanism for the ANT1-Mediated Fatty Acid Anion Transport in Lipid Bilayers

Mitochondrial adenine nucleotide translocase (ANT) exchanges ADP for ATP to maintain energy production in the cell. Its protonophoric function in the presence of long-chain fatty acids (FA) is also recognized. Our previous results imply that proton/FA transport can be best described with the FA cycling model, in which protonated FA transports the proton to the mitochondrial matrix. The mechanism by which ANT1 transports FA anions back to the intermembrane space remains unclear. Using a combined approach involving measurements of the current through the planar lipid bilayers reconstituted with ANT1, site-directed mutagenesis and molecular dynamics simulations, we show that the FA anion is first attracted by positively charged arginines or lysines on the matrix side of ANT1 before moving along the positively charged protein-lipid interface and binding to R79, where it is protonated. We show that R79 is also critical for the competitive binding of ANT1 substrates (ADP and ATP) and inhibitors (carboxyatractyloside and bongkrekic acid). The binding sites are well conserved in mitochondrial SLC25 members, suggesting a general mechanism for transporting FA anions across the inner mitochondrial membrane.

Host trehalose metabolism disruption by validamycin A results in reduced fitness of parasitoid offspring

The general cutworm, Spodoptera litura (Lepidoptera: Noctuidae) is a worldwide destructive omnivorous pest and the endoparasitoid wasp Meteorus pulchricornis (Hymenoptera: Braconidae) is the dominant endoparasitoid of S. litura larvae. Trehalase is a key enzyme in insect trehalose metabolism and plays an important role in the growth and development of insects. However, the specific function of trehalase in parasitoid and host associations has been less reported. In this study, we obtained two trehalase genes (SlTre1 and SlTre2) from our previously constructed S. litura transcriptome database; they were highly expressed in 3rd instar larvae. SlTre1 was mainly expressed in the midgut, and SlTre2 was expressed highest in the head. SlTre1 and SlTre2 were highly expressed 5 days after parasitization by M. pulchricornis. Treatment with the trehalase inhibitor validamycin A significantly inhibited the expression levels of SlTre1 and SlTre2, and the trehalase activity. Besides, the content of trehalose was increased but the content of glucose was decreased 24 h after validamycin A treatment in parasitized S. litura larvae. In addition, the immune-related genes in phenoloxidase (PO) pathway and fatty acid synthesis-related genes in lipid metabolism were upregulated in parasitized host larvae after validamycin A treatment. Importantly, the emergence rate, proportion of normal adults, and body size of parasitoid offspring was decreased in parasitized S. litura larvae after validamycin A treatment, indicating that validamycin A disrupts the trehalose metabolism of parasitized host and thus reduces the fitness of parasitoid offspring. The present study provides a novel perspective for coordinating the application of biocontrol and antibiotics in agroecosystem.

Perspectives of Proteomics in Respiratory Allergic Diseases

Proteomics in respiratory allergic diseases has such a battery of techniques and programs that one would almost think there is nothing impossible to find, invent or mold. All the resources that we document here are involved in solving problems in allergic diseases, both diagnostic and prognostic treatment, and immunotherapy development. The main perspectives, according to this version, are in three strands and/or a lockout immunological system: (1) Blocking the diapedesis of the cells involved, (2) Modifications and blocking of paratopes and epitopes being understood by modifications to antibodies, antagonisms, or blocking them, and (3) Blocking FcεRI high-affinity receptors to prevent specific IgEs from sticking to mast cells and basophils. These tools and targets in the allergic landscape are, in our view, the prospects in the field. However, there are still many allergens to identify, including some homologies between allergens and cross-reactions, through the identification of structures and epitopes. The current vision of using proteomics for this purpose remains a constant; this is also true for the basis of diagnostic and controlled systems for immunotherapy. Ours is an open proposal to use this vision for treatment.

Membrane mediated mechanical stimuli produces distinct active-like states in the AT1 receptor

The Angiotensin II Type 1 (AT1) receptor is one of the most widely studied GPCRs within the context of biased signaling. While the AT1 receptor is activated by agonists such as the peptide AngII, it can also be activated by mechanical stimuli such as membrane stretch or shear in the absence of a ligand. Despite the importance of mechanical activation of the AT1 receptor in biological processes such as vasoconstriction, little is known about the structural changes induced by external physical stimuli mediated by the surrounding lipid membrane. Here, we present a systematic simulation study that characterizes the activation of the AT1 receptor under various membrane environments and mechanical stimuli. We show that stability of the active state is highly sensitive to membrane thickness and tension. Structural comparison of membrane-mediated vs. agonist-induced activation shows that the AT1 receptor has distinct active conformations. This is supported by multi-microsecond free energy calculations that show unique landscapes for the inactive and various active states. Our modeling results provide structural insights into the mechanical activation of the AT1 receptor and how it may produce different functional outcomes within the framework of biased agonism.

DNA facilitates heterodimerization between human transcription factors FoxP1 and FoxP2 by increasing their conformational flexibility

Transcription factors regulate gene expression by binding to DNA. They have disordered regions and specific DNA-binding domains. Binding to DNA causes structural changes, including folding and interactions with other molecules. The FoxP subfamily of transcription factors in humans is unique because they can form heterotypic interactions without DNA. However, it is unclear how they form heterodimers and how DNA binding affects their function. We used computational and experimental methods to study the structural changes in FoxP1's DNA-binding domain when it forms a heterodimer with FoxP2. We found that FoxP1 has complex and diverse conformational dynamics, transitioning between compact and extended states. Surprisingly, DNA binding increases the flexibility of FoxP1, contrary to the typical folding-upon-binding mechanism. In addition, we observed a 3-fold increase in the rate of heterodimerization after FoxP1 binds to DNA. These findings emphasize the importance of structural flexibility in promoting heterodimerization to form transcriptional complexes.

Association of a single amino acid replacement with dorsal pigmentation in a lizard from the Qinghai-Tibetan Plateau

Reptiles can evolve adaptive colors in different environments, but relatively little is known about the genetic mechanisms. Here, we identified the MC1R gene and its association with intraspecific color variation in the lizard Phrynocephalus erythrurus. Analysis of the MC1R sequence in 143 individuals from dark South Qiangtang Plateau (SQP) and light North Qiangtang plateau (NQP) populations, revealed two amino acid sites that showed significant differences in frequency between two areas. One SNP, corresponding to Glu183Lys residue, was found to be a highly significant outlier and differentially fixed for SQP and NQP populations. This residue is located in an extracellular area in the second small extracellular loop within the secondary structure of MC1R, which represents an "attachment pocket" part of the 3D structure. Cytological expression of MC1R alleles with the Glu183Lys replacement showed a 39 % increase in intracellular agonist-induced cyclic AMP levels and a 23.18 % greater cell surface expression of MC1R protein in the SQP relative to the NQP allele. Further in silico 3D modeling and in vitro binding experiments indicated a higher MC1R-α-MSH binding for the SQP allele, and elevated melanin synthesis. We provide an overview of how a single amino acid replacement leads to fundamental changes in MC1R function, and hence shapes variation in dorsal pigmentation in lizards from different environments.

Investigation of the Inhibition Mechanism of Xanthine Oxidoreductase by Oxipurinol: A Computational Study

Xanthine oxidoreductase (XOR) is an enzyme found in various organisms. It converts hypoxanthine to xanthine and urate, which are crucial steps in purine elimination in humans. Elevated uric acid levels can lead to conditions like gout and hyperuricemia. Therefore, there is significant interest in developing drugs that target XOR for treating these conditions and other diseases. Oxipurinol, an analogue of xanthine, is a well-known inhibitor of XOR. Crystallographic studies have revealed that oxipurinol directly binds to the molybdenum cofactor (MoCo) in XOR. However, the precise details of the inhibition mechanism are still unclear, which would be valuable for designing more effective drugs with similar inhibitory functions. In this study, molecular dynamics and quantum mechanics/molecular mechanics calculations are employed to investigate the inhibition mechanism of XOR by oxipurinol. The study examines the structural and dynamic effects of oxipurinol on the pre-catalytic structure of the metabolite-bound system. Our results provide insights on the reaction mechanism catalyzed by the MoCo center in the active site, which aligns well with experimental findings. Furthermore, the results provide insights into the residues surrounding the active site and propose an alternative mechanism for developing alternative covalent inhibitors.

A novel domain-duplicated SlitFAR3 gene involved in sex pheromone biosynthesis in Spodoptera litura

Fatty acyl reductases (FARs) are key enzymes that participate in sex pheromone biosynthesis by reducing fatty acids to fatty alcohols. Lepidoptera typically harbor numerous FAR gene family members. Although FAR genes are involved in the biosynthesis of sex pheromones in moths, the key FAR gene of Spodoptera litura remains unclear. In this work, we predicted 30 FAR genes from the S. litura genome and identified a domain duplication within gene SlitFAR3, which exhibited high and preferential expression in the sexually mature female pheromone glands (PGs) and a rhythmic expression pattern during the scotophase of sex pheromone production. The molecular docking of SlitFAR3, as predicted using a 3D model, revealed a co-factor NADPH binding cavity and 2 substrate binding cavities. Functional expression in yeast cells combined with comprehensive gas chromatography indicated that the SlitFAR3 gene could produce fatty alcohol products. This study is the first to focus on the special phenomenon of FAR domain duplication, which will advance our understanding of biosynthesis-related genes from the perspective of evolutionary biology.

Role of hinge motion and ATP dynamics in factors for inversion stimulation FIS protein deduced while targeting drug resistant Orientia tsutsugamushi

Orientia tsutsugamushi, the causative agent of scrub typhus has been found resistant to various classes of antibiotics such as penicillins, gentamycin and cephalosporins. Review of current literature suggests that the prevalence of scrub typhus has increased globally. Therefore, the current study has aimed at exploring the genome of O. tsutsugamushi to identify potential drug target proteins that can be used for developing novel antibiotics against the pathogen. Subtractive proteomics approach has revealed FIS as a potential drug target protein involved in two component system (TCS), a signaling pathway crucial for bacteria to survive and adjust in changing environmental conditions. Molecular docking studies have revealed compound-356 (CHEMBRIDGE-10040641-3710.356) as a potential inhibitor in both chains A and B of the FIS protein. Simulation results suggest that the docked complex has remained stable and compact throughout the 200 ns run. Significant conformational changes including the hinge motion was observed in the DNA binding domain. Furthermore, the presence of salt bridge between GLU910 and ARG417, rearrangement of interaction residues and displacement of ATP in the central AAA + domain upon binding to the inhibitor were also observed playing a role in stabilizing the protein structure.

Phytopathogenic bacteria utilize host glucose as a signal to stimulate virulence through LuxR homologues

Chemical signal-mediated biological communication is common within bacteria and between bacteria and their hosts. Many plant-associated bacteria respond to unknown plant compounds to regulate bacterial gene expression. However, the nature of the plant compounds that mediate such interkingdom communication and the underlying mechanisms remain poorly characterized. Xanthomonas campestris pv. campestris (Xcc) causes black rot disease on brassica vegetables. Xcc contains an orphan LuxR regulator (XccR) which senses a plant signal that was validated to be glucose by HPLC-MS. The glucose concentration increases in apoplast fluid after Xcc infection, which is caused by the enhanced activity of plant sugar transporters translocating sugar and cell-wall invertases releasing glucose from sucrose. XccR recruits glucose, but not fructose, sucrose, glucose 6-phosphate, and UDP-glucose, to activate pip expression. Deletion of the bacterial glucose transporter gene sglT impaired pathogen virulence and pip expression. Structural prediction showed that the N-terminal domain of XccR forms an alternative pocket neighbouring the AHL-binding pocket for glucose docking. Substitution of three residues affecting structural stability abolished the ability of XccR to bind to the luxXc box in the pip promoter. Several other XccR homologues from plant-associated bacteria can also form stable complexes with glucose, indicating that glucose may function as a common signal molecule for pathogen-plant interactions. The conservation of a glucose/XccR/pip-like system in plant-associated bacteria suggests that some phytopathogens have evolved the ability to utilize host compounds as virulence signals, indicating that LuxRs mediate an interkingdom signalling circuit.

Odorant-Binding Protein 6 Contributes High Binding Affinity to Insecticides in a Parasitic Wasp Meteorus pulchricornis (Hymenoptera: Braconidae)

Meteorus pulchricornis is a preponderant parasitic wasp of various lepidopteran pests. The extensive application of broad-spectrum insecticides usually causes serious threats to the olfactory recognition of nontarget insects such as parasitoid wasps. However, the binding mechanism of odorant-binding proteins (OBPs) to insecticides in parasitoid wasps remains unknown. Herein, we find that the MpulOBP6 protein had a strong binding affinity to three insecticides (phoxim, chlorpyrifos, and chlorfenapyr). Results of computational simulations revealed that the hydrophobic interaction contributed by a mass of nonpolar amino acid residues was the primary driving force in the formation and stabilization of MpulOBP6-insecticide complexes. Among them, four residues (Met75, Val84, Phe121, and Pro122) and two residues (Val84 and Phe111) play an essential role in the binding of MpulOBP6 to phoxim and chlorfenapyr, respectively. Our findings could be instrumental to elucidate the effects of insecticide application toward the olfactory recognition of nontarget insects in the processes of agricultural production.

Protein structure and folding pathway prediction based on remote homologs recognition using PAthreader

Recognition of remote homologous structures is a necessary module in AlphaFold2 and is also essential for the exploration of protein folding pathways. Here, we propose a method, PAthreader, to recognize remote templates and explore folding pathways. Firstly, we design a three-track alignment between predicted distance profiles and structure profiles extracted from PDB and AlphaFold DB, to improve the recognition accuracy of remote templates. Secondly, we improve the performance of AlphaFold2 using the templates identified by PAthreader. Thirdly, we explore protein folding pathways based on our conjecture that dynamic folding information of protein is implicitly contained in its remote homologs. The results show that the average accuracy of PAthreader templates is 11.6% higher than that of HHsearch. In terms of structure modelling, PAthreader outperform AlphaFold2 and ranks first on the CAMEO blind test for the latest three months. Furthermore, we predict protein folding pathways for 37 proteins, in which the results of 7 proteins are almost consistent with those of biological experiments, and the other 30 human proteins have yet to be verified by biological experiments, revealing that folding information can be exploited from remote homologous structures.

Engineering of flavonoid 3'-O-methyltransferase for improved biosynthesis of chrysoeriol in Escherichia coli

O-Methylation catalyzed by O-methyltransferases (OMTs) is an important modification of flavonoids for improving the transport efficiency across membranes and metabolic stability in mammalian cells. Chrysoeriol, also known as 3'-O-methylated luteolin, is a methylated flavonoid compound with health-promoting activities. The generation of chrysoeriol from luteolin can be catalyzed by a rice-derived 3'-OMT named ROMT-9, which has a high regiospecificity and activity toward flavonoids in vitro. Herein, we explored the potential of ROMT-9 for in vivo biosynthesis of chrysoeriol in Escherichia coli and adopted semi-rational enzyme engineering guided by homology modeling and molecular docking to improve the bio-production. Two positive variants including L34Q and W284A were obtained which promoted chrysoeriol formation to more than 85 mg/L from 200 mg/L of luteolin in 24 h compared with a titer of 55 mg/L for the strain expressing the native enzyme. Further biochemical analysis confirmed that such improvement in production stemmed from a higher enzyme expression level for the L34Q variant and higher efficiency in substrate binding and catalysis for the W284A variant. This study provides some insights into the engineering of other flavonoid OMTs and will facilitate high-level biosynthesis of methylated flavonoids in engineered microorganisms. KEY POINTS: • Biosynthesis of chrysoeriol from luteolin in E. coli using ROMT-9 • Engineering of ROMT-9 for better bio-production • ROMT-9 variants promote production via better expression or better catalysis.

Two O-Methyltransferases from Phylogenetically Unrelated Cow Parsley (Anthriscus sylvestris) and Hinoki-Asunaro (Thujopsis dolabrata var. hondae) as a Signature of Lineage-Specific Evolution in Lignan Biosynthesis

O-Methyltransferases (OMTs) play important roles in antitumor lignan biosynthesis. To date, six OMTs catalyzing the methylation of dibenzylbutyrolactone lignans as biosynthetic precursors of antitumor lignans have been identified. However, there is still no systematic understanding of the diversity and regularity of the biosynthetic mechanisms among various plant lineages. Herein, we report the characterization of two OMTs from Anthriscus sylvestris and Thujopsis dolabrata var. hondae [designated as AsSecoNorYatein (SNY) OMT and TdSNYOMT] together with the six known OMTs to evaluate their diversity and regularity. Although A. sylvestris 5-O-methylthujaplicatin (SecoNorYatein) and 4-O-demethylyatein (NorYatein) OMT (AsSNYOMT) and TdSNYOMT accept 5-O-methylthujaplicatin and 4-O-demethylyatein as substrates, phylogenetic analysis indicated that these two OMTs shared low amino acid sequence identity, 33.8%, indicating a signature of parallel evolution. The OMTs and the six previously identified OMTs were found to be diverse in terms of their substrate specificity, regioselectivity and amino acid sequence identity, indicating independent evolution in each plant species. Meanwhile, two-entropy analysis detected four amino acid residues as being specifically acquired by dibenzylbutyrolactone lignan OMTs. Site-directed mutation of AsSNYOMT indicated that two of them contributed specifically to 5-O-methylthujaplicatin methylation. The results provide a new example of parallel evolution and the diversity and regularity of OMTs in plant secondary (specialized) metabolism.

Thermostabilizing mechanisms of canonical single amino acid substitutions at a GH1 β-glucosidase probed by multiple MD and computational approaches

β-glucosidases play a pivotal role in second-generation biofuel (2G-biofuel) production. For this application, thermostable enzymes are essential due to the denaturing conditions on the bioreactors. Random amino acid substitutions have originated new thermostable β-glucosidases, but without a clear understanding of their molecular mechanisms. Here, we probe by different molecular dynamics simulation approaches with distinct force fields and submitting the results to various computational analyses, the molecular bases of the thermostabilization of the Paenibacillus polymyxa GH1 β-glucosidase by two-point mutations E96K (TR1) and M416I (TR2). Equilibrium molecular dynamic simulations (eMD) at different temperatures, principal component analysis (PCA), virtual docking, metadynamics (MetaDy), accelerated molecular dynamics (aMD), Poisson-Boltzmann surface analysis, grid inhomogeneous solvation theory and colony method estimation of conformational entropy allow to converge to the idea that the stabilization carried by both substitutions depend on different contributions of three classic mechanisms: (i) electrostatic surface stabilization; (ii) efficient isolation of the hydrophobic core from the solvent, with energetic advantages at the solvation cap; (iii) higher distribution of the protein dynamics at the mobile active site loops than at the protein core, with functional and entropic advantages. Mechanisms i and ii predominate for TR1, while in TR2, mechanism iii is dominant. Loop A integrity and loops A, C, D, and E dynamics play critical roles in such mechanisms. Comparison of the dynamic and topological changes observed between the thermostable mutants and the wildtype protein with amino acid co-evolutive networks and thermostabilizing hotspots from the literature allow inferring that the mechanisms here recovered can be related to the thermostability obtained by different substitutions along the whole family GH1. We hope the results and insights discussed here can be helpful for future rational approaches to the engineering of optimized β-glucosidases for 2G-biofuel production for industry, biotechnology, and science.

The Isoleucine at Position 118 in Transmembrane 2 Is Responsible for the Selectivity of Xamoterol, Nebivolol, and ICI89406 for the Human β1-Adrenoceptor

Known off-target interactions frequently cause predictable drug side-effects (e.g., β1-antagonists used for heart disease, risk β2-mediated bronchospasm). Computer-aided drug design would improve if the structural basis of existing drug selectivity was understood. A mutagenesis approach determined the ligand-amino acid interactions required for β1-selective affinity of xamoterol and nebivolol, followed by computer-based modeling to provide possible structural explanations. ³H-CGP12177 whole cell binding was conducted in Chinese hamster ovary cells stably expressing human β1, β2, and chimeric β1/β2-adrenoceptors (ARs). Single point mutations were investigated in transiently transfected cells. Modeling studies involved docking ligands into three-dimensional receptor structures and performing molecular dynamics simulations, comparing interaction frequencies between apo and holo structures of β1 and β2-ARs. From these observations, an ICI89406 derivative was investigated that gave further insights into selectivity. Stable cell line studies determined that transmembrane 2 was crucial for the β1-selective affinity of xamoterol and nebivolol. Single point mutations determined that the β1-AR isoleucine (I118) rather than the β2 histidine (H93) explained selectivity. Studies of other β1-ligands found I118 was important for ICI89406 selective affinity but not that for betaxolol, bisoprolol, or esmolol. Modeling studies suggested that the interaction energies and solvation of β1-I118 and β2-H93 are factors determining selectivity of xamoterol and ICI89406. ICI89406 without its phenyl group loses its high β1-AR affinity, resulting in the same affinity as for the β2-AR. The human β1-AR residue I118 is crucial for the β1-selective affinity of xamoterol, nebivolol, and ICI89406 but not all β1-selective compounds. SIGNIFICANCE STATEMENT: Some ligands have selective binding affinity for the human β1 versus the β2-adrenoceptor; however, the molecular/structural reason for this is not known. The transmembrane 2 residue isoleucine I118 is responsible for the selective β1-binding of xamoterol, nebivolol, and ICI89406 but does not explain the selective β1-binding of betaxolol, bisoprolol, or esmolol. Understanding the structural basis of selectivity is important to improve computer-aided ligand design, and targeting I118 in β1-adrenoceptors is likely to increase β1-selectivity of drugs.

In Silico Analysis of a Drosophila Parasitoid Venom Peptide Reveals Prevalence of the Cation-Polar-Cation Clip Motif in Knottin Proteins

As generalist parasitoid wasps, Leptopilina heterotoma are highly successful on many species of fruit flies of the genus Drosophila. The parasitoids produce specialized multi-strategy extracellular vesicle (EV)-like structures in their venom. Proteomic analysis identified several immunity-associated proteins, including the knottin peptide, LhKNOT, containing the structurally conserved inhibitor cysteine knot (ICK) fold, which is present in proteins from diverse taxa. Our structural and docking analysis of LhKNOT's 36-residue core knottin fold revealed that in addition to the knottin motif itself, it also possesses a Cation-Polar-Cation (CPC) clip. The CPC clip motif is thought to facilitate antimicrobial activity in heparin-binding proteins. Surprisingly, a majority of ICKs tested also possess the CPC clip motif, including 75 bona fide plant and arthropod knottin proteins that share high sequence and/or structural similarity with LhKNOT. Like LhKNOT and these other 75 knottin proteins, even the Drosophila Drosomycin antifungal peptide, a canonical target gene of the fly's Toll-NF-kappa B immune pathway, contains this CPC clip motif. Together, our results suggest a possible defensive function for the parasitoid LhKNOT. The prevalence of the CPC clip motif, intrinsic to the cysteine knot within the knottin proteins examined here, suggests that the resultant 3D topology is important for their biochemical functions. The CPC clip is likely a highly conserved structural motif found in many diverse proteins with reported heparin binding capacity, including amyloid proteins. Knottins are targets for therapeutic drug development, and insights into their structure-function relationships will advance novel drug design.

Amino Acid Polymorphisms on the Brazilian Strain of Yellow Fever Virus Methyltransferase Are Related to the Host's Immune Evasion Mediated by Type I Interferon

Since late 2016, a yellow fever virus (YFV) variant carrying a set of nine amino acid variations has circulated in South America. Three of them were mapped on the methyltransferase (MTase) domain of viral NS5 protein. To assess whether these changes affected viral infectivity, we synthesized YFV carrying the MTase of circulating lineage as well as its isoform with the residues of the previous strains (NS5 K101R, NS5 V138I, and NS5 G173S). We observed a slight difference in viral growth properties and plaque phenotype between the two synthetic YFVs. However, the MTase polymorphisms associated with the Brazilian strain of YFV (2016-2019) confer more susceptibility to the IFN-I. In addition, in vitro MTase assay revealed that the interaction between the YFV MTase and the methyl donor molecule (SAM) is altered in the Brazilian MTase variant. Altogether, the results reported here describe that the MTase carrying the molecular signature of the Brazilian YFV circulating since 2016 might display a slight decrease in its catalytic activity but virtually no effect on viral fitness in the parameters comprised in this study. The most marked influence of these residues stands in the immune escape against the antiviral response mediated by IFN-I.

Dynamic conformational states of apo, ATP and cabozantinib bound TAM kinases to differentiate active-inactive kinetic models

The dynamically active and inactive conformations of kinases play a crucial role in the activation of intracellular downstream signaling pathways. The all-atom molecular dynamics (MD) simulations at microsecond (µs) timescale and longer provide robust insights into the structural details of conformational alterations in kinases that contribute to their cellular metabolic activities and signaling pathways. Tyro3, Axl and Mer (TAM) receptor tyrosine kinases (RTKs) are overexpressed in several types of human cancers. Cabozantinib, a small molecule inhibitor constrains the activity of TAM kinases at nanomolar concentrations. The apo, complexes of ATP (active state) and cabozantinib (active and inactive states) with TAM RTKs were studied by 1 µs MD simulations followed by trajectory analyses. The dynamic mechanistic pathways intrinsic to the kinase activity and protein conformational landscape in the cabozantinib bound TAM kinases are revealed due to the alterations in the P-loop, α-helix and activation loop that result in breaking the regulatory (R) and catalytic (C) spines, while the active states of ATP bound TAM kinases are retained. The co-existence of dynamical states when bound to cabozantinib was observed and the long-lived kinetic transition states of distinct active and inactive structural models were deciphered from MD simulation trajectories that have not been revealed so far.Communicated by Ramaswamy H. Sarma.

Molecular dynamics studies of CED-4/CED-9/EGL-1 ternary complex reveal CED-4 release mechanism in the linear apoptotic pathway of Caenorhabditis elegans

Many steps in programmed cell death are evolutionarily conserved across different species. The Caenorhabditis elegans proteins CED-9, CED-4 and EGL-1 involved in apoptosis are respectively homologous to anti-apoptotic Bcl-2 proteins, Apaf-1 and the "BH3-only" pro-apototic proteins in mammals. In the linear apoptotic pathway of C. elegans, EGL-1 binding to CED-9 leads to the release of CED-4 from CED-9/CED-4 complex. The molecular events leading to this process are not clearly elucidated. While the structures of CED-9 apo, CED-9/EGL-1 and CED-9/CED-4 complexes are known, the CED-9/CED-4/EGL-1 ternary complex structure is not yet determined. In this work, we modeled this ternary complex and performed molecular dynamics simulations of six different systems involving CED-9. CED-9 displays differential dynamics depending upon whether it is bound to CED-4 and/or EGL-1. CED-4 exists as an asymmetric dimer (CED4a and CED4b) in CED-9/CED-4 complex. CED-4a exhibits higher conformational flexibility when simulated without CED-4b. Principal Component Analysis revealed that the direction of CED-4a's winged-helix domain motion differs in the ternary complex. Upon EGL-1 binding, majority of non-covalent interactions involving CARD domain in the CED-4a-CED-9 interface have weakened and only half of the contacts found in the crystal structure between α/β domain of CED4a and CED-9 are found to be stable. Additional stable contacts in the ternary complex and differential dynamics indicate that winged-helix domain may play a key role in CED-4a's dissociation from CED-9. This study has provided a molecular level understanding of potential intermediate states that are likely to occur when CED-4a is released from CED-9.

Molecular evolution and structural analyses of proteins involved in metabolic pathways of volatile organic compounds in Petunia hybrida (Solanaceae)

The association between plants and their pollinators is essential for increasing the diversity in angiosperms. Morphological and physiological traits, mainly floral scent, can influence the pollination dynamics and select pollinators for each plant species. In this work, we studied two proteins involved in producing volatile organic compounds in plants, conyferyl alcohol acyltransferase (CFAT) and benzoyl-CoA:benzyl alcohol/phenyl ethanol benzoyl transferase (BPBT) genes. We aimed to understand these proteins with respect to evolutionary and structural aspects and functions in Solanaceae using phylogenetic methods and comparative molecular modeling. We used Bayesian inference to describe the proteins' evolutionary history using Petunia x hybrida as a query to search for homologs in the Solanaceae family. Theoretical 3D models were obtained for both proteins using Panicum virgatum as a template. The phylogenetic tree included several different enzymes with diverse biological roles in Solanaceae, displaying the transferase domain. We identified only one sequence of CFAT in the databases, which belongs to Petunia x hybrida, and found several BPBT sequences from the genera Nicotiana, Solanum, and Capsicum. The 3D structures of CFAT and BPBT have two different domains, and we have identified the amino acid residues essential for the enzymatic activity and interaction with substrates.

Funnel metadynamics and behavioral studies reveal complex effect of D2AAK1 ligand on anxiety-like processes

Anxiety is a troublesome symptom for many patients, especially those suffering from schizophrenia. Its regulation involves serotonin receptors, targeted e.g. by antipsychotics or psychedelics such as LSD. 5-HT_2A receptors are known for an extremely long LSD residence time, enabling minute doses to exert a long-lasting effect. In this work, we explore the changes in anxiety-like processes induced by the previously reported antipsychotic, D2AAK1. In vivo studies revealed that the effect of D2AAK1 on the anxiety is mediated through serotonin 5-HT_1A and 5-HT_2A receptors, and that it is time-dependent (anxiogenic after 30 min, anxiolytic after 60 min) and dose-dependent. The funnel metadynamics simulations suggest complicated ligand-5HT_2AR interactions, involving an allosteric site located under the third extracellular loop, which is a possible explanation of the time-dependency. The binding of D2AAK1 at the allosteric site results in a broader opening of the extracellular receptor entry, possibly altering the binding kinetics of orthosteric ligands.

A rare loss-of-function genetic mutation suggest a role of dermcidin deficiency in hidradenitis suppurativa pathogenesis

Hidradenitis suppurativa (HS) is a chronic inflammatory skin disease with a multifactorial aetiology that involves a strict interplay between genetic factors, immune dysregulation and lifestyle. Familial forms represent around 40% of total HS cases and show an autosomal dominant mode of inheritance of the disease. In this study, we conducted a whole-exome sequence analysis on an Italian family of 4 members encompassing a vertical transmission of HS. Focusing on rare damaging variants, we identified a rare insertion of one nucleotide (c.225dupA:p.A76Sfs*21) in the DCD gene encoding for the antimicrobial peptide dermcidin (DCD) that was shared by the proband, his affected father and his 11-years old daughter. Since several transcriptome studies have shown a significantly decreased expression of DCD in HS skin, we hypothesised that the identified frameshift insertion was a loss-of-function mutation that might be associated with HS susceptibility in this family. We thus confirmed by mass spectrometry that DCD levels were diminished in the affected members and showed that the antimicrobial activity of a synthetic DCD peptide resulting from the frameshift mutation was impaired. In order to define the consequences related to a decrease in DCD activity, skin microbiome analyses of different body sites were performed by comparing DCD mutant and wild type samples, and results highlighted significant differences between the groins of mutated and wild type groups. Starting from genetic analysis conducted on an HS family, our findings showed, confirming previous transcriptome results, the potential role of the antimicrobial DCD peptide as an actor playing a crucial part in the etio-pathogenesis of HS and in the maintenance of the skin's physiological microbiome composition; so, we can hypothesise that DCD could be used as a novel target for personalised therapeutic approach.

Evaluating Allosteric Perturbations in Cannabinoid Receptor 1 by In Silico Single-Point Mutation

Cannabinoid receptor 1 (CB1) is a promising drug target involved in many physiological processes. Using atomistic molecular dynamics (MD) simulations, we examined the structural effect of F237L mutation on CB1, a mutation that has qualitatively similar effects to allosteric ligand ORG27569 binding. This mutation showed a global effect on CB1 conformations. Among the observed effects, TM6 outward movement and the conformational change of the NPxxY motif upon receptor activation by CB1 agonist CP55940 were hindered compared to wt CB1. Within the orthosteric binding site, CP55940 interactions with CB1 were altered. Our results revealed that allosteric perturbations introduced by the mutation had a global impact on receptor conformations, suggesting that the mutation site is a key region for allosteric modulation in CB1.

Serial Passaging of Seasonal H3N2 Influenza A/Singapore/G2-31.1/2014 Virus in MDCK-SIAT1 Cells and Primary Chick Embryo Cells Generates HA D457G Mutation and Other Variants in HA, NA, PB1, PB1-F2, and NS1

Influenza remains one of the most prevalent viruses circulating amongst humans and has resulted in several pandemics. The prevention and control of H3N2 influenza is complicated by its propensity for evolution, which leads to vaccine mismatch and reduced vaccine efficacies. This study employed the strategy of serial passaging to compare the evolution of the human seasonal influenza strain A/Singapore/G2-31.1/2014(H3N2) in MDCK-SIAT1 versus primary chick embryo fibroblast (CEF) cells. Genetic analysis of the HA, NS1, NA, and PB1 gene segments by Sanger sequencing revealed the presence of specific mutations and a repertoire of viral quasispecies following serial passaging. Most quasispecies were also found in PB1, which exhibited consistently high transversion-to-transition ratios in all five MDCK-SIAT1 passages. Most notably, passage 5 virus harbored the D457G substitution in the HA2 subunit, while passage 3 virus acquired K53Q and Q69H mutations in PB1-F2. An A971 variant leading to a non-synonymous R316Q substitution in PB1 was also identified in MDCK-SIAT1 passages 2 and 4. With an increasing number of passages, the proportion of D457G mutations progressively increased and was associated with larger virus plaque sizes. However, microneutralization assays revealed no significant differences in the neutralizing antibody profiles of human-influenza-immune serum samples against pre-passaged virus and passage 5 virus. In contrast, viable virus was only detected in passage 1 of CEF cells, which gave rise to multiple viral RNA quasispecies. Our findings highlight that serial passaging is able to drive differential adaptation of H3N2 influenza in different host species and may alter viral virulence. More studies are warranted to elucidate the complex relationships between H3N2 virus evolution, viral virulence changes, and low vaccine efficacy.

Comparison of herbicide specificity of CYP81A cytochrome P450s from rice and a multiple-herbicide resistant weed, Echinochloa phyllopogon

BACKGROUND

CYP81A cytochrome P450s (CYP81As) play a key role in herbicide detoxification in Poaceae plants. Crop CYP81As confer natural tolerance to multiple herbicides, whereas CYP81As in weeds disrupt herbicide action. Identifying differences in CYP81A herbicide specificity between crops and weeds could provide valuable information for controlling weeds. In this study, we quantitatively compared herbicide specificity between CYP81A6 from rice (Oryza sativa) and CYP81A12 and CYP81A21 from a weed, Echinochloa phyllopogon, using transgenic Escherichia coli and Arabidopsis.

RESULTS

All three CYP81As metabolized the five tested herbicides and formed similar metabolites, with the highest relative activities of 400 to 580% toward bentazone compared to their activity on bensulfuron-methyl (defined as 100%). However, they showed differing activity toward propyrisulfuron. The relative activities of Echinochloa phyllopogon CYP81A12 (12.2%) and CYP81A21 (34.4%) toward propyrisulfuron were lower than that of rice CYP81A6 (98.5%). Additionally, rice CYP81A6 produced O-demethylated propyrisulfuron and hydroxylated products, whereas Echinochloa phyllopogon CYP81As produced only hydroxylated products. Arabidopsis expressing CYP81A12 and CYP81A21 exhibited lower levels of resistance against propyrisulfuron compared to that in Arabidopsis expressing CYP81A6. Homology modeling and in silico docking revealed that bensulfuron-methyl docked well into the active centers of all three CYP81As, whereas propyrisulfuron docked into rice CYP81A6 but not into Echinochloa phyllopogon CYP81As.

CONCLUSION

The differing herbicide specificity displayed by rice CYP81A6 and Echinochloa phyllopogon CYP81A12 and CYP81A21 will help design inhibitors (synergists) of weed CYP81As, as well as develop novel herbicide ingredients that are selectively metabolized by crop CYP81As, to overcome the problem of herbicide resistance. © 2022 Society of Chemical Industry.

Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

Characterizations of a novel peptide encoded by a circular RNA using in-silico analyses

CircRNAs have gained importance in recent times due to their involvement in gene regulation and also in the prognosis of cancer. Generally, the circRNA directly interact with miRNA or RNA binding proteins to exert their action, but some of them can be translated. These translated peptides often participate in the regulation of cellular processes. The circPPP1R12A translated peptide has been shown to influence the functioning of the Mst pathway. The Mst signaling is noteworthy for its role in the process of development, but it also has a function as a regulator of apoptosis, which is significant for regulation of cancer. Overexpression of this novel peptide deactivates the Mst signaling to induce the expression of the proliferative oncogene, Yap. Its molecular interaction with the molecules in the Mst pathway is hitherto unknown. In this short report we present our findings from in-silico studies the plausible structure of the peptide through bioinformatics and dynamics simulation studies. This is the first such report on the structure of the novel peptide encoded by circPPP1R12A, which could be important to predict in future its molecular interactions to understand its functionality.

Quaternary organization of the human eEF1B complex reveals unique multi-GEF domain assembly

Protein synthesis in eukaryotic cell is spatially and structurally compartmentalized that ensures high efficiency of this process. One of the distinctive features of higher eukaryotes is the existence of stable multi-protein complexes of aminoacyl-tRNA synthetases and translation elongation factors. Here, we report a quaternary organization of the human guanine-nucleotide exchange factor (GEF) complex, eEF1B, comprising α, β and γ subunits that specifically associate into a heterotrimeric form eEF1B(αβγ)3. As both the eEF1Bα and eEF1Bβ proteins have structurally conserved GEF domains, their total number within the complex is equal to six. Such, so far, unique structural assembly of the guanine-nucleotide exchange factors within a stable complex may be considered as a 'GEF hub' that ensures efficient maintenance of the translationally active GTP-bound conformation of eEF1A in higher eukaryotes.

MAIA, Fc receptor-like 3, supersedes JUNO as IZUMO1 receptor during human fertilization

Gamete fusion is a critical event of mammalian fertilization. A random one-bead one-compound combinatorial peptide library represented synthetic human egg mimics and identified a previously unidentified ligand as Fc receptor-like 3, named MAIA after the mythological goddess intertwined with JUNO. This immunoglobulin super family receptor was expressed on human oolemma and played a major role during sperm-egg adhesion and fusion. MAIA forms a highly stable interaction with the known IZUMO1/JUNO sperm-egg complex, permitting specific gamete fusion. The complexity of the MAIA isotype may offer a cryptic sexual selection mechanism to avoid genetic incompatibility and achieve favorable fitness outcomes.

Delineating specific regions of N- terminal domain of T3SS ATPase YsaN of Yersinia enterocolitica governing its different oligomerization states

Oligomerization of YsaN, a putative T3SS-ATPase is a necessary and crucial event for T3SS functioning in Y. enterocolitica. Different oligomeric states have been proposed for similar ATPases, yet, the true nature of its activation and formation of different oligomers is still poorly understood. In-vitro studies of YsaN reveal that its activation and oligomerization depend on its N-terminal region and occur as a result of active catalysis of ATP in an ATP concentration-dependent manner following two-step cooperative kinetics. Also, the N-terminal 83 amino acid residues of YsaN are crucial for higher-order oligomer formation while YsaN∆83 is capable of hexamer formation upon oligomerization. Enzyme kinetics study shows reduced ATPase activity of YsaN∆83 (3.19 ± 0.09 μmol/min/mg) in comparison to YsaN (9.076 ± 0.72 μmol/min/mg). Negative-TEM study of YsaN and YsaN∆83 oligomer suggests that the formation of higher-order oligomer (probably dodecamer) occurs by stacking of two hexamers through their N-terminal faces involving N-terminal 83 amino acid residues which have been further supported by the docking of two hexamers during the in-silico study. These results suggest that YsaN is an oligomerization-activated T3SS ATPase, where distinct regions of its N-terminal domain regulate its different oligomeric nature and is essential for its activation.

Structure and dynamics of Toll immunoreceptor activation in the mosquito Aedes aegypti

Aedes aegypti has evolved to become an efficient vector for arboviruses but the mechanisms of host-pathogen tolerance are unknown. Immunoreceptor Toll and its ligand Spaetzle have undergone duplication which may allow neofunctionalization and adaptation. Here we present cryo-EM structures and biophysical characterisation of low affinity Toll5A complexes that display transient but specific interactions with Spaetzle1C, forming asymmetric complexes, with only one ligand clearly resolved. Loop structures of Spaetzle1C and Toll5A intercalate, temporarily bridging the receptor C-termini to promote signalling. By contrast unbound receptors form head-to-head homodimers that keep the juxtamembrane regions far apart in an inactive conformation. Interestingly the transcriptional signature of Spaetzle1C differs from other Spaetzle cytokines and controls genes involved in innate immunity, metabolism and tissue regeneration. Taken together our results explain how upregulation of Spaetzle1C in the midgut and Toll5A in the salivary gland shape the concomitant immune response.

Connective tissue growth factor-targeting DNA aptamer suppresses pannus formation as diagnostics and therapeutics for rheumatoid arthritis

Connective tissue growth factor (CTGF) has been recently acknowledged as an ideal biomarker in the early disease course, participating in the pathogenesis of pannus formation in rheumatoid arthritis (RA). However, existing approaches for the detection of or antagonist targeting CTGF are either lacking or unsatisfactory in the diagnosis and treatment of RA. To address this, we synthesized and screened high-affinity single-stranded DNA aptamers targeting CTGF through a protein-based SELEX procedure. The structurally optimized variant AptW2-1-39-PEG was characterized thoroughly for its high-affinity (KD 7.86 nM), sensitivity (minimum protein binding concentration, 2 ng), specificity (negative binding to other biomarkers of RA), and stability (viability-maintaining duration in human serum, 48 h) properties using various biochemical and biophysical assays. Importantly, we showed the antiproliferative and antiangiogenic activities of the aptamers obtained using functional experiments and further verified the therapeutic effect of the aptamers on joint injury and inflammatory response in collagen-induced arthritis (CIA) mice, thus advancing this study into actual therapeutic application. Furthermore, we revealed that the binding within AptW2-1-39-PEG/CTGF was mediated by the thrombospondin 1 (TSP1) domain of CTGF using robust bioinformatics tools together with immunofluorescence. In conclusion, our results revealed a novel aptamer that holds promise as an additive or alternative approach for CTGF-targeting diagnostics and therapeutics for RA.

Expression of a novel mycobacterial phosphodiesterase successfully lowers cAMP levels resulting in reduced tolerance to cell wall-targeting antimicrobials

cAMP and antimicrobial susceptibility in mycobacteriaAntimicrobial tolerance, the ability to survive exposure to antimicrobials via transient nonspecific means, promotes the development of antimicrobial resistance (AMR). The study of the molecular mechanisms that result in antimicrobial tolerance is therefore essential for the understanding of AMR. In gram-negative bacteria, the second messenger molecule 3'',5''-cAMP has been previously shown to be involved in AMR. In mycobacteria, however, the role of cAMP in antimicrobial tolerance has been difficult to probe due to its particular complexity. In order to address this difficulty, here, through unbiased biochemical approaches consisting in the fractionation of clear protein lysate from a mycobacterial strain deleted for the known cAMP phosphodiesterase (Rv0805c) combined with mass spectrometry techniques, we identified a novel cyclic nucleotide-degrading phosphodiesterase enzyme (Rv1339) and developed a system to significantly decrease intracellular cAMP levels through plasmid expression of Rv1339 using the constitutive expression system, pVV16. In Mycobacterium smegmatis mc²155, we demonstrate that recombinant expression of Rv1339 reduced cAMP levels threefold and resulted in altered gene expression, impaired bioenergetics, and a disruption in peptidoglycan biosynthesis leading to decreased tolerance to antimicrobials that target cell wall synthesis such as ethambutol, D-cycloserine, and vancomycin. This work increases our understanding of the role of cAMP in mycobacterial antimicrobial tolerance, and our observations suggest that nucleotide signaling may represent a new target for the development of antimicrobial therapies.

RET c.1901G>A and Novel SLC12A3 Mutations in Familial Pheochromocytomas

Familial PHEOs (pheochromocytomas) are inherited as an autosomal dominant trait, and inherited PHEOs can be one clinical phenotype of clinical syndromes, such as multiple endocrine neoplasia type 2A (MEN2A). In recent years, there has been a lot of controversy about the factors affecting the penetrance of PHEOs in MEN2A, of which the effects of RET (rearranged during transfection) proto-oncogene mutations are the primary concern. In this report, we performed genetic screening of patients in one family presenting with PHEOs and found they carried a RET c.1901G>A mutation. They were ultimately diagnosed with familial MEN2A. We found that MEN2A patients with the RET c.1901G>A mutation tended to have bilateral PHEOs that appeared earlier than medullary thyroid carcinoma. Genetic analysis showed that the patients also carried novel SLC12A3 (solute carrier family 12 member 3) variants, which are highly associated with Giteman syndrome. The results of protein structure prediction models suggest this SLC12A3 mutant has altered both the protein structure and the interaction with surrounding amino acids. Further studies of the phenotypes and related mechanisms of the gene mutations are required to guide individual assessment and treatment.

Protein interactions with metallothionein-3 promote vectorial active transport in human proximal tubular cells

Metallothionein 3 (MT-3) is a small, cysteine-rich protein that binds to essential metals required for homeostasis, as well as to heavy metals that have the potential to exert toxic effects on cells. MT-3 is expressed by epithelial cells of the human kidney, including the cells of the proximal tubule. Our laboratory has previously shown that mortal cultures of human proximal tubular (HPT) cells express MT-3 and form domes in the cell monolayer, a morphological feature indicative of vectorial active transport, an essential function of the proximal tubule. However, an immortalized proximal tubular cell line HK-2 lacks the expression of MT-3 and fails to form domes in the monolayer. Transfection of HK-2 cells with the MT-3 gene restores dome formation in these cells suggesting that MT-3 is required for vectorial active transport. In order to determine how MT-3 imparts this essential feature to the proximal tubule, we sought to identify proteins that interact either directly or indirectly with MT-3. Using a combination of pulldowns, co-immunoprecipitations, and mass spectrometry analysis, putative protein interactants were identified and subsequently confirmed by Western analysis and confocal microscopy, following which proteins with direct physical interactions were investigated through molecular docking. Our data shows that MT-3 interacts with myosin-9, aldolase A, enolase 1, β-actin, and tropomyosin 3 and that these interactions are maximized at the periphery of the apical membrane of doming proximal tubule cells. Together these observations reveal that MT-3 interacts with proteins involved in cytoskeletal organization and energy metabolism, and these interactions at the apical membrane support vectorial active transport and cell differentiation in proximal tubule cultures.

Migraine-Associated Mutation in the Na,K-ATPase Leads to Disturbances in Cardiac Metabolism and Reduced Cardiac Function

Background

Mutations in ATP1A2 gene encoding the Na,K‐ATPase α₂ isoform are associated with familial hemiplegic migraine type 2. Migraine with aura is a known risk factor for heart disease. The Na,K‐ATPase is important for cardiac function, but its role for heart disease remains unknown. We hypothesized that ATP1A2 is a susceptibility gene for heart disease and aimed to assess the underlying disease mechanism.

Methods and Results

Mice heterozygous for the familial hemiplegic migraine type 2–associated G301R mutation in the Atp1a2 gene (α₂^+/G301R mice) and matching wild‐type controls were compared. Reduced expression of the Na,K‐ATPase α₂ isoform and increased expression of the α₁ isoform were observed in hearts from α₂^+/G301R mice (Western blot). Left ventricular dilation and reduced ejection fraction were shown in hearts from 8‐month‐old α₂^+/G301R mice (cardiac magnetic resonance imaging), and this was associated with reduced nocturnal blood pressure (radiotelemetry). Cardiac function and blood pressure of 3‐month‐old α₂^+/G301R mice were similar to wild‐type mice. Amplified Na,K‐ATPase–dependent Src kinase/Ras/Erk1/2 (p44/42 mitogen‐activated protein kinase) signaling was observed in hearts from 8‐month‐old α₂^+/G301R mice, and this was associated with mitochondrial uncoupling (respirometry), increased oxidative stress (malondialdehyde measurements), and a heart failure–associated metabolic shift (hyperpolarized magnetic resonance). Mitochondrial membrane potential (5,5´,6,6´‐tetrachloro‐1,1´,3,3´‐tetraethylbenzimidazolocarbocyanine iodide dye assay) and mitochondrial ultrastructure (transmission electron microscopy) were similar between the groups. Proteomics of heart tissue further suggested amplified Src/Ras/Erk1/2 signaling and increased oxidative stress and provided the molecular basis for systolic dysfunction in 8‐month‐old α₂^+/G301R mice.

Conclusions

Our findings suggest that ATP1A2 mutation leads to disturbed cardiac metabolism and reduced cardiac function mediated via Na,K‐ATPase–dependent reactive oxygen species signaling through the Src/Ras/Erk1/2 pathway.

Pregnane-based steroids are novel positive NMDA receptor modulators that may compensate for the effect of loss-of-function disease-associated GRIN mutations

BACKGROUND AND PURPOSE

N-methyl-D-aspartate receptors (NMDARs) play a critical role in synaptic plasticity, and mutations in human genes encoding NMDAR subunits have been described in individuals with various neuropsychiatric disorders. Compounds with a positive allosteric effect are thought to compensate for reduced receptor function.

EXPERIMENTAL APPROACH

We have used whole-cell patch-clamp electrophysiology on recombinant rat NMDARs and human variants found in individuals with neuropsychiatric disorders, in combination with in silico modelling, to explore the site of action of novel epipregnanolone-based NMDAR modulators.

KEY RESULTS

Analysis of the action of 4-(20-oxo-5β-pregnan-3β-yl) butanoic acid (EPA-But) at the NMDAR indicates that the effect of this steroid with a "bent" structure is different from that of cholesterol and oxysterols and shares a disuse-dependent mechanism of NMDAR potentiation with the "planar" steroid 20-oxo-pregn-5-en-3β-yl sulfate (PE-S). The potentiating effects of EPA-But and PE-S are additive. Alanine scan mutagenesis identified residues that reduce the potentiating effect of EPA-But. No correlation was found between the effects of EPA-But and PE-S at mutated receptors that were less sensitive to either steroid. The relative degree of potentiation induced by the two steroids also differed in human NMDARs carrying rare variants of hGluN1 or hGluN2B subunits found in individuals with neuropsychiatric disorders, including intellectual disability, epilepsy, developmental delay, and autism spectrum disorder.

CONCLUSION AND IMPLICATIONS

Our results show novel sites of action for pregnanolones at the NMDAR and provide an opportunity for the development of new therapeutic neurosteroid-based ligands to treat diseases associated with glutamatergic system hypofunction.

Binding Affinity Characterization of Four Antennae-Enriched Odorant-Binding Proteins From Harmonia axyridis (Coleoptera: Coccinellidae)

Harmonia axyridis is an important natural enemy that consumes many agricultural and forestry pests. It relies on a sensitive olfactory system to find prey and mates. Odorant-binding proteins (OBPs) as the first-step of recognizing volatiles, transport odors through sensillum lymph to odorant receptors (ORs). However, little is known about the molecular mechanisms of H. axyridis olfaction. In this study, four H. axyridis antenna specific OBP genes, HaxyOBP3, 5, 12, and 15, were bacterially expressed and the binding features of the four recombinant proteins to 40 substances were investigated using fluorescence competitive binding assays. Three-dimensional structure modeling and molecular docking analysis predicted the binding sites between HaxyOBPs and candidate volatiles. Developmental expression analyses showed that the four HaxyOBP genes displayed a variety of expression patterns at different development stages. The expression levels of HaxyOBP3 and HaxyOBP15 were higher in the adult stage than in the other developmental stages, and HaxyOBP15 was significantly transcriptionally enriched in adult stage. Ligand-binding analysis demonstrated that HaxyOBP3 and HaxyOBP12 only combined with two compounds, β-ionone and p-anisaldehyde. HaxyOBP5 protein displayed binding affinities with methyl salicylate, β-ionone, and p-anisaldehyde (K_i = 18.15, 11.71, and 13.45 μM). HaxyOBP15 protein had a broad binding profile with (E)-β-farnesene, β-ionone, α-ionone, geranyl acetate, nonyl aldehyde, dihydro-β-ionone, and linalyl acetate (K_i = 4.33–31.01 μM), and hydrophobic interactions played a key role in the binding of HaxyOBP15 to these substances according to molecular docking. Taken together, HaxyOBP15 exhibited a broader ligand-binding spectrum and a higher expression in adult stage than HaxyOBP3, 5, and 12, indicating HaxyOBP15 may play a greater role in binding volatiles than other three HaxyOBPs. The results will increase our understanding of the molecular mechanism of H. axyridis olfaction and may also result in new management strategies (attractants/repellents) that increase the biological control efficacy of H. axyridis.

An in-frame deletion mutation in the degron tail of auxin coreceptor IAA2 confers resistance to the herbicide 2,4-D in Sisymbrium orientale

Synthetic auxin herbicides intersect basic plant developmental biology and applied weed management. We investigated resistance to 2,4-D in the Australian weed Sisymbrium orientale (Indian hedge mustard). We identified a mechanism involving an in-frame 27-bp deletion in the degron tail of auxin coreceptor IAA2, one member of the gene family of Aux/IAA auxin co-receptors. We show that this deletion in IAA2 is a gain-of-function mutation that confers synthetic auxin resistance. This field-evolved mechanism of resistance to synthetic auxin herbicides confirms previous biochemical studies showing the role of the Aux/IAA degron tail in regulating Aux/IAA protein degradation upon auxin perception. The deletion mutation could be generated in crops using gene-editing approaches for cross-resistance to multiple synthetic auxin herbicides.

The natural auxin indole-3-acetic acid (IAA) is a key regulator of many aspects of plant growth and development. Synthetic auxin herbicides such as 2,4-D mimic the effects of IAA by inducing strong auxinic-signaling responses in plants. To determine the mechanism of 2,4-D resistance in a Sisymbrium orientale (Indian hedge mustard) weed population, we performed a transcriptome analysis of 2,4-D-resistant (R) and -susceptible (S) genotypes that revealed an in-frame 27-nucleotide deletion removing nine amino acids in the degron tail (DT) of the auxin coreceptor Aux/IAA2 (SoIAA2). The deletion allele cosegregated with 2,4-D resistance in recombinant inbred lines. Further, this deletion was also detected in several 2,4-D-resistant field populations of this species. Arabidopsis transgenic lines expressing the SoIAA2 mutant allele were resistant to 2,4-D and dicamba. The IAA2-DT deletion reduced binding to TIR1 in vitro with both natural and synthetic auxins, causing reduced association and increased dissociation rates. This mechanism of synthetic auxin herbicide resistance assigns an in planta function to the DT region of this Aux/IAA coreceptor for its role in synthetic auxin binding kinetics and reveals a potential biotechnological approach to produce synthetic auxin-resistant crops using gene-editing.