Elucidating the role of diverse mineralisation paradigms on bone biomechanics - a coarse-grained molecular dynamics investigation.

Bone as a hierarchical composite structure plays a myriad of roles in vertebrate skeletons including providing the structural stability of the body. Despite this critical role, the mechanical behaviour at the sub-micron levels of bone's hierarchy remains poorly understood. At this scale, bone is composed of Mineralised Collagen Fibrils (MCF) embedded within an extra-fibrillar matrix that consists of hydroxyapatite minerals and non-collagenous proteins. Recent experimental studies hint at the significance of the extra-fibrillar matrix in providing the bone with the stiffness and ductility needed to serve its structural roles. However, due to limited resolution of experimental tools, it is not clear how the arrangement of minerals, and in particular their relative distribution between the intra- and extra-fibrillar space contribute to bone's remarkable mechanical properties. In this study, a Coarse Grained Molecular Dynamics (CGMD) framework was developed to study the mechanical properties of MCFs embedded within an extra-fibrillar mineral matrix and the precise roles extra- and intra-fibrillar mineralisation on the load-deformation response was investigated. It was found that the presence of extra-fibrillar mineral resulted in the development of substantial residual stress in the system, by limiting MCF shortening that took place during intra-fibrillar mineralisation, resulting in substantial compressive residual stresses in the extra-fibrillar mineral phase. The simulation results also revealed the crucial role of extra-fibrillar mineralisation in determining the elastic response of the Extrafibrillar mineralised MCF (EFM-MCF) system up to the yield point, while the fibrillar collagen affected the post-yield behaviour. When physiological levels of mineralisation were considered, the mechanical response of the EFM-MCF systems was characterised by high ductility and toughness, with micro-cracks being distributed across the extra-fibrillar matrix, and MCFs effectively bridging these cracks leading to an excellent combination of strength and toughness. Together, these results provide novel insight into the deformation mechanisms of an EFM-MCF system and highlight that this universal building block, which forms the basis for lamellar bone, can provide an excellent balance of stiffness, strength and toughness, achieving mechanical properties that are far beyond the capabilities of the individual constituents acting alone.

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.

Ions Adsorbed at Amorphous Solid/Solution Interfaces Form Wigner Crystal-like Structures.

When a surface is immersed in a solution, it usually acquires a charge, which attracts counterions and repels co-ions to form an electrical double layer. The ions directly adsorbed to the surface are referred to as the Stern layer. The structure of the Stern layer normal to the interface was described decades ago, but the lateral organization within the Stern layer has received scant attention. This is because instrumental limitations have prevented visualization of the ion arrangements except for atypical, model, crystalline surfaces. Here, we use high-resolution amplitude modulated atomic force microscopy (AFM) to visualize in situ the lateral structure of Stern layer ions adsorbed to polycrystalline gold, and amorphous silica and gallium nitride (GaN). For all three substrates, when the density of ions in the layer exceeds a system-dependent threshold, correlation effects induce the formation of close packed structures akin to Wigner crystals. Depending on the surface and the ions, the Wigner crystal-like structure can be hexagonally close packed, cubic, or worm-like. The influence of the electrolyte concentration, species, and valence, as well as the surface type and charge, on the Stern layer structures is described. When the system parameters are changed to reduce the Stern layer ion surface excess below the threshold value, Wigner crystal-like structures do not form and the Stern layer is unstructured. For gold surfaces, molecular dynamics (MD) simulations reveal that when sufficient potential is applied to the surface, ion clusters form with dimensions similar to the Wigner crystal-like structures in the AFM images. The lateral Stern layer structures presented, and in particular the Wigner crystal-like structures, will influence diverse applications in chemistry, energy storage, environmental science, nanotechnology, biology, and medicine.

Synthesis of Chiral Triazine Frameworks for Enantiodiscrimination.

Manipulation of the structure of covalent organic frameworks at the molecular level is an efficient strategy to shift their biological, physicochemical, optical, and electrical properties in the desired windows. In this work, we report on a new method to construct chiral triazine frameworks using metal-driven polymerization for enantiodiscrimination. The nucleophilic substitution reaction between melamine and cyanuric chloride was performed in the presence of PdCl2, ZnCl2, and CuCl2 as chirality-directing agents. Palladium, with the ability of planar complex formation, was able to assemble monomers in two-dimensions and drive the reaction in two directions, leading to a two-dimensional triazine network with several micrometers lateral size. Nonplanar arrangements of monomers in the presence of ZnCl2 and CuCl2, however, resulted in calix and bouquet structures, respectively. While 2D and bouquet structures showed strong negative and positive bands in the CD spectra, respectively, their calix counterparts displayed long-range weak negative bands. In spite of the ability of both calix and bouquet networks to load l-histidine 35 and 50% more than d-histidine from pure enantiomers, respectively, only calix counterparts were able to take up this enantiomer (78%) from the racemic mixture. The two-dimensional polytriazine network did not show any specific interactions with pure enantiomers or their racemic mixtures.

Competition of opsonins and dysopsonins on the nanoparticle surface.

The biological behavior and fate of nanoparticles are dependent on their retention time in the blood circulation system. The protein corona components, especially opsonins, and dysopsonins, adsorbed on the nanoparticle surface determine their blood circulation time. The protein corona formation is a dynamic process that involves the competition between different proteins to be adsorbed on the nanoparticles. Therefore, studying how proteins compete and are oriented on the nanoparticle surface is essential. We hypothesized that the presence of opsonins (immunoglobulin (IgG)) might affect the adsorption of dysopsonins (human serum albumin (HSA)) and vice versa. Using the molecular dynamics simulations, we showed that the adsorption of HSA on the GO surface after the IgG adsorption is more probable than the opposite order of adsorption. It was also observed that the higher lateral diffusion of the HSA compared to the IgG helped the system reach a more stable configuration while the initial adsorption of the HSA limits the lateral diffusion of IgG. Therefore, replacing IgG adsorbed on the GO surface with HSA is plausible while the reverse process is less likely to occur. This study revealed that albumin might extend the blood circulation time of GO by replacing opsonins (IgG).

A coarse-grained molecular dynamics investigation of the role of mineral arrangement on the mechanical properties of mineralized collagen fibrils.

Mineralized collagen fibrils (MCFs) comprise collagen molecules and hydroxyapatite (HAp) crystals and are considered universal building blocks of bone tissue, across different bone types and species. In this study, we developed a coarse-grained molecular dynamics (CGMD) framework to investigate the role of mineral arrangement on the load-deformation behaviour of MCFs. Despite the common belief that the collagen molecules are responsible for flexibility and HAp minerals are responsible for stiffness, our results showed that the mineral phase was responsible for limiting collagen sliding in the large deformation regime, which helped the collagen molecules themselves undergo high tensile loading, providing a substantial contribution to the ultimate tensile strength of MCFs. This study also highlights different roles for the mineralized and non-mineralized protofibrils within the MCF, with the mineralized groups being primarily responsible for load carrying due to the presence of the mineral phase, while the non-mineralized groups are responsible for crack deflection. These results provide novel insight into the load-deformation behaviour of MCFs and highlight the intricate role that both collagen and mineral components have in dictating higher scale bone biomechanics.

Energy dissipation of osteopontin at a HAp mineral interface: Implications for bone biomechanics.

Osteopontin (OPN) is a one of the most abundant non-collagenous proteins in the bone's organic matrix. OPN is responsible for mediating bonding at mineral interfaces in the extrafibrillar space and recent evidence shows that it is a major contributor to bone's fracture resistance. While several experimental studies have identified an important role for calcium ions in mediating energy dissipation in OPN protein networks, the underlying molecular mechanisms remain largely unknown. In the current study, the role of calcium ions on energy dissipation at OPN interface with hydroxyapatite (HAp) as the main bone mineral was investigated. For the first time, the three-dimensional structure of OPN proteins were predicted, and it was found that calcium ions greatly influenced the final protein configuration and energy dissipation performance. Under small deformation, the compact cOPN structure, resulting from calcium ions presence, facilitated greater energy dissipation through sacrificial bond breaking and mechanisms mediated by the surface-bound calcium. At larger deformation, the compact structure also enabled cOPN to dissipate higher energy. Moreover, it was found that phosphorylation of OPN played an important role in energy dissipation. While previous studies have shown that OPN dissipated energy by forming aggregate networks, this study also showed that network formation is not necessary and that individual OPN proteins can dissipate large amounts of energy at HAp interfaces.

The structural role of osteocalcin in bone biomechanics and its alteration in Type-2 Diabetes.

This study presents an investigation into the role of Osteocalcin (OC) on bone biomechanics, with the results demonstrating that the protein's α-helix structures play a critical role in energy dissipation behavior in healthy conditions. In the first instance, α-helix structures have high affinity with the Hydroxyapatite (HAp) mineral surface and provide favorable conditions for adsorption of OC proteins onto the mineral surface. Using steered molecular dynamics simulation, several key energy dissipation mechanisms associated with α-helix structures were observed, which included stick-slip behavior, a sacrificial bond mechanism and a favorable binding feature provided by the Ca2+ motif on the OC protein. In the case of Type-2 Diabetes, this study demonstrated that possible glycation of the OC protein can occur through covalent crosslinking between Arginine and N-terminus regions, causing disruption of α-helices leading to a lower protein affinity to the HAp surface. Furthermore, the loss of α-helix structures allowed protein deformation to occur more easily during pulling and key energy dissipation mechanisms observed in the healthy configuration were no longer present. This study has significant implications for our understanding of bone biomechanics, revealing several novel mechanisms in OC's involvement in energy dissipation. Furthermore, these mechanisms can be disrupted following the onset of Type-2 Diabetes, implying that glycation of OC could have a substantial contribution to the increased bone fragility observed during this disease state.

Mechanistic Understanding of the Interactions between Nano-Objects with Different Surface Properties and α-Synuclein.

Aggregation of the natively unfolded protein α-synuclein (α-syn) is key to the development of Parkinson's disease (PD). Some nanoparticles (NPs) can inhibit this process and in turn be used for treatment of PD. Using simulation strategies, we show here that α-syn self-assembly is electrostatically driven. Dimerization by head-to-head monomer contact is triggered by dipole-dipole interactions and subsequently stabilized by van der Waals interactions and hydrogen bonds. Therefore, we hypothesized that charged nano-objects could interfere with this process and thus prevent α-syn fibrillation. In our simulations, positively and negatively charged graphene sheets or superparamagnetic iron oxide NPs first interacted with α-syn's N/C terminally charged residues and then with hydrophobic residues in the non-amyloid-ß component (61-95) region. In the experimental setup, we demonstrated that the charged nano-objects have the capacity not only to strongly inhibit α-syn fibrillation (both nucleation and elongation) but also to disaggregate the mature fibrils. Through the α-syn fibrillation process, the charged nano-objects induced the formation of off-pathway oligomers.

Disease-related metabolites affect protein-nanoparticle interactions.

Once in biological fluids, the surface of nanoparticles (NPs) is rapidly covered with a layer of biomolecules (i.e., the "protein corona") whose composition strongly determines their biological identity, regulates interactions with biological entities including cells and the immune system, and consequently directs the biological fate and pharmacokinetics of nanoparticles. We recently introduced the concept of a "personalized protein corona" which refers to the formation of different biological identities of the exact same type of NP after being exposed to extract plasmas from individuals who have various types of diseases. As different diseases have distinct metabolomic profiles and metabolites can interact with proteins, it is legitimate to hypothesize that metabolomic profiles in plasma may have the capacity to, at least partially, drive the formation of a personalized protein corona. To test this hypothesis, we employed a multi-scale approach composed of coarse-grained (CG) and all atom (AA) molecular dynamics (MD) simulations to probe the role of glucose and cholesterol (model metabolites in diabetes and hypercholesterolemia patients) in the interaction of fibrinogen protein and polystyrene NPs. Our results revealed that glucose and cholesterol had the capacity to induce substantial changes in the binding site of fibrinogen to the surface of NPs. More specifically, the simulation results demonstrated that increasing the metabolite amount could change the profiles of fibrinogen adsorption and replacement, what is known as the Vroman effect, on the NP surface. In addition, we also found out that metabolites can substantially determine the immune triggering potency of the fibrinogen-NP complex. Our proof-of-concept outcomes further emphasize the need for the development of patient-specific NPs in a disease type-specific manner for high yielding and safe clinical applications.