Kinetics of nucleotide-dependent structural transitions in the kinesin-1 hydrolysis cycle

To dissect the kinetics of structural transitions underlying the stepping cycle of kinesin-1 at physiological ATP, we used interferometric scattering microscopy to track the position of gold nanoparticles attached to individual motor domains in processively stepping dimers. Labeled heads resided stably at positions 16.4 nm apart, corresponding to a microtubule-bound state, and at a previously unseen intermediate position, corresponding to a tethered state. The chemical transitions underlying these structural transitions were identified by varying nucleotide conditions and carrying out parallel stopped-flow kinetics assays. At saturating ATP, kinesin-1 spends half of each stepping cycle with one head bound, specifying a structural state for each of two rate-limiting transitions. Analysis of stepping kinetics in varying nucleotides shows that ATP binding is required to properly enter the one-head–bound state, and hydrolysis is necessary to exit it at a physiological rate. These transitions differ from the standard model in which ATP binding drives full docking of the flexible neck linker domain of the motor. Thus, this work defines a consensus sequence of mechanochemical transitions that can be used to understand functional diversity across the kinesin superfamily.Kinesin-1 is a motor protein that steps processively toward microtubule plus-ends, tracking single protofilaments and hydrolyzing one ATP molecule per step (1–6). Step sizes corresponding to the tubulin dimer spacing of 8.2 nm are observed when the molecule is labeled by its C-terminal tail (7–10) and to a two-dimer spacing of 16.4 nm when a single motor domain is labeled (4, 11, 12), consistent with the motor walking in a hand-over-hand fashion. Kinesin has served as an important model system for advancing single-molecule techniques (7–10) and is clinically relevant for its role in neurodegenerative diseases (13), making dissection of its step a popular ongoing target of study.Despite decades of work, many essential components of the mechanochemical cycle remain disputed, including (i) how much time kinesin-1 spends in a one-head–bound (1HB) state when stepping at physiological ATP concentrations, (ii) whether the motor waits for ATP in a 1HB or two-heads–bound (2HB) state, and (iii) whether ATP hydrolysis occurs before or after tethered head attachment (4, 11, 14–20). These questions are important because they are fundamental to the mechanism by which kinesins harness nucleotide-dependent structural changes to generate mechanical force in a manner optimized for their specific cellular tasks. Addressing these questions requires characterizing a transient 1HB state in the stepping cycle in which the unattached head is located between successive binding sites on the microtubule. This 1HB intermediate is associated with the force-generating powerstroke of the motor and underlies the detachment pathway that limits motor processivity. Optical trapping (7, 19, 21, 22) and single-molecule tracking studies (4, 8–11) have failed to detect this 1HB state during stepping. Single-molecule fluorescence approaches have detected a 1HB intermediate at limiting ATP concentrations (11, 12, 14, 15), but apart from one study that used autocorrelation analysis to detect a 3-ms intermediate (17), the 1HB state has been undetectable at physiological ATP concentrations.Single-molecule microscopy is a powerful tool for studying the kinetics of structural changes in macromolecules (23). Tracking steps and potential substeps for kinesin-1 at saturating ATP has until now been hampered by the high stepping rates of the motor (up to 100 s⁻¹), which necessitates high frame rates, and the small step size (8.2 nm), which necessitates high spatial precision (7). Here, we apply interferometric scattering microscopy (iSCAT), a recently established single-molecule tool with high spatiotemporal resolution (24–27) to directly visualize the structural changes underlying kinesin stepping. By labeling one motor domain in a dimeric motor, we detect a 1HB intermediate state in which the tethered head resides over the bound head for half the duration of the stepping cycle at saturating ATP. We further show that at physiological stepping rates, ATP binding is required to enter this 1HB state and that ATP hydrolysis is required to exit it. This work leads to a significant revision of the sequence and kinetics of mechanochemical transitions that make up the kinesin-1 stepping cycle and provides a framework for understanding functional diversity across the kinesin superfamily.     

Quantum mechanical calculations suggest that lytic polysaccharide monooxygenases use a copper-oxyl,oxygen-rebound mechanism

Lytic polysaccharide monooxygenases (LPMOs) exhibit a mononuclear copper-containing active site and use dioxygen and a reducing agent to oxidatively cleave glycosidic linkages in polysaccharides. LPMOs represent a unique paradigm in carbohydrate turnover and exhibit synergy with hydrolytic enzymes in biomass depolymerization. To date, several features of copper binding to LPMOs have been elucidated, but the identity of the reactive oxygen species and the key steps in the oxidative mechanism have not been elucidated. Here, density functional theory calculations are used with an enzyme active site model to identify the reactive oxygen species and compare two hypothesized reaction pathways in LPMOs for hydrogen abstraction and polysaccharide hydroxylation; namely, a mechanism that employs a η¹-superoxo intermediate, which abstracts a substrate hydrogen and a hydroperoxo species is responsible for substrate hydroxylation, and a mechanism wherein a copper-oxyl radical abstracts a hydrogen and subsequently hydroxylates the substrate via an oxygen-rebound mechanism. The results predict that oxygen binds end-on (η¹) to copper, and that a copper-oxyl–mediated, oxygen-rebound mechanism is energetically preferred. The N-terminal histidine methylation is also examined, which is thought to modify the structure and reactivity of the enzyme. Density functional theory calculations suggest that this posttranslational modification has only a minor effect on the LPMO active site structure or reactivity for the examined steps. Overall, this study suggests the steps in the LPMO mechanism for oxidative cleavage of glycosidic bonds.Carbohydrates are the most diverse set of biomolecules, and thus, many enzyme classes have evolved to assemble, modify, and depolymerize carbohydrates, including glycosyltransferases, glycoside hydrolases, carbohydrate esterases, and polysaccharide lyases (1). Recently, a new enzymatic paradigm was discovered that employs copper-dependent oxidation to cleave glycosidic bonds in polysaccharides (2–13). These newly classified enzymes, termed lytic polysaccharide monooxygenases (LPMOs), broadly resemble other copper monooxygenases and some hydroxylation catalysts (14–21).The discovery that LPMOs use an oxidative mechanism has attracted interest both because it is a unique paradigm for carbohydrate modification that employs a powerful C–H activation mechanism, and also because LPMOs synergize with hydrolytic enzymes in biomass conversion to sugars because they act directly on the crystalline polysaccharide surface without the requirement for depolymerization (4, 22, 23), making them of interest in biofuels production. LPMOs were originally characterized as Family 61 glycoside hydrolases (GH61s, reclassified as auxiliary activity 9, AA9) or Family 33 carbohydrate-binding modules (CBM33s, reclassified as AA10), which are structurally similar enzymes found in fungi and nonfungal organisms (22), respectively. In 2005, Vaaje-Kolstad et al. described the synergism (24) of a chitin-active CBM33 (chitin-binding protein, CBP21) with hydrolases, but the mechanism was not apparent. Harris et al. demonstrated that a GH61 boosts hydrolytic enzyme activity on lignocellulosic biomass (2). Vaaje-Kolstad et al. subsequently showed that CBP21 employs an oxidative mechanism to cleave glycosidic linkages in chitin (4).Following these initial discoveries, multiple features of LPMOs have been elucidated. LPMOs use copper (5–7) and produce either aldonic acids or 4-keto sugars at oxidized chain ends, believed to result from hydroxylation at the C1 or C4 carbon, respectively. Hydroxylation at the C1 carbon is proposed to spontaneously undergo elimination to a lactone followed by hydrolytic ring opening to an aldonic acid, whereas hydroxylation and elimination at C4 yields a 4-keto sugar at the nonreducing end (5–12). The active site is a mononuclear type(II) copper center ligated by a “histidine brace” (5, 12), comprising a bidentate N-terminal histidine ligand via the amino terminus and an imidazole ring nitrogen atom and another histidine residue also via a ring nitrogen atom. Hemsworth et al. reported a bacterial LPMO structure wherein the active site copper ion was photoreduced to Cu(I) (12), and Aachmann et al. demonstrated that Cu(I) binds with higher affinity than Cu(II) in CBP21 (13). A structural study of a fungal LPMO revealed an N-terminal methylation on a nitrogen atom in the imidazole ring of unknown function (5), but some LPMOs are active without this modification (6, 11). LPMOs require reducing agents for activity such as ascorbate (2–8, 10–12), and cellobiose dehydrogenase (CDH), a common fungal secretome component, can potentiate LPMO activity in lieu of a small-molecule reducing agent (7, 8).Overall, many structural and mechanistic insights have been reported since the discoveries that LPMOs are oxidative enzymes (4–10). However, many questions remain regarding LPMO function (22, 25). Here, we examine the LPMO catalytic mechanism with density functional theory (DFT) calculations on an active site model (ASM) of a fungal LPMO. We seek to (i) understand the identity of the reactive oxygen species (ROS), (ii) compare two hypothesized catalytic mechanisms, and (iii) examine the role of N-terminal methylation in catalysis.     

Convergent evolution toward an improved growth rate and a reduced resistance range in Prochlorococcus strains resistant to phage

Prochlorococcus is an abundant marine cyanobacterium that grows rapidly in the environment and contributes significantly to global primary production. This cyanobacterium coexists with many cyanophages in the oceans, likely aided by resistance to numerous co-occurring phages. Spontaneous resistance occurs frequently in Prochlorococcus and is often accompanied by a pleiotropic fitness cost manifested as either a reduced growth rate or enhanced infection by other phages. Here, we assessed the fate of a number of phage-resistant Prochlorococcus strains, focusing on those with a high fitness cost. We found that phage-resistant strains continued evolving toward an improved growth rate and a narrower resistance range, resulting in lineages with phenotypes intermediate between those of ancestral susceptible wild-type and initial resistant substrains. Changes in growth rate and resistance range often occurred in independent events, leading to a decoupling of the selection pressures acting on these phenotypes. These changes were largely the result of additional, compensatory mutations in noncore genes located in genomic islands, although genetic reversions were also observed. Additionally, a mutator strain was identified. The similarity of the evolutionary pathway followed by multiple independent resistant cultures and clones suggests they undergo a predictable evolutionary pathway. This process serves to increase both genetic diversity and infection permutations in Prochlorococcus populations, further augmenting the complexity of the interaction network between Prochlorococcus and its phages in nature. Last, our findings provide an explanation for the apparent paradox of a multitude of resistant Prochlorococcus cells in nature that are growing close to their maximal intrinsic growth rates.Large bacterial populations are present in the oceans, playing important roles in primary production and the biogeochemical cycling of matter. These bacterial communities are highly diverse (1–4) yet form stable and reproducible bacterial assemblages under similar environmental conditions (5–7).These bacteria are present together with high abundances of viruses (phages) that have the potential to infect and kill them (8–11). Although studied only rarely in marine organisms (12–16), this coexistence is likely to be the result of millions of years of coevolution between these antagonistic interacting partners, as has been well documented for other systems (17–20). From the perspective of the bacteria, survival entails the selection of cells that are resistant to infection, preventing viral production and enabling the continuation of the cell lineage. Resistance mechanisms include passively acquired spontaneous mutations in cell surface molecules that prevent phage entry into the cell and other mechanisms that actively terminate phage infection intracellularly, such as restriction–modification systems and acquired resistance by CRISPR-Cas systems (21, 22). Mutations in the phage can also occur that circumvent these host defenses and enable the phage to infect the recently emerged resistant bacterium (23).Acquisition of resistance by bacteria is often associated with a fitness cost. This cost is frequently, but not always, manifested as a reduction in growth rate (24–27). Recently, an additional type of cost of resistance was identified, that of enhanced infection whereby resistance to one phage leads to greater susceptibility to other phages (14, 15, 28).Over the years, a number of models have been developed to explain coexistence in terms of the above coevolutionary processes and their costs (16, 29–32). In the arms race model, repeated cycles of host mutation and virus countermutation occur, leading to increasing breadths of host resistance and viral infectivity. However, experimental evidence generally indicates that such directional arms race dynamics do not continue indefinitely (25, 33, 34). Therefore, models of negative density-dependent fluctuations due to selective trade-offs, such as kill-the-winner, are often invoked (20, 33, 35, 36). In these models, fluctuations are generally considered to occur between rapidly growing competition specialists that are susceptible to infection and more slowly growing resistant strains that are considered defense specialists. Such negative density-dependent fluctuations are also likely to occur between strains that have differences in viral susceptibility ranges, such as those that would result from enhanced infection (30).The above coevolutionary processes are considered to be among the major mechanisms that have led to and maintain diversity within bacterial communities (32, 35, 37–39). These processes also influence genetic microdiversity within populations of closely related bacteria. This is especially the case for cell surface-related genes that are often localized to genomic islands (14, 40, 41), regions of high gene content, and gene sequence variability among members of a population. As such, populations in nature display an enormous degree of microdiversity in phage susceptibility regions, potentially leading to an assortment of subpopulations with different ranges of susceptibility to coexisting phages (4, 14, 30, 40).Prochlorococcus is a unicellular cyanobacterium that is the numerically dominant photosynthetic organism in vast oligotrophic expanses of the open oceans, where it contributes significantly to primary production (42, 43). Prochlorococcus consists of a number of distinct ecotypes (44–46) that form stable and reproducible population structures (7). These populations coexist in the oceans with tailed double-stranded DNA phage populations that infect them (47–49).Previously, we found that resistance to phage infection occurs frequently in two high-light–adapted Prochlorococcus ecotypes through spontaneous mutations in cell surface-related genes (14). These genes are primarily localized to genomic island 4 (ISL4) that displays a high degree of genetic diversity in environmental populations (14, 40). Although about a third of Prochlorococcus-resistant strains had no detectable associated cost, the others came with a cost manifested as either a slower growth rate or enhanced infection by other phages (14). In nature, Prochlorococcus seems to be growing close to its intrinsic maximal growth rate (50–52). This raises the question as to the fate of emergent resistant Prochlorococcus lineages in the environment, especially when resistance is accompanied with a high growth rate fitness cost.To begin addressing this question, we investigated the phenotype of Prochlorococcus strains with time after the acquisition of resistance. We found that resistant strains evolved toward an improved growth rate and a reduced resistance range. Whole-genome sequencing and PCR screening of many of these strains revealed that these phenotypic changes were largely due to additional, compensatory mutations, leading to increased genetic diversity. These findings suggest that the oceans are populated with rapidly growing Prochlorococcus cells with varying degrees of resistance and provide an explanation for how a multitude of presumably resistant Prochlorococcus cells are growing close to their maximal known growth rate in nature.     

Mechanisms of NDV-3 vaccine efficacy in MRSA skin versus invasive infection

Increasing rates of life-threatening infections and decreasing susceptibility to antibiotics urge development of an effective vaccine targeting Staphylococcus aureus. This study evaluated the efficacy and immunologic mechanisms of a vaccine containing a recombinant glycoprotein antigen (NDV-3) in mouse skin and skin structure infection (SSSI) due to methicillin-resistant S. aureus (MRSA). Compared with adjuvant alone, NDV-3 reduced abscess progression, severity, and MRSA density in skin, as well as hematogenous dissemination to kidney. NDV-3 induced increases in CD3+ T-cell and neutrophil infiltration and IL-17A, IL-22, and host defense peptide expression in local settings of SSSI abscesses. Vaccine induction of IL-22 was necessary for protective mitigation of cutaneous infection. By comparison, protection against hematogenous dissemination required the induction of IL-17A and IL-22 by NDV-3. These findings demonstrate that NDV-3 protective efficacy against MRSA in SSSI involves a robust and complementary response integrating innate and adaptive immune mechanisms. These results support further evaluation of the NDV-3 vaccine to address disease due to S. aureus in humans.The bacterium Staphylococcus aureus is the leading cause of skin and skin structure infections (SSSIs), including cellulitis, furunculosis, and folliculitis (1–4), and a common etiologic agent of impetigo (5), erysipelas (6), and superinfection in atopic dermatitis (7). This bacterium is a significant cause of surgical or traumatic wound infections (8, 9), as well as decuibitus and diabetic skin lesions (10). Moreover, SSSI is an important risk factor for systemic infection. The skin is a key portal of entry for hematogenous dissemination, particularly in association with i.v. catheters. S. aureus is now the second most common bloodstream isolate in healthcare settings (11), and SSSI is a frequent source of invasive infections such as pneumonia or endocarditis (12, 13). Despite a recent modest decline in rates of methicillin-resistant S. aureus (MRSA) infection in some cohorts (13), infections due to S. aureus remain a significant problem (14, 15). Even with appropriate therapy, up to one-third of patients diagnosed with S. aureus bacteremia succumb—accounting for more attributable annual deaths than HIV, tuberculosis, and viral hepatitis combined (16).The empiric use of antibiotics in healthcare-associated and community-acquired settings has increased S. aureus exposure to these agents, accelerating selection of resistant strains. As a result, resistance to even the most recently developed agents is emerging at an alarming pace (17, 18). The impact of this trend is of special concern in light of high rates of mortality associated with invasive MRSA infection (e.g., 15–40% in bacteremia or endocarditis), even with the most recently developed antistaphylococcal therapeutics (19, 20). Moreover, patients who experience SSSI due to MRSA exhibit high 1-y recurrence rates, often prompting surgical debridement (21) and protracted antibiotic treatment.Infections due to MRSA are a special concern in immune-vulnerable populations, including hemodialysis (22), neutropenic (23, 24), transplantation (25), and otherwise immunosuppressed patients (26, 27), and in patients with inherited immune dysfunctions (28–31) or cystic fibrosis (32). Patients having deficient interleukin 17 (IL-17) or IL-22 responses (e.g., signal transduction mediators STAT3, DOCK8, or CARD9 deficiencies) exhibit chronic or “cold” abscesses, despite high densities of pathogens such as S. aureus (33, 34). For example, patients with Chronic Granulomatous Disease (CGD; deficient Th1 and oxidative burst response) have increased risk of disseminated S. aureus infection. In contrast, patients with Job’s Syndrome (deficient Th17 response) typically have increased risk to SSSI and lung infections, but less so for systemic S. aureus bacteremia (35, 36). This pattern contrasts that observed in neutropenic or CGD patients (37). These themes suggest efficacious host defenses against MRSA skin and invasive infections involve complementary but distinct molecular and cellular immune responses.From these perspectives, vaccines or immunotherapeutics that prevent or lessen severity of MRSA infections, or that enhance antibiotic efficacy, would be significant advances in patient care and public health. However, to date, there are no licensed prophylactic or therapeutic vaccine immunotherapies for S. aureus or MRSA infection. Unfortunately, efforts to develop vaccines targeting S. aureus capsular polysaccharide type 5 or 8 conjugates, or the iron-regulated surface determinant B protein, have not been successful thus far (38, 39). Likewise, passive immunization using monoclonal antibodies targeting the S. aureus adhesin clumping factor A (ClfA, tefibazumab) (40) or lipoteichoic acid (pagibaximab) (41) have not shown efficacy against invasive infections in human clinical studies to date. Moreover, the striking recurrence rates of SSSI due to MRSA imply that natural exposure does not induce optimal preventive immunity or durable anamnestic response to infection or reinfection. Thus, significant challenges exist in the development of an efficacious vaccine targeting diseases caused by S. aureus (42) that are perhaps not optimally addressed by conventional approaches.The NDV-3 vaccine reflects a new strategy to induce durable immunity targeting S. aureus. Its immunogen is engineered from the agglutinin-like sequence 3 (Als3) adhesin/invasin of Candida albicans, which we discovered to be a structural homolog of S. aureus adhesins (43). NDV-3 is believed to cross-protect against S. aureus and C. albicans due to sequence (T-cell) and conformational (B-cell) epitopes paralleled in both organisms (44). Our prior data have shown that NDV-3 is efficacious in murine models of hematogenous and mucosal candidiasis (45), as well as S. aureus bacteremia (46–48). Recently completed phase I clinical trials demonstrate the safety, tolerability, and immunogenicity of NDV-3 in humans (49).     

Novel serologic biomarkers provide accurate estimates of recent Plasmodium falciparum exposure for individuals and communities

Tools to reliably measure Plasmodium falciparum (Pf) exposure in individuals and communities are needed to guide and evaluate malaria control interventions. Serologic assays can potentially produce precise exposure estimates at low cost; however, current approaches based on responses to a few characterized antigens are not designed to estimate exposure in individuals. Pf-specific antibody responses differ by antigen, suggesting that selection of antigens with defined kinetic profiles will improve estimates of Pf exposure. To identify novel serologic biomarkers of malaria exposure, we evaluated responses to 856 Pf antigens by protein microarray in 186 Ugandan children, for whom detailed Pf exposure data were available. Using data-adaptive statistical methods, we identified combinations of antibody responses that maximized information on an individual’s recent exposure. Responses to three novel Pf antigens accurately classified whether an individual had been infected within the last 30, 90, or 365 d (cross-validated area under the curve = 0.86–0.93), whereas responses to six antigens accurately estimated an individual’s malaria incidence in the prior year. Cross-validated incidence predictions for individuals in different communities provided accurate stratification of exposure between populations and suggest that precise estimates of community exposure can be obtained from sampling a small subset of that community. In addition, serologic incidence predictions from cross-sectional samples characterized heterogeneity within a community similarly to 1 y of continuous passive surveillance. Development of simple ELISA-based assays derived from the successful selection strategy outlined here offers the potential to generate rich epidemiologic surveillance data that will be widely accessible to malaria control programs.Many countries have extensive programs to reduce the burden of Plasmodium falciparum (Pf), the parasite responsible for most malaria morbidity and mortality (1). Effectively using limited resources for malaria control or elimination and evaluating interventions require accurate measurements of the risk of being infected with Pf (2–15). To reflect the rate at which individuals are infected with Pf in a useful way, metrics used to estimate exposure in a community need to account for dynamic changes over space and time, especially in response to control interventions (16–18).A variety of metrics can be used to estimate Pf exposure, but tools that are more precise and low cost are needed for population surveillance. Existing metrics have varying intrinsic levels of precision and accuracy and are subject to a variety of extrinsic factors, such as cost, time, and availability of trained personnel (19). For example, entomological measurements provide information on mosquito to human transmission for a community but are expensive, require specially trained staff, and lack standardized procedures, all of which reduce precision and/or make interpretation difficult (19–22). Parasite prevalence can be measured by detecting parasites in the blood of individuals from a cross-sectional sample of a community and is, therefore, relatively simple and inexpensive to perform, but results may be imprecise, especially in areas of low transmission (19, 23), and biased by a number of factors, including immunity and access to antimalarial treatment (5, 6, 19, 23–25). The burden of symptomatic disease in a community can be estimated from routine health systems data; however, such data are frequently unreliable (5, 26–28) and generally underestimate the prevalence of Pf infection in areas of intense transmission. Precise and quantitative information about exposure at an individual level can be reliably obtained from cohort studies by measuring the incidence of asymptomatic and/or symptomatic Pf infection (i.e., by measuring the molecular force of infection) (29–35). Unfortunately, the expense of cohort studies limits their use to research settings. The end result is that most malaria-endemic regions lack reliable, timely data on Pf exposure, limiting the capabilities of malaria control programs to guide and evaluate interventions.Serologic assays offer the potential to provide incidence estimates for symptomatic and asymptomatic Pf infection, which are currently obtained from cohort studies, at the cost of cross-sectional studies (36–38). Although Pf infections are transient, a record of infection remains detectable in an individual’s antibody profile. Thus, appropriately chosen antibody measurements integrated with age can provide information about an individual’s exposure history. Antibodies can be measured by simple ELISAs and obtained from dried blood spots, which are easy to collect, transport, and store (39–41). Serologic responses to Pf antigens have been explored as potential epidemiological tools (42–45), and estimated rates of seroconversion to well-characterized Pf antigens accurately reflect stable rates of exposure in a community, whereas distinct changes in these rates are obtained from successful interventions (22, 39, 41, 46–53). However, current serologic assays are not designed to detect short-term or gradual changes in Pf exposure or measure exposure to infection at an individual level. The ability to calibrate antibody responses to estimates of exposure in individuals could allow for more flexible sampling of a population (e.g., not requiring age stratification), improve accuracy of exposure estimates from small sample sizes, and better characterize heterogeneity in exposure within a community.Different Pf antigens elicit antibody responses with different magnitudes and kinetics, providing a large and diverse set of potential biomarkers for exposure (38, 54–58). We hypothesized that new and more highly informative serologic biomarkers better able to characterize an individual’s recent exposure history could be identified by analyzing antibody responses to a large number of candidate Pf antigens in participants with well-characterized exposure histories. To test this hypothesis, we probed plasma from participants in two cohort studies in Uganda against a protein microarray containing 856 Pf antigens. The primary aim of this analysis was to identify responses to select antigens that were most informative of recent exposure using robust, data-adaptive statistical methods. Each participant’s responses to these selected antigens were used as predictors for two primary outcomes of their recent exposure to Pf: (i) days since last Pf infection and (ii) the incidence of symptomatic malaria in the last year. These individual-level estimates were then aggregated across a population to assess community-level malaria exposure. The selection strategy presented here identified accurate biomarkers of exposure for children living in areas of moderate to high Pf exposure and illustrates the utility of this flexible and broadly applicable approach.     

Structural interactions of a voltage sensor toxin with lipid membranes

Protein toxins from tarantula venom alter the activity of diverse ion channel proteins, including voltage, stretch, and ligand-activated cation channels. Although tarantula toxins have been shown to partition into membranes, and the membrane is thought to play an important role in their activity, the structural interactions between these toxins and lipid membranes are poorly understood. Here, we use solid-state NMR and neutron diffraction to investigate the interactions between a voltage sensor toxin (VSTx1) and lipid membranes, with the goal of localizing the toxin in the membrane and determining its influence on membrane structure. Our results demonstrate that VSTx1 localizes to the headgroup region of lipid membranes and produces a thinning of the bilayer. The toxin orients such that many basic residues are in the aqueous phase, all three Trp residues adopt interfacial positions, and several hydrophobic residues are within the membrane interior. One remarkable feature of this preferred orientation is that the surface of the toxin that mediates binding to voltage sensors is ideally positioned within the lipid bilayer to favor complex formation between the toxin and the voltage sensor.Protein toxins from venomous organisms have been invaluable tools for studying the ion channel proteins they target. For example, in the case of voltage-activated potassium (Kv) channels, pore-blocking scorpion toxins were used to identify the pore-forming region of the channel (1, 2), and gating modifier tarantula toxins that bind to S1–S4 voltage-sensing domains have helped to identify structural motifs that move at the protein–lipid interface (3–5). In many instances, these toxin–channel interactions are highly specific, allowing them to be used in target validation and drug development (6–8).Tarantula toxins are a particularly interesting class of protein toxins that have been found to target all three families of voltage-activated cation channels (3, 9–12), stretch-activated cation channels (13–15), as well as ligand-gated ion channels as diverse as acid-sensing ion channels (ASIC) (16–21) and transient receptor potential (TRP) channels (22, 23). The tarantula toxins targeting these ion channels belong to the inhibitor cystine knot (ICK) family of venom toxins that are stabilized by three disulfide bonds at the core of the molecule (16, 17, 24–31). Although conventional tarantula toxins vary in length from 30 to 40 aa and contain one ICK motif, the recently discovered double-knot toxin (DkTx) that specifically targets TRPV1 channels contains two separable lobes, each containing its own ICK motif (22, 23).One unifying feature of all tarantula toxins studied thus far is that they act on ion channels by modifying the gating properties of the channel. The best studied of these are the tarantula toxins targeting voltage-activated cation channels, where the toxins bind to the S3b–S4 voltage sensor paddle motif (5, 32–36), a helix-turn-helix motif within S1–S4 voltage-sensing domains that moves in response to changes in membrane voltage (37–41). Toxins binding to S3b–S4 motifs can influence voltage sensor activation, opening and closing of the pore, or the process of inactivation (4, 5, 36, 42–46). The tarantula toxin PcTx1 can promote opening of ASIC channels at neutral pH (16, 18), and DkTx opens TRPV1 in the absence of other stimuli (22, 23), suggesting that these toxin stabilize open states of their target channels.For many of these tarantula toxins, the lipid membrane plays a key role in the mechanism of inhibition. Strong membrane partitioning has been demonstrated for a range of toxins targeting S1–S4 domains in voltage-activated channels (27, 44, 47–50), and for GsMTx4 (14, 50), a tarantula toxin that inhibits opening of stretch-activated cation channels in astrocytes, as well as the cloned stretch-activated Piezo1 channel (13, 15). In experiments on stretch-activated channels, both the d- and l-enantiomers of GsMTx4 are active (14, 50), implying that the toxin may not bind directly to the channel. In addition, both forms of the toxin alter the conductance and lifetimes of gramicidin channels (14), suggesting that the toxin inhibits stretch-activated channels by perturbing the interface between the membrane and the channel. In the case of Kv channels, the S1–S4 domains are embedded in the lipid bilayer and interact intimately with lipids (48, 51, 52) and modification in the lipid composition can dramatically alter gating of the channel (48, 53–56). In one study on the gating of the Kv2.1/Kv1.2 paddle chimera (53), the tarantula toxin VSTx1 was proposed to inhibit Kv channels by modifying the forces acting between the channel and the membrane. Although these studies implicate a key role for the membrane in the activity of Kv and stretch-activated channels, and for the action of tarantula toxins, the influence of the toxin on membrane structure and dynamics have not been directly examined. The goal of the present study was to localize a tarantula toxin in membranes using structural approaches and to investigate the influence of the toxin on the structure of the lipid bilayer.     

Dissecting neural pathways for forgetting in Drosophila olfactory aversive memory

Recent studies have identified molecular pathways driving forgetting and supported the notion that forgetting is a biologically active process. The circuit mechanisms of forgetting, however, remain largely unknown. Here we report two sets of Drosophila neurons that account for the rapid forgetting of early olfactory aversive memory. We show that inactivating these neurons inhibits memory decay without altering learning, whereas activating them promotes forgetting. These neurons, including a cluster of dopaminergic neurons (PAM-β′1) and a pair of glutamatergic neurons (MBON-γ4>γ1γ2), terminate in distinct subdomains in the mushroom body and represent parallel neural pathways for regulating forgetting. Interestingly, although activity of these neurons is required for memory decay over time, they are not required for acute forgetting during reversal learning. Our results thus not only establish the presence of multiple neural pathways for forgetting in Drosophila but also suggest the existence of diverse circuit mechanisms of forgetting in different contexts.Although forgetting commonly has a negative connotation, it is a functional process that shapes memory and cognition (1–4). Recent studies, including work in relatively simple invertebrate models, have started to reveal basic biological mechanisms underlying forgetting (5–15). In Drosophila, single-session Pavlovian conditioning by pairing an odor (conditioned stimulus, CS) with electric shock (unconditioned stimulus, US) induces aversive memories that are short-lasting (16). The memory performance of fruit flies is observed to drop to a negligible level within 24 h, decaying rapidly early after training and slowing down thereafter (17). Memory decay or forgetting requires the activation of the small G protein Rac, a signaling protein involved in actin remodeling, in the mushroom body (MB) intrinsic neurons (6). These so-called Kenyon cells (KCs) are the neurons that integrate CS–US information (18, 19) and support aversive memory formation and retrieval (20–22). In addition to Rac, forgetting also requires the DAMB dopamine receptor (7), which has highly enriched expression in the MB (23). Evidence suggests that the dopamine-mediated forgetting signal is conveyed to the MB by dopamine neurons (DANs) in the protocerebral posterior lateral 1 (PPL1) cluster (7, 24). Therefore, forgetting of olfactory aversive memory in Drosophila depends on a particular set of intracellular molecular pathways within KCs, involving Rac, DAMB, and possibly others (25), and also receives modulation from extrinsic neurons. Although important cellular evidence supporting the hypothesis that memory traces are erased under these circumstances is still lacking, these findings lend support to the notion that forgetting is an active, biologically regulated process (17, 26).Although existing studies point to the MB circuit as essential for forgetting, several questions remain to be answered. First, whereas the molecular pathways for learning and forgetting of olfactory aversive memory are distinct and separable (6, 7), the neural circuits seem to overlap. Rac-mediated forgetting has been localized to a large population of KCs (6), including the γ-subset, which is also critical for initial memory formation (21, 27). The site of action of DAMB for forgetting has yet to be established; however, the subgroups of PPL1-DANs implicated in forgetting are the same as those that signal aversive reinforcement and are required for learning (28–30). It leaves open the question of whether the brain circuitry underlying forgetting and learning is dissociable, or whether forgetting and learning share the same circuit but are driven by distinct activity patterns and molecular machinery (26). Second, shock reinforcement elicits multiple memory traces through at least three dopamine pathways to different subdomains in the MB lobes (28, 29). Functional imaging studies have also revealed Ca²⁺-based memory traces in different KC populations (31). It is poorly understood how forgetting of these memory traces differs, and it remains unknown whether there are multiple regulatory neural pathways. Notably, when PPL1-DANs are inactivated, forgetting still occurs, albeit at a lower rate (7). This incomplete block suggests the existence of an additional pathway(s) that conveys forgetting signals to the MB. Third, other than memory decay over time, forgetting is also observed through interference (32, 33), when new learning or reversal learning is introduced after training (6, 34, 35). Time-based and interference-based forgetting shares a similar dependence on Rac and DAMB (6, 7). However, it is not known whether distinct circuits underlie forgetting in these different contexts.In the current study, we focus on the diverse set of MB extrinsic neurons (MBENs) that interconnect the MB lobes with other brain regions, which include 34 MB output neurons (MBONs) of 21 types and ∼130 dopaminergic neurons of 20 types in the PPL1 and protocerebral anterior medial (PAM) clusters (36, 37). These neurons have been intensively studied in olfactory memory formation, consolidation, and retrieval in recent years (e.g., 24, 28–30, 38–48); however, their roles in forgetting have not been characterized except for the aforementioned PPL1-DANs. In a functional screen, we unexpectedly found that several Gal4 driver lines of MBENs showed significantly better 3-h memory retention when the Gal4-expressing cells were inactivated. The screen has thus led us to identify two types of MBENs that are not involved in initial learning but play important and additive roles in mediating memory decay. Furthermore, neither of these MBEN types is required for reversal learning, supporting the notion that there is a diversity of neural circuits that drive different forms of forgetting.     

UVB radiation generates sunburn pain and affects skin by activating epidermal TRPV4 ion channels and triggering endothelin-1 signaling

At our body surface, the epidermis absorbs UV radiation. UV overexposure leads to sunburn with tissue injury and pain. To understand how, we focus on TRPV4, a nonselective cation channel highly expressed in epithelial skin cells and known to function in sensory transduction, a property shared with other transient receptor potential channels. We show that following UVB exposure mice with induced Trpv4 deletions, specifically in keratinocytes, are less sensitive to noxious thermal and mechanical stimuli than control animals. Exploring the mechanism, we find that epidermal TRPV4 orchestrates UVB-evoked skin tissue damage and increased expression of the proalgesic/algogenic mediator endothelin-1. In culture, UVB causes a direct, TRPV4-dependent Ca²⁺ response in keratinocytes. In mice, topical treatment with a TRPV4-selective inhibitor decreases UVB-evoked pain behavior, epidermal tissue damage, and endothelin-1 expression. In humans, sunburn enhances epidermal expression of TRPV4 and endothelin-1, underscoring the potential of keratinocyte-derived TRPV4 as a therapeutic target for UVB-induced sunburn, in particular pain.The surface epithelium (epidermis) of skin provides barrier protection against dehydration and the potentially harmful external environment (1). Accordingly, skin is the site of first interaction between ambient environment and immunologically competent organismal structures, and also the site for sentient responses (2). Sensory neurons in the dorsal root ganglia (DRG) and trigeminal ganglia (TG) are endowed with sensory transduction capacity for heat, cold, mechanical cues, itch, and pain, and their axons directly interface with skin epithelium (2–4).Against a background of suggestive findings (2, 5–7), we wondered whether the epidermis as a “forefront” of sensory signaling may function in sensitizing pain transduction in response to naturally occurring irritating cues. To elucidate mechanisms, we used a mouse sunburn model and induced a state of lowered sensory thresholds associated with tissue injury caused by UV radiation (8–10). UV-sunburn-evoked lowering of sensory thresholds shares major hallmarks of pathological pain, a valuable feature of this model. Skin tissue injury caused by UVB has been elucidated to be mediated by cytokines and chemokines, known from immunological responses, such as IL-1β and IL-6, which are also known to cause and facilitate pain (11–19). Another more recent study identified a proinflammatory chemokine, CXCL5, as proalgesic in response to UVB overexposure of rat and human skin (20). An exciting new arena pertaining to molecular mechanisms of the skin’s response to noxious UV was recently opened by an elegant study that reported the role of UVB-mediated damage to noncoding RNA molecules in the skin (21). Unraveling a molecular mechanism, the Toll-like receptor 3 gene was found critical in signaling the proinflammatory actions of the UVB-damaged noncoding RNA molecules. However, this study focused on molecular mechanisms of acute inflammation in the skin.We intended to identify pain mechanisms that mediate the pain associated with UVB-mediated tissue injury. Pain in response to external environmental cues has been understood better because of scientific progress in the field of transient receptor potential (TRP) ion channels that have been found responsive to such cues, and which were found expressed in DRG and TG peripheral sensory neurons, which are the cells believed to be the primary transducers. Indeed, TRPV1, one of the founding members of the TRPV channel subfamily, has been identified as relevant for pain, including pathological pain, response to thermal cues, and most recently for itch (22–31). However, TRPA1 (transient receptor potential ion channel, ankyrin subfamily, family member #1) and TRPM8 seem to be involved in transduction of pain-inducing stimuli as well (32–36).Also a family member of the TRPV subfamily, TRPV4 is a multimodally activated, nonselective cation channel that is involved in physiological pain evoked by osmotic and mechanical, but not thermal, cues (37–40). For pathological pain, it is relevant for inflammation- and nerve-damage-induced pain sensitization (41–43). Of note, Trpv4^−/− mice exhibit impaired skin-barrier function (44, 45). That said, TRPV4 is expressed in a number of different cell types, including robust expression in epidermal keratinocytes and also is detectable in skin-innervating sensory neurons. This “dual-location expression” of TRPV4 leaves the cellular mechanisms involved in the channel’s function and the functional contribution of environment-exposed keratinocytes vs. skin-innervating sensory neurons unclear.Against this background of dual-location TRPV4 expression and the role of TRPV4 in inflammatory and neuropathic pain, we now address whether epidermally derived TRPV4 is pathophysiologically relevant in sunburn pain and tissue damage. Using Trpv4 gene-targeted mice, selectively inducing targeting in postnatal keratinocytes, and topically applying selective TRPV4 inhibitors, we demonstrate that epidermal TRPV4 plays a prominent, hitherto unrecognized role in UVB-evoked skin tissue damage and pain of sunburn.     

Cytonuclear interactions affect adaptive traits of the annual plant Arabidopsis thaliana in the field

Although the contribution of cytonuclear interactions to plant fitness variation is relatively well documented at the interspecific level, the prevalence of cytonuclear interactions at the intraspecific level remains poorly investigated. In this study, we set up a field experiment to explore the range of effects that cytonuclear interactions have on fitness-related traits in Arabidopsis thaliana. To do so, we created a unique series of 56 cytolines resulting from cytoplasmic substitutions among eight natural accessions reflecting within-species genetic diversity. An assessment of these cytolines and their parental lines scored for 28 adaptive whole-organism phenotypes showed that a large proportion of phenotypic traits (23 of 28) were affected by cytonuclear interactions. The effects of these interactions varied from slight but frequent across cytolines to strong in some specific parental pairs. Two parental pairs accounted for half of the significant pairwise interactions. In one parental pair, Ct-1/Sha, we observed symmetrical phenotypic responses between the two nuclear backgrounds when combined with specific cytoplasms, suggesting nuclear differentiation at loci involved in cytonuclear epistasis. In contrast, asymmetrical phenotypic responses were observed in another parental pair, Cvi-0/Sha. In the Cvi-0 nuclear background, fecundity and phenology-related traits were strongly affected by the Sha cytoplasm, leading to a modified reproductive strategy without penalizing total seed production. These results indicate that natural variation in cytoplasmic and nuclear genomes interact to shape integrative traits that contribute to adaptation, thereby suggesting that cytonuclear interactions can play a major role in the evolutionary dynamics of A. thaliana.The genomes of eukaryotes originate from ancient endosymbiotic associations that eventually led to energy-harnessing organelles: mitochondria, common to all eukaryotes, and chloroplasts in the “green” lineage. The evolution of endosymbionts into cellular organelles was accompanied by massive gene loss, with a large proportion being transferred to the nucleus (1, 2). Nevertheless, mitochondria and chloroplasts retained a few (30–80) protein-encoding genes that play crucial roles in energy metabolism (respiration and photosynthesis). Mitochondrion and chloroplast metabolisms rely on the proper interaction of nuclear-encoded proteins and their counterparts encoded in the organelle genome. Consequently, the genes in nuclear and organellar compartments are expected to be coadapted (3).Cytonuclear coadaptation has been demonstrated by altered phenotypes observed on interspecific exchanges of cytoplasm between related species in mammals (4), yeast (5), arthropods (6), and plants, whose interspecific crosses are frequently successful (7). These alterations affect organelle function and even the organism phenotype, indicating epistasis between nuclear and cytoplasmic genes. Although cytonuclear coadaptation is generally studied at the interspecific level, the existence of intraspecific genetic diversity in organelle genomes suggests a potential for genomic coadaptation within species. A few studies have reported phenotypic effects of intraspecific cytonuclear epistasis in nonplant species (8–11). In plants, many studies have focused on cytoplasmic male sterility (CMS), an impairment of pollen production governed by nucleo-mitochondrial interactions in some hermaphroditic species (12), in particular in crops and their relatives (13). The phenotypic effects of intraspecific cytonuclear epistasis other than CMS have been reported in only a limited number of plant systems (14–17), with evidence that cytoplasmic variation contributes to local adaptation (18, 19).In recent years, several studies using reciprocal segregating populations of the model plant Arabidopsis thaliana have investigated the effect of cytonuclear epistasis on a number of laboratory-measured phenotypes such as the metabolome, defense chemistry and growth (17, 20, 21), water-use efficiency (22, 23), and seed germination (24, 25). Although some studies have reported significant effects of cytonuclear epistasis (17, 20, 21, 23, 25), others have found additive cytoplasmic effects but with weak or no cytonuclear epistasis (22). Each of these studies (with the exception of ref. 25) was, however, based on a single reciprocal cross between two natural accessions, thereby preventing the estimation of the prevalence of cytonuclear epistasis in this species. In addition, although these reports involve adaptive traits (26–30), the investigation of the effect of cytonuclear epistasis on adaptive phenotypes in field conditions is, at best, scarce in A. thaliana.Here, following the modern standards of ecological genomics (31), we explored the prevalence of cytonuclear interactions on adaptive whole-organism traits in the model plant A. thaliana in a field experiment. To do so, based on eight natural accessions of a core collection that covers a significant part of the species’ cytoplasmic and nuclear genetic diversity in A. thaliana (25, 32), we created eight series of seven cytolines. Cytolines are genotypes that combine the nuclear genome from one parent with the organelle genomes of another (33). We examined the cytolines and their parental accessions for effects of cytonuclear interactions on 28 field-measured traits related to germination, phenology, resource acquisition, plant architecture and seed dispersal, fecundity, and survival.     

Microscopic identification of the order parameter governing liquid–liquid transition in a molecular liquid

A liquid–liquid transition (LLT) in a single-component substance is an unconventional phase transition from one liquid to another. LLT has recently attracted considerable attention because of its fundamental importance in our understanding of the liquid state. To access the order parameter governing LLT from a microscopic viewpoint, here we follow the structural evolution during the LLT of an organic molecular liquid, triphenyl phosphite (TPP), by time-resolved small- and wide-angle X-ray scattering measurements. We find that locally favored clusters, whose characteristic size is a few nanometers, are spontaneously formed and their number density monotonically increases during LLT. This strongly suggests that the order parameter of LLT is the number density of locally favored structures and of nonconserved nature. We also show that the locally favored structures are distinct from the crystal structure and these two types of orderings compete with each other. Thus, our study not only experimentally identifies the structural order parameter governing LLT, but also may settle a long-standing debate on the nature of the transition in TPP, i.e., whether the transition is LLT or merely microcrystal formation.Liquid-liquid transition (LLT) is an intriguing phenomenon in which a liquid transforms into another one via a first-order transition. This means that there can be more than two liquid states for a single-component substance. Despite its counterintuitive nature, there have recently been many pieces of experimental and numerical evidence for the existence of LLT, for various liquids such as water (1–5), aqueous solutions (6–8), triphenyl phosphite (9–12), l-butanol (13), phosphorus (14), silicon (15, 16), germanium (17), and Y₂O₃–Al₂O₃ (18, 19). This suggests that the LLT may be rather universally observed for various types of liquids. However, none of the LLTs reported so far is free from criticisms (20, 21), mainly because these LLTs take place under experimentally difficult conditions [e.g., at high temperature and pressure (14, 15, 17–19)] or in a supercooled state below the melting point (1–3, 5–7, 9, 10), where the transition is inevitably contaminated by microcrystal formation. The latter is not limited to experiments but arises in numerical simulations, often causing many controversies [LLT (22–25) vs. crystallization (26–28)]. For ST2 water, however, this issue has recently been settled by an extensive simulation study by Palmer et al. (4).One of the hottest and long-standing debates is on the nature of the transition found in a molecular liquid, triphenyl phosphite (TPP), by Kivelson and his coworkers (29). The transition is very easy to access experimentally, because it takes place at ambient pressure and at a temperature range between 230 and 210 K and the transformation speed is slow enough to follow the kinetics. Since the finding of this transition (29, 30), many researchers thus have been interested in this intriguing phenomenon and there have been hot discussions on the nature of the transition (20, 21). Some people interpreted this as a liquid-associated phenomenon (9, 10, 31, 32), but others interpret it differently. All of the controversies come from the fact that this transition accompanies microcrystal formation and thus the final state, which is called “glacial phase,” often contains microcrystallites. This led many researchers to explain the transition by non-LLT scenarios, which include a defect-ordered phase scenario predicted by a frustration limited domain theory (29, 30, 33, 34), a microcrystallization scenario (35–38), and a liquid-crystal or plastic-crystal phase scenario (39). Each scenario captures a certain feature of the glacial phase, but fails in explaining all of the experimental results in a consistent manner. Similar situations are often seen in other candidates of LLTs, such as l-butanol [LLT (13) vs. microcrystallization (40–43)], confined water [LLT (5) vs. other phenomena (44–46)], and aqueous solutions [LLT (6, 7) vs. microcrystallization (8, 28, 47, 48)]. For TPP, however, some pieces of experimental evidence supportive of the LLT scenario rather than the microcrystallization scenario have recently been reported (11, 12).We propose a two-order-parameter (TOP) model of a liquid to explain LLT (20, 49). The main point of this model is that it is necessary to consider the spatiotemporal hierarchical nature of a liquid to understand LLT. More specifically, we argue that in addition to density order parameter ρ describing a gas–liquid transition, we need an additional scalar order parameter S, which is the number density of locally favored structures (LFS). In this model, LLT is a consequence of the cooperative ordering of the scalar nonconserved order parameter S, i.e., the cooperative formation of LFS. In other words, LLT is regarded as a gas–liquid-like transition of LFS: one liquid is a gas state of LFS (low-S state), and the other is its liquid state (high-S state). Recently, it was proposed by Anisimov and coworkers (50, 51) that the thermodynamic ordering field conjugate to the order parameter is the conversion equilibrium constant, which further characterizes the nature of LLT. We explained our experimental observation of LLT in TPP in terms of this model (9, 10). We also studied the phase transition dynamics and the physical and chemical properties of the second liquid state (liquid II), which were also explained by the model (20, 21).However, we have not had any direct experimental evidence for the formation of such LFS up to now; thus, an open question is, what is the relevant order parameter governing LLT, although the link of the order parameter to the enthalpy (9, 10), the refractive index (or, density) (9, 10, 29, 30), and the polarity associated with local molecular ordering (12) has been suggested for LLT in TPP. There have been structural studies on LLT by X-ray and neutron scattering measurements, focusing on local liquid structures at an inter- and intramolecular scale (36, 38, 52–54) and mesoscopic structures (34, 55). However, there has been no experimental evidence for the presence of locally favored structures, which characterize the liquid state uniquely, or the order parameter has still not been identified from a microscopic viewpoint.Here we study the structural change of TPP during LLT by time-resolved small- and wide-angle X-ray scattering measurements, which cover a length scale from a single molecule size ( ～  1 nm) to more than tens of nanometers. We show, to our knowledge, the first direct evidence for the presence of LFS and the temporal increase upon the liquid I-to-liquid II transformation. Furthermore, we also find an indication of the formation of microcrystallites during LLT. However, we reveal that LFS and microcrystallites have different sizes and growth kinetics, indicating that although they sometimes appear simultaneously during the process of LLT, LLT itself is driven by the formation of LFS and not by that of microcrystallites. We also discover that LFS are destroyed upon crystallization, clearly indicating not only that these two types of orderings are competing with each other but also that LFS is a structure unique to the liquid state. Our findings provide a comprehensive view on the long-standing controversy on the origin of the glacial phase, which was discovered by Kivelson and his coworkers (29, 30), and show that the fraction of LFS may be the relevant order parameter of LLT. This suggests that a liquid can have a spatiotemporal hierarchical structure at a low temperature, contrary to the common picture of a high-temperature liquid where the structure is random and homogeneous beyond the molecular size.     

Fuxianhuiid ventral nerve cord and early nervous system evolution in Panarthropoda

Panarthropods are typified by disparate grades of neurological organization reflecting a complex evolutionary history. The fossil record offers a unique opportunity to reconstruct early character evolution of the nervous system via exceptional preservation in extinct representatives. Here we describe the neurological architecture of the ventral nerve cord (VNC) in the upper-stem group euarthropod Chengjiangocaris kunmingensis from the early Cambrian Xiaoshiba Lagerstätte (South China). The VNC of C. kunmingensis comprises a homonymous series of condensed ganglia that extend throughout the body, each associated with a pair of biramous limbs. Submillimetric preservation reveals numerous segmental and intersegmental nerve roots emerging from both sides of the VNC, which correspond topologically to the peripheral nerves of extant Priapulida and Onychophora. The fuxianhuiid VNC indicates that ancestral neurological features of Ecdysozoa persisted into derived members of stem-group Euarthropoda but were later lost in crown-group representatives. These findings illuminate the VNC ground pattern in Panarthropoda and suggest the independent secondary loss of cycloneuralian-like neurological characters in Tardigrada and Euarthropoda.The nervous system represents a critical source of phylogenetic information and has been used extensively for exploring the evolutionary relationships of extant Panarthropoda (i.e., Onychophora, Tardigrada, Euarthropoda) (1–7). Identification of fossilized nervous tissues has provided a unique perspective on early euarthropod brain neuroanatomy and suggests that broad patterns of extant neurological diversity were already in place by the Cambrian (8–11). The ventral nerve cord (VNC) reflects fundamental aspects of panarthropod body organization that complement the organization of the brain and together illuminate the evolution of the CNS (1–3, 5, 7, 12–16). The early evolutionary history of the panarthropod postcephalic CNS, however, remains obscure due to the exclusive preservation of brains in most available fossils (8, 10, 11). Moreover, the unresolved phylogenetic relationships within Panarthropoda complicate accurate reconstruction of the CNS ground pattern (16–22). In this study, we demonstrate the exceptional preservation of postcephalic neurological features in the early Cambrian fuxianhuiid Chengjiangocaris kunmingensis, an upper stem-group euarthropod (17) from the Xiaoshiba Lagerstätte, South China (23). These fossils clarify the neurological organization of the VNC in early euarthropod ancestors, thereby polarizing the evolution of the panarthropod CNS.     

Kinetics of protein–ligand unbinding: Predicting pathways,rates, and rate-limiting steps

The ability to predict the mechanisms and the associated rate constants of protein–ligand unbinding is of great practical importance in drug design. In this work we demonstrate how a recently introduced metadynamics-based approach allows exploration of the unbinding pathways, estimation of the rates, and determination of the rate-limiting steps in the paradigmatic case of the trypsin–benzamidine system. Protein, ligand, and solvent are described with full atomic resolution. Using metadynamics, multiple unbinding trajectories that start with the ligand in the crystallographic binding pose and end with the ligand in the fully solvated state are generated. The unbinding rate k_off is computed from the mean residence time of the ligand. Using our previously computed binding affinity we also obtain the binding rate k_on. Both rates are in agreement with reported experimental values. We uncover the complex pathways of unbinding trajectories and describe the critical rate-limiting steps with unprecedented detail. Our findings illuminate the role played by the coupling between subtle protein backbone fluctuations and the solvation by water molecules that enter the binding pocket and assist in the breaking of the shielded hydrogen bonds. We expect our approach to be useful in calculating rates for general protein–ligand systems and a valid support for drug design.Understanding the thermodynamics and kinetics of protein–ligand interactions is of paramount relevance in the early stages of drug discovery (1–3). So far the major emphasis has been placed on predicting the most likely binding pose as determined by the highest binding affinity (4, 5). In contrast, it has not been possible to predict the pathways for unbinding and the associated rates. However, it is by now well-recognized that one of the most pertinent factors for sustained drug efficacy and safety is not just its affinity, but possibly even more so, the mean lifetime of the protein–ligand complex (1–3). The latter property is strictly related to the time during which the ligand remains in the binding site (1, 2), and is typically expressed by its inverse, the dissociation rate k_off (2). In principle k_off should be amenable to calculations through all-atom molecular dynamics (MD) simulations. These simulations could give detailed and useful insights into the atomic interactions at work during unbinding, especially in the ephemeral but kinetically most relevant transition state ensemble (TSE) (6, 7). Such information is of great value in designing modifications of the ligand that might improve its pharmaceutical properties.However, despite the potential of MD simulations no such calculation has yet been reported. This is a consequence of the limited timescales of MD simulations. Even with the most modern purpose-built supercomputers or massive distributed computing, one can barely reach the timescale of milliseconds (3). Unfortunately most of the reported ligand–protein dissociation times far exceed this timescale (2). These timescales can be reached either by transition path sampling methods (8, 9), quasi-classical approximations (10), by the construction of Markov state models (11, 12), or through carefully designed enhanced sampling methods (8, 13–30) that make accessible the timescale of seconds and beyond in a controlled and accurate way. The enhanced sampling method we use in this work is based on metadynamics (13–15), which has been widely and successfully applied to a variety of systems including complex protein–ligand systems (25–30), and has been rigorously proven to converge to the correct free-energy surface (31, 32).Recently, we have extended the scope of metadynamics by showing that it can also be used to recover kinetic information (15). Furthermore, we showed that by using an a posteriori statistical analysis (33) one can also establish the reliability of the kinetics thus generated. The use of metadynamics for obtaining kinetic information is still in its infancy, however its usefulness has been tested by us and other groups in a range of systems (15, 33–36).In this work, we demonstrate that the scope of the method reported in ref. 15 can be extended to study protein–ligand dissociation pathways and to determine in an accurate way the ligand unbinding rates. We reach well into the hundreds of milliseconds regime and longer, maintaining at the same time full atomic resolution for protein, ligand, and solvent. Specifically, we study the unbinding of the inhibitor benzamidine from trypsin, a serine protease protein (27, 37, 38) using classical force fields (39, 40). Using our acceleration method (15, 33) we are able to harness 21 independent successful unbinding trajectories in which the ligand goes from the bound to the fully unbound state. We find that one of the most distinctive features of the unbinding process is the role played by the water molecules (41, 42). In particular, the solvent promotes unbinding by assisting in the breakage of shielded hydrogen bonds through the formation of water bridge interactions (41).From the analysis of the unbinding trajectories we find that along the unbinding pathways the ligand rests for times ranging from nanoseconds to milliseconds in a number of intermediate structures. We calculate the rates for all possible transitions between these intermediates and construct a Markov model for the unbinding process (11, 43, 44). The overall escape rate computed from this Markov model is in good agreement with the direct estimation of the mean unbinding time that comes from the metadynamics runs. Reassured by this agreement we use the Markov model to determine the dominant unbinding pathways and rate-limiting steps. To this end, starting from the metadynamics reactive trajectories, we perform a committor analysis and determine the TSE (6). Using the recently computed value of the binding affinity (27) we also estimate the binding rate constant k_on. Our calculated unbinding and binding rates compare reasonably well with the known experimental measurement (37), especially taking into account the margin of error in the experiment and the inaccuracy of the force field used in the simulations (42). Unprecedented structural features of the target are also disclosed. In particular, we find that in its apo state trypsin can exist in two forms. In the first form, loop Val207–Tyr224 (hereafter labeled loop L) oscillates around the crystallographic state. In the other form, a small distortion of this loop is stabilized. The mean lifetime of this distorted state is nearly 0.7 ms and during this time the ligand cannot reach the binding site.We believe that this metadynamics-based strategy is, to our knowledge, the first direct approach for calculating k_off from MD simulations of unbinding. Previous studies have focused on the calculation of k_on and the magnitude of k_off was only indirectly obtained (12, 38). Our strategy should be easily applicable for calculating unbinding pathways and rates for generic protein–ligand systems, thus complementing and extending the role of enhanced sampling-based simulations in drug discovery.     

Hyperpolarization-activated,cyclic nucleotide-gated cation channels in Aplysia: Contribution to classical conditioning

Hyperpolarization-activated, cyclic nucleotide-gated cation (HCN) channels are critical regulators of neuronal excitability, but less is known about their possible roles in synaptic plasticity and memory circuits. Here, we characterized the HCN gene organization, channel properties, distribution, and involvement in associative and nonassociative forms of learning in Aplysia californica. Aplysia has only one HCN gene, which codes for a channel that has many similarities to the mammalian HCN channel. The cloned acHCN gene was expressed in Xenopus oocytes, which displayed a hyperpolarization-induced inward current that was enhanced by cGMP as well as cAMP. Similarly to its homologs in other animals, acHCN is permeable to K⁺ and Na⁺ ions, and is selectively blocked by Cs⁺ and ZD7288. We found that acHCN is predominantly expressed in inter- and motor neurons, including LFS siphon motor neurons, and therefore tested whether HCN channels are involved in simple forms of learning of the siphon-withdrawal reflex in a semiintact preparation. ZD7288 (100 μM) significantly reduced an associative form of learning (classical conditioning) but had no effect on two nonassociative forms of learning (intermediate-term sensitization and unpaired training) or baseline responses. The HCN current is enhanced by nitric oxide (NO), which may explain the postsynaptic role of NO during conditioning. HCN current in turn enhances the NMDA-like current in the motor neurons, suggesting that HCN channels contribute to conditioning through this pathway.Hyperpolarization-activated, cyclic nucleotide-gated (HCN), cation nonselective ion channels generate hyperpolarization-activated inward currents (I_h) and thus tend to stabilize membrane potential (1–3). In addition, binding of cyclic nucleotides (cAMP and cGMP) to the C-terminal cyclic nucleotide binding domain (CNBD) enhances I_h and thus couples membrane excitability with intracellular signaling pathways (2, 4). HCN channels are widely important for numerous systemic functions such as hormonal regulation, heart contractility, epilepsy, pain, central pattern generation, sensory perception (4–15), and learning and memory (16–24).However, in previous studies it has been difficult to relate the cellular effects of HCN channels directly to their behavioral effects, because of the immense complexity of the mammalian brain. We have therefore investigated the role of HCN channels in Aplysia, which has a numerically simpler nervous system (25). We first identified and characterized an HCN gene in Aplysia, and showed that it codes for a channel that has many similarities to the mammalian HCN channel. We found that the Aplysia HCN channel is predominantly expressed in motor neurons including LFS neurons in the siphon withdrawal reflex circuit (26, 27). We therefore investigated simple forms of learning of that reflex in a semiintact preparation (28–30) and found that HCN current is involved in classical conditioning and enhances the NMDA-like current in the motor neurons. These results provide a direct connection between HCN channels and behavioral learning and suggest a postsynaptic mechanism of that effect. HCN current in turn is enhanced by nitric oxide (NO), a transmitter of facilitatory interneurons, and thus may contribute to the postsynaptic role of NO during conditioning.