Structural Reconstruction of E. coli Ubi Metabolon Using an AlphaFold2-Based Computational Framework.

Ubiquinone (UQ) is a redox polyisoprenoid lipid found in the membranes of bacteria and eukaryotes that has important roles, notably one in respiratory metabolism, which sustains cellular bioenergetics. In Escherichia coli, several steps of the UQ biosynthesis take place in the cytosol. To perform these reactions, a supramolecular assembly called Ubi metabolon is involved. This latter is composed of seven proteins (UbiE, UbiG, UbiF, UbiH, UbiI, UbiJ, and UbiK), and its structural organization is unknown as well as its protein stoichiometry. In this study, a computational framework has been designed to predict the structure of this macromolecular assembly. In several successive steps, we explored the possible protein interactions as well as the protein stoichiometry, to finally obtain a structural organization of the complex. The use of AlphaFold2-based methods combined with evolutionary information enabled us to predict several models whose quality and confidence were further analyzed using different metrics and scores. Our work led to the identification of a "core assembly" that will guide functional and structural characterization of the Ubi metabolon.

Towards Molecular Understanding of the Functional Role of UbiJ-UbiK2 Complex in Ubiquinone Biosynthesis by Multiscale Molecular Modelling Studies.

Ubiquinone (UQ) is a polyisoprenoid lipid found in the membranes of bacteria and eukaryotes. UQ has important roles, notably in respiratory metabolisms which sustain cellular bioenergetics. Most steps of UQ biosynthesis take place in the cytosol of E. coli within a multiprotein complex called the Ubi metabolon, that contains five enzymes and two accessory proteins, UbiJ and UbiK. The SCP2 domain of UbiJ was proposed to bind the hydrophobic polyisoprenoid tail of UQ biosynthetic intermediates in the Ubi metabolon. How the newly synthesised UQ might be released in the membrane is currently unknown. In this paper, we focused on better understanding the role of the UbiJ-UbiK2 heterotrimer forming part of the metabolon. Given the difficulties to gain functional insights using biophysical techniques, we applied a multiscale molecular modelling approach to study the UbiJ-UbiK2 heterotrimer. Our data show that UbiJ-UbiK2 interacts closely with the membrane and suggests possible pathways to enable the release of UQ into the membrane. This study highlights the UbiJ-UbiK2 complex as the likely interface between the membrane and the enzymes of the Ubi metabolon and supports that the heterotrimer is key to the biosynthesis of UQ8 and its release into the membrane of E. coli.

Synthesis of α-l-Araf and ß-d-Galf series furanobiosides using mutants of a GH51 α-l-arabinofuranosidase.

The GH-51 α-l-arabinofuranosidase from Thermobacillus xylanilyticus (TxAbf) possesses versatile catalytic properties, displaying not only the ability to hydrolyze glycosidic linkages but also to synthesize furanobiosides in α-l-Araf and ß-d-Galf series. Herein, mutants are investigated to evaluate their ability to perform self-condensation, assessing both yield improvements and changes in regioselectivity. Overall yields of oligo-α-l-arabino- and oligo-ß-d-galactofuranosides were increased up to 4.8-fold compared to the wild-type enzyme. In depth characterization revealed that the mutants exhibit increased transfer rates and thus a hydrolysis/self-condensation ratio in favor of synthesis. The consequence of the substitution N216W is the creation of an additional binding subsite that provides the basis for an alternative acceptor substrate binding mode. As a result, mutants bearing N216W synthesize not only (1,2)-linked furanobiosides, but also (1,3)- and even (1,5)-linked furanobiosides. Since the self-condensation is under kinetic control, the yield of homo-disaccharides was maximized using higher substrate concentrations. In this way, the mutant R69H-N216W produced oligo-ß-d-galactofuranosides in > 70% yield. Overall, this study further demonstrates the potential usefulness of TxAbf mutants for glycosynthesis and shows how these might be used to synthesize biologically-relevant glycoconjugates.

A computational study of the molecular basis of antibiotic resistance in a DXR mutant.

Proteins of the independent mevalonate pathway for isoprenoid biosynthesis are important targets for the development of new antibacterial compounds as this pathway is present in most pathogenic organisms such as Mycobacterium tuberculosis, DPlasmodium falciparum and Escherichia coli, but is not present in mammalian species, including humans. Deoxy-D-xylulose 5-phosphate reductoisomerase (DXR) is an important target in this pathway and the most effective DXR inhibitor to date is fosmidomycin, which is used to treat malaria and, more recently, tuberculosis. Recently, Armstrong C. M. et al. showed that a mutant of DXR, S222T, induces a loss of the fosmidomycin inhibition efficiency, even though the bacteria culture is still viable and able to produce isoprenoids. As this represents a potential fosmidomycin-resistant mutation, it is important to understand the mechanism of this apparent mutation-induced resistance to fosmidomycin. Here, we used molecular dynamics simulations and Molecular Mechanics/Poisson Boltzmann Surface Area analysis to understand the structural and energetic basis of the resistance. Our results suggest that the point mutation results in changes to the structural dynamics of an active site loop that probably protects the active site and facilitates enzymatic reaction. From the simulation analysis, we also showed that the mutation results in changes in the interaction energy profiles in a way that can explain the observed activity of the mutant protein toward the natural inhibitor deoxy-D-xylulose 5-phosphate. These results should be taken into consideration in future efforts to develop new therapeutic antibiotic compounds that target DXR.

Synthesis and biological evaluation of phosphate isosters of fosmidomycin and analogs as inhibitors of Escherichia coli and Mycobacterium smegmatis 1-deoxyxylulose 5-phosphate reductoisomerases.

Hydroxamate analogs of fosfoxacin, the phosphate homolog of fosmidomycin, have been synthesized and their activity tested on Escherichia coli and Mycobacterium smegmatis DXRs. Except for compound 4b, the IC50 values of phosphate derivatives are approximately 10-fold higher than those of the corresponding phosphonates. Although their inhibitory activity on Escherichia coli DXR is less efficient than their phosphonate analogs, we report the ability of phosphate compounds to inhibit the growth of Escherichia coli. This work points out that the uptake of fosfoxacin and its analogs is taking place via the GlpT and UhpT transporters. As expected, these compounds are inefficient to inhibit the growth of M. smegmatis growth inhibition probably due to a lack of uptake.

Sequence-structure relationship study in all-α transmembrane proteins using an unsupervised learning approach.

Transmembrane proteins (TMPs) are major drug targets, but the knowledge of their precise topology structure remains highly limited compared with globular proteins. In spite of the difficulties in obtaining their structures, an important effort has been made these last years to increase their number from an experimental and computational point of view. In view of this emerging challenge, the development of computational methods to extract knowledge from these data is crucial for the better understanding of their functions and in improving the quality of structural models. Here, we revisit an efficient unsupervised learning procedure, called Hybrid Protein Model (HPM), which is applied to the analysis of transmembrane proteins belonging to the all-α structural class. HPM method is an original classification procedure that efficiently combines sequence and structure learning. The procedure was initially applied to the analysis of globular proteins. In the present case, HPM classifies a set of overlapping protein fragments, extracted from a non-redundant databank of TMP 3D structure. After fine-tuning of the learning parameters, the optimal classification results in 65 clusters. They represent at best similar relationships between sequence and local structure properties of TMPs. Interestingly, HPM distinguishes among the resulting clusters two helical regions with distinct hydrophobic patterns. This underlines the complexity of the topology of these proteins. The HPM classification enlightens unusual relationship between amino acids in TMP fragments, which can be useful to elaborate new amino acids substitution matrices. Finally, two challenging applications are described: the first one aims at annotating protein functions (channel or not), the second one intends to assess the quality of the structures (X-ray or models) via a new scoring function deduced from the HPM classification.

VLDP web server: a powerful geometric tool for analysing protein structures in their environment.

Protein structures are an ensemble of atoms determined experimentally mostly by X-ray crystallography or Nuclear Magnetic Resonance. Studying 3D protein structures is a key point for better understanding protein function at a molecular level. We propose a set of accurate tools, for analysing protein structures, based on the reliable method of Voronoi-Laguerre tessellations. The Voronoi Laguerre Delaunay Protein web server (VLDPws) computes the Laguerre tessellation on a whole given system first embedded in solvent. Through this fine description, VLDPws gives the following data: (i) Amino acid volumes evaluated with high precision, as confirmed by good correlations with experimental data. (ii) A novel definition of inter-residue contacts within the given protein. (iii) A measure of the residue exposure to solvent that significantly improves the standard notion of accessibility in some cases. At present, no equivalent web server is available. VLDPws provides output in two complementary forms: direct visualization of the Laguerre tessellation, mostly its polygonal molecular surfaces; files of volumes; and areas, contacts and similar data for each residue and each atom. These files are available for download for further analysis. VLDPws can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/vldp.

A functional mini-GDE transgene corrects impairment in models of glycogen storage disease type III.

Glycogen storage disease type III (GSDIII) is a rare inborn error of metabolism affecting liver, skeletal muscle, and heart due to mutations of the AGL gene encoding for the glycogen debranching enzyme (GDE). No curative treatment exists for GSDIII. The 4.6 kb GDE cDNA represents the major technical challenge toward the development of a single recombinant adeno-associated virus-derived (rAAV-derived) vector gene therapy strategy. Using information on GDE structure and molecular modeling, we generated multiple truncated GDEs. Among them, an N-terminal-truncated mutant, ΔNter2-GDE, had a similar efficacy in vivo compared with the full-size enzyme. A rAAV vector expressing ΔNter2-GDE allowed significant glycogen reduction in heart and muscle of Agl-/- mice 3 months after i.v. injection, as well as normalization of histology features and restoration of muscle strength. Similarly, glycogen accumulation and histological features were corrected in a recently generated Agl-/- rat model. Finally, transduction with rAAV vectors encoding ΔNter2-GDE corrected glycogen accumulation in an in vitro human skeletal muscle cellular model of GSDIII. In conclusion, our results demonstrated the ability of a single rAAV vector expressing a functional mini-GDE transgene to correct the muscle and heart phenotype in multiple models of GSDIII, supporting its clinical translation to patients with GSDIII.

Modeling Protein Complexes and Molecular Assemblies Using Computational Methods.

Many biological molecules are assembled into supramolecular complexes that are necessary to perform functions in the cell. Better understanding and characterization of these molecular assemblies are thus essential to further elucidate molecular mechanisms and key protein-protein interactions that could be targeted to modulate the protein binding affinity or develop new binders. Experimental access to structural information on these supramolecular assemblies is often hampered by the size of these systems that make their recombinant production and characterization rather difficult. Computational methods combining both structural data, molecular modeling techniques, and sequence coevolution information can thus offer a good alternative to gain access to the structural organization of protein complexes and assemblies. Herein, we present some computational methods to predict structural models of the protein partners, to search for interacting regions using coevolution information, and to build molecular assemblies. The approach is exemplified using a case study to model the succinate-quinone oxidoreductase heterocomplex.

Mutagenesis of Dimer Interfacial Residues Improves the Activity and Specificity of Methyltransferase for cis-α-Irone Biosynthesis.

Promiscuous enzymes show great potential to establish new-to-nature pathways and expand chemical diversity. Enzyme engineering strategies are often employed to tailor such enzymes to improve their activity or specificity. It is paramount to identify the target residues to be mutated. Here, by exploring the inactivation mechanism with the aid of mass spectrometry, we have identified and mutated critical residues at the dimer interface region of the promiscuous methyltransferase (pMT) that converts psi-ionone to irone. The optimized pMT12 mutant showed ∼1.6-4.8-fold higher kcat than the previously reported best mutant, pMT10, and increased the cis-α-irone percentage from ∼70 to ∼83%. By one-step biotransformation, ∼121.8 mg L-1 cis-α-irone was produced from psi-ionone by the pMT12 mutant. The study offers new opportunities to engineer enzymes with enhanced activity and specificity.

Total enzymatic synthesis of cis-α-irone from a simple carbon source.

Metabolic engineering has become an attractive method for the efficient production of natural products. However, one important pre-requisite is to establish the biosynthetic pathways. Many commercially interesting molecules cannot be biosynthesized as their native biochemical pathways are not fully elucidated. Cis-α-irone, a top-end perfumery molecule, is an example. Retrobiosynthetic pathway design by employing promiscuous enzymes provides an alternative solution to this challenge. In this work, we design a synthetic pathway to produce cis-α-irone with a promiscuous methyltransferase (pMT). Using structure-guided enzyme engineering strategies, we improve pMT activity and specificity towards cis-α-irone by >10,000-fold and >1000-fold, respectively. By incorporating the optimized methyltransferase into our engineered microbial cells, ~86 mg l-1 cis-α-irone is produced from glucose in a 5 l bioreactor. Our work illustrates that integrated retrobiosynthetic pathway design and enzyme engineering can offer opportunities to expand the scope of natural molecules that can be biosynthesized.

An integrated approach reveals how lipo-chitooligosaccharides interact with the lysin motif receptor-like kinase MtLYR3.

N-acetylglucosamine containing compounds acting as pathogenic or symbiotic signals are perceived by plant-specific Lysin Motif Receptor-Like Kinases (LysM-RLKs). The molecular mechanisms of this perception are not fully understood, notably those of lipo-chitooligosaccharides (LCOs) produced during root endosymbioses with nitrogen-fixing bacteria or arbuscular mycorrhizal fungi. In Medicago truncatula, we previously identified the LysM-RLK LYR3 (MtLYR3) as a specific LCO-binding protein. We also showed that the absence of LCO binding to LYR3 of the non-mycorrhizal Lupinus angustifolius, (LanLYR3), was related to LysM3, which differs from that of MtLYR3 by several amino acids and, particularly, by a critical tyrosine residue absent in LanLYR3. Here, we aimed to define the LCO binding site of MtLYR3 by using molecular modelling and simulation approaches, combined with site-directed mutagenesis and LCO binding experiments. 3D models of MtLYR3 and LanLYR3 ectodomains were built, and homology modelling and molecular dynamics (MD) simulations were performed. Molecular docking and MD simulation on the LysM3 identified potential key residues for LCO binding. We highlighted by steered MD simulations that in addition to the critical tyrosine, two other residues were important for LCO binding in MtLYR3. Substitution of these residues in LanLYR3-LysM3 by those of MtLYR3-LysM3 allowed the recovery of high-affinity LCO binding in experimental radioligand-binding assays. An analysis of selective constraints revealed that the critical tyrosine has experienced positive selection pressure and is absent in some LYR3 proteins. These findings now pave the way to uncover the functional significance of this specific evolutionary pattern.

Comparative analysis of threshold and tessellation methods for determining protein contacts.

The 3D structure of a protein is the main physical support of a protein's biological function; 3D protein folds are primarily maintained through interactions between amino acids. Inter-residue contacts are essential for the stability of protein folds. Therefore, many methodologies in the fields of structure analysis, structure prediction, and structure-function relationships are based on residue contacts. The present study provides a comparative analysis of two approaches for determining contacts: the classical distance-threshold method and an application of Laguerre, or weighted Voronoi tessellation. First, we examined mean contact distributions and their dependence on residue volumes, accessibility and hydrophobicity. In general, the different methods gave concordant results, although the method based on Cα distances showed significant discrepancies with the all-atom tessellation method. We also analyzed preferential contacts between all amino acid species and studied the influence of protein chain length, the proximity of the residues along the sequence, and the secondary structure environment. Interestingly, the discrepancies between methods were occasionally large enough to substantially change the relative preferences of some contacts. Finally, a case study on disulfide bridges demonstrated the importance of the structural environment in determining contacts from tessellation. In conclusion, the tessellation method is more accurate because of its fine adaptation to local protein topology, with far-reaching implications for most contact-based prediction methods of protein folding.

Surface charge distribution: a key parameter for understanding protein behavior in chromatographic processes.

Multi-component adsorption of proteins still requires a better understanding of local phenomena to improve the development of predictive models. In this work, all-atom Molecular Dynamics (MD) simulations were used to investigate the influence of protein charge distribution on the adsorption capacity. The simultaneous adsorption of α-chymotrypsin and lysozyme on a cation exchanger, SP Sepharose FF, was studied through MD simulations and compared to macroscopic isotherm experiments. It appears that the charge distribution is a relevant information to better understand specific phenomena, such as a multilayer adsorption caused by the particular electrostatic profile of α-chymotrypsin. Therefore, MD simulations seem to be an interesting way to visualize and highlight these behaviors.

A tripartite carbohydrate-binding module to functionalize cellulose nanocrystals.

The development of protein and microorganism engineering have led to rising expectations of biotechnology in the design of emerging biomaterials, putatively of high interest to reduce our dependence on fossil carbon resources. In this way, cellulose, a renewable carbon based polysaccharide and derived products, displays unique properties used in many industrial applications. Although the functionalization of cellulose is common, it is however limited in terms of number and type of functions. In this work, a Carbohydrate-Binding Module (CBM) was used as a central core to provide a versatile strategy to bring a large diversity of functions to cellulose surfaces. CBM3a from Clostridium thermocellum, which has a high affinity for crystalline cellulose, was flanked through linkers with a streptavidin domain and an azide group introduced through a non-canonical amino acid. Each of these two extra domains was effectively produced and functionalized with a variety of biological and chemical molecules. Structural properties of the resulting tripartite chimeric protein were investigated using molecular modelling approaches, and its potential for the multi-functionalization of cellulose was confirmed experimentally. As a proof of concept, we show that cellulose can be labelled with a fluorescent version of the tripartite protein grafted to magnetic beads and captured using a magnet.

A novel evaluation of residue and protein volumes by means of Laguerre tessellation.

Amino acids control the protein folding process and maintain its functional fold. This study underlines the interest of the Laguerre tessellation to determine relevant amino acid volumes in proteins. Previous studies used a limited number of proteins and only buried residues. The present computations improve the method and results on three main points: (i) a large, high-quality updated and refined data bank of proteins is used; (ii) all residues are taken into account, including those at the protein surface, thanks to (iii) the addition of a realistic solvent. The new values of the average and standard deviation of amino acid volumes show significant corrections with respect to previous studies. Another issue of the method is the polyhedral protein/water interface area (PIA) which quantifies the exposure of atoms or residues to the solvent. We propose this PIA as a new, parameter-free, alternative for measuring accessibility. The comparison with NACCESS is satisfactory; however, the methods disagree in pointing out buried residues: where NACCESS evaluates to zero, the exposure given by PIA ranges from 0 to 20%. Variations of average residue volumes have been analyzed under several conditions, e.g., how they depend on protein size and on secondary structure environments. As it is based on strong mathematical grounds and on numerous high-quality protein structures, our work gives a reliable methodology and up-to-date values of amino acid volumes and surface accessibility.

Understanding adsorption behavior of α-chymotrypsin onto cation exchanger using all-atom molecular dynamics simulations.

The interest for a better understanding of ion-exchange mechanisms at the atomic level has strongly increased over the past decades. Indeed, molecular-level information about physico-chemical mechanisms could help optimizing chromatographic processes for protein purification, which are sub-optimized due to the lack of predictive models. A promising approach is based on the use of Molecular Dynamics (MD) simulations to study local phenomena inside adsorbents which can then be challenged against experimental results. In this work, macroscopic experimental data, consisting in the ion-exchange uptake of α-chymotrypsin onto SP Sepharose FF, have been compared to the adsorption behavior predicted by MD simulations. The chromatographic surface, represented as a uniform distribution of ligands with a counterion layer, in the presence of the protein was modeled using all-atom representation. The SMA formalism was used to describe single adsorption isotherms and to relate macroscopic observations with molecular simulations. Two SMA parameters based on physical principles, the characteristic charge n and the steric factor σ, have been estimated by both experiments and MD simulations. At pH 5 and NaCl concentration of 100 mM, our study shows a fairly good agreement between both results, especially for the characteristic charge. It is shown that the steric factor calculation is strongly dependent on the ligand density on the adsorbent surface, whose value must be carefully determined in order to obtain reliable predictions. In addition, four binding patches were identified as being involved in the adsorption, which have been confirmed through binding free energy calculations.

Creation of (R)-Amine Transaminase Activity within an α-Amino Acid Transaminase Scaffold.

The enzymatic transamination of ketones into (R)-amines represents an important route for accessing a range of pharmaceuticals or building blocks. Although many publications have dealt with enzyme discovery, protein engineering, and the application of (R)-selective amine transaminases [(R)-ATA] in biocatalysis, little is known about the actual in vivo role and how these enzymes have evolved from the ubiquitous α-amino acid transaminases (α-AATs). Here, we show the successful introduction of an (R)-transaminase activity in an α-amino acid aminotransferase with one to six amino acid substitutions in the enzyme's active site. Bioinformatic analysis combined with computational redesign of the d-amino acid aminotransferase (DATA) led to the identification of a sextuple variant having a specific activity of 326 milliunits mg-1 in the conversion of (R)-phenylethylamine and pyruvate to acetophenone and d-alanine. This value is similar to those of natural (R)-ATAs, which typically are in the range of 250 milliunits mg-1. These results demonstrate that (R)-ATAs can evolve from α-AAT as shown here for the DATA scaffold.

Analysis of the diversity of the glycoside hydrolase family 130 in mammal gut microbiomes reveals a novel mannoside-phosphorylase function.

Mannoside phosphorylases are involved in the intracellular metabolization of mannooligosaccharides, and are also useful enzymes for the in vitro synthesis of oligosaccharides. They are found in glycoside hydrolase family GH130. Here we report on an analysis of 6308 GH130 sequences, including 4714 from the human, bovine, porcine and murine microbiomes. Using sequence similarity networks, we divided the diversity of sequences into 15 mostly isofunctional meta-nodes; of these, 9 contained no experimentally characterized member. By examining the multiple sequence alignments in each meta-node, we predicted the determinants of the phosphorolytic mechanism and linkage specificity. We thus hypothesized that eight uncharacterized meta-nodes would be phosphorylases. These sequences are characterized by the absence of signal peptides and of the catalytic base. Those sequences with the conserved E/K, E/R and Y/R pairs of residues involved in substrate binding would target ß-1,2-, ß-1,3- and ß-1,4-linked mannosyl residues, respectively. These predictions were tested by characterizing members of three of the uncharacterized meta-nodes from gut bacteria. We discovered the first known ß-1,4-mannosyl-glucuronic acid phosphorylase, which targets a motif of the Shigella lipopolysaccharide O-antigen. This work uncovers a reliable strategy for the discovery of novel mannoside-phosphorylases, reveals possible interactions between gut bacteria, and identifies a biotechnological tool for the synthesis of antigenic oligosaccharides.

α-Galactosidase and Sucrose-Kinase Relationships in a Bi-functional AgaSK Enzyme Produced by the Human Gut Symbiont Ruminococcus gnavus E1.

Plant α-galactosides belonging to the raffinose family oligosaccharides (RFOs) and considered as prebiotics, are commonly degraded by α-galactosidases produced by the human gut microbiome. In this environment, the Ruminococcus gnavus E1 symbiont-well-known for various benefit-is able to produce an original RgAgaSK bifunctional enzyme. This enzyme contains an hydrolytic α-galactosidase domain linked to an ATP dependent extra-domain, specifically involved in the α-galactoside hydrolysis and the phosphorylation of the glucose, respectively. However, the multi-modular relationships between both catalytic domains remained hitherto unexplored and has been, consequently, herein investigated. Biochemical characterization of heterologously expressed enzymes either in full-form or in separated domains revealed similar kinetic parameters. These results were supported by molecular modeling studies performed on the whole enzyme in complex with different RFOs. Further enzymatic analysis associated with kinetic degradation of various substrates followed by high pressure anionic exchange chromatography revealed that catalytic efficiency decreased as the number of D-galactosyl moieties branched onto the oligosaccharide increased, suggesting a preference of RgAgaSK for RFO's short chains. A wide prevalence and abundance study on a human metagenomic library showed a high prevalence of the RgAgaSK encoding gene whatever the health status of the individuals. Finally, phylogeny and synteny studies suggested a limited spread by horizontal transfer of the clusters' containing RgAgaSK to only few species of Firmicutes, highlighting the importance of these undispersed tandem activities in the human gut microbiome.