Effects of Cell-Adhesive Ligand Presentation on Pentapeptide Supramolecular Assembly and Gelation: Simulations and Experiments.

The extracellular matrix (ECM) is a complex, hierarchical material containing structural and bioactive components. This complexity makes decoupling the effects of biomechanical properties and cell-matrix interactions difficult, especially when studying cellular processes in a 3D environment. Matrix mechanics and cell adhesion are both known regulators of specific cellular processes such as stem cell proliferation and differentiation. However, more information is required about how such variables impact various neural lineages that could, upon transplantation, therapeutically improve neural function after a central nervous system injury or disease. Rapidly Assembling Pentapeptides for Injectable Delivery (RAPID) hydrogels are one biomaterial approach to meet these goals, consisting of a family of peptide sequences that assemble into physical hydrogels in physiological media. In this study, we studied our previously reported supramolecularly-assembling RAPID hydrogels functionalized with the ECM-derived cell-adhesive peptide ligands RGD, IKVAV, and YIGSR. Using molecular dynamics simulations and experimental rheology, we demonstrated that these integrin-binding ligands at physiological concentrations (3-12 mm) did not impact the assembly of the KYFIL peptide system. In simulations, molecular measures of assembly such as hydrogen bonding and pi-pi interactions appeared unaffected by cell-adhesion sequence or concentration. Visualizations of clustering and analysis of solvent-accessible surface area indicated that the integrin-binding domains remained exposed. KYFIL or AYFIL hydrogels containing 3 mm of integrin-binding domains resulted in mechanical properties consistent with their non-functionalized equivalents. This strategy of doping RAPID gels with cell-adhesion sequences allows for the precise tuning of peptide ligand concentration, independent of the rheological properties. The controllability of the RAPID hydrogel system provides an opportunity to investigate the effect of integrin-binding interactions on encapsulated neural cells to discern how hydrogel microenvironment impacts growth, maturation, or differentiation.

Effects of solvent conditions on the self-assembly of heterotrimeric collagen-like peptide (CLP) triple helices: a coarse-grained simulation study.

We perform coarse-grained (CG) molecular dynamics (MD) simulations to investigate the self-assembly of collagen-like peptide (CLP) triple helices into fibrillar structures and percolated networks as a function of solvent quality. The focus of this study is on CLP triple helices whose strands are different lengths (i.e., heterotrimers), leading to dangling 'sticky ends'. These 'sticky ends' are segments of the CLP strands that have unbonded hydrogen-bonding donor/acceptor sites that drive heterotrimeric CLP triple helices to physically associate with one another, leading to assembly into higher-order structures. We use a validated CG model for CLP in implicit solvent and capture varying solvent quality through changing strength of attraction between CG beads representing the amino acids in the CLP strands. Our CG MD simulations show that, at lower CLP concentrations, CLP heterotrimers assemble into fibrils and, at higher CLP concentrations, into percolated networks. At higher concentrations, decreasing solvent quality causes (i) the formation of heterogeneous network structures with a lower degree of branching at network junctions and (ii) increases in the diameter of network strands and pore sizes. We also observe a nonmonotonic effect of solvent quality on distances between network junctions due to the balance between heterotrimer end-end associations driven by hydrogen bonding and side-side associations driven by worsening solvent quality. Below the percolation threshold, we observe that decreasing solvent quality leads to the formation of fibrils composed of multiple aligned CLP triple helices, while the number of 'sticky ends' governs the spatial extent (radius of gyration) of the assembled fibrils.

Explicit solvent machine-learned coarse-grained model of sodium polystyrene sulfonate to capture polymer structure and dynamics.

Strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties.

Impact of collagen-like peptide (CLP) heterotrimeric triple helix design on helical thermal stability and hierarchical assembly: a coarse-grained molecular dynamics simulation study.

Collagen-like peptides (CLP) are multifunctional materials garnering a lot of recent interest from the biomaterials community due to their hierarchical assembly and tunable physicochemical properties. In this work, we present a computational study that links the design of CLP heterotrimers to the thermal stability of the triple helix and their self-assembly into fibrillar aggregates and percolated networks. Unlike homotrimeric helices, the CLP heterotrimeric triple helices in this study are made of CLP strands of different chain lengths that result in 'sticky' ends with available hydrogen bonding groups. These 'sticky' ends at one end or both ends of the CLP heterotrimer then facilitate inter-helix hydrogen bonding leading to self-assembly into fibrils (clusters) and percolated networks. We consider the cases of three sticky end lengths - two, four, and six repeat units - present entirely on one end or split between two ends of the CLP heterotrimer. We observe in CLP heterotrimer melting curves generated using coarse grained Langevin dynamics simulations at low CLP concentration that increasing sticky end length results in lower melting temperatures for both one and two sticky ended CLP designs. At higher CLP concentrations, we observe non-monotonic trends in cluster sizes with increasing sticky end length with one sticky end but not for two sticky ends with the same number of available hydrogen bonding groups as the one sticky end; this nonmonotonicity stems from the formation of turn structures stabilized by hydrogen bonds at the single, sticky end for sticky end lengths greater than four repeat units. With increasing CLP concentration, heterotrimers also form percolated networks with increasing sticky end length with a minimum sticky end length of four repeat units required to observe percolation. Overall, this work informs the design of thermoresponsive, peptide-based biomaterials with desired morphologies using strand length and dispersity as a handle for tuning thermal stability and formation of supramolecular structures.

Combining simulations and experiments for the molecular engineering of multifunctional collagen mimetic peptide-based materials.

Assembling peptides allow the creation of structurally complex materials, where amino acid selection influences resulting properties. We present a synergistic approach of experiments and simulations for examining the influence of natural and non-natural amino acid substitutions via incorporation of charged residues and a reactive handle on the thermal stability and assembly of multifunctional collagen mimetic peptides (CMPs). Experimentally, we observed inclusion of charged residues significantly decreased the melting temperature of CMP triple helices with further destabilization upon inclusion of the reactive handle. Atomistic simulations of a single CMP triple helix in explicit water showed increased residue-level and helical structural fluctuations caused by the inclusion of the reactive handle; however, these atomistic simulations cannot be used to predict changes in CMP melting transition. Coarse-grained (CG) simulations of CMPs at experimentally relevant solution conditions, showed, qualitatively, the same trends as experiments in CMP melting transition temperature with CMP design. These simulations show that when charged residues are included electrostatic repulsions significantly destabilize the CMP triple helix and that an additional inclusion of a reactive handle does not significantly change the melting transition. Based on findings from both experiments and simulations, the sequence design was refined for increased CMP triple helix thermal stability, and the reactive handle was utilized for the incorporation of the assembled CMPs within covalently crosslinked hydrogels. Overall, a unique approach was established for predicting stability of CMP triple helices for various sequences prior to synthesis, providing molecular insights for sequence design towards the creation of bulk nanostructured soft biomaterials.

Effect of Peptide Sequence on the LCST-Like Transition of Elastin-Like Peptides and Elastin-Like Peptide-Collagen-Like Peptide Conjugates: Simulations and Experiments.

Elastin-like polypeptides (ELPs) are thermoresponsive biopolymers that undergo an LCST-like phase transition in aqueous solutions. The temperature of this LCST-like transition, Tt , can be tuned by varying the number of repeat units in the ELP, sequence and composition of the repeat units, the solution conditions, and via conjugation to other biomacromolecules. In this study, we show how and why the choice of guest (X) residue in the VPGXG pentad repeat tunes the Tt of short ELPs, (VPGXG)4, in the free state and when conjugated to collagen-like peptides (CLPs). In experiments, the (VPGWG)4 chain (in short, WWWW) has a Tt < 278 K, while (VPGFG)4 or FFFF has a Tt > 353 K in both free ELP and ELP-CLP systems. The Tt for the FWWF ELP sequence decreases from being >353 K for free ELP to <278 K for the corresponding ELP-CLP system. The decrease in Tt upon conjugation to CLP has been shown to be due to the crowding of ELP chains that decreases the entropic loss upon ELP aggregation. Even though the net hydrophobicity of ELP has been reasoned to drive the Tt , the origins of lower Tt of WWWW compared to FFFF are unclear, as there is disagreement in hydrophobicity scales in how phenylalanine (F) compares to tryptophan (W). Motivated by these experimental observations, we use a combination of atomistic and coarse-grained (CG) molecular dynamics simulations. Atomistic simulations of free and tethered ELPs show that WWWW are more prone to acquire ß-turn structures than FFFF at lower temperatures. Also, the atomistically informed CG simulations show that the increased local stiffness in W than F due to the bulkier side chain in W compared to F, alone does not cause the shift in the transition of WWWW versus FFFF. The experimentally observed lower Tt of WWWW than FFFF is achieved in CG simulations only when the CG model incorporates both the atomistically informed local stiffness and stronger effective attractions localized at the W position versus the F position. The effective interactions localized at the guest residue in the CG model is guided by our atomistically observed increased propensity for ß-turn structure in WWWW versus FFFF and by past experimental work of Urry et al. quantifying hydrophobic differences through enthalpy of association for W versus F.

Proximal humeral fracture fixation: a biomechanical comparison of two constructs.

BACKGROUND: Different options exist for stabilizing proximal humeral fractures. This study compared the mechanical stability of 2 common proximal humeral fixation plates in bending and torsion. METHODS: Tests were conducted on 40 synthetic and 10 matched pairs of cadaveric humeri (evenly fixed with DePuy S3 proximal humeral plating system [DePuy Orthopaedics, Warsaw, IN, USA] and Synthes proximal humerus locking compression plate [Synthes, Paoli, PA, USA]). Half of the humeri were tested by cantilevered bending in flexion, extension, varus, and valgus for 100 cycles of ±5 mm of displacement at 1 mm/s before loading to failure in varus. The other half were tested in torsion for 100 cycles of ±8° of rotational displacement at 1°/s before loading to failure in external rotation. RESULTS: Peak cyclic loads for synthetic constructs were higher for DePuy plates than Synthes plates in varus and valgus (P < .0001), but a difference was not detected in extension (P > .40) or flexion (P = .0675). Peak cyclic loads for cadaveric constructs showed a significant difference in extension and flexion (Synthes > DePuy, P < .0001) and in varus (DePuy > Synthes, P < .05) but not in valgus (P > .10). Bending stiffness during varus failure testing was higher for DePuy plates than Synthes plates (P < .0001) for synthetic constructs. Regarding torsion of synthetic and cadaveric constructs, DePuy plates experienced higher peak cyclic torques over all cycles in both directions (P < .0001). For synthetic constructs, DePuy plates showed higher torsional stiffness in external failure than Synthes plates (P < .0001). CONCLUSIONS: The DePuy plate was stiffer than the Synthes plate with varus and valgus bending, as well as in torsion. The Synthes plate tended to be stiffer in flexion and extension.

Molecular Modeling and Simulations of Peptide-Polymer Conjugates.

Peptide-polymer conjugates are a class of soft materials composed of covalently linked blocks of protein/polypeptides and synthetic/natural polymers. These materials are practically useful in biological applications, such as drug delivery, DNA/gene delivery, and antimicrobial coatings, as well as nonbiological applications, such as electronics, separations, optics, and sensing. Given their broad applicability, there is motivation to understand the molecular and macroscale structure, dynamics, and thermodynamic behavior exhibited by such materials. We focus on the past and ongoing molecular simulation studies aimed at obtaining such fundamental understanding and predicting molecular design rules for the target function. We describe briefly the experimental work in this field that validates or motivates these computational studies. We also describe the various models (e.g., atomistic, coarse-grained, or hybrid) and simulation methods (e.g., stochastic versus deterministic, enhanced sampling) that have been used and the types of questions that have been answered using these computational approaches.

Placement of Tyrosine Residues as a Design Element for Tuning the Phase Transition of Elastin-peptide-containing Conjugates: Experiments and Simulations.

Elastin-like polypeptides (ELP) have been widely used in the biomaterials community due to their controllable, thermoresponsive properties and biocompatibility. Motivated by our previous work on the effect of tryptophan (W) substitutions on the LCST-like transitions of short ELPs, we studied a series of short ELPs containing tyrosine (Y) and/or phenylalanine (F) guest residues with only 5 or 6 pentapeptide repeat units. A combination of experiments and molecular dynamics (MD) simulations illustrated that the substitution of F with Y guest residues impacted the transition temperature (Tt) of short ELPs when conjugated to collagen-like-peptides (CLP), with a reduction in the transition temperature observed only after substitution of at least two residues. Placement of the Y residues near the N-terminal end of the ELP, away from the tethering point to the CLP, resulted in a lower Tt than that observed for peptides with the Y residues near the tethering point. Atomistic and coarse-grained MD simulations indicated an increase in intra- and inter- peptide hydrogen bonds in systems containing Y guest residues that are suggested to enhance the ability of the peptides to coacervate, with a concomitantly lower Tt. Simulations also revealed that the placement of Y-containing pentads near the N-terminus (i.e., away from CLP tethering point) versus C-terminus of the ELP led to more π-π stacking interactions at low temperatures, in agreement with our experimental observations of a lower Tt. Overall, this study provides mechanistic insights into the driving forces for the LCST-like transitions of ELPs and offers additional means for tuning the Tt of short ELPs for biomedical applications such as on-demand drug delivery and tissue engineering.

Parallel Plating of Simulated Distal Humerus Fractures Demonstrates Increased Stiffness Relative to Orthogonal Plating With a Distal Humerus Locking Plate System.

OBJECTIVES: This study compared the stiffness of precontoured parallel and orthogonal locking plate configurations in cyclic torsion and bending, and then extension to failure. METHODS: Tests were conducted on 9 matched pairs of cadaveric humeri. A 10 mm block of bone was excised from the distal humerus metaphysis to simulate comminution, and fractures were repaired in matched fashion using parallel or orthogonal Biomet ALPS distal humerus locking plates (Biomet, Inc, Warsaw, IN). Specimens were tested at 0.5 degree per second up to ±2 N-m in internal and external rotation for 20 cycles, then in flexion and extension cantilevered bending at 1 mm/s to ±50 N, followed by bending at 1 mm/s to failure in extension. RESULTS: Torsional stiffness of the parallel configuration group was greater than the orthogonal configuration in both internal and external rotation (P < 0.0001). Also, stiffness in bending was significantly greater in the parallel configuration group in both flexion and extension (P < 0.0001). In extension to failure testing, the parallel plate construct stiffness was significantly greater than the orthogonal configuration (P < 0.005). CONCLUSIONS: The parallel plate configuration demonstrated significantly greater stiffness than the orthogonal plate configuration in torsion and bending using locked distal humerus plates. This greater stiffness may prove desirable in the postoperative management of patients with comminuted distal humerus fractures, providing a stable anatomic reconstruction of the joint to allow early range of motion.