1
|
Naleem N, Abreu CRA, Warmuz K, Tong M, Kirmizialtin S, Tuckerman ME. An exploration of machine learning models for the determination of reaction coordinates associated with conformational transitions. J Chem Phys 2023; 159:034102. [PMID: 37458344 DOI: 10.1063/5.0147597] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/23/2023] [Accepted: 06/23/2023] [Indexed: 07/20/2023] Open
Abstract
Determining collective variables (CVs) for conformational transitions is crucial to understanding their dynamics and targeting them in enhanced sampling simulations. Often, CVs are proposed based on intuition or prior knowledge of a system. However, the problem of systematically determining a proper reaction coordinate (RC) for a specific process in terms of a set of putative CVs can be achieved using committor analysis (CA). Identifying essential degrees of freedom that govern such transitions using CA remains elusive because of the high dimensionality of the conformational space. Various schemes exist to leverage the power of machine learning (ML) to extract an RC from CA. Here, we extend these studies and compare the ability of 17 different ML schemes to identify accurate RCs associated with conformational transitions. We tested these methods on an alanine dipeptide in vacuum and on a sarcosine dipeptoid in an implicit solvent. Our comparison revealed that the light gradient boosting machine method outperforms other methods. In order to extract key features from the models, we employed Shapley Additive exPlanations analysis and compared its interpretation with the "feature importance" approach. For the alanine dipeptide, our methodology identifies ϕ and θ dihedrals as essential degrees of freedom in the C7ax to C7eq transition. For the sarcosine dipeptoid system, the dihedrals ψ and ω are the most important for the cisαD to transαD transition. We further argue that analysis of the full dynamical pathway, and not just endpoint states, is essential for identifying key degrees of freedom governing transitions.
Collapse
Affiliation(s)
- Nawavi Naleem
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
| | - Charlles R A Abreu
- Chemical Engineering Department, Escola de Química, Universidade Federal do Rio de Janeiro, 21941-909 Rio de Janeiro, RJ, Brazil
| | - Krzysztof Warmuz
- Computer Science Program, Science Division, New York University, Abu Dhabi, UAE
| | - Muchen Tong
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
| | - Serdal Kirmizialtin
- Chemistry Program, Science Division, New York University, Abu Dhabi, UAE
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Center for Smart Engineering Materials, New York University, Abu Dhabi, UAE
| | - Mark E Tuckerman
- Department of Chemistry, New York University (NYU), New York, New York 10003, USA
- Courant Institute of Mathematical Sciences, New York University, New York, New York 10012, USA
- NYU-ECNU Center for Computational Chemistry at NYU Shanghai, 3663 Zhongshan Rd. North, Shanghai 200062, China
- Simons Center for Computational Physical Chemistry at New York University, New York, New York 10003, USA
| |
Collapse
|
2
|
Crapster JA, Guzei IA, Blackwell HE. A peptoid ribbon secondary structure. Angew Chem Int Ed Engl 2013; 52:5079-84. [PMID: 23576308 DOI: 10.1002/anie.201208630] [Citation(s) in RCA: 89] [Impact Index Per Article: 8.1] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/26/2012] [Revised: 01/02/2013] [Indexed: 01/02/2023]
Affiliation(s)
- J Aaron Crapster
- Department of Chemistry, University of Wisconsin-Madison, 1101 University Avenue, Madison, WI 53706-1322, USA
| | | | | |
Collapse
|
3
|
Park SH, Szleifer I. Structural and dynamical characteristics of peptoid oligomers with achiral aliphatic side chains studied by molecular dynamics simulation. J Phys Chem B 2011; 115:10967-75. [PMID: 21819112 PMCID: PMC3177600 DOI: 10.1021/jp2025957] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
Abstract
All-atom molecular dynamics simulations of N-substituted glycine peptoid oligomers with methyl and methoxyethyl side chains have been carried out for chain lengths of 5, 10, 20, and 50 residues in aqueous phase at room temperature. The (ϕ, ψ) backbone dihedral angle distributions in the Ramachandran plots show that helical structures, similar to polyproline type I and type II helices, are the most favorable conformations in most peptoid oligomers studied. The left-handed helical structures are shown to be increasingly favored as the oligomer chain length grows. A significant population of cis amide bond configurations has been identified in the peptoid oligomers. By combining the analysis of ϕ and ω backbone dihedral angles, we determined the relative composition of the four major conformations favored by the backbone dihedral angles. The trans α(D) conformation is found to be most favored for all peptoid oligomers studies. The time correlation functions of the end-to-end distance highlight a rigid backbone structure relative to side chains for peptoid oligomers. The transition between right-handed and left-handed helical conformations is found to be very rare and between cis and trans isomerism in the amide bond completely absent in the simulation time scale. The radii of gyration for all peptoid oligomers have been found to be consistently larger in comparison to the peptide counterparts, suggesting slightly open structures for peptoids relative to peptides, whereas the fluctuations in the radius of gyration support a rigid backbone structure of peptoids.
Collapse
Affiliation(s)
- Sung Hyun Park
- Department of Biomedical Engineering and Chemistry of Life Processes, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, U.S.A
| | - Igal Szleifer
- Department of Biomedical Engineering and Chemistry of Life Processes, Northwestern University, 2145 Sheridan Road, Evanston, Illinois 60208, U.S.A
| |
Collapse
|
4
|
Crapster JA, Stringer JR, Guzei IA, Blackwell HE. Design and conformational analysis of peptoids containing N-hydroxy amides reveals a unique sheet-like secondary structure. Biopolymers 2011; 96:604-16. [PMID: 22180908 PMCID: PMC3448284 DOI: 10.1002/bip.21599] [Citation(s) in RCA: 44] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/30/2022]
Abstract
N-hydroxy amides can be found in many naturally occurring and synthetic compounds and are known to act as both strong proton donors and chelators of metal cations. We have initiated studies of peptoids, or N-substituted glycines which contain N-hydroxy amide side chains to investigate the potential effects of these functional groups on peptoid backbone amide rotamer equilibria and local conformations. We reasoned that the propensity of these functional groups to participate in hydrogen bonding could be exploited to enforce intramolecular or intermolecular interactions that yield new peptoid structures. Here, we report the design, synthesis, and detailed conformational analysis of a series of model N-hydroxy peptoids. These peptoids were readily synthesized, and their structures were analyzed in solution by 1D and 2D NMR and in the solid-state by X-ray crystallography. The N-hydroxy amides were found to strongly favor trans conformations with respect to the peptoid backbone in chloroform. More notably, unique sheet-like structures held together via intermolecular hydrogen bonds were observed in the X-ray crystal structures of an N-hydroxy amide peptoid dimer, which to our knowledge represent the first structure of this type reported for peptoids. These results suggest that the N-hydroxy amide can be utilized to control both local backbone geometries and longer-range intermolecular interactions in peptoids, and represents a new functional group in the peptoid design toolbox.
Collapse
Affiliation(s)
- J. Aaron Crapster
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, Wisconsin 53706-1322
| | - Joseph R. Stringer
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, Wisconsin 53706-1322
| | - Ilia A. Guzei
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, Wisconsin 53706-1322
| | - Helen E. Blackwell
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, Wisconsin 53706-1322
| |
Collapse
|
5
|
Stringer JR, Crapster JA, Guzei IA, Blackwell HE. Construction of peptoids with all trans-amide backbones and peptoid reverse turns via the tactical incorporation of N-aryl side chains capable of hydrogen bonding. J Org Chem 2010; 75:6068-78. [PMID: 20722367 PMCID: PMC3175427 DOI: 10.1021/jo101075a] [Citation(s) in RCA: 74] [Impact Index Per Article: 5.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/30/2022]
Abstract
The ability to design foldamers that mimic the defined structural motifs of natural biopolymers is critical for the continued development of functional biomimetic molecules. Peptoids, or oligomers of N-substituted glycine, represent a versatile class of foldamers capable of folding into defined secondary and tertiary structures. However, the rational design of discretely folded polypeptoids remains a challenging task, due in part to an incomplete understanding of the covalent and noncovalent interactions that direct local peptoid folding. We have found that simple, peptoid monomer model systems allow for the effective isolation of individual interactions within the peptoid backbone and side chains and can facilitate the study of the role of these interactions in restricting local peptoid conformation. Herein, we present an analysis of a set of peptoid monomers and an oligomer containing N-aryl side chains capable of hydrogen bonding with the peptoid backbone. These model peptoids were found to exhibit well-defined local conformational preferences, allowing for control of the ω, ϕ, and ψ dihedral angles adopted by the systems. Fundamental studies of the peptoid monomers enabled the design and synthesis of an acyclic peptoid reverse-turn structure, in which N-aryl side chains outfitted with ortho-hydrogen bond donors were hypothesized to play a critical role in the stabilization of the turn. This trimeric peptoid was characterized by X-ray crystallography and 2D NMR spectroscopy and was shown to adopt a unique acyclic peptoid reverse-turn conformation. Further analysis of this turn revealed an n→π*(C═O) interaction within the peptoid backbone, which represents the first reported example of this type of stereoelectronic interaction occurring exclusively within a polypeptoid backbone. The installation of N-aryl side chains capable of hydrogen bonding into peptoids is straightforward and entirely compatible with current solid-phase peptoid synthesis methodologies. As such, we anticipate that the strategic incorporation of these N-aryl side chains should facilitate the construction of peptoids capable of adopting discrete structural motifs, both turnlike and beyond, and will facilitate the continued development of well-folded peptoids.
Collapse
Affiliation(s)
- Joseph R. Stringer
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, WI 53706-1322
| | - J. Aaron Crapster
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, WI 53706-1322
| | - Ilia A. Guzei
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, WI 53706-1322
| | - Helen E. Blackwell
- Department of Chemistry, University of Wisconsin–Madison, 1101 University Avenue, Madison, WI 53706-1322
| |
Collapse
|
6
|
Armand P, Kirshenbaum K, Goldsmith RA, Farr-Jones S, Barron AE, Truong KT, Dill KA, Mierke DF, Cohen FE, Zuckermann RN, Bradley EK. NMR determination of the major solution conformation of a peptoid pentamer with chiral side chains. Proc Natl Acad Sci U S A 1998; 95:4309-14. [PMID: 9539733 PMCID: PMC22485 DOI: 10.1073/pnas.95.8.4309] [Citation(s) in RCA: 268] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
Polymers of N-substituted glycines ("peptoids") containing chiral centers at the alpha position of their side chains can form stable structures in solution. We studied a prototypical peptoid, consisting of five para-substituted (S)-N-(1-phenylethyl)glycine residues, by NMR spectroscopy. Multiple configurational isomers were observed, but because of extensive signal overlap, only the major isomer containing all cis-amide bonds was examined in detail. The NMR data for this molecule, in conjunction with previous CD spectroscopic results, indicate that the major species in methanol is a right-handed helix with cis-amide bonds. The periodicity of the helix is three residues per turn, with a pitch of approximately 6 A. This conformation is similar to that anticipated by computational studies of a chiral peptoid octamer. The helical repeat orients the amide bond chromophores in a manner consistent with the intensity of the CD signal exhibited by this molecule. Many other chiral polypeptoids have similar CD spectra, suggesting that a whole family of peptoids containing chiral side chains is capable of adopting this secondary structure motif. Taken together, our experimental and theoretical studies of the structural properties of chiral peptoids lay the groundwork for the rational design of more complex polypeptoid molecules, with a variety of applications, ranging from nanostructures to nonviral gene delivery systems.
Collapse
Affiliation(s)
- P Armand
- Chiron Technologies, Chiron Corporation, 4560 Horton Street, Emeryville, CA 94608, USA
| | | | | | | | | | | | | | | | | | | | | |
Collapse
|
7
|
Kirshenbaum K, Barron AE, Goldsmith RA, Armand P, Bradley EK, Truong KT, Dill KA, Cohen FE, Zuckermann RN. Sequence-specific polypeptoids: a diverse family of heteropolymers with stable secondary structure. Proc Natl Acad Sci U S A 1998; 95:4303-8. [PMID: 9539732 PMCID: PMC22484 DOI: 10.1073/pnas.95.8.4303] [Citation(s) in RCA: 379] [Impact Index Per Article: 14.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/07/2023] Open
Abstract
We have synthesized and characterized a family of structured oligo-N-substituted-glycines (peptoids) up to 36 residues in length by using an efficient solid-phase protocol to incorporate chemically diverse side chains in a sequence-specific fashion. We investigated polypeptoids containing side chains with a chiral center adjacent to the main chain nitrogen. Some of these sequences have stable secondary structure, despite the achirality of the polymer backbone and its lack of hydrogen bond donors. In both aqueous and organic solvents, peptoid oligomers as short as five residues give rise to CD spectra that strongly resemble those of peptide alpha-helices. Differential scanning calorimetry and CD measurements show that polypeptoid secondary structure is highly stable and that unfolding is reversible and cooperative. Thermodynamic parameters obtained for unfolding are similar to those obtained for the alpha-helix to coil transitions of peptides. This class of biomimetic polymers may enable the design of self-assembling macromolecules with novel structures and functions.
Collapse
Affiliation(s)
- K Kirshenbaum
- Department of Pharmaceutical Chemistry, University of California, San Francisco, CA 94143
| | | | | | | | | | | | | | | | | |
Collapse
|