Interaction energy between neuronal cell receptors and peptoid ligands

Development of a Systematic and Extensible Force Field for Peptoids (STEPs)

Peptoids (N-substituted glycines) are a class of biomimetic polymers that have attracted significant attention due to their accessible synthesis and enzymatic and thermal stability relative to their naturally occurring counterparts (polypeptides). While these polymers provide the promise of more robust functional materials via hierarchical approaches, they present a new challenge for computational structure prediction for material design. The reliability of calculations hinges on the accuracy of interactions represented in the force field used to model peptoids. For proteins, structure prediction based on sequence and de novo design has made dramatic progress in recent years; however, these models are not readily transferable for peptoids. Current efforts to develop and implement peptoid-specific force fields are spread out, leading to replicated efforts and a fragmented collection of parameterized sidechains. Here, we developed a peptoid-specific force field containing 70 different side chains, using GAFF2 as starting point. The new model is validated based on the generation of Ramachandran-like plots from DFT optimization compared against force field reproduced potential energy and free energy surfaces as well as the reproduction of equilibrium cis/trans values for some residues experimentally known to form helical structures. Equilibrium cis/trans distributions (Kct) are estimated for all parameterized residues to identify which residues have an intrinsic propensity for cis or trans states in the monomeric state.

A hybrid approach for coarse-graining helical peptoids: Solvation, secondary structure, and assembly

Protein mimics such as peptoids form self-assembled nanostructures whose shape and function are governed by the side chain chemistry and secondary structure. Experiments have shown that a peptoid sequence with a helical secondary structure assembles into microspheres that are stable under various conditions. The conformation and organization of the peptoids within the assemblies remains unknown and is elucidated in this study via a hybrid, bottom-up coarse-graining approach. The resultant coarse-grained (CG) model preserves the chemical and structural details that are critical for capturing the secondary structure of the peptoid. The CG model accurately captures the overall conformation and solvation of the peptoids in an aqueous solution. Furthermore, the model resolves the assembly of multiple peptoids into a hemispherical aggregate that is in qualitative agreement with the corresponding results from experiments. The mildly hydrophilic peptoid residues are placed along the curved interface of the aggregate. The composition of the residues on the exterior of the aggregate is determined by two conformations adopted by the peptoid chains. Hence, the CG model simultaneously captures sequence-specific features and the assembly of a large number of peptoids. This multiscale, multiresolution coarse-graining approach could help in predicting the organization and packing of other tunable oligomeric sequences of relevance to biomedicine and electronics.

Predicting neurite extension for varying extracellular matrix stiffness and topography

Neurite extension is a dynamic process and is dependent on the microenvironment. The mechanical properties of the extracellular matrix (ECM), such as stiffness and topography influence the microenvironment and affects neurite extension; however, the mechanistic basis for this dynamic response of neurite extension remains elusive. In this study, we develop a computational model that predicts neurite extension dynamics process as the stiffness and patterned topography of ECM changes. The model includes the contribution of receptors integrin and neural cellular adhesion molecule toward the growth of neurite tip. We use non-linear finite element analysis (FEA) to model the neuronal cell, neurite, and the ECM, which is then coupled to the force-deformation receptor properties obtained from molecular dynamics simulations. Using an empirical relation, we develop a neurite extension algorithm that simulates the dynamic process of growth cone induced by growth cone extension, receptor density, and rupture. We investigate the dependence of neurite extension on ECM stiffness using three distinct materials, the effect of width and spacing of continuous (cylindrical) and discontinuous (pillar) patterned topography, as well as the topography steepness and stiffness gradient. We find that an increasing stiffness and width of patterned topography results in increased neurite extension, but the magnitude of the increase differs depending on the growth cone extension and receptor density between them. These findings will aid in vitro studies in determining an ECM with appropriate mechanical properties, such as stiffness and topography that will improve neurite extension, thus resulting in the formation of functional neurons.