A senescence restriction point acting on chromatin integrates oncogenic signals.

We identify a senescence restriction point (SeRP) as a critical event for cells to commit to senescence. The SeRP integrates the intensity and duration of oncogenic stress, keeps a memory of previous stresses, and combines oncogenic signals acting on different pathways by modulating chromatin accessibility. Chromatin regions opened upon commitment to senescence are enriched in nucleolar-associated domains, which are gene-poor regions enriched in repeated sequences. Once committed to senescence, cells no longer depend on the initial stress signal and exhibit a characteristic transcriptome regulated by a transcription factor network that includes ETV4, RUNX1, OCT1, and MAFB. Consistent with a tumor suppressor role for this network, the levels of ETV4 and RUNX1 are very high in benign lesions of the pancreas but decrease dramatically in pancreatic ductal adenocarcinomas. The discovery of senescence commitment and its chromatin-linked regulation suggests potential strategies for reinstating tumor suppression in human cancers.

The Combination of IFN ß and TNF Induces an Antiviral and Immunoregulatory Program via Non-Canonical Pathways Involving STAT2 and IRF9.

Interferon (IFN) ß and Tumor Necrosis Factor (TNF) are key players in immunity against viruses. Compelling evidence has shown that the antiviral and inflammatory transcriptional response induced by IFNß is reprogrammed by crosstalk with TNF. IFNß mainly induces interferon-stimulated genes by the Janus kinase (JAK)/signal transducer and activator of transcription (STAT) pathway involving the canonical ISGF3 transcriptional complex, composed of STAT1, STAT2, and IRF9. The signaling pathways engaged downstream of the combination of IFNß and TNF remain elusive, but previous observations suggested the existence of a response independent of STAT1. Here, using genome-wide transcriptional analysis by RNASeq, we observed a broad antiviral and immunoregulatory response initiated in the absence of STAT1 upon IFNß and TNF costimulation. Additional stratification of this transcriptional response revealed that STAT2 and IRF9 mediate the expression of a wide spectrum of genes. While a subset of genes was regulated by the concerted action of STAT2 and IRF9, other gene sets were independently regulated by STAT2 or IRF9. Collectively, our data supports a model in which STAT2 and IRF9 act through non-canonical parallel pathways to regulate distinct pool of antiviral and immunoregulatory genes in conditions with elevated levels of both IFNß and TNF.

Isolation and Analysis of Rare Norovirus Recombinants from Coinfected Mice Using Drop-Based Microfluidics.

UNLABELLED: Human noroviruses (HuNoVs) are positive-sense RNA viruses that can cause severe, highly infectious gastroenteritis. HuNoV outbreaks are frequently associated with recombination between circulating strains. Strain genotyping and phylogenetic analyses show that noroviruses often recombine in a highly conserved region near the junction of the viral polyprotein (open reading frame 1 [ORF1]) and capsid (ORF2) genes and occasionally within the RNA-dependent RNA polymerase (RdRP) gene. Although genotyping methods are useful for tracking changes in circulating viral populations, they report only the dominant recombinant strains and do not elucidate the frequency or range of recombination events. Furthermore, the relatively low frequency of recombination in RNA viruses has limited studies to cell culture or in vitro systems, which do not reflect the complexities and selective pressures present in an infected organism. Using two murine norovirus (MNV) strains to model coinfection, we developed a microfluidic platform to amplify, detect, and recover individual recombinants following in vitro and in vivo coinfection. One-step reverse transcriptase PCR (RT-PCR) was performed in picoliter drops with primers that identified the wild-type and recombinant progenies and scanned for recombination breakpoints at ∼1-kb intervals. We detected recombination between MNV strains at multiple loci spanning the viral protease, RdRP, and capsid ORFs and isolated individual recombinant RNA genomes that were present at a frequency of 1/300,000 or higher. This study is the first to examine norovirus recombination following coinfection of an animal and suggests that the exchange of RNA among viral genomes in an infected host occurs in multiple locations and is an important driver of genetic diversity. IMPORTANCE: RNA viruses increase diversity and escape host immune barriers by genomic recombination. Studies using a number of viral systems indicate that recombination occurs via template switching by the virus-encoded RNA-dependent RNA polymerase (RdRP). However, factors that govern the frequency and positions of recombination in an infected organism remain largely unknown. This work leverages advances in the applied physics of drop-based microfluidics to isolate and sequence rare recombinants arising from the coinfection of mice with two distinct strains of murine norovirus. This study is the first to detect and analyze norovirus recombination in an animal model.

The Streptomyces-produced antibiotic fosfomycin is a promiscuous substrate for archaeal isopentenyl phosphate kinase.

Isopentenyl phosphate kinase (IPK) catalyzes the phosphorylation of isopentenyl phosphate to form the isoprenoid precursor isopentenyl diphosphate in the archaeal mevalonate pathway. This enzyme is highly homologous to fosfomycin kinase (FomA), an antibiotic resistance enzyme found in a few strains of Streptomyces and Pseudomonas whose mode of action is inactivation by phosphorylation. Superposition of Thermoplasma acidophilum (THA) IPK and FomA structures aligns their respective substrates and catalytic residues, including H50 and K14 in THA IPK and H58 and K18 in Streptomyces wedmorensis FomA. These residues are conserved only in the IPK and FomA members of the phosphate subdivision of the amino acid kinase family. We measured the fosfomycin kinase activity of THA IPK [K(m) = 15.1 ± 1.0 mM, and k(cat) = (4.0 ± 0.1) × 10⁻² s⁻¹], resulting in a catalytic efficiency (k(cat)/K(m) = 2.6 M⁻¹ s⁻¹) that is 5 orders of magnitude lower than that of the native reaction. Fosfomycin is a competitive inhibitor of IPK (K(i) = 3.6 ± 0.2 mM). Molecular dynamics simulation of the IPK·fosfomycin·MgATP complex identified two binding poses for fosfomycin in the IP binding site, one of which results in a complex analogous to the native IPK·IP·ATP complex that engages H50 and the lysine triangle formed by K5, K14, and K205. The other binding pose leads to a dead-end complex that engages K204 near the IP binding site to bind fosfomycin. Our findings suggest a mechanism for acquisition of FomA-based antibiotic resistance in fosfomycin-producing organisms.

Structural basis for µ-opioid receptor binding and activation.

Opioids that stimulate the µ-opioid receptor (MOR1) are the most frequently prescribed and effective analgesics. Here we present a structural model of MOR1. Molecular dynamics simulations show a ligand-dependent increase in the conformational flexibility of the third intracellular loop that couples with the G protein complex. These simulations likewise identified residues that form frequent contacts with ligands. We validated the binding residues using site-directed mutagenesis coupled with radioligand binding and functional assays. The model was used to blindly screen a library of ∼1.2 million compounds. From the 34 compounds predicted to be strong binders, the top three candidates were examined using biochemical assays. One compound showed high efficacy and potency. Post hoc testing revealed this compound to be nalmefene, a potent clinically used antagonist, thus further validating the model. In summary, the MOR1 model provides a tool for elucidating the structural mechanism of ligand-initiated cell signaling and for screening novel analgesics.

A physical model reveals the mechanochemistry responsible for dynein's processive motion.

The molecular motor dynein is associated with various cellular activities, such as directed transport along microtubules and the rhythmic beating of the axoneme. Because of the size and complexity of the protein, a detailed understanding of the mechanochemistry that drives dynein's processive motion is lacking. To overcome this deficiency, we developed the first (to our knowledge) computational model for two-headed dynein that couples conformational changes of the motor's subunits to the biochemical steps involved in ATP hydrolysis. Analysis of the model provides what we believe are several novel insights into how the protein functions: 1), structural constraints limit the motion of the free microtubule binding domain to one dimension, increasing the efficiency with which this domain finds a binding site; 2), in addition to the power stroke of the bound head, recovery of the free head to a pre-power-stroke conformation is required for this head to reach a forward binding site; 3), the order in which the power stroke and recovery transitions occur affects the probability of back-stepping; and 4), the existence of multiple equilibria in the motor's bending energy provides a mechanism for processive back-stepping. To the best of our knowledge, our computational model provides the first complete mechanochemical description of the motor protein dynein, and the findings presented here should motivate new experimental investigations to test its predictions.

Molecular modeling tools and approaches for CFTR and cystic fibrosis.

Cystic fibrosis is a multi-faceted disease resulting from the dysfunction of the CFTR channel. Understanding the structural basis of channel function and the structural origin of the defect is imperative in the development of therapeutic strategies. Here, we describe molecular modeling tools that, in conjunction with complementary experimental tools, lead to significant findings on CFTR channel function and on the effect of the pathogenic mutant F508del.

Multiscale approaches for studying energy transduction in dynein.

Cytoplasmic dynein is an important motor that drives all minus-end directed movement along microtubules. Dynein is a complex motor whose processive motion is driven by ATP-hydrolysis. Dynein's run length has been measured to be several millimetres with typical velocities in the order of a few nanometres per second. Therefore, the average time between steps is a fraction of a second. When this time scale is compared with typical time scales for protein side chain and backbone movements (approximately 10(-9) s and approximately 10(-5) s, respectively), it becomes clear that a multi-timescale modelling approach is required to understand energy transduction in this protein. Here, we review recent efforts to use computational and mathematical modelling to understand various aspects of dynein's chemomechanical cycle. First, we describe a structural model of dynein's motor unit showing a heptameric organization of the motor subunits. Second, we describe our molecular dynamics simulations of the motor unit that are used to investigate the dynamics of the various motor domains. Third, we present a kinetic model of the coordination between the two dynein heads. Lastly, we investigate the various potential geometries of the dimer during its hydrolytic and stepping cycle.

Kinetic models for the coordinated stepping of cytoplasmic dynein.

To generate processive motion along a polymer track requires that motor proteins couple their ATP hydrolysis cycle with conformational changes in their structural subunits. Numerous experimental and theoretical efforts have been devoted to establishing how this chemomechanical coupling occurs. However, most processive motors function as dimers. Therefore a full understanding of the motor's performance also requires knowledge of the coordination between the chemomechanical cycles of the two heads. We consider a general two-headed model for cytoplasmic dynein that is built from experimental measurements on the chemomechanical states of monomeric dynein. We explore different possible scenarios of coordination that simultaneously satisfy two main requirements of the dimeric protein: high processivity (long run length) and high motor velocity (fast ATP turnover). To demonstrate the interplay between these requirements and the necessity for coordination, we first develop and analyze a simple mechanical model for the force-induced stepping in the absence of ATP. Next we use a simplified model of dimeric dynein's chemomechanical cycle to establish the kinetic rules that must be satisfied for the model to be consistent with recent data for the motor's performance from single molecule experiments. Finally, we use the results of these investigations to develop a full model for dimeric dynein's chemomechanical cycle and analyze this model to make experimentally testable predictions.

Identification and rational redesign of peptide ligands to CRIP1, a novel biomarker for cancers.

Cysteine-rich intestinal protein 1 (CRIP1) has been identified as a novel marker for early detection of cancers. Here we report on the use of phage display in combination with molecular modeling to identify a high-affinity ligand for CRIP1. Panning experiments using a circularized C7C phage library yielded several consensus sequences with modest binding affinities to purified CRIP1. Two sequence motifs, A1 and B5, having the highest affinities for CRIP1, were chosen for further study. With peptide structure information and the NMR structure of CRIP1, the higher-affinity A1 peptide was computationally redesigned, yielding a novel peptide, A1M, whose affinity was predicted to be much improved. Synthesis of the peptide and saturation and competitive binding studies demonstrated approximately a 10-28-fold improvement in the affinity of A1M compared to that of either A1 or B5 peptide. These techniques have broad application to the design of novel ligand peptides.

Multiple membrane-cytoplasmic domain contacts in the cystic fibrosis transmembrane conductance regulator (CFTR) mediate regulation of channel gating.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a unique ATP-binding cassette (ABC) ion channel mutated in patients with cystic fibrosis. The most common mutation, deletion of phenylalanine 508 (DeltaF508) and many other disease-associated mutations occur in the nucleotide binding domains (NBD) and the cytoplasmic loops (CL) of the membrane-spanning domains (MSD). A recently constructed computational model of the CFTR three-dimensional structure, supported by experimental data (Serohijos, A. W., Hegedus, T., Aleksandrov, A. A., He, L., Cui, L., Dokholyan, N. V., and Riordan, J. R. (2008) Proc. Natl. Acad. Sci. U. S. A. 105, 3256-3261) revealed that several of these mutations including DeltaF508 disrupted interfaces between these domains. Here we have used cysteine cross-linking experiments to verify all NBD/CL interfaces predicted by the structural model and observed that their cross-linking has a variety of different effects on channel gating. The interdomain contacts comprise aromatic clusters important for stabilization of the interfaces and also involve the Q-loops and X-loops that are in close proximity to the ATP binding sites. Cross-linking of all domain-swapping contacts between NBDs and MSD cytoplasmic loops in opposite halves of the protein rapidly and reversibly arrest single channel gating while those in the same halves have lesser impact. These results reinforce the idea that mediation of regulatory signals between cytoplasmic- and membrane-integrated domains of the CFTR channel apparently relies on an array of precise but highly dynamic interdomain structural joints.

Diminished self-chaperoning activity of the DeltaF508 mutant of CFTR results in protein misfolding.

The absence of a functional ATP Binding Cassette (ABC) protein called the Cystic Fibrosis Transmembrane Conductance Regulator (CFTR) from apical membranes of epithelial cells is responsible for cystic fibrosis (CF). Over 90% of CF patients carry at least one mutant allele with deletion of phenylalanine at position 508 located in the N-terminal nucleotide binding domain (NBD1). Biochemical and cell biological studies show that the DeltaF508 mutant exhibits inefficient biosynthetic maturation and susceptibility to degradation probably due to misfolding of NBD1 and the resultant misassembly of other domains. However, little is known about the direct effect of the Phe508 deletion on the NBD1 folding, which is essential for rational design strategies of cystic fibrosis treatment. Here we show that the deletion of Phe508 alters the folding dynamics and kinetics of NBD1, thus possibly affecting the assembly of the complete CFTR. Using molecular dynamics simulations, we find that meta-stable intermediate states appearing on wild type and mutant folding pathways are populated differently and that their kinetic accessibilities are distinct. The structural basis of the increased misfolding propensity of the DeltaF508 NBD1 mutant is the perturbation of interactions in residue pairs Q493/P574 and F575/F578 found in loop S7-H6. As a proof-of-principle that the S7-H6 loop conformation can modulate the folding kinetics of NBD1, we virtually design rescue mutations in the identified critical interactions to force the S7-H6 loop into the wild type conformation. Two redesigned NBD1-DeltaF508 variants exhibited significantly higher folding probabilities than the original NBD1-DeltaF508, thereby partially rescuing folding ability of the NBD1-DeltaF508 mutant. We propose that these observed defects in folding kinetics of mutant NBD1 may also be modulated by structures separate from the 508 site. The identified structural determinants of increased misfolding propensity of NBD1-DeltaF508 are essential information in correcting this pathogenic mutant.

Computational studies reveal phosphorylation-dependent changes in the unstructured R domain of CFTR.

The cystic fibrosis transmembrane conductance regulator (CFTR) is a cAMP-dependent chloride channel that is mutated in cystic fibrosis, an inherited disease of high morbidity and mortality. The phosphorylation of its approximately 200 amino acid R domain by protein kinase A is obligatory for channel gating under normal conditions. The R domain contains more than ten PKA phosphorylation sites. No individual site is essential but phosphorylation of increasing numbers of sites enables progressively greater channel activity. In spite of numerous studies of the role of the R domain in CFTR regulation, its mechanism of action remains largely unknown. This is because neither its structure nor its interactions with other parts of CFTR have been completely elucidated. Studies have shown that the R domain lacks well-defined secondary structural elements and is an intrinsically disordered region of the channel protein. Here, we have analyzed the disorder pattern and employed computational methods to explore low-energy conformations of the R domain. The specific disorder and secondary structure patterns detected suggest the presence of molecular recognition elements (MoREs) that may mediate phosphorylation-regulated intra- and inter-domain interactions. Simulations were performed to generate an ensemble of accessible R domain conformations. Although the calculated structures may represent more compact conformers than occur in vivo, their secondary structure propensities are consistent with predictions and published experimental data. Equilibrium simulations of a mimic of a phosphorylated R domain showed that it exhibited an increased radius of gyration. In one possible interpretation of these findings, by changing its size, the globally unstructured R domain may act as an entropic spring to perturb the packing of membrane-spanning sequences that constitute the ion permeability pathway and thereby activate channel gating.

Phenylalanine-508 mediates a cytoplasmic-membrane domain contact in the CFTR 3D structure crucial to assembly and channel function.

Deletion of phenylalanine-508 (Phe-508) from the N-terminal nucleotide-binding domain (NBD1) of the cystic fibrosis transmembrane conductance regulator (CFTR), a member of the ATP-binding cassette (ABC) transporter family, disrupts both its folding and function and causes most cystic fibrosis. Most mutant nascent chains do not pass quality control in the ER, and those that do remain thermally unstable, only partially functional, and are rapidly endocytosed and degraded. Although the lack of the Phe-508 peptide backbone diminishes the NBD1 folding yield, the absence of the aromatic side chain is primarily responsible for defective CFTR assembly and channel gating. However, the site of interdomain contact by the side chain is unknown as is the high-resolution 3D structure of the complete protein. Here we present a 3D structure of CFTR, constructed by molecular modeling and supported biochemically, in which Phe-508 mediates a tertiary interaction between the surface of NBD1 and a cytoplasmic loop (CL4) in the C-terminal membrane-spanning domain (MSD2). This crucial cytoplasmic membrane interface, which is dynamically involved in regulation of channel gating, explains the known sensitivity of CFTR assembly to many disease-associated mutations in CL4 as well as NBD1 and provides a sharply focused target for small molecules to treat CF. In addition to identifying a key intramolecular site to be repaired therapeutically, our findings advance understanding of CFTR structure and function and provide a platform for focused biochemical studies of other features of this unique ABC ion channel.

A structural model reveals energy transduction in dynein.

Intracellular active transport is driven by ATP-hydrolyzing motor proteins that move along cytoskeletal filaments. In particular, the microtubule-associated dynein motor is involved in the transport of organelles and vesicles, the maintenance of the Golgi, and mitosis. However, unlike kinesin and myosin, the mechanism by which dynein converts chemical energy into mechanical force remains largely a mystery, due primarily to the lack of a high-resolution molecular structure. Using homology modeling and normal mode analysis, we propose a complete atomic structure and a mechanism for force generation by the motor protein dynein. In agreement with very recent electron microscopy (EM) reconstructions showing dynein as a ring-shaped heptamer, our model consists of six ATPases of the AAA (ATPases associated with various cellular activities) superfamily and a C-terminal domain, which is experimentally known to control motor function. Our model shows a coiled coil spanning the diameter of the motor that accounts for previously unidentified structures in EM studies and provides a potential mechanism for long-range communication between the AAA domains. Furthermore, normal mode analysis reveals that the subunits of the motor that contain the nucleotide binding sites exhibit minimal movement, whereas the rest of the motor is very mobile. Our analysis suggests the likely domain rearrangements of the motor unit that generate its power stroke. This study provides insights into the structure and function of dynein that can guide further experimental investigations into energy transduction in dynein.