ABSTRACT
Enzyme engineering approaches have allowed to extend the collection of enzymatic tools available for synthetic purposes. However, controlling the regioselectivity of the reaction remains challenging, in particular when dealing with carbohydrates bearing numerous reactive hydroxyl groups as substrates. Here, we used a computer-aided design framework to engineer the active site of a sucrose-active [Formula: see text]-transglucosylase for the 1,2-cis-glucosylation of a lightly protected chemically synthesized tetrasaccharide, a common precursor for the synthesis of serotype-specific S. flexneri O-antigen fragments. By targeting 27 amino acid positions of the acceptor binding subsites of a GH70 branching sucrase, we used a RosettaDesign-based approach to propose 49 mutants containing up to 15 mutations scattered over the active site. Upon experimental evaluation, these mutants were found to produce up to six distinct pentasaccharides, whereas only two were synthesized by the parental enzyme. Interestingly, we showed that by introducing specific mutations in the active site of a same enzyme scaffold, it is possible to control the regiospecificity of the 1,2-cis glucosylation of the tetrasaccharide acceptor and produce a unique diversity of pentasaccharide bricks. This work offers novel opportunities for the development of highly convergent chemo-enzymatic routes toward S. flexneri haptens.
Subject(s)
Glucose/analysis , Glucose/chemistry , Oligosaccharides/chemistry , Polysaccharides/chemistry , Sucrase/chemistry , Biotechnology , Carbohydrates/chemistry , Catalytic Domain , Chromatography, High Pressure Liquid , Computer-Aided Design , Enzymes/chemistry , Glycosylation , Haptens , Hydrolases/metabolism , Molecular Biology , Mutation , O Antigens , Protein Engineering/methods , Shigella flexneri , Spectrometry, Mass, Electrospray Ionization , Tandem Mass SpectrometryABSTRACT
The (chemo-)enzymatic synthesis of oligosaccharides has been hampered by the lack of appropriate enzymatic tools with requisite regio- and stereo-specificities. Engineering of carbohydrate-active enzymes, in particular targeting the enzyme active site, has notably led to catalysts with altered regioselectivity of the glycosylation reaction thereby enabling to extend the repertoire of enzymes for carbohydrate synthesis. Using a collection of 22 mutants of ΔN123-GBD-CD2 branching sucrase, an enzyme from the Glycoside Hydrolase family 70, containing between one and three mutations in the active site, and a lightly protected chemically synthesized tetrasaccharide as an acceptor substrate, we showed that altered glycosylation product specificities could be achieved compared to the parental enzyme. Six mutants were selected for further characterization as they produce higher amounts of two favored pentasaccharides compared to the parental enzyme and/or new products. The produced pentasaccharides were shown to be of high interest as they are precursors of representative haptens of Shigella flexneri serotypes 3a, 4a and 4b. Furthermore, their synthesis was shown to be controlled by the mutations introduced in the active site, driving the glucosylation toward one extremity or the other of the tetrasaccharide acceptor. To identify the molecular determinants involved in the change of ΔN123-GBD-CD2 regioselectivity, extensive molecular dynamics simulations were carried out in combination with in-depth analyses of amino acid residue networks. Our findings help to understand the inter-relationships between the enzyme structure, conformational flexibility and activity. They also provide new insight to further engineer this class of enzymes for the synthesis of carbohydrate components of bacterial haptens.
Subject(s)
Bacterial Proteins , Haptens/biosynthesis , Oligosaccharides/biosynthesis , Protein Engineering , Shigella flexneri/metabolism , Sucrase , Bacterial Proteins/genetics , Bacterial Proteins/metabolism , Haptens/genetics , Oligosaccharides/genetics , Shigella flexneri/genetics , Sucrase/genetics , Sucrase/metabolismABSTRACT
Progress in glycoscience is strongly dependent on the availability of broadly diverse tailor-made, well-defined, and often complex oligosaccharides. Herein, going beyond natural resources and aiming to circumvent chemical boundaries in glycochemistry, we tackle the development of an in vitro chemoenzymatic strategy holding great potential to answer the need for molecular diversity characterizing microbial cell-surface carbohydrates. The concept is exemplified in the context of Shigella flexneri, a major cause of diarrhoeal disease. Aiming at a broad serotype coverage S. flexneri glycoconjugate vaccine, a non-natural lightly protected tetrasaccharide was designed for compatibility with (i) serotype-specific glucosylations and O-acetylations defining S. flexneri O-antigens, (ii) recognition by suitable α-transglucosylases, and (iii) programmed oligomerization following enzymatic α-d-glucosylation. The tetrasaccharide core was chemically synthesized from two crystalline monosaccharide precursors. Six α-transglucosylases found in the glycoside hydrolase family 70 were shown to transfer glucosyl residues on the non-natural acceptor. The successful proof of concept is achieved for a pentasaccharide featuring the glucosylation pattern from the S. flexneri type IV O-antigen. It demonstrates the potential of appropriately planned chemoenzymatic pathways involving non-natural acceptors and low-cost donor/transglucosylase systems to achieve the demanding regioselective α-d-glucosylation of large substrates, paving the way to microbial oligosaccharides of vaccinal interest.
Subject(s)
O Antigens , Shigella flexneri , Carbohydrate Sequence , Oligosaccharides , SerogroupABSTRACT
Among α-transglucosylases from Glycoside-Hydrolase family 70, the ΔN123-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 ΔN123-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.
Subject(s)
Leuconostoc/enzymology , Sucrase/chemistry , Catalytic Domain , Crystallography, X-Ray , Leuconostoc/chemistry , Leuconostoc/metabolism , Molecular Dynamics Simulation , Protein Binding , Protein Conformation , Sucrase/metabolism , Sucrose/metabolismABSTRACT
Recently, a 3-hydroxychromone based nucleoside 3HCnt has been developed as a highly environment-sensitive nucleoside surrogate to investigate protein-DNA interactions. When it is incorporated in DNA, the probe is up to 50-fold brighter than 2-aminopurine, the reference fluorescent nucleoside. Although the insertion of 3HCnt in DNA was previously shown to not alter the overall DNA structure, the possibility of the probe inducing local effects cannot be ruled out. Hence, a systematic structural and dynamic study is required to unveil the 3HCnt's limitations and to properly interpret the data obtained with this universal probe. Here, we investigated by NMR a 12-mer duplex, in which a central adenine was replaced by 3HCnt. The chemical shifts variations and nOe contacts revealed that the 3HCnt is well inserted in the DNA double helix with extensive stacking interactions with the neighbor base pairs. These observations are in excellent agreement with the steady-state and time-resolved fluorescence properties indicating that the 3HCnt fluorophore is protected from the solvent and does not exhibit rotational motion. The 3HCnt insertion in DNA is accompanied by the extrusion of the opposite nucleobase from the double helix. Molecular dynamics simulations using NMR-restraints demonstrated that 3HCnt fluorophore exhibits only translational dynamics. Taken together, our data showed an excellent intercalation of 3HCnt in the DNA double helix, which is accompanied by localized perturbations. This confirms 3HCnt as a highly promising tool for nucleic acid labeling and sensing.
Subject(s)
Chromones/chemistry , DNA/chemistry , Fluorescence , Molecular Dynamics Simulation , Nucleic Acid ConformationABSTRACT
The accurate prediction of the structure and dynamics of DNA remains a major challenge in computational biology due to the dearth of precise experimental information on DNA free in solution and limitations in the DNA force-fields underpinning the simulations. A new generation of force-fields has been developed to better represent the sequence-dependent B-DNA intrinsic mechanics, in particular with respect to the BI â BII backbone equilibrium, which is essential to understand the B-DNA properties. Here, the performance of MD simulations with the newly updated force-fields Parmbsc0εζOLI and CHARMM36 was tested against a large ensemble of recent NMR data collected on four DNA dodecamers involved in nucleosome positioning. We find impressive progress towards a coherent, realistic representation of B-DNA in solution, despite residual shortcomings. This improved representation allows new and deeper interpretation of the experimental observables, including regarding the behavior of facing phosphate groups in complementary dinucleotides, and their modulation by the sequence. It also provides the opportunity to extensively revisit and refine the coupling between backbone states and inter base pair parameters, which emerges as a common theme across all the complementary dinucleotides. In sum, the global agreement between simulations and experiment reveals new aspects of intrinsic DNA mechanics, a key component of DNA-protein recognition.
Subject(s)
Models, Chemical , Molecular Dynamics Simulation , Nucleic Acid Conformation , Sequence Analysis, Protein/methods , Base Sequence , Molecular Sequence Data , Structure-Activity RelationshipABSTRACT
We investigated how the intrinsic sequence-dependent properties probed via the phosphate linkages (BI â BII equilibrium) influence the preferred shape of free DNA, and how this affects the nucleosome formation. First, this exploits NMR solution studies of four B-DNA dodecamers that together cover 39 base pairs of the 5' half of the sequence 601, of special interest for nucleosome formation. The results validate our previous prediction of a systematic, general sequence effect on the intrinsic backbone BII propensities. NMR provides new evidence that the backbone behavior is intimately coupled to the minor groove width. Second, application of the backbone behavior predictions to the full sequence 601 and other relevant sequences demonstrates that alternation of intrinsic low and high BII propensities, coupled to intrinsic narrow and wide minor grooves, largely coincides with the sinusoidal variations of the DNA minor groove width observed in crystallographic structures of the nucleosome. This correspondence is much poorer with low affinity sequences. Overall, the results indicate that nucleosome formation involves an indirect readout process implicating pre-existing DNA minor groove conformations. It also illustrates how the prediction of the intrinsic structural DNA behavior offers a powerful framework to gain explanatory insight on how proteins read DNA.