Metal Cation-Binding Mechanisms of Q-Proline Peptoid Macrocycles in Solution

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.

Minimal Peptoid Dynamics Inform Self-Assembly Propensity

Peptoids are structural isomers of natural peptides, with side chain attachment at the amide nitrogen, conferring this class of compounds with the ability to access both cis and trans ω torsions as well as an increased diversity of ψ/φ states with respect to peptides. Sampling within these dimensions is controlled through side chain selection, and an expansive set of viable peptoid residues exists. It has been shown recently that "minimal" di- and tripeptoids with aromatic side chains can self-assemble into highly ordered structures, with size and morphological definition varying as a function of sequence pattern (e.g., XFF and FXF, where X = a nonaromatic peptoid monomer). Aromatic groups, such as phenylalanine, are regularly used in the design of minimal peptide assemblers. In recognition of this, and to draw parallels between these compounds classes, we have developed a series of descriptors for intramolecular dynamics of aromatic side chains to discern whether these dynamics, in a preassembly condition, can be related to experimentally observed nanoscale assemblies. To do this, we have built on the atomistic peptoid force field reported by Weiser and Santiso (CGenFF-WS) through the rigorous fitting of partial charges and the collation of Charmm General Force Field (CGenFF) parameters relevant to these systems. Our study finds that the intramolecular dynamics of side chains, for a given sequence, is dependent on the specific combination of backbone ω torsions and that homogeneity of sampling across these states correlates well with the experimentally observed ability to assemble into nanomorphologies with long-range order. Sequence patterning is also shown to affect sampling, in a manner consistent for both tripeptoids and tripeptides. Additionally, sampling similarities between the nanofiber forming tripeptoid, Nf-Nke-Nf in the cc state, and the nanotube forming dipeptide FF, highlight a structural motif which may be relevant to the emergence of extended linear assemblies. To assess these properties, a variety of computational approaches have been employed.

Addressing Challenges of Macrocyclic Conformational Sampling in Polar and Apolar Solvents: Lessons for Chameleonicity

We evaluated a workflow to reliably sample the conformational space of a set of 47 peptidic macrocycles. Starting from SMILES strings, we use accelerated molecular dynamics simulations to overcome high energy barriers, in particular, the cis-trans isomerization of peptide bonds. We find that our approach performs very well in polar solvents like water and dimethyl sulfoxide. Interestingly, the protonation state of a secondary amine in the ring only slightly influences the conformational ensembles of our test systems. For several of the macrocycles, determining the conformational distribution in chloroform turns out to be considerably more challenging. Especially, the choice of partial charges crucially influences the ensembles in chloroform. We address these challenges by modifying initial structures and the choice of partial charges. Our results suggest that special care has to be taken to understand the configurational distribution in apolar solvents, which is a key step toward a reliable prediction of membrane permeation of macrocycles and their chameleonic properties.

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.

Evaluating the Conformations and Dynamics of Peptoid Macrocycles

Peptoid macrocycles are versatile and chemically diverse peptidomimetic oligomers. However, the conformations and dynamics of these macrocycles have not been evaluated comprehensively and require extensive further investigation. Recent studies indicate that two degrees of freedom, and four distinct conformations, adequately describe the behavior of each monomer backbone unit in most peptoid oligomers. On the basis of this insight, we conducted molecular dynamics simulations of model macrocycles using an exhaustive set of idealized possible starting conformations. Simulations of various sizes of peptoid macrocycles yielded a limited set of populated conformations. In addition to reproducing all relevant experimentally determined conformations, the simulations accurately predicted a cyclo-octamer conformation for which we now present the first experimental observation. Sets of three adjacent dihedral angles (ϕ_i, ψ_i, ω_i+1) exhibited correlated crankshaft motions over the course of simulation for peptoid macrocycles of six residues and larger. These correlated motions may occur in the form of an inversion of one amide bond and the concerted rotation of the preceding ϕ and ψ angles to their mirror-image conformation, a variation on "crankshaft flip" motions studied in polymers and peptides. The energy landscape of these peptoid macrocycles can be described as a network of conformations interconnected by transformations of individual crankshaft flips. For macrocycles of up to eight residues, our mapping of the landscape is essentially complete.