SwissParam: a fast force field generation tool for small organic molecules.

The drug discovery process has been deeply transformed recently by the use of computational ligand-based or structure-based methods, helping the lead compounds identification and optimization, and finally the delivery of new drug candidates more quickly and at lower cost. Structure-based computational methods for drug discovery mainly involve ligand-protein docking and rapid binding free energy estimation, both of which require force field parameterization for many drug candidates. Here, we present a fast force field generation tool, called SwissParam, able to generate, for arbitrary small organic molecule, topologies, and parameters based on the Merck molecular force field, but in a functional form that is compatible with the CHARMM force field. Output files can be used with CHARMM or GROMACS. The topologies and parameters generated by SwissParam are used by the docking software EADock2 and EADock DSS to describe the small molecules to be docked, whereas the protein is described by the CHARMM force field, and allow them to reach success rates ranging from 56 to 78%. We have also developed a rapid binding free energy estimation approach, using SwissParam for ligands and CHARMM22/27 for proteins, which requires only a short minimization to reproduce the experimental binding free energy of 214 ligand-protein complexes involving 62 different proteins, with a standard error of 2.0 kcal mol(-1), and a correlation coefficient of 0.74. Together, these results demonstrate the relevance of using SwissParam topologies and parameters to describe small organic molecules in computer-aided drug design applications, together with a CHARMM22/27 description of the target protein. SwissParam is available free of charge for academic users at www.swissparam.ch.

Functional analysis of human POT1 and RPA : insights into the role of single stranded DNA binding proteins at telomeres

Abstract Telomeres, the natural ends of chromosomes, need to be protected from chromosome end fusions, aberrant homologous recombination and degradation. In humans, chromosome ends are specified through arrays of tandemly repeated 5'-TTAGGG-3' hexamers, ending in a 3' overhang. A complex formed by the six proteins TRF1, TRF2, hRap1, TIN2, TPP1 and POT1 specifically assocìates with and protects telomeres. Telomeres are maintained by semiconservative DNA replication and by a specialized reverse transcriptase, telomerase, that carries an RNA subunit which templates new telomeric repeat synthesis. The telomeric single stranded (ss) DNA binding protein POT1 protects the telomeric 3' overhang and modulates telomerase-mediated telomere elongation. It is possible that POT1 also influences DNA synthesis during semiconservative DNA replication, which is initiated by the DNA polymerase alpha-primase complex. The heterotrimeric ss DNA-binding protein RPA plays essential roles during DNA replication. RPA binds to ss DNA with high affinity in order to stabilize ss DNA and facilitate nascent strand synthesis at the replication fork. Here we investigate how the two proteins RPA and POT1 contribute to telomere maintenance by regulating semi-conservative DNA replication and telomerase. Using chromatin immunoprecipitation experiments, we show that RPA associates with telomeres during S-phase. Analysis of telomere structure in cells shRNA-depleted for RPA and POT1 reveals that loss of RPA and POT1 causes exposure of single-stranded DNA at telomeres, suggestive of incomplete DNA replication. Biochemical experiments using purified recombinant POT1 and RPA show that saturating telomeric oligonucleotides with POT1 or RPA reduces the primase activity of the DNA polymerase alpha-primase complex and the overall activity of telomerase. POT1 and RPA also increase the primer extension by DNA polymerase alpha-primase complex and the processivity of telomerase under certain conditions, although POT1 increases the activities to a greater extent than RPA. We propose that POT1 is required for proper replication of the lagging strand of telomeres and that some phenotypes observed in POT1-depleted cells may stern from incomplete DNA replication rather than de-protection of the single-stranded overhang. Résumé Les télomères, les extrémités normales des chromosomes linéaires, doivent être protégés des fusions chromosomiques, d'événements de recombinaison homologue aberrants et de phénomènes de dégradation. Chez l'Homme, les extrémités des chromosomes sont constitués d'ADN double brin répétitif de séquence 5'-TTAGGG-3', d'une extension simple brin 3' sortante et d'un complexe protéique formé des six facteurs TRF1, TRF2, hRap1, TIN2, TPP1 et POT1 qui, s'associant à cette séquence, protègent l'ADN télomèrique. Les télomères sont maintenus par la télomérase, une transcriptase inverse capable d'allonger l'extension 3' sortante télomérique. POT1 lie l'ADN simple brin télomérique et module l'élongation des télomères par la télomérase. POT1 pourrait en théorie également influencer la réplication semi-conservative de l'ADN. L'ADN-polymérase Pal alpha-primase amorce et initie la synthèse d'ADN. Pendant la réplication, l'ADN simple brin est stabilisé par RPA, un complexe hétérotrimèrique qui lie l'ADN simple brin. RPA facilite la synthèse du brin naissant à la fourche de réplication. Ici nous avons étudié comment ces deux protéines qui lient l'ADN simple brin, RPA et POT1, régulent la réplication des télomères par la télomérase et la machinerie classique de réplication de l'ADN. Par immunoprécipitation de chromatine (ChIP), nous montrons que RPA est localisé aux télomères lors de la phase S du cycle cellulaire. De plus, l'analyse de la structure des télomeres indique que !a perte de RPA ou de POT1 conduit à l'apparition d'ADN simple brin télomérique, suggérant une réplication incomplète de l'ADN télomérique in vivo. Par une approche complémentaire biochimique utilisant les protéines POT1 et RPA recombinantes purifiées, nous montrons également que la liaison de POT1 ou de RPA à des oligonucléotides télomériques bloque l'activité primase du complexe polymérase alpha/primase et réduit l'activité télomérase sur ces substrats. En revanche, leur liaison augmente l'activité ADN-polymérase du complexe polymérase alpha/primase, ainsi que fa processivité de la télomérase dans certaines conditions, POT1 étant le plus efficace des deux facteurs. Nous proposons que POT1 est nécessaire à la réplication du brin retardé au niveau des télomères, ce qui suggère que certains phénotypes des cellules déplétés en POT1 puissent résulter d'une réplication incomplète de l'ADN télémétrique plutôt que d'une déprotection de l'extrémité sortante des télomères.

Synthesis of a radioiodinated photoreactive MAGE-1 peptide derivative and photoaffinity labeling of cell-associated human leukocyte antigen-A1 molecules.

The synthesis of a photoreactive derivative of the human leukocyte antigen-A1 (HLA-A1)-restricted MAGE-1 peptide 161-169 (EADPTGHSY) is described. Using conventional automated solid-phase peptide synthesis, a photoreactive derivative of this peptide was synthesized by replacing histidine-167 with photo-reactive N-beta-4-azidosalicyloyl-L-2,3-diaminopropionic acid. The C-terminal tyrosine was incorporated as phosphotyrosine. This peptide derivative was radioiodinated in the presence of chloramine T. This iodination took place selectively at the photoreactive group, because the phosphate ester prevented tyrosine iodination. Following dephosphorylation with alkaline phosphatase and chromatographic purification, the radiolabeled peptide derivative was incubated with cells expressing HLA-A1 or other HLA molecules. Photoactivation resulted in efficient photoaffinity labeling of HLA-A1. Other HLA molecules or other cellular components were not detectably labeled. This labeling was inhibited by HLA-A1 but not by HLA-A2-binding peptides. This synthesis is generally applicable and can also be adapted to the synthesis of well-defined radiolabeled nonphotoreactive peptide derivatives.

A deep dig-hindsight on Holocene plant composition from ancient environmental DNA

Want a glimpse at past vegetation? Studying pollen and other plant remains, which are preserved for example in lake sediments or mires for thousands of years, allows us to document regional occurrences of plant species over radiocarbon-dated time series. Such vegetation reconstructions derived from optical analyses of fossil samples are inherently incomplete because they only comprise taxa that contribute sufficient amounts of pollen, spores, macrofossil or other evidences. To complement optical analyses for paleoecological inference, molecular markers applied to ancient DNA (aDNA) may help in disclosing information hitherto inaccessible to biologists. Parducci etal. (2013) targeted aDNA from sediment cores of two lakes in the Scandes Mountains with generic primers in a meta-barcoding approach. When compared to palynological records from the same cores, respective taxon lists show remarkable differences in their compositions, but also in quantitative representation and in taxonomic resolution similar to a previous study (JOrgensen etal. 2012). While not free of assumptions that need critical and robust testing, notably the question of possible contamination, this study provides thrilling prospects to improve our knowledge about past vegetation composition, but also other organismic groups, stored as a biological treasure in the ground.

Homologous DNA pairing domain peptides of RecA protein: intrinsic propensity to form beta-structures and filaments.

The 20 amino acid residue peptides derived from RecA loop L2 have been shown to be the pairing domain of RecA. The peptides bind to ss- and dsDNA, unstack ssDNA, and pair the ssDNA to its homologous target in a duplex DNA. As shown by circular dichroism, upon binding to DNA the disordered peptides adopt a beta-structure conformation. Here we show that the conformational change of the peptide from random coil to beta-structure is important in binding ss- and dsDNA. The beta-structure in the DNA pairing peptides can be induced by many environmental conditions such as high pH, high concentration, and non-micellar sodium dodecyl sulfate (6 mM). This behavior indicates an intrinsic property of these peptides to form a beta-structure. A beta-structure model for the loop L2 of RecA protein when bound to DNA is thus proposed. The fact that aromatic residues at the central position 203 strongly modulate the peptide binding to DNA and subsequent biochemical activities can be accounted for by the direct effect of the aromatic amino acids on the peptide conformational change. The DNA-pairing domain of RecA visualized by electron microscopy self-assembles into a filamentous structure like RecA. The relevance of such a peptide filamentous structure to the structure of RecA when bound to DNA is discussed.

Molecular characterization of HLA-C incompatibilities in HLA-ABDR-matched unrelated bone marrow donor-recipient pairs. Sequence of two new Cw alleles (Cw*02023 and Cw*0707) and recognition by cytotoxic T lymphocytes.

While the influence of HLA-AB and -DRB1 matching on the outcome of bone marrow transplantation (BMT) with unrelated donors is clear, the evaluation of HLA-C has been hampered by its poor serological definition. Because the low resolution of standard HLA-C typing could explain the significant number of positive cytotoxic T lymphocyte precursor frequency (CTLpf) tests found among HLA-AB-subtype, DRB1/B3/B5-subtype matched patient/donor pairs, we have identified by sequencing the incompatibilities recognized by CD8+ CTL clones obtained from such positive CTLpf tests. In most cases the target molecules were HLA-C antigens that had escaped detection by serology (e.g. Cw*1601, 1502 or 0702). Direct recognition of HLA-C by a CTL clone was demonstrated by lysis of the HLA class I-negative 721.221 cell line transfected with Cw*1601 cDNA. Because of the functional importance of Cw polymorphism, a PCR-SSO oligotyping procedure was set up allowing the resolution of 29 Cw alleles. Oligotyping of a panel of 382 individuals (including 101 patients and their 272 potential unrelated donors, 5 related donors and 4 platelet donors) allowed to determine HLA-C and HLA A-B-Cw-DRB1 allelic frequencies, as well as a number of A-Cw, B-Cw, and DRB1-Cw associations. Two new HLA-Cw alleles (Cw*02023 and Cw*0707) were identified by DNA sequencing of PCR-amplified exon 2-intron 2-exon 3 amplicons. Furthermore, we determined the degree of HLA-C compatibility in 287 matched pairs that could be formed from 73 patients and their 184 potential unrelated donors compatible for HLA-AB by serology and for HLA-DRB1/ B3/B5 by oligotyping. Cw mismatches were identified in 42.1% of these pairs, and AB-subtype oligotyping showed that 30% of these Cw-incompatible pairs were also mismatched for A or B-locus subtype. The degree of HLA-C incompatibility was strongly influenced by the linkage with B alleles and by the ABDR haplotypes. Cw alleles linked with B*4403, B*5101, B18, and B62 haplotypes were frequently mismatched. Apparently high resolution DNA typing for HLA-AB does not result in full matching at locus C. Since HLA-C polymorphism is recognized by alloreactive CTLs, such incompatibilities might be as relevant as AB-subtype mismatches in clinical transplantation.

Electron microscopy of RecA-DNA complexes: Two different states, their functional significance and relation to the solved crystal structure

Seven different electron microscopy techniques habe been employed to study the RecA protein of E. coli. This review provides a summary of the conclusions that have been drawn from these studies, and attempts to relate these observations to models for the role of RecA protein in homologous recombination.

Effect of early antiretroviral therapy during primary HIV-1 infection on cell-associated HIV-1 DNA and plasma HIV-1 RNA.

Background: Early initiation of combination antiretroviral therapy (ART) during primary HIV-1 infection may prevent the establishment of large viral reservoirs, possibly resulting in improved control of plasma viraemia rebound after ART cessation.Methods: Levels of cell-associated HIV-1 DNA and plasma HIV-1 RNA were measured longitudinally in 32 acutely and recently infected patients, who started ART <= 120 days after the estimated date of infection, and interrupted ART after 18 months (median) of continuous therapy. Averages of HIV-1 DNA and RNA concentrations present in blood 30-365 days after therapy interruption (median duration 300 days, range 195-358) were compared between patients who started ART <= 60 days after the estimated date of infection (early starters), those who started between 61 and 120 days (later starters), and, for HIV-1 RNA only, with 89 untreated participants of the Swiss HIV Cohort Study with documented sero-conversion and longitudinal measurements collected 90-455 days after the first positive HIV test.Results: In early ART starters, average levels of plasma HIV-1 RNA and cell-associated HIV-1 DNA after treatment interruption were 1 log(10) (P=0.008) and 0.4 log(10) (P=0.03) lower compared with later starters. Average post-treatment plasma HIV-1 RNA levels in early starters were significantly lower, respectively, compared with untreated controls (-1.2 log(10); P<0.0004).Conclusions: Early treatment initiation within 2 months after HIV infection compared with later therapy initiation resulted in reduced levels of plasma viraemia and proviral HIV-1 DNA for >= 1 year after subsequent ART cessation. Plasma HIV-1 RNA levels in early starters were also significantly lower than in untreated controls.

The efficiency of a novel delivery system (PEI600-Tat) in development of potent DNA vaccine using HPV16 E7 as a model antigen.

DNA vaccination is a promising approach for inducing both humoral and cellular immune responses. The mode of plasmid DNA delivery is critical to make progress in DNA vaccination. Using human papillomavirus type 16 E7 as a model antigen, this study evaluated the effect of peptide-polymer hybrid including PEI600-Tat conjugate as a novel gene delivery system on the potency of antigen-specific immunity in mice model. At ratio of 10:50 PEI-Tat/E7DNA (w/w), both humoral and cellular immune responses were significantly enhanced as compared with E7DNA construct and induced Th1 response. Therefore, this new delivery system could have promising applications in gene therapy.

DNA structure-specific priming of ATR activation by DNA-PKcs.

Three phosphatidylinositol-3-kinase-related protein kinases implement cellular responses to DNA damage. DNA-dependent protein kinase catalytic subunit (DNA-PKcs) and ataxia-telangiectasia mutated respond primarily to DNA double-strand breaks (DSBs). Ataxia-telangiectasia and RAD3-related (ATR) signals the accumulation of replication protein A (RPA)-covered single-stranded DNA (ssDNA), which is caused by replication obstacles. Stalled replication intermediates can further degenerate and yield replication-associated DSBs. In this paper, we show that the juxtaposition of a double-stranded DNA end and a short ssDNA gap triggered robust activation of endogenous ATR and Chk1 in human cell-free extracts. This DNA damage signal depended on DNA-PKcs and ATR, which congregated onto gapped linear duplex DNA. DNA-PKcs primed ATR/Chk1 activation through DNA structure-specific phosphorylation of RPA32 and TopBP1. The synergistic activation of DNA-PKcs and ATR suggests that the two kinases combine to mount a prompt and specific response to replication-born DSBs.

Thiamine-auxotrophic mutants of Pseudomonas fluorescens CHA0 are defective in cell-cell signaling and biocontrol factor expression.

In the biocontrol strain Pseudomonas fluorescens CHA0, the Gac/Rsm signal transduction pathway positively controls the synthesis of antifungal secondary metabolites and exoenzymes. In this way, the GacS/GacA two-component system determines the expression of three small regulatory RNAs (RsmX, RsmY, and RsmZ) in a process activated by the strain's own signal molecules, which are not related to N-acyl-homoserine lactones. Transposon Tn5 was used to isolate P. fluorescens CHA0 insertion mutants that expressed an rsmZ-gfp fusion at reduced levels. Five of these mutants were gacS negative, and in them the gacS mutation could be complemented for exoproduct and signal synthesis by the gacS wild-type allele. Furthermore, two thiamine-auxotrophic (thiC) mutants that exhibited decreased signal synthesis in the presence of 5 x 10(-8) M thiamine were found. Under these conditions, a thiC mutant grew normally but showed reduced expression of the three small RNAs, the exoprotease AprA, and the antibiotic 2,4-diacetylphloroglucinol. In a gnotobiotic system, a thiC mutant was impaired for biological control of Pythium ultimum on cress. Addition of excess exogenous thiamine restored all deficiencies of the mutant. Thus, thiamine appears to be an important factor in the expression of biological control by P. fluorescens.

DNA binding properties of peroxisome proliferator-activated receptor subtypes on various natural peroxisome proliferator response elements. Importance of the 5'-flanking region.

The three subtypes of the peroxisome proliferator-activated receptors (PPARalpha, beta/delta, and gamma) form heterodimers with the 9-cis-retinoic acid receptor (RXR) and bind to a common consensus response element, which consists of a direct repeat of two hexanucleotides spaced by one nucleotide (DR1). As a first step toward understanding the molecular mechanisms determining PPAR subtype specificity, we evaluated by electrophoretic mobility shift assays the binding properties of the three PPAR subtypes, in association with either RXRalpha or RXRgamma, on 16 natural PPAR response elements (PPREs). The main results are as follows. (i) PPARgamma in combination with either RXRalpha or RXRgamma binds more strongly than PPARalpha or PPARbeta to all natural PPREs tested. (ii) The binding of PPAR to strong elements is reinforced if the heterodimerization partner is RXRgamma. In contrast, weak elements favor RXRalpha as heterodimerization partner. (iii) The ordering of the 16 natural PPREs from strong to weak elements does not depend on the core DR1 sequence, which has a relatively uniform degree of conservation, but correlates with the number of identities of the 5'-flanking nucleotides with respect to a consensus element. This 5'-flanking sequence is essential for PPARalpha binding and thus contributes to subtype specificity. As a demonstration of this, the PPARgamma-specific element ARE6 PPRE is able to bind PPARalpha only if its 5'-flanking region is exchanged with that of the more promiscuous HMG PPRE.

Respective roles of calcitonin receptor-like receptor (CRLR) and receptor activity-modifying proteins (RAMP) in cell surface expression of CRLR/RAMP heterodimeric receptors.

Receptor activity modifying proteins RAMP1, RAMP2, and RAMP3 are responsible for defining affinity to ligands of the calcitonin receptor-like receptor (CRLR). It has also been proposed that receptor activity-modifying proteins (RAMP) are molecular chaperones required for CRLR transport to the cell surface. Here, we have studied the respective roles of CRLR and RAMP in transporting CRLR/RAMP heterodimers to the plasma membrane by using a highly specific binding assay that allows quantitative detection of cell surface-expressed CRLR or RAMP in the Xenopus oocytes expression system. We show that: (i) heterodimer assembly is not a prerequisite for efficient cell surface expression of CRLR, (ii) N-glycosylated RAMP2 and RAMP3 are expressed at the cell surface and their transport to the plasma membrane requires N-glycans, (iii) RAMP1 is not N-glycosylated and is transported to the plasma membrane only upon formation of heterodimers with CRLR, and (iv) introduction of N-glycosylation sites in the RAMP1 sequence (D58N/G60S, Y71N, and K103N/P105S) allows cell surface expression of these mutants at levels similar to that of wild-type RAMP1 co-expressed with CRLR. Our data argue against a chaperone function for RAMP and identify the role of N-glycosylation in targeting these molecules to the cell surface.

Dissociation of thymic positive and negative selection in transgenic mice expressing major histocompatibility complex class I molecules exclusively on thymic cortical epithelial cells.

Thymic positive and negative selection of developing T lymphocytes confronts us with a paradox: How can a T-cell antigen receptor (TCR)-major histocompatibility complex (MHC)/peptide interaction in the former process lead to transduction of signals allowing for cell survival and in the latter induce programmed cell death or a hyporesponsive state known as anergy? One of the hypotheses put forward states that the outcome of a TCR-MHC/peptide interaction depends on the cell type presenting the selecting ligand to the developing thymocyte. Here we describe the development and lack of self-tolerance of CD8(+) T lymphocytes in transgenic mice expressing MHC class I molecules in the thymus exclusively on cortical epithelial cells. Despite the absence of MHC class I expression on professional antigen-presenting cells, normal numbers of CD8(+) cells were observed in the periphery. Upon specific activation, transgenic CD8(+) T cells efficiently lysed syngeneic MHC class I(+) targets in vitro and in vivo, indicating that thymic cortical epithelium (in contrast to medullary epithelium and antigen-presenting cells of hematopoietic origin) is incapable of tolerance induction. Thus, compartmentalization of the antigen-presenting cells involved in thymic positive selection and tolerance induction can (at least in part) explain the positive/negative selection paradox.