Gel phase in hydrated calcium dipicolinate

The mineralization of dipicolinic acid (DPA) molecules in bacterial spore cores with Ca2+ ions to form Ca-DPA is critical to the wet-heat resistance of spores. This resistance to “wet-heat” also depends on the physical properties of water and DPA in the hydrated Ca-DPA-rich protoplasm. Using reactive molecular dynamics simulations, we have determined the phase diagram of hydrated Ca-DPA as a function of temperature and water concentration, which shows the existence of a gel phase along with distinct solid-gel and gel-liquid phase transitions. Simulations reveal monotonically decreasing solid-gel-liquid transition temperatures with increasing hydration, which explains the experimental trend of wet-heat resistance of bacterial spores. Our observation of different phases of water also reconciles previous conflicting experimental findings on the state of water in bacterial spores. Further comparison with an unmineralized hydrated DPA system allows us to quantify the importance of Ca mineralization in decreasi...

6][7] These spores have a wide range of lifespans ranging from a couple of days to millions of years. 8The physical and chemical resistance of these spores, their ability to germinate after sub-lethal damage, and their suitability for genetic manipulations are of great practical and commercial importance. 9acterial spores have a complex, layered structure, comprising a large and loose-fitting outermost exosporium, followed by an inner spore coat containing over 50 proteins. 2,10urther inward lies the outer membrane, which has no major contribution to spore's resistance to heat or permeability. 11he outer membrane is followed by the cortex which is made up of peptidoglycan (PG), which is responsible for the decrease in the water content inside the core.The inner most part of the spore is the spore core, and it consists of metabolic enzymes, proteins, DNA, and Dipicolinic Acid (DPA), and its mineralized forms contain divalent cations such as Ca 2þ , Mg 2þ , and Mn 2þ . 11arious techniques, including heat, chemical, and radioactive treatments, have been proposed for the sterilization of harmful bacterial spores in the defense and food industries.Recently, nitrogen cold atmospheric plasma has been used as a biocide because it causes changes in surface morphology and loss of spore integrity. 12High pressure treatment of spores has been applied, resulting in their loss of viability.This is primarily attributed to the loss of DPA or due to the possible oxidative losses in the absence of DPA. 13 Chemical treatments involving decontaminating agents like dimethyldioxirane (DMDO) have also been studied and found to initiate possible oxidative losses due to a variety of oxidation reactions. 14However, these sterilization techniques carry disadvantages such as increased energy consumption and treatment under extreme conditions associated with plasma production.Considering these factors, wet-heat treatments are considered a better alternative for killing spores. 15,16Application of wet-heat in the range of 80-95 C leads to the loss of viability in the spore core, resulting in compromised germination and no outgrowth leading to spore death.Furthermore, it has been observed that heat resistance of bacterial spore decreases with the increasing hydration content during wet-heat treatment. 17lthough the role of DPA in its mineralized form (chelated with Ca 2þ to form Ca-DPA) in affecting the wet-heat resistance has been quantified, 16,[18][19][20][21][22] the mechanism responsible for this improved heat resistance is not well understood.A leading hypothesis posits that the heat resistance is an outcome of the suppressed enzyme activity and protein agglomeration due to the extremely low diffusivity in the Ca-DPA rich spore core. 23Further, the structure of water molecules in this material system also remains unresolved, with prior experiments suggesting both mobile liquid-like and immobile glassy configurations of water molecules in the spore core.Recent nuclear magnetic resonance experiments indicate that water is immobilized in the spore cores of Bacillus subtilis and Bacillus pumilus, i.e., the spore core is glass-like with no significant diffusivity of water molecules, 24 rather than being gel-like. 15,16On the contrary, measurements of water 2 H and 17 O spin relaxation rates in D 2 O-exchanged Bacillus subtilis indicate high water mobility, indicating a gel-like configuration of water molecules. 23Here, we reconcile these findings using reactive molecular dynamics (RMD) simulations to quantify the structure and diffusivity of water molecules 25,26 as a function of temperature and water content.Our simulation results clearly distinguish between solid, gel-like, and molten states in the temperaturewater content phase diagram.Understanding molecular mobility and phase transformations in the spore core may help explain the role of the Ca-DPA rich core matrix in arresting protein conformational changes, agglomeration, and denaturation and may shed light on the wet-heat resistance of bacterial spores.
We use RMD simulations based on first principlesderived reactive force fields (ReaxFF) [27][28][29][30] to calculate the mean square displacement (MSD) and coefficient of selfdiffusion of water and DPA molecules in hydrated Ca-DPA simulation cells as a function of temperature and level of hydration.We have developed a force field for Ca-DPAÁ3H 2 O that reproduces experimental lattice constants and important bond-lengths for this system.Our main motivation to use the reactive force field was for hydrated systems.The calculated MSD values are also used to distinguish between different material phases since the MSD for a system in a solid state is constant in time for large times, while a linear dependence of MSD with time indicates a fluid phase.In RMD simulations, the time evolution of atomic trajectories is determined by an environment-dependent force field based on the concept of bond order.The ReaxFF describes chemical bond breakage and formation by introducing bond orders and charge transfer between atoms based on a charge-equilibration approach.We use an in-house scalable parallel reactive molecular dynamics (RMD) code, in which the time-consuming iterative determination of atomic charges is replaced by a computationally efficient extended Lagrangian scheme. 31e start our RMD simulations with equilibrated zeropressure simulation cells of the hydrated Ca-DPA crystal, where the water concentration ranges between 20.8 and 33.54 wt.%.Zero-pressure configurations are obtained for all hydrated Ca-DPA systems at various temperatures using the following procedure.The initial structures of Ca-DPA are first relaxed by setting atomic velocities to zero at every 100 steps for a total of 20,000 steps, where the unit time step is 0.3 fs.Subsequently, the systems are heated to À173 C and thermalized for 200,000 time steps in canonical (NVT) ensemble.These relaxed configurations are then expanded in steps of 1.5% by volume, and each of these expanded systems is heated up to 400 C in steps of 50 C.The length of each of these runs is again 200,000 steps.After that, the zero-pressure configuration for each of these expanded systems is obtained from their pressure-temperature plot.These zero-pressure systems are further relaxed for 1000 ps for the MSD calculation.It can be observed that the volume of the Ca-DPA simulation cell expands linearly with the water content, increasing by 33.5% as the water content increases from 20.The length of our MD simulations is 1000 ps.The length of the correlation function hR(t) 2 i is 60 ps; however, 600 initial conditions are used for better averages (Fig. 3).Selfdiffusion coefficients, D, for the various molecules are computed from the derivative of the MSD with time t using the following equation: (1) This quantification of constant of self-diffusion can be used to distinguish between liquid-like, gel-like, and crystallike behavior of water and DPA molecules in the spore.In our analysis, we consider Ca-DPA and water diffusivities greater than 3 Â 10 À8 cm 2 /s and 1.2 Â 10 À7 cm 2 /s to represent the onset of liquid-like mobility, based on the magnitude to measure the diffusivity of confined water molecules in silica gels. 32o better understand the physical state of the water in hydrated Ca-DPA, we monitor the trajectories of different types of water molecules in Ca-DPA at 25 C and above 100 C. The trajectories of water molecules in Figs.4(a ) in the temperature range of 250 C to 92 C leads to a spore core more amenable to processes that precede spore death.We refer to this phase as the "gel" phase, denoted by light orange color in Fig. 5.A further increase in the level of hydration leads to liquid-like diffusion in both the Ca-DPA and water bound with Ca 2þ ions and gives rise to a fully molten phase.
Most importantly for the moist-heat treatment of bacterial spores, we demonstrate that the temperature required for the solid-gel and gel-molten phase transformations decreases with increasing hydration and saturates at 92 C for the water content beyond 31 wt.% H 2 O.This observation is consistent with experimental measurements of heat resistance of spore cores in contact with wet-heat [Fig.5(b)], which shows that the heat resistance of the spore core starts decreasing when the water content in the protoplast becomes more than 30%. 17he phase diagram in Fig. 5 is dominated by phase transitions involving Ca-coordinated water molecules, H 2 O 1 and H 2 O 2 .In order to understand the effect of Ca mineralization on the phase transformations, we compare Fig. 5 to the phase diagram of the Ca-free hydrated DPA crystal (Fig. 6) at different temperatures and different water contents, specifically 0 wt.% (anhydrous DPA), 33 9.7 wt.% (DPA monohydrate, DPAÁH 2 O), 34 and 17.7 wt.% (DPA dihydrate, DPAÁ2H 2 O).
Similar to the hydrated Ca-DPA phase diagram in Fig. 5, the transition temperatures for the solid-gel and gel-liquid phase transitions decrease with increasing hydration from pure anhydrous DPA to DPA monohydrate to DPA dihydrate.Figure 5 shows the temperature for the onset of long-range diffusion of DPA and water molecules dropping from 400 C in anhydrous DPA to 260 C in DPAÁH 2 O (9.7 wt.% H 2 O) and subsequently to 180 C in DPAÁ2H 2 O (17.7 wt.% H 2 O).For a more direct illustration of the impact of Ca mineralization on the heat resistance, we compare the melting point (i.e., gel-molten phase transition temperature) of DPAÁ2H  In conclusion, we used RMD simulations to identify the structure and calculate the diffusivity of water, DPA, and Ca-DPA molecules as a function of temperature and hydration.The analysis described here has three main conclusions.First, it helps to reconcile outstanding differences in experimental characterization of water molecules in bacterial spore cores.Second, we show a direct relationship between molecular diffusivities in the spore core and its extent of hydration, which helps to explain the greater heat resistance of dry bacterial spores.Further, we also calculate the impact of Ca mineralization on heat resistance of spores by quantifying the melting points of unmineralized and Ca-mineralized hydrated DPA.This research was supported by the Defense Threat Reduction Agency Grant No. HDTRA1-14-1-0074.We thank Dr. Douglas Allen Dalton and Dr. Suhithi Peiris for their encouragement and continued support throughout this investigation.The computations were performed at the Center for High Performance Computing of the University of Southern California.
Figure 1 shows the unit cell of the Ca-DPA crystal containing four Ca ions and four DPA molecules along with twelve water molecules that are bound to Ca ions.In addition, the figure also contains 11 water molecules that represent excess added water found in Ca-DPA systems at higher hydration levels.The crystal structure of the native form of calcium dipicolinate [Ca-DPAÁ3H 2 O] 22 contains two types of H 2 O molecules, distinguishable by the number of calcium ions in their coordination shell.We refer to H 2 O molecules containing one and two calcium ions in their coordination shells as H 2 O 1 and H 2 O 2 , respectively.Wet-heat treatments lead to increased hydration of the Ca-DPA matrix, leading to a composition of [Ca-DPAÁ3H 2 O] 4 þ nH 2 O-containing n added water molecules per unit cell, denoted as H 2 O 0 , that are not bonded to any of the Ca ions in the Ca-DPA crystal.These nominally unbound water molecules occupy interstices within the hydrated CA-DPA structure.For the rest of this paper, we refer to the oxygen atoms found in H 2 O 1 , H 2 O 2 , and H 2 O 0 as O 1 , O 2 , and O 0 , respectively, and the oxygen in the DPA molecule as O DPA .We will use their RMD trajectories to monitor the mean square displacement (MSD) and determine the coefficient of self-diffusion of the three types of H 2 O and DPA molecules.

Figure 2 (
FIG. 1. Configuration of [Ca-DPAÁ3H 2 O] 4 þ 11H 2 O (water content of 33.5 wt.% H 2 O), prepared from [Ca-DPAÁ3H 2 O] 4 by adding 11 additional water molecules to the unit cell (outlined in black).O 1 (in yellow) and O 2 (in magenta) are O atoms of singly and doubly coordinated water molecules with calcium of Ca-DPA, respectively, and O 0 (in green) is the oxygen atom of the added water molecules.

FIG. 2 .
FIG. 2. (a) Volume expansion (%, relative to the volume of [Ca-DPAÁ3H 2 O] 4 ) required to obtain the zero pressure configuration in the [Ca-DPAÁ3H 2 O] 4 þ nH 2 O crystal as a function of water content at 25 C. (b) Mean-square displacement of three different types of water molecules present in a representative molten system, [Ca-DPAÁ3H 2 O] 4 þ 11H 2 O (will become 1 ns above 100 C).

Figure 3
Figure 3 shows computed diffusivities for different types of H 2 O molecules and DPA for temperatures ranging from 25 C to 320 C at different levels of hydration.Comparison of Figs.3(b)-3(d) shows that the threshold temperature for the onset of diffusion is the lowest for the mobile O 0 species followed by O 1 and O 2 species consistent with the relative degree of binding to the Ca-DPA matrix.In contrast to the relatively mobile water molecules, computed diffusivities of DPA molecules are approximately one order of magnitude lower [Fig.5(a)], reflecting their greater size, weight, and higher melting point.This quantification of constant of self-diffusion can be used to distinguish between liquid-like, gel-like, and crystallike behavior of water and DPA molecules in the spore.In FIG. 3. Diffusion constant with the increasing value for (a) O DPA , (b) O 2 , (c) O 1 , and (d) O 0 present in the [Ca-DPAÁ3H 2 O] 4 þ nH 2 O system.Dotted lines denote the value above which the system is considered to be liquid.The diffusion threshold for DPA and water is 3 Â 10 À8 and 1.2 Â 10 À7 cm 2 /s, respectively.

FIG. 4 .
FIG. 4. Snapshots of RMD simulations at (a) 25 C and (b) 400 C. Trace of the trajectory of each type of water molecules in Ca-DPA for 2.4 ps is drawn as a line.Light yellow lines denote the O 1 atoms, green lines denote the O 2 atoms, and dark yellow lines denote the O 0 atoms.Multimedia view: https:// doi.org/10.1063/1.5000394.1;https://doi.org/10.1063/1.5000394.2 2 O with 18 wt.% H 2 O (180 C) to that of [Ca-DPAÁ3H 2 O] 4 with 21 wt.% H 2 O (310 C).It is notable that the Ca-mineralized [Ca-DPAÁ3H 2 O] has a significantly higher melting point even though it has a higher water content.