Design of New Dispersants Using Machine Learning and Visual Analytics.

Artificial intelligence (AI) is an emerging technology that is revolutionizing the discovery of new materials. One key application of AI is virtual screening of chemical libraries, which enables the accelerated discovery of materials with desired properties. In this study, we developed computational models to predict the dispersancy efficiency of oil and lubricant additives, a critical property in their design that can be estimated through a quantity named blotter spot. We propose a comprehensive approach that combines machine learning techniques with visual analytics strategies in an interactive tool that supports domain experts' decision-making. We evaluated the proposed models quantitatively and illustrated their benefits through a case study. Specifically, we analyzed a series of virtual polyisobutylene succinimide (PIBSI) molecules derived from a known reference substrate. Our best-performing probabilistic model was Bayesian Additive Regression Trees (BART), which achieved a mean absolute error of 5.50±0.34 and a root mean square error of 7.56±0.47, as estimated through 5-fold cross-validation. To facilitate future research, we have made the dataset, including the potential dispersants used for modeling, publicly available. Our approach can help accelerate the discovery of new oil and lubricant additives, and our interactive tool can aid domain experts in making informed decisions based on blotter spot and other key properties.

Dissolving Hydroxyolite: A DNA Molecule into Its Hydroxyapatite Mold.

In spite of the clinical importance of hydroxyapatite (HAp), the mechanism that controls its dissolution in acidic environments remains unclear. Knowledge of such a process is highly desirable to provide better understanding of different pathologies, as for example osteoporosis, and of the HAp potential as vehicle for gene delivery to replace damaged DNA. In this work, the mechanism of dissolution in acid conditions of HAp nanoparticles encapsulating double-stranded DNA has been investigated at the atomistic level using computer simulations. For this purpose, four consecutive (multi-step) molecular dynamics simulations, involving different temperatures and proton transfer processes, have been carried out. Results are consistent with a polynuclear decalcification mechanism in which proton transfer processes, from the surface to the internal regions of the particle, play a crucial role. In addition, the DNA remains protected by the mineral mold and transferred proton from both temperature and chemicals. These results, which indicate that biomineralization imparts very effective protection to DNA, also have important implications in other biomedical fields, as for example in the design of artificial bones or in the fight against osteoporosis by promoting the fixation of Ca(2+) ions.

An experimental-computer modeling study of inorganic phosphates surface adsorption on hydroxyapatite particles.

The adsorption of orthophosphate, pyrophosphate, triphosphate and a trisphosphonate onto hydroxyapatite has been examined using experiments and quantum mechanical calculations. Adsorption studies with FTIR and X-ray photoelectron spectroscopies have been performed considering both crystalline hydroxyapatite (HAp) and amorphous calcium phosphate particles, which were specifically prepared and characterized for this purpose. Density functional theory (DFT) calculations have been carried out considering the (100) and (001) surfaces of HAp, which were represented using 1 × 2 × 2 and 3 × 3 × 1 slab models, respectively. The adsorption of phosphate onto the two crystallographic surfaces is very much favored from an energetic point of view, which is fully consistent with current interpretations of the HAp growing process. The structures calculated for the adsorption of pyrophosphate and triphosphate evidence that this process is easier for the latter than for the former. Thus, the adsorption of pyrophosphate is severely limited by the surface geometry while the flexibility of triphosphate allows transforming repulsive electrostatic interactions into molecular strain. On the other hand, calculations predict that the trisphosphonate only adsorbs onto the (001) surface of HAp. Theoretical predictions are fully consistent with experimental data. Thus, comparison of DFT results and spectroscopic data suggests that the experimental conditions used to prepare HAp particles promote the predominance of the (100) surface. Accordingly, experimental identification of the adsorption of trisphosphonate onto such crystalline particles is unclear while the adsorption of pyrophosphate and triphosphate is clearly observed.

Synergistic approach to elucidate the incorporation of magnesium ions into hydroxyapatite.

Although the content of Mg(2+) in hard tissues is very low (typically ≤1.5 wt %), its incorporation into synthetic hydroxyapatite (HAp) particles and its role in the mineral's properties are still subject of intensive debate. A combined experimental-computational approach is used to answer many of the open questions. Mg(2+) -enriched HAp particles are prepared using different synthetic approaches and considering different concentrations of Mg(2+) in the reaction medium. The composition, morphology and structure of the resulting particles are investigated using X-ray photoelectron spectroscopy, energy dispersive X-ray spectroscopy, scanning and transmission electron microscopies, FTIR, and wide-angle X-ray diffraction. After this scrutiny, the role of the Mg(2+) in the first nucleation stages, before HAp formation, is investigated using atomistic molecular dynamics simulations. Saturated solutions are simulated with and without the presence of DNA, which has been recently used as a soft template in the biomineralization process. This synergistic investigation provides a complete picture of how Mg(2+) ions affect the mineralization from the first stages onwards.

Mineralization of DNA into nanoparticles of hydroxyapatite.

Encapsulation of DNA into hydroxyapatite (HAp) has been investigated using a rational approach that involves computer simulation and experimental techniques. The temporal evolution of the radial distribution functions derived from atomistic molecular dynamics simulations of Ca(2+), PO4(3-) and OH(-)-containing aqueous solutions in the presence and absence of B-DNA has been used to conclude that the backbone of the double helix acts as a template for HAp growth. More specifically, results reveal the formation of calcium phosphate clusters at the first stages of the simulations, which subsequently re-organize to nucleate HAp. This effect is produced in the absence and, especially, presence, of DNA indicating that the biomolecules do not inhibit but even promote mineral growth. Furthermore, computer simulations suggest that the diffusion of the OH(-) anions through the inorganic solution is the limiting step for the nucleation of the biomineral. Nanocapsules and crystalline nanorods of HAp containing DNA molecules inside have been prepared by mixing solutions containing Ca(2+) and PO4(3-) ions with fish sperm DNA at high pH. The dimensions and morphology of such nanostructures have been examined by transmission electron microscopy, while the characterization of the biomineral has been focused on the identification of DNA inside HAp using infrared, X-ray photoelectron and UV-vis spectroscopies, as well as gel electrophoresis. The biominerals reported in this work are important for biomedical applications requiring the protection of DNA from aggressive environmental conditions.

DNA adsorbed on hydroxyapatite surfaces.

Hydroxyapatite (HAp) particles with very different surface charges and compositions (i.e. different Ca/P and CO3 2-/PO4 3- ratios) have been obtained by varying the experimental conditions used during the chemical precipitation process. The DNA adsorption capacity and protection imparted against the attack of nucleases of HAp particles have been proved to depend on the surface charge while the buffering capacity is affected by the chemical composition. On the basis of both the surface charge and the crystallinity, the predominant planes at the surfaces of HAp particles have been identified. Atomistic molecular dynamics simulations of surfaces constructed with these planes (i.e. (001) and the two terminations of (010)) with the adsorbed B-DNA double helix have been performed to get microscopic understanding of the influence of the mineral in the biomolecule structure and the interaction energies. The results indicate that the DNA secondary structure is perfectly preserved on the (001) surface, this stability being accompanied by an attractive binding energy. In contrast, the (010) surface with PO4 3-, OH- and Ca2+ ions in the termination induces significant local and global deformations in the double helix, repulsive OH-(HAp)PO4 3- (DNA) interactions provoking the desorption of the biomolecule. Finally, although the termination of the (010) surface with PO4 3- and Ca2+ ions also deforms the double helix, it forms very strong attractive interactions with the biomolecule. These binding characteristics are in excellent agreement with the DNA adsorption and protection abilities experimentally determined for the HAp samples. Finally, the surface charge has been found less decisive than the chemical composition in the efficacy of the transfection process.

Modeling biominerals formed by apatites and DNA.

Different aspects of biominerals formed by apatite and DNA have been investigated using computer modeling tools. Firstly, the structure and stability of biominerals in which DNA molecules are embedded into hydroxyapatite and fluoroapatite nanopores have been examined by combining different molecular mechanics methods. After this, the early processes in the nucleation of hydroxyapatite at a DNA template have been investigated using molecular dynamics simulations. Results indicate that duplexes of DNA adopting a B double helix can be encapsulated inside nanopores of hydroxyapatite without undergoing significant distortions in the inter-strand hydrogen bonds and the intra-strand stacking. This ability of hydroxyapatite is practically independent of the DNA sequence, which has been attributed to the stabilizing role of the interactions between the calcium atoms of the mineral and the phosphate groups of the biomolecule. In contrast, the fluorine atoms of fluoroapatite induce pronounced structural distortions in the double helix when embedded in a pore of the same dimensions, resulting in the loss of its most relevant characteristics. On the other hand, molecular dynamics simulations have allowed us to observe the formation of calcium phosphate clusters at the surface of the B-DNA template. Electrostatic interactions between the phosphate groups of DNA and Ca(2+) have been found to essential for the formation of stable ion complexes, which were the starting point of calcium phosphate clusters by incorporating PO3(4) from the solution.

Factors governing the conformational tendencies of C(α)-ethylated α-amino acids: chirality and side-chain size effects.

The intrinsic conformational properties of the N-acetyl-N'-methylamide derivatives of d-C(α)-ethylglycine (Abu), d-C(α)-methyl-C(α)-ethylglycine (Iva), and C(α,α)-diethylglycine (Deg) have been investigated using quantum mechanical calculations in the gas phase and in chloroform, dichloromethane and aqueous solutions. Although the large number of flexible dihedral angles results in many minimum energy conformations, only a few of them are energetically representative because of the repulsive interactions between the ethyl groups and the backbone atoms. The conformational restrictions imposed by such repulsions increase as follows: Abu < Iva < Deg. The most important characteristics of the investigated residues are as follows: (i) the conformation tendencies of Abu resemble those of standard residues with similar constitution (e.g., Val or Leu); (ii) the properties of Iva are relatively similar to those of C(α,α)-dimethylglycine, which is the simplest C(α,α)-dialkylated α-amino acid, even though the former is more restricted than the latter; (iii) the conformational flexibility of Deg is even lower than that of C(α,α)-dialkylated residues with bulkier side groups (e.g., C(α,α)-diphenylglycine and C(α,α)-dibenzylglycine).

Intermolecular interactions in electron transfer through stretched helical peptides.

The helical peptide Cys-Ala-Lys-(Glu-Ala-Ala-Ala-Lys)(2)-Ala-NH-(CH(2))(2)-SH has been organized forming a self-assembled monolayer on gold (0.602 peptides per nm(2)), its conductance behavior under stretching conditions being studied using scanning tunnelling microscopy and current sensing atomic force microscopy. The helical conformation of the peptide has been found to play a fundamental role in the conductance. Moreover, variation of the current upon molecular stretching indicates that peptides can be significantly elongated before the conductance drops to zero, the critical elongation being 1.22 ± 0.47 nm. Molecular dynamics simulations of a single peptide in the free state and of a variable number of peptides tethered to a gold surface (i.e. densities ranging from 0.026 to 1.295 peptides per nm(2)) have indicated that the helical conformation is intrinsically favored in solvated environments while in desolvated environments it is retained because of the fundamental role played by peptide-peptide intermolecular interactions. The structure obtained for the system with 24 tethered peptides, with a density of 0.634 peptides per nm(2) closest to the experimental one, is in excellent agreement with experimental observations. On the other hand, simulations in which a single molecule is submitted to different compression and stretching processes while the rest remain in the equilibrium have been used to mimic the variation of the tip-substrate distance in experimental measures. Results allowed us to identify the existence, and in some cases coexistence, of intermolecular and intramolecular ionic ladders, suggesting that peptide-mediated electron transfer occurs through the hopping mechanism. Finally, quantum mechanical calculations have been used to investigate the variation of the electronic structure upon compression and stretching deformations.

Conformational profile of a proline-arginine hybrid.

The intrinsic conformational preferences of a new nonproteinogenic amino acid have been explored by computational methods. This tailored molecule, named ((ß)Pro)Arg, is conceived as a replacement for arginine in bioactive peptides when the stabilization of folded turn-like conformations is required. The new residue features a proline skeleton that bears the guanidilated side chain of arginine at the C(ß) position of the five-membered pyrrolidine ring, in either a cis or a trans orientation with respect to the carboxylic acid. The conformational profiles of the N-acetyl-N'-methylamide derivatives of the cis and trans isomers of ((ß)Pro)Arg have been examined in the gas phase and in solution by B3LYP/6-31+G(d,p) calculations and molecular dynamics simulations. The main conformational features of both isomers represent a balance between geometric restrictions imposed by the five-membered pyrrolidine ring and the ability of the guanidilated side chain to interact with the backbone through hydrogen bonds. Thus, both cis- and trans-((ß)Pro)Arg exhibit a preference for the α(L) conformation as a consequence of the interactions established between the guanidinium moiety and the main-chain amide groups.