Extensive changes in DNA methylation are common in cancer and may contribute to oncogenesis through transcriptional silencing of tumor-suppressor genes1. Genome-scale studies have yielded important insights into these changes2, 3, 4, 5 but have focused on CpG islands or gene promoters. We used whole-genome bisulfite sequencing (bisulfite-seq) to comprehensively profile a primary human colorectal tumor and adjacent normal colon tissue at single-basepair resolution. Regions of focal hypermethylation in the tumor were located primarily at CpG islands and were concentrated within regions of long-range (>100 kb) hypomethylation. These hypomethylated domains covered nearly half of the genome and coincided with late replication and attachment to the nuclear lamina in human cell lines. We confirmed the confluence of hypermethylation and hypomethylation within these domains in 25 diverse colorectal tumors and matched adjacent tissue. We propose that widespread DNA methylation changes in cancer are linked to silencing programs orchestrated by the three-dimensional organization of chromatin within the nucleus.

DNA methylation loss occurs frequently in cancer genomes, primarily within lamina-associated, late-replicating regions termed Partially Methylated Domains (PMDs). We profiled 39 diverse primary tumors and 8 matched adjacent tissues using Whole-Genome Bisulfite Sequencing (WGBS), and analyzed them alongside 343 additional human and 206 mouse WGBS datasets. We identified a local CpG sequence context associated with preferential hypomethylation in PMDs. Analysis of CpGs in this context (“Solo-WCGWs”) revealed previously undetected PMD hypomethylation in almost all healthy tissue types. PMD hypomethylation increased with age, beginning during fetal development, and appeared to track the accumulation of cell divisions. In cancer, PMD hypomethylation depth correlated with somatic mutation density and cell-cycle gene expression, consistent with its reflection of mitotic history, and suggesting its application as a mitotic clock. We propose that late replication leads to lifelong progressive methylation loss, which acts as a biomarker for cellular aging and which may contribute to oncogenesis.

Recent discoveries suggest cysteine-stabilized toxins and antimicrobial peptides have structure-activity parallels derived by common ancestry. Here, human antimicrobial peptide hBD-2 and rattlesnake venom-toxin crotamine were compared in phylogeny, 3D structure, target cell specificity, and mechanisms of action. Results indicate a striking degree of structural and phylogenetic congruence. Importantly, these polypeptides also exhibited functional reciprocity: (i) they exerted highly similar antimicrobial pH optima and spectra; (ii) both altered membrane potential consistent with ion channel-perturbing activities; and (iii) both peptides induced phosphatidylserine accessibility in eukaryotic cells. However, the Nav channel-inhibitor tetrodotoxin antagonized hBD-2 mechanisms, but not those of crotamine. As crotamine targets eukaryotic ion channels, computational docking was used to compare hBD-2 versus crotamine interactions with prototypic bacterial, fungal, or mammalian Kv channels. Models support direct interactions of each peptide with Kv channels. However, while crotamine localized to occlude Kv channels in eukaryotic but not prokaryotic cells, hBD-2 interacted with prokaryotic and eukaryotic Kv channels but did not occlude either. Together, these results support the hypothesis that antimicrobial and cytotoxic polypeptides have ancestral structurefunction homology, but evolved to preferentially target respective microbial versus mammalian ion channels via residue-specific interactions. These insights may accelerate development of antiinfective or therapeutic peptides that selectively target microbial or abnormal host cells.channel ͉ defensin ͉ host defense ͉ toxin

ABSTRACT:Topologically, platelet factor-4 kinocidins consist of distinct N-terminal extended, C-terminal helical, and interposing g-core structural domains. The C-terminal a-helices autonomously confer direct microbicidal activity, and the synthetic antimicrobial peptide RP-1 is modeled upon these domains. In this study, the structure of RP-1 was assessed using several complementary techniques. The highresolution structure of RP-1 was determined by NMR in anionic sodium dodecyl sulfate (SDS) and zwitterionic dodecylphosphocholine (DPC) micelles, which approximate prokaryotic and eukaryotic membranes, respectively. NMR data indicate the peptide assumes an amphipathic a-helical backbone conformation in both micelle environments. However, small differences were observed in the side-chain orientations of lysine, tyrosine, and phenylalanine residues in SDS versus DPC environments.NMR experiments with a paramagnetic probe indicated differences in positioning of the peptide within the two micelle types. Molecular dynamics (MD) simulations of the peptide in both micelle types were also performed to add insight into the peptide/micelle interactions and to assess the validity of this technique to predict the structure of peptides in complex with micelles. MD independently predicted RP-1 to interact only peripherally with the DPC micelle, leaving its spherical shape intact. In contrast, RP-1 entered deeply into and significantly distorted the SDS micelle. Overall, the experimental and MD results support a preferential specificity of RP-1 for anionic membranes over zwitterionic membranes. This specificity likely derives from differences in RP-1 interaction with distinct lipid systems, including subtle differences in side chain orientations, rather than gross changes in RP-1 structure in the two lipid environments. #

In this communication we introduce an efficient implementation of adaptive biasing that greatly improves the speed of free energy computation in molecular dynamics simulations. We investigated the use of accelerated simulations to inform on compound design using a recently reported and clinically relevant inhibitor of the chromatin regulator BRD4 (bromodomain-containing protein 4). Benchmarking on our local compute cluster, our implementation achieves up to 2.5 times more force calls per day than plumed2. Results of five 1 µs-long simulations are presented, which reveal a conformational switch in the BRD4 inhibitor between a binding competent and incompetent state. Stabilization of the switch led to a −3 kcal/mol improvement of absolute binding free energy. These studies suggest an unexplored ligand design principle and offer new actionable hypotheses for medicinal chemistry efforts against this druggable epigenetic target class. Published by AIP Publishing. [http://dx