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

Molecular Driving Forces in the Self-Association of Silaffin Peptide R5 from MD Simulations

The 19-residue silaffin-R5 peptide has been widely studied for its ability to precipitate uniform SiO2 particles through mild temperature and pH pathways, in the absence of any organic solvents. There is consensus that post-translational modification (PTM) of side chains has a large impact on the biomineralization process. Thus, it is imperative to understand the precise mechanisms that dictate the formation of SiO2 from R5 peptide, including the effects of PTM on peptide aggregation and peptide-surface adsorption. In this work, we use molecular dynamics (MD) simulations to study the aggregation of R5 dimer with multiple PTMs, with the presence of different ions in solution. Since this system has strong interactions with deep metastable states, we use parallel bias metadynamics with partitioned families to efficiently sample the different states of the system. We find that peptide aggregation is a prerequisite for biomineralization. We observe that the electrostatic interactions are essential in the R5 dimer aggregation; for wild type R5 that only has positively charged residues, phosphate ions HPO4 2- in the solution form a bridge between two peptides and are essential for peptide aggregation.

High-Throughput Computational Screening of Solid-Binding Peptides

Inspired by biomineralization, a naturally occurring, protein-facilitated process, solid-binding peptides (SBPs) have gained much attention for their potential to fabricate various shaped nanocrystals and hierarchical nanostructures. The advantage of SBPs over other traditionally used synthetic polymers or short ligands is their tunable interaction with the solid material surface via carefully programmed sequence and being solution-dependent simultaneously. However, designing a sequence with targeted binding affinity or selectivity often involves intensive processes such as phage display, and only a limited number of sequences can be identified. Other computational efforts have also been introduced, but the validation process remains prohibitively expensive once a suitable sequence has been identified. In this paper, we present a new model to rapidly estimate the binding free energy of any given sequence to a solid surface. We show how the overall binding of a polypeptide can be estimated from the free energy contribution of each residue based on the statistics of the thermodynamically stable structure ensemble. We validated our model using five silica-binding peptides of different binding affinities and lengths and showed that the model is accurate and robust across a wider range of chemistries and binding strengths. The computational cost of this method can be as low as 3% of the commonly used enhanced sampling scheme for similar studies and has a great potential to be used in high-throughput algorithms to obtain larger training data sets for machine learning SBP screening.

Nucleation Domains in Biomineralization: Biomolecular Sequence and Conformational Features

Biomolecules play a vital role in the regulation of biomineralization. However, the characteristics of practical nucleation domains are still sketchy. Herein, the effects of the representative biomolecular sequence and conformations on calcium phosphate (Ca-P) nucleation and mineralization are investigated. The results of computer simulations and experiments prove that the line in the arrangement of dual acidic/essential amino acids with a single interval (Bc (Basic) -N (Neutral) -Bc-N-Ac (Acidic)- NN-Ac-N) is most conducive to the nucleation. 2α-helix conformation can best induce Ca-P ion cluster formation and nucleation. "Ac- × × × -Bc" sequences with α-helix are found to be the features of efficient nucleation domains, in which process, molecular recognition plays a non-negligible role. It further indicates that the sequence determines the potential of nucleation/mineralization of biomolecules, and conformation determines the ability of that during functional execution. The findings will guide the synthesis of biomimetic mineralized materials with improved performance for bone repair.

Towards Polycaprolactone-Based Scaffolds for Alveolar Bone Tissue Engineering: A Biomimetic Approach in a 3D Printing Technique

The alveolar bone is a unique type of bone, and the goal of bone tissue engineering (BTE) is to develop methods to facilitate its regeneration. Currently, an emerging trend involves the fabrication of polycaprolactone (PCL)-based scaffolds using a three-dimensional (3D) printing technique to enhance an osteoconductive architecture. These scaffolds are further modified with hydroxyapatite (HA), type I collagen (CGI), or chitosan (CS) to impart high osteoinductive potential. In conjunction with cell therapy, these scaffolds may serve as an appealing alternative to bone autografts. This review discusses research gaps in the designing of 3D-printed PCL-based scaffolds from a biomimetic perspective. The article begins with a systematic analysis of biological mineralisation (biomineralisation) and ossification to optimise the scaffold's structural, mechanical, degradation, and surface properties. This scaffold-designing strategy lays the groundwork for developing a research pathway that spans fundamental principles such as molecular dynamics (MD) simulations and fabrication techniques. Ultimately, this paves the way for systematic in vitro and in vivo studies, leading to potential clinical applications.

Exploration of high dimensional free energy landscapes by a combination of temperature-accelerated sliced sampling and parallel biasing

Temperature-accelerated sliced sampling (TASS) is an enhanced sampling method for achieving accelerated and controlled exploration of high-dimensional free energy landscapes in molecular dynamics simulations. With the aid of umbrella bias potentials, the TASS method realizes a controlled exploration and divide-and-conquer strategy for computing high-dimensional free energy surfaces. In TASS, diffusion of the system in the collective variable (CV) space is enhanced with the help of metadynamics bias and elevated-temperature of the auxiliary degrees of freedom (DOF) that are coupled to the CVs. Usually, a low-dimensional metadynamics bias is applied in TASS. In order to further improve the performance of TASS, we propose here to use a highdimensional metadynamics bias, in the same form as in a parallel bias metadynamics scheme. Here, a modified reweighting scheme, in combination with artificial neural network is used for computing unbiased probability distribution of CVs and projections of high-dimensional free energy surfaces. We first validate the accuracy and efficiency of our method in computing the four-dimensional free energy landscape for alanine tripeptide in vacuo. Subsequently, we employ the approach to calculate the eight-dimensional free energy landscape of alanine pentapeptide in vacuo. Finally, the method is applied to a more realistic problem wherein we compute the broad four-dimensional free energy surface corresponding to the deacylation of a drug molecule which is covalently complexed with a β-lactamase enzyme. We demonstrate that using parallel bias in TASS improves the efficiency of exploration of high-dimensional free energy landscapes.

Developments and Ongoing Challenges for Analysis of Surface-Bound Proteins

Proteins at surfaces and interfaces play important roles in the function and performance of materials in applications ranging from diagnostic assays to biomedical devices. To improve the performance of these materials, detailed molecular structure (conformation and orientation) along with the identity and concentrations of the surface-bound proteins on those materials must be determined. This article describes radiolabeling, surface plasmon resonance, quartz crystal microbalance with dissipation, X-ray photoelectron spectroscopy, secondary ion mass spectrometry, sum frequency generation spectroscopy, and computational techniques along with the information each technique provides for characterizing protein films. A multitechnique approach using both experimental and computation methods is required for these investigations. Although it is now possible to gain much insight into the structure of surface-bound proteins, it is still not possible to obtain the same level of structural detail about proteins on surfaces as can be obtained about proteins in crystals and solutions, especially for large, complex proteins. However, recent results have shown it is possible to obtain detailed structural information (e.g., backbone and side chain orientation) about small peptides (5-20 amino sequences) on surfaces. Current studies are extending these investigations to small proteins such as protein G B1 (∼6 kDa). Approaches for furthering the capabilities for characterizing the molecular structure of surface-bound proteins are proposed.

Role of Water in CaCO3 Biomineralization

Biomineralization occurs in aqueous environments. Despite the ubiquity and relevance of CaCO₃ biomineralization, the role of water in the biomineralization process has remained elusive. Here, we demonstrate that water reorganization accompanies CaCO₃ biomineralization for sea urchin spine generation in a model system. Using surface-specific vibrational spectroscopy, we probe the water at the interface of the spine-associated protein during CaCO₃ mineralization. Our results show that, while the protein structure remains unchanged, the structure of interfacial water is perturbed differently in the presence of both Ca²⁺ and CO₃^2– compared to the addition of only Ca²⁺. This difference is attributed to the condensation of prenucleation mineral species. Our findings are consistent with a nonclassical mineralization pathway for sea urchin spine generation and highlight the importance of protein hydration in biomineralization.

Probing the binding modes and dynamics of histidine on fumed silica surfaces by solid-state NMR

Silica nanoparticles can be designed to exhibit a diverse range of morphologies (e.g. non-porous, mesoporous), physical properties (e.g. hydrophobic, hydrophilic) and a wide range of chemical and biomolecular surface functionalizations. In the present work, the adsorption complex of histidine (His) and fumed silica nanoparticles (FSN) is probed using thermal analysis (TGA/DTG) and a battery of solid-state (SS) NMR methods supported by DFT chemical shift calculations. Multinuclear (1H/13C/15N) one- and two-dimensional magic angle spinning (MAS) SSNMR experiments were applied to determine site-specific interactions between His and FSN surfaces as a function of adsorption solution concentration, pH and hydration state. By directly comparing SSNMR observables (linewidth, chemical shift and relaxation parameters) for His-FSN adsorption complexes to various crystalline, amorphous and aqueous His forms, the His structural and dynamic environment on FSN surfaces could be determined at an atomic level. The observed 13C and 15N MAS NMR chemical shifts, linewidths and relaxation parameters show that the His surface layer on FSN has a significant dependence on pH and hydration state. His is highly dynamic on FSN surfaces under acidic conditions (pH 4) as evidenced by sharp resonances with near isotropic chemical shifts regardless of hydration level indicating a non-specific binding arrangement while, a considerably more rigid His environment with defined protonation states is observed at near neutral pH with subtle variations between hydrated and anhydrous complexes. At near neutral pH, less charge repulsion occurs on the FSN surface and His is more tightly bound as evidenced by considerable line broadening likely due to chemical shift heterogeneity and a distribution in hydrogen-bonding strengths on the FSN surface. Multiple His sites exchange with a tightly bound water layer in hydrated samples while, direct interaction with the FSN surface and significant chemical shift perturbations for imidazole ring nitrogen sites and some carbon resonances are observed after drying. The SSNMR data was used to propose an interfacial molecular binding model between His and FSN surfaces under varying conditions setting the stage for future multinuclear, multidimensional SSNMR studies of His-containing peptides on silica nanoparticles and other nanomaterials of interest.

Data-Driven Molecular Dynamics: A Multifaceted Challenge

The big data concept is currently revolutionizing several fields of science including drug discovery and development. While opening up new perspectives for better drug design and related strategies, big data analysis strongly challenges our current ability to manage and exploit an extraordinarily large and possibly diverse amount of information. The recent renewal of machine learning (ML)-based algorithms is key in providing the proper framework for addressing this issue. In this respect, the impact on the exploitation of molecular dynamics (MD) simulations, which have recently reached mainstream status in computational drug discovery, can be remarkable. Here, we review the recent progress in the use of ML methods coupled to biomolecular simulations with potentially relevant implications for drug design. Specifically, we show how different ML-based strategies can be applied to the outcome of MD simulations for gaining knowledge and enhancing sampling. Finally, we discuss how intrinsic limitations of MD in accurately modeling biomolecular systems can be alleviated by including information coming from experimental data.

Studies of Dynamic Binding of Amino Acids to TiO2 Nanoparticle Surfaces by Solution NMR and Molecular Dynamics Simulations

Adsorption of biomolecules onto material surfaces involves a potentially complex mechanism where molecular species interact to varying degrees with a heterogeneous material surface. Surface adsorption studies by atomic force microscopy, sum frequency generation spectroscopy, and solid-state NMR detect the structures and interactions of biomolecular species that are bound to material surfaces, which, in the absence of a solid-liquid interface, do not exchange rapidly between surface-bound forms and free molecular species in bulk solution. Solution NMR has the potential to complement these techniques by detecting and studying transiently bound biomolecules at the liquid-solid interface. Herein, we show that dark-state exchange saturation transfer (DEST) NMR experiments on gel-stabilized TiO₂ nanoparticle (NP) samples detect several forms of biomolecular adsorption onto titanium(IV) oxide surfaces. Specifically, we use the DEST approach to study the interaction of amino acids arginine (Arg), lysine (Lys), leucine (Leu), alanine (Ala), and aspartic acid (Asp) with TiO₂ rutile NP surfaces. Whereas Leu, Ala, and Asp display only a single weakly interacting form in the presence of TiO₂ NPs, Arg and Lys displayed at least two distinct bound forms: a species that is surface bound and retains a degree of reorientational motion and a second more tightly bound form characterized by broadened DEST profiles upon the addition of TiO₂ NPs. Molecular dynamics simulations indicate different surface bound states for both Lys and Arg depending on the degree of TiO₂ surface hydroxylation but only a single bound state for Asp regardless of the degree of surface hydroxylation, in agreement with results obtained from the analysis of DEST profiles.

Structural Biology of Calcium Phosphate Nanoclusters Sequestered by Phosphoproteins

Biofluids that contain stable calcium phosphate nanoclusters sequestered by phosphopeptides make it possible for soft and hard tissues to co-exist in the same organism with relative ease. The stability diagram of a solution of nanocluster complexes shows how the minimum concentration of phosphopeptide needed for stability increases with pH. In the stable region, amorphous calcium phosphate cannot precipitate. Nevertheless, if the solution is brought into contact with hydroxyapatite, the crystalline phase will grow at the expense of the nanocluster complexes. The physico-chemical principles governing the formation, composition, size, structure, and stability of the complexes are described. Examples are given of complexes formed by casein, osteopontin, and recombinant phosphopeptides. Application of these principles and properties to blood serum, milk, urine, and resting saliva is described to show that under physiological conditions they are in the stable region of their stability diagram and so cannot cause soft tissue calcification. Stimulated saliva, however, is in the metastable region, consistent with its role in tooth remineralization. Destabilization of biofluids, with consequential ill-effects, can occur when there is a failure of homeostasis, such as an increase in pH without a balancing increase in the concentration of sequestering phosphopeptides.

Impact of Glutamate Carboxylation in the Adsorption of the α-1 Domain of Osteocalcin to Hydroxyapatite and Titania

One proposed mechanism of implant fouling is attributed to the nonspecific adsorption of non-collagenous bone matrix proteins (NCPs) onto a newly implanted interface. With the goal of capturing the fundamental mechanistic and thermodynamic forces that govern changes in these NCP recognition domains as a function of γ-carboxyglutamic acid (Gla) post-translational modification and surface chemistry, we probe the adsorption process of the most commonly occurring NCP, osteocalcin, onto a mineral and metal oxide surface. Here, we apply two enhanced sampling methods to independently probe the effects of post-translational modification and peptide structure on adsorption. First, well-tempered metadynamics was used to capture the binding of acetyl and N-methylamide capped glutamic acid and Gla single amino acids onto crystalline hydroxyapatite and titania model surfaces at physiological pH. Following this, parallel tempering metadynamics in the well-tempered ensemble (PTMetaD-WTE) was used to study adsorption of the α-1 domain of osteocalcin onto hydroxyapatite and titania. Simulations were performed for the α-1 domain of osteocalcin in both its fully decarboxylated (dOC) and fully carboxylated (OC) form. Our simulations find that increased charge density due to carboxylation results in increased interactions at the interface, and stronger adsorption of the single amino acids to both surfaces. Interestingly, the role of Gla in promoting compact and helical structure in the α-1 domain resulted in disparate binding modes at the two surfaces, which is attributed to differences in interfacial water behavior. Overall, this work provides a benchmark for understanding the mechanisms that drive adsorption of Gla-containing mineralizing proteins onto different surface chemistries.

Anisotropy in Stable Conformations of Hydroxylate Ions between the {001} and {110} Planes of TiO2 Rutile Crystals for Glycolate, Lactate, and 2-Hydroxybutyrate Ions Studied by Metadynamics Method

Control over TiO2 rutile crystal growth and morphology using additives is essential for the development of functional materials. Computer simulation studies on the thermodynamically stable conformations of additives at the surfaces of rutile crystals contribute to understanding the mechanisms underlying this control. In this study, a metadynamics method was combined with molecular dynamics simulations to investigate the thermodynamically stable conformations of glycolate, lactate, and 2-hydroxybutyrate ions at the {001} and {110} planes of rutile crystals. Two simple atom-atom distances were selected as collective variables for the metadynamics method. At the {001} plane, a conformation in which the COO- group was oriented toward the surface was found to be the most stable for the lactate and 2-hydroxybutyrate ions, whereas a conformation in which the COO- group was oriented toward water was the most stable for the glycolate ion. At the {110} plane, a conformation in which the COO- group was oriented toward the surface was the most stable for all three hydroxylate ions, and a second most stable conformation was also observed for the lactate ion at positions close to the {110} plane. For all three hydroxylate ions (α-hydroxycarboxylate ions), the stability of the most stable conformation was higher for the {110} plane than for the {001} plane. At both planes, the stability of the most stable conformation was highest for the 2-hydroxybutyrate ion and lowest for the glycolate ion. Supposing that all three hydroxylate ions serve to decrease the surface free energy at the rutile surface and that a more stable conformation at the rutile surface leads to a greater decrease in the surface free energy, the present results partially explain experimentally observed differences in the changes in growth rate and morphology of rutile crystals in the presence of glycolic, lactic, and 2-hydroxybutyric acids.

Computer simulations of the adsorption of an N-terminal peptide of statherin, SN15, and its mutants on hydroxyapatite surfaces

Salt-bridge adsorption of the SN15 peptide and its mutants on the HAP(001) surface.