Self-assembly of amphiphilic Janus dendrimers into uniform onion-like dendrimersomes with predictable size and number of bilayers

A constitutional isomeric library synthesized by a modular approach has been used to discover six amphiphilic Janus dendrimer primary structures, which self-assemble into uniform onion-like vesicles with predictable dimensions and number of internal bilayers. These vesicles, denoted onion-like dendrimersomes, are assembled by simple injection of a solution of Janus dendrimer in a water-miscible solvent into water or buffer. These dendrimersomes provide mimics of double-bilayer and multibilayer biological membranes with dimensions and number of bilayers predicted by the Janus compound concentration in water. The simple injection method of preparation is accessible without any special equipment, generating uniform vesicles, and thus provides a promising tool for fundamental studies as well as technological applications in nanomedicine and other fields.Most living organisms contain single-bilayer membranes composed of lipids, glycolipids, cholesterol, transmembrane proteins, and glycoproteins (1). Gram-negative bacteria (2, 3) and the cell nucleus (4), however, exhibit a strikingly special envelope that consists of a concentric double-bilayer membrane. More complex membranes are also encountered in cells and their various organelles, such as multivesicular structures of eukaryotic cells (5) and endosomes (6), and multibilayer structures of endoplasmic reticulum (7, 8), myelin (9, 10), and multilamellar bodies (11, 12). This diversity of biological membranes inspired corresponding biological mimics. Liposomes (Fig. 1) self-assembled from phospholipids are the first mimics of single-bilayer biological membranes (13–16), but they are polydisperse, unstable, and permeable (14). Stealth liposomes coassembled from phospholipids, cholesterol, and phospholipids conjugated with poly(ethylene glycol) exhibit improved stability, permeability, and mechanical properties (17–20). Polymersomes (21–24) assembled from amphiphilic block copolymers exhibit better mechanical properties and permeability, but are not always biocompatible and are polydisperse. Dendrimersomes (25–28) self-assembled from amphiphilic Janus dendrimers and minidendrimers (26–28) have also been elaborated to mimic single-bilayer biological membranes. Amphiphilic Janus dendrimers take advantage of multivalency both in their hydrophobic and hydrophilic parts (23, 29–32). Dendrimersomes are assembled by simple injection (33) of a solution of an amphiphilic Janus dendrimer (26) in a water-soluble solvent into water or buffer and produce uniform (34), impermeable, and stable vesicles with excellent mechanical properties. In addition, their size and properties can be predicted by their primary structure (27). Amphiphilic Janus glycodendrimers self-assemble into glycodendrimersomes that mimic the glycan ligands of biological membranes (35). They have been demonstrated to be bioactive toward biomedically relevant bacterial, plant, and human lectins, and could have numerous applications in nanomedicine (20).Open in a separate windowFig. 1.Strategies for the preparation of single-bilayer vesicles and multibilayer onion-like vesicles.More complex and functional cell mimics such as multivesicular vesicles (36, 37) and multibilayer onion-like vesicles (38–40) have also been discovered. Multivesicular vesicles compartmentalize a larger vesicle (37) whereas multibilayer onion-like vesicles consist of concentric alternating bilayers (40). Currently multibilayer vesicles are obtained by very complex and time-consuming methods that do not control their size (39) and size distribution (40) in a precise way. Here we report the discovery of “single–single” (28) amphiphilic Janus dendrimer primary structures that self-assemble into uniform multibilayer onion-like dendrimersomes (Fig. 1) with predictable size and number of bilayers by simple injection of their solution into water or buffer.     

Solid-phase peptide synthesis and solid-phase fragment coupling mediated by isonitriles

The synthesis of polypeptides on solid phase via mediation by isonitriles is described. The acyl donor is a thioacid, which presumably reacts with the isonitrile to generate a thio-formimidate carboxylate mixed anhydride intermediate. Applications of this chemistry to reiterative solid-phase peptide synthesis as well as solid-phase fragment coupling are described.Amide bond formations are arguably among the most important constructions in organic chemistry (1, 2). The centrality of the amide linkage, as found in polypeptides and proteins, in the maintenance of life hardly needs restatement. Numerous strategies, resulting in a vast array of protocols to synthesize biologically active polypeptides and proteins, have been demonstrated (3, 4). Central to reiterative polypeptide bond formations was the discovery and remarkable development of solid-phase peptide synthesis (SPPS) (5, 6). The extraordinary impact of SPPS in fostering enhanced access to homogeneous polypeptides is clear to everyone in the field.As we have described elsewhere, by classical, mechanistic reasoning, we were led to conjecture about some hitherto-unexplored possibilities relevant to the chemistry of isonitriles (7–14). It was anticipated that isonitriles might be able to mediate the acylation of amines, thus giving rise to amides (15). Early experiments focused on free carboxylic acids as the acylating agents. As our studies progressed, it was found that the combination of thioacids, amines, and isonitriles leads to the efficient formation of amide bonds under stoichiometric or near-stoichiometric conditions (7–13, 16, 17). Although there remain unresolved issues of detail and nuance, the governing mechanism for amide formation under these conditions involves reaction of the thioacid, 1, with an isonitrile, 2, to generate a thio-formimidate carboxylate mixed anhydride (thio-FCMA), 3, which is intercepted by the “acyl-accepting” amine to generate amide, 5, and thioformamide, 6 (Fig. 1). The efficiency of the amidation was further improved through the use of hydroxybenzotriazole (HOBt) (18), which could well give rise to HOBt ester 7, although this pathway has not been mechanistically proven.Open in a separate windowFig. 1.Isonitrile-mediated amidation; structure of OT.The potentialities of the isonitrile-mediated amidation method were foreshadowed via its application to the synthesis of cyclosporine (19). The power of the method was particularly well demonstrated in the context of our recent total synthesis of oxytocin (OT) (20), wherein isonitrile mediation was used in each of the peptide bond constructions, leading to the synthesis of the hormone in high yield and excellent purity. This nonapeptide is involved in a range of biological functions including parturition and lactation (21, 22). Signaling of OT to its receptor (OTR) is apparently an important factor in quality maintenance of various CNS functions (23). The ability to synthesize such modestly sized, but bio-impactful peptides in both native (wild-type) form, and as strategically modified variants, is one of the current missions of our laboratory, with the objective of possible applications to the very serious problem of autism (24–26).     

Contingency and entrenchment in protein evolution under purifying selection

The phenotypic effect of an allele at one genetic site may depend on alleles at other sites, a phenomenon known as epistasis. Epistasis can profoundly influence the process of evolution in populations and shape the patterns of protein divergence across species. Whereas epistasis between adaptive substitutions has been studied extensively, relatively little is known about epistasis under purifying selection. Here we use computational models of thermodynamic stability in a ligand-binding protein to explore the structure of epistasis in simulations of protein sequence evolution. Even though the predicted effects on stability of random mutations are almost completely additive, the mutations that fix under purifying selection are enriched for epistasis. In particular, the mutations that fix are contingent on previous substitutions: Although nearly neutral at their time of fixation, these mutations would be deleterious in the absence of preceding substitutions. Conversely, substitutions under purifying selection are subsequently entrenched by epistasis with later substitutions: They become increasingly deleterious to revert over time. Our results imply that, even under purifying selection, protein sequence evolution is often contingent on history and so it cannot be predicted by the phenotypic effects of mutations assayed in the ancestral background.Whether a heritable mutation is advantageous or deleterious to an organism often depends on the evolutionary history of the population. A mutation that is beneficial at the time of its introduction may confer its beneficial effect only in the presence of other potentiating or permissive mutations (1–9). Thus, the fate of a mutation arising in a population may be contingent on previous mutations (10–13). Conversely, once a mutation has fixed in a population, the mutation becomes part of the genetic background onto which subsequent modifications are introduced. Because the beneficial effects of the subsequent modifications may depend on the focal mutation, as time passes reversion of the focal mutation may become increasingly deleterious, leading to a type of evolutionary conservatism, or entrenchment (14–18).In the context of protein evolution, the effects of contingency and entrenchment are most easily studied by considering a sequence of single amino acid changes (19) that extends both forward and backward in time from some focal substitution. To assess the roles of contingency and entrenchment we can study the degree to which each focal substitution was facilitated by previous substitutions, and the degree to which the focal substitution influences the subsequent course of evolution (Fig. 1A).Open in a separate windowFig. 1.(A) A schematic model indicating how a focal substitution may be contingent on prior substitutions and may constrain future substitutions along an evolutionary trajectory, owing to epistasis. (B) A model of protein evolution under weak mutation and purifying selection for thermodynamic stability. Starting from the wild-type sequence of argT we propose 10 random 1-aa point mutations. For each of the proposed mutants we compute its predicted stability (ΔG) using FoldX, and its associated fitness. The fitness function is assumed to be either Gaussian or semi-Gaussian, with a maximum at the wild-type stability. One of the proposed mutants fixes in the population, based on its relative fixation probability under the Moran model with effective population size N_e. This process is iterated for 30 consecutive substitutions to produce an evolutionary trajectory. We simulate 100 replicate trajectories, each initiated at the wild-type argT sequence.Dependencies within a sequence of substitutions are closely connected to the concept of epistasis—that is, the idea that the phenotypic effect of a mutation at a particular genetic site may depend on the genetic background in which it arises (20–24). In the absence of epistasis, a mutation has the same effect regardless of its context and therefore regardless of any prior history or subsequent evolution. By contrast, in the presence of epistasis, each substitution may be contingent on the entire prior history of the protein, and it may constrain all subsequent evolution.The potential for epistasis to play an important role in evolution, including protein evolution, has not been overlooked by researchers (1, 8, 25–34), nor have the concepts of contingency (3, 4, 9, 12, 35–38) and, more recently, entrenchment (18, 39, 40). However, most studies have addressed the role of epistasis in the context of adaptive evolution (1–9, 27, 30, 31, 36, 38), whereas the consequences of epistasis under purifying selection have received less attention (18, 41–44). Indeed, although some more sophisticated models have been proposed (e.g., refs. 45–50), all commonly used phylogenetic models of long-term protein evolution assume that epistasis is absent so that sites evolve independently (51–56).Here we explore the relationships between epistasis, contingency, and entrenchment under long-term purifying selection on protein stability. Our analysis combines computational models for protein structures with population-genetic models for evolutionary dynamics. We use a force-field-based model, FoldX (57), to characterize the effects of point mutations on a protein’s stability and fitness. This approach allows us to simulate evolutionary trajectories of protein sequences under purifying selection, by the sequential fixation of nearly neutral mutations. We can then dissect the epistatic relationships between these substitutions by systematically inserting or reverting particular substitutions at various time points along the evolutionary trajectory.Our analysis considers epistasis both at the level of protein stability and at the level of fitness. Whereas empirical studies in diverse proteins have demonstrated that the stability effects of point mutations are typically additive across sites (58, 59), in this study we are specifically interested in epistasis for stability among the mutations that fix during evolution. Even if most random mutations are virtually additive in their effects on stability, the mutations that fix under purifying selection are highly nonrandom, and so there is reason to suspect that epistasis for stability may be enriched among such mutations. Moreover, because the mapping from stability to fitness is itself nonlinear (18, 26, 60, 61) and because selection is sensitive to selection coefficients as small as the inverse of the population size (62), even slight variation in the stability effects of mutations across different genetic backgrounds may be sufficient to influence the course of evolution.Using the computational approach summarized above, we will demonstrate that the nearly neutral mutations that fix under purifying selection are, indeed, often epistatic with each other for both stability and fitness. In particular, we find that each mutation that fixes is typically permitted to fix by the presence of preceding substitutions—that is, most substitutions would be too deleterious to fix were it not for epistasis with preceding substitutions. Conversely, we also find that mutations that fix typically become entrenched over time by epistasis—so that a substitution that was nearly neutral when it fixed becomes increasingly deleterious to revert as subsequent substitutions accumulate (18, 39). These results imply an important role for epistasis in shaping the course of sequence evolution in a protein under selection to maintain thermodynamic stability.     

The evolution of self-control

Cognition presents evolutionary research with one of its greatest challenges. Cognitive evolution has been explained at the proximate level by shifts in absolute and relative brain volume and at the ultimate level by differences in social and dietary complexity. However, no study has integrated the experimental and phylogenetic approach at the scale required to rigorously test these explanations. Instead, previous research has largely relied on various measures of brain size as proxies for cognitive abilities. We experimentally evaluated these major evolutionary explanations by quantitatively comparing the cognitive performance of 567 individuals representing 36 species on two problem-solving tasks measuring self-control. Phylogenetic analysis revealed that absolute brain volume best predicted performance across species and accounted for considerably more variance than brain volume controlling for body mass. This result corroborates recent advances in evolutionary neurobiology and illustrates the cognitive consequences of cortical reorganization through increases in brain volume. Within primates, dietary breadth but not social group size was a strong predictor of species differences in self-control. Our results implicate robust evolutionary relationships between dietary breadth, absolute brain volume, and self-control. These findings provide a significant first step toward quantifying the primate cognitive phenome and explaining the process of cognitive evolution.Since Darwin, understanding the evolution of cognition has been widely regarded as one of the greatest challenges for evolutionary research (1). Although researchers have identified surprising cognitive flexibility in a range of species (2–40) and potentially derived features of human psychology (41–61), we know much less about the major forces shaping cognitive evolution (62–71). With the notable exception of Bitterman’s landmark studies conducted several decades ago (63, 72–74), most research comparing cognition across species has been limited to small taxonomic samples (70, 75). With limited comparable experimental data on how cognition varies across species, previous research has largely relied on proxies for cognition (e.g., brain size) or metaanalyses when testing hypotheses about cognitive evolution (76–92). The lack of cognitive data collected with similar methods across large samples of species precludes meaningful species comparisons that can reveal the major forces shaping cognitive evolution across species, including humans (48, 70, 89, 93–98).To address these challenges we measured cognitive skills for self-control in 36 species of mammals and birds (Fig. 1 and Tables S1–S4) tested using the same experimental procedures, and evaluated the leading hypotheses for the neuroanatomical underpinnings and ecological drivers of variance in animal cognition. At the proximate level, both absolute (77, 99–107) and relative brain size (108–112) have been proposed as mechanisms supporting cognitive evolution. Evolutionary increases in brain size (both absolute and relative) and cortical reorganization are hallmarks of the human lineage and are believed to index commensurate changes in cognitive abilities (52, 105, 113–115). Further, given the high metabolic costs of brain tissue (116–121) and remarkable variance in brain size across species (108, 122), it is expected that the energetic costs of large brains are offset by the advantages of improved cognition. The cortical reorganization hypothesis suggests that selection for absolutely larger brains—and concomitant cortical reorganization—was the predominant mechanism supporting cognitive evolution (77, 91, 100–106, 120). In contrast, the encephalization hypothesis argues that an increase in brain volume relative to body size was of primary importance (108, 110, 111, 123). Both of these hypotheses have received support through analyses aggregating data from published studies of primate cognition and reports of “intelligent” behavior in nature—both of which correlate with measures of brain size (76, 77, 84, 92, 110, 124).Open in a separate windowFig. 1.A phylogeny of the species included in this study. Branch lengths are proportional to time except where long branches have been truncated by parallel diagonal lines (split between mammals and birds ∼292 Mya).With respect to selective pressures, both social and dietary complexities have been proposed as ultimate causes of cognitive evolution. The social intelligence hypothesis proposes that increased social complexity (frequently indexed by social group size) was the major selective pressure in primate cognitive evolution (6, 44, 48, 50, 87, 115, 120, 125–141). This hypothesis is supported by studies showing a positive correlation between a species’ typical group size and the neocortex ratio (80, 81, 85–87, 129, 142–145), cognitive differences between closely related species with different group sizes (130, 137, 146, 147), and evidence for cognitive convergence between highly social species (26, 31, 148–150). The foraging hypothesis posits that dietary complexity, indexed by field reports of dietary breadth and reliance on fruit (a spatiotemporally distributed resource), was the primary driver of primate cognitive evolution (151–154). This hypothesis is supported by studies linking diet quality and brain size in primates (79, 81, 86, 142, 155), and experimental studies documenting species differences in cognition that relate to feeding ecology (94, 156–166).Although each of these hypotheses has received empirical support, a comparison of the relative contributions of the different proximate and ultimate explanations requires (i) a cognitive dataset covering a large number of species tested using comparable experimental procedures; (ii) cognitive tasks that allow valid measurement across a range of species with differing morphology, perception, and temperament; (iii) a representative sample within each species to obtain accurate estimates of species-typical cognition; (iv) phylogenetic comparative methods appropriate for testing evolutionary hypotheses; and (v) unprecedented collaboration to collect these data from populations of animals around the world (70).Here, we present, to our knowledge, the first large-scale collaborative dataset and comparative analysis of this kind, focusing on the evolution of self-control. We chose to measure self-control—the ability to inhibit a prepotent but ultimately counterproductive behavior—because it is a crucial and well-studied component of executive function and is involved in diverse decision-making processes (167–169). For example, animals require self-control when avoiding feeding or mating in view of a higher-ranking individual, sharing food with kin, or searching for food in a new area rather than a previously rewarding foraging site. In humans, self-control has been linked to health, economic, social, and academic achievement, and is known to be heritable (170–172). In song sparrows, a study using one of the tasks reported here found a correlation between self-control and song repertoire size, a predictor of fitness in this species (173). In primates, performance on a series of nonsocial self-control control tasks was related to variability in social systems (174), illustrating the potential link between these skills and socioecology. Thus, tasks that quantify self-control are ideal for comparison across taxa given its robust behavioral correlates, heritable basis, and potential impact on reproductive success.In this study we tested subjects on two previously implemented self-control tasks. In the A-not-B task (27 species, n = 344), subjects were first familiarized with finding food in one location (container A) for three consecutive trials. In the test trial, subjects initially saw the food hidden in the same location (container A), but then moved to a new location (container B) before they were allowed to search (Movie S1). In the cylinder task (32 species, n = 439), subjects were first familiarized with finding a piece of food hidden inside an opaque cylinder. In the following 10 test trials, a transparent cylinder was substituted for the opaque cylinder. To successfully retrieve the food, subjects needed to inhibit the impulse to reach for the food directly (bumping into the cylinder) in favor of the detour response they had used during the familiarization phase (Movie S2).Thus, the test trials in both tasks required subjects to inhibit a prepotent motor response (searching in the previously rewarded location or reaching directly for the visible food), but the nature of the correct response varied between tasks. Specifically, in the A-not-B task subjects were required to inhibit the response that was previously successful (searching in location A) whereas in the cylinder task subjects were required to perform the same response as in familiarization trials (detour response), but in the context of novel task demands (visible food directly in front of the subject).     

A family of starch-active polysaccharide monooxygenases

The recently discovered fungal and bacterial polysaccharide monooxygenases (PMOs) are capable of oxidatively cleaving chitin, cellulose, and hemicelluloses that contain β(1→4) linkages between glucose or substituted glucose units. They are also known collectively as lytic PMOs, or LPMOs, and individually as AA9 (formerly GH61), AA10 (formerly CBM33), and AA11 enzymes. PMOs share several conserved features, including a monocopper center coordinated by a bidentate N-terminal histidine residue and another histidine ligand. A bioinformatic analysis using these conserved features suggested several potential new PMO families in the fungus Neurospora crassa that are likely to be active on novel substrates. Herein, we report on NCU08746 that contains a C-terminal starch-binding domain and an N-terminal domain of previously unknown function. Biochemical studies showed that NCU08746 requires copper, oxygen, and a source of electrons to oxidize the C1 position of glycosidic bonds in starch substrates, but not in cellulose or chitin. Starch contains α(1→4) and α(1→6) linkages and exhibits higher order structures compared with chitin and cellulose. Cellobiose dehydrogenase, the biological redox partner of cellulose-active PMOs, can serve as the electron donor for NCU08746. NCU08746 contains one copper atom per protein molecule, which is likely coordinated by two histidine ligands as shown by X-ray absorption spectroscopy and sequence analysis. Results indicate that NCU08746 and homologs are starch-active PMOs, supporting the existence of a PMO superfamily with a much broader range of substrates. Starch-active PMOs provide an expanded perspective on studies of starch metabolism and may have potential in the food and starch-based biofuel industries.Polysaccharide monooxygenases (PMOs) are enzymes secreted by a variety of fungal and bacterial species (1–5). They have recently been found to oxidatively degrade chitin (6–8) and cellulose (8–14). PMOs have been shown to oxidize either the C1 or C4 atom of the β(1→4) glycosidic bond on the surface of chitin (6, 7) or cellulose (10–12, 14), resulting in the cleavage of this bond and the creation of new chain ends that can be subsequently processed by hydrolytic chitinases and cellulases. Several fungal PMOs were shown to significantly enhance the degradation of cellulose by hydrolytic cellulases (9), indicating that these enzymes can be used in the conversion of plant biomass into biofuels and other renewable chemicals.There are three families of PMOs characterized thus far: fungal PMOs that oxidize cellulose (9–12) (also known as GH61 and AA9); bacterial PMOs that are active either on chitin (6, 8) or cellulose (8, 13) (also known as CBM33 and AA10); and fungal PMOs that oxidize chitin (AA11) (7). Sequence homology between these three families is very low. Nevertheless, the available structures of PMOs from all three families reveal a conserved fold, including an antiparallel β-sandwich core and a highly conserved monocopper active site on a flat protein surface (Fig. 1A) (2, 6, 7, 9, 10, 15–17). Two histidine residues in a motif termed the histidine brace coordinate the copper center. The N-terminal histidine ligand binds in a bidentate mode, and its imidazole ring is methylated at the N_ε position in fungal PMOs (Fig. 1A).Open in a separate windowFig. 1.(A) Representative overall and active site structures of fungal PMOs (PDB ID code 2YET) (10). (B) Structure of cellulose (18, 19). Chitin also contains β(1→4) linkages and has similar crystalline higher order structure to cellulose. (C) Model structure of amylopectin (23–25). Hydrogen bonds are shown with green dashed lines.Considering the conserved structural features, it is not surprising that the currently known PMOs act on substrates with similar structures. Cellulose and chitin contain long linear chains of β(1→4) linked glucose units and N-acetylglucosamine units, respectively (Fig. 1B). The polymer chains form extensive hydrogen bonding networks, which result in insoluble and very stable crystalline structures (18–21). PMOs are thought to bind to the substrate with their flat active site surface, which orients the copper center for selective oxidation at the C1 or C4 position (6, 16, 22). Some bacterial chitin-binding proteins are cellulose-active PMOs (8, 13, 14), further suggesting that the set of PMO substrates is restricted to β(1→4) linked polymers of glucose and glucose derivatives.Here, we report on the identification of new families of PMOs that contain several key features of previously characterized PMOs, but act on substrates different from cellulose or chitin. A member of one of these novel families of PMOs, NCU08746, was shown to oxidatively cleave amylose, amylopectin, and starch. We designate the NCU08746 family as starch-active PMOs. Both amylose and amylopectin contain linear chains of α(1→4) linked glucose, whereas the latter also contains α(1→6) glycosidic linkages at branch points in the otherwise α(1→4) linked polymer. Unlike cellulose and chitin, amylose and amylopectin do not form microcrystals; instead, they exist in disordered, single helical, and double helical forms (23–27) (see Fig. 1C for example). Starch exists partially in nanocrystalline form, but lacks the flat molecular surfaces as those found in chitin and cellulose. The discovery of starch-active PMOs shows that this oxidative mechanism of glycosidic bond cleavage is more widespread than initially expected.     

Exchange protein directly activated by cAMP plays a critical role in bacterial invasion during fatal rickettsioses

Rickettsiae are responsible for some of the most devastating human infections. A high infectivity and severe illness after inhalation make some rickettsiae bioterrorism threats. We report that deletion of the exchange protein directly activated by cAMP (Epac) gene, Epac1, in mice protects them from an ordinarily lethal dose of rickettsiae. Inhibition of Epac1 suppresses bacterial adhesion and invasion. Most importantly, pharmacological inhibition of Epac1 in vivo using an Epac-specific small-molecule inhibitor, ESI-09, completely recapitulates the Epac1 knockout phenotype. ESI-09 treatment dramatically decreases the morbidity and mortality associated with fatal spotted fever rickettsiosis. Our results demonstrate that Epac1-mediated signaling represents a mechanism for host–pathogen interactions and that Epac1 is a potential target for the prevention and treatment of fatal rickettsioses.Rickettsiae are responsible for some of the most devastating human infections (1–4). It has been forecasted that temperature increases attributable to global climate change will lead to more widespread distribution of rickettsioses (5). These tick-borne diseases are caused by obligately intracellular bacteria of the genus Rickettsia, including Rickettsia rickettsii, the causative agent of Rocky Mountain spotted fever (RMSF) in the United States and Latin America (2, 3), and Rickettsia conorii, the causative agent of Mediterranean spotted fever endemic to southern Europe, North Africa, and India (6). A high infectivity and severe illness after inhalation make some rickettsiae (including Rickettsia prowazekii, R. rickettsii, Rickettsia typhi, and R. conorii) bioterrorism threats (7). Although the majority of rickettsial infections can be controlled by appropriate broad-spectrum antibiotic therapy if diagnosed early, up to 20% of misdiagnosed or untreated (1, 3) and 5% of treated RMSF cases (8) result in a fatal outcome caused by acute disseminated vascular endothelial infection and damage (9). Fatality rates as high as 32% have been reported in hospitalized patients diagnosed with Mediterranean spotted fever (10). In addition, strains of R. prowazekii resistant to tetracycline and chloramphenicol have been developed in laboratories (11). Disseminated endothelial infection and endothelial barrier disruption with increased microvascular permeability are the central features of SFG rickettsioses (1, 2, 9). The molecular mechanisms involved in rickettsial infection remain incompletely elucidated (9, 12). A comprehensive understanding of rickettsial pathogenesis and the development of novel mechanism-based treatment are urgently needed.Living organisms use intricate signaling networks for sensing and responding to changes in the external environment. cAMP, a ubiquitous second messenger, is an important molecular switch that translates environmental signals into regulatory effects in cells (13). As such, a number of microbial pathogens have evolved a set of diverse virulence-enhancing strategies that exploit the cAMP-signaling pathways of their hosts (14). The intracellular functions of cAMP are predominantly mediated by the classic cAMP receptor, protein kinase A (PKA), and the more recently discovered exchange protein directly activated by cAMP (Epac) (15). Thus, far, two isoforms, Epac1 and Epac2, have been identified in humans (16, 17). Epac proteins function by responding to increased intracellular cAMP levels and activating the Ras superfamily small GTPases Ras-proximate 1 and 2 (Rap1 and Rap2). Accumulating evidence demonstrates that the cAMP/Epac1 signaling axis plays key regulatory roles in controlling various cellular functions in endothelial cells in vitro, including cell adhesion (18–21), exocytosis (22), tissue plasminogen activator expression (23), suppressor of cytokine signaling 3 (SOCS-3) induction (24–27), microtubule dynamics (28, 29), cell–cell junctions, and permeability and barrier functions (30–37). Considering the critical importance of endothelial cells in rickettsioses, we examined the functional roles of Epac1 in rickettsial pathogenesis in vivo, taking advantage of the recently generated Epac1 knockout mouse (38) and Epac-specific inhibitors (39, 40) generated from our laboratory. Our studies demonstrate that Epac1 plays a key role in rickettsial infection and represents a therapeutic target for fatal rickettsioses.     

Asymmetrical effects of unilateral right or left amygdala damage on auditory cortical processing of vocal emotions

We tested whether human amygdala lesions impair vocal processing in intact cortical networks. In two functional MRI experiments, patients with unilateral amygdala resection either listened to voices and nonvocal sounds or heard binaural vocalizations with attention directed toward or away from emotional information on one side. In experiment 1, all patients showed reduced activation to voices in the ipsilesional auditory cortex. In experiment 2, emotional voices evoked increased activity in both the auditory cortex and the intact amygdala for right-damaged patients, whereas no such effects were found for left-damaged amygdala patients. Furthermore, the left inferior frontal cortex was functionally connected with the intact amygdala in right-damaged patients, but only with homologous right frontal areas and not with the amygdala in left-damaged patients. Thus, unilateral amygdala damage leads to globally reduced ipsilesional cortical voice processing, but only left amygdala lesions are sufficient to suppress the enhanced auditory cortical processing of vocal emotions.Socially relevant and emotionally charged stimuli evoke increased activation in sensory cortices, both during the visual processing of emotional pictures or facial expressions (1, 2) and during the auditory processing of vocally expressed emotions (3–5). Such increases are assumed to be remotely driven by the amygdala, which is critically involved in decoding the emotional value of stimuli (6–8). Moreover, these effects seem to be predominantly (although not exclusively) mediated by ipsilateral anatomical (9) and functional connections between amygdala and sensory areas (10, 11).In line with this view, recent studies conducted in patients with amygdala lesions reported impairments in the recognition of facial expressions (12), emotional words (13), or vocal emotions (14, 15). Furthermore, studies in both human patients (16, 17) and monkeys (18) showed significant changes in visual cortical activations to facial expressions following lesions of the amygdala. These changes in cortical processing are assumed to be remotely driven by the impaired emotional processing in the amygdala (10, 16). Distant effects of amygdala damage have also been observed for visual stimuli in cats (19) and for auditory stimuli in rats (8). However, other results have challenged this view, with some studies reporting no impairment in recognition (20–25) or changes in cortical processing for emotional stimuli in patients with amygdala lesions (26). Notably, Edmiston et al. (26) observed normal visual increases in response to emotional scenes for patients with unilateral amygdala resection, arguing against a direct role for the amygdala in modulating activity in sensory cortical areas. However, in that study (26), such increases could be related to attentional effects driven by greater interest or complexity of emotional scenes (27, 28).Thus, evidence for impaired cortical responses to emotional stimuli after unilateral amygdala damage in humans remains inconsistent. In addition, unlike in rodents (8), to date, no study has investigated how the cortical processing of emotionally salient auditory stimuli might be affected by amygdala lesions in humans. Here, we tested for the first time, to our knowledge, whether unilateral amygdala damage in patients with left or right medial temporal lobe (MTL) lesions would modify auditory responses in intact cortical areas to voices and vocally expressed emotions. Previous studies consistently found differential activity in several subregions of auditory cortex in response to vocal emotions (29, 30), as well as in the amygdala (4, 5, 29, 31, 32), especially for angry voices (3–5). These auditory effects predominate in the superior temporal gyrus (STG) and superior temporal sulcus (STS), attributed to the processing of emotional valence in the amygdala (10) and presumably mediated by direct anatomical connections between the latter and auditory cortex (9, 33). Previous studies also consistently reported a response to emotional voices in the inferior frontal cortex, which may support higher level categorization processes (34) and thus constitutes an important component of the distributed network involved in detecting and decoding vocal emotions (29, 35).We therefore hypothesized that cortical processing of human vocalizations in general, and of vocal emotions in particular, might be impaired in patients with lesions to the amygdala. This impairment is thought to result from a reduced emotional decoding of affective vocal cues in the amygdala, which is generally sensitive to emotional cues in voices (4, 5, 29, 31, 36) and usually is assumed to enhance cortical processing remotely (10, 16, 19). We also hypothesized that left and right amygdala lesions might have different effects. Whereas a right MTL lesion may strongly impair the processing of facial expressions due to well-known hemispheric asymmetries in face processing (12, 17, 37), the left amygdala seems to be more strongly involved in the decoding of emotional cues expressed in speech (13) or speech-like material (3, 4, 29). In two experiments, we tested brain responses to human vocalizations in general (experiment 1) and to emotional vocalizations embedded in pseudolanguage (experiment 2) (Fig. 1) while 10 patients with unilateral left amygdala lesions and 10 patients with unilateral right amygdala lesions (SI Results, Fig. S1A, and Table S1) underwent functional MRI (fMRI) scanning. In experiment 2, emotional voices were presented in either the attended or unattended ear during a dichotic listening task (38). We expected, first, that vocalizations, as socially salient stimuli, would generally produce weaker cortical processing in interconnected regions due to unilateral amygdala damage (experiment 1) (39). Second, in keeping with predominant left amygdala activity in healthy individuals during the processing of vocal emotions, we expected more severe impairment in cortical processing of emotional cues in patients with left amygdala lesions compared with right amygdala lesions (experiment 2) (29, 40).Open in a separate windowFig. 1.Experiment 2 included three emotion conditions, with angry voices presented in the left or right ear or neither. On an trials, an angry voice was heard on the task-relevant side, whereas on na trials, an angry voice was heard on the task-irrelevant side. On nn trials, neutral voices were presented to each ear. (A) Examples show all three conditions when attention was focused on the right ear. The same trials were also performed when attention was focused on the left ear (not shown here). (B) RTs and error rates for gender decisions on the attended voice revealed a main effect for the factor condition, indicating increased RTs and error rates during the an condition, as indicated by the asterisks.     

INAUGURAL ARTICLE by a Recently Elected Academy Member:Lagging adaptation to warming climate in Arabidopsis thaliana

If climate change outpaces the rate of adaptive evolution within a site, populations previously well adapted to local conditions may decline or disappear, and banked seeds from those populations will be unsuitable for restoring them. However, if such adaptational lag has occurred, immigrants from historically warmer climates will outperform natives and may provide genetic potential for evolutionary rescue. We tested for lagging adaptation to warming climate using banked seeds of the annual weed Arabidopsis thaliana in common garden experiments in four sites across the species’ native European range: Valencia, Spain; Norwich, United Kingdom; Halle, Germany; and Oulu, Finland. Genotypes originating from geographic regions near the planting site had high relative fitness in each site, direct evidence for broad-scale geographic adaptation in this model species. However, genotypes originating in sites historically warmer than the planting site had higher average relative fitness than local genotypes in every site, especially at the northern range limit in Finland. This result suggests that local adaptive optima have shifted rapidly with recent warming across the species’ native range. Climatic optima also differed among seasonal germination cohorts within the Norwich site, suggesting that populations occurring where summer germination is common may have greater evolutionary potential to persist under future warming. If adaptational lag has occurred over just a few decades in banked seeds of an annual species, it may be an important consideration for managing longer-lived species, as well as for attempts to conserve threatened populations through ex situ preservation.Rapid climate change has already caused species range shifts and local extinctions (1) and is predicted to have greater future impacts (2). As the suitable climate space for a species shifts poleward (3), populations previously well adapted to the historical climate in a particular region may experience strong selection to adapt to rapidly warming local temperatures (4–10). Rapid evolutionary response to climate change has already been observed (11, 12), but it remains unclear whether evolutionary response can keep pace with rapidly changing local adaptive optima (6, 8, 13–15). If local adaptation is slower than the rate of climate change, the average fitness of local populations may decline over time (7, 14, 16, 17), possibly resulting in local extinctions and range collapse at the warmer margin. Where such lag exists, we expect that local seeds banked for conservation may no longer be well adapted to their sites of origin (18). However, such adaptational lag may be mitigated by migration or gene flow from populations in historically warmer sites if those populations are better adapted to current conditions in a site than local populations (8, 13, 19, 20). Although adaptational lag has been predicted (4–6, 8, 14, 15, 19, 21, 22), the distinctive signature of mismatch between local population performance and current climate optima has not yet been explicitly demonstrated in nature.Despite evidence for local adaptation in many organisms (23), there have been few explicit tests for the role of specific climate factors in shaping local fitness optima (4, 9, 13). Such tests require growing many genotypes from populations spanning a range of climates in common gardens across a species’ range to decouple climate of origin from geographic variation in other selective factors (4, 6, 14). If adaptation to local climate has occurred, then genotypes from climates similar to each planting site are expected to have high fitness in that site relative to genotypes from dissimilar climates (6). However, if local adaptive optima have shifted with rapid warming trends over the last 50 y, we expect that banked genotypes from historically warmer climates will have higher fitness within a site than banked genotypes of local origin (6, 21, 22).We tested for lagging adaptation to climate using Arabidopsis thaliana, a naturally inbreeding annual species that inhabits a broad climate space across its native Eurasian range (24). A. thaliana exhibits strong circumstantial evidence of climate adaptation, including geographic clines in ecologically important life-history traits (25–28) and in candidate genes associated with these traits (29, 30), as well as genome-wide associations of single nucleotide polymorphisms with climatic factors (31–34). To test explicitly for local adaptation to climate we measured the lifetime fitness of more than 230 accessions from banked seeds originating from a broad range of climates in replicated field experiments in four sites across the species’ native climate range (Fig. 1). We observed that genotypes originating in historically warmer climates outperformed local genotypes, particularly at the northern range limit.Open in a separate windowFig. 1.Map of common garden sites and sites of origin of the 241 native A. thaliana accessions represented in our experiments.     

Detection of genomic variations and DNA polymorphisms and impact on analysis of meiotic recombination and genetic mapping

DNA polymorphisms are important markers in genetic analyses and are increasingly detected by using genome resequencing. However, the presence of repetitive sequences and structural variants can lead to false positives in the identification of polymorphic alleles. Here, we describe an analysis strategy that minimizes false positives in allelic detection and present analyses of recently published resequencing data from Arabidopsis meiotic products and individual humans. Our analysis enables the accurate detection of sequencing errors, small insertions and deletions (indels), and structural variants, including large reciprocal indels and copy number variants, from comparisons between the resequenced and reference genomes. We offer an alternative interpretation of the sequencing data of meiotic products, including the number and type of recombination events, to illustrate the potential for mistakes in single-nucleotide polymorphism calling. Using these examples, we propose that the detection of DNA polymorphisms using resequencing data needs to account for nonallelic homologous sequences.DNA polymorphisms are ubiquitous genetic variations among individuals and include single nucleotide polymorphisms (SNPs), insertions and deletions (indels), and other larger rearrangements (1–3) (Fig. 1 A and B). They can have phenotypic consequences and also serve as molecular markers for genetic analyses, facilitating linkage and association studies of genetic diseases, and other traits in humans (4–6), animals, plants, (7–10) and other organisms. Using DNA polymorphisms for modern genetic applications requires low-error, high-throughput analytical strategies. Here, we illustrate the use of short-read next-generation sequencing (NGS) data to detect DNA polymorphisms in the context of whole-genome analysis of meiotic products.Open in a separate windowFig. 1.(A) SNPs and small indels between two ecotype genomes. (B) Possible types of SVs. Col genotypes are marked in blue and Ler in red. Arrows indicate DNA segments involved in SVs between the two ecotypes. (C) Meiotic recombination events including a CO and a GC (NCO). Centromeres are denoted by yellow dots.There are many methods for detecting SNPs (11–14) and structural variants (SVs) (15–25), including NGS, which can capture nearly all DNA polymorphisms (26–28). This approach has been widely used to analyze markers in crop species such as rice (29), genes associated with diseases (6, 26), and meiotic recombination in yeast and plants (30, 31). However, accurate identification of DNA polymorphisms can be challenging, in part because short-read sequencing data have limited information for inferring chromosomal context.Genomes usually contain repetitive sequences that can differ in copy number between individuals (26–28, 31); therefore, resequencing analyses must account for chromosomal context to avoid mistaking highly similar paralogous sequences for polymorphisms. Here, we use recently published datasets to describe several DNA sequence features that can be mistaken as allelic (32, 33) and describe a strategy for differentiating between repetitive sequences and polymorphic alleles. We illustrate the effectiveness of these analyses by examining the reported polymorphisms from the published datasets.Meiotic recombination is initiated by DNA double-strand breaks (DSBs) catalyzed by the topoisomerase-like SPORULATION 11 (SPO11). DSBs are repaired as either crossovers (COs) between chromosomes (Fig. 1C), or noncrossovers (NCOs). Both COs and NCOs can be accompanied by gene conversion (GC) events, which are the nonreciprocal transfer of sequence information due to the repair of heteroduplex DNA during meiotic recombination. Understanding the control of frequency and distribution of CO and NCO (including GC) events has important implications for human health (including cancer and aneuploidy), crop breeding, and the potential for use in genome engineering. COs can be detected relatively easily by using polymorphic markers in the flanking sequences, but NCO products can only be detected if they are accompanied by a GC event. Because GCs associated with NCO result in allelic changes at polymorphic sites without exchange of flanking sequences, they are more difficult to detect. Recent advances in DNA sequencing have made the analysis of meiotic NCOs more feasible (30–32, 34); however, SVs present a challenge in these analyses. We recommend a set of guidelines for detection of DNA polymorphisms by using genomic resequencing short-read datasets. These measures improve the accuracy of a wide range of analyses by using genomic resequencing, including estimation of COs, NCOs, and GCs.     

De novo production of the plant-derived alkaloid strictosidine in yeast

The monoterpene indole alkaloids are a large group of plant-derived specialized metabolites, many of which have valuable pharmaceutical or biological activity. There are ∼3,000 monoterpene indole alkaloids produced by thousands of plant species in numerous families. The diverse chemical structures found in this metabolite class originate from strictosidine, which is the last common biosynthetic intermediate for all monoterpene indole alkaloid enzymatic pathways. Reconstitution of biosynthetic pathways in a heterologous host is a promising strategy for rapid and inexpensive production of complex molecules that are found in plants. Here, we demonstrate how strictosidine can be produced de novo in a Saccharomyces cerevisiae host from 14 known monoterpene indole alkaloid pathway genes, along with an additional seven genes and three gene deletions that enhance secondary metabolism. This system provides an important resource for developing the production of more complex plant-derived alkaloids, engineering of nonnatural derivatives, identification of bottlenecks in monoterpene indole alkaloid biosynthesis, and discovery of new pathway genes in a convenient yeast host.Monoterpene indole alkaloids (MIAs) are a diverse family of complex nitrogen-containing plant-derived metabolites (1, 2). This metabolite class is found in thousands of plant species from the Apocynaceae, Loganiaceae, Rubiaceae, Icacinaceae, Nyssaceae, and Alangiaceae plant families (2, 3). Many MIAs and MIA derivatives have medicinal properties; for example, vinblastine, vincristine, and vinflunine are approved anticancer therapeutics (4, 5). These structurally complex compounds can be difficult to chemically synthesize (6, 7). Consequently, industrial production relies on extraction from the plant, but these compounds are often produced in small quantities as complex mixtures, making isolation challenging, laborious, and expensive (8–10). Reconstitution of plant pathways in microbial hosts is proving to be a promising approach to access plant-derived compounds as evidenced by the successful production of terpenes, flavonoids, and benzylisoquinoline alkaloids in microorganisms (11–19). Microbial hosts can also be used to construct hybrid biosynthetic pathways to generate modified natural products with potentially enhanced bioactivities (8, 20, 21). Across numerous plant species, strictosidine is believed to be the core scaffold from which all 3,000 known MIAs are derived (1, 2). Strictosidine undergoes a variety of redox reactions and rearrangements to form the thousands of compounds that comprise the MIA natural product family (Fig. 1) (1, 2). Due to the importance of strictosidine, the last common biosynthetic intermediate for all known MIAs, we chose to focus on heterologous production of this complex molecule (1). Therefore, strictosidine reconstitution represents the necessary first step for heterologous production of high-value MIAs.Open in a separate windowFig. 1.Strictosidine, the central intermediate in monoterpene indole alkaloid (MIA) biosynthesis, undergoes a series of reactions to produce over 3,000 known MIAs such as vincristine, quinine, and strychnine.     

From the Cover: Why the proportion of transmission during early-stage HIV infection does not predict the long-term impact of treatment on HIV incidence

Antiretroviral therapy (ART) reduces the infectiousness of HIV-infected persons, but only after testing, linkage to care, and successful viral suppression. Thus, a large proportion of HIV transmission during a period of high infectiousness in the first few months after infection (“early transmission”) is perceived as a threat to the impact of HIV “treatment-as-prevention” strategies. We created a mathematical model of a heterosexual HIV epidemic to investigate how the proportion of early transmission affects the impact of ART on reducing HIV incidence. The model includes stages of HIV infection, flexible sexual mixing, and changes in risk behavior over the epidemic. The model was calibrated to HIV prevalence data from South Africa using a Bayesian framework. Immediately after ART was introduced, more early transmission was associated with a smaller reduction in HIV incidence rate—consistent with the concern that a large amount of early transmission reduces the impact of treatment on incidence. However, the proportion of early transmission was not strongly related to the long-term reduction in incidence. This was because more early transmission resulted in a shorter generation time, in which case lower values for the basic reproductive number (R₀) are consistent with observed epidemic growth, and R₀ was negatively correlated with long-term intervention impact. The fraction of early transmission depends on biological factors, behavioral patterns, and epidemic stage and alone does not predict long-term intervention impacts. However, early transmission may be an important determinant in the outcome of short-term trials and evaluation of programs.Recent studies have confirmed that effective antiretroviral therapy (ART) reduces the transmission of HIV among stable heterosexual couples (1–3). This finding has generated interest in understanding the population-level impact of HIV treatment on reducing the rate of new HIV infections in generalized epidemic settings (4). Research, including mathematical modeling (5–10), implementation research (11), and major randomized controlled trials (12–14), are focused on how ART provision might be expanded strategically to maximize its public health benefits (15, 16).One concern is that if a large fraction of HIV transmission occurs shortly after a person becomes infected, before the person can be diagnosed and initiated on ART, this will limit the potential impact of HIV treatment on reducing HIV incidence (9, 17, 18). Data suggest that persons are more infectious during a short period of “early infection” after becoming infected with HIV (19–22), although there is debate about the extent, duration, and determinants of elevated infectiousness (18, 23). The amount of transmission that occurs also will depend on patterns of sexual behavior and sexual networks (17, 24–27). There have been estimates for the contribution of early infection to transmission from mathematical models (7, 17, 21, 24–26) and phylogenetic analyses (28–31), but these vary widely, from 5% to above 50% (23).In this study, we use a mathematical model to quantify how the proportion of transmission that comes from persons who have been infected recently affects the impact of treatment scale-up on HIV incidence. The model is calibrated to longitudinal HIV prevalence data from South Africa using a Bayesian framework. Thus, the model accounts for not only the early epidemic growth rate highlighted in previous research (5, 9, 18), but also the heterogeneity and sexual behavior change to explain the peak and decline in HIV incidence observed in sub-Saharan African HIV epidemics (32, 33).The model calibration allows uncertainty about factors that determine the amount of early transmission, including the relative infectiousness during early infection, heterogeneity in propensity for sexual risk behavior, assortativity in sexual partner selection, reduction in risk propensity over the life course, and population-wide reductions in risk behavior in response to the epidemic (32, 33). This results in multiple combinations of parameter values that are consistent with the observed epidemic and variation in the amount of early transmission. We simulated the impact of a treatment intervention and report how the proportion of early transmission correlates with the reduction in HIV incidence from the intervention over the short- and long-term.     

Self-assembly of alkynylplatinum(II) terpyridine amphiphiles into nanostructures via steric control and metal–metal interactions

A series of mono- and dinuclear alkynylplatinum(II) terpyridine complexes containing the hydrophilic oligo(para-phenylene ethynylene) with two 3,6,9-trioxadec-1-yloxy chains was designed and synthesized. The mononuclear alkynylplatinum(II) terpyridine complex was found to display a very strong tendency toward the formation of supramolecular structures. Interestingly, additional end-capping with another platinum(II) terpyridine moiety of various steric bulk at the terminal alkyne would lead to the formation of nanotubes or helical ribbons. These desirable nanostructures were found to be governed by the steric bulk on the platinum(II) terpyridine moieties, which modulates the directional metal−metal interactions and controls the formation of nanotubes or helical ribbons. Detailed analysis of temperature-dependent UV-visible absorption spectra of the nanostructured tubular aggregates also provided insights into the assembly mechanism and showed the role of metal−metal interactions in the cooperative supramolecular polymerization of the amphiphilic platinum(II) complexes.Square-planar d⁸ platinum(II) polypyridine complexes have long been known to exhibit intriguing spectroscopic and luminescence properties (1–54) as well as interesting solid-state polymorphism associated with metal−metal and π−π stacking interactions (1–14, 25). Earlier work by our group showed the first example, to our knowledge, of an alkynylplatinum(II) terpyridine system [Pt(tpy)(C ≡ CR)]⁺ that incorporates σ-donating and solubilizing alkynyl ligands together with the formation of Pt···Pt interactions to exhibit notable color changes and luminescence enhancements on solvent composition change (25) and polyelectrolyte addition (26). This approach has provided access to the alkynylplatinum(II) terpyridine and other related cyclometalated platinum(II) complexes, with functionalities that can self-assemble into metallogels (27–31), liquid crystals (32, 33), and other different molecular architectures, such as hairpin conformation (34), helices (35–38), nanostructures (39–45), and molecular tweezers (46, 47), as well as having a wide range of applications in molecular recognition (48–52), biomolecular labeling (48–52), and materials science (53, 54). Recently, metal-containing amphiphiles have also emerged as a building block for supramolecular architectures (42–44, 55–59). Their self-assembly has always been found to yield different molecular architectures with unprecedented complexity through the multiple noncovalent interactions on the introduction of external stimuli (42–44, 55–59).Helical architecture is one of the most exciting self-assembled morphologies because of the uniqueness for the functional and topological properties (60–69). Helical ribbons composed of amphiphiles, such as diacetylenic lipids, glutamates, and peptide-based amphiphiles, are often precursors for the growth of tubular structures on an increase in the width or the merging of the edges of ribbons (64, 65). Recently, the optimization of nanotube formation vs. helical nanostructures has aroused considerable interests and can be achieved through a fine interplay of the influence on the amphiphilic property of molecules (66), choice of counteranions (67, 68), or pH values of the media (69), which would govern the self-assembly of molecules into desirable aggregates of helical ribbons or nanotube scaffolds. However, a precise control of supramolecular morphology between helical ribbons and nanotubes remains challenging, particularly for the polycyclic aromatics in the field of molecular assembly (64–69). Oligo(para-phenylene ethynylene)s (OPEs) with solely π−π stacking interactions are well-recognized to self-assemble into supramolecular system of various nanostructures but rarely result in the formation of tubular scaffolds (70–73). In view of the rich photophysical properties of square-planar d⁸ platinum(II) systems and their propensity toward formation of directional Pt···Pt interactions in distinctive morphologies (27–31, 39–45), it is anticipated that such directional and noncovalent metal−metal interactions might be capable of directing or dictating molecular ordering and alignment to give desirable nanostructures of helical ribbons or nanotubes in a precise and controllable manner.Herein, we report the design and synthesis of mono- and dinuclear alkynylplatinum(II) terpyridine complexes containing hydrophilic OPEs with two 3,6,9-trioxadec-1-yloxy chains. The mononuclear alkynylplatinum(II) terpyridine complex with amphiphilic property is found to show a strong tendency toward the formation of supramolecular structures on diffusion of diethyl ether in dichloromethane or dimethyl sulfoxide (DMSO) solution. Interestingly, additional end-capping with another platinum(II) terpyridine moiety of various steric bulk at the terminal alkyne would result in nanotubes or helical ribbons in the self-assembly process. To the best of our knowledge, this finding represents the first example of the utilization of the steric bulk of the moieties, which modulates the formation of directional metal−metal interactions to precisely control the formation of nanotubes or helical ribbons in the self-assembly process. Application of the nucleation–elongation model into this assembly process by UV-visible (UV-vis) absorption spectroscopic studies has elucidated the nature of the molecular self-assembly, and more importantly, it has revealed the role of metal−metal interactions in the formation of these two types of nanostructures.     

Understanding the kinetic mechanism of RNA single base pair formation

RNA functions are intrinsically tied to folding kinetics. The most elementary step in RNA folding is the closing and opening of a base pair. Understanding this elementary rate process is the basis for RNA folding kinetics studies. Previous studies mostly focused on the unfolding of base pairs. Here, based on a hybrid approach, we investigate the folding process at level of single base pairing/stacking. The study, which integrates molecular dynamics simulation, kinetic Monte Carlo simulation, and master equation methods, uncovers two alternative dominant pathways: Starting from the unfolded state, the nucleotide backbone first folds to the native conformation, followed by subsequent adjustment of the base conformation. During the base conformational rearrangement, the backbone either retains the native conformation or switches to nonnative conformations in order to lower the kinetic barrier for base rearrangement. The method enables quantification of kinetic partitioning among the different pathways. Moreover, the simulation reveals several intriguing ion binding/dissociation signatures for the conformational changes. Our approach may be useful for developing a base pair opening/closing rate model.RNAs perform critical cellular functions at the level of gene expression and regulation (1–4). RNA functions are determined not only by RNA structure or structure motifs [e.g., tetraloop hairpins (5, 6)] but also by conformational distributions and dynamics and kinetics of conformational changes. For example, riboswitches can adopt different conformations in response to specific conditions of the cellular environment (7, 8). Understanding the kinetics, such as the rate and pathways for the conformational changes, is critical for deciphering the mechanism of RNA function (9–19). Extensive experimental and theoretical studies on RNA folding kinetics have provided significant insights into the kinetic mechanism of RNA functions (19–36). However, due to the complexity of the RNA folding energy landscape (37–46) and the limitations of experimental tools (47–55), many fundamental problems, including single base flipping and base pair formation and fraying, remain unresolved. These unsolved fundamental problems have hampered our ability to resolve other important issues, such as RNA hairpin and larger structure folding kinetics. Several key questions remain unanswered, such as whether the hairpin folding is rate-limited by the conformational search of the native base pairs, whose formation leads to fast downhill folding of the whole structure, or by the breaking of misfolded base pairs before refolding to the native structure (18, 19, 54–73).Motivated by the need to understand the basic steps of nucleic acids folding, Hagan et al. (74) performed forty-three 200-ps unfolding trajectories at 400 K and identified both on- and off-pathway intermediates and two dominant unfolding pathways for a terminal C-G base pair in a DNA duplex. In one of the pathways, base pairing and stacking interactions are broken concomitantly, whereas in the other pathway, base stacking is broken after base pairing is disrupted. Furthermore, the unfolding requires that the Cyt diffuse away from the pairing Gua to a distance such that the C-G hydrogen bond cannot reform easily. More recently, Colizzi and Bussi (75) performed molecular dynamics (MD) pulling simulations for an RNA duplex and construct free energy landscape from the pulling simulation. The simulation showed that the base pair opening reaction starts with the unbinding of the 5′-base, followed by the unbinding of the 3′-base (i.e., the 5′-base is less stable than the 3′-base). These previous unfolding simulations offered significant insights into the pathways and transition states. However, as shown below, several important issues remain.One intriguing problem is the rate model for base pairing. There are currently three main types of models. In the first type of model, the barrier ΔG+‡ for closing a base pair is dominated by the entropic cost ΔS for positioning the nucleotides to the base-paired configuration and the barrier ΔG−‡ for opening a base pair is the enthalpic cost ΔH for disrupting the hydrogen bonds and base stacking interactions (18, 59, 60). In the second type of model, ΔG+‡ is the net free energy change for base pairing ΔG = ΔH ? TΔS and ΔG−‡ is zero (76, 77). In the third type of model, ΔG±‡=±ΔG/2 is used (78). In addition to the above three main types, other models, such as more sophisticated hybrid rate models, have been proposed (29).In this paper, we report a hybrid method (see Fig. 1) to investigate the single base pairing process. In contrast to the previous simulations for temperature- or force-induced unfolding reactions, we directly model the folding process here (i.e., the base pair closing process). Specifically, we use MD simulations to identify the conformational clusters. Based on the network of the conformational clusters as a reduced conformational ensemble, we apply kinetic Monte Carlo (KMC) and master equation (ME) methods to elucidate the detailed roles of base pairing and stacking interactions, as well as the roles of water and ions (79–82). The study reveals previously unidentified kinetics pathways, misfolded states, and rate-limiting steps. A clear understanding of the microscopic details of the elementary kinetic move is a prerequisite for further rigorous study of large-scale RNA kinetic studies. The method described here may provide a feasible way to develop a rate model for the base pair/stack-based kinetic move set. Furthermore, the mechanism of RNA single base folding may provide useful insights into many biologically significant processes, such as nucleotide flipping (83) in helicases and base pair fraying (84) (as the possible first step for nucleic duplex melting in nucleic acid enzymatic processes).Open in a separate windowFig. 1.(A) Folding of a single nucleotide (G1, red) from the unfolded (Left) to the native folded (Right) state. (B) Exhaustive sampling for the (discrete) conformations of the G1 nucleotide (Right) through enumeration of the torsion angles (formed by the blue bonds). (C) Schematic plot shows the trajectories on the energy landscape (depicted with two reaction coordinates for clarity) explored by the MD simulations. The lines, open circles, and hexagons denote the trajectories; the initial states; and the (centroid structures of the) clusters, respectively. (D) Conformational network based on six clusters. (E) The rmsds to the different clusters provide information about the structural changes in a MD trajectory.     

Hippocampal molecular mechanisms involved in the enhancement of fear extinction caused by exposure to novelty

Exposure to a novel environment enhances the extinction of contextual fear. This has been explained by tagging of the hippocampal synapses used in extinction, followed by capture of proteins from the synapses that process novelty. The effect is blocked by the inhibition of hippocampal protein synthesis following the novelty or the extinction. Here, we show that it can also be blocked by the postextinction or postnovelty intrahippocampal infusion of the NMDA receptor antagonist 2-amino-5-phosphono pentanoic acid; the inhibitor of calcium/calmodulin-dependent protein kinase II (CaMKII), autocamtide-2–related inhibitory peptide; or the blocker of L-voltage–dependent calcium channels (L-VDCCs), nifedipine. Inhibition of proteasomal protein degradation by β-lactacystin has no effect of its own on extinction or on the influence of novelty thereon but blocks the inhibitory effects of all the other substances except that of rapamycin on extinction, suggesting that their action depends on concomitant synaptic protein turnover. Thus, the tagging-and-capture mechanism through which novelty enhances fear extinction involves more molecular processes than hitherto thought: NMDA receptors, L-VDCCs, CaMKII, and synaptic protein turnover.Frey and Morris (1, 2) and their collaborators (3–7) proposed a mechanism whereby relatively “weak” hippocampal long-term potentiation (LTP) or long-term depression (LTD) lasting only a few minutes can nevertheless “tag” the synapses involved with proteins synthesized ad hoc, so that other plasticity-related proteins (PRPs) produced at other sets of synapses by other LTPs or LTDs can be captured by the tagged synapses and strengthen their activity to “long” LTPs or LTDs lasting hours or days (8). LTDs and LTPs can “cross”-tag each other; that is, LTDs can enhance both LTDs and LTPs, and vice versa (6, 8). Because many learned behaviors rely on hippocampal LTP or LTD (7–9), among them the processing of novelty (9, 10) and the making of extinction (11–13), interactions between consecutive learnings can also be explained by the “tagging-and-capture” hypothesis (9, 10, 13), whose application to behavior became known as “behavioral tagging and capture” (5, 7, 9, 13). Typically, exposure to a novel environment [e.g., a nonanxiogenic 50 × 50 × 40-cm open field (OF) (5, 7, 9, 10, 14)] is interpolated before testing for another task, which becomes enhanced (4–10, 13). The usual reaction to novelty is orienting and exploration (14), followed by habituation of this response (14–16). Habituation is perhaps the simplest form of learning, and it consists of inhibition of the orienting/exploratory response (14, 16).We recently showed that the brief exposure of rats to a novel environment (the OF) within a limited time window enhances the extinction of contextual fear conditioning (CFC) through a mechanism of synaptic tagging and capture (13), which is a previously unidentified example of behavioral tagging of inhibitory learning. Fear extinction is most probably due to LTD in the hippocampus (11, 12), although the possibility that it may also involve LTP is not discarded (13). The enhancement of extinction by novelty probably relies on the habituation to the novel environment, which is also probably due to LTD (15, 16). The enhancement of extinction by the exposure to novelty depends on hippocampal gene expression and ribosomal protein synthesis following extinction training and on both ribosomal and nonribosomal protein synthesis caused by the novel experience (13). Nonribosomal protein synthesis that can be blocked by rapamycin is believed to be dendritic (13, 17), so it would be strategically located for tagging-and-capture processes, but it has not been studied in synaptic tagging to date (3–8) or in other forms of behavioral tagging (7–10). As occurs with the interactions between LTPs and/or LTDs (4), the enhancement of extinction by novelty relies on hippocampal but not amygdalar processes (13).Recent findings indicate that several hippocampal processes related to learning and memory, such as the reconsolidation of spatial learning, are highly dependent on NMDA glutamate receptors, calcium/calmodulin protein kinase II (CaMKII), and long-term voltage channel blockers (L-VDCCs), which, in turn, rely on the proteasomal degradation of proteins (18). Here, we study the effects of an NMDA blocker, 2-amino-5-phosphono pentanoic acid (AP5); the L-VDCC blocker nifedipine (Nife); a CaMKII inhibitor, the autocamtide-2–related inhibitory peptide (AIP); and the irreversible proteasome blocker β-lactacystin (12, 13) on the interaction between novelty and extinction (11). As will be seen, we found that both the setting up of tags by extinction and the presumable production of PRPs by the processing of novelty are dependent on NMDA receptors, CaMKII, and L-VDCCs. This endorses and expands the hypothesis that the novelty–extinction interaction relies on synaptic tagging and capture (13).     

Potentiation of cytotoxic chemotherapy by growth hormone-releasing hormone agonists

The dismal prognosis of malignant brain tumors drives the development of new treatment modalities. In view of the multiple activities of growth hormone-releasing hormone (GHRH), we hypothesized that pretreatment with a GHRH agonist, JI-34, might increase the susceptibility of U-87 MG glioblastoma multiforme (GBM) cells to subsequent treatment with the cytotoxic drug, doxorubicin (DOX). This concept was corroborated by our findings, in vivo, showing that the combination of the GHRH agonist, JI-34, and DOX inhibited the growth of GBM tumors, transplanted into nude mice, more than DOX alone. In vitro, the pretreatment of GBM cells with JI-34 potentiated inhibitory effects of DOX on cell proliferation, diminished cell size and viability, and promoted apoptotic processes, as shown by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide proliferation assay, ApoLive-Glo multiplex assay, and cell volumetric assay. Proteomic studies further revealed that the pretreatment with GHRH agonist evoked differentiation decreasing the expression of the neuroectodermal stem cell antigen, nestin, and up-regulating the glial maturation marker, GFAP. The GHRH agonist also reduced the release of humoral regulators of glial growth, such as FGF basic and TGFβ. Proteomic and gene-expression (RT-PCR) studies confirmed the strong proapoptotic activity (increase in p53, decrease in v-myc and Bcl-2) and anti-invasive potential (decrease in integrin α3) of the combination of GHRH agonist and DOX. These findings indicate that the GHRH agonists can potentiate the anticancer activity of the traditional chemotherapeutic drug, DOX, by multiple mechanisms including the induction of differentiation of cancer cells.Glioblastoma multiforme (GBM) is one of the most aggressive human cancers, and the afflicted patients inevitably succumb. The dismal outcome of this malignancy demands great efforts to find improved methods of treatment (1). Many compounds have been synthesized in our laboratory in the past few years that have proven to be effective against diverse malignant tumors (2–14). These are peptide analogs of hypothalamic hormones: luteinizing hormone-releasing hormone (LHRH), growth hormone-releasing hormone (GHRH), somatostatin, and analogs of other neuropeptides such as bombesin and gastrin-releasing peptide. The receptors for these peptides have been found to be widely distributed in the human body, including in many types of cancers (2–14). The regulatory functions of these hypothalamic hormones and other neuropeptides are not confined to the hypothalamo–hypophyseal system or, even more broadly, to the central nervous system (CNS). In particular, GHRH can induce the differentiation of ovarian granulosa cells and other cells in the reproductive system and function as a growth factor in various normal tissues, benign tumors, and malignancies (2–4, 6, 11, 14–18). Previously, we also reported that antagonistic cytototoxic derivatives of some of these neuropeptides are able to inhibit the growth of several malignant cell lines (2–14).Our earlier studies showed that treatment with antagonists of LHRH or GHRH rarely effects complete regression of glioblastoma-derived tumors (5, 7, 10, 11). Previous studies also suggested that growth factors such as EGF or agonistic analogs of LHRH serving as carriers for cytotoxic analogs and functioning as growth factors may sensitize cancer cells to cytotoxic treatments (10, 19) through the activation of maturation processes. We therefore hypothesized that pretreatment with one of our GHRH agonists, such as JI-34 (20), which has shown effects on growth and differentiation in other cell lines (17, 18, 21, 22), might decrease the pluripotency and the adaptability of GBM cells and thereby increase their susceptibility to cytotoxic treatment.In vivo, tumor cells were implanted into athymic nude mice, tumor growth was recorded weekly, and final tumor mass was measured upon autopsy. In vitro, proliferation assays were used for the determination of neoplastic proliferation and cell growth. Changes in stem (nestin) and maturation (GFAP) antigen expression was evaluated with Western blot studies in vivo and with immunocytochemistry in vitro. The production of glial growth factors (FGF basic, TGFβ) was verified by ELISA. Further, using the Human Cancer Pathway Finder real-time quantitative PCR, numerous genes that play a role in the development of cancer were evaluated. We placed particular emphasis on the measurement of apoptosis, using the ApoLive-Glo Multiplex Assay kit and by detection of the expression of the proapoptotic p53 protein. This overall approach permitted the evaluation of the effect of GHRH agonist, JI-34, on the response to chemotherapy with doxorubicin.