skip to main content


This content will become publicly available on May 1, 2024

Title: Crystal structure of MbnF: an NADPH-dependent flavin monooxygenase from Methylocystis strain SB2
Methanobactins (MBs) are ribosomally produced and post-translationally modified peptides (RiPPs) that are used by methanotrophs for copper acquisition. The signature post-translational modification of MBs is the formation of two heterocyclic groups, either an oxazolone, pyrazinedione or imidazolone group, with an associated thioamide from an X -Cys dipeptide. The precursor peptide (MbnA) for MB formation is found in a gene cluster of MB-associated genes. The exact biosynthetic pathway of MB formation is not yet fully understood, and there are still uncharacterized proteins in some MB gene clusters, particularly those that produce pyrazinedione or imidazolone rings. One such protein is MbnF, which is proposed to be a flavin monooxygenase (FMO) based on homology. To help to elucidate its possible function, MbnF from Methylocystis sp. strain SB2 was recombinantly produced in Escherichia coli and its X-ray crystal structure was resolved to 2.6 Å resolution. Based on its structural features, MbnF appears to be a type A FMO, most of which catalyze hydroxylation reactions. Preliminary functional characterization shows that MbnF preferentially oxidizes NADPH over NADH, supporting NAD(P)H-mediated flavin reduction, which is the initial step in the reaction cycle of several type A FMO enzymes. It is also shown that MbnF binds the precursor peptide for MB, with subsequent loss of the leader peptide sequence as well as the last three C-terminal amino acids, suggesting that MbnF might be needed for this process to occur. Finally, molecular-dynamics simulations revealed a channel in MbnF that is capable of accommodating the core MbnA fragment minus the three C-terminal amino acids.  more » « less
Award ID(s):
1912482
NSF-PAR ID:
10433420
Author(s) / Creator(s):
; ; ; ; ; ; ; ;
Date Published:
Journal Name:
Acta Crystallographica Section F Structural Biology Communications
Volume:
79
Issue:
5
ISSN:
2053-230X
Page Range / eLocation ID:
111 to 118
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Kelly, Robert M. (Ed.)
    ABSTRACT Methanobactins (MBs) are ribosomally synthesized and posttranslationally modified peptides (RiPPs) produced by methanotrophs for copper uptake. The posttranslational modification that defines MBs is the formation of two heterocyclic groups with associated thioamines from X-Cys dipeptide sequences. Both heterocyclic groups in the MB from Methylosinus trichosporium OB3b (MB-OB3b) are oxazolone groups. The precursor gene for MB-OB3b is mbnA , which is part of a gene cluster that contains both annotated and unannotated genes. One of those unannotated genes, mbnC , is found in all MB operons and, in conjunction with mbnB , is reported to be involved in the formation of both heterocyclic groups in all MBs. To determine the function of mbnC , a deletion mutation was constructed in M. trichosporium OB3b, and the MB produced from the Δ mbnC mutant was purified and structurally characterized by UV-visible absorption spectroscopy, mass spectrometry, and solution nuclear magnetic resonance (NMR) spectroscopy. MB-OB3b from the Δ mbnC mutant was missing the C-terminal Met and was also found to contain a Pro and a Cys in place of the pyrrolidinyl-oxazolone-thioamide group. These results demonstrate MbnC is required for the formation of the C-terminal pyrrolidinyl-oxazolone-thioamide group from the Pro-Cys dipeptide, but not for the formation of the N-terminal 3-methylbutanol-oxazolone-thioamide group from the N-terminal dipeptide Leu-Cys. IMPORTANCE A number of environmental and medical applications have been proposed for MBs, including bioremediation of toxic metals and nanoparticle formation, as well as the treatment of copper- and iron-related diseases. However, before MBs can be modified and optimized for any specific application, the biosynthetic pathway for MB production must be defined. The discovery that mbnC is involved in the formation of the C-terminal oxazolone group with associated thioamide but not for the formation of the N-terminal oxazolone group with associated thioamide in M. trichosporium OB3b suggests the enzymes responsible for posttranslational modification(s) of the two oxazolone groups are not identical. 
    more » « less
  2. Federle, Michael J. ; Dunny, Gary M. (Ed.)
    ABSTRACT Streptococcus gallolyticus subsp. gallolyticus is an emerging opportunistic pathogen responsible for septicemia and endocarditis in the elderly. Invasive infections by S. gallolyticus subsp. gallolyticus are strongly linked to the occurrence of colorectal cancer (CRC). It was previously shown that increased secondary bile salts under CRC conditions enhance the bactericidal activity of gallocin, a bacteriocin produced by S. gallolyticus subsp. gallolyticus , enabling it to colonize the mouse colon by outcompeting resident enterococci (L. Aymeric, F. Donnadieu, C. Mulet, L. du Merle, et al., Proc Natl Acad Sci U S A 115:E283–E291, 2018, https://doi.org/10.1073/pnas.1715112115 ). In a separate study, we showed that S. gallolyticus subsp. gallolyticus produces and secretes a 21-mer peptide that activates bacteriocin production (A. Proutière, L. du Merle, B. Périchon, H. Varet, et al., mBio 11:e03187-20, 2020, https://doi.org/10.1128/mBio.03187-20 ). This peptide was named CSP because of its sequence similarity with competence-stimulating peptides found in other streptococci. Here, we demonstrate that CSP is a bona fide quorum sensing peptide involved in activation of gallocin gene transcription. We therefore refer to CSP as GSP (gallocin-stimulating peptide). GSP displays some unique features, since its N-terminal amino acid lies three residues after the double glycine leader sequence. Here, we set out to investigate the processing and export pathway that leads to mature GSP. Heterologous expression in Lactococcus lactis of the genes encoding GSP and the BlpAB transporter is sufficient to produce the 21-mer form of GSP in the supernatant, indicating that S. gallolyticus subsp. gallolyticus BlpAB displays an atypical cleavage site. We also conducted the first comprehensive structure-activity relationship (SAR) analysis of S. gallolyticus subsp. gallolyticus GSP to identify its key structural features and found that unlike many other similar streptococci signaling peptides (such as CSPs), nearly half of the mature GSP sequence can be removed (residues 1 to 9) without significantly impacting the peptide activity. IMPORTANCE Streptococcus gallolyticus subsp. gallolyticus is an opportunistic pathogen associated with colorectal cancer (CRC) and endocarditis. S. gallolyticus subsp. gallolyticus utilizes quorum sensing (QS) to regulate the production of a bacteriocin (gallocin) and gain a selective advantage in colonizing the colon. In this article, we report (i) the first structure-activity relationship study of the S. gallolyticus subsp. gallolyticus QS pheromone that regulates gallocin production, (ii) evidence that the active QS pheromone is processed to its mature form by a unique ABC transporter and not processed by an extracellular protease, and (iii) supporting evidence of interspecies interactions between streptococcal pheromones. Our results revealed the minimal pheromone scaffold needed for gallocin activation and uncovered unique interactions between two streptococcal QS signals that warrant further study. 
    more » « less
  3. Johnson, Karyn N. (Ed.)
    ABSTRACT A pervasive pest of stored leguminous products, the bean beetle Callosobruchus maculatus (Coleoptera: Chrysomelidae) associates with a simple bacterial community during adulthood. Despite its economic importance, little is known about the compositional stability, heritability, localization, and metabolic potential of the bacterial symbionts of C. maculatus . In this study, we applied community profiling using 16S rRNA gene sequencing to reveal a highly conserved bacterial assembly shared between larvae and adults. Dominated by Firmicutes and Proteobacteria , this community is localized extracellularly along the epithelial lining of the bean beetle’s digestive tract. Our analysis revealed that only one species, Staphylococcus gallinarum (phylum Firmicutes ), is shared across all developmental stages. Isolation and whole-genome sequencing of S. gallinarum from the beetle gut yielded a circular chromosome (2.8 Mb) and one plasmid (45 kb). The strain encodes complete biosynthetic pathways for the production of B vitamins and amino acids, including tyrosine, which is increasingly recognized as an important symbiont-supplemented precursor for cuticle biosynthesis in beetles. A carbohydrate-active enzyme search revealed that the genome codes for a number of digestive enzymes, reflecting the nutritional ecology of C. maculatus . The ontogenic conservation of the gut microbiota in the bean beetle, featuring a “core” community composed of S. gallinarum , may be indicative of an adaptive role for the host. In clarifying symbiont localization and metabolic potential, we further our understanding and study of a costly pest of stored products. IMPORTANCE From supplementing essential nutrients to detoxifying plant secondary metabolites and insecticides, bacterial symbionts are a key source of adaptations for herbivorous insect pests. Despite the pervasiveness and geographical range of the bean beetle Callosobruchus maculatus , the role of microbial symbioses in its natural history remains understudied. Here, we demonstrate that the bean beetle harbors a simple gut bacterial community that is stable throughout development. This community localizes along the insect’s digestive tract and is largely dominated by Staphylococcus gallinarum . In elucidating symbiont metabolic potential, we highlight its possible adaptive significance for a widespread agricultural pest. 
    more » « less
  4. This data set for the manuscript entitled "Design of Peptides that Fold and Self-Assemble on Graphite" includes all files needed to run and analyze the simulations described in the this manuscript in the molecular dynamics software NAMD, as well as the output of the simulations. The files are organized into directories corresponding to the figures of the main text and supporting information. They include molecular model structure files (NAMD psf or Amber prmtop format), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD configuration files, Colvars configuration files, NAMD log files, and NAMD output including restart files (in binary NAMD format) and trajectories in dcd format (downsampled to 10 ns per frame). Analysis is controlled by shell scripts (Bash-compatible) that call VMD Tcl scripts or python scripts. These scripts and their output are also included.

    Version: 2.0

    Changes versus version 1.0 are the addition of the free energy of folding, adsorption, and pairing calculations (Sim_Figure-7) and shifting of the figure numbers to accommodate this addition.


    Conventions Used in These Files
    ===============================

    Structure Files
    ----------------
    - graph_*.psf or sol_*.psf (original NAMD (XPLOR?) format psf file including atom details (type, charge, mass), as well as definitions of bonds, angles, dihedrals, and impropers for each dipeptide.)

    - graph_*.pdb or sol_*.pdb (initial coordinates before equilibration)
    - repart_*.psf (same as the above psf files, but the masses of non-water hydrogen atoms have been repartitioned by VMD script repartitionMass.tcl)
    - freeTop_*.pdb (same as the above pdb files, but the carbons of the lower graphene layer have been placed at a single z value and marked for restraints in NAMD)
    - amber_*.prmtop (combined topology and parameter files for Amber force field simulations)
    - repart_amber_*.prmtop (same as the above prmtop files, but the masses of non-water hydrogen atoms have been repartitioned by ParmEd)

    Force Field Parameters
    ----------------------
    CHARMM format parameter files:
    - par_all36m_prot.prm (CHARMM36m FF for proteins)
    - par_all36_cgenff_no_nbfix.prm (CGenFF v4.4 for graphene) The NBFIX parameters are commented out since they are only needed for aromatic halogens and we use only the CG2R61 type for graphene.
    - toppar_water_ions_prot_cgenff.str (CHARMM water and ions with NBFIX parameters needed for protein and CGenFF included and others commented out)

    Template NAMD Configuration Files
    ---------------------------------
    These contain the most commonly used simulation parameters. They are called by the other NAMD configuration files (which are in the namd/ subdirectory):
    - template_min.namd (minimization)
    - template_eq.namd (NPT equilibration with lower graphene fixed)
    - template_abf.namd (for adaptive biasing force)

    Minimization
    -------------
    - namd/min_*.0.namd

    Equilibration
    -------------
    - namd/eq_*.0.namd

    Adaptive biasing force calculations
    -----------------------------------
    - namd/eabfZRest7_graph_chp1404.0.namd
    - namd/eabfZRest7_graph_chp1404.1.namd (continuation of eabfZRest7_graph_chp1404.0.namd)

    Log Files
    ---------
    For each NAMD configuration file given in the last two sections, there is a log file with the same prefix, which gives the text output of NAMD. For instance, the output of namd/eabfZRest7_graph_chp1404.0.namd is eabfZRest7_graph_chp1404.0.log.

    Simulation Output
    -----------------
    The simulation output files (which match the names of the NAMD configuration files) are in the output/ directory. Files with the extensions .coor, .vel, and .xsc are coordinates in NAMD binary format, velocities in NAMD binary format, and extended system information (including cell size) in text format. Files with the extension .dcd give the trajectory of the atomic coorinates over time (and also include system cell information). Due to storage limitations, large DCD files have been omitted or replaced with new DCD files having the prefix stride50_ including only every 50 frames. The time between frames in these files is 50 * 50000 steps/frame * 4 fs/step = 10 ns. The system cell trajectory is also included for the NPT runs are output/eq_*.xst.

    Scripts
    -------
    Files with the .sh extension can be found throughout. These usually provide the highest level control for submission of simulations and analysis. Look to these as a guide to what is happening. If there are scripts with step1_*.sh and step2_*.sh, they are intended to be run in order, with step1_*.sh first.


    CONTENTS
    ========

    The directory contents are as follows. The directories Sim_Figure-1 and Sim_Figure-8 include README.txt files that describe the files and naming conventions used throughout this data set.

    Sim_Figure-1: Simulations of N-acetylated C-amidated amino acids (Ac-X-NHMe) at the graphite–water interface.

    Sim_Figure-2: Simulations of different peptide designs (including acyclic, disulfide cyclized, and N-to-C cyclized) at the graphite–water interface.

    Sim_Figure-3: MM-GBSA calculations of different peptide sequences for a folded conformation and 5 misfolded/unfolded conformations.

    Sim_Figure-4: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.

    Sim_Figure-5: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 295 K.

    Sim_Figure-5_replica: Temperature replica exchange molecular dynamics simulations for the peptide cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) with 20 replicas for temperatures from 295 to 454 K.

    Sim_Figure-6: Simulation of the peptide molecule cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) in free solution (no graphite).

    Sim_Figure-7: Free energy calculations for folding, adsorption, and pairing for the peptide CHP1404 (sequence: cyc(GTGSGTG-GPGG-GCGTGTG-SGPG)). For folding, we calculate the PMF as function of RMSD by replica-exchange umbrella sampling (in the subdirectory Folding_CHP1404_Graphene/). We make the same calculation in solution, which required 3 seperate replica-exchange umbrella sampling calculations (in the subdirectory Folding_CHP1404_Solution/). Both PMF of RMSD calculations for the scrambled peptide are in Folding_scram1404/. For adsorption, calculation of the PMF for the orientational restraints and the calculation of the PMF along z (the distance between the graphene sheet and the center of mass of the peptide) are in Adsorption_CHP1404/ and Adsorption_scram1404/. The actual calculation of the free energy is done by a shell script ("doRestraintEnergyError.sh") in the 1_free_energy/ subsubdirectory. Processing of the PMFs must be done first in the 0_pmf/ subsubdirectory. Finally, files for free energy calculations of pair formation for CHP1404 are found in the Pair/ subdirectory.

    Sim_Figure-8: Simulation of four peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) where the peptides are far above the graphene–water interface in the initial configuration.

    Sim_Figure-9: Two replicates of a simulation of nine peptide molecules with the sequence cyc(GTGSGTG-GPGG-GCGTGTG-SGPG) at the graphite–water interface at 370 K.

    Sim_Figure-9_scrambled: Two replicates of a simulation of nine peptide molecules with the control sequence cyc(GGTPTTGGGGGGSGGPSGTGGC) at the graphite–water interface at 370 K.

    Sim_Figure-10: Adaptive biasing for calculation of the free energy of the folded peptide as a function of the angle between its long axis and the zigzag directions of the underlying graphene sheet.

     

    This material is based upon work supported by the US National Science Foundation under grant no. DMR-1945589. A majority of the computing for this project was performed on the Beocat Research Cluster at Kansas State University, which is funded in part by NSF grants CHE-1726332, CNS-1006860, EPS-1006860, and EPS-0919443. This work used the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1548562, through allocation BIO200030. 
    more » « less
  5. Abstract Background The explosive radiation and diversification of the advanced snakes (superfamily Colubroidea) was associated with changes in all aspects of the shared venom system. Morphological changes included the partitioning of the mixed ancestral glands into two discrete glands devoted for production of venom or mucous respectively, as well as changes in the location, size and structural elements of the venom-delivering teeth. Evidence also exists for homology among venom gland toxins expressed across the advanced snakes. However, despite the evolutionary novelty of snake venoms, in-depth toxin molecular evolutionary history reconstructions have been mostly limited to those types present in only two front-fanged snake families, Elapidae and Viperidae. To have a broader understanding of toxins shared among extant snakes, here we first sequenced the transcriptomes of eight taxonomically diverse rear-fanged species and four key viperid species and analysed major toxin types shared across the advanced snakes. Results Transcriptomes were constructed for the following families and species: Colubridae - Helicops leopardinus , Heterodon nasicus , Rhabdophis subminiatus ; Homalopsidae – Homalopsis buccata ; Lamprophiidae - Malpolon monspessulanus , Psammophis schokari , Psammophis subtaeniatus , Rhamphiophis oxyrhynchus ; and Viperidae – Bitis atropos , Pseudocerastes urarachnoides , Tropidolaeumus subannulatus , Vipera transcaucasiana . These sequences were combined with those from available databases of other species in order to facilitate a robust reconstruction of the molecular evolutionary history of the key toxin classes present in the venom of the last common ancestor of the advanced snakes, and thus present across the full diversity of colubroid snake venoms. In addition to differential rates of evolution in toxin classes between the snake lineages, these analyses revealed multiple instances of previously unknown instances of structural and functional convergences. Structural convergences included: the evolution of new cysteines to form heteromeric complexes, such as within kunitz peptides (the beta-bungarotoxin trait evolving on at least two occasions) and within SVMP enzymes (the P-IIId trait evolving on at least three occasions); and the C-terminal tail evolving on two separate occasions within the C-type natriuretic peptides, to create structural and functional analogues of the ANP/BNP tailed condition. Also shown was that the de novo evolution of new post-translationally liberated toxin families within the natriuretic peptide gene propeptide region occurred on at least five occasions, with novel functions ranging from induction of hypotension to post-synaptic neurotoxicity. Functional convergences included the following: multiple occasions of SVMP neofunctionalised in procoagulant venoms into activators of the clotting factors prothrombin and Factor X; multiple instances in procoagulant venoms where kunitz peptides were neofunctionalised into inhibitors of the clot destroying enzyme plasmin, thereby prolonging the half-life of the clots formed by the clotting activating enzymatic toxins; and multiple occasions of kunitz peptides neofunctionalised into neurotoxins acting on presynaptic targets, including twice just within Bungarus venoms. Conclusions We found novel convergences in both structural and functional evolution of snake toxins. These results provide a detailed roadmap for future work to elucidate predator–prey evolutionary arms races, ascertain differential clinical pathologies, as well as documenting rich biodiscovery resources for lead compounds in the drug design and discovery pipeline. 
    more » « less