Macrocyclization of Peptide Side Chains by Ugi Reaction: Achieving Peptide Folding and Exocyclic N-Functionalization in One Shot Aldrin V. Vasco,, Carlos S. Pérez, Fidel E. Morales, Hilda E. Garay, ǁ Dimitar Vasilev, José A. Gavín, ǂ Ludger A. Wessjohann,, * and Daniel G. Rivera,, * Center for Natural Products Research, Faculty of Chemistry, University of Havana, Zapata y G, 10400, La Habana, Cuba. Facultad de Ingeniería Química. Instituto Superior Politécnico José Antonio Echeverría, CUJAE, Calle 114 # 11901, 11500, La Habana, Cuba. ǁ Synthetic Peptides Group, Center for Genetic Engineering and Biotechnology, P.O. Box 6162, La Habana, Cuba. ǂ Instituto Universitario de Bioorgánica Antonio González and Departamento de Química Orgánica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain. Department of Bioorganic Chemistry, Leibniz Institute of Plant Biochemistry, Weinberg 3, D- 06120, Halle/Saale, Germany. Corresponding Authors: *Emails: dgr@fq.uh.cu; wessjohann@ipb-halle.de S1
Table of contents NMR Structure determination Molecular Dynamics Simulations Figure 1. HPLC chromatogram of crude peptide 3. Figure 2. ESI-MS of peptide 3. Figure 3. HPLC Chromatogram of crude peptide 5. Figure 4. ESI-MS of peptide 5 Figure 5. HPLC Chromatogram of pure cyclic peptide 2. Figure 6. ESI-MS of cyclic peptide 2. Figure 7. HPLC Chromatogram of pure cyclic peptide 4. Figure 8. Mass spectrum of cyclic peptide 4. Figure 9. HPLC Chromatogram of pure cyclic peptide 6. S3 S4 S5 S5 S6 S6 S7 S7 S8 S8 S9 Figure 10. Mass spectrum of cyclic peptide 6. S9 Figure 11. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 2. S10 Figure 12. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 2. S11 Figure 13. TOCSY spectrum in DMSO-d6 of cyclic peptide 2. S12 Figure 14. HSQC spectrum of in DMSO-d6 cyclic peptide 2. S13 Figure 15. HMBC spectrum in DMSO-d6 of cyclic peptide 2. S14 Figure 16. ROESY spectrum in DMSO-d6 of cyclic peptide 2. S15 Figure 17. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 4. S17 Figure 18. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 4. S18 Figure 19. TOCSY spectrum in DMSO-d6 of cyclic peptide 4. S19 Figure 20. HSQC spectrum in DMSO-d6 of cyclic peptide 4. S20 S2
Figure 21. HMBC spectrum in DMSO-d6 of cyclic peptide 4. S21 Figure 22. ROESY spectrum in DMSO-d6 of cyclic peptide 4. S22 Figure 23. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 6. S23 Figure 24. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 6. S24 Figure 25. TOCSY spectrum of in DMSO-d6 cyclic peptide 6. S25 Figure 26. HSQC spectrum in DMSO-d6 of cyclic peptide 6. S26 Figure 27. HMBC spectrum in DMSO-d6 of cyclic peptide 6. S27 Figure 28. ROESY spectrum in DMSO-d6 of cyclic peptide 6. S28 Figure 29. 600 MHz 1 H-NMR spectrum in acetone-d6 of cyclic peptide 8. S29 Figure 30. 150 MHz 13 C-NMR spectrum acetone-d6 of cyclic peptide 8. S30 Figure 31. 600 MHz 1 H-NMR spectrum in CDCl 3 of cyclic peptide 10. S31 Figure 32. 150 MHz 13 C NMR spectrum in CDCl 3 of cyclic peptide 10. S32 References S33 NMR Structure determination NMR structure determination was performed in Xplor-NIH 1,2 through simulated annealing regularization and refinement in torsion angle space, using experimental data as inter-proton distances and dihedral angles restraints. Simulations were based on examples scripts distributed in the /eginput folder within Xplor-NIH package. For cyclic peptide 2, φ angles from the 3 J NHCHα of residues Leu1, Lys2, Phe3, Val4 and Glu5 as well as χ 1 angle from the Phe3 residue were used as dihedral constraints. For cyclic peptide 4, φ angles from the 3 J NHCHα of residues Val1, Ile2, Lys6 and Phe7 as well as χ 1 angle from the Phe7 residue were used as dihedral constraints. For cyclic peptide 6, φ angles from the 3 J NHCHα of residues Phe1 and Lys2 were used as dihedral constraints. For simulated annealing regularization 200 starting structures were randomly generated. Geometric bond, angle and improper energy contributions were scaled to 1.0, 0.4, and 0.1 respectively. Non-bonded Van der Waals contributions were set to 0.004 and experimental NOE S3
and dihedral constraints were scaled to 2 and 10 respectively. An initial 100 steps conjugated gradient Powell minimization was performed. A 100 ps molecular dynamics simulation at 3500 K was performed with a timestep of 3 fs. Before the cooling stage, the energy function due to experimental dihedral restraints was scaled to 200. The system was cooled from 3500 to 25 K, with a temperature step of 12.5 K. At each temperature step, 0.2 ps of molecular dynamics simulation was performed. During this stage the Van der Waals energy term was progressively scaled from 0.004 to 4, and the experimental NOE potential was climbed to 30. Angle and improper terms were both progressively scaled to 1.0. A 500 steps torsion angle minimization was performed and finally the system was optimized by means of 500 steps conjugated gradient Powell cartesian minimization. The refinement protocol consisted in a slow cooling simulated annealing from the regularized structures. The initial weights of the energy functions was 0.4, 0.1 and 1.0 for the angle, improper, and bond terms; 10 and 2 respectively for the NOE and dihedral experimental restraints; and 0.004 for the non-bonded Van der Waals term. A 10 ps molecular dynamics simulation at 3000 K was achieved with a time-step of 3 fs, afterwards dihedral restraints contribution was set to 200. The system was cooled with a temperature step of 12.5 K and a simulation time of 0.2 ps at each temperature. During this simulation the Van der Waals energy term was scaled from 0.003 to 4. During this stage the Van der Waals energy term was progressively scaled from 0.004 to 4, and the experimental NOE potential was climbed to 30. Angle and improper terms were both progressively scaled to 1.0. A 500 torsion angle minimization was performed afterwards a second 500 steps minimization was achieved in cartesian coordinates. A finally 1000 steps Powell minimization with an energy function nondependent of experimental restraints was executed. Final RMSD were calculated with the plugin RMSD calculator within VMD 1.9.2, 3 excluding N- and C-terminal residues in the calculation. Molecular Dynamics Simulations Solvent boxes were generated with the Solvate plugin within VMD 1.9.2. 3 Trajectories analysis of interatomic distances were performed in VEGA ZZ 4 while time distribution of ϕ and ψ angles (Ramachandran Plots) were generated within VMD. 3 The analysis of hydrogen bonds in time was performed within VMD using a cut-off criteria of 20 0 and a donor-acceptor distance criteria lesser than 3Å. S4
Figure 1. HPLC chromatogram of crude peptide 3. fidel-1 5 (0.108) AM (Cen,2, 80.00, Ht,5000.0,0.00,1.00); Sm (Mn, 2x3.00); Cm (1:53) 466.77 100 TOF MS ES+ 2.00e4 467.27 % 458.26 678.37 549.33 932.54 429.75 455.27 430.25 467.78 550.34 712.43 479.30 583.33 790.39 792.46 415.75 430.76 680.38 495.27 551.34 860.50 934.56 532.30 584.33 610.86 645.32 661.36 728.45 733.85 774.44 780.87 799.38 820.87 840.40 874.48 903.51 914.33 935.55 0 m/z 400 420 440 460 480 500 520 540 560 580 600 620 640 660 680 700 720 740 760 780 800 820 840 860 880 900 920 940 679.38 711.43 791.45 Figure 2. ESI-MS of peptide 3. 858.49 859.50 933.55 S5
Figure 3. HPLC Chromatogram of crude peptide 5. fidel-2 85 (1.643) AM (Cen,2, 80.00, Ht,5000.0,0.00,1.00); Sm (Mn, 2x3.00); Cm (2:86) 746.41 100 % 747.42 0 730.04 429.12 445.16 557.33 640.32 658.35 748.43 768.41 817.45898.35 1010.67 1119.571138.59 1213.00 1244.22 Figure 4. ESI-MS of peptide 5 S6
Figure 5. HPLC Chromatogram of pure cyclic peptide 2. Figure 6. ESI-MS of cyclic peptide 2. S7
Figure 7. HPLC Chromatogram of pure cyclic peptide 4. favv361-pico2 17 (0.326) AM (Cen,2, 80.00, Ht,5000.0,0.00,1.00); Sm (Mn, 2x3.00); Cm (2:23) A 100 1027.59 % A2 514.30 1028.59 505.79 514.80 0 440.73 515.30 533.27 665.40 733.39 806.48 807.48 880.47 939.03 953.55 1049.58 1050.59 1066.56 1179.57 Figure 8. Mass spectrum of cyclic peptide 4. S8
Figure 9. HPLC Chromatogram of pure cyclic peptide 6. avv037 110 (2.100) AM (Cen,2, 80.00, Ht,5000.0,0.00,1.00); Sm (Mn, 2x3.00); Cm (2:119) 100 885.44 % 886.45 0 443.23 462.20 532.29562.27 608.31 652.32 696.37 887.45 717.36 797.39 769.39 826.45 868.43 908.43 989.50 Figure 10. Mass spectrum of cyclic peptide 6. S9
Figure 11. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 2. S10
Figure 12. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 2. S11
Figure 13. TOCSY spectrum in DMSO-d6 of cyclic peptide 2. S12
Figure 14. HSQC spectrum of in DMSO-d6 cyclic peptide 2. S13
Figure 15. HMBC spectrum in DMSO-d6 of cyclic peptide 2. S14
Figure 16. ROESY spectrum in DMSO-d6 of cyclic peptide 2. S15
S16
Figure 17. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 4. S17
Figure 18. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 4. S18
Figure 19. TOCSY spectrum in DMSO-d6 of cyclic peptide 4. S19
Figure 20. HSQC spectrum in DMSO-d6 of cyclic peptide 4. S20
Figure 21. HMBC spectrum in DMSO-d6 of cyclic peptide 4. S21
Figure 22. ROESY spectrum in DMSO-d6 of cyclic peptide 4. S22
Figure 23. 600 MHz 1 H-NMR spectrum in DMSO-d6 of cyclic peptide 6. S23
Figure 24. 150 MHz 13 C-NMR spectrum in DMSO-d6 of cyclic peptide 6. S24
Figure 25. TOCSY spectrum of in DMSO-d6 cyclic peptide 6. S25
Figure 26. HSQC spectrum in DMSO-d6 of cyclic peptide 6. S26
Figure 27. HMBC spectrum in DMSO-d6 of cyclic peptide 6. S27
Figure 28. ROESY spectrum in DMSO-d6 of cyclic peptide 6. S28
Figure 29. 600 MHz 1 H-NMR spectrum in acetone-d6 of cyclic peptide 8. S29
Figure 30. 150 MHz 13 C-NMR spectrum acetone-d6 of cyclic peptide 8. S30
Figure 31. 600 MHz 1 H-NMR spectrum in CDCl 3 of cyclic peptide 10. S31
Figure 32. 150 MHz 13 C NMR spectrum in CDCl 3 of cyclic peptide 10. S32
References (1) Schwieters, C. D.; Kuszewski, J. J.; Clore, G. M. Prog. Nucl. Magn. Reson. Spectrosc. 2006, 48, 47 (2) Schwieters, C. D.; Kuszewski, J. J.; Tjandra, N.; G. M. Clore J. Magn. Reson. 2003, 160, 65. (3) Humphrey, W.; Dalke, A.; Schulten, K. J. Molec. Graphics 1996, 14, 33. (4) Pedretti, A.; Villa, L.; Vistoli, G. J. Mol. Graph. 2002, 21, 47. S33