Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Docking”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,153 records · Page 64

Detecting and Characterizing Fracture Zones Using a Convolutional Neural Network

This project directly supports the Geothermal Technologies Office (GTO) objectives outlined in the Multi-Year Program Plan (MYPP) by advancing two key research areas: “Exploration and Characterization” and “Data, Modeling, and Analysis.” This project has successfully demonstrated a pre-drilling ability to image and characterize the distribution and connectivity of subsurface faults and fractures, key parameters for identifying permeable pathways that enable geothermal fluids to circulate and produce energy. Specifically, we developed and implemented innovative machine learning methodologies to enhance geothermal exploration. Large-scale faults were detected using a Convolutional Neural Network (CNN), while small-scale fractures were characterized using a novel Double-Beam Neural Network (DBNN). These tools have proven both technically effective and cost-efficient by reducing reliance on expensive exploratory drilling. Through collaboration with our geothermal industry partner, this research has significantly advanced techniques for identifying hidden geothermal systems and extending the productive lifespan of existing geothermal fields. We applied our methods to two geothermal fields—Soda Lake (Nevada) and Lightning Dock (New Mexico)—to identify shallow steam-charged fracture zones and characterize deep faults at depths of 1.5-2 km. The steam zone identified at the Soda Lake geothermal field showed excellent agreement with prior drilling data, validating the effectiveness of our approaches. In addition, the analysis revealed three new prospective drilling targets for further development and verification. The outcomes of this project improve our scientific understanding of geothermal reservoir behavior, enhance exploration efficiency, extend the economic life of existing geothermal plants. Ultimately, these advancements contribute to GTO’s goal of achieving more sustainable, affordable, and data-driven geothermal energy development across the United States.

15 GEOTHERMAL ENERGY↗

Where do the fish go in winter? A year of observing seasonal changes in Sequim Bay’s nearshore fish community

The use of environmental DNA (eDNA) sampling has been proposed as a complementary method to monitor fish species in marine environments. eDNA offers a non-invasive, cost-effective, and scalable way to detect aquatic species. It is appealing in environments where traditional methods are limited by access or visibility, especially in complex or sensitive habitats such as tidal channels and other marine energy sites. Before eDNA can be fully relied upon, we must verify its accuracy against established methods, like underwater photography. In this study, we collected eDNA samples and concurrently deployed a 360-degree camera near the floating dock of PNNL-Sequim in the tidal channel of Sequim Bay once a month for twelve consecutive months. During deployments, the camera remained on the seafloor for several hours overlapping slack tide and captured time-lapse photographs at ten second intervals. Counts and identifications of fish species observed in the images were used to calculate monthly Shannon diversity and Pielou evenness indices. These values were compared across the months using a Kruskal-Wallis test paired with a Conover-Iman post-hoc test. Cliff’s Delta was also calculated to quantify the effect size of the monthly differences. No fish were observed from December through April, likely due to seasonal behavior changes within the local fish community. Fish returned to the shoreline in May, with the greatest diversity and evenness recorded in September. These findings reveal substantial seasonal variation in nearshore fish communities. Many monthly comparisons were found to contain statistically significant differences within the diversity and evenness, and even more were found to have large effect sizes, signifying large ecological changes throughout the seasons. The absence of fish observed during the winter months is a key outcome of this survey that will hopefully be reflected in the eDNA results still to come, which would help validate the eDNA approach.

59 BASIC BIOLOGICAL SCIENCES↗

Models and Processes to Extract Drug-like Molecules From Natural Language Text

Researchers worldwide are seeking to repurpose existing drugs or discover new drugs to counter the disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). A promising source of candidates for such studies is molecules that have been reported in the scientific literature to be drug-like in the context of viral research. However, this literature is too large for human review and features unusual vocabularies for which existing named entity recognition (NER) models are ineffective. We report here on a project that leverages both human and artificial intelligence to detect references to such molecules in free text. We present 1) a iterative model-in-the-loop method that makes judicious use of scarce human expertise in generating training data for a NER model, and 2) the application and evaluation of this method to the problem of identifying drug-like molecules in the COVID-19 Open Research Dataset Challenge (CORD-19) corpus of 198,875 papers. We show that by repeatedly presenting human labelers only with samples for which an evolving NER model is uncertain, our human-machine hybrid pipeline requires only modest amounts of non-expert human labeling time (tens of hours to label 1778 samples) to generate an NER model with an F-1 score of 80.5%—on par with that of non-expert humans—and when applied to CORD’19, identifies 10,912 putative drug-like molecules. This enriched the computational screening team’s targets by 3,591 molecules, of which 18 ranked in the top 0.1% of all 6.6 million molecules screened for docking against the 3CLPro protein.

60 APPLIED LIFE SCIENCES↗

Discovery of Small-Molecule Inhibitors of SARS-CoV-2 Proteins Using a Computational and Experimental Pipeline

A rapid response is necessary to contain emergent biological outbreaks before they can become pandemics. The novel coronavirus (SARS-CoV-2) that causes COVID-19 was first reported in December of 2019 in Wuhan, China and reached most corners of the globe in less than two months. In just over a year since the initial infections, COVID-19 infected almost 100 million people worldwide. Although similar to SARS-CoV and MERS-CoV, SARS-CoV-2 has resisted treatments that are effective against other coronaviruses. Crystal structures of two SARS-CoV-2 proteins, spike protein and main protease, have been reported and can serve as targets for studies in neutralizing this threat. We have employed molecular docking, molecular dynamics simulations, and machine learning to identify from a library of 26 million molecules possible candidate compounds that may attenuate or neutralize the effects of this virus. The viability of selected candidate compounds against SARS-CoV-2 was determined experimentally by biolayer interferometry and FRET-based activity protein assays along with virus-based assays. In the pseudovirus assay, imatinib and lapatinib had IC 50 values below 10 μM, while candesartan cilexetil had an IC 50 value of approximately 67 µM against M pro in a FRET-based activity assay. Comparatively, candesartan cilexetil had the highest selectivity index of all compounds tested as its half-maximal cytotoxicity concentration 50 (CC 50 ) value was the only one greater than the limit of the assay (>100 μM).

59 BASIC BIOLOGICAL SCIENCES↗

Characterization of a Haematococcus pluvialis Diacylglycerol Acyltransferase 1 and Its Potential in Unsaturated Fatty Acid-Rich Triacylglycerol Production

The unicellular green alga Haematococcus pluvialis has been recognized as an industry strain to produce simultaneously esterified astaxanthin (EAST) and triacylglycerol (TAG) under stress induction. It is necessary to identify the key enzymes involving in synergistic accumulation of EAST and TAG in H. pluvialis . In this study, a novel diacylglycerol acyltransferase 1 was systematically characterized by in vivo and in silico assays. The upregulated expression of HpDGAT1 gene was positively associated with the significant increase of TAG and EAST contents under stress conditions. Functional complementation by overexpressing HpDGAT1 in a TAG-deficient yeast strain H1246 revealed that HpDGAT1 could restore TAG biosynthesis and exhibited a high substrate preference for monounsaturated fatty acyl-CoAs (MUFAs) and polyunsaturated fatty acyl-CoAs (PUFAs). Notably, heterogeneous expression of HpDGAT1 in Chlamydomonas reinhardtii and Arabidopsis thaliana resulted in a significant enhancement of total oils and concurrently a high accumulation of MUFAs- and PUFAs-rich TAGs. Furthermore, molecular docking analysis indicated that HpDGAT1 contained AST-binding sites. These findings evidence a possible dual-function role for HpDGAT1 involving in TAG and EAST synthesis, demonstrating that it is a potential target gene to enrich AST accumulation in this alga and to design oil production in both commercial algae and oil crops.

Cui, Hongli↗

Functional Characterization and Structural Insights Into Stereoselectivity of Pulegone Reductase in Menthol Biosynthesis

Monoterpenoids are the main components of plant essential oils and the active components of some traditional Chinese medicinal herbs like Mentha haplocalyx Briq ., Nepeta tenuifolia Briq ., Perilla frutescens (L.) Britt and Pogostemin cablin (Blanco) Benth. Pulegone reductase is the key enzyme in the biosynthesis of menthol and is required for the stereoselective reduction of the Δ 2,8 double bond of pulegone to produce the major intermediate menthone, thus determining the stereochemistry of menthol. However, the structural basis and mechanism underlying the stereoselectivity of pulegone reductase remain poorly understood. In this study, we characterized a novel (-)-pulegone reductase from Nepeta tenuifolia ( Nt PR), which can catalyze (-)-pulegone to (+)-menthone and (-)-isomenthone through our RNA-seq, bioinformatic analysis in combination with in vitro enzyme activity assay, and determined the structure of (+)-pulegone reductase from ( Mp PR) by using X-ray crystallography, molecular modeling and docking, site-directed mutagenesis, molecular dynamics simulations, and biochemical analysis. We identified and validated the critical residues in the crystal structure of Mp PR involved in the binding of the substrate pulegone. We also further identified that residues Leu56, Val282, and Val284 determine the stereoselectivity of the substrate pulegone, and mainly contributes to the product stereoselectivity. This work not only provides a starting point for the understanding of stereoselectivity of pulegone reductases, but also offers a basis for the engineering of menthone/menthol biosynthetic enzymes to achieve high-titer, industrial-scale production of enantiomerically pure products.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗

Indoline-6-Sulfonamide Inhibitors of the Bacterial Enzyme DapE

Inhibitors of the bacterial enzyme dapE-encoded N-succinyl-l,l-diaminopimelic acid desuccinylase (DapE; EC 3.5.1.18) hold promise as antibiotics with a new mechanism of action. Herein we describe the discovery of a new series of indoline sulfonamide DapE inhibitors from a high-throughput screen and the synthesis of a series of analogs. Inhibitory potency was measured by a ninhydrin-based DapE assay recently developed by our group. Molecular docking experiments suggest active site binding with the sulfonamide acting as a zinc-binding group (ZBG).

59 BASIC BIOLOGICAL SCIENCES↗

Improving Protein–Ligand Interaction Modeling with cryo-EM Data, Templates, and Deep Learning in 2021 Ligand Model Challenge

Elucidating protein–ligand interaction is crucial for studying the function of proteins and compounds in an organism and critical for drug discovery and design. The problem of protein–ligand interaction is traditionally tackled by molecular docking and simulation, which is based on physical forces and statistical potentials and cannot effectively leverage cryo-EM data and existing protein structural information in the protein–ligand modeling process. In this work, we developed a deep learning bioinformatics pipeline (DeepProLigand) to predict protein–ligand interactions from cryo-EM density maps of proteins and ligands. DeepProLigand first uses a deep learning method to predict the structure of proteins from cryo-EM maps, which is averaged with a reference (template) structure of the proteins to produce a combined structure to add ligands. The ligands are then identified and added into the structure to generate a protein–ligand complex structure, which is further refined. The method based on the deep learning prediction and template-based modeling was blindly tested in the 2021 EMDataResource Ligand Challenge and was ranked first in fitting ligands to cryo-EM density maps. These results demonstrate that the deep learning bioinformatics approach is a promising direction for modeling protein–ligand interactions on cryo-EM data using prior structural information.

59 BASIC BIOLOGICAL SCIENCES↗

Molecular Dynamics Simulations of a Catalytic Multivalent Peptide–Nanoparticle Complex

Molecular modeling of a supramolecular catalytic system is conducted resulting from the assembling between a small peptide and the surface of cationic self-assembled monolayers on gold nanoparticles, through a multiscale iterative approach including atomistic force field development, flexible docking with Brownian Dynamics and µs-long Molecular Dynamics simulations. Self-assembly is a prerequisite for the catalysis, since the catalytic peptides do not display any activity in the absence of the gold nanocluster. Atomistic simulations reveal details of the association dynamics as regulated by defined conformational changes of the peptide due to peptide length and sequence. Our results show the importance of a rational design of the peptide to enhance the catalytic activity of peptide–nanoparticle conjugates and present a viable computational approach toward the design of enzyme mimics having a complex structure–function relationship, for technological and nanomedical applications.

59 BASIC BIOLOGICAL SCIENCES↗

Enhancing Paraoxon Binding to Organophosphorus Hydrolase Active Site

Organophosphorus hydrolase (OPH) is a metalloenzyme that can hydrolyze organophosphorus agents resulting in products that are generally of reduced toxicity. The best OPH substrate found to date is diethyl p-nitrophenyl phosphate (paraoxon). Most structural and kinetic studies assume that the binding orientation of paraoxon is identical to that of diethyl 4-methylbenzylphosphonate, which is the only substrate analog co-crystallized with OPH. In the current work, we used a combined docking and molecular dynamics (MD) approach to predict the likely binding mode of paraoxon. Then, we used the predicted binding mode to run MD simulations on the wild type (WT) OPH complexed with paraoxon, and OPH mutants complexed with paraoxon. Additionally, we identified three hot-spot residues (D253, H254, and I255) involved in the stability of the OPH active site. We then experimentally assayed single and double mutants involving these residues for paraoxon binding affinity. The binding free energy calculations and the experimental kinetics of the reactions between each OPH mutant and paraoxon show that mutated forms D253E, D253E-H254R, and D253E-I255G exhibit enhanced substrate binding affinity over WT OPH. Interestingly, our experimental results show that the substrate binding affinity of the double mutant D253E-H254R increased by 19-fold compared to WT OPH.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Revising the Role of Cortical Cytoskeleton during Secretion: Actin and Myosin XI Function in Vesicle Tethering

In plants, secretion of cell wall components and membrane proteins plays a fundamental role in growth and development as well as survival in diverse environments. Exocytosis, as the last step of the secretory trafficking pathway, is a highly ordered and precisely controlled process involving tethering, docking, and fusion of vesicles at the plasma membrane (PM) for cargo delivery. Although the exocytic process and machinery are well characterized in yeast and animal models, the molecular players and specific molecular events that underpin late stages of exocytosis in plant cells remain largely unknown. Here, by using the delivery of functional, fluorescent-tagged cellulose synthase (CESA) complexes (CSCs) to the PM as a model system for secretion, as well as single-particle tracking in living cells, we describe a quantitative approach for measuring the frequency of vesicle tethering events. Genetic and pharmacological inhibition of cytoskeletal function, reveal that the initial vesicle tethering step of exocytosis is dependent on actin and myosin XI. In contrast, treatments with the microtubule inhibitor, oryzalin, did not significantly affect vesicle tethering or fusion during CSC exocytosis but caused a minor increase in transient or aborted tethering events. With data from this new quantitative approach and improved spatiotemporal resolution of single particle events during secretion, we generate a revised model for the role of the cortical cytoskeleton in CSC trafficking.

Zhang, Weiwei↗

Characterization of Interactions between CTX-M-15 and Clavulanic Acid, Desfuroylceftiofur, Ceftiofur, Ampicillin, and Nitrocefin

Cefotaximase-Munich (CTX-M) extended-spectrum beta-lactamases (ESBLs) are commonly associated with Gram-negative, hospital-acquired infections worldwide. Several beta-lactamase inhibitors, such as clavulanate, are used to inhibit the activity of these enzymes. To understand the mechanism of CTX-M-15 activity, we have determined the crystal structures of CTX-M-15 in complex with two specific classes of beta-lactam compounds, desfuroylceftiofur (DFC) and ampicillin, and an inhibitor, clavulanic acid. The crystal structures revealed that Ser70 and five other residues (Lys73, Tyr105, Glu166, Ser130, and Ser237) participate in catalysis and binding of those compounds. Based on analysis of steady-state kinetics, thermodynamic data, and molecular docking to both wild-type and S70A mutant structures, we determined that CTX-M-15 has a similar affinity for all beta-lactam compounds (ceftiofur, nitrocefin, DFC, and ampicillin), but with lower affinity for clavulanic acid. A catalytic mechanism for tested β-lactams and two-step inhibition mechanism of clavulanic acid were proposed. CTX-M-15 showed a higher activity toward DFC and nitrocefin, but significantly lower activity toward ampicillin and ceftiofur. The interaction between CTX-M-15 and both ampicillin and ceftiofur displayed a higher entropic but lower enthalpic effect, compared with DFC and nitrocefin. DFC, a metabolite of ceftiofur, displayed lower entropy and higher enthalpy than ceftiofur. This finding suggests that compounds containing amine moiety (e.g., ampicillin) and the furfural moiety (e.g., ceftiofur) could hinder the hydrolytic activity of CTX-M-15.

59 BASIC BIOLOGICAL SCIENCES↗

Crystal Structure of CYP3A4 Complexed with Fluorol Identifies the Substrate Access Channel as a High-Affinity Ligand Binding Site

Cytochrome P450 3A4 (CYP3A4) is a major human drug-metabolizing enzyme, notoriously known for its extreme substrate promiscuity, allosteric behavior, and implications in drug–drug interactions. Despite extensive investigations, the mechanism of ligand binding to CYP3A4 is not fully understood. We determined the crystal structure of CYP3A4 complexed with fluorol, a small fluorescent dye that can undergo hydroxylation. In the structure, fluorol associates to the substrate channel, well suited for the binding of planar polyaromatic molecules bearing polar groups, through which stabilizing H-bonds with the polar channel residues, such as Thr224 and Arg372, can be established. Mutagenesis, spectral, kinetic, and functional data confirmed the involvement but not strict requirement of Thr224 for the association of fluorol. Collectively, our data identify the substrate channel as a high-affinity ligand binding site and support the notion that hydrophobic ligands first dock to the nearby peripheral surface, before migrating to the channel and, subsequently, into the active site.

59 BASIC BIOLOGICAL SCIENCES↗

Structural Similarities and Overlapping Activities among Dihydroflavonol 4-Reductase, Flavanone 4-Reductase, and Anthocyanidin Reductase Offer Metabolic Flexibility in the Flavonoid Pathway

Flavonoids are potent antioxidants that play a role in defense against pathogens, UV-radiation, and the detoxification of reactive oxygen species. Dihydroflavonol 4-reductase (DFR) and flavanone 4-reductase (FNR) reduce dihydroflavonols and flavanones, respectively, using NAD(P)H to produce flavan-(3)-4-(di)ols in flavonoid biosynthesis. Anthocyanidin reductase (ANR) reduces anthocyanidins to flavan-3-ols. In addition to their sequences, the 3D structures of recombinant DFR, FNR and ANR from sorghum and switchgrass showed a high level of similarity. The catalytic mechanism, substrate-specificity and key residues of three reductases were deduced from crystal structures, site-directed mutagenesis, molecular docking, kinetics, and thermodynamic ana-lyses. Although DFR displayed its highest activity against dihydroflavonols, it also showed activity against flavanones and anthocyanidins. It was inhibited by the flavonol quercetin and high concentrations of dihydroflavonols/flavonones. SbFNR1 and SbFNR2 did not show any activity against dihydroflavonols. However, SbFNR1 displayed activity against flavanones and ANR activity against two anthocyanidins, cyanidin and pelargonidin. Therefore, SbFNR1 and SbFNR2 could be specific ANR isozymes without delphinidin activity. Sorghum has high concentrations of 3-deoxyanthocyanidins in vivo, supporting the observed high activity of SbDFR against flavonols. Mining of expression data indicated substantial induction of these three reductase genes in both switchgrass and sorghum in response to biotic stress. Key signature sequences for proper DFR/ANR classification are proposed and could form the basis for future metabolic engineering of flavonoid metabolism.

3-deoxyanthocyanidin↗

Structural and Interactional Analysis of the Flavonoid Pathway Proteins: Chalcone Synthase, Chalcone Isomerase and Chalcone Isomerase-like Protein

Chalcone synthase (CHS) and chalcone isomerase (CHI) catalyze the first two committed steps of the flavonoid pathway that plays a pivotal role in the growth and reproduction of land plants, including UV protection, pigmentation, symbiotic nitrogen fixation, and pathogen resistance. Based on the obtained X-ray crystal structures of CHS, CHI, and chalcone isomerase-like protein (CHIL) from the same monocotyledon, Panicum virgatum, along with the results of the steady-state kinetics, spectroscopic/thermodynamic analyses, intermolecular interactions, and their effect on each catalytic step are proposed. In addition, PvCHI’s unique activity for both naringenin chalcone and isoliquiritigenin was analyzed, and the observed hierarchical activity for those type-I and -II substrates was explained with the intrinsic characteristics of the enzyme and two substrates. The structure of PvCHS complexed with naringenin supports uncompetitive inhibition. PvCHS displays intrinsic catalytic promiscuity, evident from the formation of p-coumaroyltriacetic acid lactone (CTAL) in addition to naringenin chalcone. In the presence of PvCHIL, conversion of p-coumaroyl-CoA to naringenin through PvCHS and PvCHI displayed ~400-fold increased Vmax with reduced formation of CTAL by 70%. Supporting this model, molecular docking, ITC (Isothermal Titration Calorimetry), and FRET (Fluorescence Resonance Energy Transfer) indicated that both PvCHI and PvCHIL interact with PvCHS in a non-competitive manner, indicating the plausible allosteric effect of naringenin on CHS. Significantly, the presence of naringenin increased the affinity between PvCHS and PvCHIL, whereas naringenin chalcone decreased the affinity, indicating a plausible feedback mechanism to minimize spontaneous incorrect stereoisomers. These are the first findings from a three-body system from the same species, indicating the importance of the macromolecular assembly of CHS-CHI-CHIL in determining the amount and type of flavonoids produced in plant cells.

59 BASIC BIOLOGICAL SCIENCES↗

Prospecting Biotechnologically-Relevant Monooxygenases from Cold Sediment Metagenomes: An In Silico Approach

The goal of this work was to identify sequences encoding monooxygenase biocatalysts with novel features by in silico mining an assembled metagenomic dataset of polar and subpolar marine sediments. The targeted enzyme sequences were Baeyer–Villiger and bacterial cytochrome P450 monooxygenases (CYP153). These enzymes have wide-ranging applications, from the synthesis of steroids, antibiotics, mycotoxins and pheromones to the synthesis of monomers for polymerization and anticancer precursors, due to their extraordinary enantio-, regio-, and chemo- selectivity that are valuable features for organic synthesis. Phylogenetic analyses were used to select the most divergent sequences affiliated to these enzyme families among the 264 putative monooxygenases recovered from the ~14 million protein-coding sequences in the assembled metagenome dataset. Three-dimensional structure modeling and docking analysis suggested features useful in biotechnological applications in five metagenomic sequences, such as wide substrate range, novel substrate specificity or regioselectivity. Further analysis revealed structural features associated with psychrophilic enzymes, such as broader substrate accessibility, larger catalytic pockets or low domain interactions, suggesting that they could be applied in biooxidations at room or low temperatures, saving costs inherent to energy consumption. This work allowed the identification of putative enzyme candidates with promising features from metagenomes, providing a suitable starting point for further developments.

59 BASIC BIOLOGICAL SCIENCES↗

Snails In Silico: A Review of Computational Studies on the Conopeptides

Marine cone snails are carnivorous gastropods that use peptide toxins called conopeptides both as a defense mechanism and as a means to immobilize and kill their prey. These peptide toxins exhibit a large chemical diversity that enables exquisite specificity and potency for target receptor proteins. This diversity arises in terms of variations both in amino acid sequence and length, and in posttranslational modifications, particularly the formation of multiple disulfide linkages. Most of the functionally characterized conopeptides target ion channels of animal nervous systems, which has led to research on their therapeutic applications. Many facets of the underlying molecular mechanisms responsible for the specificity and virulence of conopeptides, however, remain poorly understood. In this review, we will explore the chemical diversity of conopeptides from a computational perspective. First, we discuss current approaches used for classifying conopeptides. Next, we review different computational strategies that have been applied to understanding and predicting their structure and function, from machine learning techniques for predictive classification to docking studies and molecular dynamics simulations for molecular-level understanding. We then review recent novel computational approaches for rapid high-throughput screening and chemical design of conopeptides for particular applications. We close with an assessment of the state of the field, emphasizing important questions for future lines of inquiry.

59 BASIC BIOLOGICAL SCIENCES↗

Evolutionary Diversification of Host-Targeted Bartonella Effectors Proteins Derived from a Conserved FicTA Toxin-Antitoxin Module

Proteins containing a FIC domain catalyze AMPylation and other post-translational modifications (PTMs). In bacteria, they are typically part of FicTA toxin-antitoxin modules that control conserved biochemical processes such as topoisomerase activity, but they have also repeatedly diversified into host-targeted virulence factors. Among these, Bartonella effector proteins (Beps) comprise a particularly diverse ensemble of FIC domains that subvert various host cellular functions. However, no comprehensive comparative analysis has been performed to infer molecular mechanisms underlying the biochemical and functional diversification of FIC domains in the vast Bep family. Here, we used X-ray crystallography, structural modelling, and phylogenetic analyses to unravel the expansion and diversification of Bep repertoires that evolved in parallel in three Bartonella lineages from a single ancestral FicTA toxin-antitoxin module. Our analysis is based on 99 non-redundant Bep sequences and nine crystal structures. Inferred from the conservation of the FIC signature motif that comprises the catalytic histidine and residues involved in substrate binding, about half of them represent AMP transferases. A quarter of Beps show a glutamate in a strategic position in the putative substrate binding pocket that would interfere with triphosphate-nucleotide binding but may allow binding of an AMPylated target for deAMPylation or another substrate to catalyze a distinct PTM. The β-hairpin flap that registers the modifiable target segment to the active site exhibits remarkable structural variability. The corresponding sequences form few well-defined groups that may recognize distinct target proteins. The binding of Beps to promiscuous FicA antitoxins is well conserved, indicating a role of the antitoxin to inhibit enzymatic activity or to serve as a chaperone for the FIC domain before translocation of the Bep into host cells. Taken together, our analysis indicates a remarkable functional plasticity of Beps that is mostly brought about by structural changes in the substrate pocket and the target dock. These findings may guide future structure–function analyses of the highly versatile FIC domains.

59 BASIC BIOLOGICAL SCIENCES↗