ramachandran plots ( 1 ) of protein backbone torsions / are frequently used to validate 3d structures of proteins ( 2,3 ) .
good when ( preferably all ) amino acids have / values located in  allowed regions of the plot . for carbohydrates
the glycosidic torsions / (  ) are the main determinants of the 3d structure and it is straightforward to validate the quality of an experimentally determined carbohydrate structure in a similar way to protein structures .
in contrast to proteins , however , the  allowed regions on a conformational map for a given glycosidic linkage not only depend on the linked monosaccharide types , but also on the linkage type and  a feature completely different from proteins  the degree of branching of the glycan .
the number of available high quality crystal structures of carbohydrates is too limited to serve as a basis to determine the 
allowed regions for all linkage types so the accessible conformational space of carbohydrate linkages has therefore to be estimated using computational methods .
consequently , over the last 30 years a considerable effort has been put into the development of force fields that are able to predict accurately the local minima and flexibility of glycosidic torsions ( 4 ) .
the force field most frequently used to calculate conformational maps of glycosidic linkages is mm3 ( 5 ) .
a variety of methods exist to calculate the energy of a carbohydrate as a function of the glycosidic torsions / ( or // for 16 linkages ) .
a  relaxed map is obtained by systematically changing  and  in small intervals ( normally 10 ) and minimizing all degrees of freedom while restraining / using an external force .
relaxed conformational maps depend on the orientation of the exocyclic torsions in the starting conformation of the carbohydrate and therefore the calculation of an 
adiabatic map is advisable . to generate a fully adiabatic map the calculation of 3 relaxed maps
the computational cost is reduced by taking into account only gg and gt ( or gt and tg for monosaccharides that have the oh4 group in an axial orientation ) conformations for the hydroxymethyl groups and a clock- and anticlockwise orientation for the hydroxyl groups .
this reduces the number of relaxed maps to be calculated to eight resulting in a total of 10 368 conformations to minimize for each 
it is obvious that this approach is limited to disaccharides since the number of conformations to be minimized for a trisaccharide would already be 3  10 ( assuming the second linkage is searched in intervals of 30 ) , far too many for a large - scale project where thousands of maps need to be calculated . high temperature molecular dynamics simulation is a robust and efficient method to explore the accessible conformational space of carbohydrates ( 8) .
conformational free energy maps can be derived from population analysis by applying the boltzmann equation .
this approach has several advantages compared to the systematic search methods : 
 it is directly applicable to branched oligosaccharides , so that the same method can be used for disaccharides and larger oligosaccharides , the low energy conformational space only is explored and no computational time is wasted to calculate unrealistic high energy conformations , the required computational cost increases therefore only moderately with the number of atoms of the oligosaccharide , the data - flow can be easily optimised in such a way that the conformational maps are generated automatically and only minimal human interaction is required .
it is directly applicable to branched oligosaccharides , so that the same method can be used for disaccharides and larger oligosaccharides , the low energy conformational space only is explored and no computational time is wasted to calculate unrealistic high energy conformations , the required computational cost increases therefore only moderately with the number of atoms of the oligosaccharide , the data - flow can be easily optimised in such a way that the conformational maps are generated automatically and only minimal human interaction is required .
it was therefore decided that free energy maps derived from md simulations were a good first set of conformational maps with which to populate the glycomaps database .
an in - house library of carbohydrate fragments ( up to pentasaccharides ) derived from structures described in the carbbank ( 9 ) and built using the sweet - ii program ( 10 ) served as input .
the mm3 force field as implemented in the tinker suite ( ) [ for a comparision with the original mm3 implementation see ( 11 ) ] was used to calculate the trajectories at 1000 k. the length of the md simulation was 10 ns for disaccharides and 30 ns for larger oligosaccharides .
the conformational analysis tools ( cat ) software ( ) was used for data processing and analysis .
conformational maps of carbohydrate linkages can be retrieved from the database by entering the disaccharide fragment in iupac form into the search database web interface of glycomapsdb ( figure 1 ) .
maps found in the database are displayed as a list , with each entry showing a preview picture and the full structure of the carbohydrate for which the map has been calculated in extended iupac form ( figure 2 ) .
difference maps can be calculated e.g. to evaluate the influence of branching on the accessible conformational space of a linkage ( figure 3 ) .
ii ( 10)can be displayed using jmol ( ) . if experimental data for the disaccharide fragment is available in the protein data bank ( pdb ) ( 12 ) , a link is displayed , which leads to a page where / torsion values retrieved from pdb data [ using the glytorsion tool ( 13 ) ] are overlaid onto the calculated conformational map ( figure 4 ) .
some statistical data regarding the fit of the experimental data to the calculated data is also displayed .
if further information such as literature references or nmr data for the carbohydrate is available in glycosciencesdb ( 14 ) a link to the corresponding entry in that database is displayed .
monosaccharide residues can be selected from the pull down menus or entered directly into the appropriate fields .
the complete carbohydrate structure used for the calculation of the conformational map is displayed in iupac extended nomenclature . the disaccharide fragment that defines
the method used to calculate the maps is shown as well as a preview picture of the map that serves as an active link to a more detailed display of the database entry .
shown the conformational maps of the disaccharide fragment -d - glcpnac-(1 - 4)--d - manp either as part of a highly branched oligosaccharide ( b ) or as part of a linear chain ( c ) are compared .
the reduction of accessible conformational space caused by neighbouring residues in the branched oligosaccharide can clearly be seen in the difference map ( a , deep red areas represent regions where the neighbouring residue causes a strong energy penalty mainly by steric conflicts ) .
comparison of calculated conformational maps with glycosidic torsion values derived from the pdb using the glytorsion tool ( 13 ) .
the pdb database entries can be explored using the  show details link . if available
visualization of carbohydrate structures is performed using the java applet jmol or the plugin chime ( ) . diagrams and plots are generated in scalable vector graphics ( svg ) format . for browsers that can not display svg files ,
the service is hosted at and maintained by the german cancer research centre in heidelberg , germany .
currently , the glycomapsdb contains 2500 conformational maps of carbohydrate fragments originally described in the carbbank .
the direct crosslink between calculated maps and pdb data opens a very efficient route to crosscheck the quality of experimental structures as well as the quality of the maps .
in general the amount of available high quality experimental data currently available for carbohydrates in the pdb database is rather limited compared to that for proteins and for some linkages there is no experimental data available at all , so there is a clear lack of experimental reference data . in this respect , the conformational maps contained in glycomapsdb might help crystallographers to crosscheck their data before submission similar to the ramachandran plot analysis for proteins . for this purpose ,
maps from the glycomapsdb are also accessed by the carp ( carbohydrate ramachandran plot ) software ( 13 ) , where users can upload a structure in pdb file format and retrieve plots comparing the torsions present in the structure with the conformational maps .
glycomapsdb indirectly offers an interface for  data mining in the pdb database , e.g. to find carbohydrate entries with unusual glycosidic torsion values . for -d - glcpnac-(1 - 4)--d - glcpnac ,
a frequent fragment contained in n - glycans , the agreement between experimental and calculated data is remarkably good ( figure 4 ) .
more than 80% of the crystal structures have values in low energy areas of the conformational map .
but there are also some outliers that are located in  not allowed areas of the maps .
of special interest are also those structures that have  or  values of 180 ( anti conformation ) .
the  show details link would help to find the corresponding entries in the pdb . in the near future
we will provide also conformational maps calculated with force fields other than mm3 , so that the influence of different parameter sets can be investigated . an upload feature for maps will be added so that users can upload and compare their own calculated maps with the maps stored in the database .
in addition , a functionality that will allow users to calculate maps not already contained in the database is planned .
 conformational energy maps of the glycosidic linkages are a valuable resource to gain information about preferred conformations and flexibility of carbohydrates . here 
 we present glycomapsdb , a new database containing more than 2500 calculated conformational maps for a variety of di- to pentasaccharide fragments contained in n- and o - glycans . 
 oligosaccharides representing branchpoints of n - glycans are included in the set of fragments , thus the influence of neighbouring residues is reflected in the conformational maps . during refinement of new crystal structures , maps contained in glycomapsdb 
 can serve as a valuable resource to check whether the torsion values of a glycosidic linkage are located in an  
 allowed region similar to the ramachandran plot analysis for proteins . 
 this might help to improve the structural quality of the glycan data contained in the protein data bank ( pdb ) . 
 a link between glycomapsdb and the pdb has been established so that the glycosidic torsions of all glycans contained in the pdb can be retrieved and compared to calculated data . 
 the service is available at . 
INTRODUCTION
CALCULATION OF CONFORMATIONAL MAPS
ACCESSING AND ANALYSING CONFORMATIONAL MAPS
IMPLEMENTATION
DISCUSSION AND OUTLOOK