Assembly and Analysis of Cell-Scale Membrane Envelopes

Molecular dynamics simulations of biological membranes and membrane-associated phenomena across scales

Membranes are fundamental components of cells that are involved in a wide variety of cellular functions. They are inherently complex, being composed of hugely diverse collections of lipids and proteins, and their various functions arise directly from the intricate interplay between their components. To investigate the interactions between these components in detail, molecular dynamics (MD) simulations have proven to be an invaluable tool. In this mini-review, we highlight several recent studies that illustrate the current state of the art in using MD to study membrane systems. In particular, we first examine how MD is being used to characterize membrane binding of peripheral membrane proteins, we next describe how interactions between lipids and integral membrane proteins are being probed with MD, and we conclude by discussing new tools that have recently been developed to address the significant challenge of building simulateable models of large-scale membranes with complex curvature.

All-atom virus simulations to tackle airborne disease

We briefly review the latest computational studies focused on modeling viruses with classical all-atom (AA) molecular dynamics. We report on the challenges, current solutions, and ongoing developments in constructing and simulating whole viruses, and discuss unique insights derived from AA mesoscale simulations that cannot be achieved by other means. Finally, we present new opportunities in computational virology to understand viral aerostability within the context of respiratory disease transmission. Overall, we highlight the value of large-scale AA simulation and champion the need for increased interdisciplinary collaboration to generate novel insights and guide future research in respiratory disease.

Advancing Molecular Simulations: Merging Physical Models, Experiments, and AI to Tackle Multiscale Complexity

Proteins and protein complexes form adaptable networks that regulate essential biochemical pathways and define cell phenotypes through dynamic mechanisms and interactions. Advances in structural biology and molecular simulations have revealed how protein systems respond to changes in their environments, such as ligand binding, stress conditions, or perturbations like mutations and post-translational modifications, influencing signal transduction and cellular phenotypes. Here, we discuss how computational approaches, ranging from molecular dynamics (MD) simulations to AI-driven methods, are instrumental in studying protein dynamics from isolated molecules to large assemblies. These techniques elucidate conformational landscapes, ligand-binding mechanisms, and protein-protein interactions and are starting to support the construction of multiscale realistic representations of highly complex systems, ranging up to whole cell models. With cryo-electron microscopy, cryo-electron tomography, and AlphaFold accelerating the structural characterization of protein networks, we suggest that integrating AI and Machine Learning with multiscale MD methods will enhance fundamental understating for systems of ever-increasing complexity, usher in exciting possibilities for predictive modeling of the behavior of cell compartments or even whole cells. These advances are indeed transforming biophysics and chemical biology, offering new opportunities to study biomolecular mechanisms at atomic resolution.

Enabling Atomistic Modeling and Simulation of Complex Curved Cellular Membranes with xMAS Builder

As more powerful high performance computing resources are becoming available, there is a new opportunity to bring the unique capabilities of molecular dynamics (MD) simulations to cell-scale systems. Membranes are ubiquitous within cells and are responsible for a diverse set of essential biological functions, but building atomistic models of cell-scale membranes for MD simulations is immensely challenging because of their vast sizes, complex geometries, and complex compositions. To meet this challenge, we have developed xMAS Builder (Experimentally-Derived Membranes of Arbitrary Shape Builder), which is designed to take experimental lipidomics and structural (e.g., electron microscopy and tomography) data as input and use them to build MD-ready models of cellular membrane systems. To test xMAS Builder's capabilities, we have used it to build two models (one ~12.0 million atoms and the other ~11.6 million atoms) of a test system with a representative complex lipid composition and geometry. The two models, which differed only in their lipid packing densities, both maintained their membrane integrity during an extended MD simulation (250 ns and 386 ns), but their highly divergent relaxation dynamics indicate that the proper packing density of curved membranes is determined by leaflet volume rather than surface area. These results suggest that xMAS Builder's algorithms produce high quality models and that simulation of these models will provide profound biophysical insights into the behavior of cellular membranes.

Recent Progress in Modeling and Simulation of Biomolecular Crowding and Condensation Inside Cells

Macromolecular crowding in the cellular cytoplasm can potentially impact diffusion rates of proteins, their intrinsic structural stability, binding of proteins to their corresponding partners as well as biomolecular organization and phase separation. While such intracellular crowding can have a large impact on biomolecular structure and function, the molecular mechanisms and driving forces that determine the effect of crowding on dynamics and conformations of macromolecules are so far not well understood. At a molecular level, computational methods can provide a unique lens to investigate the effect of macromolecular crowding on biomolecular behavior, providing us with a resolution that is challenging to reach with experimental techniques alone. In this review, we focus on the various physics-based and data-driven computational methods developed in the past few years to investigate macromolecular crowding and intracellular protein condensation. We review recent progress in modeling and simulation of biomolecular systems of varying sizes, ranging from single protein molecules to the entire cellular cytoplasm. We further discuss the effects of macromolecular crowding on different phenomena, such as diffusion, protein-ligand binding, and mechanical and viscoelastic properties, such as surface tension of condensates. Finally, we discuss some of the outstanding challenges that we anticipate the community addressing in the next few years in order to investigate biological phenomena in model cellular environments by reproducing in vivo conditions as accurately as possible.

Modeling membranes in situ

Molecular dynamics simulations of cellular membranes have come a long way-from simple model lipid bilayers to multicomponent systems capturing the crowded and complex nature of real cell membranes. In this opinionated minireview, we discuss the current challenge to simulate the dynamics of membranes in their native environment, in situ, with the prospect of reaching the level of whole cells and cell organelles using an integrative modeling framework.

There and back again: bridging meso- and nano-scales to understand lipid vesicle patterning

We describe a complete methodology to bridge the scales between nanoscale molecular dynamics and (micrometer) mesoscale Monte Carlo simulations in lipid membranes and vesicles undergoing phase separation, in which curving molecular species are furthermore embedded. To go from the molecular to the mesoscale, we notably appeal to physical renormalization arguments enabling us to rigorously infer the mesoscale interaction parameters from its molecular counterpart. We also explain how to deal with the physical timescales at stake at the mesoscale. Simulating the as-obtained mesoscale system enables us to equilibrate the long wavelengths of the vesicles of interest, up to the vesicle size. Conversely, we then backmap from the meso- to the nano-scale, which enables us to equilibrate in turn the short wavelengths down to the molecular length-scales. By applying our approach to the specific situation of patterning a vesicle membrane, we show that macroscopic membranes can thus be equilibrated at all length-scales in achievable computational time offering an original strategy to address the fundamental challenge of timescale in simulations of large bio-membrane systems.

VMD as a Platform for Interactive Small Molecule Preparation and Visualization in Quantum and Classical Simulations

Modeling and simulation of small molecules such as drugs and biological cofactors have been both a major focus of computational chemistry for decades and a growing need among computational biophysicists who seek to investigate the interaction of different types of ligands with biomolecules. Of particular interest in this regard are quantum mechanical (QM) calculations that are used to more accurately describe such small molecules, which can be of heterogeneous structures and chemistry, either in purely QM calculations or in hybrid QM/molecular mechanics (MM) simulations. QM programs are also used to develop MM force field parameters for small molecules to be used along with established force fields for biomolecules in classical simulations. With this growing need in mind, here we report a set of software tools developed and closely integrated within the broadly used molecular visualization/analysis program, VMD, that allow the user to construct, modify, and parametrize small molecules and prepare them for QM, hybrid QM/MM, or classical simulations. The tools also provide interactive analysis and visualization capabilities in an easy-to-use and integrated environment. In this paper, we briefly report on these tools and their major features and capabilities, along with examples of how they can facilitate molecular research in computational biophysics that might be otherwise prohibitively complex.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

Mesoscale simulations: An indispensable approach to understand biomembranes

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.

Molecular dynamics simulation of an entire cell

The ultimate microscope, directed at a cell, would reveal the dynamics of all the cell's components with atomic resolution. In contrast to their real-world counterparts, computational microscopes are currently on the brink of meeting this challenge. In this perspective, we show how an integrative approach can be employed to model an entire cell, the minimal cell, JCVI-syn3A, at full complexity. This step opens the way to interrogate the cell's spatio-temporal evolution with molecular dynamics simulations, an approach that can be extended to other cell types in the near future.

LongBondEliminator: A Molecular Simulation Tool to Remove Ring Penetrations in Biomolecular Simulation Systems

We develop a workflow, implemented as a plugin to the molecular visualization program VMD, that can fix ring penetrations with minimal user input. LongBondEliminator, detects ring piercing artifacts by the long, strained bonds that are the local minimum energy conformation during minimization for some assembled simulation system. The LongBondEliminator tool then automatically treats regions near these long bonds using multiple biases applied through NAMD. By combining biases implemented through the collective variables module, density-based forces, and alchemical techniques in NAMD, LongBondEliminator will iteratively alleviate long bonds found within molecular simulation systems. Through three concrete examples with increasing complexity, a lignin polymer, an viral capsid assembly, and a large, highly glycosylated protein aggrecan, we demonstrate the utility for this method in eliminating ring penetrations from classical MD simulation systems. The tool is available via gitlab as a VMD plugin, and has been developed to be generically useful across a variety of biomolecular simulations.

A scoring function for the prediction of protein complex interfaces based on the neighborhood preferences of amino acids

Proteins often assemble into functional complexes, the structures of which are more difficult to obtain than those of the individual protein molecules. Given the structures of the subunits, it is possible to predict plausible complex models via computational methods such as molecular docking. Assessing the quality of the predicted models is crucial to obtain correct complex structures. Here, an energy-scoring function was developed based on the interfacial residues of structures in the Protein Data Bank. The statistically derived energy function (Nepre) imitates the neighborhood preferences of amino acids, including the types and relative positions of neighboring residues. Based on the preference statistics, a program iNepre was implemented and its performance was evaluated with several benchmarking decoy data sets. The results show that iNepre scores are powerful in model ranking to select the best protein complex structures.

A brief history of visualizing membrane systems in molecular dynamics simulations

Understanding lipid dynamics and function, from the level of single, isolated molecules to large assemblies, is more than ever an intensive area of research. The interactions of lipids with other molecules, particularly membrane proteins, are now extensively studied. With advances in the development of force fields for molecular dynamics simulations (MD) and increases in computational resources, the creation of realistic and complex membrane systems is now common. In this perspective, we will review four decades of the history of molecular dynamics simulations applied to membranes and lipids through the prism of molecular graphics.

Activation energy and force fields during topological transitions of fluid lipid vesicles

Topological transitions of fluid lipid membranes are fundamental processes for cell life. For example, they are required for endo- and exocytosis or to enable neurotransmitters to cross the neural synapses. Here, inspired by the idea that fusion and fission proteins could have evolved in Nature in order to carry out a minimal work expenditure, we evaluate the minimal free energy pathway for the transition between two spherical large unilamellar vesicles and a dumbbell-shaped one. To address the problem, we propose and successfully use a Ginzburg-Landau type of free energy, which allows us to uniquely describe without interruption the whole, full-scale topological change. We also compute the force fields needed to overcome the involved energy barriers. The obtained forces are in excellent agreement, in terms of intensity, scale, and spatial localization with experimental data on typical fission protein systems, whereas they suggest the presence of additional features in fusion proteins.

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Organelle Interaction and Drug Discovery: Towards Correlative Nanoscopy and Molecular Dynamics Simulation

The inter-organelle interactions, including the cytomembrane, endoplasmic reticulum, mitochondrion, lysosome, dictyosome, and nucleus, play the important roles in maintaining the normal function and homeostasis of cells. Organelle dysfunction can lead to a range of diseases (e.g., Alzheimer's disease (AD), Parkinson's disease (PD), and cancer), and provide a new perspective for drug discovery. With the development of imaging techniques and functional fluorescent probes, a variety of algorithms and strategies have been developed for the ever-improving estimation of subcellular structures, organelle interaction, and organelle-related drug discovery with accounting for the dynamic structures of organelles, such as the nanoscopy technology and molecular dynamics (MD) simulations. Accordingly, this work summarizes a series of state-of-the-art examples of the recent progress in this rapidly changing field and uncovering the drug screening based on the structures and interactions of organelles. Finally, we propose the future outlook for exciting applications of organelle-related drug discovery, with the cooperation of nanoscopy and MD simulations.

Integrating cellular and molecular structures and dynamics into whole-cell models

A complete description of the state of the cell requires knowledge of its size, shape, components, intracellular reactions, and interactions with its environment-all of these as a function of time and cell growth. Adding to this list is the need for theoretical models and simulations that integrate and help to interpret this daunting amount of experimental data. It seems like an overwhelming list of requirements, but progress is being made on many fronts. In this review, we discuss the current challenges and problems in obtaining sufficient information about each aspect of a dynamical whole-cell model (DWCM) for simple and well-studied bacterial systems.