In pursuit of an accurate spatial and temporal model of biomolecules at the atomistic level: a perspective on computer simulation

Atomistic Molecular Dynamics Simulations of DNA in Complex 3D Arrangements for Comparison with Lower Resolution Structural Experiments

Atomic-level computer simulations are a very useful tool for describing the structure and dynamics of complex biomolecules such as DNA and for providing detail at a resolution where experimental techniques cannot arrive. Molecular dynamics (MD) simulations of mechanically distorted DNA caused by agents like supercoiling and protein binding are computationally challenging due to the large size of the associated systems and timescales. However, nowadays they are achievable thanks to the efficient usage of GPU and to the improvements of continuum solvation models. This together with the concurrent improvements in the resolution of single-molecule experiments, such as atomic force microscopy (AFM), makes possible the convergence between the two. Here we present detailed protocols for doing so: for performing molecular dynamics (MD) simulations of DNA adopting complex three-dimensional arrangements and for comparing the outcome of the calculations with single-molecule experimental data with a lower resolution than atomic.

Exploring the dynamics of flagellar dynein within the axoneme with Fluctuating Finite Element Analysis

Flagellar dyneins are the molecular motors responsible for producing the propagating bending motions of cilia and flagella. They are located within a densely packed and highly organised super-macromolecular cytoskeletal structure known as the axoneme. Using the mesoscale simulation technique Fluctuating Finite Element Analysis (FFEA), which represents proteins as viscoelastic continuum objects subject to explicit thermal noise, we have quantified the constraints on the range of molecular conformations that can be explored by dynein-c within the crowded architecture of the axoneme. We subsequently assess the influence of crowding on the 3D exploration of microtubule-binding sites, and specifically on the axial step length. Our calculations combine experimental information on the shape, flexibility and environment of dynein-c from three distinct sources; negative stain electron microscopy, cryo-electron microscopy (cryo-EM) and cryo-electron tomography (cryo-ET). Our FFEA simulations show that the super-macromolecular organisation of multiple protein complexes into higher-order structures can have a significant influence on the effective flexibility of the individual molecular components, and may, therefore, play an important role in the physical mechanisms underlying their biological function.

Securing the future of research computing in the biosciences

Improvements in technology often drive scientific discovery. Therefore, research requires sustained investment in the latest equipment and training for the researchers who are going to use it. Prioritising and administering infrastructure investment is challenging because future needs are difficult to predict. In the past, highly computationally demanding research was associated primarily with particle physics and astronomy experiments. However, as biology becomes more quantitative and bioscientists generate more and more data, their computational requirements may ultimately exceed those of physical scientists. Computation has always been central to bioinformatics, but now imaging experiments have rapidly growing data processing and storage requirements. There is also an urgent need for new modelling and simulation tools to provide insight and understanding of these biophysical experiments. Bioscience communities must work together to provide the software and skills training needed in their areas. Research-active institutions need to recognise that computation is now vital in many more areas of discovery and create an environment where it can be embraced. The public must also become aware of both the power and limitations of computing, particularly with respect to their health and personal data.

Focused conformational sampling in proteins

A detailed understanding of the conformational dynamics of biological molecules is difficult to obtain by experimental techniques due to resolution limitations in both time and space. Computer simulations avoid these in theory but are often too short to sample rare events reliably. Here we show that the progress index-guided sampling (PIGS) protocol can be used to enhance the sampling of rare events in selected parts of biomolecules without perturbing the remainder of the system. The method is very easy to use as it only requires as essential input a set of several features representing the parts of interest sufficiently. In this feature space, new states are discovered by spontaneous fluctuations alone and in unsupervised fashion. Because there are no energetic biases acting on phase space variables or projections thereof, the trajectories PIGS generates can be analyzed directly in the framework of transition networks. We demonstrate the possibility and usefulness of such focused explorations of biomolecules with two loops that are part of the binding sites of bromodomains, a family of epigenetic "reader" modules. This real-life application uncovers states that are structurally and kinetically far away from the initial crystallographic structures and are also metastable. Representative conformations are intended to be used in future high-throughput virtual screening campaigns.

Multiscale Methods in Drug Design Bridge Chemical and Biological Complexity in the Search for Cures

Drug action is inherently multiscale: it connects molecular interactions to emergent properties at cellular and larger scales. Simulation techniques at each of these different scales are already central to drug design and development, but methods capable of connecting across these scales will extend understanding of complex mechanisms and the ability to predict biological effects. Improved algorithms, ever-more-powerful computing architectures and the accelerating growth of rich datasets are driving advances in multiscale modeling methods capable of bridging chemical and biological complexity from the atom to the cell. Particularly exciting is the development of highly detailed, structure-based, physical simulations of biochemical systems, which are now able to access experimentally relevant timescales for large systems and, at the same time, achieve unprecedented accuracy. In this Perspective, we discuss how emerging data-rich, physics-based multiscale approaches are of the cusp of realizing long-promised impact in the discovery, design and development of novel therapeutics. We highlight emerging methods and applications in this growing field, and outline how different scales can be combined in practical modelling and simulation strategies.

Fluctuating Finite Element Analysis (FFEA): A continuum mechanics software tool for mesoscale simulation of biomolecules

Fluctuating Finite Element Analysis (FFEA) is a software package designed to perform continuum mechanics simulations of proteins and other globular macromolecules. It combines conventional finite element methods with stochastic thermal noise, and is appropriate for simulations of large proteins and protein complexes at the mesoscale (length-scales in the range of 5 nm to 1 μm), where there is currently a paucity of modelling tools. It requires 3D volumetric information as input, which can be low resolution structural information such as cryo-electron tomography (cryo-ET) maps or much higher resolution atomistic co-ordinates from which volumetric information can be extracted. In this article we introduce our open source software package for performing FFEA simulations which we have released under a GPLv3 license. The software package includes a C ++ implementation of FFEA, together with tools to assist the user to set up the system from Electron Microscopy Data Bank (EMDB) or Protein Data Bank (PDB) data files. We also provide a PyMOL plugin to perform basic visualisation and additional Python tools for the analysis of FFEA simulation trajectories. This manuscript provides a basic background to the FFEA method, describing the implementation of the core mechanical model and how intermolecular interactions and the solvent environment are included within this framework. We provide prospective FFEA users with a practical overview of how to set up an FFEA simulation with reference to our publicly available online tutorials and manuals that accompany this first release of the package.

Atomistic Molecular Dynamics Simulations of DNA Minicircle Topoisomers: A Practical Guide to Setup, Performance, and Analysis

While DNA supercoiling is ubiquitous in vivo, the structure of supercoiled DNA is more challenging to study experimentally than simple linear sequences because the DNA must have a closed topology in order to sustain superhelical stress. DNA minicircles, which are closed circular double-stranded DNA sequences typically containing between 60 and 500 base pairs, have proven to be useful biochemical tools for the study of supercoiled DNA mechanics. We present detailed protocols for constructing models of DNA minicircles in silico, for performing atomistic molecular dynamics (MD) simulations of supercoiled minicircle DNA, and for analyzing the results of the calculations. These simulations are computationally challenging due to the large system sizes. However, improvements in parallel computing software and hardware promise access to improve conformational sampling and simulation timescales. Given the concurrent improvements in the resolution of experimental techniques such as atomic force microscopy (AFM) and cryo-electron microscopy, the study of DNA minicircles will provide a more complete understanding of both the structure and the mechanics of supercoiled DNA.

Encapsulated membrane proteins: A simplified system for molecular simulation

Over the past 50years there has been considerable progress in our understanding of biomolecular interactions at an atomic level. This in turn has allowed molecular simulation methods employing full atomistic modelling at ever larger scales to develop. However, some challenging areas still remain where there is either a lack of atomic resolution structures or where the simulation system is inherently complex. An area where both challenges are present is that of membranes containing membrane proteins. In this review we analyse a new practical approach to membrane protein study that offers a potential new route to high resolution structures and the possibility to simplify simulations. These new approaches collectively recognise that preservation of the interaction between the membrane protein and the lipid bilayer is often essential to maintain structure and function. The new methods preserve these interactions by producing nano-scale disc shaped particles that include bilayer and the chosen protein. Currently two approaches lead in this area: the MSP system that relies on peptides to stabilise the discs, and SMALPs where an amphipathic styrene maleic acid copolymer is used. Both methods greatly enable protein production and hence have the potential to accelerate atomic resolution structure determination as well as providing a simplified format for simulations of membrane protein dynamics. This article is part of a Special Issue entitled: Biosimulations edited by Ilpo Vattulainen and Tomasz Róg.