Force Field Parameterization of Actinyl Molecular Cations Using the 12-6-4 Model

Fast and Robust Modeling of Lanthanide and Actinide Complexes, Biomolecules, and Molecular Crystals with the Extended GFN-FF Model

Lanthanides (Ln) and actinides (An) have recently become important tools in biomedical and materials science. However, the development of computational methods able to describe such elements in various environments has not kept up with the pace of the field. Addressing this challenge, this work introduces and showcases an extension of the GFN-FF to An alongside a reparameterization for Ln. This development fills a gap for fast computational methods that are out-of-the-box applicable to large f-element-containing systems with thousands of atoms. We discuss the reparameterization of the charge model and the covalent topology setup and showcase the model through various applications: Molecular dynamics simulations, optimization of Ln-containing biomolecules, and optimizations of several periodic structures. With the presented improvements, GFN-FF is a powerful method that routinely delivers robust and accurate geometries for large Ln/An systems with thousands of atoms.

Intelligent Design of Antithrombotic Peptide Targeting Collagen

Binding of blood components to collagen was proved to be a key step in thrombus formation. Intelligent Design of Protein Matcher (IDProMat), a neural network model, was then developed based on the principle of seq2seq to design an antithrombotic peptide targeting collagen. The encoding and decoding of peptide sequence data and the interaction patterns of peptide chains at the interface were studied, and then, IDProMat was applied to the design of peptides to cover collagen. The 99.3% decrease in seq2seq loss and 58.3% decrease in MLP loss demonstrated that IDProMat learned the interaction patterns between residues at the binding interface. An efficient peptide, LRWNSYY, was then designed using this model. Validations on its binding on collagen and its inhibition of platelet adhesion were obtained using docking, MD simulations, and experimental approaches.

Towards tailoring hydrophobic interaction with uranyl(VI) oxygen for C-H activation

Bovine serum albumin (BSA) has a uranyl(VI) binding hotspot where uranium is tightly bound by three carboxylates. Uranyl oxygen is "soaked" into the hydrophobic core of BSA. Isopropyl hydrogen of Val is trapped near UO22+ and upon photoexcitation, C-H bond cleavage is initiated. A unique hydrophobic contact with "yl"-oxygen, as observed here, can be used to induce C-H activation.

Modeling Lanthanide Ions in Solution: A Versatile Force Field in Aqueous and Organic Solvents

In this paper, we propose a new nonpolarizable force field for describing the Ln3+ (Ln = lanthanide) series based on a 12-6-4 Lennard-Jones potential. The development of the force field was performed in pure water by adjusting both the ion-oxygen distance and the hydration free energy. This force field accurately reproduces the Ln3+ hydration properties through the series, especially the coordination number that is hardly accessible using a nonpolarizable force field. Then, the validity and the transferability of the current force field were evaluated for two different systems containing Ln3+ in various solvents, namely, 0.1 mol L-1 La(NO3)3 salts in methanol and Eu(NO3)3 salts in solvent organic phases composed of DMDOHEMA molecules in n-heptane. The good agreement between our simulations and the data available in the literature confirms the accuracy of the force field for describing the lanthanide cations in both aqueous and nonaqueous media.

Nature and coordination geometry of geologically relevant aqueous Uranium(VI) complexes up to 400 ºC: A review and new data

The structure of the uranyl aqua ion (UO₂²⁺) and a number of its inorganic complexes (specifically, UO₂Cl⁺, UO₂Cl₂⁰, UO₂SO₄⁰, [Formula: see text] , [Formula: see text] and UO₂OH₄^2-) have been characterised using X-Ray absorption spectroscopy/extended X-Ray absorption fine structure (XAS/EXAFS) at temperatures ranging from 25 to 326 ºC. Results of ab initio molecular dynamics (MD) calculations are also reported for uranyl in chloride and sulfate-bearing fluids from 25 to 400 ºC and 600 bar to 20 kilobar (kb). These results are reported alongside a comprehensive review of prior structural characterisation work with particular focus given to EXAFS works to provide a consistent and up-to-date view of the structure of these complexes under conditions relevant to U mobility in ore-forming systems and around high-grade nuclear waste repositories. Regarding reported EXAFS results, average equatorial coordination was found to decrease in uranyl and its sulfate and chloride complexes as temperature rose - the extent of this decrease differed between species and solution compositions but typically resulted in an equatorial coordination number of ∼3-4 at temperatures above 200 ºC. The [Formula: see text] complex was observed at temperatures from 25 to 247 ºC and exhibited no major structural changes over this temperature range. UO₂OH₄^2- exhibited only minor structural changes over a temperature range from 88 to 326 ºC and was suggested to manifest fivefold coordination with four hydroxyl molecules and one water molecule around its equator. Average coordination values derived from fits of the reported EXAFS data were compared to average coordination values calculated using the experimentally derived thermodynamic data for chloride complexes reported by Dargent et al. (2013) and Migdisov et al. (2018b), and for sulfate complexes reported by Alcorn et al. (2019) and Kalintsev et al. (2019). Sulfate EXAFS data were well described by available thermodynamic data, and chloride EXAFS data were described well by the thermodynamic data of Migdisov et al. (2018b), but not by the data of Dargent et al. (2013). The ab initio molecular dynamics calculations confirmed the trends in equatorial coordination observed with EXAFS and were also able to provide an insight into the effect of pressure in equatorial water coordination - for a given temperature, higher pressures appear to lead to a greater number of equatorially bound waters counteracting the temperature effect.