Single-molecule spectroscopy reveals how calmodulin activates NO synthase by controlling its conformational fluctuation dynamics

Mechanisms that regulate the nitric oxide synthase enzymes (NOS) are of interest in biology and medicine. Although NOS catalysis relies on domain motions, and is activated by calmodulin binding, the relationships are unclear. We used single-molecule fluorescence resonance energy transfer (FRET) spectroscopy to elucidate the conformational states distribution and associated conformational fluctuation dynamics of the two electron transfer domains in a FRET dye-labeled neuronal NOS reductase domain, and to understand how calmodulin affects the dynamics to regulate catalysis. We found that calmodulin alters NOS conformational behaviors in several ways: It changes the distance distribution between the NOS domains, shortens the lifetimes of the individual conformational states, and instills conformational discipline by greatly narrowing the distributions of the conformational states and fluctuation rates. This information was specifically obtainable only by single-molecule spectroscopic measurements, and reveals how calmodulin promotes catalysis by shaping the physical and temporal conformational behaviors of NOS.Although proteins adopt structures determined by their amino acid sequences, they are not static objects and fluctuate among ensembles of conformations (1). Transitions between these states can occur on a variety of length scales (Å to nm) and time scales (ps to s) and have been linked to functionally relevant phenomena such as allosteric signaling, enzyme catalysis, and protein–protein interactions (2–4). Indeed, protein conformational fluctuations and dynamics, often associated with static and dynamic inhomogeneity, are thought to play a crucial role in biomolecular functions (5–11). It is difficult to characterize such spatially and temporally inhomogeneous dynamics in bulk solution by an ensemble-averaged measurement, especially in proteins that undergo multiple-conformation transformations. In such cases, single-molecule spectroscopy is a powerful approach to analyze protein conformational states and dynamics under physiological conditions, and can provide a molecular-level perspective on how a protein’s structural dynamics link to its functional mechanisms (12–21).A case in point is the nitric oxide synthase (NOS) enzymes (22–24), whose nitric oxide (NO) biosynthesis involves electron transfer reactions that are associated with relatively large-scale movement (tens of Å) of the enzyme domains (Fig. 1A). During catalysis, NADPH-derived electrons first transfer into an FAD domain and an FMN domain in NOS that together comprise the NOS reductase domain (NOSr), and then transfer from the FMN domain to a heme group that is bound in a separate attached “oxygenase” domain, which then enables NO synthesis to begin (22, 25–27). The electron transfers into and out of the FMN domain are the key steps for catalysis, and they appear to rely on the FMN domain cycling between electron-accepting and electron-donating conformational states (28, 29) (Fig. 1B). In this model, the FMN domain is suggested to be highly dynamic and flexible due to a connecting hinge that allows it to alternate between its electron-accepting (FAD→FMN) or closed conformation and electron-donating (FMN→heme) or open conformation (Fig. 1 A and B) (28, 30–36). In the electron-accepting closed conformation, the FMN domain interacts with the NADPH/FAD domain (FNR domain) to receive electrons, whereas in the electron donating open conformation the FMN domain has moved away to expose the bound FMN cofactor so that it may transfer electrons to a protein acceptor like the NOS oxygenase domain, or to a generic protein acceptor like cytochrome c. In this way, the reductase domain structure cycles between closed and open conformations to deliver electrons, according to a conformational equilibrium that determines the movements and thus the electron flux capacity of the FMN domain (25, 28, 32, 34, 35, 37). A similar conformational switching mechanism is thought to enable electron transfer through the FMN domain in the related flavoproteins NADPH-cytochrome P450 reductase and methionine synthase reductase (38–42).Open in a separate windowFig. 1.(A) The nNOSr ribbon structure (from PDB: 1TLL) showing bound FAD (yellow) in FNR domain (green), FMN (orange) in FMN domain (yellow), connecting hinge (blue), and the Cy3–Cy5 label positions (pink) and distance (42 Å, dashed line). (B) Cartoon of an equilibrium between the FMN-closed and FMN-open states, with Cy dye label positions indicated. (C) Cytochrome c reductase activity of nNOSr proteins in their CaM-bound and CaM-free states. Color scheme of bar graphs: Black, WT nNOSr unlabeled; Red, Cys-lite (CL) nNOSr unlabeled; Blue, E827C/Q1268C CL nNOSr unlabeled; and Dark cyan, E827C/Q1268C CL nNOSr labeled.NOS enzymes also contain a calmodulin (CaM) binding domain that is located just before the N terminus of the FMN domain (Fig. 1B), and this provides an important layer of regulation (25, 27). CaM binding to NOS enzymes increases electron transfer from NADPH through the reductase domain and also triggers electron transfer from the FMN domain to the NOS heme as is required for NO synthesis (31, 32). The ability of CaM, or similar signaling proteins, to regulate electron transfer reactions in enzymes is unusual, and the mechanism is a topic of interest and intensive study. It has long been known that CaM binding alters NOSr structure such that, on average, it populates a more open conformation (43, 44). Recent equilibrium studies have detected a buildup of between two to four discreet conformational populations in NOS enzymes and in related flavoproteins, and in some cases, have also estimated the distances between the bound FAD and FMN cofactors in the different species (26, 36, 37, 39, 40), and furthermore, have confirmed that CaM shifts the NOS population distribution toward more open conformations (34, 36, 45). Although valuable, such ensemble-averaged results about conformational states cannot explain how electrons transfer through these enzymes, or how CaM increases the electron flux in NOS, because answering these questions requires a coordinate understanding of the dynamics of the conformational fluctuations. Indeed, computer modeling has indicated that a shift toward more open conformations as is induced by CaM binding to nNOS should, on its own, actually diminish electron flux through nNOS and through certain related flavoproteins (38). Despite its importance, measuring enzyme conformational fluctuation dynamics is highly challenging, and as far as we know, there have been no direct measures on the NOS enzymes or on related flavoproteins, nor studies on how CaM binding might influence the conformational fluctuation dynamics in NOS.To address this gap, we used single-molecule fluorescence energy resonance transfer (FRET) spectroscopy to characterize individual molecules of nNOSr that had been labeled at two specific positions with Cyanine 3 (Cy3) donor and Cyanine 5 (Cy5) acceptor dye molecules, regarding their conformational states distribution and the associated conformational fluctuation dynamics, which in turn enabled us to determine how CaM binding impacts both parameters. This work provides a unique perspective and a novel study of the NOS enzymes and within the broader flavoprotein family, which includes the mammalian enzymes methionine synthase reductase (MSR) and cytochrome P450 reductase (CPR), and reveals how CaM’s control of the conformational behaviors may regulate the electron transfer reactions of NOS catalysis.     

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).     

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.     

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.     

Human genome–guided identification of memory-modulating drugs

In the last decade there has been an exponential increase in knowledge about the genetic basis of complex human traits, including neuropsychiatric disorders. It is not clear, however, to what extent this knowledge can be used as a starting point for drug identification, one of the central hopes of the human genome project. The aim of the present study was to identify memory-modulating compounds through the use of human genetic information. We performed a multinational collaborative study, which included assessment of aversive memory—a trait central to posttraumatic stress disorder—and a gene-set analysis in healthy individuals. We identified 20 potential drug target genes in two genomewide-corrected gene sets: the neuroactive ligand–receptor interaction and the long-term depression gene set. In a subsequent double-blind, placebo-controlled study in healthy volunteers, we aimed at providing a proof of concept for the genome-guided identification of memory modulating compounds. Pharmacological intervention at the neuroactive ligand–receptor interaction gene set led to significant reduction of aversive memory. The findings demonstrate that genome information, along with appropriate data mining methodology, can be used as a starting point for the identification of memory-modulating compounds.Recent advances in human genetics have led to an unprecedented rate of discovery of genes related to complex human disease, including neuropsychiatric disorders (1–3). The human genome–based gain of knowledge is certainly expected to have a large impact on drug discovery in complex human disease (4–6). It is, however, still not clear to what extent this knowledge can be used as a starting point for the identification of druggable molecular pathways of complex traits (7), including mental disorders (8).Genomewide association studies (GWASs) using single-marker statistics have been very successful in identifying trait-associated single-gene loci (9). It is, however, widely accepted that single marker–based analyses have limited power to identify the genetic basis of a given trait, as for example, many loci will fail to reach stringent genomewide significance threshold, despite the fact that they may be genuinely associated with the trait. Triggered by statistical approaches for the analysis of gene expression, gene set–based analytical methods have recently become available. These methods aim at identifying biologically meaningful sets of genes associated with a certain trait, rather than focusing on a single GWAS gene locus (10). By taking into account prior biological knowledge, gene set–based approaches examine whether test statistics for a group of related genes have consistent deviation from chance (10). As shown recently in studies on autism (11), bipolar disorder (12, 13), attention deficit hyperactivity disorder (ADHD) (14), and schizophrenia (15), such approaches can convincingly identify convergent molecular pathways relevant to neuropsychiatry. Importantly, the identification of groups of functionally related genes is likely to facilitate drug discovery, because the most significant single loci from a GWAS might not be the best candidates for therapeutic intervention (7, 10).In the present study, we focused on emotionally aversive memory—a trait central to anxiety disorders such as posttraumatic stress disorder (PTSD) (16–23). Strong memory for emotionally arousing events can be seen as a primarily adaptive phenomenon, which helps us to remember vital information (e.g., dangerous situations). In case of an extremely aversive event, however, this mechanism can also lead to intrusive and persistent traumatic memories, thereby contributing to the development and symptoms of PTSD (18–22). Symptoms related to aversive memory include intrusive daytime recollections, traumatic nightmares, and flashbacks in which components of the event are relived. Aversive memory is a genetically complex trait as shown both in healthy subjects and in traumatized individuals (17, 23). Furthermore, we recently reported evidence suggesting a genetic link between the predisposition to build strong aversive memories and the risk for PTSD (16).Based on these observations, we developed a process (Fig. 1) aimed at identifying gene sets related to human aversive memory, followed by a pharmacological intervention study as proof-of-concept for the genome-guided identification of memory-modulating drugs.Open in a separate windowFig. 1.Drug discovery process.     

Tautomerism provides a molecular explanation for the mutagenic properties of the anti-HIV nucleoside 5-aza-5,6-dihydro-2′-deoxycytidine

Viral lethal mutagenesis is a strategy whereby the innate immune system or mutagenic pool nucleotides increase the error rate of viral replication above the error catastrophe limit. Lethal mutagenesis has been proposed as a mechanism for several antiviral compounds, including the drug candidate 5-aza-5,6-dihydro-2′-deoxycytidine (KP1212), which causes A-to-G and G-to-A mutations in the HIV genome, both in tissue culture and in HIV positive patients undergoing KP1212 monotherapy. This work explored the molecular mechanism(s) underlying the mutagenicity of KP1212, and specifically whether tautomerism, a previously proposed hypothesis, could explain the biological consequences of this nucleoside analog. Establishing tautomerism of nucleic acid bases under physiological conditions has been challenging because of the lack of sensitive methods. This study investigated tautomerism using an array of spectroscopic, theoretical, and chemical biology approaches. Variable temperature NMR and 2D infrared spectroscopic methods demonstrated that KP1212 existed as a broad ensemble of interconverting tautomers, among which enolic forms dominated. The mutagenic properties of KP1212 were determined empirically by in vitro and in vivo replication of a single-stranded vector containing a single KP1212. It was found that KP1212 paired with both A (10%) and G (90%), which is in accord with clinical observations. Moreover, this mutation frequency is sufficient for pushing a viral population over its error catastrophe limit, as observed before in cell culture studies. Finally, a model is proposed that correlates the mutagenicity of KP1212 with its tautomeric distribution in solution.Many viruses exhibit a high mutation rate when replicating their genomes, enabling quick adaptation to both changing cellular environments and therapeutics (1–5). Mammalian innate immune systems have developed a mechanism to exploit this high mutation rate against the virus; in a phenomenon termed “lethal mutagenesis,” (6–14) the immune system employs nucleic acid-modifying enzymes (e.g., APOBEC and ADAR) to increase the viral mutation rate sharply, stressing the functional gene product repertoire of the virus to the point that the viral population collapses (15–17). Several antiviral agents are proposed to work at least in part by a chemical version of lethal mutagenesis [e.g., ribavirin against hepatitis C virus (18–22), 5-hydroxy-2′-deoxycytidine against HIV (7), and T-705 against influenza viruses (23)]. When a sufficient number of these mutagenic nucleoside analogs is incorporated into viral genomes, the analogs increase the viral mutation rate above the error catastrophe limit, the rate above which no viable progeny are produced (6, 24–27). This work aimed to understand the molecular basis underlying the biological phenomenon of lethal mutagenesis induced by mutagenic nucleotides.The nucleoside analog 5-aza-5,6-dihydro-2′-deoxycytidine (KP1212) (Fig. 1A) is specifically designed to induce lethal mutagenesis in HIV (28–30). KP1212, the only anti-HIV drug candidate in clinical trials to use this mechanism, has been shown to increase the mutation rate of HIV both in cell culture and in isolates from humans undergoing monotherapy (28, 29). The mutagenic properties of KP1212 in cell culture reveal that it is likely to base pair promiscuously with A and G, and that the progressive acquisition of mutations (primarily A-to-G and G-to-A transitions) precedes population collapse (Fig. 1B) (29). These data are supported by biochemical experiments performed using purified polymerases that establish the ability of KP1212 to pair with either A or G, both when the modified base enters DNA from the nucleotide pool and when it acts as a template base (30). Understanding the chemical and structural basis of mutagenesis of this drug candidate is critical for both its future clinical progress and the development of new therapeutic agents that work by the principle of lethal mutagenesis.Open in a separate windowFig. 1.Schematic presentation of KP1212''s mutagenic effect on viruses. (A) KP1212 exists as an array of different tautomeric forms, whereas cytosine almost exclusively exists as one form, the canonical keto-amino tautomer. (B) The deoxynucleotide analog of KP1212 is incorporated by viral polymerases, causing G-to-A and A-to-G mutations during viral replication. KP1212 is a poor substrate for human polymerases, which provides selectivity in its action against the virus. The progressive acquisition of mutations in the viral genome leads to viral population collapse.Tautomerism of KP1212 leading to viral mutagenesis has been proposed by us and others (29, 30) to be the basis for the clinical activity of this drug candidate. There are, however, no direct data to support that view. Tautomerism as the basis of mutagenesis of natural bases has long been proposed (31–35), and substantiated in part by experimental evidence of minor tautomeric forms of both canonical bases (36–38) and certain base analogs (e.g., 5-hydroxy-2′-deoxycytidine) (39). In a search for a chemical rationale to explain the ambiguous pairing of KP1212 during replication, the present study revealed that the compound readily adopts multiple tautomeric forms, some of which were unexpected. Previously, spectroscopic methods (e.g., UV, Raman, NMR) have been used to study tautomerism of nucleobases (37, 39, 40). In the current work, we also used a battery of spectroscopic tools (1D, 2D, and variable temperature NMR; FTIR; and 2D IR) (41) to quantify and structurally characterize the array of tautomers exhibited by KP1212. Tautomer interconversion equilibria deconvoluted from NMR spectra provided data on the relative levels of tautomers in solution. In parallel with the spectroscopic studies, the qualitative and quantitative features of KP1212 mutagenesis were directly determined by inserting the KP1212 base into a single-stranded viral vector and measuring the intrinsic mutagenic properties of the base, both in vitro and in vivo. Finally, a model is proposed that correlates the mutagenic and clinical properties of KP1212 with its ability to exist as multiple tautomers.     

Multiaddressable molecular rectangles with reversible host–guest interactions: Modulation of pH-controlled guest release and capture

A series of multiaddressable platinum(II) molecular rectangles with different rigidities and cavity sizes has been synthesized by endcapping the U-shaped diplatinum(II) terpyridine moiety with various bis-alkynyl ligands. The studies of the host–guest association with various square planar platinum(II), palladium(II), and gold(III) complexes and the related low-dimensional gold(I) complexes, most of which are potential anticancer therapeutics, have been performed. Excellent guest confinement and selectivity of the rectangular architecture have been shown. Introduction of pH-responsive functionalities to the ligand backbone generates multifunctional molecular rectangles that exhibit reversible guest release and capture on the addition of acids and bases, indicating their potential in controlled therapeutics delivery on pH modulation. The reversible host–guest interactions are found to be strongly perturbed by metal–metal and π–π interactions and to a certain extent, electrostatic interactions, giving rise to various spectroscopic changes depending on the nature of the guest molecules. Their binding mode and thermodynamic parameters have been determined by 2D NMR and van’t Hoff analysis and supported by computational study.The study of metal–metal interactions has drawn enormous attention since the past two decades because of the intriguing spectroscopic and photophysical properties arising from the close proximity of the metal centers (1, 2). Square planar d⁸ platinum(II) complexes with coordination unsaturation are one of the important classes of metal complexes that have been extensively explored because of their capability to exhibit metal–metal interactions and display rich photophysical properties (3–26). Platinum(II) terpyridine complexes have been found to exhibit rich polymorphism in the solid state (16–20) owing to their square planar coordination geometry, which permits facile access to Pt(II)···Pt(II) interactions as well as π–π interactions between the chromophores. It was not until 2001 that the first successful synthesis of platinum(II) terpyridine alkynyl complexes, which possess enhanced solubility and luminescence compared with the chloro counterpart, was reported (16). Additional efforts have been devoted to the use of the system to respond to external stimuli, such as variation in solvent composition (17, 18), pH (19, 20), temperature (21, 22), addition of ionic (24–26), and polymeric species (27, 28), in which spectral changes induced by strong Pt(II)···Pt(II) and π−π interactions have been displayed.In the past few decades, enormous efforts have been devoted to the construction of molecular architectures by fusing the organic framework to the transition metal centers through self-assembly processes (29–57). There has been continuous interest in the construction of stimuli-responsive metallosupramolecular architectures with diverse sizes, shapes, and symmetries to rationalize the criteria for molecular recognition and impart them on unique areas of applications, such as stereoselective guest encapsulation and molecular transporting devices (45–65). Although such a variety of metal–organic macrocyclic architectures has been reported, those involving the use of noncovalent interactions other than those of hydrogen bonding, donor–acceptor, electrostatic, and hydrophobic–hydrophobic interactions as well as luminescence changes that depend on the nature of the guests, which would be attractive for chemo- and biosensing, have been rare and are rather underexplored. Examples of such systems that can exhibit reversible host–guest association are also limited.Since the discovery of anticancer properties of cisplatin in 1969 (58), the coordination chemistry and the development of related species with enhanced properties and reduced cytotoxicity have received enormous attention. Although the potency and cytotoxicity studies are important, the availability of the drugs and their transport and release to the site of action are equally important. Thus, the design of smart drug delivery systems has been an area of growing interest. The first phosphorescent molecular tweezers making use of the alkynylplatinum(II) terpyridine moiety have been reported by our group to show their host–guest interactions with transition metal complexes (57). However, the opened structures of the tweezers have limited their selectivity and functionality. To accomplish the controlled drug delivery functionalities, the first main strategy is to rigidify the molecular architecture of the host from tweezers to a rectangle, so that the guest molecules would be better accommodated within the cavity, which may lead to a more selective encapsulation of guests within a definite size and steric environment. The possibility of introducing responsive functionalities into the molecular rectangles, which may serve as models for the study of on-demand controlled guest capture and release systems, has also been explored. pH-sensitive pyridine moieties have, therefore, been incorporated into the backbone of the rectangle to modulate the reversible host–guest interaction within the constrained rectangle environment on protonation/deprotonation of the pyridine nitrogen atom to achieve multiaddressable functions that would not have been readily achievable with the molecular tweezers structure. Additionally, the use of various platinum and gold complexes as guest molecules, which have been shown to display anticancer therapeutic behavior (58–65), may lead to the design of a smart multiaddressable molecular rectangle system that could capture and release specific guest molecules under different pH conditions to achieve proof-of-principle on-demand controlled drug delivery. Herein, the design and synthesis of a series of alkynylplatinum(II) terpyridine molecular rectangles (Fig. 1) with different geometries, topologies and electronic properties are reported. Moreover, the encapsulation of various guest molecules is also investigated in detail to provide a proof-of-principle model for the design of artificial drug delivery systems with the modulation of drug release by pH.Open in a separate windowFig. 1.Molecular structures of rectangles 1−4.     

Shape control and compartmentalization in active colloidal cells

Small autonomous machines like biological cells or soft robots can convert energy input into control of function and form. It is desired that this behavior emerges spontaneously and can be easily switched over time. For this purpose we introduce an active matter system that is loosely inspired by biology and which we term an active colloidal cell. The active colloidal cell consists of a boundary and a fluid interior, both of which are built from identical rotating spinners whose activity creates convective flows. Similarly to biological cell motility, which is driven by cytoskeletal components spread throughout the entire volume of the cell, active colloidal cells are characterized by highly distributed energy conversion. We demonstrate that we can control the shape of the active colloidal cell and drive compartmentalization by varying the details of the boundary (hard vs. flexible) and the character of the spinners (passive vs. active). We report buckling of the boundary controlled by the pattern of boundary activity, as well as formation of core–shell and inverted Janus phase-separated configurations within the active cell interior. As the cell size is increased, the inverted Janus configuration spontaneously breaks its mirror symmetry. The result is a bubble–crescent configuration, which alternates between two degenerate states over time and exhibits collective migration of the fluid along the boundary. Our results are obtained using microscopic, non–momentum-conserving Langevin dynamics simulations and verified via a phase-field continuum model coupled to a Navier–Stokes equation.Active matter describes particulate systems with the characteristic that each “particle” (agent) converts energy into motion (1, 2). Active matter covers a range of length scales that include molecular motors in the cytoskeleton (3–5), swimming bacteria (6–8), driven colloids (9, 10), flocks of birds and fish (11–14), and people and vehicles in motion (15). Over the last decade, studies of active matter have demonstrated behavior not seen in equilibrium systems, including giant number fluctuations (16, 17), emergent attraction and superdiffusion (18–20), clustering (21, 22), swarming (23–27), and self-assembled motifs (28, 29). These systems provide interesting theoretical and engineering challenges as well as opportunities to explore and target novel behaviors that proceed outside of thermodynamic equilibrium.Of particular interest are systems found in nature or inspired by natural phenomena. Biological systems usually operate in confined regions of space––think of intracellular space, interfaces and membranes, and the crowding of cells near surfaces. The role of hydrodynamics in confinement has been studied for biological swimmers, such as bacteria and sperm, showing accumulation at the walls (30–32) and upstream swimming along surfaces (33) or in a spiral vortex (34–36). Attraction to walls has also been reported in the absence of hydrodynamics for disks (37, 38), spheres (39), and dumbbell swimmers (40). But, whereas these examples study the behavior under the influence of hard boundaries, biological swimmers typically interact with soft boundaries, such as membranes and biofilms. Another design variable is the possibility that the boundary itself is active, as in the surface of a bacterium covered with flagellae or, as demonstrated recently, active nematic vesicles (41).In this work, we propose and investigate an active matter system under flexible, active confinement. We call this system an active colloidal cell. Our realization of an active colloidal cell consists of independent particles, called spinners (42), that translate and rotate in two dimensions and are constrained within a finite area by a flexible boundary that is also built from spinners. Each spinner has a gear-like geometry, which consists of a large central disk and four smaller satellite disks (Fig. 1A). Similar gear-shaped rigid aggregates of self-propelled particles have been formed experimentally (43). Spinners are freely mobile in the cell interior. On the cellular boundary, spinners are connected to one another by a flexible chain of beads attached by finitely extensible springs. Both the interior and the boundary spinners can be subject to a clockwise or counterclockwise driving torque, which makes them active.Open in a separate windowFig. 1.Schematic of the confined spinner models. (A) The active colloidal cell is made up of spinners driven counterclockwise (blue) or clockwise (yellow). Boundary spinners are connected by a flexible bead–spring chain (gray). We compare the behavior of a continuum model (B) to a microscopic model (C). The compartmentalization of interior spinners is visualized by coloring the Voronoi tessellation in the microscopic model.Rotationally driven particles can synchronize and self-organize (44, 45) in the absence (42) and in the presence (46–48) of hydrodynamic interactions. Crystallization has recently been observed in rotating magnetic Janus colloids (49) and fast-moving bacteria (50). Spinners in the interior of the cell resemble molecular motors that push themselves forward on their neighbors and, thus, sustain convective dynamics. The effect of the boundary spinners is similar to that found in the cilia of living tissues, which stir nearby fluid. Our results demonstrate that a natural consequence of the activity present in the colloidal cell is control over both its external shape and internal structure. We report compartmentalization into regions of clockwise and counterclockwise spinners––a behavior which is affected by, and can be controlled via, properties of the enclosing boundary configuration as previously suggested (51). Transitions in the internal structure of the colloidal cell occur as its radius increases, and as the composition of the interior spinners and the patterning of the boundary are varied.A previous study of spinners in bulk (42) showed phase separation into clockwise- and counterclockwise domains. Cates and collaborators (6, 52, 53) have suggested that phase separation is a generic consequence of local energy input in an otherwise equilibrium system. Here and in the study of bulk spinners we demonstrate phase separation due to local rotational, rather than translational, energy input. We obtain our results using a particulate, microscopic model (Fig. 1C) as well as a continuum model (Fig. 1B). This allows us to conclude that the phenomena we observe are robust with respect to details of the model.In this study we use two models to study the behavior of an active colloidal cell, illustrated in Fig. 1. The microscopic model describes spinners as individual particles and simulates their motion using Langevin dynamics. It resolves the behavior of individual spinners but does not include hydrodynamic effects. In contrast, the continuum model describes the spinner system as a viscous binary fluid, which is governed by an incompressible Navier–Stokes equation coupled to a Cahn–Hilliard equation. Both models are described in detail in Materials and Methods below. Note that the microscopic model was introduced in earlier work using Brownian dynamics (42) and is extended here to include boundaries.     

Decision-related pupil dilation reflects upcoming choice and individual bias

A number of studies have shown that pupil size increases transiently during effortful decisions. These decision-related changes in pupil size are mediated by central neuromodulatory systems, which also influence the internal state of brain regions engaged in decision making. It has been proposed that pupil-linked neuromodulatory systems are activated by the termination of decision processes, and, consequently, that these systems primarily affect the postdecisional brain state. Here, we present pupil results that run contrary to this proposal, suggesting an important intradecisional role. We measured pupil size while subjects formed protracted decisions about the presence or absence (“yes” vs. “no”) of a visual contrast signal embedded in dynamic noise. Linear systems analysis revealed that the pupil was significantly driven by a sustained input throughout the course of the decision formation. This sustained component was larger than the transient component during the final choice (indicated by button press). The overall amplitude of pupil dilation during decision formation was bigger before yes than no choices, irrespective of the physical presence of the target signal. Remarkably, the magnitude of this pupil choice effect (yes > no) reflected the individual criterion: it was strongest in conservative subjects choosing yes against their bias. We conclude that the central neuromodulatory systems controlling pupil size are continuously engaged during decision formation in a way that reveals how the upcoming choice relates to the decision maker’s attitude. Changes in brain state seem to interact with biased decision making in the face of uncertainty.Changes in pupil size at constant luminance have long been used as a marker of central autonomic processes linked to cognition (1–4). Many studies over the past decades reported that the pupil dilates while subjects engage in demanding perceptual, cognitive, or economic decision tasks (1–3, 5–17). This decision-related pupil dilation has commonly been linked to the final choice terminating the decision process (6, 14, 16) and the consolidation of the committed decision (6, 16).Changes in pupil size are also linked to changes in brain state. It has been proposed that the decision-related pupil dilation tracks the activity of certain neuromodulatory systems of the brainstem—in particular, the noradrenergic locus coeruleus (5, 7–9, 18) and, possibly, the cholinergic basal forebrain (19) systems. These neuromodulatory systems also activate briefly (“phasically”) during perceptual decisions, such as visual target detection (5, 20–24), likely mediated via feedback connections from the prefrontal cortex (5, 25). The modulatory neurotransmitters released from the projections of these brainstem systems, in turn, shape the internal state of cortical networks, for instance, by boosting the gain of neural interactions (5, 7, 26). Thus, these brainstem systems might also shape decision computations in cortical networks—provided that they are activated already during decision formation. If so, these systems might affect the decision process, over and above shortening the time to respond. For instance, they might govern the decision maker’s ability to overcome his or her intrinsic bias.Here, we addressed these issues noninvasively in humans by linking decision-related pupil dilation to the time course, outcome, and bias of a protracted perceptual decision process. Many perceptual decisions are not transient events but evolve gradually over several hundreds of milliseconds, due to the slow accumulation of noisy sensory information (27–33). Further, perceptual decisions are, like economic decisions (34), prone to strong biases that are not due to external asymmetries in the magnitude or probability of payoffs for certain choices. In particular “yes” vs. “no” detection decisions depend on the idiosyncratic (liberal or conservative) attitude of the decision maker with respect to saying “yes” or “no” (35, 36).We thus measured pupil size in subjects performing a challenging yes–no visual contrast detection task at constant luminance (Fig. 1A). A general linear model (GLM) (37) allowed us to disentangle different temporal components of the neural input to the sluggish system controlling pupil size. This approach revealed that decision-related pupil dilation was not only driven by subjects’ final choice and the concomitant motor response, but also by a (stronger) sustained component throughout the preceding decision process. Further, the dilation amplitude was bigger for yes than for no choices. This pupil choice effect was due to the conservative subjects who decided yes against their bias. Taken together, our findings point to an intricate interplay between changes in internal brain state and biased decision making in the face of uncertainty.Open in a separate windowFig. 1.Task and behavioral results. (A) Sequence of events during a single trial. Dynamic noise is continuously present in a circular aperture around fixation. During the decision interval (onset cued by a tone), the subject searches for a faint grating signal superimposed onto the noise and indicates the yes or no choice by button press. The signal is shown at high contrast for illustration purposes only. In the actual experiment, its contrast was titrated to each individual’s detection threshold. (B) Stimulus types during the decision interval and possible choices of the subject, yielding the four trial categories of signal-detection theory. (C) Distribution of trial types, pooled across all subjects. (D) Reaction-time distributions for each trial type, pooled across all subjects. (E) Normalized reaction times, sorted by trial type and averaged across the group. RT, reaction time. Error bars, SEM.     

Robust efficiency and actuator saturation explain healthy heart rate control and variability

The correlation of healthy states with heart rate variability (HRV) using time series analyses is well documented. Whereas these studies note the accepted proximal role of autonomic nervous system balance in HRV patterns, the responsible deeper physiological, clinically relevant mechanisms have not been fully explained. Using mathematical tools from control theory, we combine mechanistic models of basic physiology with experimental exercise data from healthy human subjects to explain causal relationships among states of stress vs. health, HR control, and HRV, and more importantly, the physiologic requirements and constraints underlying these relationships. Nonlinear dynamics play an important explanatory role––most fundamentally in the actuator saturations arising from unavoidable tradeoffs in robust homeostasis and metabolic efficiency. These results are grounded in domain-specific mechanisms, tradeoffs, and constraints, but they also illustrate important, universal properties of complex systems. We show that the study of complex biological phenomena like HRV requires a framework which facilitates inclusion of diverse domain specifics (e.g., due to physiology, evolution, and measurement technology) in addition to general theories of efficiency, robustness, feedback, dynamics, and supporting mathematical tools.Biological systems display a variety of well-known rhythms in physiological signals (1–6), with particular patterns of variability associated with a healthy state (2–6). Decades of research demonstrate that heart rate (HR) in healthy humans has high variability, and loss of this high HR variability (HRV) is correlated with adverse states such as stress, fatigue, physiologic senescence, or disease (6–13). The dominant approach to analysis of HRV has been to focus on statistics and patterns in HR time series that have been interpreted as fractal, chaotic, scale-free, critical, etc. (6–17). The appeal of time series analysis is understandable as it puts HRV in the context of a broad and popular approach to complex systems (5, 18), all while requiring minimal attention to domain-specific (e.g., physiological) details. However, despite intense research activity in this area, there is limited consensus regarding causation or mechanism and minimal clinical application of the observed phenomena (10). This paper takes a completely different approach, aiming for more fundamental rigor (19–24) and methods that have the potential for clinical relevance. Here we use and model data from experimental studies of exercising healthy athletes, to add simple physiological explanations for the largest source of HRV and its changes during exercise. We also present methods that can be used to systematically pursue further explanations about HRV that can generalize to less healthy subjects.Fig. 1 shows the type of HR data analyzed, collected from healthy young athletes (n = 5). The data display responses to changes in muscle work rate on a stationary bicycle during mostly aerobic exercise. Fig. 1A shows three separate exercise sessions with identical workload fluctuations about three different means. With proper sleep, hydration, nutrition, and prevention from overheating, trained athletes can maintain the highest workload in Fig. 1 for hours and the lower and middle levels almost indefinitely. This ability requires robust efficiency: High workloads are sustained while robustly maintaining metabolic homeostasis, a particularly challenging goal in the case of the relatively large, metabolically demanding, and fragile human brain.Open in a separate windowFig. 1.HR responses to simple changes in muscle work rate on a stationary bicycle: Each experimental subject performed separate stationary cycle exercises of ∼10 min for each workload profile, with different means but nearly identical square wave fluctuations around the mean. A typical result is shown from subject 1 for three workload profiles with time on the horizontal axis (zoomed in to focus on a 6-min window). (A) HR (red) and workload (blue); linear local piecewise static fits (black) with different parameters for each exercise. The workload units (most strenuous exercise on top of graph) are shifted and scaled so that the blue curves are also the best global linear fit. (B) Corresponding dynamics fits, either local piecewise linear (black) or global linear (blue). Note that, on all time scales, mean HR increases and variability (HRV) goes down with the increasing workload. Breathing was spontaneous (not controlled).Whereas mean HR in Fig. 1A increases monotonically with workloads, both slow and fast fluctuations (i.e., HRV) in HR are saturating nonlinear functions of workloads, meaning that both high- and low-frequency HRV component goes down. Results from all subjects showed qualitatively similar nonlinearities (SI Appendix). We will argue that this saturating nonlinearity is the simplest and most fundamental example of change in HRV in response to stressors (11, 12, 25) [exercise in the experimental case, but in general also fatigue, dehydration, trauma, infection, even fear and anxiety (6–9, 11, 12, 25)].Physiologists have correlated HRV and autonomic tone (7, 11, 12, 14), and the (im)balance between sympathetic stimulation and parasympathetic withdrawal (12, 26–28). The alternation in autonomic control of HR (more sympathetic and less parasympathetic tone during exercise) serves as an obvious proximate cause for how the HRV changes as shown in Fig. 1, but the ultimate question remains as to why the system is implemented this way. It could be an evolutionary accident, or could follow from hard physiologic tradeoff requirements on cardiovascular control, as work in other systems suggests (1). Here, the explanation of HRV similarly involves hard physiological tradeoffs in robust efficiency and employs the mathematical tools necessary to make this explanation rigorous in the context of large measurement and modeling uncertainties.     

Aptamer photoregulation in vivo

The in vivo application of aptamers as therapeutics could be improved by enhancing target-specific accumulation while minimizing off-target uptake. We designed a light-triggered system that permits spatiotemporal regulation of aptamer activity in vitro and in vivo. Cell binding by the aptamer was prevented by hybridizing the aptamer to a photo-labile complementary oligonucleotide. Upon irradiation at the tumor site, the aptamer was liberated, leading to prolonged intratumoral retention. The relative distribution of the aptamer to the liver and kidney was also significantly decreased, compared to that of the free aptamer.Aptamers are single-stranded nucleic acids that have emerged as a promising class of therapeutics owing to their relative ease of synthesis and high affinity and selectivity toward a range of targets including small molecules, proteins, viral particles, and living cells (1–6). Aptamers can fold into well-defined conformations and are more resistant to enzymatic degradation than other oligonucleotides (7–9). Aptamers have been suggested for imaging applications because their relatively small size and molecular mass (∼10 kDa) allow fast tissue penetration and clearance from blood (10, 11). The same characteristics make aptamers promising for effective delivery of diagnostic and therapeutic agents to tissues or organs. However, nonspecific accumulation of aptamers in normal tissues is undesirable (12–15) because it diminishes the proportion of aptamer that targets the desired tissue. This can adversely affect the therapeutic index of the aptamer; this may be particularly true if the aptamer is conjugated to a drug or drug delivery device. Moreover, aptamers themselves can have nonspecific toxic effects (16, 17). Ideally aptamers would achieve a high concentration in a pathological tissue of interest while maintaining low levels elsewhere. The activity of aptamers can be modulated in vivo by binding to polymers or complementary oligonucleotide sequences (18–20), but spatiotemporal regulation of aptamer activity in vivo has not been achieved, whereby activity would be enhanced in target tissues and not others. Here, we report a strategy to provide light-triggered control of aptamer function and distribution in vivo.Light is an excellent means of providing external spatiotemporal control of biological systems (21–26). Many strategies have been developed to incorporate photosensitive groups in nucleotides that can control cellular function or affect biological pathways or gene expression by light (24–29). Of particular interest, such approaches can be used to provide spatiotemporal control of gene activation (24). Here we hypothesized that light triggering can be used to achieve spatiotemporal control of binding of an aptamer injected systemically to its target tissue in vivo, which would have implications for control of delivery of therapeutic aptamers and/or conjugated drugs or drug delivery systems. We designed a photo-triggerable system whereby the aptamer of interest is inactivated by hybridization to a photo-labile complementary oligonucleotide. Upon irradiation, the complementary sequence breaks down, releasing the functional aptamer (Fig. 1). The aptamer of interest is the single-stranded DNA 26-mer aptamer AS1411 (A₁₄₁₁; sequence: 5′-GGT GGT GGT GGT TGT GGT GGT GGT GG-3′) that binds with high affinity and selectivity to nucleolin (30–32), which is overexpressed on the cell membrane of several types of cancer cells, including the 4T1 breast cancer cells used here (33–35). A₁₄₁₁ has been used for cancer targeting in vitro and in vivo (35, 36). A complementary photo-triggerable inhibitory oligonucleotide (OliP) was designed [sequence: 5′-CCA CCA//CCA CCA//CAA CCA C-3′, where // indicates photo-labile 1-(2-nitrophenyl)ethyl bonds (37); Scheme S1].Open in a separate windowFig. 1.The AS1411 aptamer (A₁₄₁₁, red DNA strand) is hybridized to a complementary oligonucleotide (OliP, green DNA strand) containing photo-cleavable bonds (black dots). The hybridized complex cannot bind cells. The A₁₄₁₁ can be released by light-triggered breakage of the OliP, allowing binding to cell surfaces.     

Aggregation propensities of superoxide dismutase G93 hotspot mutants mirror ALS clinical phenotypes

Protein framework alterations in heritable Cu, Zn superoxide dismutase (SOD) mutants cause misassembly and aggregation in cells affected by the motor neuron disease ALS. However, the mechanistic relationship between superoxide dismutase 1 (SOD1) mutations and human disease is controversial, with many hypotheses postulated for the propensity of specific SOD mutants to cause ALS. Here, we experimentally identify distinguishing attributes of ALS mutant SOD proteins that correlate with clinical severity by applying solution biophysical techniques to six ALS mutants at human SOD hotspot glycine 93. A small-angle X-ray scattering (SAXS) assay and other structural methods assessed aggregation propensity by defining the size and shape of fibrillar SOD aggregates after mild biochemical perturbations. Inductively coupled plasma MS quantified metal ion binding stoichiometry, and pulsed dipolar ESR spectroscopy evaluated the Cu²⁺ binding site and defined cross-dimer copper–copper distance distributions. Importantly, we find that copper deficiency in these mutants promotes aggregation in a manner strikingly consistent with their clinical severities. G93 mutants seem to properly incorporate metal ions under physiological conditions when assisted by the copper chaperone but release copper under destabilizing conditions more readily than the WT enzyme. Altered intradimer flexibility in ALS mutants may cause differential metal retention and promote distinct aggregation trends observed for mutant proteins in vitro and in ALS patients. Combined biophysical and structural results test and link copper retention to the framework destabilization hypothesis as a unifying general mechanism for both SOD aggregation and ALS disease progression, with implications for disease severity and therapeutic intervention strategies.ALS is a lethal degenerative disease of the human motor system (1). Opportunities for improved understanding and clinical intervention arose from the discovery that up to 23.5% of familial ALS cases and 7% of spontaneous cases are caused by mutations in the superoxide dismutase 1 (SOD1) gene encoding human Cu, Zn SOD (2–4). SOD is a highly conserved (5), dimeric, antioxidant metalloenzyme that detoxifies superoxide radicals (6, 7), but overexpression of SOD1 ALS mutants is sufficient to cause disease in mice (8). Misfolded and/or aggregated SOD species are deposited within mouse neuronal and glial inclusions (9, 10), even before symptoms appear (11, 12). Although human familial ALS has a symptomatic phenotype indistinguishable from sporadic cases (13), individual SOD1 mutations can result in highly variable disease progression and penetrance (14, 15).Many nongeneral mechanisms, including loss of activity or gain of function, were postulated to explain the roles of SOD mutants in ALS (3, 16–19). Recently, however, an initial hypothesis proposing that SOD manifests disease symptoms by framework destabilization (protein instability caused by structural defects) and consequent protein misassembly and aggregation has gained renewed support (2, 10, 14, 20–23). Ironically, WT SOD is an unusually stable protein (7, 24–26), and precisely how SOD mutations cause disease remains unclear. For instance, human SOD free cysteine residues C6 and C111 have been implicated in protein aggregation by promoting cross-linking (27, 28) and/or stability changes associated with oxidative modifications (29–33). Mutation of the chemically reactive thiols significantly decreases the irreversible denaturation rate for human and bovine SOD (24, 34). However, ALS mutants in a C6A/C111S SOD (AS-SOD) background (35, 36) maintain the native C57–C146 disulfide bond but can still undergo aggregation, and mutations of the free cysteines can cause ALS (37, 38). These results imply that free cysteines are not strictly required but rather, may alter aggregation kinetics (20). SOD also contains two metal ion cofactors in each subunit: a catalytic copper ion (6) and a structurally stabilizing zinc ion (34, 39, 40) (Fig. 1A). In higher eukaryotes, a copper chaperone for SOD (CCS) plays an important role in catalyzing both the copper incorporation and native disulfide bond formation (41). Structural analyses of apo WT SOD point to greater flexibility or increased solvent accessibility of C6 otherwise buried in the stable dimer interface (42, 43), and molecular dynamics simulations also suggest a critical role for metal ions in protein structure, because SOD’s β-sheet propensity decreases in the absence of metals (44). As a result, apo SOD readily forms protein aggregates (45, 46), but the molecular structures of SOD aggregates are likely polymorphic and represent a controversial topic (23, 47–51). The intertwined effects of the aggregation-enhancing free cysteines, dimer-stabilizing metal ions, and CCS maturation of SOD complicate the study of the ALS-causing SOD mutations themselves, and therefore, a clear cause-and-effect relationship remains obscure and requires deconvolution.Open in a separate windowFig. 1.Comparison of crystallographic and solution structures of WT and G93A SOD. (A) Overall architecture of the WT SOD dimer is displayed in 90° rotated views. G93 (small red spheres) resides on a surface-exposed interstrand loop between the fifth and sixth sequential β-strands of SOD and is expected to be innocuous in facilitating protein stability; however, this site harbors the most substitutions observed to result in ALS. G93 is also distant from both (Upper) the dimer interface and (Lower Left) the SOD active site (gold and silver spheres), which are generally implicated as the major determinants for SOD stability. Small blue spheres denote free cysteines. (Lower Right) The close-up view of the mutation site (boxed region in Lower Left tilted forward) shows high similarity between WT (purple) and G93A (red) SOD crystal structures [Protein Data Bank ID codes 1PU0 (WT) and 2ZKY (G93A)]. Hydrogen bonds characteristic of a β-bulge motif are indicated, whereby G93 (or A93) represents position 1. The main chain carbonyl group of β-barrel cork residue L38 is adjacent to the G93 site. (B) SAXS-derived electron pair P(r) distributions from WT (purple) and G93A (red) SOD samples in solution are compared with the theoretical curve for 1PU0. P(r) plots are normalized to peak height. Ab initio models of WT SOD derived from P(r) data are depicted in purple, with crystal structure docked into mesh envelope. Contributions to major and minor peaks from subunit and dimer dimensions are indicated.To better understand the structural effects of ALS mutations on SOD architecture, we coupled the wealth of crystallographic knowledge on SOD structure (7, 52, 53) with small-angle X-ray scattering (SAXS) experiments to characterize misassembly aggregates of ALS mutant SODs in solution. Over 20 y ago, we solved the first atomic structure of the human WT SOD protein (Fig. 1A) (20, 34) and proposed the framework destabilization hypothesis to explain how diverse mutations located throughout the 153-residue β-barrel enzyme might produce a similar disease phenotype (2), albeit with distinctions in the progression trajectory. Since that time, a staggering number of ALS mutations has been documented in patients [178 (mostly missense) (54)], with a similar phenotype in dogs (55, 56). Solution-based techniques are increasingly being applied to connect structure to biological outcome, for instance, through examination of intermolecular interactions within stress-activated pathways, for instance (57, 58). SAXS, which can probe structures for a wide size range of species, also provides higher resolution insights (59), for instance, over visible light-scattering techniques, readily distinguishing unfolded from folded proteins (60).Here, we monitor the initial events of protein aggregation in a subset of ALS mutants localized to a mutational hotspot site at glycine 93. Specifically, we wished to test a possible structural basis for how G93 mutations (to A, C, D, R, S, or V) modulate age of onset and clinical severity in ALS patients (14, 15). The G93 substitution occurs in a β-bulge region (61) between sequential β-strands of the protein (Fig. 1A) on a protruding loop roughly ∼20 Å from T54, the nearest residue of the opposing subunit, and the metal-containing active site (Fig. S1). A priori, mutation of this outer loop position would not be expected to interfere with active site chemistry or buried molecular interfaces. However, we discovered correlations of aggregation nucleation kinetics of SOD proteins with ALS mutations at this site, the stabilizing effects of metal ion retention, and available data for clinical phenotypes in patients with the same mutation. Furthermore, by measuring and exploiting the dimer geometry to observe intrinsic SOD conformers, we show that G93 mutant proteins natively reveal increased intradimer conformational flexibility in the absence of aggregation, which may reflect an increased tendency for ALS mutants to become metal-deficient and misfolding-prone and further explain the correlation to disease severity. Collective results on G93 mutants, thus, support and extend the framework destabilization hypothesis.     

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.