Free energy computation by controlled Langevin processes

We propose a nonequilibrium sampling method for computing free energy profiles along a given reaction coordinate. The method consists of two parts: a controlled Langevin sampler that generates nonequilibrium bridge paths conditioned by the reaction coordinate, and Jarzynski’s formula for reweighting the paths. Our derivation of the equations of motion of the sampler is based on stochastic perturbation of a controlled dissipative Hamiltonian system, for which we prove Jarzynski’s identity as a special case of the Feynman-Kac formula. We illustrate our method by means of a suitable numerical example and briefly discuss issues of optimally choosing the control protocol for the reaction coordinate

Sequential Change Point Detection in Molecular Dynamics Trajectories

Motivated from a molecular dynamics context we propose a sequential change point detection algorithm for vector-valued autoregressive models based on Bayesian model selection. The algorithm does not rely on any sampling procedure or assumptions underlying the dynamics of the transitions and is designed to cope with high-dimensional data. We show the applicability of the algorithm on a time series obtained from numerical simulation of a penta-peptide molecule

A Structure-preserving numerical discretization of reversible diffusions

We propose a robust and efficient numerical discretization scheme for the infinitesimal generator of a diffusion process based on a finite volume approximation. The resulting discrete-space operator can be interpreted as a jump process on the mesh whose invariant measure is precisely the cell approximation of the Boltzmann distribution of the original process. Moreover the resulting jump process preserves the detailed balance property of the original stochastic process

Genomic and SNP Analyses Demonstrate a Distant Separation of the Hospital and Community-Associated Clades of Enterococcus faecium

Recent studies have pointed to the existence of two subpopulations of Enterococcus faecium, one containing primarily commensal/community-associated (CA) strains and one that contains most clinical or hospital-associated (HA) strains, including those classified by multi-locus sequence typing (MLST) as belonging to the CC17 group. The HA subpopulation more frequently has IS16, pathogenicity island(s), and plasmids or genes associated with antibiotic resistance, colonization, and/or virulence. Supporting the two clades concept, we previously found a 3–10% difference between four genes from HA-clade strains vs. CA-clade strains, including 5% difference between pbp5-R of ampicillin-resistant, HA strains and pbp5-S of ampicillin-sensitive, CA strains. To further investigate the core genome of these subpopulations, we studied 100 genes from 21 E. faecium genome sequences; our analyses of concatenated sequences, SNPs, and individual genes all identified two distinct groups. With the concatenated sequence, HA-clade strains differed by 0–1% from one another while CA clade strains differed from each other by 0–1.1%, with 3.5–4.2% difference between the two clades. While many strains had a few genes that grouped in one clade with most of their genes in the other clade, one strain had 28% of its genes in the CA clade and 72% in the HA clade, consistent with the predicted role of recombination in the evolution of E. faecium. Using estimates for Escherichia coli, molecular clock calculations using sSNP analysis indicate that these two clades may have diverged ≥1 million years ago or, using the higher mutation rate for Bacillus anthracis, ∼300,000 years ago. These data confirm the existence of two clades of E. faecium and show that the differences between the HA and CA clades occur at the core genomic level and long preceded the modern antibiotic era

Matrix Metalloproteinase-2 and -9 Secreted by Leukemic Cells Increase the Permeability of Blood-Brain Barrier by Disrupting Tight Junction Proteins

Central nervous system (CNS) involvement remains an important cause of morbidity and mortality in acute leukemia, the mechanisms of leukemic cell infiltration into the CNS have not yet been elucidated. The blood-brain barrier (BBB) makes CNS become a refugee to leukemic cells and serves as a resource of cells that seed extraneural sites. How can the leukemic cells disrupt this barrier and invasive the CNS, even if many of the currently available chemotherapies can not cross the BBB? Tight junction in endothelial cells occupies a central role in the function of the BBB. Except the well known role of degrading extracellular matrix in metastasis of cancer cells, here we show matrix metalloproteinase (MMP)-2 and -9, secreted by leukemic cells, mediate the BBB opening by disrupting tight junction proteins in the CNS leukemia. We demonstrated that leukemic cells impaired tight junction proteins ZO-1, claudin-5 and occludin resulting in increased permeability of the BBB. However, these alterations reduced when MMP-2 and -9 activities were inhibited by RNA interference strategy or by MMP inhibitor GM6001 in an in vitro BBB model. We also found that the disruption of the BBB in company with the down-regulation of ZO-1, claudin-5 and occludin and the up-regulation of MMP-2 and -9 in mouse brain tissues with leukemic cell infiltration by confocal imaging and the assay of in situ gelatin zymography. Besides, GM6001 protected all mice against CNS leukemia. Our findings suggest that the degradation of tight junction proteins ZO-1, claudin-5 and occludin by MMP-2 and -9 secreted by leukemic cells constitutes an important mechanism in the BBB breakdown which contributes to the invasion of leukemic cells to the CNS in acute leukemia

Genome Sequence of a Lancefield Group C Streptococcus zooepidemicus Strain Causing Epidemic Nephritis: New Information about an Old Disease

Outbreaks of disease attributable to human error or natural causes can provide unique opportunities to gain new information about host-pathogen interactions and new leads for pathogenesis research. Poststreptococcal glomerulonephritis (PSGN), a sequela of infection with pathogenic streptococci, is a common cause of preventable kidney disease worldwide. Although PSGN usually occurs after infection with group A streptococci, organisms of Lancefield group C and G also can be responsible. Despite decades of study, the molecular pathogenesis of PSGN is poorly understood. As a first step toward gaining new information about PSGN pathogenesis, we sequenced the genome of Streptococcus equi subsp. zooepidemicus strain MGCS10565, a group C organism that caused a very large and unusually severe epidemic of nephritis in Brazil. The genome is a circular chromosome of 2,024,171 bp. The genome shares extensive gene content, including many virulence factors, with genetically related group A streptococci, but unexpectedly lacks prophages. The genome contains many apparently foreign genes interspersed around the chromosome, consistent with the presence of a full array of genes required for natural competence. An inordinately large family of genes encodes secreted extracellular collagen-like proteins with multiple integrin-binding motifs. The absence of a gene related to speB rules out the long-held belief that streptococcal pyrogenic exotoxin B or antibodies reacting with it singularly cause PSGN. Many proteins previously implicated in GAS PSGN, such as streptokinase, are either highly divergent in strain MGCS10565 or are not more closely related between these species than to orthologs present in other streptococci that do not commonly cause PSGN. Our analysis provides a comparative genomics framework for renewed appraisal of molecular events underlying APSGN pathogenesis

Cellular binding partners of the human papillomavirus E6 protein

The high-risk strains of human papillomavirus (HR-HPV) are known to be causative agents of cervical cancer and have recently also been implicated in cancers of the oropharynx. E6 is a potent oncogene of HR-HPVs, and its role in the progression to malignancy has been and continues to be explored. E6 is known to interact with and subsequently inactivate numerous cellular proteins pivotal in the mediation of apoptosis, transcription of tumor suppressor genes, maintenance of epithelial organization, and control of cell proliferation. Binding of E6 to these proteins cumulatively contributes to the oncogenic potential of HPV. This paper provides an overview of these cellular protein partners of HR-E6, the motifs known to mediate oncoprotein binding, and the agents that have the potential to interfere with E6 expression and activity and thus prevent the subsequent progression to oncogenesis