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.

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.

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 k_cat 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.

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.

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.

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.

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.

α-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.

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.

Analyzing protein topology based on Laguerre tessellation of a pore-traversing water network

Given the tight relation between protein structure and function, we present a set of methods to analyze protein topology, implemented in the VLDP program, relying on Laguerre space partitions built from series of molecular dynamics snapshots. The Laguerre partition specifies inter-atomic contacts, formalized in graphs. The deduced properties are the existence and count of water aggregates, possible passage ways and constrictions, the structure, connectivity, stability and depth of the water network. As a test-case, the membrane protein FepA is investigated in its full environment, yielding a more precise description of the protein surface. Inside FepA, the solvent splits into isolated clusters and an intricate network connecting both sides of the lipid bilayer. The network is dynamic, connections set on and off, occasionally substantially relocating traversing paths. Subtle differences are detected between two forms of FepA, ligand-free and complexed with its natural iron carrier, the enterobactin. The complexed form has more constricted and more centered openings in the upper part whereas, in the lower part, constriction is released: two main channels between the plug and barrel lead directly to the periplasm. Reliability, precision and the variety of topological features are the main interest of the method.

Combining multi-scale modelling methods to decipher molecular motions of a branching sucrase from glycoside-hydrolase family 70

Among α-transglucosylases from Glycoside-Hydrolase family 70, the ΔN₁₂₃-GB-CD2 enzyme derived from the bifunctional DSR-E from L. citreum NRRL B-1299 is particularly interesting as it was the first described engineered Branching Sucrase, not able to elongate glucan polymers from sucrose substrate. The previously reported overall structural organization of this multi-domain enzyme is an intricate U-shape fold conserved among GH70 enzymes which showed a certain conformational variability of the so-called domain V, assumed to play a role in the control of product structures, in available X-ray structures. Understanding the role of functional dynamics on enzyme reaction and substrate recognition is of utmost interest although it remains a challenge for biophysical methods. By combining long molecular dynamics simulation (1μs) and multiple analyses (NMA, PCA, Morelet Continuous Wavelet Transform and Cross Correlations Dynamics), we investigated here the dynamics of ΔN₁₂₃-GB-CD2 alone and in interaction with sucrose substrate. Overall, our results provide the detailed picture at atomic level of the hierarchy of motions occurring along different timescales and how they are correlated, in agreement with experimental structural data. In particular, detailed analysis of the different structural domains revealed cooperative dynamic behaviors such as twisting, bending and wobbling through anti- and correlated motions, and also two structural hinge regions, of which one was unreported. Several highly flexible loops surrounding the catalytic pocket were also highlighted, suggesting a potential role in the acceptor promiscuity of ΔN123-GBD-CD2. Normal modes and essential dynamics underlined an interesting two-fold dynamic of the catalytic domain A, pivoting about an axis splitting the catalytic gorge in two parts. The comparison of the conformational free energy landscapes using principal component analysis of the enzyme in absence or in presence of sucrose, also revealed a more harmonic basin when sucrose is bound with a shift population of the bending mode, consistent with the substrate binding event.

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.