Development of a Coarse-Grained Model of Collagen-Like Peptide (CLP) for Studies of CLP Triple Helix Melting

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.

Investigating the Hydrogen Bond-Induced Self-Assembly of Polysulfamides Using Molecular Simulations and Experiments

In this paper, we present a synergistic, experimental, and computational study of the self-assembly of N,N'-disubstituted polysulfamides driven by hydrogen bonds (H-bonds) between the H-bonding donor and acceptor groups present in repeating sulfamides as a function of the structural design of the polysulfamide backbone. We developed a coarse-grained (CG) polysulfamide model that captures the directionality of H-bonds between the sulfamide groups and used this model in molecular dynamics (MD) simulations to study the self-assembly of these polymers in implicit solvent. The CGMD approach was validated by reproducing experimentally observed trends in the extent of crystallinity for three polysulfamides synthesized with aliphatic and/or aromatic repeating units. After validation of our CGMD approach, we computationally predicted the effect of repeat unit bulkiness, length, and uniformity of segment lengths in the polymers on the extent of orientational and positional order among the self-assembled polysulfamide chains, providing key design principles for tuning the extent of crystallinity in polysulfamides in experiments. Those computational predictions were then experimentally tested through the synthesis and characterization of polysulfamide architectures.

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.

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.

FEK self-assembled peptide hydrogels facilitate primary hepatocytes culture and pharmacokinetics screening

A FEFEFKFK (FEK, F, phenylalaninyl; E, glutamyl; K, lysinyl)-based self-assembling peptide hydrogel (FEK-SAPH) was developed to replace sandwich culture (SC) for improved culture of primary hepatocytes in vitro. Under neutral conditions, FEK self-assembles to form β-sheet nanofibers, which in turn form FEK-SAPH. For the culture of rat primary hepatocytes (RPH), the use of FEK-SAPH simplified operation steps and promoted excellent cell-cell interactions while maintaining the SC-related RPH polarity trend. Compared with SC, FEK-SAPH cultured RPH for 14 days, the bile duct network was formed, the secretion of albumin and urea was improved, and the metabolic clearance rate based on cytochrome P450 (CYPs) was comparable. In FEK-SAPH culture, the expression level of the biliary efflux transporter bile salt export pump increased by 230.7%, while the biliary excretion index value of deuterium-labeled sodium taurocholate (d8-TCA) differed slightly from the SC value (72% and 77%, respectively, p = .0195). The inhibitory effect of cholestasis drugs on FEK-SAPH was significantly higher than that of SC. In FEK-SAPH, hepatoprotective drugs were more effective in antagonizing hepatotoxicity induced by lithocholic acid (LCA). FEK-SAPH cultured RPH with hepatoprotective drugs can better recover from LCA-induced damage. In summary, FEK-SAPH can be used as a substitute for SC for pharmacokinetic screening to evaluate the drug absorption, disposition, metabolism, excretion, and toxicity (ADMET) in hepatocytes.

Biomaterials via peptide assembly: Design, characterization, and application in tissue engineering

A core challenge in biomaterials, with both fundamental significance and technological relevance, concerns the rational design of bioactive microenvironments. Designed properly, peptides can undergo supramolecular assembly into dynamic, physical hydrogels that mimic the mechanical, topological, and biochemical features of native tissue microenvironments. The relatively facile, inexpensive, and automatable preparation of peptides, coupled with low batch-to-batch variability, motivates the expanded use of assembling peptide hydrogels for biomedical applications. Integral to realizing dynamic peptide assemblies as functional biomaterials for tissue engineering is an understanding of the molecular and macroscopic features that govern assembly, morphology, and biological interactions. In this review, we first discuss the design of assembling peptides, including primary structure (sequence), secondary structure (e.g., α-helix and β-sheets), and molecular interactions that facilitate assembly into multiscale materials with desired properties. Next, we describe characterization tools for elucidating molecular structure and interactions, morphology, bulk properties, and biological functionality. Understanding of these characterization methods enables researchers to access a variety of approaches in this ever-expanding field. Finally, we discuss the biological properties and applications of peptide-based biomaterials for engineering several important tissues. By connecting molecular features and mechanisms of assembling peptides to the material and biological properties, we aim to guide the design and characterization of peptide-based biomaterials for tissue engineering and regenerative medicine. STATEMENT OF SIGNIFICANCE: Engineering peptide-based biomaterials that mimic the topological and mechanical properties of natural extracellular matrices provide excellent opportunities to direct cell behavior for regenerative medicine and tissue engineering. Here we review the molecular-scale features of assembling peptides that result in biomaterials that exhibit a variety of relevant extracellular matrix-mimetic properties and promote beneficial cell-biomaterial interactions. Aiming to inspire and guide researchers approaching this challenge from both the peptide biomaterial design and tissue engineering perspectives, we also present characterization tools for understanding the connection between peptide structure and properties and highlight the use of peptide-based biomaterials in neural, orthopedic, cardiac, muscular, and immune engineering applications.

Computational Reverse Engineering Analysis for Scattering Experiments (CREASE) on Vesicles Assembled from Amphiphilic Macromolecular Solutions

In this paper we present the development and validation of the "Computational Reverse-Engineering Analysis for Scattering Experiments" (CREASE) method for analyzing scattering results from vesicle structures that are commonly found upon assembly of synthetic, biomimetic, or bioderived amphiphilic copolymers in solution. The two-step CREASE method takes the amphiphilic polymer chemistry and small-angle scattering intensity profile, I exp(q), as input and determines the vesicles' structural features on multiple length scales ranging from assembled vesicle wall's individual layer thicknesses to the monomer-level packing and distribution of polymer conformations. In the first step of CREASE, a genetic algorithm (GA) is used to determine the relevant vesicle dimensions from the input macromolecular solution information and I exp(q) by identifying the structure whose computed scattering profile best matches the input I exp(q). Then in the second step, the GA-determined dimensions are used for molecular reconstruction of the vesicle structure. To validate CREASE for vesicles, we test CREASE on input scattering intensity profiles generated mathematically (termed as in silico I exp(q) vs q) from a variety of vesicle sizes with known dimensions. We also test CREASE on in silico I exp(q) vs q generated from vesicles with dispersity in all relevant dimensions, resembling real experiments. After successful validation of CREASE, we compare the CREASE-determined dimensions against those obtained from the traditional approach of fitting the scattering intensity profile to relevant analytical model in SASVIEW package. We show that CREASE performs better than or as well as the core-multishell analytical model's fitting in SASVIEW in determining vesicle dimensions with dispersity. We also show that CREASE provides structural information beyond those possible from traditional scattering analysis using the core-multishell model, such as the distribution of solvophilic monomers between the vesicle wall's inner and outer layers in the vesicle wall and the chain-level packing within each vesicle layer.

Rapid Production of Multifunctional Self-Assembling Peptides for Incorporation and Visualization within Hydrogel Biomaterials

Peptides are of continued interest for therapeutic applications, from soluble and immobilized ligands that promote desired binding or uptake to self-assembled supramolecular structures that serve as scaffolds in vitro and in vivo. These applications require efficient and scalable synthetic approaches because of the large amounts of material that often are needed for studies of bulk material properties and their translation. In this work, we establish new methods for the synthesis, purification, and visualization of assembling peptides, with a focus on multifunctional collagen mimetic peptides (mfCMPs) relevant for formation and integration within hydrogel-based biomaterials. First, a methodical approach useful for the microwave-assisted synthesis of assembling peptide sequences prone to deletions was established, beginning with the identification of the deleted residues and their locations and followed by targeted use of dual chemistry couplings for those specific residues. Second, purification techniques that integrate the principles of heating and ion displacement with traditional chromatography and dialysis were implemented to improve separation and isolation of the desired multifunctional peptide product, which contained blocks for thermoresponsiveness and ionic interactions. Third, an approach for fluorescent labeling of these mfCMPs, which is orthogonal to their assembly and their covalent incorporation into a bulk hydrogel material, was established, allowing visualization of the resulting hierarchical fibrillar structures in three dimensions within hydrogels using confocal microscopy. The methods presented in this work allow the production of multifunctional peptides in scalable quantities and with minimal deletions, enabling future studies for better understanding of composition-structure-property relationships and for translating these biomaterials into a range of applications. Although mfCMPs are the focus of this work, the methods demonstrated could prove useful for other assembling peptide systems and for the production of peptides more broadly for therapeutic applications.

Molecular Modeling and Simulation of Polymer Nanocomposites with Nanorod Fillers

We present a coarse-grained (CG) molecular dynamics (MD) simulation study of polymer nanocomposites (PNCs) containing nanorods with homogeneous and patchy surface chemistry/functionalization, modeled with isotropic and directional nanorod-nanorod attraction, respectively. We show how the PNC morphology is impacted by the nanorod design (i.e., aspect ratio, homogeneous or patchy surface chemistry/functionalization) for nanorods with a diameter equal to the Kuhn length of the polymer in the matrix. For PNCs with 10 vol % nanorods that have an aspect ratio ≤5, we observe percolated morphology with directional nanorod-nanorod attraction and phase-separated (i.e., nanorod aggregation) morphology with isotropic nanorod-nanorod attraction. In contrast, for nanorods with higher aspect ratios, both types of attractions result in aggregated nanorods morphology due to the dominance of entropic driving forces that cause long nanorods to form orientationally ordered aggregates. For most PNCs with isotropic or directional nanorod-nanorod attractions, the average matrix polymer conformation is not perturbed by the inclusion of up to 20 vol % nanorods. The polymer chains in contact with nanorods (i.e., interfacial chains) are on average extended and statistically different from the conformations the matrix chains adopt in the pure melt state (with no nanorods); in contrast, the polymer chains far from nanorods (i.e., bulk chains) adopt the same conformations as the matrix chains adopt in the pure melt state. We also study the effect of other parameters, such as attraction strength, nanorod volume fraction, and matrix chain length, for PNCs with isotropic or directional nanorod-nanorod attractions. Collectively, our results provide valuable design rules to achieve specific PNC morphologies (i.e., dispersed, aggregated, percolated, and orientationally aligned nanorods) for various potential applications.

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.

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.

Monitoring the Hierarchical Evolution from a Double-Stranded Helix to a Well-Defined Microscopic Morphology Based on a Turbine-like Aromatic Molecule

1H-Indazolo[1,2-b]phthalazine-5,10-dione IPDD with an approximate turbine-like spatial structure, primary assembled double-stranded helices at the first level, was predicted by quantum chemical calculations and confirmed by atomic force microscopy. The higher-dimensional hierarchical architectures including fibrils, helical fibers, spherical shells, and porous prismatic structures were observed in sequence by the scanning electron microscopy technique. The final porous prismatic structures sensitive to NH₃ vapors have the potential to be applied in gas sensing and absorbing materials.

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.

100th Anniversary of Macromolecular Science Viewpoint: Modeling and Simulation of Macromolecules with Hydrogen Bonds: Challenges, Successes, and Opportunities

Macromolecular materials with directional interactions such as hydrogen bonds exhibit numerous attractive features in terms of structure, thermodynamics, and dynamics. Besides enabling precise tuning of desirable geometries in the assembled state (e.g., programmable coordination numbers depending on the valency of the directional interaction), mixing in a blend/composite through stabilization via hydrogen bonds between the various components, hydrogen bonds can also impart responsiveness to external stimuli (e.g., temperature, pH). In biomacromolecules (e.g., proteins, DNA, polysaccharides), hydrogen bonds play a key role in stabilizing secondary and tertiary structures, which in turn define the function of these macromolecules. In this Viewpoint, I present the challenges, successes, and opportunities for molecular modeling and simulations to conduct fundamental and application-focused research on macromolecular materials with hydrogen bonding interactions. The past successes and limitations of atomistic simulations are discussed first, followed by highlights from recent developments in coarse-grained modeling and their use in studies of (synthetic and biologically relevant) macromolecular materials. Model development focused on polynucleotides (e.g., DNA, RNA, etc.), polypeptides, polysaccharides, and synthetic polymers at experimentally relevant conditions are highlighted. This viewpoint ends with potential future directions for macromolecular modeling and simulations with other types of directional interactions beyond hydrogen bonding.

Discovery and design of soft polymeric bio-inspired materials with multiscale simulations and artificial intelligence

Establishing the “Materials 4.0” paradigm requires intimate knowledge of the virtual space in materials design.

Self-assemble peptide biomaterials and their biomedical applications

Inspired by self-assembling peptides found in native proteins, deliberately designed engineered peptides have shown outstanding biocompatibility, biodegradability, and extracellular matrix-mimicking microenvironments. Assembly of the peptides can be triggered by external stimuli, such as electrolytes, temperature, and pH. The formation of nanostructures and subsequent nanocomposite materials often occur under physiological conditions. The respective properties of side chains in each amino acids provide numerous sites for chemical modification and conjugation choices of the peptides, enabling various resulting supramolecular nanostructures and hydrogels with adjustable mechanical and physicochemical properties. Moreover, additional functionalities can be easily induced into the hydrogels, including shear-thinning, bioactivity, self-healing, and shape memory. It further broaden the scope of application of self-assemble peptide materials. This review outlines designs of self-assembly peptide (β-sheet, α-helix, collagen-like peptides, elastin-like polypeptides, and peptide amphiphiles) with potential additional functionalities and their biomedical applications in bioprinting, tissue engineering, and drug delivery.

Coarse-grained molecular dynamics simulations of α-1,3-glucan

In this paper we present a computational study of aggregation in aqueous solutions of α-1,3-glucan captured using a coarse-grained (CG) model that can be extended to other polysaccharides. This CG model captures atomistic geometry (i.e., relative placement of the hydrogen bonding donors and acceptors within the monomer) of the α-1,3-glucan monomer, the directional interactions due to the donor-acceptor hydrogen bonds, and their effect on aggregation of multiple α-1,3-glucan chains without the extensive computational resources needed for simulations with atomistic models. Using this CG model, we conduct molecular dynamics simulations to assess the effect of varying α-1,3-glucan chain length and hydrogen bond interaction strengths on the aggregation of multiple chains at finite concentrations in implicit solvent. We quantify the hydrogen bonding strength needed for multiple chains to aggregate, the distribution of inter- and intra-chain hydrogen bonds within the aggregate and in some cases, the shapes of the aggregate. We also explore the effect of substitution/silencing of some randomly selected or specific hydrogen bonding sites in the chain on the aggregation and aggregate structure. In the unmodified α-1,3-glucan solution, the inter-chain hydrogen bonds cause the chains to aggregate into sheets. Random silencing of hydrogen bonding donor sites only increases the hydrogen bond strength needed for aggregation but retains the same aggregate structure as the unmodified chains. Specific silencing of the hydrogen-bonding site on the C6 carbon leads to the chains aggregating into planar sheets that then fold over to form hollow cylinders at intermediate hydrogen bond strength - 4.7 to 5.3 kcal mol-1. These cylindrical aggregates assemble end-to-end to form larger aggregates at higher hydrogen bond strengths.

Conformational preferences and phase behavior of intrinsically disordered low complexity sequences: insights from multiscale simulations

While many proteins and protein regions utilize a complex repertoire of amino acids to achieve their biological function, a subset of protein sequences are enriched in a reduced set of amino acids. These so-called low complexity (LC) sequences, specifically intrinsically disordered variants of LC sequences, have been the focus of recent investigations owing to their roles in a range of biological functions, specifically phase separation. Computational studies of LC sequences have provided rich insights into their behavior both as individual proteins in dilute solutions and as the drivers and modulators of higher-order assemblies. Here, we review how simulations performed across distinct resolutions have provided different types of insights into the biological role of LC sequences.