Thermal Deformations of Crystal Structures in the L-Aspartic Acid/L-Glutamic Acid System and DL-Aspartic Acid

: The method of temperature-resolved powder X-ray diffraction (TRPXRD) was used to determine the elevated temperature behavior of L-aspartic acid (L-asp), DL-aspartic acid (DL-asp), L-glutamic acid (L-glu), and an L-asp 0.25 ,L-glu 0.75 solid solution. These amino acids were not found to undergo any solid-phase (polymorph) transformations. When heated, they all experienced only thermal deformations. The corresponding parameters of the monoclinic cells of L-asp and DL-asp, and the orthorhombic cells of L-glu and L-asp 0.25 ,L-glu 0.75 , were calculated for the entire range of studied temperatures (up to 220 ◦ C). The data obtained were used to calculate the parameters of the thermal deformation tensors, and to plot the ﬁgures of their thermal expansion coefﬁcients. A correlation between the maximum and minimum values of thermal expansion coefﬁcients and the length, type, direction, and number of hydrogen bonds in the crystal structures of the investigated amino acids was established. The observed negative thermal expansion (contraction) of crystal structures of L-asp and DL-asp along the ac plane can be explained as a result of shear deformations occurring in monoclinic crystals with a non-ﬁxed angle β . The studies were related to the presence of amino acids in various natural and technological processes occurring at different temperatures.


Introduction
Amino acids of the general formula NH 2 -CHR-COOH play an important role in natural biochemical processes. They are typical representatives of organic crystals composed of chiral molecules. Chirality is the ability of a substance to occur in two or more configurations of molecules, which are referred to as enantiomers and diastereomers. Molecules of enantiomers are mirror images of one another. Relative configurations of enantiomers are designated as L (levorotary) and D (dextrorotary). All of the naturally occurring protein molecules are composed of L-enantiomers of amino acids. The respective D-amino acids are less frequently found in nature; for example, they are present in the cellular walls of some bacteria, in some antibiotic drugs, etc. Amino acids are involved in the regulation of metabolic processes; they also stimulate oxidation-reduction reactions in the brain and perform neurotransmission functions in the central nervous system [1]. Moreover, amino acids are widely used in the pharmaceutical industry, and participate in many technological processes occurring at different temperatures [2][3][4]. The chirality of amino acids and their wide distribution in the geological media make these compounds useful tools for dating sediments and paleontological objects [5][6][7]. Considering all of the above-mentioned, investigation of amino acids' behavior at elevated temperatures is a subject of high relevance.
Thermal deformations in organic compounds [8,9] are considerably less studied than those in inorganic substances [10][11][12], while thermal deformations of chiral organic Thermal deformations in organic compounds [8,9] are considerably less studied than those in inorganic substances [10][11][12], while thermal deformations of chiral organic compounds are hardly investigated at all. Surprisingly enough, there are only a few reported works discussing thermal deformations of amino acids and presenting figures of their thermal expansion coefficients (CTEs). Here, it is worth mentioning the studies of K. Nakata et al. [13], S.J. Coles et al. [14], and C.H. Görbitz et al. [15,16], which discuss the high-temperature polymorph transformations of L-and DL-methionine, L-norvaline, Land DL-2-aminobutyric acid, and L-and DL-norleucine.
In our previous work, we investigated the thermal deformations of amino acid crystal structures and plotted the figures of CTEs for the following substances: (1) both components and two solid solutions formed in the L-threonine/L-allo-threonine diastereomer system [17], (2) both components and the non-equimolar discrete compound formed in the L-valine/L-isoleucine system [18], and (3) the components formed in the L-alanine/L-serine amino acids system [19].
The present work is seen as continuation of the previous investigations of thermal deformations in amino acid crystal structures, aimed at a better understanding of the structure-property relationships of this important class of compounds. In the system of levorotary enantiomers L-aspartic acid/β-L-glutamic acid, both components and the solid solution L-asp0.25,L-glu0.75 (L-asp/L-glu = 25/75 mol.%) were studied. In addition, racemic aspartic acid (DL-asp) was examined.

Materials
Aspartic acid (C4H7NO4) and glutamic acid (C5H9NO4) are proteinogenic aliphatic amino acids. They are characterized by having two carboxylic COOH groups in their molecules, and differ only in the absence or presence of a CH2 group (Figure 1). Both acids have only one chiral center, i.e., they can exist in two different configurations, referred to as L and D enantiomers. In the event that the enantiomers form a binary equimolar compound, the D and L molecules contribute equally to the formation of the resulting DL crystalline structure. The crystals of L-asp and DL-asp show monoclinic syngony with the space groups P21 and C2/c, respectively [20,21]. The L-glu enantiomer possesses orthorhombic syngony with the space group P212121 [22]. The reactants used were as follows: (1) L-and DL-aspartic acid of 99% purity, obtained from AppliChem GmbH, Germany, and (2) β-L-glutamic acid of 99% purity, obtained from Acros Organics-Thermo Fisher Scientific, Belgium. The samples of the L-asp and L-glu enantiomers, DL-asp, and the solid solution L-asp0.25,L-glu0.75, applied in the The reactants used were as follows: (1) L-and DL-aspartic acid of 99% purity, obtained from AppliChem GmbH, Germany, and (2) β-L-glutamic acid of 99% purity, obtained from Acros Organics-Thermo Fisher Scientific, Belgium. The samples of the L-asp and L-glu enantiomers, DL-asp, and the solid solution L-asp 0.25 ,L-glu 0.75 , applied in the experiments, were obtained by dissolving the corresponding compounds in distilled water at~70 • C, followed by crystallization via isothermal evaporation of the solvent at 55 • C for 48 h.

Methods
The above-mentioned samples were pre-investigated by means of the PXRD method (Rigaku MiniFlex II, Tokyo, Japan). The behavior of the samples at elevated temperatures was studied using the temperature-resolved powder X-ray diffraction (TRPXRD) technique. The experiments were performed in atmospheric air using a Rigaku Ultima IV (Tokyo, Japan) diffractometer equipped with a high-temperature accessory. The experimental conditions were as follows: Co Kα irradiation, 2θ range of 5-50 • , and temperature range of 23-220 • C, with a temperature step of 10 • C. The unit cell parameters at each temperature were calculated via the Rietveld method using Topaz software (Bruker AXS GmbH, Karlsruhe, Germany).
Temperature dependencies of the unit cell parameters were applied to calculate the thermal expansion coefficients of the crystal structure (α) along the three mutually orthogonal axes of the thermal deformation tensor (α 11 , α 22 , α 33 ), and along the crystallographic axes a, b, and c. The resulting values were utilized to plot projections of the CTE figures onto the ab, ac, and bc planes of the orthorhombic (L-glu and L-asp 0.25 ,L-glu 0.75 ) and monoclinic (L-asp and DL-asp) unit cells. In the case of crystal structures with orthorhombic syngony, the tensor axes were found to be co-directional with the crystallographic axes. The tensor calculations and the plotting of the CTE figures were performed using TEV V1.0.1 software (Thomas Langreiter and Volker Kahlenberg, Institute of Mineralogy and Petrography, Innsbruck, Austria).           As the temperature increases, the peaks are shifted, usually towards the lower 2θ angles. Since the shift values of peaks with different hkl are not the same, some peaks merge.

Temperature-Resolved Powder X-ray Diffraction (TRPXRD) Data
The sample of the L-glu enantiomer ( Figure 4) is already a mixture of two phases at the starting temperature. It is related to the polymorph modifications of β-L-glu (~86 wt.%) and α-L-glu (~14 wt.%). The intensity of the peaks corresponding to the α-L-glu component persists up to the melting point. This verifies that under the conditions of the TRPXRD experiment the β-L-glu and α-L-glu phases do not undergo solid-phase polymorph transformations to change into one another.
The X-ray patterns of the L-asp 0.25 ,L-glu 0.75 solid solution ( Figure 5) correspond to those of the β-L-glu enantiomer. It is worth mentioning that the solid solution sample did not contain the α-L-glu admixture phase. This may be due to the fact that the conditions of preparing the solutions of L-glu and L-asp 0.25 ,L-glu 0.75 (the temperature, degree of supersaturation, and mode of stirring) were slightly different. Several attempts to crystallize a pure α-L-glu phase using the method proposed in [22][23][24] were unsuccessful.
The TRPXRD results were used to calculate the parameters and volume of monoclinic cells of L-asp and DL-asp and orthorhombic cells of L-glu and L-asp 0.25 ,L-glu 0.75 within the whole range of studied temperatures. Figure 6 shows changes in the a, b, c, and β parameters, and in the volume V of the monoclinic cells of L-and DL-asp versus temperature. In all cases, the errors in calculating the parameters and volume did not exceed the marker size. For purposes of comparison, the corresponding dependencies of the parameters and the volume of L-asp ( Figure 6a) and DL-asp ( Figure 6b) are presented in the same scale. It can be seen that the linear parameters a, b, and c, as well as the volumes V, grow linearly with the increase in temperature for both L-and DL-asp. However, the monoclinic angle β behaves differently: it decreases linearly for L-asp and increases linearly for DL-asp. The CTE values α V of the volumes calculated for the monoclinic cells of the said compounds are close to one another (Table 1).   The temperature-derived changes in the orthorhombic cell parameters a, b, c, and β and volume V of L-glu and the L-asp0.25,L-glu0.75 solid solution are presented in Figure  7a,b, respectively-again at the same scale. For both the enantiomer and the solid solution, increasing the temperature caused a linear increase in all of the parameters and in the volume of the respective orthorhombic cells. This confirms that heating the samples results only in the thermal expansion of their crystal structures. As seen in Figure 7, elevating the temperature leads to a greater increase in the parameters and the volume of the unit cell for L-glu than for the L-asp0.25,L-glu0.75 solid solution. Accordingly, the CTE values αV calculated for these samples appear to be rather different ( Table 2).  The temperature-derived changes in the orthorhombic cell parameters a, b, c, and β and volume V of L-glu and the L-asp 0.25 ,L-glu 0.75 solid solution are presented in Figure 7a,b, respectively-again at the same scale. For both the enantiomer and the solid solution, increasing the temperature caused a linear increase in all of the parameters and in the volume of the respective orthorhombic cells. This confirms that heating the samples results only in the thermal expansion of their crystal structures. As seen in Figure 7, elevating the temperature leads to a greater increase in the parameters and the volume of the unit cell for L-glu than for the L-asp 0.25 ,L-glu 0.75 solid solution. Accordingly, the CTE values α V calculated for these samples appear to be rather different ( Table 2). 79 (2) All of the temperature dependencies of the parameters and the volumes of the unit cells of L-asp, DL-asp, L-glu, and L-asp0.25,L-glu0.75 were approximated by polynomials of the first degree. The calculated data served as the basis for estimation of the CTE values of the crystal structures along the axes of the thermal deformation tensor and the crystallographic axes, given for L-asp and DL-asp in Table 3, and for L-glu and the L-asp0.25,L-glu0.75 solid solution in Table 4. Table 3. Thermal expansion coefficients (α * 10 −6 °C −1 ) for the monoclinic crystal structures of L-asp and DL-asp along the axes of the thermal deformation tensor, and along the crystallographic axes.  All of the temperature dependencies of the parameters and the volumes of the unit cells of L-asp, DL-asp, L-glu, and L-asp 0.25 ,L-glu 0.75 were approximated by polynomials of the first degree. The calculated data served as the basis for estimation of the CTE values of the crystal structures along the axes of the thermal deformation tensor and the crystallographic axes, given for L-asp and DL-asp in Table 3, and for L-glu and the L-asp 0.25 ,L-glu 0.75 solid solution in Table 4. Table 3. Thermal expansion coefficients (α * 10 −6 • C −1 ) for the monoclinic crystal structures of L-asp and DL-asp along the axes of the thermal deformation tensor, and along the crystallographic axes.   The data shown in Table 3 indicate the presence of negative (or anomalous) thermal expansion-that is, contraction-of the L-asp and DL-asp monoclinic crystal structures. In both substances, the contraction occurs in the crystallographic plane ac, but in L-asp it acts along the α 11 axis of the tensor, while in DL-asp it acts along the α 33 axis.

Sample
The data presented in Tables 3 and 4 were used to plot projections of the CTE figures onto the crystallographic planes ab, ac, and bc of L-and DL-asp (Figures 8 and 9, right), and L-glu and L-asp 0.25 ,L-glu 0.75 (Figure 10, right, solid black and red dashed lines, respectively). Figures 8 and 9 (left) also show projections of the monoclinic structures of L-and DL-asp onto the a*b, ac, and bc* planes, which are orthogonal to the c, b, and a axes, respectively. The analogous Figure 10 Table 4. Thermal expansion coefficients (α * 10 −6 °C −1 ) for the orthorhombic crystal structures of L-glu and L-asp0.25,L-glu0.75 along the axes of the thermal deformation tensor, and the crystallographic axes.
The data presented in Tables 3 and 4 were used to plot projections of the CTE figures onto the crystallographic planes ab, ac, and bc of L-and DL-asp (Figures 8 and 9, right), and L-glu and L-asp0.25,L-glu0.75 (Figure 10, right, solid black and red dashed lines, respectively). Figures 8 and 9 (left) also show projections of the monoclinic structures of L-and DL-asp onto the a*b, ac, and bc* planes, which are orthogonal to the c, b, and a axes, respectively. The analogous Figure 10     green dashed lines. The crystal structure projections were plotted using the structural data from CSD (LASPRT) [25].

Discussion
In the following section, first, the results for L-asp and DL-asp crystallizing in a monoclinic lattice are discussed (Figures 8 and 9), followed by the findings for L-glu and the L-asp/L-glu solid solution crystallizing in an orthorhombic syngony ( Figure 10).

Monoclinic Amino Acids L-asp and DL-asp
The molecules in the L-asp crystal structure (Figure 8) are arranged in zigzag chains extended along the direction close to that of the bisector of the acute angle a^c. The chain

Discussion
In the following section, first, the results for L-asp and DL-asp crystallizing in a monoclinic lattice are discussed (Figures 8 and 9), followed by the findings for L-glu and the L-asp/L-glu solid solution crystallizing in an orthorhombic syngony ( Figure 10).

Monoclinic Amino Acids L-asp and DL-asp
The molecules in the L-asp crystal structure (Figure 8) are arranged in zigzag chains extended along the direction close to that of the bisector of the acute angle aˆc. The chain units are interconnected via the strongest hydrogen bond in the structure-namely, O 4 -H 1 . . . O 1 = 2577 Å-which is highlighted in the projection of the L-asp crystal structure onto the ac plane. This bond is essentially co-directional with the direction of the negative thermal expansion (α 33 ). There are two other hydrogen bonds-N-H 5 . . . O 2 and N-H 6 . . . O 3formed in approximately the same direction and within the same plane. The maximum thermal expansion is directed along the α 33 axis of the thermal deformation tensor, which approximates the bisector of the obtuse angle aˆc, i.e., the β angle; only one hydrogen bondnamely, N-H 4 . . . O 2 -is extended in that direction (Figure 8, projection ac). A sufficiently great value of the CTE is observed along the α 22 axis (which coincides with axis b), which can be explained by the relatively low concentration of hydrogen bonds along this direction compared to the crystallographic axis a (Figure 8, projection ab).
Therefore, at elevated temperatures, the L-asp crystal structure naturally contracts along the direction in which several hydrogen bonds are grouped, including the strongest hydrogen bond: O 4 -H 1 . . . O 1 . This contraction is accompanied by a diminution in the β angle (see Figure 6a), and is characterized by maximum anisotropy of the thermal deformations within the ac plane ( Figure 8).
The crystal structure of DL-asp is also characterized by negative thermal expansion in the ac plane ( Figure 9). In contrast to L-asp, however, the racemate structure contracts along the direction of the bisector of the obtuse angle aˆc, i.e., the β angle (Figure 9, projection ac). This may be due to the fact that this is also the direction of the strongest hydrogen bond: O 4 -H 7 . . . O 2 = 2558 Å. The closely located hydrogen bonds of the N-H . . . O type may also facilitate the contraction. Therefore, in this case, contraction with increase in temperature is accompanied by enlargement of the β angle.
The greatest thermal expansion of the DL-asp crystal structure is observed in the direction of crystallographic axis b (α b = α 22 ). This may be the result of molecular layers formed from interconnected chains of molecules that are arranged in the direction orthogonal to the b axis (within the ac plane). Each D-or L-asp molecule located within the layer is bonded to two molecules located in the previous and subsequent layers via N-H 2 . . . O 1 and N-H 3 . . . O 1 hydrogen bonds. This interlayer interaction is not very strong (Figure 9, left, projection bc*) and, therefore, the maximum expansion of the DL-asp crystal structure is directed along the b axis.
In the crystal structures of the aspartic acid enantiomer and the racemate, the negative thermal expansion in the ac plane is co-directed with the direction of the strongest hydrogen The effect of changing the angle β in monoclinic crystals with elevation of temperature is described in detail by S.K. Filatov et al. [26,27]. Such thermal behavior is typical of crystalline structures, where either one angular parameter (monoclinic cell) or three angular parameters (triclinic cell) are not symmetrically fixed. In addition to thermal expansion, such crystalline structures can also undergo shear deformations characterized by shifts in the molecular rows relative to one another. Figure 11 illustrates simulated shear deformations in monoclinic crystal models projected for the ac plane. Increasing the β 1 angle of the initial orthorhombic structure (Figure 11a, solid line) to reach the β 2 value (dashed line) is accompanied by contraction of the crystal structure along the bisector of this angle, and by its simultaneous expansion in the orthogonal direction.
Similarly, decreasing the initial β 1 angle (Figure 11b, solid line) to the smaller value β 2 (dashed line) also leads to contraction of the crystal structure; however, this contraction occurs along the direction orthogonal to the β angle bisector, and at the same time the structure is expanded in the direction parallel to this bisector. Figure 11a,b also show the CTE figures and the regions corresponding to the positive (axis α 11 ) and negative (axis α 33 ) thermal expansion.
Crystals 2021, 11, x FOR PEER REVIEW 12 of 14 and at the same time the structure is expanded in the direction parallel to this bisector. Figure 11a,b also show the CTE figures and the regions corresponding to the positive (axis α11) and negative (axis α33) thermal expansion. Figure 11. Schematic representation of shear deformation in relation to the figures of the thermal expansion coefficients (α11 and α33) for increased (a) and decreased (b) angle β. Angles β1 and β2 are the monoclinic angles before and after the shear deformation, respectively. The image was plotted with partial use of the data reported in [26].
A decrease in the β angle with increase in temperature was also observed for enantiomers of two other monoclinic amino acids-namely, L-valine and L-isoleucine [9]. However, in these cases the changes in the crystal structures at elevated temperatures were not accompanied by negative thermal expansion, as revealed for L-and DL-aspartic acid.

Orthorhombic Amino Acids L-glu and L-asp0.25,L-glu0.75
The strongest hydrogen bond in the orthorhombic crystal structure of L-glu is O4-H1…O2 = 2519 Å ( Figure 10, left, highlighted), which bonds the molecules into zigzagged chains extended along the crystallographic axis с. This explains the minimal CTE value observed along the с axis compared to the other crystallographic directions (Figure 10, right, solid black line). The chains are joined to one another by means of a hydrogen bond (N-H5…O1) to form layers arranged orthogonally to the b axis. The layers, in turn, are bonded via two other hydrogen bonds, i.e., N-H6…O3 and N-H9…O1. The maximal thermal expansion observed along the b axis can be explained by the fact that the direction of the hydrogen bonds uniting the molecular layers deviates to a considerable extent from direction of the b axis. The most significant anisotropy of the thermal deformation was noticed within the ас and bc planes (Figure 10  . Schematic representation of shear deformation in relation to the figures of the thermal expansion coefficients (α 11 and α 33 ) for increased (a) and decreased (b) angle β. Angles β 1 and β 2 are the monoclinic angles before and after the shear deformation, respectively. The image was plotted with partial use of the data reported in [26].
A decrease in the β angle with increase in temperature was also observed for enantiomers of two other monoclinic amino acids-namely, L-valine and L-isoleucine [9]. However, in these cases the changes in the crystal structures at elevated temperatures were not accompanied by negative thermal expansion, as revealed for L-and DL-aspartic acid. in the thermal changes to their orthorhombic cell parameters and, accordingly, different CTE values, can be caused by a different degree of imperfection of their crystal structure. Substitution of some of the L-glu molecules with smaller L-asp molecules results in the formation of defects, such as voids, in the crystal structure of the solid solution. Expansion of its crystal structure at elevated temperatures can be partially compensated for by the presence of these voids. For this reason, the intensity of changes in the L-asp 0.25 ,L-glu 0.75 orthorhombic parameters is considerably lower than that for the L-glu enantiomer.
The above results were compared with similar studies of thermal deformations for (1) the orthorhombic components in the L-threonine/L-allo-threonine (L-thr/L-athr) diastereomer system, and two solid solutions of both-namely, L-thr 0.34 ,L-athr 0.66 and Lthr 0.10 ,L-athr 0.90 [17]-and (2) the orthorhombic components in the L-alanine (L-ala) and L-serine (L-ser) enantiomer system [19]. It appears that heating resulted in increases in all orthorhombic parameters of L-glu and those of the solid solutions L-asp 0.25 ,L-glu 0.75 , L-thr 0.34 ,L-athr 0.66 , and L-thr 0.10 ,L-athr 0.90 , and in diminution of one of the parameters of the diastereomers L-thr and L-athr and the enantiomers L-ala and L-ser. Increasing all of the parameters indicates the predominantly thermal expansion of the crystal structure. Decreasing one of the parameters allows, in addition, the use of the concept of the hinge mechanism of thermal deformations [26].

Conclusions
It was established that, when heated, the L-enantiomers of aspartic and glutamic acid, the aspartic acid racemate, and the L-asp 0.25 ,L-glu 0.75 solid solution do not undergo polymorphic transitions, but experience thermal deformations. The mode of the thermal deformations observed in all of the above crystal structures was consistent with the distance, type, geometry, and concentration of the intermolecular hydrogen bonds.
Thermal deformations of L-asp and DL-asp are characterized by the manifestation of negative thermal expansion-i.e., contraction-of their crystal structures in the monoclinicity plane ac. Explanation of this effect of contraction is based on the contribution of shear deformations to thermal deformations of the crystal structure.
Thermal deformations of L-glu and the L-asp 0.25 ,L-glu 0.75 solid solution are characterized by increases in all of the parameters of their orthorhombic cells. The CTE values obtained for the solid solution were lower than those calculated for the L-glu enantiomer. The reason for this difference seems to be imperfection of the solid solution crystal structure due to isomorphic substitution of some L-glu molecules by the smaller-sized L-asp molecules.
The results of the present study were compared with data for other amino acids-such as L-valine, L-isoleucine, L-threonine, L-allo-threonine, and solid solutions of the latter two-that we had investigated previously [9,17,19].
The results obtained made it possible to reveal the individual specific features of the amino acids' thermal behavior, which is important for their efficient application.