Structural and mutational analysis of the nonribosomal peptide synthetase heterocyclization domain provides insight into catalysis

Kristjan Bloudoff, Christopher D. Fage, Mohamed A. Marahiel, and T. Martin Schmeing
PNAS January 3, 2017 114 (1) 95-100; published ahead of print December 19, 2016 https://doi.org/10.1073/pnas.1614191114

  1. Edited by Peter B. Moore, Yale University, New Haven, CT, and approved November 30, 2016 (received for review August 31, 2016)

Significance

Nonribosomal peptide synthetases produce peptides with wide varieties of therapeutic and biological activities. Monomer substrates are typically linked by a condensation domain. However, in many modules, a heterocyclization (Cy) domain takes its place and performs both condensation and cyclodehydration of a cysteine, serine, or threonine to form a five-membered ring in the peptide backbone. Although studied for decades, the mechanisms of condensation and cyclodehydration by Cy domains were previously unknown. The crystal structure of a Cy domain, and accompanying mutagenic and bioinformatics analyses, uncover the importance of an aspartate and a threonine for the cyclodehydration reaction. This study provides insight into the catalysis of condensation by the Cy domain and enables the proposal of a reaction mechanism for cyclodehydration.

Abstract

Nonribosomal peptide synthetases (NRPSs) are a family of multidomain, multimodule enzymes that synthesize structurally and functionally diverse peptides, many of which are of great therapeutic or commercial value. The central chemical step of peptide synthesis is amide bond formation, which is typically catalyzed by the condensation (C) domain. In many NRPS modules, the C domain is replaced by the heterocyclization (Cy) domain, a homologous domain that performs two consecutive reactions by using hitherto unknown catalytic mechanisms. It first catalyzes amide bond formation, and then the intramolecular cyclodehydration between a Cys, Ser, or Thr side chain and the backbone carbonyl carbon to form a thiazoline, oxazoline, or methyloxazoline ring. The rings are important for the form and function of the peptide product. We present the crystal structure of an NRPS Cy domain, Cy2 of bacillamide synthetase, at a resolution of 2.3 Å. Despite sharing the same fold, the active sites of C and Cy domains have important differences. The structure allowed us to probe the roles of active-site residues by using mutational analyses in a peptide synthesis assay with intact bacillamide synthetase. The drastically different effects of these mutants, interpreted by using our structural and bioinformatic results, provide insight into the catalytic mechanisms of the Cy domain and implicate a previously unexamined Asp-Thr dyad in catalysis of the cyclodehydration reaction.

Nonribosomal peptide synthetases (NRPSs) are a family of large, multimodular enzymes that produce a range of important bioactive secondary metabolites (1). NRPS products have great diversity because they can use more than 500 different acyl monomer substrates, including l-amino acids, d-amino acids, aryl acids, fatty acids, hydroxy acids, and keto acids, and they can subsequently modify these moieties during peptide synthesis.

NRPSs function in a modular, assembly-line fashion. A typical elongation module consists of a condensation (C), an adenylation (A), and a peptidyl-carrier protein (PCP) domain. The A domain specifically recognizes and activates a monomer acyl substrate through adenylation, then transfers it onto the prosthetic phosphopantetheine arm of the PCP domain. This acyl-PCP then travels to the C domain acceptor site for condensation with the upstream module’s acyl-PCP at the C domain donor site (2⇓⇓–5). The PCP domain brings the elongated peptide chain to the downstream module, where it is passed off and further elongated in the next condensation reaction. This process is repeated in each module until the growing peptide reaches the termination module, where it is elongated and then released from the NRPS, often by a thioesterase domain. However, most NRPSs, along with their C, A, and PCP domains, also include tailoring and/or alternative domains, which cosynthetically modify the nonribosomal peptide.

One important modification that can occur during peptide synthesis is cyclodehydration of Cys, Ser, or Thr residues into thiazoline, oxazoline, or methyloxazoline rings, respectively, by the NRPS heterocyclization (Cy) domain (6⇓⇓⇓⇓⇓–12). These heterocyclic rings are found in many peptides with important clinical and research utility, such as the antibiotics bacitracin A (6) and zelkovamycin (13), the antitumor agents bleomycin (14) and epothilone (8), the immunosuppressant argyrin (15), the siderophores mycobactin (16) and yersiniabactin (17), and the microbiome genotoxin colibactin (18). Introduction of the five-membered heterocyclic ring makes the peptide resistant to proteolytic cleavage and can induce conformations in the peptide that favor interaction with biological targets (19).

In NRPSs that synthesize these heterocycle-containing products, the module specific for the Cys, Ser, or Thr monomer substrate contains a Cy domain in place of the C domain. Cy domains are evolutionarily and structurally related to C domains (20). The Cy domain first catalyzes nucleophilic attack on the thioester of a PCP-linked donor substrate by the α-amino group of a Cys-, Ser-, or Thr-PCP substrate, presumably in a manner similar to C domains (6, 7, 10⇓–12, 21⇓–23) (Fig. 1). In the two-step cyclodehydration reaction that follows, the thiol of the Cys side chain or hydroxyl of the Ser or Thr side chain first attacks the carbonyl carbon of the newly formed amide to create the heterocycle (10, 11), and then the former carbonyl oxygen is removed in a dehydration reaction, which introduces the carbon-nitrogen double bond of the thiazoline or (methyl)oxazoline ring. The nascent heterocyclic peptidyl-PCP can be used as the donor substrate by the next module’s C domain, or is first oxidized or reduced by discrete oxidase or reductase domains (10, 24).

Fig. 1.
Fig. 1.

Schematic representation of BmdB and the bacillamide E biosynthesis cycle. The A domains (orange) adenylate alanine and cysteine and transfer them onto the phosphopantetheine arm of the PCP domains (blue). The Cy2 domain (dark green) first catalyzes amide bond formation between the PCP-linked Ala and Cys residues, then catalyzes the intramolecular cyclodehydration reaction. The C3 domain (light green) catalyzes amide bond formation between the PCP-linked Ala-thiazoline moiety and free tryptamine, which releases the bacillamide E from the NRPS.

All C domain superfamily domains share the same protein fold, so the overall form of the Cy domain is not in doubt (25, 26). However, the features that allow the Cy domain to catalyze two separate and different reactions are not known. Cy domains contain a conserved Asp motif, DxxxxD, which directly replaces the catalytic His motif, HHxxxD, of C domains (6). Mutation of the aspartate residues of the Asp motif diminishes or abolishes catalytic activity of the protein (12, 23, 27, 28). Furthermore, other mutations have differential effects on the condensation reaction and the cyclodehydration reaction, suggesting that the reactions are not catalyzed by a completely overlapping set of residues (10, 22). Because Cy domains are not larger than C domains, the added function in Cy domains must occur in a confined sequence space.

Here, we present the crystal structure of Cy2 of bacillamide synthetase at 2.3 Å resolution. Bacillamide synthetase (Fig. 1) is a trimodular NRPS that produces bacillamide E, whose derivatives bacillamide A–D exhibit algicidal activity against dinoflagellates, raphidophytes, and cyanobacteria (SI Appendix, Fig. S1) (29, 30). Structural determination, along with mutagenic analysis of this Cy domain active site in the context of the full, intact, and active bacillamide synthetase, and bioinformatic investigation allowed us to identify D1226 and T1196 as important residues for cyclodehydration, providing a better understanding of the structure and mechanism of the Cy domain.

Results and Discussion

The Crystal Structure of an NRPS Heterocyclization Domain.

We have solved the crystal structure of an NRPS heterocyclization domain by X-ray crystallography to a resolution of 2.3 Å. Like all known C domain superfamily proteins, BmdB-Cy2 adopts two chloramphenicol acetyltransferase folds (2), with the N- and C-terminal lobes forming a pseudodimer (Fig. 2). In C domains, these two lobes assume a range of relative conformations (more “open” or “closed”; ref. 5). Superimposition of BmdB-Cy2 with each C domain shows that it fits within the observed range and is quite similar to AB3403 (31) (SI Appendix, Fig. S2). The latch formed by a crossover between N- and C-terminal subdomains in C domains is present, meaning that the Cy domain active site is also situated in the center of a tunnel, between donor and acceptor PCP binding sites (see SI Appendix, SI Results and Discussion for more detailed description of the overall Cy domain structure).

Fig. 2.
Fig. 2.

The crystal structure of BmdB-Cy2. (A) Ribbon representation of BmdB-Cy2. The cyclization domain adopts two chloramphenicol acetyltransferase folds, similar to C domains. (B) Close-up view of the BmdB-Cy2 DxxxxD motif (orange). D980 makes a bifurcated hydrogen bond with the hydroxyl of S873 and a hydrogen bond with the backbone nitrogen of L982. D985 makes hydrogen bonds with S988 and the amides of A987 and F1134, and water-mediated interaction with R1120.

BmdB-Cy2 Active Site.

The main established signature that differentiates Cy domains from C domains is an active site motif: C domains have a catalytic His motif of HHxxxD, and Cy domains have an Asp motif of DxxxxD (6). The first aspartate of the Asp motif is essential for catalytic activity in only some Cy domains. When the first aspartate is mutated to alanine in HMWP2-Cy1, AngN-Cy1, and AngN-Cy2, condensation and cyclodehydration are completely eliminated (23, 27), but in VibF-Cy1, VibF-Cy2, and EpoB-Cy1, both condensation and cyclodehydration occur, although they are significantly diminished (12, 23, 28) (SI Appendix, Table S1 lists current and previous Cy domain mutations). The second aspartate is critical for activity (12, 23, 27, 28), but whether its role is catalytic or structural (like the aspartate at that position in C domains) (32) was not definitively determined. Furthermore, a triple mutant of the Asp motif of HMWP2-Cy1 to introduce a C domain His motif results in a catalytically inactive protein (23), demonstrating that a straight swap of the motifs is not sufficient to interconvert catalytic activities.

In the structure of BmdB-Cy2, both aspartate residues in DxxxxD (D980 and D985) are oriented away from the active site (Fig. 2B). The side chain of D980 makes a bifurcated hydrogen bond with the hydroxyl of S873 and a hydrogen bond with the backbone amide of L982. D985 hydrogen bonds with the sidechain of S988 and the amides of A987 and F1134. The interaction this aspartate makes with a nearby arginine in C domains is present, but is water-mediated for D985-R1120 of BmdB-Cy2. Thus, these aspartate residues directly replace the first histidine and the aspartate in the C domain HHxxxD motif, occupying the same position in their respective domains and making the same or similar interactions with the rest of the protein (SI Appendix, Fig. S2B). Indeed, the whole DxxxxD motif is essentially in the same conformation as the HHxxxD motif. The position of the second histidine of the HHxxxD motif, which is the most important residue for condensation in C domains (2⇓⇓–5), is occupied by A981 in BmdB-Cy2 and, thus, unable to contribute to catalysis. The HHxxxD motif in C domains does not reorient upon substrate binding (5), and there is no indication that the Cy domain DxxxxD motif would do so. Overall, the structure strongly suggests that both of these aspartate residues play structural rather than catalytic roles, as has been established for the corresponding residues in C domains (2⇓⇓–5).

Bacillamide Synthesis Assay and Mutational Analysis.

To determine the importance of other active site residues in condensation and cyclization, we established a bacillamide synthesis assay. Bacillamide synthetase is a one-protein, six-domain, 265-kDa NRPS (Fig. 1) (33). It adenylates and thiolates its Ala and Cys substrates, after which the Cy2 domain performs condensation and cyclodehydration. The final domain in this NRPS is not a thioesterase domain, but a specialized C domain that condenses the thiazoline-containing intermediate with free tryptamine (Tpm) to release bacillamide E. Therefore, this terminal C domain plays all of the roles typically associated with a four-domain termination module—substrate selection, peptide bond formation, and release of the peptide product—which is rare, but not unprecedented in NRPSs (12, 27).

The Tpm in the product makes the reaction convenient to follow by HPLC at 280 nm (Fig. 3 and SI Appendix, Fig. S3). Wild-type bacillamide synthetase showed a large peak at 15 min that corresponded by high resolution mass spectrometry, fragmentation mass spectrometry, and NMR to bacillamide E (1) (SI Appendix, Fig. S3C). A minor peak at ∼15.5 min, with 10% of the intensity of the main peak, corresponded to the uncyclized tripeptide Ala-Cys-Tpm (2). This product shows that, at least in vitro, BmdB-Cy2 is not completely efficient at cyclizing all of the intermediates it condenses, and that the terminal C3 domain is not completely selective for the cyclized substrate. Furthermore, a small peak for the dipeptide Cys-Tpm is also evident (3; SI Appendix, Fig. S3), indicating that the terminal C3 domain also uses Cys-PCP2 as a donor substrate to some extent. Note that bacillamide E contains a thiazoline ring, and not a thiazole ring as in bacillamide A–D extracted from natural sources (29, 34, 35). The adjacent oxidase BmdC likely performs the oxidation postsynthetically.

Fig. 3.
Fig. 3.

Effects of structure-guided mutations in BmdB-Cy2 on bacillamide E synthesis. (A) Representative HPLC trace of a BmdB activity assay. Compound 1 is bacillamide E, compound 2 is linear Ala-Cys-Tpm, and compound 3 is Cys-Tpm. (B) Representative mass spectra of 1 (blue) and 2 (red). A representative mass spectrum of 3 is shown in SI Appendix, Fig. S3A. (C) Quantification of the relative production of 1 (blue) and 2 (red) in reactions with mutant BmdB. All reactions were measured in triplicate, except F1118G, which was measured in duplicate. Error bars represent SD.

We next produced a series of bacillamide synthetase constructs with mutations in the Cy domain. We targeted residues in spatial proximity to the DxxxxD motif in the active site tunnel of BmdB-Cy2, with side chains that could play an important role in condensation or cyclization, and which are largely conserved in alignments of characterized Cy domains (SI Appendix, Fig. S4). Therefore, bacillamide synthetase with mutations Y859F, T1116A, F1118G, T1135A, T1196A, D1226G, and D1226N, as well as N1114A and S1197A (two residues shown previously to selectively effect cyclodehydration; ref. 10) were purified and assayed.

The bacillamide synthetase mutants had radically different effects on condensation and cyclodehydration (Fig. 3 and SI Appendix, Fig. S3 and SI Results and Discussion). Mutant S1197A had almost no effect on production of linear or heterocyclized tripeptide. With Y859F, bacillamide E production decreased to 79% of wild type, but linear Ala-Cys-Tpm production doubled. Although the hydroxyl of Y859 points directly into the heart of the active site, the relatively minor effect of its removal suggests it is unlikely to act chemically. Mutants N1114A and T1116A each drastically reduced synthesis of both the linear and heterocyclic product. T1116A maintained approximately the same ratio of heterocyclic to linear product as wild type, whereas the little product made by N1114A was successfully heterocyclized, consistent with previously reported decoupling of condensation and heterocyclization by this residue (10). Both N1114 and T1116 thus appear important for condensation, and N1114 is also important for cyclodehydration.

For insight into the cyclization reaction, the most interesting mutations are those that selectively affect cyclodehydration. Mutation of residues F1118, T1135, T1196, and D1226 resulted in substantially more linear Ala-Cys-Tpm and less bacillamide E than the wild-type enzyme. F1118 is positioned along the tunnel that the phosphopantetheine arm of donor PCP1 occupies when presenting the Ala substrate to the active site (SI Appendix, Fig. S5). In the structure of BmdB-Cy2, the F1118 side chain completely blocks that tunnel. This residue must move to allow alanyl-PCP1 to bind and participate in the condensation reaction. However, for cyclization, PCP1 likely departs, leaving the cyclization substrate Ala-Cys-PCP2 bound at the acceptor site. It is possible that F1118 aids cyclization by blocking the donor side of the tunnel to create a single-entry active site more optimal for cyclodehydration. This feature is reminiscent of the (permanent) blocking of the acceptor site by a tryptophan side chain to form a single-entry active site in related epimerization domains, observed in a recently published crystal structure (25). Even more striking are the effects of the T1196A mutation, which allowed only minimal cyclization, and mutations of D1226 to glycine or asparagine, which obliterated cyclization while retaining robust condensation. These two residues are adjacent to one another, oriented toward the active site cavity, and hydrogen bond to one another. Their position and importance for cyclization activity make them the most likely to play a direct role in catalysis (see below).

Bioinformatic Analysis of Cy Domains.

We undertook a bioinformatic analysis of Cy domains. We retrieved 36,853 C domain superfamily sequences that had a maximum pairwise sequence identity of 90% and sorted them by using the Natural Products Domain Seeker web server (36). Multiple sequence alignment of 1,790 Cy domains showed two motifs in the C-terminal region that stand out as highly conserved in Cy domains but absent from the other C domain superfamily proteins (SI Appendix, Fig. S6). These motifs bear the consensus sequences PVVFTS and SQTPQVxLD (Fig. 4A) and are part of what had been recognized by Konz et al. (6) as conserved signature sequences 6 and 7. They exist in BmdB-Cy2 as PIVFTS (1192–1197) and ARTPQVYLD (1218–1226). These motifs are as close as 9 Å in space to the DxxxxD, but they are separated from it in sequence by more than 200 amino acids. The motifs form a surface on the acceptor substrate side of the active site (Fig. 4B and SI Appendix, Fig. S6A). Interestingly, this region is the putative location of the Ala-Cys-phosphopantetheinyl substrate in the cyclization reaction and includes the residues that displayed the most drastic effect on cyclization when mutated, namely T1196 of the PVVFTS motif and D1226 of the SQTPQVxLD motif. D1226 is one of the most conserved residues in Cy domains, present in 96% of the 1,790 sequences. T1196 is somewhat more variable: It is a threonine or serine in 88% of the sequences.

Fig. 4.
Fig. 4.

Model of the cyclodehydration intermediate and critical residues for cyclodehydration reaction. (A) WebLogo3 (57) of Cy domain motifs PVVFTS (core Cy6) and SQTPQVxLD (part of core Cy7) compared with the sequence of BmdB-Cy2. Putative catalytic residues T1196 and D1226 are labeled with yellow asterisks. (B) BmdB-Cy2 (mostly green, with DxxxxD in orange, PVVFTS in purple, and SQTPQVxLD in brown) with a model of the cyclodehydration intermediate (blue). Putative catalytic residues T1196 and D1226 are shown in sticks, as is V1228, a position occupied by glutamine in most Cy domains. (C) Possible mechanism of the dehydration step. We suggest that deprotonation of the amino hydrogen occurs after dehydration and double bond formation. The full putative reaction pathway is diagramed in SI Appendix, Fig. S7.

Putative functions had not been assigned to these motifs. S1197 was the only residue in either motif previously targeted by mutation, and it had differential effects on Cy domain reactivity in EpoB, BacA2, and BmdB, consistent with its moderate conservation (Figs. 3 and 4A and SI Appendix, Table S1) (10, 12). Notably, a portion of Cy domain sequence that largely overlaps with the two new motifs is assigned PFam08415. It is annotated only as being found in NRPSs with C and PCP domains, and ends partway through the SQTPQVxLD motif, before the critical D1226. We propose that this motif be extended to residue 1241 (to incorporate the whole SQTPQVxLD and a conserved WD; SI Appendix, Fig. S6A) and annotated as a signature for Cy domains.

A Trend for Tandem Cy Domains in (Methyl)oxazoline-Forming Modules.

Cy domains can be subdivided by their use of thiol (Cys) or hydroxyl (Ser or Thr) groups as the nucleophile in the cyclization reaction. Each subset contains the highly conserved PVVFTS and SQTPQVxLD motifs, and there are no discernable differences between Cy domain sequences in each subset. However, one clear trend did emerge: In general, modules that incorporate Cys substrates contain only one Cy domain per module, whereas those using Thr or Ser have tandem Cy-Cy (or Cy-C) domains. Of a set of 505 full-length Cy domain-containing proteins, 454 are predicted to use Cys (440 with single Cy domains, 14 with Cy-Cy domains), whereas 51 are predicted to use Ser or Thr (2 with a single Cy domain, 49 with Cy-Cy domains). This trend is maintained after the proteins have been filtered for unique architecture and/or unique predicted product: Of 200 remaining proteins, 187 use Cys (177 with a single Cy domain, 10 with Cy-Cy domains) and 13 use Ser or Thr (1 with single a Cy domain, 12 with Cy-Cy domains).

The trend of tandem Cy domains for hydroxyl-containing substrates and single Cy domains for thiol-containing substrates holds for most characterized NRPS systems that feature Cy domains (6⇓–8, 10, 12, 22, 23, 37⇓⇓⇓–41). The tandem Cy domain arrangement is exemplified by the well-characterized VibF protein. In VibF, Cy2 is primarily responsible for condensation and Cy1 for cyclodehydration of the Thr substrate (22). Furthermore, the tandem Cy arrangement can be maintained even when the module is split between two proteins, as in serratiochelin synthetase SchF1 and SchF2, which also cyclodehydrates Thr (40, 41). Exceptions to the trend are mycobactin synthetase (42) (Thr-specific module with a single Cy domain in MbtB) and anguibactin synthetase (Cys-specific module with Cy-Cy domains in AngN) (27). The tandem Cy domains in anguibactin synthetase could be a remnant of evolution from a Thr-using ancestor, because, other than containing a Cys in place of a Thr, actinomycin is identical to acinetobactin, and the two synthetases share identical domain configurations (43). Furthermore, the Cy domains in AngN can each perform both condensation and cyclodehydration reactions and are largely redundant with one another (27).

The Cys thiol is a better nucleophile than the Ser or Thr hydroxyl (44). This difference in reactivity may be why thiazoline-forming NRPS modules are ∼10× more prevalent than (methyl)oxazoline-forming modules and may explain the trend to dedicate two domains in tandem in (methyl)oxazoline-forming modules. Tandem domains may increase the probability of cyclodehydration before the peptide is donated in the downstream module’s condensation reaction. The downstream C domains are likely not completely selective for the heterocycle-containing peptide, as shown with BmdB (Fig. 3) and VibF (28), so efficiency of cyclodehydration would be important for cognate heterocycle production.

Insight into Catalysis and Model of the Cyclodehydration Intermediate.

The first reaction performed by Cy domains is condensation. The N1114A and T1116A mutants of BmdB-Cy2 exhibited the greatest effects on condensation without fully abolishing it (Fig. 3), and neither residue is highly conserved (SI Appendix, Fig. S6A). Likewise, all mutations that were reported to decrease condensation in Cy domains are shown by the BmdB-Cy2 structure to be too distal to (residues 900, 988, 1089, 1114, 1120) or turned away from (residues 980, 985 of the DxxxxD motif) the atoms participating in condensation (SI Appendix, Table S1) (10, 12, 23, 27, 28). Unless the DxxxxD motif radically reorients (which is possible, although there is no supporting data), there do not appear to be any critical, conserved residues that could act to abstract or donate a proton in the condensation reaction in Cy domains. We have argued that substrate orientation is the principal source of catalytic power for condensation in C domains (21), and this hypothesis appears to be even more likely in Cy domains. Comparison of the kinetic parameters of uncatalyzed peptide bond formation and ribosome-catalyzed peptide bond formation shows that the rate enhancement provided by the ribosome comes predominantly from the lowering of activation entropy via positioning of ester substrates (and perhaps water molecules) for reaction (45). Thus, Cy domains should be able to rely principally on substrate positioning to catalyze condensation of their more reactive thioester substrates.

The second reaction performed by Cy domains, cyclodehydration, is a nonsimplistic, two-step reaction. Cy domains and at least two other types of enzymes perform this reaction to form five-membered rings in independent natural product biosynthesis systems. YcaO- and TruD-type enzymes processively modify ribosomally synthesized and posttranslationally modified peptides by using cyclodehydration reactions to introduce thiazoline and (methyl)oxazoline rings in the synthesis of peptides such as microcin and the trunkamides (46). YcaO and TruD, which are structurally unrelated to Cy domains, covalently attach a phosphate or adenylate from ATP to the carbonyl oxygen to promote both cyclization and dehydration (46). In Cy domains, which do not use ATP, the energetics of cyclodehydration is presumably linked to the high-energy thioester bond that is broken in the subsequent step of bacillamide synthesis. Cy domains share this characteristic with the newly described specialty C domain, NocB-C5, that catalyzes β-lactam ring formation (47). NocB-C5 is proposed to dehydrate a serine to β-alanine via base catalysis by a histidine directly upstream of its HHxxxD motif, prior to cyclization. Dehydration before cyclization is the reaction order shared with lantibiotic cyclodehydratases, but not with Cy domains (48). The position of the upstream histidine is occupied by a hydrophobic residue in Cy domains, and this V979 in BmdB-Cy2 is somewhat recessed from the active site. Thus, Cy domains rely on a completely different mechanism for catalysis than YcaO, TruD, NcoB-C5, or the lantibiotic cyclodehydratases.

To help integrate our structural, mutagenic, and conservation data, we created a model of the intermediate of the BmdB-catalyzed cyclodehydration reaction (Fig. 4B). The configuration of the model and the present data are consistent with the absolutely critical catalytic action of D1226 in the first step of the cyclodehydration reaction, orienting and abstracting a proton from the thiol (or hydroxyl in Ser- or Thr-specific Cy domains) (SI Appendix, Fig. S7). T1196, 3.0 Å from D1226, appears well positioned to aid in catalysis by interacting with D1226, and may donate a proton to the former carbonyl oxygen to form the hydroxyl-thiazolidine. D1226 may protonate the same oxygen again to allow it to leave as water in the dehydration reaction (Fig. 4 B and C and SI Appendix, Fig. S7). Upon dehydration, the pKa of the amino proton is lowered such that it can be facilely lost to solvent to produce the final thiazoline moiety.

Our mutational results indicate that D1226 is critical for cyclodehydration, and T1196 is nearly as important (Fig. 3A). However, T1196 is only a threonine or serine in 88% of Cy sequences, and appears as an alanine in 7% of cases, including in the functionally characterized PchE-Cy1 (39). How can its apparent importance in the reaction and its relative lack of conservation be reconciled? Homology modeling of PchE-Cy1 onto BmdB-Cy2 provides a potential answer: The residue at the nearby position equivalent to 1228 of PchE-Cy1 is a glutamine, which could compensate for the missing T1196. A double mutant of T1196A/V1228Q did not restore cyclization activity to BmdB (SI Appendix, Fig. S8), but is not surprising given the BmdB-Cy has evolved to perform cyclodehydration with T1196 and not Q1228. In the analyzed Cy domain sequences that have an alanine at position 1196, all but two have a glutamine at 1228 (with the remaining two having glutamate or asparagine). In the majority of Cy domains, all three residues are present to form a T-D-Q triad. A similar S/T-D-Q catalytic triad in Escherichia coli thioesterase II activates water for nucleophilic attack (49). We suggest that mutation of the threonine or glutamine in the Cy domain triad may be tolerated in some instances, because chemistry is not rate-limiting in the overall slow synthetic cycle of NRPSs, and at least in thiazoline-forming Cy domains, it would not be difficult to deprotonate the Cys side chain (pKa of ∼8.2 in aqueous solution). Deprotonation of Ser and Thr side chains is more difficult and could require the full catalytic triad. Variation in catalytic residues is not abnormal in nature, as exemplified by serine protease active sites (50). Analogous to the observed variation of T1196, thioesterase II enzymes feature hydrophobic residues in place of the catalytic serine/threonine in more than 5% of proteins.

Finally, parallels can be drawn between the role of D1226 in deprotonating the substrate side chain for nucleophilic attack with the mechanisms of microbial transglutaminase, cytosolic phospholipase A2, patatin, and TEM-1 β-lactamase, in which the aspartate of a C-D or S-D dyad deprotonates the cysteine or serine for nucleophilic attack on a carbonyl carbon (51⇓–53).

During the review process for this study, a structure and mutational analysis of the heterocyclization domain of epothilone synthase (EpoB-Cy) was reported by Dowling et al. (54). The structures of EpoB-Cy and BmdB-Cy2 are similar. Their mutational data are consistent with those presented here, in particular the identification of D1226 (D449 in EpoB-Cy) as important for catalysis, although without differentiation between defects in condensation and heterocyclization (SI Appendix, SI Results and Discussion and Table S1). Their study and ours complement one another and provide a greater understanding of NRPS heterocyclization domains.

Conclusion

In summary, we have presented the crystal structure of an NRPS heterocyclization domain, solved to a resolution of 2.3 Å. We cloned, expressed, and purified the entire bacillamide synthetase containing this Cy domain, and used it to assay effects of Cy domain mutations on peptide production. We were able to identify two residues, T1196 and D1226, which are important for the cyclodehydration reaction. Finally, we presented a putative mechanism for cyclodehydration in which D1226 acts as general acid/base catalyst.

Methods

Bacillamide Synthetase Cy2 Crystallography.

The Thermoactinomyces vulgaris BmdB-Cy2 construct, with an N-terminal octahistidine tag, tobacco etch virus (TEV) protease cleavage site, and BmdB residues 844–1287, was synthesized by DNA 2.0. BmdB-Cy2 was heterologously expressed in E. coli and purified to homogeneity (SI Appendix, SI Methods). BmdB-Cy2 was crystallized by using 0.88% Tween 20, 1.62 M ammonium sulfate, 0.1 M Hepes pH 7.5, 2.67% (wt/vol) PEG 400, and 3% (wt/vol) 6-aminohexanoic acid, with 20% (vol/vol) ethylene glycol, added before flash-cooling. Diffraction data were collected at CLS beamline 08ID-1 (SI Appendix, Table S2) and the structure of BmdB-Cy2 was determined by molecular replacement.

Peptide Synthesis Assay for Full-Length BmdB.

bmdB was cloned from T. vulgaris F-5595 genomic DNA into a pET21-derived vector containing an N-terminal TEV-cleavable calmodulin binding peptide tag and a C-terminal TEV-cleavable octa-histidine tag. Mutations were introduced via site-directed mutagenesis (SI Appendix, Table S3). Wild-type and mutant protein was expressed in E. coli cells and purified to homogeneity. In a 1-mL reaction, 50 mM Tris⋅HCl pH 7.5, 100 mM NaCl, and 10 mM MgCl2, 1 mM Ala, 1 mM Cys, 10 mM tryptamine, 2 mM ATP, and 100 nM BmdB were incubated for 2 h at 37 °C. Reactions were analyzed by HPLC using C18 media and a gradient of water/0.1% TFA to acetonitrile/0.1% TFA. High-resolution MS and fragmentation MS analysis was performed at the Proteomics Platform at the Research Institute of the McGill University Health Centre, and NMR analysis at Quebec/Eastern Canada High Field NMR Facility (QANUC), McGill University.

Bioinformatic Analysis.

Sequences of known C domains were queried against the nr90 database (55), using low threshold search. Matched sequences were themselves filtered for maximum pairwise sequence identity of 90%, which gave 36,853 C domain superfamily sequences. These sequences were classified into subtypes by using NaPDoS (36), and included 1,790 Cy domains. To look for trends in Cy domains by substrates, all full-length protein sequences from which the 1,790 Cy domains originate were retrieved. Sequences were discarded if there was not an A domain adjacent to the Cy domain, or if they started as a single Cy domain leaving a set of 505 proteins. The Cy domain acceptor/cyclodehydration substrates were inferred from substrate of the A domains in the same module, as predicted by the program ANTISMASH (56).

Acknowledgments

We thank Chaitan Khosla and James Kuo for the kind gift of BAP1 cells; Vikram Alva and Johannes Soeding for providing the current nr90 database; Christopher Boddy and Adrian Keatinge-Clay for conversations; Asfarul Haque for performing the ITC experiments; Janice Reimer for schematic figure design; the Proteomics Platform at the RI-MUHC (Montreal, Canada) for LC-MS analysis; members of the T.M.S. laboratory for helpful discussions; Alexander Wahba for MS analysis; Varoujan Yaylayan for fragmentation MS assignment; and Robin Stein for NMR analysis.

Footnotes

  • ↵1To whom correspondence should be addressed. Email: martin.schmeing{at}mcgill.ca.

  • Author contributions: K.B. and T.M.S. designed research; K.B. and T.M.S. performed research; K.B., C.D.F., M.A.M., and T.M.S. analyzed data; and K.B. and T.M.S. wrote the paper with comments from C.D.F. and M.A.M.

  • The authors declare no conflict of interest.

  • This article is a PNAS Direct Submission.

  • Data deposition: The atomic coordinates and structure factors have been deposited in the Protein Data Bank, www.pdb.org (PDB ID code 5T3E).

  • This article contains supporting information online at www.pnas.org/lookup/suppl/doi:10.1073/pnas.1614191114/-/DCSupplemental.

References

  1. ↵
    1. Weissman KJ
    (2015) The structural biology of biosynthetic megaenzymes. Nat Chem Biol 11(9):660–670.
    .
    OpenUrlCrossRefPubMed
  2. ↵
    1. Keating TA,
    2. Marshall CG,
    3. Walsh CT,
    4. Keating AE
    (2002) The structure of VibH represents nonribosomal peptide synthetase condensation, cyclization and epimerization domains. Nat Struct Biol 9(7):522–526.
    .
    OpenUrlCrossRefPubMed
  3. ↵
    1. Samel SA,
    2. Schoenafinger G,
    3. Knappe TA,
    4. Marahiel MA,
    5. Essen LO
    (2007) Structural and functional insights into a peptide bond-forming bidomain from a nonribosomal peptide synthetase. Structure 15(7):781–792.
    .
    OpenUrlCrossRefPubMed
  4. ↵
    1. Bergendahl V,
    2. Linne U,
    3. Marahiel MA
    (2002) Mutational analysis of the C-domain in nonribosomal peptide synthesis. Eur J Biochem 269(2):620–629.
    .
    OpenUrlPubMed
  5. ↵
    1. Bloudoff K,
    2. Rodionov D,
    3. Schmeing TM
    (2013) Crystal structures of the first condensation domain of CDA synthetase suggest conformational changes during the synthetic cycle of nonribosomal peptide synthetases. J Mol Biol 425(17):3137–3150.
    .
    OpenUrlCrossRefPubMed
  6. ↵
    1. Konz D,
    2. Klens A,
    3. Schörgendorfer K,
    4. Marahiel MA
    (1997) The bacitracin biosynthesis operon of Bacillus licheniformis ATCC 10716: Molecular characterization of three multi-modular peptide synthetases. Chem Biol 4(12):927–937.
    .
    OpenUrlCrossRefPubMed
  7. ↵
    1. Marshall CG,
    2. Burkart MD,
    3. Keating TA,
    4. Walsh CT
    (2001) Heterocycle formation in vibriobactin biosynthesis: Alternative substrate utilization and identification of a condensed intermediate. Biochemistry 40(35):10655–10663.
    .
    OpenUrl
  8. ↵
    1. Chen H,
    2. O’Connor S,
    3. Cane DE,
    4. Walsh CT
    (2001) Epothilone biosynthesis: Assembly of the methylthiazolylcarboxy starter unit on the EpoB subunit. Chem Biol 8(9):899–912.
    .
    OpenUrlCrossRefPubMed
  9. ↵
    1. Keating TA,
    2. Marshall CG,
    3. Walsh CT
    (2000) Reconstitution and characterization of the Vibrio cholerae vibriobactin synthetase from VibB, VibE, VibF, and VibH. Biochemistry 39(50):15522–15530.
    .
    OpenUrlCrossRefPubMed
  10. ↵
    1. Duerfahrt T,
    2. Eppelmann K,
    3. Müller R,
    4. Marahiel MA
    (2004) Rational design of a bimodular model system for the investigation of heterocyclization in nonribosomal peptide biosynthesis. Chem Biol 11(2):261–271.
    .
    OpenUrlCrossRefPubMed
  11. ↵
    1. Gehring AM,
    2. Mori I,
    3. Perry RD,
    4. Walsh CT
    (1998) The nonribosomal peptide synthetase HMWP2 forms a thiazoline ring during biogenesis of yersiniabactin, an iron-chelating virulence factor of yersinia pestis. Biochemistry 37(48):17104.
    .
    OpenUrlCrossRefPubMed
  12. ↵
    1. Kelly WL,
    2. Hillson NJ,
    3. Walsh CT
    (2005) Excision of the epothilone synthetase B cyclization domain and demonstration of in trans condensation/cyclodehydration activity. Biochemistry 44(40):13385–13393.
    .
    OpenUrlCrossRefPubMed
  13. ↵
    1. Tabata N,
    2. Tomoda H,
    3. Zhang H,
    4. Uchida R,
    5. Omura S
    (1999) Zelkovamycin, a new cyclic peptide antibiotic from Streptomyces sp. K96-0670. II. Structure elucidation. J Antibiot (Tokyo) 52(1):34–39.
    .
    OpenUrlPubMed
  14. ↵
    1. Shen B, et al.
    (2002) Cloning and characterization of the bleomycin biosynthetic gene cluster from Streptomyces verticillus ATCC15003. J Nat Prod 65(3):422–431.
    .
    OpenUrlCrossRefPubMed
  15. ↵
    1. Vollbrecht L, et al.
    (2002) Argyrins, immunosuppressive cyclic peptides from myxobacteria. II. Structure elucidation and stereochemistry. J Antibiot (Tokyo) 55(8):715–721.
    .
    OpenUrlCrossRefPubMed
  16. ↵
    1. Snow GA
    (1970) Mycobactins: Iron-chelating growth factors from mycobacteria. Bacteriol Rev 34(2):99–125.
    .
    OpenUrlFREE Full Text
  17. ↵
    1. Drechsel H, et al.
    (1995) Structure elucidation of Yersiniabactin, a siderophore from highly virulent yersinia strains. Liebigs Ann 1995(10):1727–1733.
    .
    OpenUrl
  18. ↵
    1. Vizcaino MI,
    2. Crawford JM
    (2015) The colibactin warhead crosslinks DNA. Nat Chem 7(5):411–417.
    .
    OpenUrlCrossRefPubMed
  19. ↵
    1. Roy RS,
    2. Gehring AM,
    3. Milne JC,
    4. Belshaw PJ,
    5. Walsh CT
    (1999) Thiazole and oxazole peptides: Biosynthesis and molecular machinery. Nat Prod Rep 16(2):249–263.
    .
    OpenUrlCrossRefPubMed
  20. ↵
    1. Rausch C,
    2. Hoof I,
    3. Weber T,
    4. Wohlleben W,
    5. Huson DH
    (2007) Phylogenetic analysis of condensation domains in NRPS sheds light on their functional evolution. BMC Evol Biol 7:78.
    .
    OpenUrlCrossRefPubMed
  21. ↵
    1. Bloudoff K,
    2. Alonzo DA,
    3. Schmeing TM
    (2016) Chemical probes allow structural insight into the condensation reaction of nonribosomal peptide synthetases. Cell Chem Biol 23(3):331–339.
    .
    OpenUrl
  22. ↵
    1. Marshall CG,
    2. Hillson NJ,
    3. Walsh CT
    (2002) Catalytic mapping of the vibriobactin biosynthetic enzyme VibF. Biochemistry 41(1):244–250.
    .
    OpenUrlCrossRefPubMed
  23. ↵
    1. Keating TA,
    2. Miller DA,
    3. Walsh CT
    (2000) Expression, purification, and characterization of HMWP2, a 229 kDa, six domain protein subunit of Yersiniabactin synthetase. Biochemistry 39(16):4729–4739.
    .
    OpenUrlCrossRefPubMed
  24. ↵
    1. Patel HM,
    2. Walsh CT
    (2001) In vitro reconstitution of the Pseudomonas aeruginosa nonribosomal peptide synthesis of pyochelin: Characterization of backbone tailoring thiazoline reductase and N-methyltransferase activities. Biochemistry 40(30):9023–9031.
    .
    OpenUrlCrossRefPubMed
  25. ↵
    1. Chen WH,
    2. Li K,
    3. Guntaka NS,
    4. Bruner SD
    (2016) Interdomain and intermodule organization in epimerization domain containing nonribosomal peptide synthetases. ACS Chem Biol 11(8):2293–2303.
    .
    OpenUrlCrossRefPubMed
  26. ↵
    1. Samel SA,
    2. Czodrowski P,
    3. Essen LO
    (2014) Structure of the epimerization domain of tyrocidine synthetase A. Acta Crystallogr D Biol Crystallogr 70(Pt 5):1442–1452.
    .
    OpenUrlCrossRef
  27. ↵
    1. Di Lorenzo M,
    2. Stork M,
    3. Naka H,
    4. Tolmasky ME,
    5. Crosa JH
    (2008) Tandem heterocyclization domains in a nonribosomal peptide synthetase essential for siderophore biosynthesis in Vibrio anguillarum. Biometals 21(6):635–648.
    .
    OpenUrlCrossRefPubMed
  28. ↵
    1. Marshall CG,
    2. Burkart MD,
    3. Meray RK,
    4. Walsh CT
    (2002) Carrier protein recognition in siderophore-producing nonribosomal peptide synthetases. Biochemistry 41(26):8429–8437.
    .
    OpenUrl
  29. ↵
    1. Jeong SY,
    2. Ishida K,
    3. Ito Y,
    4. Okada S,
    5. Murakami M
    (2003) Bacillamide, a novel algicide from the marine bacterium, Bacillus sp SY-1, against the harmful dinoflagellate, Cochlodinium polykrikoides. Tetrahedron Lett 44(43):8005–8007.
    .
    OpenUrlCrossRef
  30. ↵
    1. Churro C, et al.
    (2009) Effects of bacillamide and newly synthesized derivatives on the growth of cyanobacteria and microalgae cultures. J Appl Phycol 21(4):429–442.
    .
    OpenUrlCrossRef
  31. ↵
    1. Drake EJ, et al.
    (2016) Structures of two distinct conformations of holo-non-ribosomal peptide synthetases. Nature 529(7585):235–238.
    .
    OpenUrlCrossRefPubMed
  32. ↵
    1. Roche ED,
    2. Walsh CT
    (2003) Dissection of the EntF condensation domain boundary and active site residues in nonribosomal peptide synthesis. Biochemistry 42(5):1334–1344.
    .
    OpenUrlCrossRefPubMed
  33. ↵
    1. Yuwen L, et al.
    (2013) The role of aromatic L-amino acid decarboxylase in bacillamide C biosynthesis by Bacillus atrophaeus C89. Sci Rep 3:1753.
    .
    OpenUrl
  34. ↵
    1. Ivanova V, et al.
    (2007) Microbiaeratin, a new natural indole alkaloid from a Microbispora aerata strain, isolated from Livingston Island, Antarctica. Prep Biochem Biotechnol 37(2):161–168.
    .
    OpenUrlCrossRefPubMed
  35. ↵
    1. Socha AM,
    2. Long RA,
    3. Rowley DC
    (2007) Bacillamides from a hypersaline microbial mat bacterium. J Nat Prod 70(11):1793–1795.
    .
    OpenUrlCrossRefPubMed
  36. ↵
    1. Ziemert N, et al.
    (2012) The natural product domain seeker NaPDoS: A phylogeny based bioinformatic tool to classify secondary metabolite gene diversity. PLoS One 7(3):e34064.
    .
    OpenUrlCrossRefPubMed
  37. ↵
    1. Du L,
    2. Chen M,
    3. Zhang Y,
    4. Shen B
    (2003) BlmIII and BlmIV nonribosomal peptide synthetase-catalyzed biosynthesis of the bleomycin bithiazole moiety involving both in cis and in trans aminoacylation. Biochemistry 42(32):9731–9740.
    .
    OpenUrl
  38. ↵
    1. Silakowski B, et al.
    (1999) New lessons for combinatorial biosynthesis from myxobacteria. The myxothiazol biosynthetic gene cluster of Stigmatella aurantiaca DW4/3-1. J Biol Chem 274(52):37391–37399.
    .
    OpenUrlAbstract/FREE Full Text
  39. ↵
    1. Quadri LE,
    2. Keating TA,
    3. Patel HM,
    4. Walsh CT
    (1999) Assembly of the Pseudomonas aeruginosa nonribosomal peptide siderophore pyochelin: In vitro reconstitution of aryl-4, 2-bisthiazoline synthetase activity from PchD, PchE, and PchF. Biochemistry 38(45):14941–14954.
    .
    OpenUrlCrossRefPubMed
  40. ↵
    1. Seyedsayamdost MR, et al.
    (2012) Mixing and matching siderophore clusters: Structure and biosynthesis of serratiochelins from Serratia sp. V4. J Am Chem Soc 134(33):13550–13553.
    .
    OpenUrlCrossRefPubMed
  41. ↵
    1. Wang H,
    2. Fewer DP,
    3. Holm L,
    4. Rouhiainen L,
    5. Sivonen K
    (2014) Atlas of nonribosomal peptide and polyketide biosynthetic pathways reveals common occurrence of nonmodular enzymes. Proc Natl Acad Sci USA 111(25):9259–9264.
    .
    OpenUrlAbstract/FREE Full Text
  42. ↵
    1. McMahon MD,
    2. Rush JS,
    3. Thomas MG
    (2012) Analyses of MbtB, MbtE, and MbtF suggest revisions to the mycobactin biosynthesis pathway in Mycobacterium tuberculosis. J Bacteriol 194(11):2809–2818.
    .
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Dorsey CW, et al.
    (2004) The siderophore-mediated iron acquisition systems of Acinetobacter baumannii ATCC 19606 and Vibrio anguillarum 775 are structurally and functionally related. Microbiology 150(Pt 11):3657–3667.
    .
    OpenUrlCrossRefPubMed
  44. ↵
    1. Belshaw PJ,
    2. Roy RS,
    3. Kelleher NL,
    4. Walsh CT
    (1998) Kinetics and regioselectivity of peptide-to-heterocycle conversions by microcin B17 synthetase. Chem Biol 5(7):373–384.
    .
    OpenUrlCrossRefPubMed
  45. ↵
    1. Beringer M, et al.
    (2005) Essential mechanisms in the catalysis of peptide bond formation on the ribosome. J Biol Chem 280(43):36065–36072.
    .
    OpenUrlAbstract/FREE Full Text
  46. ↵
    1. Dunbar KL,
    2. Melby JO,
    3. Mitchell DA
    (2012) YcaO domains use ATP to activate amide backbones during peptide cyclodehydrations. Nat Chem Biol 8(6):569–575.
    .
    OpenUrlCrossRefPubMed
  47. ↵
    1. Gaudelli NM,
    2. Long DH,
    3. Townsend CA
    (2015) β-Lactam formation by a non-ribosomal peptide synthetase during antibiotic biosynthesis. Nature 520(7547):383–387.
    .
    OpenUrlCrossRefPubMed
  48. ↵
    1. Zhang Q,
    2. Yu Y,
    3. Vélasquez JE,
    4. van der Donk WA
    (2012) Evolution of lanthipeptide synthetases. Proc Natl Acad Sci USA 109(45):18361–18366.
    .
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Li J,
    2. Derewenda U,
    3. Dauter Z,
    4. Smith S,
    5. Derewenda ZS
    (2000) Crystal structure of the Escherichia coli thioesterase II, a homolog of the human Nef binding enzyme. Nat Struct Biol 7(7):555–559.
    .
    OpenUrlCrossRefPubMed
  50. ↵
    1. Ekici OD,
    2. Paetzel M,
    3. Dalbey RE
    (2008) Unconventional serine proteases: Variations on the catalytic Ser/His/Asp triad configuration. Protein Sci 17(12):2023–2037.
    .
    OpenUrlCrossRefPubMed
  51. ↵
    1. Dessen A, et al.
    (1999) Crystal structure of human cytosolic phospholipase A2 reveals a novel topology and catalytic mechanism. Cell 97(3):349–360.
    .
    OpenUrlCrossRefPubMed
  52. ↵
    1. Rydel TJ, et al.
    (2003) The crystal structure, mutagenesis, and activity studies reveal that patatin is a lipid acyl hydrolase with a Ser-Asp catalytic dyad. Biochemistry 42(22):6696–6708.
    .
    OpenUrlCrossRefPubMed
  53. ↵
    1. Minasov G,
    2. Wang X,
    3. Shoichet BK
    (2002) An ultrahigh resolution structure of TEM-1 beta-lactamase suggests a role for Glu166 as the general base in acylation. J Am Chem Soc 124(19):5333–5340.
    .
    OpenUrlCrossRefPubMed
  54. ↵
    1. Dowling DP, et al.
    (2016) Structural elements of an NRPS cyclization domain and its intermodule docking domain. Proc Natl Acad Sci USA 113(44):12432–12437.
    .
    OpenUrlAbstract/FREE Full Text
  55. ↵
    1. Biegert A,
    2. Söding J
    (2009) Sequence context-specific profiles for homology searching. Proc Natl Acad Sci USA 106(10):3770–3775.
    .
    OpenUrlAbstract/FREE Full Text
  56. ↵
    1. Weber T, et al.
    (2015) antiSMASH 3.0-a comprehensive resource for the genome mining of biosynthetic gene clusters. Nucleic Acids Res 43(W1):W237-43.
    .
    OpenUrl
  57. ↵
    1. Crooks GE,
    2. Hon G,
    3. Chandonia JM,
    4. Brenner SE
    (2004) WebLogo: A sequence logo generator. Genome Res 14(6):1188–1190.
    .
    OpenUrlAbstract/FREE Full Text
View Abstract
Back to top
Article Alerts
Email Article
Citation Tools
Request Permissions
Structural insight into heterocyclization domains
Kristjan Bloudoff, Christopher D. Fage, Mohamed A. Marahiel, T. Martin Schmeing
Proceedings of the National Academy of Sciences Jan 2017, 114 (1) 95-100; DOI: 10.1073/pnas.1614191114
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Facebook logo Google logo Mendeley logo
Proceedings of the National Academy of Sciences: 116 (7)
Current Issue

Submit

Jump to section

You May Also be Interested in

Nonmonogamous strawberry poison frog (Oophaga pumilio). Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Putative signature of monogamy
A study suggests a putative gene-expression hallmark common to monogamous male vertebrates of some species, namely cichlid fishes, dendrobatid frogs, passeroid songbirds, common voles, and deer mice, and identifies 24 candidate genes potentially associated with monogamy.
Image courtesy of Yusan Yang (University of Pittsburgh, Pittsburgh).
Active lifestyles. Image courtesy of Pixabay/MabelAmber.
Meaningful life tied to healthy aging
Physical and social well-being in old age are linked to self-assessments of life worth, and a spectrum of behavioral, economic, health, and social variables may influence whether aging individuals believe they are leading meaningful lives.
Image courtesy of Pixabay/MabelAmber.