Evaluating the Conformations and Dynamics of Peptoid Macrocycles

Enhancing molecular diversity of peptoid oligomers using amino acid synthons

We report the use of unprotected amino acids as submonomer reagents in the solid-phase synthesis of N-substituted glycine peptoid oligomers. Subsequent coupling of an amine, alcohol, or thiol to the free carboxylate of the incorporated amino acid provides access to peptoids bearing amides, esters, and thioesters as side chain pendant groups and permits further elongation of the peptoid backbone. The palette of readily obtained building blocks suitable for solid-phase peptoid synthesis is substantially expanded through this protocol, further enhancing the chemical diversity and potential applications of sequence-specific peptoid oligomers.

Exploration of Tertiary Structure in Sequence-Defined Polymers Using Molecular Dynamics Simulations

Peptoids are a class of sequence-defined biomimetic polymers with peptide-like backbones and side chains located on backbone nitrogens rather than alpha carbons. These materials demonstrate a strong ability for precise control of single-chain structure, multiunit self-assembly, and macromolecular assembly through careful tuning of sequence due to the diversity of available side chains, although the driving forces behind these assemblies are often not understood. Prior experimental work has shown that linked 15mer peptoids can mimic the protein helical hairpin structure by leveraging the chirality-inducing nature of bulky side chains and hydrophobicity, but there are still gaps in our understanding of the relationship between sequence, stability, and particular secondary or tertiary structure. We present a molecular dynamics (MD) study on the folding behavior of these polymers into hairpins, discussing the differences in structure from sequences with various characteristics in water and acetonitrile, and then compare the handedness preference of common helical motifs between solvents.

Design of Triazolium-Grafted Peptidomimetic Macrocycles with Facial Amphipathicity to Target Pathogenic Bacteria

Access to 1,2,3-triazolium-grafted peptoid macrocycles was developed by macrocyclization and multivalent postmodification of linear peptoid oligomers carrying an alternance of benzylic and propargyl groups as side chains. X-ray analysis and NMR studies revealed a conformational preference for constrained hairpin-shaped structures leading to the facial amphipathic character of these macrocycles. A preliminary evaluation showed the antimicrobial activities of these new cationic amphipathic architectures.

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.

Water-Soluble Chiral Cyclic Peptoids and Their Sodium and Gadolinium Complexes: Study of Conformational and Relaxometric Properties

Cyclic peptoids are macrocyclic oligomers of N-substituted glycines with specific folding abilities and excellent metal binding properties. In this work, we show how strategic positioning of chiral (S)- and (R)-(1-carboxyethyl)glycine units influences the conformational stability of water-soluble macrocyclic peptoids as sodium complexes. The reported results are based on nuclear magnetic resonance spectroscopy, extensive computational studies, and X-ray diffraction analysis using single crystals grown from aqueous solutions. The studies include ¹H relaxometric investigations of hexameric cyclic peptoids in the presence of the Gd³⁺ ion to assess their thermodynamic stabilities and relaxivities.

A Model for the Rapid Assessment of Solution Structures for 24-Atom Macrocycles: The Impact of β-Branched Amino Acids on Conformation

Experiment and computation are used to develop a model to rapidly predict solution structures of macrocycles sharing the same Murcko framework. These 24-atom triazine macrocycles result from the quantitative dimerization of identical monomers presenting a hydrazine group and an acetal tethered to an amino acid linker. Monomers comprising glycine and the β-branched amino acids threonine, valine, and isoleucine yield macrocycles G-G, T-T, V-V, and I-I, respectively. Elements common to all members of the framework include the efficiency of macrocyclization (quantitative), the solution- and solid-state structures (folded), the site of protonation (opposite the auxiliary dimethylamine group), the geometry of the hydrazone (E), the C2 symmetry of the subunits (conserved), and the rotamer state adopted. In aggregate, the data reveal metrics predictive of the three-dimensional solution structure that derive from the fingerprint region of the 1D 1H spectrum and a network of rOes from a single resonance. The metrics also afford delineation of more nuanced structural features that allow subpopulations to be identified among the members of the framework. Well-tempered metadynamics provides free energy surfaces and population distributions of these macrocycles. The areas of the free energy surface decrease with increasing steric bulk (G-G > V-V ∼ T-T > I-I). In addition, the surfaces are increasingly isoenergetic with decreasing steric bulk (G-G > V-V ∼ T-T > I-I).

Folding Coarse-Grained Oligomer Models with PyRosetta

Non-biological foldamers are a promising class of macromolecules that share similarities to classical biopolymers such as proteins and nucleic acids. Currently, designing novel foldamers is a non-trivial process, often involving many iterations of trial synthesis and characterization until folded structures are observed. In this work, we aim to tackle these foldamer design challenges using computational modeling techniques. We developed CG PyRosetta, an extension to the popular protein folding python package, PyRosetta, which introduces coarse-grained (CG) residues into PyRosetta, enabling the folding of toy CG foldamer models. Although these models are simplified, they can help explore overarching physical hypotheses about how oligomers can form. Through systematic variation of CG parameters in these models, we can investigate various folding hypotheses at the CG scale to inform the design process of new foldamer chemistries. In this study, we demonstrate CG PyRosetta's ability to identify minimum energy structures with a diverse structural search over a range of simple models, as well as two hypothesis-driven parameter scans investigating the effects of side-chain size and internal backbone angle on secondary structures. We are able to identify several types of secondary structures from single- and double-helices to sheet-like and knot-like structures. We show how side-chain size and backbone bond angle both play an important role in the structure and energetics of these toy models. Optimal side-chain sizes promote favorable packing of side chains, while specific backbone bond angles influence the specific helix type found in folded structures.