Engineering a small diameter vascular graft with mechanical and biological properties comparable to living tissues remains challenging. Often, current devices lead to thrombosis and unsatisfactory long-term patency as a result of poor blood compatibility and a mismatch between the mechanical properties of the living tissue and the implanted biomaterial. Addressing all these requirements is essential to produce scaffolds able to survive throughout the life of the patient. For this purpose, we fabricated a novel three-layered vascular graft by combining electrospinning and braiding. Mirroring the structure of human blood vessels, the proposed device is composed of three layers: the intima, the media, and the adventitia. The intima and media layers were obtained by sequentially electrospinning silk fibroin (SF) and poly(L-lactide-co-ε-caprolactone), with ratios selected to match the mechanical properties of the native tissue. For the outer layer, the adventitia, SF yarns were braided on top of the electrospun tubes to create a structure able to withstand high pressures. Compliance, Young's modulus and deformability of the obtained scaffold were similar to that of human blood vessels. Additionally, cytocompatibility of the two layers, media and intima, was assessed in vitro by analysing cell metabolic activity and proliferation of endothelial cells and smooth muscle cells, respectively. Furthermore, heparin functionalization of the scaffolds led to improved anticoagulant properties upon incubation in whole blood. The obtained results indicate a potential application of the herewith designed three-layered construct as a vascular graft for small diameter blood vessel engineering.

A new force field has been derived for the aqueous calcium phosphate system that aims to reproduce the key thermodynamic properties of the system, including free energies of hydration of the ions and the solubility of the solid mineral phases. Interactions of three phosphate anions (PO, HPO, and HPO) with water were calibrated through comparison with the results obtained from ab initio molecular dynamics using both GGA and hybrid density functional theory with dispersion corrections. In the solid state, the force field has been evaluated by benchmarking against experiment and other existing models and is shown to reproduce the structural and mechanical properties well, despite the primary focus being on thermodynamics. To validate the force field, the thermodynamics of ion pairing for calcium phosphate species in water has been computed and shown to be in excellent agreement with experimental data.

Classical molecular dynamics simulations and free energy methods have been used to obtain a better understanding of the molecular processes occurring prior to the first nucleation event for calcium phosphate biominerals. The association constants for the formation of negatively charged complexes containing calcium and phosphate ions in aqueous solution have been computed, and these results suggest that the previously proposed calcium phosphate building unit, [Ca­(HPO4)3]4–, should only be present in small amounts under normal experimental conditions. However, the presence of an activation barrier for the removal of an HPO4 2– ion from this complex indicates that this species could be kinetically trapped. Aggregation pathways involving CaHPO4, [Ca­(HPO4)2]2–, and [Ca­(HPO4)3]4– complexes have been explored with the finding that dimerization is favorable up to a Ca/HPO4 ratio of 1:2.

Biological organisms display sophisticated control of nucleation and crystallization of minerals. In order to mimic living systems, deciphering the mechanisms by which organic molecules control the formation of mineral phases from solution is a key step. We have used computer simulations to investigate the effects of the amino acids arginine, aspartic acid, and glycine on species that form in solutions of calcium carbonate (CaCO3) at lower and higher levels of supersaturation. This provides net positive, negative, and neutral additives. In addition, we have prepared simulations containing hexapeptides of the amino acids to consider the effect of additive size on the solution species. We find that additives have limited impact on the formation of extended, liquid-like CaCO3 networks in supersaturated solutions. Additives control the amount of (bi)­carbonate in solution, but more importantly, they are able to stabilize these networks on the time scales of the simulations. This is achieved by coordinating the networks and assembled additive clusters in solutions. The association leads to subtle changes in the coordination of CaCO3 and reduced mobility of the cations. We find that the number of solute association sites and the size and topology of the additives are more important than their net charge. Our results help to understand why polymer additives are so effective at stabilizing dense liquid CaCO3 phases.

Amorphous calcium carbonate is often the first phase to precipitate from solution during the mineralization of calcium carbonate, before the formation of one of the crystalline polymorphs. In vivo, this phase is believed to be essential for the manufacture of minerals displaying nonequilibrium morphologies. The precipitation of this, usually transient, phase and its subsequent transformation into one of the crystalline polymorphs can be controlled by organic molecules. Here, we present a series of molecular dynamics simulations that explore the amorphous calcium carbonate–water interface, the attachment of amino acids onto both hydrous and anhydrous amorphous calcium carbonate, and their effect on the surface. The results show that surface ions have a different coordination number distribution from bulk ions and can diffuse up to two orders of magnitude faster than their bulk counterparts, suggesting that crystallization is much more likely to occur in this region. All the amino acids investigated bind to the amorphous calcium carbonate surfaces. However, acidic amino acids have a clear preference for the surface of amorphous CaCO3·H2O. The favored mode of interaction of the amino acids is through amine and/or guanidine moieties. The important ramifications of the results for our understanding of protein–mineral interactions are discussed.

Transdermal drug delivery is a key technology for administering drugs. However, most devices are "one-sizefits-all", even though drug diffusion through the skin varies significantly from person-to-person. For nextgeneration devices, personalization for optimal drug release would benefit from an augmented insight into the drug release and percutaneous uptake kinetics. Our objective was to quantify the changes in transdermal fentanyl uptake with regards to the patient's age and the anatomical location where the patch was placed. We also explored to which extent the drug flux from the patch could be altered by miniaturizing the contact surface area of the patch reservoir with the skin. To this end, we used validated mechanistic modeling of fentanyl diffusion, storage, and partitioning in the epidermis to quantify drug release from the patch and the uptake within the skin. A superior spatiotemporal resolution compared to experimental methods enabled in-silico identification of peak concentrations and fluxes, and the amount of stored drug and bioavailability. The patients' drug uptake showed a 36% difference between different anatomical locations after 72 h, but there was a strong interpatient variability. With aging, the drug uptake from the transdermal patch became slower and less potent. A 70-year-old patient received 26% less drug over the 72-h application period, compared to an 18-year-old patient. Additionally, a novel concept of using micron-sized drug reservoirs was explored in silico. These reservoirs induced a much higher local flux (µg cm -2 h -1 ) than conventional patches. Up to a 200-fold increase in the drug flux was obtained from these small reservoirs. This effect was mainly caused by transverse diffusion in the stratum corneum, which is not relevant for much larger conventional patches. These micronsized drug reservoirs open new ways to individualize reservoir design and thus transdermal therapy. Such computer-aided engineering tools also have great potential for in-silico design and precise control of drug delivery systems. Here, the validated mechanistic models can serve as a key building block for developing digital twins for transdermal drug delivery systems.

Transdermal drug delivery is a key technology for administering drugs. However, most devices are "one-size-fits-all", even though drug diffusion through the skin varies significantly from person-to-person. For next-generation devices, personalization for optimal drug release would benefit from an augmented insight into the drug release and percutaneous uptake kinetics. Our objective was to quantify the changes in transdermal fentanyl uptake with regards to the patient's age and the anatomical location where the patch was placed. We also explored to which extent the drug flux from the patch could be altered by miniaturizing the contact surface area of the patch reservoir with the skin. To this end, we used validated mechanistic modeling of fentanyl diffusion, storage, and partitioning in the epidermis to quantify drug release from the patch and the uptake within the skin. A superior spatiotemporal resolution compared to experimental methods enabled in-silico identification of peak concentrations and fluxes, and the amount of stored drug and bioavailability. The patients' drug uptake showed a 36% difference between different anatomical locations after 72 h, but there was a strong interpatient variability. With aging, the drug uptake from the transdermal patch became slower and less potent. A 70-year-old patient received 26% less drug over the 72-h application period, compared to an 18-year-old patient. Additionally, a novel concept of using micron-sized drug reservoirs was explored in silico. These reservoirs induced a much higher local flux (µg cm-2 h-1) than conventional patches. Up to a 200-fold increase in the drug flux was obtained from these small reservoirs. This effect was mainly caused by transverse diffusion in the stratum corneum, which is not relevant for much larger conventional patches. These micron-sized drug reservoirs open new ways to individualize reservoir design and thus transdermal therapy. Such computer-aided engineering tools also have great potential for in-silico design and precise control of drug delivery systems. Here, the validated mechanistic models can serve as a key building block for developing digital twins for transdermal drug delivery systems.

We report results from studies using four different protocols to prepare hydrated amorphous calcium carbonate, ranging from random initial structures to melting hydrated mineral structures.All protocols give good agreement with experimental X-ray structure factors. However, the thermodynamic properties, ion coordination environments, and distribution of water for the structures produced by the protocols show statistically significant variation depending on the protocols used. We discuss the diffusivity of water through the various structures and its relation to experiments. We show that one protocol (based on melting ikaite) gives a structure where the water is mobile, due to the presence of porosity in the amorphous structure. We conclude that our models of hydrated amorphous calcium carbonate do give a range of behaviour that resembles that observed experimentally, although the variation is less marked in the simulations than in experiments.3