Predicting initial dissolution rates using structural features from molecular dynamics simulations

Initial and long-term dissolution rates of aluminosilicate glasses enriched with Ti, Zr and Nd

Structural origin of high bioactivity in zirconia containing bioactive glasses.

The vitrification of radioactive waste within glass and subsequent disposal within a geological disposal facility (GDF) requires a comprehensive understanding of the effect of glass dissolution on GDF performance. This paper aims to analyse the effect of both high-level and intermediate-level waste (HLW and ILW) glass dissolution source terms on radionuclide release into the geosphere just above the disposal vault (the ‘crown’). Radionuclide migration was simulated in GoldSim for HLW in either granite or clay host rocks with a bentonite buffer using carbon steel or copper canisters, whereas ILW simulations considered either granite or clay host rocks, in either bentonite buffer or cement backfill, using concrete or cast-iron canisters. Glass dissolution source terms were varied by coupling GoldSim and MATLAB to modify the initial, residual, and resumption dissolution rates of the glass or by applying the analytical GRAAL model to glass dissolution. HLW glass results indicate no preference of granite over clay host rocks for a given canister type but that a copper canister is preferable to steel. ILW results suggest that a granite–bentonite–cast-iron environment yields lowest crown activities with cast-iron preferable to concrete as the canister, bentonite preferable to cement as the buffer/backfill, and granite preferable to clay as the host rock. Varying glass dissolution source terms (initial, residual, and resumption dissolution rates) had an understood effect on radionuclide migration, although changes were arguably insignificant considering peak crown activity for both HLW and ILW.

The initial dissolution rate of a series of multicomponent glasses is studied in order to discuss the influence of increasing magnesium content in the glass on this alteration regime and to highlight differences in behavior between calcium- and magnesium-bearing glasses. The application of MD-based topological constraint theory (TCT) is confronted to glass transition temperature (Tg) and initial dissolution rate (r0) on a glass series containing the main oxides of a French nuclear glass (AVM). In addition, a comparison between a reference magnesium-containing nuclear waste glass, AVMV4 and a proposed derived simplified composition N19M8 is performed regarding r0 values. Results indicate a similar behavior in this alteration regime for the two glasses, suggesting that this simple glass might be a good analogue to the complex one. Substituting calcium for magnesium decreases the initial dissolution rate by a factor two in the series, while an overall increase of magnesium leads to an increased dissolution rate. Analyses performed with TCT suggests that magnesium environment is better defined than calcium or sodium and may behave as an intermediate species. Finally, a correlation between the number of constraints per atom and Tg is established, while the model failed to link structural features to r0.

Dissolution rate of a basalt glass in silica-rich solutions: Implications for long-term alteration

Calcium-dependent lectin I from Pseudomonas aeruginosa (PA-IL) binds specifically to oligosaccharides presenting an alpha-galactose residue at their nonreducing end, such as the disaccharides alphaGal1-2betaGalOMe, alphaGal1-3betaGalOMe, and alphaGal1-4betaGalOMe. This provides a unique model for studying the effect of the glycosidic linkage of the ligands on structure and thermodynamics of the complexes by means of experimental and theoretical tools. The structural features of PA-IL in complex with the three disaccharides were established by docking and molecular dynamics simulations and compared with those observed in available crystal structures, including PA-IL.alphaGal1-2betaGalOMe complex, which was solved at 2.4 A resolution and reported herein. The role of a structural bridge water molecule in the binding site of PA-IL was also elucidated through molecular dynamics simulations and free energy calculations. This water molecule establishes three very stable hydrogen bonds with O6 of nonreducing galactose, oxygen from Pro-51 main chain, and nitrogen from Gln-53 main chain of the lectin binding site. Binding free energies for PA-IL in complex with the three disaccharides were investigated, and the results were compared with the experimental data determined by titration microcalorimetry. When the bridge water molecule was included in the free energy calculations, the simulations predicted the correct binding affinity trends with the 1-2-linked disaccharide presenting three times stronger affinity ligand than the other two. These results highlight the role of the water molecule in the binding site of PA-IL and indicate that it should be taken into account when designing glycoderivatives active against P. aeruginosa adhesion.

Flibanserin (FLB) is a drug used for female hypotensive sexual desire disorder approved by the FDA in August 2015. FLB exhibits extensive hepatic first-pass metabolism and low aqueous solubility, hence poor oral bioavailability. In this study, beta hydroxypropyl cyclodextrin-FLB inclusion complexes were incorporated into orally fast dissolving films. This dosage form was expected to improve FLB aqueous solubility, which would give fast onset of action and decrease presystemic metabolism, hence improving oral bioavailability. The inclusion complex at a ratio of 1:1 was prepared by the kneading method. Differential scanning calorimetry (DSC), Fourier transform infrared spectroscopy (FTIR), and powder X-ray diffractometry (XRD) were used to confirm complex formation. The Box–Behnken design (15 different formulae of FLB fast-dissolving oral films (FLBFDOFs) were utilized for the optimization of the prepared films. The Expert Design 11 program was utilized to examine the effects of three selected factors, polymer concentration (X1), plasticizer concentration (X2), and disintegrant concentration (X3) on four responses: disintegration time (DT), initial dissolution rate (IDR), dissolution efficiency (DE), and film quality (QF). Numerical optimization was performed by minimizing disintegration time (Y1), while maximizing the initial drug dissolution rate (Y2), dissolution efficiency (Y3), and the quality factor (Y4). The statistical analysis showed that X1 has a significant positive effect on the disintegration time and a significant negative effect on IDR. While X2 and X3 produced a nonsignificant negative effect on IDR. Dissolution efficiency was maximized at the middle concentration of both X2 and X3. The best film quality was observed at the middle concentration of both X1 and X2. In addition, increasing X3 leads to an improvement in film quality. The optimized film cast from an aqueous solution contains hydroxypropyl cellulose (2%) as a hydrophilic film-forming agent and propylene glycol (0.8%) as a plasticizer and cross povidone (0.2%) as a disintegrant. The prepared film released 98% of FLB after 10 min and showed good physical and mechanical properties. The optimized formula showed a disintegration time of 30 s, IDR of 16.6% per minute, DE15 of 77.7%, and QF of 90%. This dosage form is expected to partially avoid the pre-systemic metabolism with a fast onset of action, hence improving its bioavailability that favors an advantage over conventional dosage forms.

Atomic-Level Protein Structure Refinement Using Fragment-Guided Molecular Dynamics Conformation Sampling

Phospholemman (PLM) is a 72-residue bitopic cardiac transmembrane protein, which acts as a modulator of the Na(+)/K(+)-ATPase and the Na(+)/Ca(2+) exchanger and possibly forms taurine channels in nonheart tissue. This work presents a high resolution structural model obtained from a combination of site-specific infrared spectroscopy and experimentally constrained high throughput molecular dynamics (MD) simulations. Altogether, 37 experimental constraints, including nine long range orientational constraints, have been used during MD simulations in an explicit lipid bilayer/water system. The resulting tetrameric alpha-helical bundle has an average helix tilt of 7.3 degrees and a crossing angle close to 0 degrees . It does not reveal a hydrophilic pore, but instead strong interactions between various residues occlude any pore. The helix-helix packing is unusual, with Gly(19) and Gly(20) pointing to the outside of the helical bundle, facilitating potential interaction with other transmembrane proteins, thus providing a structural basis for the modulatory effect of PLM on the Na(+)/K(+)-ATPase. A two-stage model of interaction between PLM and the Na(+)/K(+)-ATPase is discussed involving PLM-ATPase interaction and subsequent formation of an unstable PLM trimer, which readily interacts with surrounding ATPase molecules. Further unconstrained MD simulations identified other packing models of PLM, one of which could potentially undergo a conformational transition to an open pore.

Use of orthophosphate complexing agents to investigate mechanisms limiting the alteration kinetics of French SON 68 nuclear glass

Quantitative structure-property relationship (QSPR) analysis is a promising approach to correlate structural features with properties of glass materials that lack long-range order and usually have complex structures. By using carefully chosen descriptors based on structural models generated from molecular dynamics (MD) simulations, correlations with properties and insights on glass behaviors can be obtained. Zirconia can significantly alter glass properties including chemical durability, even in a small amount, and hence plays an important role in vitrification of nuclear waste where long-term chemical durability is desired. In this study, borosilicate glasses with the composition of xZrO2-(61 - x)SiO2-17B2O3-18Na2O-4CaO with x = 0, 1, 2, 4, 6, and 8 were simulated using classical MD simulations with the recently developed composition-dependent potentials. Short-range (e.g., bond distance and coordination numbers) and medium-range (e.g., Qn distribution, network connectivity, and ring-size distribution) structural features altered by ZrO2 were obtained and analyzed. The use of a descriptor ( Fnet descriptor) that combines short-range structural characteristics, from MD simulations, and the cation-oxygen single bond strength was found to provide excellent linear correlations with the density and initial dissolution rate of these glasses. The results show that by combining MD simulations and QSPR analysis the composition and structural effect on the properties of complex multicomponent glasses can be elucidated, thus suggesting that this is a promising approach for future glass research and new composition design.

Structural features of ISG borosilicate nuclear waste glasses revealed from high-energy X-ray diffraction and molecular dynamics simulations

Evaluating the corrosion behavior of borosilicate glass after irradiation holds vital importance for the safe geological disposal of high‐level radioactive waste. This study examined the leaching properties of three kinds of sodium‐borosilicate (RNa2O‒B2O3‒KSiO2, denoted as “NBS”) glasses with the same R but varying K under 5 MeV Xe‐ion radiation to decipher the correlation between the corrosion mechanism and glass composition. After Xe‐ion radiation, an increase in [BO3] units was observed, indicative of the glass network depolymerization. Irradiated glasses exhibited a more significant change in [BO3] units over the leaching time compared to their pristine counterparts. There is a strong correlation between the initial dissolution rates of pristine glasses, the K value, and the proportion of Na. Additionally, the initial Si dissolution rate of the NBS glasses has been accelerated by heavy ion irradiation. For the same R value over the glass series, the initial Si dissolution rates in ternary borosilicate glasses change at similar rates after irradiation. This study concluded that the variation trend in the initial dissolution rate of Si and Na differs under ion radiation.

Vitrification is used for the immobilization and conditioning of high-level waste (HLW) arising from the reprocessing of spent nuclear fuel in the UK. Vitrification is also under consideration for the immobilization of certain intermediate-level wastes (ILW), where there may be advantages of volume reduction and removal of uncertainties in long-term waste behaviour, compared to encapsulation in a cement grout. This paper gives an overview of recent work into the chemical durability of UK vitrified wasteforms to inform the technical specification for the disposal facilities for these waste products and the treatment of their long-term behaviour in post-closure performance assessment. This has included: (1) measurements of the initial glass dissolution rates of a simulated HLW Magnox waste glass in a range of groundwater types representative of potential UK host geologies and in simulated high pH near-field porewaters relevant to potential disposal concepts, using Product Consistency Test type-B (PCT-B) at 40°C; and (2) durability testing of three simulant ILW glasses in a saturated calcium hydroxide buffered solution to simulate conditions in cement-based disposal vaults, using PCT-B tests at 50°C.The experimentally defined initial rate of HLW Magnox waste glass dissolution in a range of simulated groundwater compositions appears to be similar regardless of the ionic strength and major element composition of the solution. The release of caesium from HLW Magnox waste glass appears to be sensitive to solution composition. Caesium is selectively retained in the glass compared to other soluble components in the two low ionic strength solutions, but is released at similar rates to other soluble components in the three groundwaters and Ca(OH)2 solution. Whether this change in caesium retention is an ionic strength effect or is related to changes in the nature of the surface alteration layer formed on the glass, has yet to be established. For HLW Magnox waste glass, dissolution is accelerated at high pH in NaOH solution, however, the presence of calcium acts to mitigate the effects of high pH, at least initially. In Ca(OH)2 solution, calcium is found to react with all the glasses studied leading to the formation of calcium-containing alteration products. The initial dissolution behaviour in Ca(OH)2 solution varies with glass composition and in particular appears to be sensitive to the boron content.

Quantitative Structure Property Relationship (QSPR) analysis based on molecular dynamics (MD) simulations is a promising approach for establishing the composition‐property relationships of glasses and other materials with complex structures. A series of 20 borosilicate, aluminosilicate, and boroaluminosilicate glasses have been modeled using MD simulations with recently developed effective potentials. Short‐ and medium‐range structures of these glasses were analyzed and, based on the structural information, QSPR analysis of the initial dissolution rates ( r 0 ) of these glasses that were measured at 90°C and pH 9 by using various structural descriptors such as percentage of bridging oxygen species, network connectivity, and average ring size. The structural descriptor, F net , which contains both energetic information such as single bond strength and structural information such as cation coordination number and Q n distribution, was also used. It was found that while the overall network connectivity, average ring size and F net provide reasonable correlations with r 0 of studied glasses, F net gives the best correlation among the descriptors. For glasses that show incongruent dissolution, it was found that modification of glass compositions to account for preferential release of modifier cations is necessary to achieve best correlations. The findings were discussed with results of recent studies on evaluating the compositional dependence of glass dissolution behavior using the topological‐constraints‐based models.