Local Simulation Algorithms for Coulombic Interactions

We consider dynamically constrained Monte-Carlo dynamics and show that this leads to the generation of long ranged effective interactions. This allows us to construct a local algorithm for the simulation of charged systems without ever having to evaluate pair potentials or solve the Poisson equation. We discuss a simple implementation of a charged lattice gas as well as more elaborate off-lattice versions of the algorithm. There are analogies between our formulation of electrostatics and the bosonic Hubbard model in the phase approximation. Cluster methods developed for this model further improve the efficiency of the electrostatics algorithm.Comment: Proceedings Statphys22 10 page

Atomic Force Stiffness Imaging: Capturing differences in mechanical properties to identify and localize areas of prostate cancer tissue

Prostate cancer is now the most commonly diagnosed cancer in men in western countries. Due to the difficulty for early detection, there are an estimated 10000 deaths a year in the UK from prostate cancer alone; whereby the only curative option is interventional treatment that aims to excise all diseased cells while preserving the neurovascular bundle. To date, several studies have shown that the mechanical properties of cancer cells and tissues i.e. adhesion, stiffness, roughness and viscoelasticity are significantly different from benign cells and regions of tissue that are healthy. Building upon these results, we believe novel methods of imaging the mechanical properties of prostate cancer samples can provide new surgical intervention opportunities beyond what is possible through vision alone. In this paper, we used an Atomic Force Microscope (AFM) to measure the stiffness and topography variations correlating to regions of prostate cancer at the surface of an excised sample at a cellular level. Preliminary results show that by using an AFM we can detect structural differences in non-homogeneous tissue samples, confirming previous results that cancerous tissues appear stiffer than benign areas. Through these results, we aim to develop a stiffness imaging protocol to aid the early detection of prostate cancer, in addition to force sensing surgical tools

Targeted therapy of the XIAP/proteasome pathway overcomes TRAIL-resistance in carcinoma by switching apoptosis signaling to a Bax/Bak-independent 'type I' mode

TRAIL is a promising anticancer agent, capable of inducing apoptosis in a wide range of treatment-resistant tumor cells. In 'type II' cells, the death signal triggered by TRAIL requires amplification via the mitochondrial apoptosis pathway. Consequently, deregulation of the intrinsic apoptosis-signaling pathway, for example, by loss of Bax and Bak, confers TRAIL-resistance and limits its application. Here, we show that despite resistance of Bax/Bak double-deficient cells, TRAIL-treatment resulted in caspase-8 activation and complete processing of the caspase-3 proenzymes. However, active caspase-3 was degraded by the proteasome and not detectable unless the XIAP/proteasome pathway was inhibited. Direct or indirect inhibition of XIAP by RNAi, Mithramycin A or by the SMAC mimetic LBW-242 as well as inhibition of the proteasome by Bortezomib overcomes TRAIL-resistance of Bax/Bak double-deficient tumor cells. Moreover, activation and stabilization of caspase-3 becomes independent of mitochondrial death signaling, demonstrating that inhibition of the XIAP/proteasome pathway overcomes resistance by converting 'type II' to 'type I' cells. Our results further demonstrate that the E3 ubiquitin ligase XIAP is a gatekeeper critical for the 'type II' phenotype. Pharmacological manipulation of XIAP therefore is a promising strategy to sensitize cells for TRAIL and to overcome TRAIL-resistance in case of central defects in the intrinsic apoptosis-signaling pathway

How does conformational flexibility influence key structural features involved in activation of anaplastic lymphoma kinase?

Anaplastic Lymphoma Kinase (ALK) plays a major role in developing tumor processes and therefore has emerged as a validated therapeutic target. Applying atomistic molecular dynamics simulations on the wild type enzyme and the nine most frequently occurring and clinically important activation mutants we revealed important conformational effects on key interactions responsible for the activation of the enzyme

Tackling Exascale Software Challenges in Molecular Dynamics Simulations with GROMACS

GROMACS is a widely used package for biomolecular simulation, and over the last two decades it has evolved from small-scale efficiency to advanced heterogeneous acceleration and multi-level parallelism targeting some of the largest supercomputers in the world. Here, we describe some of the ways we have been able to realize this through the use of parallelization on all levels, combined with a constant focus on absolute performance. Release 4.6 of GROMACS uses SIMD acceleration on a wide range of architectures, GPU offloading acceleration, and both OpenMP and MPI parallelism within and between nodes, respectively. The recent work on acceleration made it necessary to revisit the fundamental algorithms of molecular simulation, including the concept of neighborsearching, and we discuss the present and future challenges we see for exascale simulation - in particular a very fine-grained task parallelism. We also discuss the software management, code peer review and continuous integration testing required for a project of this complexity.Comment: EASC 2014 conference proceedin

Local Simulation Algorithms for Coulomb Interaction

Long ranged electrostatic interactions are time consuming to calculate in molecular dynamics and Monte-Carlo simulations. We introduce an algorithmic framework for simulating charged particles which modifies the dynamics so as to allow equilibration using a local Hamiltonian. The method introduces an auxiliary field with constrained dynamics so that the equilibrium distribution is determined by the Coulomb interaction. We demonstrate the efficiency of the method by simulating a simple, charged lattice gas.Comment: Last figure changed to improve demonstration of numerical efficienc

Conformational effects on the Circular Dichroism of Human Carbonic Anhydrase II: a multilevel computational study

Circular Dichroism (CD) spectroscopy is a powerful method for investigating conformational changes in proteins and therefore has numerous applications in structural and molecular biology. Here a computational investigation of the CD spectrum of the Human Carbonic Anhydrase II (HCAII), with main focus on the near-UV CD spectra of the wild-type enzyme and it seven tryptophan mutant forms, is presented and compared to experimental studies. Multilevel computational methods (Molecular Dynamics, Semiempirical Quantum Mechanics, Time-Dependent Density Functional Theory) were applied in order to gain insight into the mechanisms of interaction between the aromatic chromophores within the protein environment and understand how the conformational flexibility of the protein influences these mechanisms. The analysis suggests that combining CD semi empirical calculations, crystal structures and molecular dynamics (MD) could help in achieving a better agreement between the computed and experimental protein spectra and provide some unique insight into the dynamic nature of the mechanisms of chromophore interactions

Fourier Acceleration of Langevin Molecular Dynamics

Fourier acceleration has been successfully applied to the simulation of lattice field theories for more than a decade. In this paper, we extend the method to the dynamics of discrete particles moving in continuum. Although our method is based on a mapping of the particles' dynamics to a regular grid so that discrete Fourier transforms may be taken, it should be emphasized that the introduction of the grid is a purely algorithmic device and that no smoothing, coarse-graining or mean-field approximations are made. The method thus can be applied to the equations of motion of molecular dynamics (MD), or its Langevin or Brownian variants. For example, in Langevin MD simulations our acceleration technique permits a straightforward spectral decomposition of forces so that the long-wavelength modes are integrated with a longer time step, thereby reducing the time required to reach equilibrium or to decorrelate the system in equilibrium. Speedup factors of up to 30 are observed relative to pure (unaccelerated) Langevin MD. As with acceleration of critical lattice models, even further gains relative to the unaccelerated method are expected for larger systems. Preliminary results for Fourier-accelerated molecular dynamics are presented in order to illustrate the basic concepts. Possible extensions of the method and further lines of research are discussed.Comment: 11 pages, two illustrations included using graphic

Stain-free identification of tissue pathology using a generative adversarial network to infer nanomechanical signatures

Intraoperative frozen section analysis can be used to improve the accuracy of tumour margin estimation during cancer resection surgery through rapid processing and pathological assessment of excised tissue. Its applicability is limited in some cases due to the additional risks associated with prolonged surgery, largely from the time-consuming staining procedure. Our work uses a measurable property of bulk tissue to bypass the staining process: as tumour cells proliferate, they influence the surrounding extra-cellular matrix, and the resulting change in elastic modulus provides a signature of the underlying pathology. In this work we accurately localise atomic force microscopy measurements of human liver tissue samples and train a generative adversarial network to infer elastic modulus from low-resolution images of unstained tissue sections. Pathology is predicted through unsupervised clustering of parameters characterizing the distributions of inferred values, achieving 89% accuracy for all samples based on the nominal assessment (n = 28), and 95% for samples that have been validated by two independent pathologists through post hoc staining (n = 20). Our results demonstrate that this technique could increase the feasibility of intraoperative frozen section analysis for use during resection surgery and improve patient outcomes