Atomistic and coarse-grained simulations of membrane proteins: A practical guide

Lipids and proteins: Insights into the dynamics of assembly, recognition, condensate formation. What is still missing?

Lipid membranes and proteins, which are part of us throughout our lives, have been studied for decades. However, every year, new discoveries show how little we know about them. In a reader-friendly manner for people not involved in the field, this paper tries to serve as a bridge between physicists and biologists and new young researchers diving into the field to show its relevance, pointing out just some of the plethora of lines of research yet to be unraveled. It illustrates how new ways, from experimental to theoretical approaches, are needed in order to understand the structures and interactions that take place in a single lipid, protein, or multicomponent system, as we are still only scratching the surface.

Transmembrane proteins-Different anchoring systems

Transmembrane proteins are active in amphipathic environments. To stabilize the protein in such surrounding the exposure of hydrophobic residues on the protein surface is required. Transmembrane proteins are responsible for the transport of various molecules. Therefore, they often represent structures in the form of channels. This analysis focused on the stability and local flexibility of transmembrane proteins, particularly those related to their biological activity. Different forms of anchorage were identified using the fuzzy oil-drop model (FOD) and its modified form, FOD-M. The mainly helical as well as β-barrel structural forms are compared with respect to the mechanism of stabilization in the cell membrane. The different anchoring system was found to stabilize protein molecules with possible local fluctuation.

Building complex membranes with Martini 3

The Martini model is a popular force field for coarse-grained simulations. Membranes have always been at the center of its development, with the latest version, Martini 3, showing great promise in capturing more and more realistic behavior. In this chapter we provide a step-by-step tutorial on how to construct starting configurations, run initial simulations and perform dedicated analysis for membrane-based systems of increasing complexity, including leaflet asymmetry, curvature gradients and embedding of membrane proteins.

Challenges in modelling dynamic processes in realistic nanostructured materials at operating conditions

The question is addressed in how far current modelling strategies are capable of modelling dynamic phenomena in realistic nanostructured materials at operating conditions. Nanostructured materials used in applications are far from perfect; they possess a broad range of heterogeneities in space and time extending over several orders of magnitude. Spatial heterogeneities from the subnanometre to the micrometre scale in crystal particles with a finite size and specific morphology, impact the material's dynamics. Furthermore, the material's functional behaviour is largely determined by the operating conditions. Currently, there exists a huge length-time scale gap between attainable theoretical length-time scales and experimentally relevant scales. Within this perspective, three key challenges are highlighted within the molecular modelling chain to bridge this length-time scale gap. Methods are needed that enable (i) building structural models for realistic crystal particles having mesoscale dimensions with isolated defects, correlated nanoregions, mesoporosity, internal and external surfaces; (ii) the evaluation of interatomic forces with quantum mechanical accuracy albeit at much lower computational cost than the currently used density functional theory methods and (iii) derivation of the kinetics of phenomena taking place in a multi-length-time scale window to obtain an overall view of the dynamics of the process. This article is part of a discussion meeting issue 'Supercomputing simulations of advanced materials'.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Multifaceted Computational Modeling in Glycoscience

Glycoscience assembles all the scientific disciplines involved in studying various molecules and macromolecules containing carbohydrates and complex glycans. Such an ensemble involves one of the most extensive sets of molecules in quantity and occurrence since they occur in all microorganisms and higher organisms. Once the compositions and sequences of these molecules are established, the determination of their three-dimensional structural and dynamical features is a step toward understanding the molecular basis underlying their properties and functions. The range of the relevant computational methods capable of addressing such issues is anchored by the specificity of stereoelectronic effects from quantum chemistry to mesoscale modeling throughout molecular dynamics and mechanics and coarse-grained and docking calculations. The Review leads the reader through the detailed presentations of the applications of computational modeling. The illustrations cover carbohydrate-carbohydrate interactions, glycolipids, and N- and O-linked glycans, emphasizing their role in SARS-CoV-2. The presentation continues with the structure of polysaccharides in solution and solid-state and lipopolysaccharides in membranes. The full range of protein-carbohydrate interactions is presented, as exemplified by carbohydrate-active enzymes, transporters, lectins, antibodies, and glycosaminoglycan binding proteins. A final section features a list of 150 tools and databases to help address the many issues of structural glycobioinformatics.

Investigating the entropic nature of membrane-mediated interactions driving the aggregation of peripheral proteins

Peripheral membrane-associated proteins are known to accumulate on the surface of biomembranes as a result of membrane-mediated interactions. For a pair of rotationally-symmetric curvature-inducing proteins, membrane mechanics at the low-temperature limit predicts pure repulsion. On the other hand, temperature-dependent entropic forces arise between pairs of stiff-binding proteins suppressing membrane fluctuations. These Casimir-like interactions have thus been suggested as candidates for attractive forces leading to aggregation. With dense assemblies of peripheral proteins on the membrane, both these abstractions encounter short-range and multi-body complications. Here, we make use of a particle-based membrane model augmented with flexible peripheral proteins to quantify purely membrane-mediated interactions and investigate their underlying nature. We introduce a continuous reaction coordinate corresponding to the progression of protein aggregation. We obtain free energy and entropy landscapes for different surface concentrations along this reaction coordinate. In parallel, we investigate time-dependent estimates of membrane entropy corresponding to membrane undulations and coarse-grained director field and how they change dynamically with protein aggregation. Congruent outcomes of the two approaches point to the conclusion that for low surface concentrations, interactions with an entropic nature may drive the aggregation. But at high concentrations, enthalpic contributions due to concerted membrane deformation by protein clusters are dominant.

Specific interactions of peripheral membrane proteins with lipids: what can molecular simulations show us?

Peripheral membrane proteins (PMPs) can reversibly and specifically bind to biological membranes to carry out functions such as cell signalling, enzymatic activity, or membrane remodelling. Structures of these proteins and of their lipid-binding domains are typically solved in a soluble form, sometimes with a lipid or lipid headgroup at the binding site. To provide a detailed molecular view of PMP interactions with the membrane, computational methods such as molecular dynamics (MD) simulations can be applied. Here, we outline recent attempts to characterise these binding interactions, focusing on both intracellular proteins, such as phosphatidylinositol phosphate (PIP)-binding domains, and extracellular proteins such as glycolipid-binding bacterial exotoxins. We compare methods used to identify and analyse lipid-binding sites from simulation data and highlight recent work characterising the energetics of these interactions using free energy calculations. We describe how improvements in methodologies and computing power will help MD simulations to continue to contribute to this field in the future.

A Journey from Structure to Function of Bacterial Lipopolysaccharides

Lipopolysaccharide (LPS) is a crucial constituent of the outer membrane of most Gram-negative bacteria, playing a fundamental role in the protection of bacteria from environmental stress factors, in drug resistance, in pathogenesis, and in symbiosis. During the last decades, LPS has been thoroughly dissected, and massive information on this fascinating biomolecule is now available. In this Review, we will give the reader a third millennium update of the current knowledge of LPS with key information on the inherent peculiar carbohydrate chemistry due to often puzzling sugar residues that are uniquely found on it. Then, we will drive the reader through the complex and multifarious immunological outcomes that any given LPS can raise, which is strictly dependent on its chemical structure. Further, we will argue about issues that still remain unresolved and that would represent the immediate future of LPS research. It is critical to address these points to complete our notions on LPS chemistry, functions, and roles, in turn leading to innovative ways to manipulate the processes involving such a still controversial and intriguing biomolecule.

Computational studies of the mitochondrial carrier family SLC25. Present status and future perspectives

The members of the mitochondrial carrier family, also known as solute carrier family 25 (SLC25), are transmembrane proteins involved in the translocation of a plethora of small molecules between the mitochondrial intermembrane space and the matrix. These transporters are characterized by three homologous domains structure and a transport mechanism that involves the transition between different conformations. Mutations in regions critical for these transporters’ function often cause several diseases, given the crucial role of these proteins in the mitochondrial homeostasis. Experimental studies can be problematic in the case of membrane proteins, in particular concerning the characterization of the structure–function relationships. For this reason, computational methods are often applied in order to develop new hypotheses or to support/explain experimental evidence. Here the computational analyses carried out on the SLC25 members are reviewed, describing the main techniques used and the outcome in terms of improved knowledge of the transport mechanism. Potential future applications on this protein family of more recent and advanced in silico methods are also suggested.

Toward a Suite of Coarse-Grained Models for Molecular Simulation of Monoclonal Antibodies and Therapeutic Proteins

A series of coarse-grained models for molecular simulation of proteins are considered, with emphasis on the application of predicting protein-protein self-interactions for monoclonal antibodies (MAbs). As an illustrative example and for quantitative comparison, the models are used to predict osmotic virial coefficients over a broad range of attractive and repulsive self-interactions and solution conditions for a series of MAbs where the second osmotic virial coefficient has been experimentally determined in prior work. The models are compared based on how well they can predict experimental behavior, their computational burdens, and scalability. An intermediate-resolution model is also introduced that can capture specific electrostatic interactions with improved efficiency and similar or improved accuracy when compared to the previously published models. Guidance is included for the selection of coarse-grained models more generally for capturing a balance of electrostatic, steric, and short-ranged nonelectrostatic interactions for proteins from low to high concentrations.