Nature of peptide wrapping onto metal nanoparticle catalysts and driving forces for size control

Osteocalcin: Promoter or Inhibitor of Hydroxyapatite Growth?

Osteocalcin is the most abundant noncollagenous bone protein and the functions in bone remineralization as well as in inhibition of bone growth have remained unclear. In this contribution, we explain the dual role of osteocalcin in the nucleation of new calcium phosphate during bone remodeling and in the inhibition of hydroxyapatite crystal growth at the molecular scale. The mechanism was derived using pH-resolved all-atom models for the protein, phosphate species, and hydroxyapatite, along with molecular dynamics simulations and experimental and clinical observations. Osteocalcin binds to (hkl) hydroxyapatite surfaces through multiple residues, identified in this work, and the fingerprint of binding residues varies as a function of the (hkl) crystal facet and pH value. On balance, the affinity of osteocalcin to hydroxyapatite slows down crystal growth. The unique tricalcium γ-carboxylglutamic acid (Gla) domain hereby rarely adsorbs to hydroxyapatite surfaces and faces instead toward the solution. The Gla domain enables prenucleation of calcium phosphate for new bone formation at a slightly acidic pH of 5. The growth of prenucleation clusters of calcium phosphate continues upon increase in pH value from 5 to 7 and is much less favorable, or not observed, on the native osteocalcin structure at and above neutral pH values of 7. The results provide mechanistic insight into the early stages of bone remodeling from the molecular scale, help inform mutations of osteocalcin to modify binding to apatites, support drug design, and guide toward potential cures for osteoporosis and hyperosteogeny.

Resolving Atomic-Scale Structure and Chemical Coordination in High-Entropy Alloy Electrocatalysts for Structure-Function Relationship Elucidation

The recent breakthrough in confining five or more atomic species in nanocatalysts, referred to as high-entropy alloy nanocatalysts (HEAs), has revealed the possibilities of multielemental interactions that can surpass the limitations of binary and ternary electrocatalysts. The wide range of potential surface configurations in HEAs, however, presents a significant challenge in resolving active structural motifs, preventing the establishment of structure-function relationships for rational catalyst design and optimization. We present a methodology for creating sub-5 nm HEAs using an aqueous-based peptide-directed route. Using a combination of pair distribution function and X-ray absorption spectroscopy, HEA structure models are constructed from reverse Monte Carlo modeling of experimental data sets and showcase a clear peptide-induced influence on atomic-structure and chemical miscibility. Coordination analysis of our structure models facilitated the construction of structure-function correlations applied to electrochemical methanol oxidation reactions, revealing the complex interplay between multiple metals that leads to improved catalytic properties. Our results showcase a viable strategy for elucidating structure-function relationships in HEAs, prospectively providing a pathway for future materials design.

Accurate Force Fields for Atomistic Simulations of Oxides, Hydroxides, and Organic Hybrid Materials up to the Micrometer Scale

The simulation of metals, oxides, and hydroxides can accelerate the design of therapeutics, alloys, catalysts, cement-based materials, ceramics, bioinspired composites, and glasses. Here we introduce the INTERFACE force field (IFF) and surface models for α-Al2O3, α-Cr2O3, α-Fe2O3, NiO, CaO, MgO, β-Ca(OH)2, β-Mg(OH)2, and β-Ni(OH)2. The force field parameters are nonbonded, including atomic charges for Coulomb interactions, Lennard-Jones (LJ) potentials for van der Waals interactions with 12-6 and 9-6 options, and harmonic bond stretching for hydroxide ions. The models outperform DFT calculations and earlier atomistic models (Pedone, ReaxFF, UFF, CLAYFF) up to 2 orders of magnitude in reliability, compatibility, and interpretability due to a quantitative representation of chemical bonding consistent with other compounds across the periodic table and curated experimental data for validation. The IFF models exhibit average deviations of 0.2% in lattice parameters, <10% in surface energies (to the extent known), and 6% in bulk moduli relative to experiments. The parameters and models can be used with existing parameters for solvents, inorganic compounds, organic compounds, biomolecules, and polymers in IFF, CHARMM, CVFF, AMBER, OPLS-AA, PCFF, and COMPASS, to simulate bulk oxides, hydroxides, electrolyte interfaces, and multiphase, biological, and organic hybrid materials at length scales from atoms to micrometers. The nonbonded character of the models also enables the analysis of mixed oxides, glasses, and certain chemical reactions, and well-performing nonbonded models for silica phases, SiO2, are introduced. Automated model building is available in the CHARMM-GUI Nanomaterial Modeler. We illustrate applications of the models to predict the structure of mixed oxides, and energy barriers of ion migration, as well as binding energies of water and organic molecules in outstanding agreement with experimental data and calculations at the CCSD(T) level. Examples of model building for hydrated, pH-sensitive oxide surfaces to simulate solid-electrolyte interfaces are discussed.

Developing a rational approach to designing recombinant proteins for peptide-directed nanoparticle synthesis

The controlled formation of nanoparticles with optimum characteristics and functional aspects has proven successful via peptide-mediated nanoparticle synthesis. However, the effects of the peptide sequence and binding motif on surface features and physicochemical properties of nanoparticles are not well-understood. In this study, we investigate in a comparative manner how a specific peptide known as Pd4 and its two known variants may form nanoparticles both in an isolated state and when attached to a green fluorescent protein (GFPuv). More importantly, we introduce a novel computational approach to predict the trend of the size and activity of the peptide-directed nanoparticles by estimating the binding affinity of the peptide to a single ion. We used molecular dynamics (MD) simulations to explore the differential behavior of the isolated and GFP-fused peptides and their mutants. Our computed palladium (Pd) binding free energies match the typical nanoparticle sizes reported from transmission electron microscope pictures. Stille coupling and Suzuki-Miyaura reaction turnover frequencies (TOFs) also correspond with computationally predicted Pd binding affinities. The results show that while using Pd4 and its two known variants (A6 and A11) in isolation produces nanoparticles of varying sizes, fusing these peptides to the GFPuv protein produces nanoparticles of similar sizes and activity. In other words, GFPuv reduces the sensitivity of the nanoparticles to the peptide sequence. This study provides a computational framework for designing free and protein-attached peptides that helps in the synthesis of nanoparticles with well-regulated properties.

Adsorption and Diffusion of Oxygen on Pure and Partially Oxidized Metal Surfaces in Ultrahigh Resolution

The interaction of gas molecules with metal and oxide surfaces plays a critical role in corrosion, catalysis, sensing, and heterogeneous materials. However, insights into the dynamics of O₂ from picoseconds to microseconds have remained unavailable to date. We obtained 3D potential energy surfaces for adsorption of O₂ on 11 common pristine and partially oxidized (hkl) surfaces of Ni and Al in picometer resolution and high accuracy of 0.1 kcal/mol, identified binding sites, and surface mobility from 25 to 300 °C. We explain relative oxidation rates and parameters for oxide growth. We employed over 150 000 molecular mechanics and molecular dynamics simulations with the interface force field (IFF) using structural data from X-ray diffraction (XRD) and low-energy electron diffraction (LEED). The methods reach 10 to 50 times higher accuracy than possible before and are suited to analyze gas interactions with metals up to the micrometer scale including defects and irregular nanostructures.

Molecular Driving Force for Facet Selectivity of Sequence-Defined Amphiphilic Peptoids at Au-Water Interfaces

Shape-controlled colloidal nanocrystal syntheses often require facet-selective solution-phase chemical additives to regulate surface free energy, atom addition/migration fluxes, or particle attachment rates. Because of their highly tunable properties and robustness to a wide range of experimental conditions, peptoids represent a very promising class of next-generation functional additives for control over nanocrystal growth. However, understanding the origin of facet selectivity at the molecular level is critical to generalizing their design. Herein we employ molecular dynamics simulations and biased sampling methods and report stronger selectivity to Au(111) than to Au(100) for Nce3Ncp6, a peptoid that has been shown to assist the formation of 5-fold twinned Au nanostars. We find that facet selectivity is achieved through synergistic effects of both peptoid-surface and solvent-surface interactions. Moreover, the amphiphilic nature of Nce3Ncp6 together with the order of peptoid-peptoid and peptoid-surface binding energies, that is, peptoid-Au(100) < peptoid-peptoid < peptoid-Au(111), further amplifies its distinct collective behavior on different Au surfaces. Our studies provide a fundamental understanding of the molecular origin of facet-selective adsorption and highlight the possibility of future designs and uses of sequence-defined peptoids for predictive syntheses of nanocrystals with designed shapes and properties.

Interpretable molecular models for molybdenum disulfide and insight into selective peptide recognition

Molybdenum disulfide (MoS2) is a layered material with outstanding electrical and optical properties. Numerous studies evaluate the performance in sensors, catalysts, batteries, and composites that can benefit from guidance by simulations in all-atom resolution. However, molecular simulations remain difficult due to lack of reliable models. We introduce an interpretable force field for MoS2 with record performance that reproduces structural, interfacial, and mechanical properties in 0.1% to 5% agreement with experiments. The model overcomes structural instability, deviations in interfacial and mechanical properties by several 100%, and empirical fitting protocols in earlier models. It is compatible with several force fields for molecular dynamics simulation, including the interface force field (IFF), CVFF, DREIDING, PCFF, COMPASS, CHARMM, AMBER, and OPLS-AA. The parameters capture polar covalent bonding, X-ray structure, cleavage energy, infrared spectra, bending stability, bulk modulus, Young's modulus, and contact angles with polar and nonpolar solvents. We utilized the models to uncover the binding mechanism of peptides to the MoS2 basal plane. The binding strength of several 7mer and 8mer peptides scales linearly with surface contact and replacement of surface-bound water molecules, and is tunable in a wide range from -86 to -6 kcal mol-1. The binding selectivity is multifactorial, including major contributions by van-der-Waals coordination and charge matching of certain side groups, orientation of hydrophilic side chains towards water, and conformation flexibility. We explain the relative attraction and role of the 20 amino acids using computational and experimental data. The force field can be used to screen and interpret the assembly of MoS2-based nanomaterials and electrolyte interfaces up to a billion atoms with high accuracy, including multiscale simulations from the quantum scale to the microscale.

Chiral Restructuring of Peptide Enantiomers on Gold Nanomaterials

The use of biomolecules has been invaluable at generating and controlling optical chirality in nanomaterials; however, the structure and properties of the chiral biotemplate are not well understood due to the complexity of peptide-nanoparticle interactions. In this study, we show that the complex interactions between d-peptides and gold nanomaterials led to a chiral restructuring of peptides as demonstrated by circular dichroism and proteolytic cleavage of d-peptides via gold-mediated inversion of peptide chirality. The gold nanoparticles synthesized using d-peptide produce a highly ordered atomic surface and restructured peptide bonds for enzyme cleavage. Differences in gold nanoparticle catalyzed reduction of 4-nitrophenol were observed on the basis of the chiral peptide used in nanoparticle synthesis. Notably, the proteolytic cleavage of d-peptides on gold provides an opportunity for designing nanoparticle based therapeutics to treat peptide venoms, access new chemistries, or modulate the catalytic activity of nanomaterials.

Assembly of Gold Nanoparticles into Chiral Superstructures Driven by Circularly Polarized Light

Photon-to-matter chirality transfer offers both simplicity and universality to chiral synthesis, but its efficiency is typically low for organic compounds. Besides the fundamental importance of this process relevant for understanding the origin of homochirality on Earth, new pathways for imposing chiral bias during chemical process are essential for a variety of technologies from medicine to informatics. The strong optical activity of inorganic nanoparticles (NPs) affords photosynthetic routes to chiral superstructures using circularly polarized photons. Although plasmonic NPs are promising candidates for such synthetic routes due to the strong rotatory power of highly delocalized plasmonic states (Ma et al. Chem. Rev. 2017, 117 (12), 8041), realization of light-driven synthesis of chiral nanostructures has been more challenging for plasmonic NPs than for the semiconductor due to the short lifetime of the plasmonic states. Here we show that illumination of gold salt solutions with circularly polarized light induces the formation of NPs and their subsequent assembly into chiral nanostructures 10-15 nm in diameter. Despite their seemingly irregular shape, the resulting nanocolloids showed circular dichroism (CD) spectra with opposite polarity after exposure to photons with left and right circular polarization. The sign and spectral position of the CD peaks of illuminated dispersions matched those calculated for nanostructures with complex geometry identified from electron tomography images. Quantification of the complex shapes of NP assemblies using chirality measures revealed a direct correlation with the experimental spectra. The light-driven assembly of chiral nanostructures originates from the asymmetric displacement of NPs in dynamic assemblies by plasmonic fields followed by particle-to-particle attachment. The ability of gold NPs to "lock" the chirality of the incident photons in assembled nanostructures can be used to create a variety of chiral nanomaterials with plasmonic resonances.

Dynamics of carbohydrate strands in water and interactions with clay minerals: influence of pH, surface chemistry, and electrolytes

Carbohydrate hydrogels are extensively used in pharmaceuticals and engineered biomaterials. Molecular conformations, assembly, and interactions of the carbohydrate strands with stabilizers such as clay minerals in aqueous solution are difficult to quantify in experiments and the hydrogel properties remain largely a result of trial-and-error studies. We analyzed the assembly of gellan gum in aqueous solution and interactions with dispersed clay minerals in all-atomic detail using molecular dynamics simulation, atomic force microscopy (AFM), and comparisons to earlier measurements. Gellan strands associate at low pH values of 2 and gradually disassemble to double strands with weak association of -0.4 kcal per mole carbohydrate ring as the pH values increases to 9. Ionization of the carbonic acid side groups in the backbone extends the chains and accelerates the conformational dynamics via rapidly changing intramolecular ion bridges. Gellan interactions with clay minerals depend on the strength of electric triple layers between clay, cations, and anionic polymer strands, as well as weaker hydrogen bonds along the edges, which are tunable as a function of the clay surface chemistry, local ionic strength, and pH values. Interaction energies range from -4 to +6 kcal per mol carbohydrate ring and were most favorable for electric triple layers with high charge mobility, which can be achieved for intermediate cation exchange capacity of the clay mineral and high pH values to increase ionization of the clay edges and of the polymer. The findings provide understanding and help control the dynamics and stabilization of carbohydrate hydrogels by clay minerals.

Atomic-Scale Structure and Stress Release Mechanism in Core-Shell Nanoparticles

Core-shell nanoparticles find applications in catalysts, sensors, and theranostics. The full internal 3D atomic structure, however, cannot be resolved by current imaging and diffraction techniques. We analyzed the atomic positions and stress-release mechanism in a cubic Au-Pd core-shell nanoparticle in approximately 1000 times higher resolution than current experimental techniques using large-scale molecular dynamics simulation to overcome these limitations. The core-shell nanocube of 73 nm size was modeled similarly to solution synthesis by random epitaxial deposition of a 4 nm thick shell of Pd atoms onto a Au core of 65 nm side length using reliable interatomic potentials. The internal structure reveals specific deformations and stress relaxation mechanisms that are caused by the +4.8% lattice mismatch of gold relative to palladium and differential confinement of extended particle facets, edges, and corners by one, two, or three Au-Pd interfaces, respectively. The three-dimensional lattice strain causes long-range, arc-like bending of atomic rows along the faces and edges of the particle, especially near the Au-Pd interface, a bulging deformation of the Pd shell, and stacking faults in the Pd shell at the corners of the particle. The strain pattern and mechanism of stress release were further characterized by profiles of the atomic layer spacing in the principal crystallographic directions. Accordingly, strain in the Pd shell is several times larger in the extended facets than near the edges and corners of the nanoparticle, which likely affects adsorption, optical, and electrochemical properties. The findings are consistent with available experimental data, including 3D reconstructions of the same cubic nanoparticle by coherent diffractive imaging (CDI) and may be verified by more powerful experimental techniques in the future. The stress release mechanisms are representative for cubic core-shell nanoparticles with fcc structure and can be explored for different shapes by the same methods.

Interactions between Metal Oxides and Biomolecules: from Fundamental Understanding to Applications

Metallo-oxide (MO)-based bioinorganic nanocomposites promise unique structures, physicochemical properties, and novel biochemical functionalities, and within the past decade, investment in research on materials such as ZnO, TiO₂, SiO₂, and GeO₂ has significantly increased. Besides traditional approaches, the synthesis, shaping, structural patterning, and postprocessing chemical functionalization of the materials surface is inspired by strategies which mimic processes in nature. Would such materials deliver new technologies? Answering this question requires the merging of historical knowledge and current research from different fields of science. Practically, we need an effective defragmentation of the research area. From our perspective, the superficial accounting of material properties, chemistry of the surfaces, and the behavior of biomolecules next to such surfaces is a problem. This is particularly of concern when we wish to bridge between technologies in vitro and biotechnologies in vivo. Further, besides the potential practical technological efficiency and advantages such materials might exhibit, we have to consider the wider long-term implications of material stability and toxicity. In this contribution, we present a critical review of recent advances in the chemistry and engineering of MO-based biocomposites, highlighting the role of interactions at the interface and the techniques by which these can be studied. At the end of the article, we outline the challenges which hamper progress in research and extrapolate to developing and promising directions including additive manufacturing and synthetic biology that could benefit from molecular level understanding of interactions occurring between inanimate (abiotic) and living (biotic) materials.

Insight into induced charges at metal surfaces and biointerfaces using a polarizable Lennard-Jones potential

Metallic nanostructures have become popular for applications in therapeutics, catalysts, imaging, and gene delivery. Molecular dynamics simulations are gaining influence to predict nanostructure assembly and performance; however, instantaneous polarization effects due to induced charges in the free electron gas are not routinely included. Here we present a simple, compatible, and accurate polarizable potential for gold that consists of a Lennard–Jones potential and a harmonically coupled core-shell charge pair for every metal atom. The model reproduces the classical image potential of adsorbed ions as well as surface, bulk, and aqueous interfacial properties in excellent agreement with experiment. Induced charges affect the adsorption of ions onto gold surfaces in the gas phase at a strength similar to chemical bonds while ions and charged peptides in solution are influenced at a strength similar to intermolecular bonds. The proposed model can be applied to complex gold interfaces, electrode processes, and extended to other metals.

Molecular dynamics models for predicting the behavior of metallic nanostructures typically do not take into account polarization effects in metals. Here, the authors introduce a polarizable Lennard–Jones potential that provides quantitative insight into the role of induced charges at metal surfaces and related complex material interfaces.

Computational approaches to cell-nanomaterial interactions: keeping balance between therapeutic efficiency and cytotoxicity

Owing to their unique properties, nanomaterials have been widely used in biomedicine since they have obvious inherent advantages over traditional ones. However, nanomaterials may also cause dysfunction in proteins, genes and cells, resulting in cytotoxic and genotoxic responses. Recently, more and more attention has been paid to these potential toxicities of nanomaterials, especially to the risks of nanomaterials to human health and safety. Therefore, when using nanomaterials for biomedical applications, it is of great importance to keep the balance between therapeutic efficiency and cytotoxicity (i.e., increase the therapeutic efficiency as well as decrease the potential toxicity). This requires a deeper understanding of the interactions between various types of nanomaterials and biological systems at the nano/bio interface. In this review, from the point of view of theoretical researchers, we will present the current status regarding the physical mechanism of cytotoxicity caused by nanomaterials, mainly based on recent simulation results. In addition, the strategies for minimizing the nanotoxicity naturally and artificially will also be discussed in detail. Furthermore, we should notice that toxicity is not always bad for clinical use since causing the death of specific cells is the main way of treating disease. Enhancing the targeting ability of nanomaterials to diseased cells and minimizing their side effects on normal cells will always be hugely challenging issues in nanomedicine. By combining the latest computational studies with some experimental verifications, we will provide special insights into recent advances regarding these problems, especially for the design of novel environment-responsive nanomaterials.