Adsorption of biomolecules and polymers on silicates, glasses, and oxides: mechanisms, predictions, and opportunities by molecular simulation

Effect of Pore Defects on Uniaxial Mechanical Properties of Bulk Hexagonal Hydroxyapatite: A Molecular Dynamics Study

Hydroxyapatite (HAP) is a calcium apatite bioceramic used in various naturally-derived and synthetic forms for bone repair and regeneration. While useful for the regrowth of osseus tissue, the poor load-bearing capacity of this material relative to other biomaterials is worsened by the propensity for pore formation during the synthetic processing of scaffolds, blocks, and granules. Here we use molecular dynamics (MD) simulations to improve the current understanding of the defect-altered uniaxial mechanical response in hexagonal HAP single crystals relative to defect-free structures. The inclusion of a central spherical pore within a repeated lattice was found to reduce both the failure stress and failure strain in uniaxial tension and compression, with up to a 30% reduction in maximum stress at the point of failure compared to a perfect crystalline structure observed when a 30 Å diameter pore was included. The Z axis ([0 0 0 1] crystalline direction) was found to be the least susceptible to pore defects in tension but the most sensitive to pore inclusion in compression. The deformation mechanisms are discussed to explain the observed mechanical responses, for which charge imbalances and geometric stress concentration factor effects caused by pore inclusion play a significant role.

Stable Binding Conformations of Polymaleic and Polyacrylic Acids at a Calcite Surface in the Presence of Countercations: A Metadynamics Study

Elucidating the stable binding conformations of additives at the surface of CaCO3 crystals is essential to biomineralization, scale inhibition, and materials technology. However, accomplishing this by experimental means is rather difficult. In this study, molecular dynamics simulations based on a metadynamics approach were conducted to elucidate the stable binding conformations of a deprotonated polymaleic acid (PMA) additive and two deprotonated poly(acrylic acid) (PAA) additives with different polymerization degrees in the presence of various countercations at a hydrated calcite (104) surface. The simulated free-energy surfaces suggested the existence of several slightly different stable binding conformations for each additive. The appearance of these distinct binding conformations is speculated to originate from different balances of interactions between the additive, the calcite surface, and the countercations. The binding conformations and binding stabilities at the calcite surface were affected by the countercations, with Ca2+ ions producing a more pronounced effect than Na+ ions. Furthermore, the simulation results suggested that the binding stability at the calcite surface was higher for the PMA additive than for the PAA additives, and the PAA additive with a polymerization degree of 10 displayed a binding stability that was similar to or lower than that of the PAA additive with a polymerization degree of 5. The present simulation method provides a new strategy for analyzing the binding conformations of complex additives at material surfaces, developing additives that stably bind to these surfaces, and designing additives to control crystal growth.

Study of nanoscale deformation mechanisms in bulk hexagonal hydroxyapatite under uniaxial loading using molecular dynamics

Hydroxyapatite (HAP) is a natural bioceramic which is currently used in scaffolds and coatings for the regrowth of osseous tissue but offers poor load-bearing capacity compared to other biomaterials. The deformation mechanisms responsible for the mechanical behavior of HAP are not well understood, although the advent of multiscale modeling offers the promise of improvements in many materials through computational materials science. This work utilizes molecular dynamics to study the nanoscale deformation mechanisms of HAP in uniaxial tension and compression. It was found that deformation mechanisms vary with loading direction in tension and compression leading to significant compression/tension asymmetry and crystal anisotropy. Bond orientation and geometry relative to the loading direction was found to be an indicator of whether a specific bond was involved in the deformation of HAP in each loading case. Tensile failure mechanisms were attributed to stretching and failure in loading case-specific ionic bond groups. The compressive failure mechanisms were attributed to coulombic repulsion in each case, although loading case-specific bond group rotation and displacement were found to affect specific failure modes. The elastic modulus was the highest for both tension and compression along the Z direction (i.e. normal to the basal plane), followed by Y and X.

An Overview of Molecular Modeling for Drug Discovery with Specific Illustrative Examples of Applications

In this paper we review the current status of high-performance computing applications in the general area of drug discovery. We provide an introduction to the methodologies applied at atomic and molecular scales, followed by three specific examples of implementation of these tools. The first example describes in silico modeling of the adsorption of small molecules to organic and inorganic surfaces, which may be applied to drug delivery issues. The second example involves DNA translocation through nanopores with major significance to DNA sequencing efforts. The final example offers an overview of computer-aided drug design, with some illustrative examples of its usefulness.

Solid-State NMR and MD Study of the Structure of the Statherin Mutant SNa15 on Mineral Surfaces

Elucidation of the structure and interactions of proteins at native mineral interfaces is key to understanding how biological systems regulate the formation of hard tissue structures. In addition, understanding how these same proteins interact with non-native mineral surfaces has important implications for the design of medical and dental implants, chromatographic supports, diagnostic tools, and a host of other applications. Here, we combine solid-state NMR spectroscopy, isotherm measurements, and molecular dynamics simulations to study how SNa15, a peptide derived from the hydroxyapatite (HAP) recognition domain of the biomineralization protein statherin, interacts with HAP, silica (SiO₂), and titania (TiO₂) mineral surfaces. Adsorption isotherms are used to characterize the binding affinity of SNa15 to HAP, SiO₂, and TiO₂. We also apply 1D ¹³C CP MAS, 1D ¹⁵N CP MAS, and 2D ¹³C-¹³C DARR experiments to SNa15 samples with uniformly ¹³C- and ¹⁵N-enriched residues to determine backbone and side-chain chemical shifts. Different computational tools, namely TALOS-N and molecular dynamics simulations, are used to deduce secondary structure from backbone and side-chain chemical shift data. Our results show that SNa15 adopts an α-helical conformation when adsorbed to HAP and TiO₂, but the helix largely unravels upon adsorption to SiO₂. Interactions with HAP are mediated in general by acidic and some basic amino acids, although the specific amino acids involved in direct surface interaction vary with surface. The integrated experimental and computational approach used in this study is able to provide high-resolution insights into adsorption of proteins on interfaces.

Adsorption Mechanisms of Nucleobases on the Hydrated Au(111) Surface

The solution environment is of fundamental importance in the adsorption of molecules on surfaces, a process that is strongly affected by the capability of the adsorbate to disrupt the hydration layer above the surface. Here we disclose how the presence of interface water influences the adsorption mechanism of DNA nucleobases on a gold surface. By means of metadynamics simulations, we describe the distinctive features of a complex free-energy landscape for each base, which manifests activation barriers for the adsorption process. We characterize the different pathways that allow each nucleobase to overcome the barriers and be adsorbed on the surface, discussing how they influence the kinetics of adsorption of single-stranded DNA oligomers with homogeneous sequences. Our findings offer a rationale as to why experimental data on the adsorption of single-stranded homo-oligonucleotides do not straightforwardly follow the thermodynamics affinity rank.

Characterization of Mechanical Degradation in Perfluoropolyether Film for Its Application to Antifingerprint Coatings

Enhancing the mechanical durability of antifingerprint films is critical for its industrial application on touch-screen devices to withstand friction damage from repeated rubbing in daily usage. Using reactive molecular dynamics simulations, we herein implement adhesion, mechanical, and deposition tests to investigate two durability-determining factors: intrachain and interchain strength, which affect the structural stability of the antifingerprint film (perfluoropolyether) on silica. From the intrachain perspective, it is found that the Si-C bond in the polymer chain is the weakest, and therefore prone to dissociation and potentially forming a C-O bond. This behavior is demonstrated consistently, regardless of the cross-linking density between polymer chains. For the interchain interaction, increasing the chain length enhances the mechanical properties of the film. Furthermore, the chain deposition test, mimicking the experimental coating process, demonstrates that placing shorter chains first to the surface of silica and then depositing longer chains is an ideal way to improve the interchain interaction in the film structure. The current study reveals a clear pathway to optimize the configuration of the polymer chain as well as its film structure to prolong the product life of the coated antifingerprint film.

A computational method for selecting short peptide sequences for inorganic material binding

Discovering or designing biofunctionalized materials with improved quality highly depends on the ability to manipulate and control the peptide-inorganic interaction. Various peptides can be used as assemblers, synthesizers, and linkers in the material syntheses. In another context, specific and selective material-binding peptides can be used as recognition blocks in mining applications. In this study, we propose a new in silico method to select short 4-mer peptides with high affinity and selectivity for a given target material. This method is illustrated with the calcite (104) surface as an example, which has been experimentally validated. A calcite binding peptide can play an important role in our understanding of biomineralization. A practical aspect of calcite is a need for it to be selectively depressed in mining sites.

Complexation of carboxylate on smectite surfaces

We report a first principles molecular dynamics (FPMD) study of carboxylate complexation on clay surfaces. By taking acetate as a model carboxylate, we investigate its inner-sphere complexes adsorbed on clay edges (including (010) and (110) surfaces) and in interlayer space. Simulations show that acetate forms stable monodentate complexes on edge surfaces and a bidentate complex with Ca²⁺ in the interlayer region. The free energy calculations indicate that the complexation on edge surfaces is slightly more stable than in interlayer space. By integrating pK_as and desorption free energies of Al coordinated water calculated previously (X. Liu, X. Lu, E. J. Meijer, R. Wang and H. Zhou, Geochim. Cosmochim. Acta, 2012, 81, 56-68; X. Liu, J. Cheng, M. Sprik, X. Lu and R. Wang, Geochim. Cosmochim. Acta, 2014, 140, 410-417), the pH dependence of acetate complexation has been revealed. It shows that acetate forms inner-sphere complexes on (110) in a very limited mildly acidic pH range while it can complex on (010) in the whole common pH range. The results presented in this study form a physical basis for understanding the geochemical processes involving clay-organics interactions.

Atomistic Modelling of Confined Polypropylene Chains between Ferric Oxide Substrates at Melt Temperature

The interactions and conformational characteristics of confined molten polypropylene (PP) chains between ferric oxide (Fe₂O₃) substrates were investigated by molecular dynamics (MD) simulations. A comparative analysis of the adsorbed amount shows strong adsorption of the chains on the high-energy surface of Fe₂O₃. Local structures formed in the polymer film were studied utilizing density profiles, orientation of bonds, and end-to-end distance of chains. At interfacial regions, the backbone carbon-carbon bonds of the chains preferably orient in the direction parallel to the surface while the carbon-carbon bonds with the side groups show a slight tendency to orient normal to the surface. Based on the conformation tensor data, the chains are compressed in the normal direction to the substrates in the interfacial regions while they tend to flatten in parallel planes with respect to the surfaces. The orientation of the bonds as well as the overall flattening of the chains in planes parallel to the solid surfaces are almost identical to that of the unconfined PP chains. Also, the local pressure tensor is anisotropic closer to the solid surfaces of Fe₂O₃ indicating the influence of the confinement on the buildup imbalance of normal and tangential pressures.

Integrated Modeling and Experimental Approaches to Control Silica Modification of Design Silk-Based Biomaterials

Mineralized polymeric biomaterials provide useful options toward mechanically robust systems for some tissue repairs. Silks as a mechanically robust protein-based material provide a starting point for biomaterial options, particularly when combined with silica toward organic-inorganic hybrid systems. To further understand the interplay between silk proteins and silica related to material properties, we systematically study the role of three key domains in bioengineered spider silk fusion proteins with respect to β-sheet formation and mineralization: (i) a core silk domain for materials assembly, (ii) a histidine tag for purification, and (iii) a silicification domain for mineralization. Computational simulations are used to identify the effect of each domain on the protein folding and accessibility of positively charged amino acids for silicification and predictions are then compared with experimental data. The results show that the addition of the silica and histidine domains reduces β-sheet structure in the materials, and increases solvent-accessible surface area to the positive charged amino acids, leading to higher levels of silica precipitation. Moreover, the simulations show that the location of the charged biomineralization domain has small effect on the protein folding and consequently surface exposure of charged amino acids. Those surfaces display correlation with the amount of silicification in experiments. The results demonstrate that the exposure of the positively charged amino acids impacts protein function related to mineralization. In addition, processing parameters (solvating agent, the method of β-sheet induction and temperature) affect protein secondary structure, folding and function. This integrated modeling and experimental approach provides insight into sequence-structure-function relationships for control of mineralized protein biomaterial structures.

Hydration in silica based mesoporous materials: a DFT model

In this work, calculable and realistic DFT models of MCM-41 material that follow temperature dependence of silanol density were developed. They can be easily applied in further studies of adsorption or as a support for catalysts.