PEPDROID: Development of a Generic DREIDING-Based Force Field for the Assessment of Peptoid Secondary Structures

Switchable supramolecular helices for asymmetric stereodivergent catalysis

Despite recent developments on the design of dynamic catalysts, none of them have been exploited for the in-situ control of multiple stereogenic centers in a single molecular scaffold. We report herein that it is possible to obtain in majority any amongst the four possible stereoisomers of an amino alcohol by means of a switchable asymmetric catalyst built on supramolecular helices. Hydrogen-bonded assemblies between a benzene-1,3,5-tricarboxamide (BTA) achiral phosphine ligand coordinated to copper and a chiral BTA comonomer are engaged in a copper-hydride catalyzed hydrosilylation and hydroamination cascade process. The nature of the product stereoisomer is related to the handedness of the helices and can thus be directed in a predictable way by changing the nature of the major enantiomer of the BTA comonomer present in the assemblies. The strategy allows all stereoisomers to be obtained one-pot with similar selectivities by conducting the cascade reaction in a concomitant manner, i.e. without inverting the handedness of the helices, or sequentially, i.e. by switching the handedness of the supramolecular helices between the hydrosilylation and hydroamination steps. Supramolecular helical catalysts appear as a unique and versatile platform to control the configuration of molecules or polymers embedding several stereogenic centers.

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.

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.

Gas-phase structure of polymer ions: Tying together theoretical approaches and ion mobility spectrometry

An increasing number of studies take advantage of ion mobility spectrometry (IMS) coupled to mass spectrometry (IMS-MS) to investigate the spatial structure of gaseous ions. Synthetic polymers occupy a unique place in the field of IMS-MS. Indeed, due to their intrinsic dispersity, they offer a broad range of homologous ions with different lengths. To help rationalize experimental data, various theoretical approaches have been described. First, the study of trend lines is proposed to derive physicochemical and structural parameters. However, the evaluation of data fitting reflects the overall behavior of the ions without reflecting specific information on their conformation. Atomistic simulations constitute another approach that provide accurate information about the ion shape. The overall scope of this review is dedicated to the synergy between IMS-MS and theoretical approaches, including computational chemistry, demonstrating the essential role they play to fully understand/interpret IMS-MS data.

Computational and Experimental Determination of the Properties, Structure, and Stability of Peptoid Nanosheets and Nanotubes

Peptoids (N-substituted glycines) are a group of highly controllable peptidomimetic polymers. Amphiphilic diblock peptoids have been engineered to assemble crystalline nanospheres, nanofibrils, nanosheets, and nanotubes with biochemical, biomedical, and bioengineering applications. The mechanical properties of peptoid nanoaggregates and their relationship to the emergent self-assembled morphologies have been relatively unexplored and are critical for the rational design of peptoid nanomaterials. In this work, we consider a family of amphiphilic diblock peptoids consisting of a prototypical tube-former (Nbrpm6Nc6, a NH₂-capped hydrophobic block of six N-((4-bromophenyl)methyl)glycine residues conjugated to a polar NH₃(CH₂)₅CO tail), a prototypical sheet-former (Nbrpe6Nc6, where the hydrophobic block comprises six N-((4-bromophenyl)ethyl)glycine residues), and an intermediate sequence that forms mixed structures ((NbrpeNbrpm)3Nc6). We combine all-atom molecular dynamics simulations and atomic force microscopy to determine the mechanical properties of the self-assembled 2D crystalline nanosheets and relate these properties to the observed self-assembled morphologies. We find good agreement between our computational predictions and experimental measurements of Young's modulus of crystalline nanosheets. A computational analysis of the bending modulus along the two axes of the planar crystalline nanosheets reveals bending to be more favorable along the axis in which the peptoids stack by interdigitation of the side chains compared to that in which they form columnar crystals with π-stacked side chains. We construct molecular models of nanotubes of the Nbrpm6Nc6 tube-forming peptoid and predict a stability optimum in good agreement with experimental measurements. A theoretical model of nanotube stability suggests that this optimum is a free energy minimum corresponding to a "Goldilocks" tube radius at which capillary wave fluctuations in the tube wall are minimized.

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.

Is an Inductive Effect Explicit Account Required for Atomic Charges Aimed at Use within the Force Fields?

Polarization and inductive effects are the concepts that have been widely used in qualitative and even quantitative descriptions of experimentally observed properties in chemistry. The polarization effect has proven to be important in cases of biomolecular modeling though still the vast majority of molecular simulations use the classical non-polarizable force fields. In the last few decades, a lot of effort has been put into promoting the polarization effect and incorporating it into modern force fields and charge calculation methods. In contrast, the inductive effect has not attracted such attention and is effectively absent in both classic and modern force fields. Thus, a question is whether this difference corresponds to the difference in the physical significance of the effects and their explicit account, or is an artifact that should be corrected in the next generation of force fields. The significance of the electronic effects is studied in this paper through the prism of performance of specific models for atomic charge calculation that take into explicit account a nested set of effects: the formal charge, the nearest neighbors, the inductive effect, and finally the model, which takes into account all effects, which are possible to account for using atomic charges. The specific choice for the methods is the following: formal charges, MMFF94 bond charge increments, Dynamic Electronegativity Relaxation (DENR), and RESP. We propose a special scheme for the separate estimation of each particular effect contribution. By pairwise comparing the residual molecular electrostatic potential (MEP) errors of those charge models (aimed at best reproducing the quantum chemical reference MEP), we sequentially revealed how the account of each effect contributes to the better-quality MEP reproduction. The following relative importance of effects was estimated; thus, the natural hierarchy of the effects was established. First, the account of formal charges is of primordial importance. Second, the nearest neighbors account is the next in significance. Third, the explicit account of inductive effect in empirical charge calculation schemes was shown to significantly─both qualitatively and quantitatively─improve the quality of MEP reproduction. Fourth, the contribution of polarization is indirectly assessed. Surprisingly, it is of the order of magnitude of the inductive effect even for the molecular systems, for which it is anticipated to be more significant. Finally, the relative importance of anisotropic effects in neutral molecules was additionally reviewed.

Simulation studies of polypeptoids using replica exchange with dynamical scaling and dihedral biasing

Polypeptoids differ from polypeptides in that the amide bond can more frequently adopt both cis and trans conformations. The transition between the two conformations requires overcoming a large energy barrier, making it difficult for conventional molecular simulations to adequately visit the cis and trans structures. A replica-exchange method is presented that allows for easy rotations of the amide bond and also an efficient linking to a high temperature replica. The method allows for just three replicas (one at the temperature and Hamiltonian of interest, a second high temperature replica with a biased dihedral potential, and a third connecting them) to overcome the amide bond sampling problem and also enhance sampling for other coordinates. The results indicate that for short peptoid oligomers, the conformations can range from all cis to all trans with an average cis/trans ratio that depends on side chain and potential model.

Hierarchical Self-Assembly Pathways of Peptoid Helices and Sheets

Peptoids (N-substituted glycines) are a class of tailorable synthetic peptidomic polymers. Amphiphilic diblock peptoids have been engineered to assemble 2D crystalline lattices with applications in catalysis and molecular separations. Assembly is induced in an organic solvent/water mixture by evaporating the organic phase, but the assembly pathways remain uncharacterized. We conduct all-atom molecular dynamics simulations of Nbrpe6Nc6 as a prototypical amphiphilic diblock peptoid comprising an NH₂-capped block of six hydrophobic N-((4-bromophenyl)ethyl)glycine residues conjugated to a polar NH₃(CH₂)₅CO tail. We identify a thermodynamically controlled assembly mechanism by which monomers assemble into disordered aggregates that self-order into 1D chiral helical rods then 2D achiral crystalline sheets. We support our computational predictions with experimental observations of 1D rods using small-angle X-ray scattering, circular dichroism, and atomic force microscopy and 2D crystalline sheets using X-ray diffraction and atomic force microscopy. This work establishes a new understanding of hierarchical peptoid assembly and principles for the design of peptoid-based nanomaterials.

Helical Peptoid Ions in the Gas Phase: Thwarting the Charge Solvation Effect by H-Bond Compensation

Folding and unfolding processes are key aspects that should be mastered for the design of foldamer molecules for targeted applications. In contrast to the solution phase, in vacuo conditions represent a well-defined environment to analyze the intramolecular interactions that largely control the folding/unfolding dynamics. Ion mobility mass spectrometry coupled to theoretical modeling represents an efficient method to decipher the spatial structures of gaseous ions, including foldamers. However, charge solvation typically compacts the ion structure in the absence of strong stabilizing secondary interactions. This is the case in peptoids that are synthetic peptide regioisomers whose side chains are connected to the nitrogen atoms of the backbone instead of α-carbon as in peptides, thus implying the absence of H-bonds among the core units of the backbone. A recent work indeed reported that helical peptoids based on Nspe units formed in solution do not retain their secondary structure when transferred to the gas phase upon electrospray ionization (ESI). In this context, we demonstrate here that the helical structure of peptoids bearing (S)-N-(1-carboxy-2-phenylethyl) bulky side chains (Nscp) is largely preserved in the gas phase by the creation of a hydrogen bond network, induced by the presence of carboxylic moieties, that compensates for the charge solvation process.

Peptoids as a Chiral Stationary Phase for Liquid Chromatography: Insights from Molecular Dynamics Simulations

Peptoids are peptide regioisomers with attractive structural tunability in terms of sequence and three-dimensional arrangement. Peptoids are foreseen to have a great potential for many diverse applications, including their utilization as a chiral stationary phase in chromatography. To achieve chiral recognition, a chiral side chain is required to allow specific interactions with a given enantiomer from a racemic mixture. One of the most studied chiral stationary phases, built with (S)-N-1-phenylethyl (Nspe) units, was shown to be successful in resolving racemic mixtures of binaphthyl derivatives. However, there is currently no description at the atomic scale of the factors favoring its enantioselectivity. Here, we take advantage of steered molecular dynamics simulations to mimic the elution process at the atomic scale and present evidence that the predominantly right-handed helical conformation of Nspe peptoids and their ability to form stronger hydrogen bonds with the (S) enantiomer are responsible for the chiral recognition of the popular chiral probe 2,2'-bihydroxy-1,1'-binaphthyl.

Modification and Validation of the DREIDING Force Field for Molecular Liquid Simulations (DREIDING-UT)

In this study, we modified Lennard-Jones (LJ) parameters and point-charge parameters of the DREIDING force field (the modified force-field model is named DREIDING-UT). While the original LJ parameters of DREIDING were derived through an analytical formula to reproduce the potential depths and the equilibrium lengths of the Buckingham potentials of DREIDING/X6, the modified LJ parameters were derived through the least square fitting of the Buckingham potentials. Because the Gasteiger-Marsili (GM) charges of the original DREIDING underestimated electrostatic interactions, we replaced it with the restrained electrostatic potential (RESP) charges calculated from the ab initio wavefunctions, taking the dynamic electron correlation and solvation effects into account. To confirm how the modified force field works, we conducted molecular dynamics (MD) simulations of typical liquids. It was found that the densities and self-diffusion coefficients of the DREIDING-UT model agreed with the experimental ones much better than those of the original model.

Peptide science: A "rule model" for new generations of peptidomimetics

Peptides have been heavily investigated for their biocompatible and bioactive properties. Though a wide array of functionalities can be introduced by varying the amino acid sequence or by structural constraints, properties such as proteolytic stability, catalytic activity, and phase behavior in solution are difficult or impossible to impart upon naturally occurring α-L-peptides. To this end, sequence-controlled peptidomimetics exhibit new folds, morphologies, and chemical modifications that create new structures and functions. The study of these new classes of polymers, especially α-peptoids, has been highly influenced by the analysis, computational, and design techniques developed for peptides. This review examines techniques to determine primary, secondary, and tertiary structure of peptides, and how they have been adapted to investigate peptoid structure. Computational models developed for peptides have been modified to predict the morphologies of peptoids and have increased in accuracy in recent years. The combination of in vitro and in silico techniques have led to secondary and tertiary structure design principles that mirror those for peptides. We then examine several important developments in peptoid applications inspired by peptides such as pharmaceuticals, catalysis, and protein-binding. A brief survey of alternative backbone structures and research investigating these peptidomimetics shows how the advancement of peptide and peptoid science has influenced the growth of numerous fields of study. As peptide, peptoid, and other peptidomimetic studies continue to advance, we will expect to see higher throughput structural analyses, greater computational accuracy and functionality, and wider application space that can improve human health, solve environmental challenges, and meet industrial needs. STATEMENT OF SIGNIFICANCE: Many historical, chemical, and functional relations draw a thread connecting peptides to their recent cognates, the "peptidomimetics". This review presents a comprehensive survey of this field by highlighting the width and relevance of these familial connections. In the first section, we examine the experimental and computational techniques originally developed for peptides and their morphing into a broader analytical and predictive toolbox. The second section presents an excursus of the structures and properties of prominent peptidomimetics, and how the expansion of the chemical and structural diversity has returned new exciting properties. The third section presents an overview of technological applications and new families of peptidomimetics. As the field grows, new compounds emerge with clear potential in medicine and advanced manufacturing.

Effects of electrospray mechanisms and structural relaxation on polylactide ion conformations in the gas phase: insights from ion mobility spectrometry and molecular dynamics simulations

Gas-phase polymer ions may retain structural features associated with their electrospray formation mechanisms.

Organogel Formation Rationalized by Hansen Solubility Parameters: Shift of the Gelation Sphere with the Gelator Structure

To rationalize how the gelation ability of a low molecular weight gelator is influenced by its molecular structure, we performed extensive solubility tests of a group of thiazole-based gelators and made use of Hansen solubility parameter formalism. We observe that the increase of a linear alkyl chain in these gelators promotes an increase of the radius of the gelation sphere as well as a gradual shift of its center to lower values of the polar (δ_P) and hydrogen bonding (δ_H) components. The molecular packing within the fibers and the crystal habit were determined by a combination of X-ray diffraction and molecular modeling. We attribute the gradual and linear shift of the gelation sphere to the fact that all of the studied gelators share the same molecular packing, so that an increasing length of the alkyl chain reduces the proportion of polar groups at the surface, resulting in a gradual increase in the contact between apolar parts of the fiber and the solvent.