M. P ꢀe rez-S aꢀ nchez et al. / Tetrahedron 67 (2011) 7708e7712
7711
column (Teknokroma) with acetonitrile/water (80:20) as a mobile
phase at a flow of 1.0 mL/min. NMR spectra were recorded on
Bruker 500 MHz spectrometers. The structure of the enzymatically
synthesized disaccharides was assigned by protoneproton shift
correlations, carboneproton shift correlation.
phenol-b-galactopiranose was initially docked into the active cen-
tre of the enzyme to determine the position of the molecule to
simulate the glycosyleenzyme intermediate. We selected the best
pose according to geometrical (close enough to the nucleophilic
2
6
residue Glu-537) and energetic criteria (favourable energy) and
set up a covalent link between the nucleophilic glutamic residue
and the galactose. In the second docking, the glycosyleenzyme was
used to perform a docking with the NAG. The best solution was
selected taking only into account energetic criteria, and it was used
as the starting structure for the subsequent molecular dynamics
studies. Docking procedures were carried out with Autodock 4
4
.2. Solvents
Glycerol-based solvents were synthesized using the same pro-
11
cedures previously described. The alcoholysis of either epichlo-
rohydrin (symmetrical derivatives) or the appropriate glycidol
ether (unsymmetrical derivatives) gave the corresponding 1,3-
dialkoxy-2-propanols. O-methylation of some selected targets
leads to the corresponding 1,2,3-trialkoxypropanes. All the solvents
were purified by distillation and used in high pure form. The
complete listing of the solvents used in this work is shown in Fig. 1.
2
7
software, using a 60Â60Â60 grid points box around nucleophilic
glutamic residue (Glu-537), to embrace all residues implicated in
ꢁ
the mechanism, and with a grid spacing of 0.375 A. A Lamarckian
genetic algorithm with the standard parameters was selected.
4.5.3. Solvent parameterisation. The glycerol derivative (5) used for
4
.3. Enzyme activity assay
the simulation was parameterised for the GROMOS 96 43a1 force
field as it was described in our recent work. The mixture of glycerol
derivate and water was previously equilibrated 300 ps at 298 K and
1 atm before solvating the system.
Protein concentration was determined by the Bradford25
method using bovine serum albumin as the standard. Hydrolytic
activity was determined spectrophotometrically by quantification
of pNP (p-nitrophenol) liberated from correspondent p-NP-
p-nitrophenyl- -galactopyranoside) in Tris/HCl buffer 10 mM pH
.3. A sample of enzyme solution (20 L) was added to 80 L of Tris/
-Gal. The reaction mixture was
incubated for 3 min at 37 C. Absorbance was measured at 410 nm.
One enzymatic unit was defined as the amount of protein that
b
-Gal
4.5.4. Molecular dynamics. Molecular dynamics simulations were
set up with the best solution obtained from docking and the gly-
cosyleenzyme intermediate. Parameters adopted in GROMOS 96
43a1 were used for the protein, while the parameters required for
the galactoseeglutamic residue were derived in a consistent
(
7
b-D
m
m
HCl buffer to containing 5 mM pNP-
b
ꢀ
2
8
manner with the force field. NAG parameters were generated by
2
9
hydrolyses 1
m
mol of pNP-
b
-Gal per minute under the conditions
the Dundee PRODRG 2.5 Server. For all systems, an initial mini-
misation was performed with 500 steps of SD followed by 1500
steps of PolakeRibiere Conjugate Gradients (CG). The minimized
complexes were solvated in different cubic boxes, fulfilling the re-
quirement of a minimum distance of 1.2 nm between any atom of
the complex and the faces of the box. In one case, the monomer was
solvated with SPC water molecules, to reproduce the buffer con-
ditions, and in the other, a mixture of glycerol derivate and SPC
water molecules was used, with the aim to reproduce the condi-
tions for the regioselectivity change. Following, each box was
minimized with 500 steps of SD and 3000 steps of CG and equili-
brated 100 ps at 298 K and 1 atm (NPT conditions). Finally, a 5 ns
production simulation for the water system was carried out while
a 20 ns production simulation for the mixture was performed. All
MD simulations were set up with GROMACS software suite
described before.
4
.4. General procedure for transglycosylation reactions using
b-galactosidase from E. coli
A solution of 100 mg (0.17 M) p-NP-b-D-Gal (donor) and 550 mg
(
1.25 M) of N-acetylglucosamine (acceptor) in 1 mL of glycerol-
ꢀ
based solvent (2 M)-buffer mixture was pre-equilibrated to 30 C.
Afterwards, 155 mmol/min (U) of b-galactosidase from E. coli were
added to the reaction mixture. Reaction was monitored by HPLC
UVevis and final products obtained by HPLC with a light scattering
detector (ELSD) were analysed. The reaction was stopped by heat-
ꢀ
ing to 100 C for 5 min. Final reaction mixture obtained was loaded
on activated carbon (50% m/m) and Celite (50% m/m) column
3
0
(
50 cmÂ2 cm, 25 cm column height) eluted with 3 volumes of
milliQ water, 3 volumes of 5% ethanol (in water) and 3 volumes of
5% ethanol (in water). Disaccharide enriched fractions were de-
(v4.0.7).
1
Acknowledgements
termined by TLC on silica gel 60 (Merck, Darmstadt, Germany), with
isopropanol/nitromethane/water (10:9:2) as eluent with detection
This work was supported by two research projects of the
by spraying the plates with 10% aq H
2
SO
4
in methanol and heating.
Spanish MICINN (Ministerio de Ciencia e Innovaci oꢀ n de Espa n~ a)
Ten, these fractions were collected in 15% ethanol; they were
pooled and analysed by HPLC-ELSD as described above. H NMR
CTQ2009-11801 and CTQ2008-05138, and one European project
(FP-62003-NMP-SMF-3, proposal 011774-2).
1
spectrum (D
2
O) was recorded to purified disaccharide on a Bruker
5
4
4
00 MHz spectrometer.
Supplementary data
.5. Molecular Modelling
.5.1. Protein model. The X-ray structure of the b-galactosidase of E.
coli was obtained from the Brookhaven Protein Data Bank (PDB ID:
PX3). The enzyme is composed by four identical units assembled
in a tetramer. To simplify the study the monomer A was selected.
The structure was relaxed with 2000 steps of Conjugative Gradient
References and notes
1
1. Anastas, P. T.; Warner, J. C. Green Chemistry: Theory and Practice; Oxford Uni-
versity: New York, NY, 1998.
2. Anastas, P. T. Chem. Rev. 2007, 107, 2167e2168.
(
CG) energy minimisation and it was used for docking and molec-
3. Horvath, I. T.; Anastas, P. T. Chem. Rev. 2007, 107, 2169e2173.
ular dynamics simulations.
4. (a) Kates, R. W.; Clark, W. C.; Corell, R.; Hall, J. M.; Jaeger, C. C.; Lowe, I.; Mc-
Carthy, J. J.; Schellnhuber, H. J.; Bolin, B.; Dickson, N. M.; Faucheux, S.; Gallopin,
G. C.; Grubler, A.; Huntley, B.; Jager, J.; Jodha, N. S.; Kasperson, R. E.; Mabogunje,
A.; Matson, P.; Mooney, H.; Moore, B.; O’Riordan, T.; Svedin, U. Science 2001,
292, 641e642; (b) Hern aꢀ iz, M. J.; Alc ꢀa ntara, R. R.; Sinisterra, J. V. Chem.dEur. J.
4.5.2. Molecular docking. In order to describe the system properly,
2
4
a two-dockings approach was followed. In a first attempt, p-nitro-