skip to main content

Title: Rational Design of an Organocatalyst for Peptide Bond Formation
Amide bonds are ubiquitous in peptides, proteins, pharmaceuticals, and polymers. The formation of amide bonds is a straightforward process: amide bonds can be synthesized with relative ease because of the availability of efficient coupling agents. However, there is a substantive need for methods that do not require excess reagents. A catalyst that condenses amino acids could have an important impact by reducing the significant waste generated during peptide synthesis. We describe the rational design of a biomimetic catalyst that can efficiently couple amino acids featuring standard protecting groups. The catalyst design combines lessons learned from enzymes, peptide biosynthesis, and organocatalysts. Under optimized conditions, 5 mol % catalyst efficiently couples Fmoc amino acids without notable racemization. Importantly, we demonstrate that the catalyst is functional for the synthesis of oligopeptides on solid phase. This result is significant because it illustrates the potential of the catalyst to function on a substrate with a multitude of amide bonds, which may be expected to inhibit a hydrogen-bonding catalyst.
Authors:
; ; ;
Award ID(s):
1807670
Publication Date:
NSF-PAR ID:
10180597
Journal Name:
Journal of the American Chemical Society
Volume:
141
Issue:
40
Page Range or eLocation-ID:
15977–15985
ISSN:
0002-7863
Sponsoring Org:
National Science Foundation
More Like this
  1. Disordered proline-rich motifs are common across the proteomes of many species and are often involved in protein-protein interactions. Proline is a unique amino acid due to the covalent bond between the backbone nitrogen and the proline side chain. The resulting five-membered ring allows proline to sample the cis state about its peptide bond, which other residues cannot do as readily. Because proline-rich disordered sequences exist as ensembles that likely include structures with the proline peptide bond in cis , a robust methodology to accurately account for these conformations in the overall ensemble is crucial. Observing the cis conformations of prolinemore »in a disordered sequence is challenging both experimentally and computationally. Nitrogen-hydrogen NMR spectroscopy cannot directly observe proline residues, which lack an amide bond, and computational methods struggle to overcome the large kinetic barrier between the cis and trans states, since isomerization usually occurs on the order of seconds. In the current work, Gaussian accelerated molecular dynamics was used to overcome this free energy barrier and simulate proline isomerization in a tetrapeptide (KPTP) and in the 12-residue proline-rich SH3 binding peptide, ArkA. We found that Gaussian accelerated molecular dynamics, when combined with a lowered peptide bond dihedral angle potential energy barrier (15 kcal/mol), allowed sufficient sampling of the proline cis and trans states on a microsecond timescale. All ArkA prolines spend a significant fraction of time in cis , leading to a more compact ensemble with less polyproline II helix structure than an ArkA ensemble with all peptide bonds in trans . The ensemble containing cis prolines also matches more closely to in vitro circular dichroism data than the all- trans ensemble. The ability of the ArkA prolines to isomerize likely affects the peptide’s ability to bind its partner SH3 domain, and should be studied further. This is the first molecular dynamics simulation study of proline isomerization in a biologically relevant proline-rich sequence that we know of, and a similar protocol could be applied to study multi-proline isomerization in other proline-containing proteins to improve conformational diversity and agreement with in vitro data.« less
  2. Abstract Short hydrogen bonds (SHBs), whose donor and acceptor heteroatoms lie within 2.7 Å, exhibit prominent quantum mechanical characters and are connected to a wide range of essential biomolecular processes. However, exact determination of the geometry and functional roles of SHBs requires a protein to be at atomic resolution. In this work, we analyze 1260 high-resolution peptide and protein structures from the Protein Data Bank and develop a boosting based machine learning model to predict the formation of SHBs between amino acids. This model, which we name as machine learning assisted prediction of short hydrogen bonds (MAPSHB), takes into accountmore »21 structural, chemical and sequence features and their interaction effects and effectively categorizes each hydrogen bond in a protein to a short or normal hydrogen bond. The MAPSHB model reveals that the type of the donor amino acid plays a major role in determining the class of a hydrogen bond and that the side chain Tyr-Asp pair demonstrates a significant probability of forming a SHB. Combining electronic structure calculations and energy decomposition analysis, we elucidate how the interplay of competing intermolecular interactions stabilizes the Tyr-Asp SHBs more than other commonly observed combinations of amino acid side chains. The MAPSHB model, which is freely available on our web server, allows one to accurately and efficiently predict the presence of SHBs given a protein structure with moderate or low resolution and will facilitate the experimental and computational refinement of protein structures.« less
  3. Peptide backbone amide substitution can dramatically alter the conformational and physiochemical properties of native sequences. Although uncommon relative to N -alkyl substituents, peptides harboring main-chain N -hydroxy groups exhibit unique conformational preferences and biological activities. Here, we describe a versatile method to prepare N -hydroxy peptide on solid support and evaluate the impact of backbone N -hydroxylation on secondary structure stability. Based on previous work demonstrating the β-sheet-stabilizing effect of α-hydrazino acids, we carried out an analogous study with N -hydroxy-α-amino acids using a model β-hairpin fold. In contrast to N -methyl substituents, backbone N -hydroxy groups are accommodated inmore »the β-strand region of the hairpin without energetic penalty. An enhancement in β-hairpin stability was observed for a di- N -hydroxylated variant. Our results facilitate access to this class of peptide derivatives and inform the use of backbone N -hydroxylation as a tool in the design of constrained peptidomimetics.« less
  4. The chemical modification of peptides is a promising approach for the design of protein-protein interaction inhibitors and peptide-based drug candidates. Among several peptidomimetic strategies, substitution of the amide backbone maintains side-chain functionality that may be important for engagement of biological targets. Backbone amide substitution has been largely limited to N-alkylation, which can promote cis amide geometry and disrupt important H-bonding interactions. In contrast, N-amination of peptides induces distinct backbone geometries and maintains H-bond donor capacity. In this chapter we discuss the conformational characteristics of designed N-amino peptides and present a detailed protocol for their synthesis on solid support. The describedmore »methods allow for backbone N-amino scanning of biologically active parent sequences.« less
  5. Abstract
    <p>This data set for the manuscript entitled &#34;Design of Peptides that Fold and Self-Assemble on Graphite&#34; includes all files needed to run and analyze the simulations described in the this manuscript in the molecular dynamics software NAMD, as well as the output of the simulations. The files are organized into directories corresponding to the figures of the main text and supporting information. They include molecular model structure files (NAMD psf or Amber prmtop format), force field parameter files (in CHARMM format), initial atomic coordinates (pdb format), NAMD configuration files, Colvars configuration files, NAMD log files, and NAMD output includingMore>>