Modeling Sequence-Specific Polymers Using Anisotropic Coarse-Grained Sites Allows Quantitative Comparison with Experiment

Uprising Unconventional Nanobiomaterials: Peptoid Nanosheets as a Multi-Modular Platform for Advanced Biological Studies

Peptoids are bio-inspired peptidomimetic polymers that can be designed to self-assemble into a variety of nanostructures. Among these different assemblies, peptoid nanosheets are the most studied. Peptoid nanosheets are 2D highly ordered nanostructures, able to free float in aqueous solutions while featuring versatile chemical displays that can be tuned to incorporate a plethora of functional units. In this review, the synthetic approach used to prepare sequence-defined oligomers and highlight their main characteristics is introduced. The ability of peptoids to fold into nanostructures is then reviewed with an extensive emphasis on peptoid nanosheets, and their physico-chemical characteristics, assembly mechanism, and stability. A particular focus is also placed on the variety of functionalization incorporated into the peptoid nanosheets to tune their properties toward specific applications, especially within the fields of biology and medicine. Finally, the comparison between peptoid nanosheets and other 2D nanomaterials is discussed to address the challenges in the current nanomaterials and underline the future development of peptoid nanosheets in the field of biology.

Martinoid: the peptoid martini force field

Many exciting innovations have been made in the development of assembling peptoid materials. Typically, these have utilised large oligomeric sequences, though elsewhere the study of peptide self-assembly has yielded numerous examples of assemblers below 6-8 residues in length, evidencing that minimal peptoid assemblers are not only feasible but expected. A productive means of discovering such materials is through the application of in silico screening methods, which often benefit from the use of coarse-grained molecular dynamics (CG-MD) simulations. At the current level of development, CG models for peptoids are insufficient and we have been motivated to develop a Martini forcefield compatible peptoid model. A dual bottom-up and top-down parameterisation approach has been adopted, in keeping with the Martini parameterisation methodology, targeting the reproduction of atomistic MD dynamics and trends in experimentally obtained log D7.4 partition coefficients, respectively. This work has yielded valuable insights into the practicalities of parameterising peptoid monomers. Additionally, we demonstrate that our model can reproduce the experimental observations of two very different peptoid assembly systems, namely peptoid nanosheets and minimal tripeptoid assembly. Further we can simulate the peptoid helix secondary structure relevant for antimicrobial sequences. To be of maximum usefulness to the peptoid research community, we have developed freely available code to generate all requisite simulation files for the application of this model with Gromacs MD software.

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.

Systematic bottom-up molecular coarse-graining via force and torque matching using anisotropic particles

We derive a systematic and general method for parametrizing coarse-grained molecular models consisting of anisotropic particles from fine-grained (e.g. all-atom) models for condensed-phase molecular dynamics simulations. The method, which we call anisotropic force-matching coarse-graining (AFM-CG), is based on rigorous statistical mechanical principles, enforcing consistency between the coarse-grained and fine-grained phase-space distributions to derive equations for the coarse-grained forces, torques, masses, and moments of inertia in terms of properties of a condensed-phase fine-grained system. We verify the accuracy and efficiency of the method by coarse-graining liquid-state systems of two different anisotropic organic molecules, benzene and perylene, and show that the parametrized coarse-grained models more accurately describe properties of these systems than previous anisotropic coarse-grained models parametrized using other methods that do not account for finite-temperature and many-body effects on the condensed-phase coarse-grained interactions. The AFM-CG method will be useful for developing accurate and efficient dynamical simulation models of condensed-phase systems of molecules consisting of large, rigid, anisotropic fragments, such as liquid crystals, organic semiconductors, and nucleic acids.

Assessment of an anisotropic coarse-grained model for cis-1,4-polybutadiene: a bottom-up approach

The spherical representation usually utilized for the coarse-grained particles of soft matter systems is an assumption and pertinent studies have shown that both structural and dynamical properties can depend on anisotropic effects. On these grounds, we develop coarse-grained equations of motion which take into account explicitly the anisotropy of the beads. As a first step, this model incorporates only conservative terms. Inclusion of the dissipative and random terms is in principle possible but is beyond the scope of this study. The translational dynamics of the beads is tracked using the position and momentum of their center of mass, while their rotational dynamics is modeled by representing their orientation through the use of quaternions, similarly to the case of rigid bodies. The associated force and torque controlling the motion are derived from atomistic molecular dynamics (MD) simulations via a bottom-up approach and define a coarse-grained potential. The assumptions of the model are clearly stated and checked for a reference system of a cis-1,4-polybutadiene melt. In particular, the choice of the angular velocity as a slow variable is justified by comparing its dynamics to atomic vibrations. The accuracy of this approach to reproduce static structural features of the polymer melt is assessed.

Self-Assembly of Minimal Peptoid Sequences

Peptoids are biofunctional N-substituted glycine peptidomimics. Their self-assembly is of fundamental interest because they demonstrate alternatives to conventional peptide structures based on backbone chirality and beta-sheet hydrogen bonding. The search for self-assembling, water-soluble "minimal" sequences, be they peptide or peptidomimic, is a further challenge. Such sequences are highly desired for their compatibility with biomacromolecules and convenient synthesis for broader application. We report the self-assembly of a set of trimeric, water-soluble α-peptoids that exhibit a relatively low critical aggregation concentration (CAC ∼ 0.3 wt %). Cryo-EM and angle-resolved DLS show different sequence-dependent morphologies, namely uniform ca. 6 nm wide nanofibers, sheets, and clusters of globular assemblies. Absorbance and fluorescence spectroscopies indicate unique phenyl environments for π-interactions in the highly ordered nanofibers. Assembly of our peptoids takes place when the sequences are fully ionized, representing a departure from superficially similar amyloid-type hydrogen-bonded peptide nanostructures and expanding the horizons of assembly for sequence-specific bio- and biomimetic macromolecules.

Generic Model for Lamellar Self-Assembly in Conjugated Polymers: Linking Mesoscopic Morphology and Charge Transport in P3HT

We develop a generic coarse-grained model of soluble conjugated polymers, capable of describing their self-assembly into a lamellar mesophase. Polymer chains are described by a hindered-rotation model, where interaction centers represent entire repeat units, including side chains. We introduce soft anisotropic nonbonded interactions to mimic the potential of mean force between atomistic repeat units. The functional form of this potential reflects the symmetry of the molecular order in a lamellar mesophase. The model can generate both nematic and lamellar (sanidic smectic) molecular arrangements. We parametrize this model for a soluble conjugated polymer poly(3-hexylthiophene) (P3HT) and demonstrate that the simulated lamellar mesophase matches morphologies of low molecular weight P3HT, experimentally observed at elevated temperatures. A qualitative charge-transport model allows us to link local chain conformations and mesoscale order to charge transport. In particular, it shows how coarsening of lamellar domains and chain extension increase the charge carrier mobility. By modeling large systems and long chains, we can capture transport between lamellar layers, which is due to rare, but thermodynamically allowed, backbone bridges between neighboring layers.

Molecular simulations of self-assembling bio-inspired supramolecular systems and their connection to experiments

In bionanotechnology, the field of creating functional materials consisting of bio-inspired molecules, the function and shape of a nanostructure only appear through the assembly of many small molecules together. The large number of building blocks required to define a nanostructure combined with the many degrees of freedom in packing small molecules has long precluded molecular simulations, but recent advances in computational hardware as well as software have made classical simulations available to this strongly expanding field. Here, we review the state of the art in simulations of self-assembling bio-inspired supramolecular systems. We will first discuss progress in force fields, simulation protocols and enhanced sampling techniques using recent examples. Secondly, we will focus on efforts to enable the comparison of experimentally accessible observables and computational results. Experimental quantities that can be measured by microscopy, spectroscopy and scattering can be linked to simulation output either directly or indirectly, via quantum mechanical or semi-empirical techniques. Overall, we aim to provide an overview of the various computational approaches to understand not only the molecular architecture of nanostructures, but also the mechanism of their formation.

Design, Synthesis, Assembly, and Engineering of Peptoid Nanosheets

Two-dimensional (2D) atomically defined organic nanomaterials are an important material class with broad applications. However, few general synthetic methods exist to produce such materials in high yields and to precisely functionalize them. One strategy to form ordered 2D organic nanomaterials is through the supramolecular assembly of sequence-defined synthetic polymers. Peptoids, one such class of polymer, are designable bioinspired heteropolymers whose main-chain length and monomer sequence can be precisely controlled. We have recently discovered that individual peptoid polymers with a simple sequence of alternating hydrophobic and ionic monomers can self-assemble into highly ordered, free-floating nanosheets. A detailed understanding of their molecular structure and supramolecular assembly dynamics provides a robust platform for the discovery of new classes of nanosheets with tunable properties and novel applications. In this Account, we discuss the discovery, characterization, assembly, molecular modeling, and functionalization of peptoid nanosheets. The fundamental properties of peptoid nanosheets, their mechanism of formation, and their application as robust scaffolds for molecular recognition and as templates for the growth of inorganic minerals have been probed by an arsenal of experimental characterization techniques (e.g., scanning probe, electron, and optical microscopy, X-ray diffraction, surface-selective vibrational spectroscopy, and surface tensiometry) and computational techniques (coarse-grained and atomistic modeling). Peptoid nanosheets are supramolecular assemblies of 16-42-mer chains that form molecular bilayers. They span tens of microns in lateral dimensions and freely float in water. Their component chains are highly ordered, with chains nearly fully extended and packed parallel to one another as a result of hydrophobic and electrostatic interactions. Nanosheets form via a novel interface-catalyzed monolayer collapse mechanism. Peptoid chains first assemble into a monolayer at either an air-water or oil-water interface, on which peptoid chains extend, order, and pack into a brick-like pattern. Upon mechanical compression of the interface, the monolayer buckles into stable bilayer structures. Recent work has focused on the design of nanosheets with tunable properties and functionality. They are readily engineerable, as functional monomers can be readily incorporated onto the nanosheet surface or into the interior. For example, functional hydrophilic "loops" have been displayed on the surfaces of nanosheets. These loops can interact with specific protein targets, serving as a potentially general platform for molecular recognition. Nanosheets can also bind metal ions and serve as 2D templates for mineral growth. Through our understanding of the formation mechanism, along with predicted features ascertained from molecular modeling, we aim to further design and synthesize nanosheets as robust protein mimetics with the potential for unprecedented functionality and stability.

Insights into Peptoid Helix Folding Cooperativity from an Improved Backbone Potential

Peptoids (N-substituted oligoglycines) are biomimetic polymers that can fold into a variety of unique structural scaffolds. Peptoid helices, which result from the incorporation of bulky chiral side chains, are a key peptoid structural motif whose formation has not yet been accurately modeled in molecular simulations. Here, we report that a simple modification of the backbone φ-angle potential in GAFF is able to produce well-folded cis-amide helices of (S)-N-(1-phenylethyl)glycine (Nspe), consistent with experiment. We validate our results against both QM calculations and NMR experiments. For this latter task, we make quantitative comparisons to sparse NOE data using the Bayesian Inference of Conformational Populations (BICePs) algorithm, a method we have recently developed for this purpose. We then performed extensive REMD simulations of Nspe oligomers as a function of chain length and temperature to probe the molecular forces driving cooperative helix formation. Analysis of simulation data by Lifson-Roig helix-coil theory show that the modified potential predicts much more cooperative folding for Nspe helices. Unlike peptides, per-residue entropy changes for helix nucleation and extension are mostly positive, suggesting that steric bulk provides the main driving force for folding. We expect these results to inform future work aimed at predicting and designing peptoid peptidomimetics and tertiary assemblies of peptoid helices.

Sequence Programmable Peptoid Polymers for Diverse Materials Applications

Polymer sequence programmability is required for the diverse structures and complex properties that are achieved by native biological polymers, but efforts towards controlling the sequence of synthetic polymers are, by comparison, still in their infancy. Traditional polymers provide robust and chemically diverse materials, but synthetic control over their monomer sequences is limited. The modular and step-wise synthesis of peptoid polymers, on the other hand, allows for precise control over the monomer sequences, affording opportunities for these chains to fold into well-defined nanostructures. Hundreds of different side chains have been incorporated into peptoid polymers using efficient reaction chemistry, allowing for a seemingly infinite variety of possible synthetically accessible polymer sequences. Combinatorial discovery techniques have allowed the identification of functional polymers within large libraries of peptoids, and newly developed theoretical modeling tools specifically adapted for peptoids enable the future design of polymers with desired functions. Work towards controlling the three-dimensional structure of peptoids, from the conformation of the amide bond to the formation of protein-like tertiary structure, has and will continue to enable the construction of tunable and innovative nanomaterials that bridge the gap between natural and synthetic polymers.

Structure-determining step in the hierarchical assembly of peptoid nanosheets

Organic two-dimensional nanomaterials are of growing importance, yet few general synthetic methods exist to produce them in high yields and to precisely functionalize them. We previously developed an efficient hierarchical supramolecular assembly route to peptoid bilayer nanosheets, where the organization of biomimetic polymer sequences is catalyzed by an air-water interface. Here we determine at which stages of assembly the nanoscale and atomic-scale order appear. We used X-ray scattering, grazing incidence X-ray scattering at the air-water interface, electron diffraction, and a recently developed computational coarse-grained peptoid model to probe the molecular ordering at various stages of assembly. We found that lateral packing and organization of the chains occurs during the formation of a peptoid monolayer, prior to its collapse into a bilayer. Identifying the structure-determining step enables strategies to influence nanosheet order, to predict and optimize production yields, and to further engineer this class of material. More generally, our results provide a guide for using fluid interfaces to catalytically assemble 2D nanomaterials.