peptide- and protein- hydration is the dominant factor in the stabilization of spatial molecular structure , in the process of protein folding by gating hydrophobic residues , and in the mechanisms of peptide and protein mediated reactions [ 14 ] .
water molecules , therefore , can be considered as an integral component of biomolecular systems with dynamic , functional , and structural roles [ 47 ] .
investigation of the structural and functional role of water molecules , bound to proteins and peptides , requires a sufficient understanding of the hydration process of their building blocks [ 1 , 2 ] .
the hydration of amino acids and their derivatives at a molecular level , therefore , is of great importance and has been extensively studied with x - ray crystallography [ 1 , 3 ] and a variety of spectroscopic techniques including multinuclear magnetic resonance spectroscopy [ 2 , 813 ] , ir and raman spectroscopy [ 1416 ] , icr mass spectrometry , and laser ablation in combination with microwave spectroscopy .
we present here , for the first time in the international literature , a comparative investigation of literature d , c , n , and o nmr and crystallographic data in order to provide a coherent hydration model of amino acids and selected derivatives at different ionization states in aqueous solution and in the crystal state .

 o nmr has received little attention in amino acid and peptide research [ 2 , 12 , 13 , 19 , 20 ] .
this neglect is due to the fact that of the three naturally occurring oxygen isotopes , only o possesses a nuclear spin ( i = 5/2 ) .
owing to its electric quadrupole moment ( qe = 2.6  10  em ) and , thus , broad line widths , and its low absolute sensitivity compared with that of h ( ~1.1  10 ) , the o- isotope is one of the more difficult to observe by nmr spectroscopy [ 12 , 13 , 21 , 22 ] .
o nmr studies , therefore , of compounds at natural abundance require high concentrations ( > 0.1  m ) and extensive signal averaging .
recording of spectra can be greatly facilitated by the use of o enriched samples [ 2327 ] .  
figure 1(a ) illustrates the natural abundance o nmr spectrum of glutamic acid , 0.1  m in o - depleted water at 40c . despite the extensive signal averaging ( number of scans ( ns ) =3  10 ) and the total experimental time of 4.2 hours
, the achievable signal - to - noise ( s / n ) ratio is very poor and practically prohibitive for the accurate determination of chemical shifts and line widths .  
figure 1(b ) illustrates the clear advantages of working with o - labelled glutamic acid ( o enrichment 1 at.% ) . 

o shieldings of various chemical functional groups are very sensitive for studying hydrogen bonding interactions because of the large chemical shift range of the o nucleus [ 12 , 13 ] .
the effect of solvent - induced hydrogen bonding interactions on (o ) of the carboxyl groups is , however , rather small compared with the substantial sensitivity of over 80  ppm to hydrogen bonding interaction of (o ) of amide and carbonyl oxygens [ 12 , 13 ] .
only a single o resonance absorption is observed for the carboxylic group since the shifts of the individual resonance absorptions (c = o ) and (oh ) are averaged out by rapid intermolecular proton transfer with protic solvents , traces of h2o , and/or through hydrogen bonding aggregates of the cooh groups in organic solvents [ 12 , 13 , 23 , 24 , 26 , 28 ] .
reuben   from dilution studies of acetic acid in 1,2-dichloroethane estimated a deshielding effect of ~12  ppm due to breaking of a hydrogen bond involving the carbonyl oxygen of the acid and a shielding effect of 6  ppm due to breaking of a oho hydrogen bond .
therefore , a total shift of only + 6  ppm is expected for the monomeric acetic acid in apolar media ( dichloroethane ) compared with the dimeric form . despite the relatively low sensitivity of the o shieldings of the carboxyl group to hydrogen bond interactions , spisni and collaborators   attempted to estimate the solvation state of the -carboxyl group of amino acids in the different ionization states . figures 2(a ) and 2(b ) show the dependence of (o ) of l - alanine and l - proline as a function of molar fraction of dmso in the ph range 7 - 8 and 12 - 13 .
since dmso can not form a hydrogen bond interaction with the carboxylate group , contrary to the case of h2o , the shielding difference of 1017  ppm between the two solvents was interpreted with the hypothesis that the carboxylate group of these amino acids is hydrated by two water molecules in aqueous solution with one hydrogen bond per carboxylate oxygen . in the acidic ph range ( figures 2(c ) , 2(d ) ) , a nonlinear behaviour of the chemical shift at high dmso molar fractions was observed . for dmso
molar fractions up to 0.6 , a linear dependence of the chemical shift was observed which , on extrapolation to 100% dmso , results in a shielding of 1517  ppm , the same as in the neutral ph .
this was interpreted with the hypothesis that two hydrogen bonds ( one to each oxygen ) are being ruptured .
when the dmso molar fraction is between 0.6 and 0.8 , it was suggested that a third molecule of water , which is hydrogen bonded to the hydroxyl hydrogen , is dissociated due to the interaction with dmso .
this might explain the deflection from linearity and the plateau - like dependence of the o shielding .
the protonated form , therefore , of the carboxyl group of the amino acids is more hydrated with an access of a bound molecule of h2o than the deprotonated form .
this conclusion is in qualitative agreement with multinuclear nmr relaxation data ( see below ) . for quadrupolar nuclei , such as d , n , and o , the longitudinal ( t1 ) and transverse ( t2 ) relaxation times are essentially due to quadrupolar interaction 
 ( 1)1t1=1t2(1+23)2f(,d ) , 
 where  is the nuclear quadrupole coupling constant .
the asymmetry parameter  varies from 0 to 1 and describes the deviation of the electric field gradient from axial symmetry , and f( , d ) is the correlation function , which depends on the rotational diffusion constant d and its relative orientation with respect to the principal axes of the field gradient tensor [ 12 , 13 ] .
when isotropic reorientation is assumed , f( , d ) reduces to a single overall correlation time mol which is given by the stokes - debye formula 
 ( 2a)mol = vmkbt , 
 where vm is the molecular volume ,  the viscosity of the solution , kb the boltzman constant , and t the absolute temperature .
vm can be estimated as 
 ( 2b)vm=0.74mwn0 , 
 where n0 is the avogadro 's number and mw and  are the molecular weight and the density of the solute ( amino acid ) , respectively . 

the removal , therefore , of these impurities is necessary in studies of t2 and t1 relaxation times .
figure 3 illustrates the ph dependence of the o line widths of 0.1  m glycine in h2o .
this o line width minimum has been previously explained by a decrease of the molecular tumbling time attributable to a reduction in hydration and , thus , intermolecular association of glycine in the zwitterionic form . in the high ph region , a broad maximum at ph 11 was observed .
addition of 2  mm ethylenediamine - n , n , n,n-tetraacetate ( edta ) to the original solution resulted in no line width variation in the neutral and high ph region .
it can , therefore , be concluded that this broad minimum at ph 11 should be attributed to the effect of paramagnetic impurities and not to a hydration change of glycine in the neutral and high ph region .  

d  t1 relaxation times of cd2 of glycine at acidic ph were shown to be shorter relative to those at neutral ph .
this shortening in t1 implies an increase in mol and , thus , in the effective molecular weight mw ( ( 1 ) , ( 2a ) , and ( 2b ) ) , which was interpreted with an increase in the hydration state in the cationic form . 

tritt goc and fiat   investigated in detail the viscosity and temperature dependence of the o nmr line width of glycine , alanine , proline , leucine , histidine , and phenylalanine at ph 2 , 7 , and 12.5 .
the experimentally observed viscosity / temperature ( /t ) dependence of the reorientation correlation time was compared with various hydrodynamic models . a model of the hydration state in the primary solvation sphere of the carboxylic group of amino acids in their cationic state
was suggested in which two water molecules are hydrogen bonded to the oxygens and one to the hydrogen of the oh group . in the zwitterionic and anionic states , the hydration model of the carboxylate group
can be presented by a structure in which one water molecule is hydrogen bonded to each of the oxygens .
the o [ 10 , 11 ] and n nmr   line widths of several protein amino acids were measured in aqueous solution to investigate the effect of molecular weight on the line widths ( table 1 ) .
the n and o line widths , under composite proton decoupling , increase with the bulk of the amino acid , and increase at low ph . assuming an isotropic molecular reorientation of a rigid sphere and , thus , a single correlation time from overall molecular reorientation ( mol ) , then
, the line width 1/2 can be expressed in the following form : 
 ( 3)1/2=1t2=0+1mw , 
 where mw is the molecular weight , 1 is the contribution to the line width of the quadrupolar coupling constant , density and temperature , ( 1 ) , ( 2a ) , and ( 2b ) , and 0 is the solvent viscosity - independent contributions to the line width due to the primary hydration sphere of the amino acids .  
the linear correlation between 1/2 and mw at ph 6 for both n and o nuclei ( figure 4 ) is in agreement with the hydrodynamic model of ( 3 ) .
furthermore , the (o ) of the amino acid is independent of both the ionization and the degree of hydration of the carboxyl group .
the increase in the o line widths at acidic ph ( ~100  31  hz ) , relative to those at neutral ph , was interpreted by a change in the rotational correlation time and , thus , effective mw of the amino acids , ( 3 ) .
this implies that the cationic form of the amino acids is more hydrated by an access of 1.3 to 2.5 molecules of water relative to that in the zwitterionic form   with lifetimes that are longer than the overall molecular rotational correlation time , presumably 210 ps . in the case of a stochastic diffusion of the amino and carboxyl groups comprising contributions from internal ( int ) and overall ( mol ) motions , the correlation time c for n or o
is given by ( 4a)c=mol[a+(b+c)(12/r)intmol+(12/r)int ] 
 with 
 ( 4b)a=34(3cos21)2 ,  b=3sin2cos2 ,  c=34sin4,where  is the angle between the rotation axis and the main field gradient ( r denotes an r - fold jump mechanism ) . since the sum of a , b , and c is equal to 1 , ( 4a )
can be rewritten as(5a)c=molamol+imol+i , 
 where 
 ( 5b)i  =  ( 12r)int.equations ( 5a ) and ( 5b ) can be rewritten as 
 ( 6)c=(1a)i+amol(1a)i2mol+i . 

since a and i can be assumed to be constant for all the amino acids , ( 4a ) and ( 4b ) can be written as 
 ( 7)1/2=0+1mw+2mw+3 , 
 where 03 are constants .
the minimization of ( 7 ) on the basis of the o experimental data gave the mean difference of 35.8  17.3 in mw between ph 0.5 and 6.0 for three different 1/2 values : 250 , 350 ( figure 4(b ) ) , and 500  hz .
the difference in the n line widths at the two ionization states ( figure 4(a ) ) should be attributed to differences in the correlation times and to a decrease in the (n ) on deprotonation of the carboxyl group . in the case of the linear model ,
the influence of variations of values of the (n ) to the line width , 1/2 , is less for small molecular weights .
therefore , for 1/2 = 70  hz ( figure 4(a ) ) , the difference in mw will be a reasonable approximation of the difference in hydration in the two states .
the calculated value was found to be 45.2  7.4 , which corresponds to an excess of 2 - 3 water molecules in the cationic form compared to that in the zwitterionic form , in reasonable agreement with the o nmr data .
more recently , takis et al .   investigated the cc longitudinal relaxation times ( t1 ) and n line widths ( 1/2 ) of amino acids and acetyl - amino acids in aqueous solutions at acid and neutral ph . both c and n values indicate that amino acids and acetyl - amino acids at acid ph interact with an access of one water molecule with respect to their deprotonated form at neutral ph . on the contrary , c and n values of betaines ( r3nch(r)coo )
crystal structure databases provide a rich source of information to extract details on the architectures and interactions of molecules .
this kind of search provides the opportunity to examine the formation of intramolecular and intermolecular hydrogen bond in small molecule crystal structures [ 33 , 34 ] .
propensities for the hydration of the -carboxylate group of amino acids and their derivatives were derived on the basis of exhaustive searches in the cambridge crystallographic database ( csd ) .
since intermolecular hydrogen bonds are preferred when five- or six - membered conjugated rings are formed , particular attention has been given to the hydrogen bond patterns in the vicinity of the carboxylate group that involves two simultaneous hydrogen acceptors .
the concept of five- and six - membered conjugated rings , along with three - center ( bifurcated ) and 4-center ( trifurcated ) hydrogen bonds , has been acknowledged and accepted widely as an important factor in determining the structure and function of molecules ranging from inorganic to organic and biological molecules [ 1 , 3539 ] .
furthermore , port and pullman   studied theoretically the formate ion - water interaction as a prototype of the carboxylate group .
three energetically favourable hydration sites were obtained , two equivalent sites on the carboxylate oxygens at the exterior of the ion and one water bridging the two oxygen atoms .
the conquest 1.13 program was used for all the statistical analysis described in this paper .
specifically , the csd version 5.32 ( november 2010 ) for small molecules was searched , with the following general search flags : r > 0.5 , 
, and  only organic  . in order to extract the number of entries present in the current database that form six - membered conjugated rings between the two oxygens of the -carboxylate and the carboxylic group with a molecule of water in the vicinity , the following geometric cut - offs were used : upper limits d = 3   for ( ow)ho = c and ( ow)ho  c , and d = 3.5   for owo = c and owo  c
44 hits were obtained for the carboxylate state ( figure 5 ) , whereas only one was derived for the protonated form .
figure 5 demonstrates that the oxygen of water , ow , is reasonably close to the carboxylate oxygens and displays a significant preference for the o1c 
there is a general correlation between hydrogen bond lengths and hydrogen bond angles ( figure 6 ) similar to that observed by jeffrey and maluszynska   in the case of water molecules in the hydrates of small biological molecules .
furthermore , crystallographic database searches were performed to identify the propensity for the formation of intramolecular hydrogen bond interaction in the carboxylate ( nho ) and the carboxylic acid ( nho = c ) state .
interestingly , 946 and 118 hits were retrieved for the carboxylate and 621 and 6 hits for the carboxylic form in the absence and presence of two molecules of water , respectively .
it is evident from figure 7 that in the presence of two bound water molecules there is a significant reduction in the number of structures with intramolecular hydrogen bond interaction for the carboxylate group and , concurrently , a significant increase in the distance ( nho ) .
it is important to note that no intramolecular nh3ooc hydrogen bonds were observed for 82 amino acid carboxylates with sp - hybridized c - atoms   in agreement with an early survey of amino acid structures determined by neutron diffraction .

 o shielding changes of amino acids as a function of molar fraction of dmso / h2o , the decrease in the longitudinal relaxation times ( t1 ) of c d and c , and the increase in line widths of n and o at acidic ph relative to those at neutral ph may be interpreted with the hypothesis that the cationic form of amino acids is more hydrated by 1 to 3 molecules of water than the zwitterionic form .
similar behaviour was also observed for acetylated derivatives of amino acids , but not for betaines , between the protonated and deprotonated carboxyl group .
although the precise hydration differences observed for various nuclei deviate somehow , it may be concluded that these hydrated complexes have lifetimes that are shorter than the nmr chemical shift time scale , but presumably longer than the overall molecular rotational correlation time of 210 ps .
an exhaustive search in the cambridge crystallographic database ( csd ) demonstrates a strong tendency of the two oxygens of the deprotonated carboxylate group to form hydrogen bonds with a single molecule of water .
even though statistical analysis of structural parameters in crystals can not be used in a straightforward way to derive quantitative structural models in solution , it is of interest to note that this mode of six - membered conjugated ring , which is absent in the case of the carboxylic group , might result in a more compact and , thus , less hydrated structure in aqueous solution , in accordance with the nmr data ( figure 8) .
furthermore , it may be concluded that the bound molecules of water alleviate the nho interaction and very probably this effect is even more pronounced in aqueous solution . from the above ,
it is evident that the reduced hydration of the carboxylate group , relative to the carboxylic group , should be attributed mainly to the strong tendency of the carboxylate group to form a six - membered conjugated ring with a single molecule of water .
constructively , the tentative models illustrated in figure 8 should be further validated by in silico and experimental approaches .
computational methods complement the experimental results by providing information on the microscope and physicochemical details on the interplay between water and the biomolecule of interest [ 4447 ] . for example , introduction of solvent effects into molecular dynamics can provide an atomic description of the folding and unfolding of a protein .
furthermore , there is an array of theoretical approaches that have been utilized for treating nmr shieldings in solution , that can be classified as continuum models [ 49 , 50 ] and molecular dynamics simulations .
experimental approaches could involve o nmr both in powders and in the crystal state   with varying degrees of hydration .
 
 2d , 13c , 14n , and 17o nmr and crystallographic data from the literature were critically evaluated in order to provide a coherent hydration model of amino acids and selected derivatives at different ionization states . 
 17o shielding variations , longitudinal relaxation times ( t1 ) of 2d and 13c and line widths ( 1/2 ) of 14n and 17o , may be interpreted with the hypothesis that the cationic form of amino acids is more hydrated by 1 to 3 molecules of water than the zwitterionic form . 
 similar behaviour was also observed for n - acetylated derivatives of amino acids . 
 an exhaustive search in crystal structure databases demonstrates the importance of six - membered hydrogen - bonded conjugated rings of both oxygens of the -carboxylate group with a molecule of water in the vicinity . 
 this type of hydrogen bond mode is absent in the case of the carboxylic groups . moreover , a considerable number of structures was identified with the propensity to form intramolecular hydrogen bond both in the carboxylic acid ( nho = c ) and in the carboxylate ( nh  o ) ionization state . in the presence of bound molecules of water 
 this interaction is significantly reduced in the case of the carboxylate group whereas it is statistically negligible in the carboxylic group . 
1. Introduction
2. Results and Discussion
3. Conclusions