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.

Recent Advances in Modeling Membrane ß-Barrel Proteins Using Molecular Dynamics Simulations: From Their Lipid Environments to Their Assemblies.

Spurred by advances in AI-driven modeling and experimental methods, molecular dynamics simulations are now acting as a platform to integrate these different approaches. This combination of methods is especially useful to understand ß-barrel proteins from the molecular level, e.g., identifying specific interactions with lipids or small molecules, up to assemblies comprised of hundreds of proteins and thousands of lipids. In this minireview, we will discuss recent advances, mainly from the last 5 years, in modeling ß-barrel proteins and their assemblies. These approaches require specific kinds of modeling and potentially different model resolutions that we will first describe in Subheading 1. We will then focus on different aspects of ß-barrel protein modeling: how different types of molecules can diffuse through ß-barrel proteins (Subheading 2); how lipids can interact with these proteins (Subheading 3); how ß-barrel proteins can interact with membrane partners (Subheading 4) or periplasmic extensions and partners (Subheading 5) to form large assemblies.

Molecular Modeling and Simulation of the Mycobacterial Cell Envelope: From Individual Components to Cell Envelope Assemblies.

Unlike typical Gram-positive bacteria, the cell envelope of mycobacteria is unique and composed of a mycobacterial outer membrane, also known as the mycomembrane, a peptidoglycan layer, and a mycobacterial inner membrane, which is analogous to that of Gram-negative bacteria. Despite its importance, however, our understanding of this complex cell envelope is rudimentary at best. Thus, molecular modeling and simulation of such an envelope can benefit the scientific community by proposing new hypotheses about the biophysical properties of its different layers. In this Perspective, we present recent advances in molecular modeling and simulation of the mycobacterial cell envelope from individual components to cell envelope assemblies. We also show how modeling other types of cell envelopes, such as that of Escherichia coli, may help modeling part of the mycobacterial envelopes. We hope that the studies presented here are just the beginning of the road and more and more new modeling and simulation studies help us to understand crucial questions related to mycobacteria such as antibiotic resistance or bacterial survival in the host.

MDverse: Shedding Light on the Dark Matter of Molecular Dynamics Simulations.

The rise of open science and the absence of a global dedicated data repository for molecular dynamics (MD) simulations has led to the accumulation of MD files in generalist data repositories, constituting the dark matter of MD - data that is technically accessible, but neither indexed, curated, or easily searchable. Leveraging an original search strategy, we found and indexed about 250,000 files and 2,000 datasets from Zenodo, Figshare and Open Science Framework. With a focus on files produced by the Gromacs MD software, we illustrate the potential offered by the mining of publicly available MD data. We identified systems with specific molecular composition and were able to characterize essential parameters of MD simulation, such as temperature and simulation length, and identify model resolution, such as all-atom and coarse-grain. Based on this analysis, we inferred metadata to propose a search engine prototype to explore collected MD data. To continue in this direction, we call on the community to pursue the effort of sharing MD data, and increase populating and standardizing metadata to reuse this valuable matter.

Supramolecular organization and dynamics of mannosylated phosphatidylinositol lipids in the mycobacterial plasma membrane.

Mycobacterium tuberculosis (Mtb) is the causative agent of tuberculosis (TB), a disease that claims ~1.6 million lives annually. The current treatment regime is long and expensive, and missed doses contribute to drug resistance. Therefore, development of new anti-TB drugs remains one of the highest public health priorities. Mtb has evolved a complex cell envelope that represents a formidable barrier to antibiotics. The Mtb cell envelop consists of four distinct layers enriched for Mtb specific lipids and glycans. Although the outer membrane, comprised of mycolic acid esters, has been extensively studied, less is known about the plasma membrane, which also plays a critical role in impacting antibiotic efficacy. The Mtb plasma membrane has a unique lipid composition, with mannosylated phosphatidylinositol lipids (phosphatidyl-myoinositol mannosides, PIMs) comprising more than 50% of the lipids. However, the role of PIMs in the structure and function of the membrane remains elusive. Here, we used multiscale molecular dynamics (MD) simulations to understand the structure-function relationship of the PIM lipid family and decipher how they self-organize to shape the biophysical properties of mycobacterial plasma membranes. We assess both symmetric and asymmetric assemblies of the Mtb plasma membrane and compare this with residue distributions of Mtb integral membrane protein structures. To further validate the model, we tested known anti-TB drugs and demonstrated that our models agree with experimental results. Thus, our work sheds new light on the organization of the mycobacterial plasma membrane. This paves the way for future studies on antibiotic development and understanding Mtb membrane protein function.

GPC3-Unc5 receptor complex structure and role in cell migration.

Neural migration is a critical step during brain development that requires the interactions of cell-surface guidance receptors. Cancer cells often hijack these mechanisms to disseminate. Here, we reveal crystal structures of Uncoordinated-5 receptor D (Unc5D) in complex with morphogen receptor glypican-3 (GPC3), forming an octameric glycoprotein complex. In the complex, four Unc5D molecules pack into an antiparallel bundle, flanked by four GPC3 molecules. Central glycan-glycan interactions are formed by N-linked glycans emanating from GPC3 (N241 in human) and C-mannosylated tryptophans of the Unc5D thrombospondin-like domains. MD simulations, mass spectrometry and structure-based mutants validate the crystallographic data. Anti-GPC3 nanobodies enhance or weaken Unc5-GPC3 binding and, together with mutant proteins, show that Unc5/GPC3 guide migrating pyramidal neurons in the mouse cortex, and cancer cells in an embryonic xenograft neuroblastoma model. The results demonstrate a conserved structural mechanism of cell guidance, where finely balanced Unc5-GPC3 interactions regulate cell migration.

Design - a new way to look at old molecules.

We discuss how design enriches molecular science, particularly structural biology and bioinformatics. We present two use cases, one in academic practice and the other to design for outreach. The first case targets the representation of ion channels and their dynamic properties. In the second, we document a transition process from a research environment to general-purpose designs. Several testimonials from practitioners are given. By describing the design process of abstracted shapes, exploded views of molecular structures, motion-averaged slices, 360-degree panoramic projections, and experiments with lit sphere shading, we document how designers help make scientific data accessible without betraying its meaning, and how a creative mind adds value over purely data-driven visualizations. A similar conclusion was drawn for public outreach, as we found that comic-book-style drawings are better suited for communicating science to a broad audience.

The guidance and adhesion protein FLRT2 dimerizes in cis via dual small-X3-small transmembrane motifs.

Fibronectin Leucine-rich Repeat Transmembrane (FLRT 1-3) proteins are a family of broadly expressed single-spanning transmembrane receptors that play key roles in development. Their extracellular domains mediate homotypic cell-cell adhesion and heterotypic protein interactions with other receptors to regulate cell adhesion and guidance. These in trans FLRT interactions determine the formation of signaling complexes of varying complexity and function. Whether FLRTs also interact at the surface of the same cell, in cis, remains unknown. Here, molecular dynamics simulations reveal two dimerization motifs in the FLRT2 transmembrane helix. Single particle tracking experiments show that these Small-X3-Small motifs synergize with a third dimerization motif encoded in the extracellular domain to permit the cis association and co-diffusion patterns of FLRT2 receptors on cells. These results may point to a competitive switching mechanism between in cis and in trans interactions, which suggests that homotypic FLRT interaction mirrors the functionalities of classic adhesion molecules.

Proteasome complexes experience profound structural and functional rearrangements throughout mammalian spermatogenesis.

During spermatogenesis, spermatogonia undergo a series of mitotic and meiotic divisions on their path to spermatozoa. To achieve this, a succession of processes requiring high proteolytic activity are in part orchestrated by the proteasome. The spermatoproteasome (s20S) is specific to the developing gametes, in which the gamete-specific α4s subunit replaces the α4 isoform found in the constitutive proteasome (c20S). Although the s20S is conserved across species and was shown to be crucial for germ cell development, its mechanism, function, and structure remain incompletely characterized. Here, we used advanced mass spectrometry (MS) methods to map the composition of proteasome complexes and their interactomes throughout spermatogenesis. We observed that the s20S becomes highly activated as germ cells enter meiosis, mainly through a particularly extensive 19S activation and, to a lesser extent, PA200 binding. Additionally, the proteasome population shifts from c20S (98%) to s20S (>82 to 92%) during differentiation, presumably due to the shift from α4 to α4s expression. We demonstrated that s20S, but not c20S, interacts with components of the meiotic synaptonemal complex, where it may localize via association with the PI31 adaptor protein. In vitro, s20S preferentially binds to 19S and displays higher trypsin- and chymotrypsin-like activities, both with and without PA200 activation. Moreover, using MS methods to monitor protein dynamics, we identified significant differences in domain flexibility between α4 and α4s. We propose that these differences induced by α4s incorporation result in significant changes in the way the s20S interacts with its partners and dictate its role in germ cell differentiation.

Visualizing protein structures - tools and trends.

Molecular visualization is fundamental in the current scientific literature, textbooks and dissemination materials. It provides an essential support for presenting results, reasoning on and formulating hypotheses related to molecular structure. Tools for visual exploration of structural data have become easily accessible on a broad variety of platforms thanks to advanced software tools that render a great service to the scientific community. These tools are often developed across disciplines bridging computer science, biology and chemistry. This mini-review was written as a short and compact overview for scientists who need to visualize protein structures and want to make an informed decision which tool they should use. Here, we first describe a few 'Swiss Army knives' geared towards protein visualization for everyday use with an existing large user base, then focus on more specialized tools for peculiar needs that are not yet as broadly known. Our selection is by no means exhaustive, but reflects a diverse snapshot of scenarios that we consider informative for the reader. We end with an account of future trends and perspectives.

Structural Basis of Teneurin-Latrophilin Interaction in Repulsive Guidance of Migrating Neurons.

Teneurins are ancient metazoan cell adhesion receptors that control brain development and neuronal wiring in higher animals. The extracellular C terminus binds the adhesion GPCR Latrophilin, forming a trans-cellular complex with synaptogenic functions. However, Teneurins, Latrophilins, and FLRT proteins are also expressed during murine cortical cell migration at earlier developmental stages. Here, we present crystal structures of Teneurin-Latrophilin complexes that reveal how the lectin and olfactomedin domains of Latrophilin bind across a spiraling beta-barrel domain of Teneurin, the YD shell. We couple structure-based protein engineering to biophysical analysis, cell migration assays, and in utero electroporation experiments to probe the importance of the interaction in cortical neuron migration. We show that binding of Latrophilins to Teneurins and FLRTs directs the migration of neurons using a contact repulsion-dependent mechanism. The effect is observed with cell bodies and small neurites rather than their processes. The results exemplify how a structure-encoded synaptogenic protein complex is also used for repulsive cell guidance.

The conical shape of DIM lipids promotes Mycobacterium tuberculosis infection of macrophages.

Phthiocerol dimycocerosate (DIM) is a major virulence factor of the pathogen Mycobacterium tuberculosis (Mtb). While this lipid promotes the entry of Mtb into macrophages, which occurs via phagocytosis, its molecular mechanism of action is unknown. Here, we combined biophysical, cell biology, and modeling approaches to reveal the molecular mechanism of DIM action on macrophage membranes leading to the first step of Mtb infection. Matrix-assisted laser desorption ionization time-of-flight (MALDI-TOF) mass spectrometry showed that DIM molecules are transferred from the Mtb envelope to macrophage membranes during infection. Multiscale molecular modeling and 31P-NMR experiments revealed that DIM adopts a conical shape in membranes and aggregates in the stalks formed between 2 opposing lipid bilayers. Infection of macrophages pretreated with lipids of various shapes uncovered a general role for conical lipids in promoting phagocytosis. Taken together, these results reveal how the molecular shape of a mycobacterial lipid can modulate the biological response of macrophages.

Molecular Graphics: Bridging Structural Biologists and Computer Scientists.

Visualization of molecular structures is one of the most common tasks carried out by structural biologists, typically using software, such as Chimera, COOT, PyMOL, or VMD. In this Perspective article, we outline how past developments in computer graphics and data visualization have expanded the understanding of biomolecular function, and we summarize recent advances that promise to further transform structural biology. We also highlight how progress in molecular graphics has been impeded by communication barriers between two communities: the computer scientists driving these advances, and the structural and computational biologists who stand to benefit. By pointing to canonical papers and explaining technical progress underlying new graphical developments in simple terms, we aim to improve communication between these communities; this, in turn, would help shape future developments in molecular graphics.

Sharing Data from Molecular Simulations.

Given the need for modern researchers to produce open, reproducible scientific output, the lack of standards and best practices for sharing data and workflows used to produce and analyze molecular dynamics (MD) simulations has become an important issue in the field. There are now multiple well-established packages to perform molecular dynamics simulations, often highly tuned for exploiting specific classes of hardware, each with strong communities surrounding them, but with very limited interoperability/transferability options. Thus, the choice of the software package often dictates the workflow for both simulation production and analysis. The level of detail in documenting the workflows and analysis code varies greatly in published work, hindering reproducibility of the reported results and the ability for other researchers to build on these studies. An increasing number of researchers are motivated to make their data available, but many challenges remain in order to effectively share and reuse simulation data. To discuss these and other issues related to best practices in the field in general, we organized a workshop in November 2018 ( https://bioexcel.eu/events/workshop-on-sharing-data-from-molecular-simulations/ ). Here, we present a brief overview of this workshop and topics discussed. We hope this effort will spark further conversation in the MD community to pave the way toward more open, interoperable, and reproducible outputs coming from research studies using MD simulations.

Cholesterol Interaction Sites on the Transmembrane Domain of the Hedgehog Signal Transducer and Class F G Protein-Coupled Receptor Smoothened.

Transduction of Hedgehog signals across the plasma membrane is facilitated by the class F G-protein-coupled-receptor (GPCR) Smoothened (SMO). Recent studies suggest that SMO is modulated via interactions of its transmembrane (TM) domain with cholesterol. We apply molecular dynamics simulations of SMO embedded in cholesterol containing lipid bilayers, revealing a direct interaction of cholesterol with the TM domain at regions distinct from those observed in class A GPCRs. In particular the extracellular tips of helices TM2 and TM3 form a well-defined cholesterol interaction site. Potential of mean force calculations yield a free energy landscape for cholesterol binding. Alongside analysis of equilibrium cholesterol occupancy, this reveals the existence of a dynamic "greasy patch" interaction with the TM domain of SMO, which may be compared with previously identified lipid interaction sites on other membrane proteins. These predictions provide molecular-level insights into cholesterol interactions with a class F GPCR, suggesting potential druggable sites.

How nanoscale protein interactions determine the mesoscale dynamic organisation of bacterial outer membrane proteins.

The spatiotemporal organisation of membranes is often characterised by the formation of large protein clusters. In Escherichia coli, outer membrane protein (OMP) clustering leads to OMP islands, the formation of which underpins OMP turnover and drives organisation across the cell envelope. Modelling how OMP islands form in order to understand their origin and outer membrane behaviour has been confounded by the inherent difficulties of simulating large numbers of OMPs over meaningful timescales. Here, we overcome these problems by training a mesoscale model incorporating thousands of OMPs on coarse-grained molecular dynamics simulations. We achieve simulations over timescales that allow direct comparison to experimental data of OMP behaviour. We show that specific interaction surfaces between OMPs are key to the formation of OMP clusters, that OMP clusters present a mesh of moving barriers that confine newly inserted proteins within islands, and that mesoscale simulations recapitulate the restricted diffusion characteristics of OMPs.

Interactions of the EphA2 Kinase Domain with PIPs in Membranes: Implications for Receptor Function.

EphA2 is a member of the receptor tyrosine kinase family. Interactions of the cytoplasmic region of EphA2 with the cell membrane are functionally important and yet remain incompletely characterized. Molecular dynamics simulations combined with biochemical studies reveal the interactions of the transmembrane, juxtamembrane (JM), and kinase domains with the membrane. We describe how the kinase domain is oriented relative to the membrane and how the JM region can modulate this interaction. We highlight the role of phosphatidylinositol phosphates (PIPs) in mediating the interaction of the kinase domain with the membrane and, conversely, how positively charged patches at the kinase surface and in the JM region induce the formation of nanoclusters of PIP molecules in the membrane. Integration of these results with those from previous studies enable computational reconstitution of a near complete EphA2 receptor within a membrane, suggesting a role for receptor-lipid interactions in modulation of EphA2.

Protein crowding and lipid complexity influence the nanoscale dynamic organization of ion channels in cell membranes.

Cell membranes are crowded and complex environments. To investigate the effect of protein-lipid interactions on dynamic organization in mammalian cell membranes, we have performed coarse-grained molecular dynamics simulations containing >100 copies of an inwardly rectifying potassium (Kir) channel which forms specific interactions with the regulatory lipid phosphatidylinositol 4,5-bisphosphate (PIP2). The tendency of protein molecules to cluster has the effect of organizing the membrane into dynamic compartments. At the same time, the diversity of lipids present has a marked effect on the clustering behavior of ion channels. Sub-diffusion of proteins and lipids is observed. Protein crowding alters the sub-diffusive behavior of proteins and lipids such as PIP2 which interact tightly with Kir channels. Protein crowding also affects bilayer properties, such as membrane undulations and bending rigidity, in a PIP2-dependent manner. This interplay between the diffusion and the dynamic organization of Kir channels may have important implications for channel function.

Lipid-Loving ANTs: Molecular Simulations of Cardiolipin Interactions and the Organization of the Adenine Nucleotide Translocase in Model Mitochondrial Membranes.

The exchange of ADP and ATP across the inner mitochondrial membrane is a fundamental cellular process. This exchange is facilitated by the adenine nucleotide translocase, the structure and function of which are critically dependent on the signature phospholipid of mitochondria, cardiolipin (CL). Here we employ multiscale molecular dynamics simulations to investigate CL interactions within a membrane environment. Using simulations at both coarse-grained and atomistic resolutions, we identify three CL binding sites on the translocase, in agreement with those seen in crystal structures and inferred from nuclear magnetic resonance measurements. Characterization of the free energy landscape for lateral lipid interaction via potential of mean force calculations demonstrates the strength of interaction compared to those of binding sites on other mitochondrial membrane proteins, as well as their selectivity for CL over other phospholipids. Extending the analysis to other members of the family, yeast Aac2p and mouse uncoupling protein 2, suggests a degree of conservation. Simulation of large patches of a model mitochondrial membrane containing multiple copies of the translocase shows that CL interactions persist in the presence of protein-protein interactions and suggests CL may mediate interactions between translocases. This study provides a key example of how computational microscopy may be used to shed light on regulatory lipid-protein interactions.

Membrane stiffness is modified by integral membrane proteins.

The ease with which a cell membrane can bend and deform is important for a wide range of biological functions. Peripheral proteins that induce curvature in membranes (e.g. BAR domains) have been studied for a number of years. Little is known, however, about the effect of integral membrane proteins on the stiffness of a membrane (characterised by the bending rigidity, Kc). We demonstrate by computer simulation that adding integral membrane proteins at physiological densities alters the stiffness of the membrane. First we establish that the coarse-grained MARTINI forcefield is able to accurately reproduce the bending rigidity of a small patch of 1500 phosphatidyl choline lipids by comparing the calculated value to both experiment and an atomistic simulation of the same system. This enables us to simulate the dynamics of large (ca. 50 000 lipids) patches of membrane using the MARTINI coarse-grained description. We find that altering the lipid composition changes the bending rigidity. Adding integral membrane proteins to lipid bilayers also changes the bending rigidity, whilst adding a simple peripheral membrane protein has no effect. Our results suggest that integral membrane proteins can have different effects, and in the case of the bacterial outer membrane protein, BtuB, the greater the density of protein, the larger the reduction in stiffness.