Modelling biomacromolecular assemblies with continuum mechanics

Computational models for studying physical instabilities in high concentration biotherapeutic formulations

Computational prediction of the behavior of concentrated protein solutions is particularly advantageous in early development stages of biotherapeutics when material availability is limited and a large set of formulation conditions needs to be explored. This review provides an overview of the different computational paradigms that have been successfully used in modeling undesirable physical behaviors of protein solutions with a particular emphasis on high-concentration drug formulations. This includes models ranging from all-atom simulations, coarse-grained representations to macro-scale mathematical descriptions used to study physical instability phenomena of protein solutions such as aggregation, elevated viscosity, and phase separation. These models are compared and summarized in the context of the physical processes and their underlying assumptions and limitations. A detailed analysis is also given for identifying protein interaction processes that are explicitly or implicitly considered in the different modeling approaches and particularly their relations to various formulation parameters. Lastly, many of the shortcomings of existing computational models are discussed, providing perspectives and possible directions toward an efficient computational framework for designing effective protein formulations.

From System Modeling to System Analysis: The Impact of Resolution Level and Resolution Distribution in the Computer-Aided Investigation of Biomolecules

The ever increasing computer power, together with the improved accuracy of atomistic force fields, enables researchers to investigate biological systems at the molecular level with remarkable detail. However, the relevant length and time scales of many processes of interest are still hardly within reach even for state-of-the-art hardware, thus leaving important questions often unanswered. The computer-aided investigation of many biological physics problems thus largely benefits from the usage of coarse-grained models, that is, simplified representations of a molecule at a level of resolution that is lower than atomistic. A plethora of coarse-grained models have been developed, which differ most notably in their granularity; this latter aspect determines one of the crucial open issues in the field, i.e. the identification of an optimal degree of coarsening, which enables the greatest simplification at the expenses of the smallest information loss. In this review, we present the problem of coarse-grained modeling in biophysics from the viewpoint of system representation and information content. In particular, we discuss two distinct yet complementary aspects of protein modeling: on the one hand, the relationship between the resolution of a model and its capacity of accurately reproducing the properties of interest; on the other hand, the possibility of employing a lower resolution description of a detailed model to extract simple, useful, and intelligible information from the latter.

Combined Force-Torque Spectroscopy of Proteins by Means of Multiscale Molecular Simulation

Assessing the structural properties of large proteins is important to gain an understanding of their function in, e.g., biological systems or biomedical applications. We propose a method to examine the mechanical properties of proteins subject to applied forces by means of multiscale simulation. Both stretching and torsional forces are considered, and these may be applied independently of each other. As a proof of principle, we apply torsional forces to a coarse-grained continuum model of the antibody protein immunoglobulin G using fluctuating finite element analysis and use it to identify the area of strongest deformation. This region is essential to the torsional properties of the molecule as a whole because it represents the softest, most deformable domain. Zooming in, this part of the molecule is subjected to torques and stretching forces using molecular dynamics simulations on an atomistically resolved level to investigate its torsional properties. We calculate the torsional resistance as a function of the rotation of the domain while subjecting it to various stretching forces. From this, we assess how the measured twist-torque profiles develop with increasing stretching force and show that they exhibit torsion stiffening, in qualitative agreement with experimental findings. We argue that combining the twist-torque profiles for various stretching forces effectively results in a combined force-torque spectroscopy analysis, which may serve as a mechanical signature for a biological macromolecule.

Nano-encapsulated Escherichia coli Divisome Anchor ZipA, and in Complex with FtsZ

The E. coli membrane protein ZipA, binds to the tubulin homologue FtsZ, in the early stage of cell division. We isolated ZipA in a Styrene Maleic Acid lipid particle (SMALP) preserving its position and integrity with native E. coli membrane lipids. Direct binding of ZipA to FtsZ is demonstrated, including FtsZ fibre bundles decorated with ZipA. Using Cryo-Electron Microscopy, small-angle X-ray and neutron scattering, we determine the encapsulated-ZipA structure in isolation, and in complex with FtsZ to a resolution of 1.6 nm. Three regions can be identified from the structure which correspond to, SMALP encapsulated membrane and ZipA transmembrane helix, a separate short compact tether, and ZipA globular head which binds FtsZ. The complex extends 12 nm from the membrane in a compact structure, supported by mesoscale modelling techniques, measuring the movement and stiffness of the regions within ZipA provides molecular scale analysis and visualisation of the early divisome.

Survey of the analysis of continuous conformational variability of biological macromolecules by electron microscopy

Single-particle analysis by electron microscopy is a well established technique for analyzing the three-dimensional structures of biological macromolecules. Besides its ability to produce high-resolution structures, it also provides insights into the dynamic behavior of the structures by elucidating their conformational variability. Here, the different image-processing methods currently available to study continuous conformational changes are reviewed.

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.