Biocatalyzed synthesis of difuranosides and their ability to trigger production of TNF-α.

Transglycosylation reactions biocatalyzed by the native arabinofuranosidase Araf51 and using d-galactosyl, d-fucosyl and 6-deoxy-6-fluoro-D-galactosyl derivatives as donors and acceptors provided di-to pentahexofuranosides. The immunostimulatory potency of these compounds, and more especially their ability to induce production of TNF-α, was evaluated on the murine macrophage cell line, Raw 264.7. The results obtained showed concentration-dependent and most importantly, structure-dependent responses. Interestingly, oligoarabinofuranosides belonging to the oligopentafuranoside family displayed concentration-, chain length and aglycon-dependent bioactivities irrespective of their fine chemical variations. Thus, neo-oligofuranosides in D-Galf series, as well as their D-Fucf and 6-fluorinated counterparts are indeed potential sources of immunostimulating agents.

Alpha-L-fucosidase isoenzyme iso2 from Paenibacillus thiaminolyticus.

BACKGROUND: α-L-Fucosidases are enzymes involved in metabolism of α-L-fucosylated molecules, compounds with a fundamental role in different life essential processes including immune response, fertilization and development, but also in some serious pathological events. According to the CAZy database, these enzymes belong to families 29 and 95. Some of them are also reported to be able to catalyze transglycosylation reactions, during which α-L-fucosylated molecules, representing compounds of interest especially for pharmaceutical industry, are formed. METHODS: Activity-based screening of a genomic library was used to isolate the gene encoding a novel α-L-fucosidase. The enzyme was expressed in E.coli and affinity chromatography was used for purification of His-tagged α-L-fucosidase. Standard activity assay was used for enzyme characterization. Thin layer chromatography and mass spectrometry were used for transglycosylation reactions evaluation. RESULTS: Using a genomic library of Paenibacillus thiaminolyticus, constructed in E.coli DH5α cells, nucleotide sequence of a new α-L-fucosidase isoenzyme was determined and submitted to the EMBL database (HE654122). However, no similarity with enzymes from CAZy database families 29 and 95 was detected. This enzyme was produced in form of histidine-tagged protein in E.coli BL21 (DE3) cells and purified by metaloaffinity chromatography. Hydrolytic and transglycosylation abilities of α-L-fucosidase iso2 were tested using different acceptor molecules. CONCLUSIONS: In this study, new enzyme α-L-fucosidase iso2 originating from Paenibacillus thiaminolyticus was described and prepared in recombinant form and its hydrolytic and transglycosylation properties were characterized. As a very low amino acid sequence similarity with known α-L-fucosidases was found, following study could be important for different biochemical disciplines involving molecular modelling.

Nonlinear vs. linear biasing in Trp-cage folding simulations.

Biased simulations have great potential for the study of slow processes, including protein folding. Atomic motions in molecules are nonlinear, which suggests that simulations with enhanced sampling of collective motions traced by nonlinear dimensionality reduction methods may perform better than linear ones. In this study, we compare an unbiased folding simulation of the Trp-cage miniprotein with metadynamics simulations using both linear (principle component analysis) and nonlinear (Isomap) low dimensional embeddings as collective variables. Folding of the mini-protein was successfully simulated in 200 ns simulation with linear biasing and non-linear motion biasing. The folded state was correctly predicted as the free energy minimum in both simulations. We found that the advantage of linear motion biasing is that it can sample a larger conformational space, whereas the advantage of nonlinear motion biasing lies in slightly better resolution of the resulting free energy surface. In terms of sampling efficiency, both methods are comparable.

The versatile enzyme Araf51 allowed efficient synthesis of rare pathogen-related ß-D-galactofuranosyl-pyranoside disaccharides.

The preparation of galactofuranosyl-containing disaccharidic parts of natural glycoconjugates was performed according to a chemo-enzymatic synthesis. Our goals were firstly to develop an alternative approach to standard chemical strategies by limiting the number of reaction and purification steps, and secondly to evaluate the scope of the Araf51 biocatalyst to transfer a galactofuranosyl moiety to a set of pyranosidic acceptors differing from each other by the series, the anomeric configuration as well as the conformation. The study of binding mode of the resulting disaccharides was also performed by molecular modeling and showed significant differences between (1→2)- and (1→6)-linked disaccharides.

α-L-fucosidase from Paenibacillus thiaminolyticus: its hydrolytic and transglycosylation abilities.

In this work, focused on possible application of α-L-fucosidases from bacterial sources in the synthesis of α-L-fucosylated glycoconjugates, several nonpathogenic aerobic bacterial strains were screened for α-L-fucosidase activity. Among them Paenibacillus thiaminolyticus was confirmed as a potent producer of enzyme with the ability to cleave the chromogenic substrate p-nitrophenyl α-L-fucopyranoside. The gene encoding α-L-fucosidase was found using the genomic library of P. thiaminolyticus constructed in the cells of Escherichia coli DH5α and sequenced (EMBL database: FN869117, carbohydrate-active enzymes database: Glycosidase family 29). The enzyme was expressed in the form of polyhistidine-tagged protein (51.2 kDa) in Escherichia coli BL21 (DE3) cells, purified using nickel-nitrilotriacetic acid agarose affinity chromatography and characterized using the chromogenic substrate p-nitrophenyl α-L-fucopyranoside (K(m) = (0.44 ± 0.02) mmol/L, K(S) = (83 ± 8) mmol/L (substrate inhibition), pH(optimum) = 8.2, t(optimum) = 48°C). By testing the ability of the enzyme to catalyze the transfer of α-L-fucosyl moiety to different types of acceptor molecules, it was confirmed that the enzyme is able to catalyze the formation of α-L-fucosylated p-nitrophenyl glycopyranosides containing α-D-galactopyranosidic, α-D-glucopyranosidic, α-D-mannopyranosidic or α-L-fucopyranosidic moiety. This enzyme is also able to catalyze α-L-fucosylation of aliphatic alcohols of different lenghs of alkyl chain and hydroxyl group positions (methanol, ethanol, 1-propanol, 2-propanol and 1-octanol) and hydroxyl group-containing amino acid derivatives (N-(tert-butoxycarbonyl)-L-serine methyl ester and N-(tert-butoxycarbonyl)-L-threonine methyl ester). These results indicate the possibility of exploiting this enzyme in the synthesis of different types of α-L-fucosylated molecules representing compounds with potential application in biotechnology and the pharmaceutical industry.

Conformational free energy modeling of druglike molecules by metadynamics in the WHIM space.

Protein-ligand affinities can be significantly influenced not only by the interaction itself but also by conformational equilibrium of both binding partners, free ligand and free protein. Identification of important conformational families of a ligand and prediction of their thermodynamics is important for efficient ligand design. Here we report conformational free energy modeling of nine small-molecule drugs in explicitly modeled water by metadynamics with a bias potential applied in the space of weighted holistic invariant molecular (WHIM) descriptors. Application of metadynamics enhances conformational sampling compared to unbiased molecular dynamics simulation and allows to predict relative free energies of key conformations. Selected free energy minima and one example of transition state were tested by a series of unbiased molecular dynamics simulation. Comparison of free energy surfaces of free and target-bound Imatinib provides an estimate of free energy penalty of conformational change induced by its binding to the target.

An enzymatic glycosylation of nucleoside analogues using ß-galactosidase from Escherichia coli.

A new enzymatic method for the synthesis of ß-galactosides of nucleosides and acyclic nucleoside analogues has been developed, using ß-galactosidase from Escherichia coli as a catalyst and lactose as a sugar donor. The method is very rapid, feasible and last but not least inexpensive. Its applicability has been proven for a broad variety of possible substrates with respect to its scaling up for preparative use. Five new compounds from a series of nucleoside and acyclic nucleoside analogues have been prepared on a scale of several hundred milligrams, in all cases revealing very good results of the method concerning the reproducibility of the reaction yields and simplicity of the purification process.

Metadynamics in the conformational space nonlinearly dimensionally reduced by Isomap.

Atomic motions in molecules are not linear. This infers that nonlinear dimensionality reduction methods can outperform linear ones in analysis of collective atomic motions. In addition, nonlinear collective motions can be used as potentially efficient guides for biased simulation techniques. Here we present a simulation with a bias potential acting in the directions of collective motions determined by a nonlinear dimensionality reduction method. Ad hoc generated conformations of trans,trans-1,2,4-trifluorocyclooctane were analyzed by Isomap method to map these 72-dimensional coordinates to three dimensions, as described by Brown and co-workers [J. Chem. Phys. 129, 064118 (2008)]. Metadynamics employing the three-dimensional embeddings as collective variables was applied to explore all relevant conformations of the studied system and to calculate its conformational free energy surface. The method sampled all relevant conformations (boat, boat-chair, and crown) and corresponding transition structures inaccessible by an unbiased simulation. This scheme allows to use essentially any parameter of the system as a collective variable in biased simulations. Moreover, the scheme we used for mapping out-of-sample conformations from the 72D to 3D space can be used as a general purpose mapping for dimensionality reduction, beyond the context of molecular modeling.

beta-D-Galactosidase from Paenibacillus thiaminolyticus catalyzing transfucosylation reactions.

A genomic library of bacterial strain Paenibacillus thiaminolyticus was constructed and the plasmid DNA of the clone, containing the gene encoding beta-d-galactosidase with beta-d-fucosidase activity, detected by 5-bromo-4-chloro-3-indoxyl beta-d-galactopyranoside, was sequenced. Cells of Escherichia coli BL21 (DE3) were used for production of the enzyme in the form of a histidine-tagged protein. This recombinant fusion protein was purified using Ni-NTA agarose affinity chromatography and characterized by using p-nitrophenyl beta-d-fucopyranoside (K(m) value of (1.18 +/- 0.06) mmol/L), p-nitrophenyl beta-d-galactopyranoside (K(m) value of (250 +/- 40) mmol/L), p-nitrophenyl beta-d-glucopyranoside (K(m) value of (77 +/- 6) mmol/L), and lactose (K(m) value of (206 +/- 5) mmol/L) as substrates. Optimal pH and temperature were estimated as 5.5 and 65 degrees C, respectively. According to the amino acid sequence, the molecular weight of the fusion protein was calculated to be 68.6 kDa and gel filtration chromatography confirmed the presence of the enzyme in a monomeric form. In the following step, its ability to catalyze transfucosylation reactions was tested. The enzyme was able to catalyze the transfer of fucosyl moiety to different p-nitrophenyl glycopyranosides (producing p-nitrophenyl beta-d-fucopyranosyl-(1,3)-beta-d-fucopyranoside, p-nitrophenyl beta-d-fucopyranosyl-(1,3)-alpha-d-glucopyranoside, p-nitrophenyl beta-d-fucopyranosyl-(1,3)-alpha-d-mannopyranoside, and p-nitrophenyl beta-d-fucopyranosyl-(1,6)-alpha-d-galactopyranoside) and alcohols (producing methyl beta-d-fucopyranoside, ethyl beta-d-fucopyranoside, 1-propyl beta-d-fucopyranoside, 2-propyl beta-d-fucopyranoside, 1-octyl beta-d-fucopyranoside, and 2-octyl beta-d-fucopyranoside). These results indicate the possibility of utilizing this enzyme as a promising tool for enzymatic synthesis of beta-d-fucosylated molecules.

Enzymatic synthesis of oligo-D-galactofuranosides and l-arabinofuranosides: from molecular dynamics to immunological assays.

D-Galactofuranosyl-containing conjugates are ubiquitous in many pathogenic microorganisms, but completely absent from mammals. As they may constitute interesting pharmacophores, recent works have been dedicated to their preparation. Besides well-reported chemical procedures, enzymatic approaches are still limited, mainly due to the lack of the corresponding biocatalysts. Based on the similarity between chemical structures, the arabinofuranosyl hydrolase Araf51 from Clostridium thermocellum was expected to recognize both the L-Araf motif and its D-Galf analogue. Molecular dynamics and STD-NMR were firstly used to confirm this hypothesis and increase our knowledge of the active site. Interestingly, this arabinofuranosidase was not only able to hydrolyze galactosyl derivatives, but was also really efficient in catalyzing oligomerisations using p-nitrophenyl furanosides as donors. The structures of the products obtained were determined using mass spectrometry and NMR. Amongst them, all the possible regioisomers of di-arabino and -galactofuranosides were synthesized, and the ratio of each regioisomer was easily tuned with respect to the reaction time. Especially, the galactofuranobioside displaying the biologically relevant sequence beta-D-Galf-(1,6)-beta-D-Galf was enzymatically prepared for the first time. All fractions going from di- to penta-arabino- and galactofuranosides were tested for their ability in eliciting the production of TNF-alpha. Interesting immunological properties were observed with arabinofuranosides as short as three sugar residues.

Metadynamics in essential coordinates: free energy simulation of conformational changes.

We propose an approach that combines an extraction of collective motions of a molecular system with a sampling of its free energy surface. A recently introduced method of metadynamics allows exploration of the free energy surface of a molecular system by means of coarse-grained dynamics with flooding of free energy minima. This free energy surface is defined as a function of a set of collective variables (e.g., interatomic distances, angles, torsions, and others). In this study, essential coordinates determined by essential dynamics (principle component analysis) were used as collective variables in metadynamics. First, dynamics of the model system (explicitly solvated alanine dipeptide, Ace-Ala-Nme) was simulated by a classical molecular dynamics simulation. The trajectory (1 ns) was then analyzed by essential dynamics to obtain essential coordinates. The free energy surface as a function of the first and second essential coordinates was then explored by metadynamics. The resulting free energy surface is in agreement with other studies of this system. We propose that a combination of these two methods (metadynamics and essential dynamics) has great potential in studies of conformational changes in peptides and proteins.

Cold-active beta-galactosidase from Arthrobacter sp. C2-2 forms compact 660 kDa hexamers: crystal structure at 1.9A resolution.

The X-ray structure of cold-active beta-galactosidase (isoenzyme C-2-2-1) from an Antarctic bacterium Arthrobacter sp. C2-2 was solved at 1.9A resolution. The enzyme forms 660 kDa hexamers with active sites opened to the central cavity of the hexamer and connected by eight channels with exterior solvent. To our best knowledge, this is the first cold-active beta-galactosidase with known structure and also the first known beta-galactosidase structure in the form of compact hexamers. The hexamer organization regulates access of substrates and ligands to six active sites and this unique packing, present also in solution, raises questions about its purpose and function. This enzyme belongs to glycosyl hydrolase family 2, similarly to Escherichia coli beta-galactosidase, forming tetramers necessary for its enzymatic function. However, we discovered significant differences between these two enzymes affecting the ability of tetramer/hexamer formation and complementation of the active site. This structure reveals new insights into the cold-adaptation mechanisms of enzymatic pathways of extremophiles.

Molecular simulations of hevein/(GlcNAc)3 complex with weakened OH/O and CH/π hydrogen bonds: implications for their role in complex stabilization.

Carbohydrate-protein complexes are often characterized by interactions via aromatic amino acid residues. Several mechanisms have been proposed to explain these stacking-like interactions between pyranose sugars and aromatic moieties. The physical basis of these interactions is being explained as either dispersion CH/π or hydrophobic. In order to elucidate the nature of these interactions, we performed a series of molecular dynamics simulation of hevein domain (HEV32) in complex with (ß-D-GlcNAc)3. Selected OH/O and CH/π hydrogen bonds involved in carbohydrate recognition were artificially weakened in 100 ns molecular dynamics simulations. Separate weakening of either OH/O or CH/π hydrogen bonds was not sufficient to destabilize the complex. This indicates that other effects, not solely CH/π dispersion interactions, contribute significantly to the stability of the complex. Significant destabilization of complexes was reached only by simultaneous weakening of OH/O and CH/π hydrogen bonds. This also shows that classical hydrogen bonds and CH/π interactions are working in concert to stabilize this carbohydrate-protein test case.

Role of CH/pi interactions in substrate binding by Escherichia coli beta-galactosidase.

Interactions between carbohydrates and aromatic amino-acid residues are often observed in structures of carbohydrate-protein complexes. They are characterized by an orientation of the pyranose or furanose ring parallel with the aromatic ring of amino-acid residues. An important role in the formation of these complexes is supposed to be played by CH/pi interactions. This paper presents an ab initio quantum chemistry study of CH/pi interactions between beta-galactosidase from E. coli and its substrates and products. The energy stabilizing the interaction between Trp999 residue and substrate bound in the shallow binding mode was calculated at the MP2/6-31+G(d) level as 5.2kcalmol(-1) for the glucose moiety of allolactose, 2.4kcalmol(-1) for the galactose moiety of allolactose and 5.0kcalmol(-1) for the glucose moiety of lactose. The energy stabilizing the interaction between Trp568 residue and galactose in the deep binding mode was calculated as 2.7kcalmol(-1). Interaction energies at the HF/6-31+G(d) and B3LYP/6-31+G(d) levels were small or repulsive; therefore, highly correlated ab initio methods were necessary to study these interactions. These unexpectedly strong interactions give a rationale for allolactose formation and illustrate the role of the Trp999 residue. In addition, this illustrates the importance of CH/pi interactions for the function of carbohydrate-binding proteins and carbohydrate-processing enzymes.

Toward an accurate conformational modeling of iduronic acid.

Iduronic acid (IdoA), unlike most other monosaccharides, can adopt different ring conformations, depending on the context of the molecular structure. Accurate modeling of this building block is essential for understanding the role of glycosaminoglycans and other glycoconjugates. Here, we use metadynamics to predict equilibria of (1)C(4), (4)C(1) and (2)S(O) conformations of α-L-IdoA-OMe and α-L-IdoA2S-OMe. Different schemes of scaling of atoms separated by three bonds (1-4 interaction) were tested. It was found that scaling (reduction) of 1-4 electrostatic interactions significantly changes conformational preferences toward the (4)C(1) conformation. More interestingly, scaling of 1-4 van der Waals interaction favors skew-boat conformations. This shows that a minor modification of noncovalent 1-4 interactions parameters can provide a good agreement between populations of conformers of iduronic acid in water from simulations and experiments.

Modelling of beta-D-glucopyranose ring distortion in different force fields: a metadynamics study.

Modelling of carbohydrate conformations is a challenging task for force field developers. Three carbohydrate force fields, namely GLYCAM06, GROMOS 45a4 and OPLS were evaluated. Free energies of different ring conformations of beta-D-glucopyranose were calculated using metadynamics in vacuum as well as in explicitly modelled water. All three force fields model the (4)C(1) conformation as the most stable by at least 6kJ/mol, as compared to other conformations. Interconversion from the (4)C(1) to any other conformation is associated with a barrier of no lower than 26kJ/mol. The free energy surface calculated in the GLYCAM06 force field is in remarkably good agreement with the recent Car-Parrinello metadynamics study. The effect of a water environment is relatively low and analogous in all tested force fields. Namely, the presence of water stabilizes the upper-left ((3,O)B) versus bottom-right (B(3,O)) area of Stoddard's plot, relative to the situation in vacuum. Comparison of free and potential surfaces is also provided for vacuum calculations.

Continuous metadynamics in essential coordinates as a tool for free energy modelling of conformational changes.

Modelling of conformational changes in biopolymers is one of the greatest challenges of molecular biophysics. Metadynamics is a recently introduced free energy modelling technique that enhances sampling of configurational (e.g. conformational) space within a molecular dynamics simulation. This enhancement is achieved by the addition of a history-dependent bias potential, which drives the system from previously visited regions. Discontinuous metadynamics in the space of essential dynamics eigenvectors (collective motions) has been proposed and tested in conformational change modelling. Here, we present an implementation of two continuous formulations of metadynamics in the essential subspace. The method was performed in a modified version of the molecular dynamics package GROMACS. These implementations were tested on conformational changes in cyclohexane, alanine dipeptide (terminally blocked alanine, Ace-Ala-Nme) and SH3 domain. The results illustrate that metadynamics in the space of essential coordinates can accurately model free energy surfaces associated with conformational changes.

Caleosin of Arabidopsis thaliana : Effect of Calcium on Functional and Structural Properties.

A non-radioactive blot binding assay has proved the capacity of a purified recombinant form of Arabidopsis thaliana caleosin (AtClo1), a key protein of this plant oil body, to bind calcium. Calcium affected recombinant caleosin aggregation state, solubility, and electrophoretic mobility on SDS-PAGE. The effect of calcium on interfacial behavior of recombinant caleosin was studied at three interfaces: air/water (A/W), purified oil/water (O/W), and air/phosholipid/water (A/PLs/W). Recombinant caleosin was able to decrease interfacial tension (IFT) at A/W and O/W interfaces as a function of concentration and calcium, whereas no interaction was detected at the A/PLs/W interface. Effect of calcium was time dependent, and its amplitude strongly varied with the interface considered. Reconstituted oil bodies were used to prove the involvement of recombinant caleosin in their calcium-driven aggregation and coalescence. Calcium ions at concentration as low as 100 nM were able to strongly modify the shape and aggregation state of purified oil bodies, as well as their behavior within a monolayer, reflecting potentially profound changes in their structure and dynamic.

Structural properties of caleosin: a MS and CD study.

We have investigated the covalent and secondary solution structure of caleosin, a 27-kDa protein also called ATS1 or AtClo1 (At4g26740) found within Arabidopsis thaliana seed lipid bodies. The native protein was partly phosphorylated at S225. Purified bacterially expressed caleosin (recClo) was not phosphorylated; cysteine residues C221 and C230 were connected by a disulfide bridge. In solution it exists as a mixture of predominant monomers and covalent dimers. We have used recClo as a model for the study of AtClo1 secondary structure. recClo is folded in aqueous solution (16% alpha-helix, 29% beta-sheet), its secondary structure being dramatically influenced by the polarity of media, as deduced from CD spectra measured in the presence of increasing concentrations of various aliphatic alcohols.