Cyclic Hexapeptides from Didemnum molle
Journal of Natural Products, 2008, Vol. 71, No. 6 945
solvent under a stream of N2, the residue was subjected to hydrolysis
and derivatization as described above with mollamide B. The D/L-
configuration for serine, proline, and leucine were assigned using the
two important NMR-derived distance restraints in case of MolB-S
(Figure 2). The distances between the R-proton of valine-1 and valine-2
were 5.95 and 6.86 Å, while those between the R-proton of valine-2
and the side chain methyl of threonine were 5.16 and 5.45 Å for MolB-R
and MolB-S, respectively. The former restraint violation appears to be
critical since the two atoms involved in the restraint flank the thiazoline
ring, and thus their positions (coordinates) are expected to be affected
by the chirality at the R-carbon of the ring. In other words, the optimized
conformation of MolB-R appeared to match the experimental distances
derived from the NMR data, while MolB-S showed two violations and
was also energetically less stable than the former. On the basis of these
observations, the chirality at the R-carbon of the thiazoline ring was
assigned to be R, corresponding to the L-amino acid configuration for
the thiazoline moiety.
C18 column, and L-isoleucine and L-allo-isoleucine were assigned using
the C5 column as described above.Then the ion at m/z 384 was extracted
for the four isoleucine isomers.
The analysis of mollamide B (1) using the C18 column established
the presence of L-Pro (34.47) [D-Pro (40.52)], L-Val (43.47), [D-Val
(54.10)], and L-Phe (52.41) [D-Phe (59.86)]. Amino acid standards using
LCTOF revealed retention times as follows: L-Thr (36.42 min) [L-allo-
Thr (34.94), D-Thr (45.46), D-allo-Thr (41.37)].
The analysis of mollamide C (2) using the C18 column established
the presence of L-Ser (21.56 min) [D-Ser (26.13)], L-Pro (34.07) [D-Pro
(40.25)], and L-Leu (51.82) [D-Leu (60.85)]. Amino acid standards using
LCTOF revealed retention times as follows: L-Ile (52.63) [L-allo-Ile
(52.93), D-Ile (70.42), D-allo-Ile (71.37)]. Due to the close retention
time of L-Ile and L-allo-Ile, mollamide C was co-injected with L-Ile
and L-allo-Ile standards, which clearly revealed that Ile has the
L-configuration. Thus, the analysis of mollamides B (1) and C (2)
established the presence of all L-amino acids.
Molecular Modeling. In order to investigate the chirality at the
R-carbon of the thiazoline moiety, two structures of molecule 1 were
drawn using the BUILDER module in Insight II (Accelrys Inc., San
Diego, CA). All regular amino acids were drawn as S, while the
R-carbon of the thiazoline moiety was drawn with both S and R
configurations to provide two initial configurations of molecule 1:
MolB-R and MolB-S (denoting the chirality at the R-carbon of the
thiazoline moiety to be L-Tzn and D-Tzn, respectively).
The structures were first optimized using a series of concurrent
minimizations employing the CVFF forcefield to a gradient of 0.001.
The optimized structures were then subjected to simulated annealing
using the DISCOVER (v 2.98) module of Insight II. All the simulations
were performed using a distance-dependent dielectric. A forcing
potential of 30 kcal/mol/rad2 was applied to the omega dihedrals for
all amino acids, except for the cis-proline, to maintain the trans
geometry. Also, NMR-derived distance restraints based on ROESY
correlations (Figure 1) were introduced during the MD simulations with
a force constant of 20 kcal/mol/rad2. The molecules were minimized
beginning with 300 steps of steepest descents and 1000 steps of
conjugate gradients to relieve any strain in the starting structures. The
structures were then “heated” gradually from 100 to 1000 K in steps
of 100 K with 1500 fs of simulation at each step. At 1000 K, the
simulation was further continued for 20 000 fs, during which 20
conformations were sampled at an interval of 1000 fs. The conforma-
tions were then gradually “cooled” to 300 K in steps of 100 K with
2500 fs simulation at each step. The structures were finally minimized
using 500 steps of steepest descents, 1500 steps of conjugate gradients,
and 500 steps of VA09A. At this stage, the NMR-derived distance
restraints were removed, and all 20 conformations were again subjected
to the same minimization protocol as given above.
The structures were analyzed, and the lowest energy conformation
of MolB-R obtained without any restraints appeared to maintain
important interproton distances in accordance with the NMR data. The
structures were then further optimized by semiempirical quantum
mechanical calculations. The AM1 Hamiltonian was used to calculate
the heat of formation (∆Hf) using the MOPAC suite of programs (v
6.0). The semiempirical calculations were carried out without any
restraints, and all degrees of freedom in the molecule were optimized.
Minimization was run with the EF algorithm to a GNORM of 0.01.
Since it is well known that most semiempirical methods underestimate
the barrier to rotation in peptides, the key word MMOK was used to
allow for a molecular mechanics correction. Stable points on the
potential energy surface (PES) of the molecule were identified by
running a FORCE calculation and were confirmed by the absence of
any imaginary (negative) frequencies.
Acknowledgment. We appreciate the assistance of F. T. Wiggers
and A. G. Shilabin in acquiring NMR and LC-TOF data. We are grateful
to the National Cancer Institute (NCI) and F. Valeriote at Wayne State
for performing cytotoxicity assays; The National Center for Natural
Products Research for performing the antimicrobial and antimalarial
assays; R. F. Schinazi’s group at Emory for HIV-1; S. G. Franzblau’s
group at UIC for Mtb; and A. M. S. Mayer for anti-inflammatory assays.
Financial support for this work was provided by grant numbers
1R01A136596 and P20 RR021929 from the NIAID and NCRR
components of the National Institutes of Health (NIH), and the contents
are solely the responsibility of the authors and do not necessarily
represent the official view of the NIH. The Egyptian Government and
the CDC are also very gratefully acknowledged for financial support.
This investigation was conducted in a facility constructed with support
from Research Facilities Improvements Program (C06 RR-14503-01)
from the National Center for Research Resources, NIH.
References and Notes
(1) Rinehart, K. L. Med. Res. ReV. 2000, 20, 1–27.
(2) Ireland, C.; Scheuer, P. J. J. Am. Chem. Soc. 1980, 102, 5688–5691.
(3) Wipf, P.; Fritch, P. C.; Geib, S. G.; Sefler, A. M. J. Am. Chem. Soc.
1998, 120, 4105–4112.
(4) Roy, R. S.; Gehring, A. M.; Milne, J. C.; Belshaw, P. J.; Walsh, C. T.
Nat. Prod. Rep. 1999, 16, 249–263.
(5) Carroll, A. R.; Bowden, B. F.; Coll, J. C.; Hockless, D. C. R.; Skelton,
B. W.; White, A. H. Aust. J. Chem. 1994, 47, 61–69.
(6) Toske, S. G.; Fenical, W. Tetrahedron Lett. 1995, 36, 8355–8358.
(7) Arrault, A.; Witczak-Legrand, A.; Gonzalez, P.; Banaigs, B. Tetra-
hedron Lett. 2002, 43, 4041–4044.
(8) Rudi, A.; Aknin, M.; Gaydou, E. M.; Kashman, Y. Tetrahedron 1998,
54, 13203–13210.
(9) Rudi, A.; Chill, L.; Aknin, M.; Kashman, Y. J. Nat. Prod. 2003, 66,
575–577.
(10) Mu¨ller, W. E. G.; Maidhof, A.; Zahn, R. K.; Conrad, J.; Rose, T.;
Uhlenbruck, G. Biol. Cell. 1984, 51, 381–388.
(11) Schmidt, E. W.; Nelson, J. T.; Rasko, D. A.; Sudek, S.; Eisen, J. A.;
Haygood, M. G.; Ravel, J. Proc. Natl. Acad. Sci. 2005, 102, 7315–
7320.
(12) Donia, M. S.; Hathaway, B. J.; Sudek, S.; Haygood, M. G.; Rosovitz,
M. J.; Ravel, J.; Schmidt, E. W. Nat. Chem. Biol. 2006, 2, 729–735.
(13) Wesson, K. J.; Hamann, M. T. J. Nat. Prod. 1996, 59, 629–631.
(14) Siemion, I. Z.; Wieland, T.; Pook, K. H. Angew. Chem., Int. Ed. Engl.
1975, 14, 702–703.
(15) Marfeys, P. Carlsberg Res. Commun. 1984, 49, 591.
(16) Yonetani, K.; Hirotsu, Y.; Shiba, T. Bull. Chem. Soc. Jpn. 1975, 48,
3302–3305.
(17) Szabo, S.; Szokan, G.; Khalafulla, A. M.; Almas, M.; Kiss, C.; Rill,
A.; Schon, I. J. Pept. Sci. 2001, 7, 316–322.
(18) Sharman, G. J.; Jones, I. C. Magn. Reson. Chem. 2001, 39, 549–554.
(19) Costantine, K. L.; Mueller, L.; Huang, S.; Abid, S.; Lam, K. S.; Li,
W.; Leet, J. E. J. Am. Chem. Soc. 2002, 124, 7284–7285.
(20) Leusch, H.; Yoshida, W. Y.; Moore, R. E.; Paul, V. J.; Corbett, T. H.
J. Am. Chem. Soc. 2001, 123, 5418–5423.
The ∆Hf values were calculated to be -175.2 kcal/mol for MolB-R
and -166.7 kcal/mol for MolB-S. Thus, MolB-R appeared to be more
stable by 8.5 kcal/mol. Also, significant violations were observed for
(21) McDonald, L. A.; Ireland, C. M. J. Nat. Prod. 1992, 55, 376–379.
NP700718P