Abstract
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.
Introduction
Biomolecular functions emerge from the molecular interactions taking place in cellular environments. Understanding each component’s role in driving cell function poses an immense challenge. For a long time, experimental techniques have been our main window into the cellular environment. By resolving biomolecular structures and probing the dynamics of biomolecular processes, both in vivo and in vitro, researchers have developed a global understanding of how a cell functions.
A limitation of these experimental techniques is the spatio-temporal resolution that they can probe, particularly within the complex and crowded environment of the cell. Molecular dynamics (MD) simulations provide a suitable alternative approach, covering the relevant length and timescales at molecular resolution, albeit over short periods of a cell cycle. Over the past decades, MD has matured into a powerful tool that functions as a computational microscope (; ). With the advances in available computer power, including the transition from using central processing units (CPUs) to graphical processing units (GPUs), the complexity and the spatio-temporal scales of MD simulations have increased remarkably. State-of-the-art simulations, containing hundreds of millions of atoms, include dynamic models of a photosynthetic chromatophore vesicle (Singharoy et al., 2019), the nuclear core complex (Mosalaganti et al., 2022), the membranes of a mitochondrion (Pezeshkian et al., 2020), the bacterial cytoplasm (Yu et al., 2016), and a virus particle embedded in a nanoscale aerosol droplet (). The natural next step is, arguably, the scale of entire cells (; ; ; Thornburg et al., 2022).
Creating a whole-cell model has long been a major goal in computational modeling. A computational cell will help us to resolve a more integral picture of how biomolecular interactions drive cell function since biomolecular processes operate on a hierarchy of interconnected scales. Thus, resolving the full cell function requires a holistic approach. The current state-of-the-art uses static representations of heterogeneous cell-scale structures such as cellPACK (; Maritan et al., 2022), genome-scale well-stirred reaction models for metabolism (; Macklin et al., 2014; ; ; Macklin et al., 2020), or mesoscale kinetic models that attempt to include both structural and chemical states of the cell such as Lattice Microbes (Roberts et al., 2013). These computational techniques, despite granting impressive insights into the complexity of cellular processes, are limited by the spatial resolution they can achieve.
Constructing whole-cell models requires the integration of a large amount of experimental data, i.e., relies on an integrative modeling approach (; ; ). Obtaining such data with high spatial and dynamic detail, particularly in living cells, is very challenging, but exciting progress is being made in elucidating the architecture and stoichiometry of cellular components at unprecedented resolutions (Reading et al., 2017; ; ; ; ; ; Narasimhan et al., 2020; Wietrzynski et al., 2020; Štefl et al., 2020), setting the stage for spatially detailed simulations of whole cells. To showcase this possibility, we consider one of the simplest cells known to date: the minimal cell created by the J. Craig Venter Institute (), a stripped-down version of a Mycoplasma bacterium. The current strain, named JCVI-syn3A (Syn3A), contains only 493 genes and is still able to replicate independently (). This cell is particularly amenable to detailed computational modeling because it is of relatively small size (measuring 400 nm in diameter), and its precise composition has largely been resolved ().
Here we present our ongoing effort to simulate Syn3A using the Martini coarse-grained (CG) force field (Marrink et al., 2022). The Martini force field employs a four-to-one mapping scheme, where up to four heavy atoms and associated hydrogens are represented by one CG bead. This reduction in the number of particles in the system, together with a smoothening of the potential energy surface, speeds up the simulations by about three orders of magnitude, enabling simulations of systems that approach the size of entire cells. In the case of Syn3A, about 550 million CG particles are required, corresponding to more than six billion atoms.
In the remainder of this work, we first describe the set of tools needed to construct a system as complex as an entire cell at the Martini level (the Martini “ecosystem”), including a description of the stepwise procedure to construct a starting model for Syn3A, from building the chromosome and modeling the cytoplasm to the addition of the cell envelope. We end with the prospects of actually simulating this model and discuss the potential future avenues of simulating entire cells. The integrative modeling workflow is schematically depicted in Figure 1.
FIGURE 1
Building cells with the Martini ecosystem
Modeling cellular environments using a coarse-grained approach requires the use of a force field that incorporates enough detail to represent all biomolecules and their interactions explicitly. In this regard, the Martini force field is an excellent candidate, as demonstrated by the wide range of application studies using Martini over the past 2 decades (Marrink et al., 2022). Additionally, parameters for a large variety of biomolecules are already available, including proteins (), lipids (Wassenaar et al., 2015), polynucleotides (Uusitalo et al., 2015; Uusitalo et al., 2017), carbohydrates (; ) and metabolites (Sousa et al., 2021; ). A curated collection of all parameters is available from the Martini Database ().
Accompanying the Martini force field is a collection of tools that compose the software side of the Martini ecosystem (Figure 1). The primary goal of this software suite is to facilitate the construction of topologies and initial coordinates for running CG Martini MD simulations. The Martini ecosystem is currently being extended around a central framework, named Vermouth. Making use of a graph-based description of molecules, Vermouth implements a unified handling of processes frequently used in Martini, such as topology and coordinate generation or resolution transformation, as a lightweight python library ().
Proteins comprise the bulk of a cell’s organic material, composing approximately 40% of the total intracellular volume (). The number of unique proteins expressed by the cell, i.e., the proteome, can range from a few hundred to several thousand. Consequently, describing realistic cellular environments requires generating topology files for proteins of varying shapes and sizes as well as packing these into a highly crowded solution. Martinize2, which is built on top of Vermouth, is designed for high-throughput generation of Martini topologies and coordinates for proteins from atomistic protein structures. The workflow used in Martinize2 additionally performs quality checks on the atomistic protein structures and alerts to potential problems making it ideal for such a high-complexity use case (). To generate dense protein solutions on a cellular scale as required for whole-cell modeling, a new software tool, called Bentopy, is currently under development. It uses an efficient collision detection scheme () to generate random packings of proteins and protein complexes within volumetric constraints. Furthermore, functional annotations of proteins can be integrated into the algorithm, biasing their spatial distribution in the cytosol based on their known biochemical function.
Constructing coordinates and input files for chromosomal DNA presents another challenge for modeling of a whole cell. Even for a comparatively small genome as that of JCVI-syn3A, approaches that rely on reading input coordinates and forward mapping such as used in Martinize2, become too inefficient due to the sheer size of the molecule. Instead the Polyply software, which was originally developed to efficiently setup general polymeric systems, can be used (). Within Polyply, a multiresolution graph-based approach is used to efficiently generate polymer topologies, in particular for DNA, only from the sequence. In addition to topologies, system coordinates can be generated using a specialized biased random walk protocol. This tool of the Martini ecosystem has successfully been applied to model dense polymer melts and simple ssDNA viral chromosomes. At the moment, the package is being extended to handle double-stranded nucleic acids, and generate more complex DNA structures such as bacterial chromosomes.
Lastly, modeling lipid membranes has historically been a leading application of the Martini force field (Marrink et al., 2019). Simulating arbitrarily complex membranes of various sizes, geometries, and lateral heterogeneities is facilitated by TS2CG (Pezeshkian et al., 2020). This tool implements a backmapping algorithm that converts triangulated surfaces into CG membrane models. As a result of the method’s high level of control, the curvature-dependent lipid concentrations in both membrane leaflets can be precisely determined by the user. In addition, proteins can be inserted into the membrane together with their characteristic lipid shells [i.e., lipid fingerprints ()], setting the stage for building cell envelopes.
In the following subsections, we describe the application of the aforementioned tools to construct a proof of principle whole-cell simulation of Syn3A, illustrating how the current Martini ecosystem enables users to study multi-component systems at the mesoscale.
Chromosome building
The minimal genome of JCVI-syn3A contains 493 genes and is encoded in a single circular chromosome of 543 kilobase pairs (kbps). Since the chromosome is contained inside the cell’s cytosol, the structural organization is heavily influenced by the crowded intracellular environment. Due to the size and near-uniform distribution of ribosomes present in the cytosol, the excluded volume interactions of these protein-RNA complexes are known to have a significant influence on the nucleoid organization (Mondal et al., 2011).
The nucleoid structure of Syn3A was previously modeled by based on the ribosome distribution and cell boundary determined by cryo-electron tomography. A Monte Carlo (MC) method grew the chromosome, modeled by a self-avoiding polygon, on a lattice inside the cell boundary. Each MC step ensured that no model constraints were violated, resulting in a circular genome without steric clashes with the ribosomes or cell membrane. The algorithm was validated by comparing the chromosome conformation capture (3C) maps of ensembles of simulated nucleoid configurations with experimental 3C maps. 3C maps show spatial correlations between chromosomal regions, which are spatially close but can be distant in the nucleotide sequence. Based on the features in the 3C maps, we infer that the chromosome is organized more like a fractal globule with little persistent supercoiling.
Whilst the previous chromosome modeling approach with a lattice polymer was tailored to be highly compatible with the whole-cell simulations using Lattice Microbes, we have subsequently developed a new method to generate circular chromosomes organized as fractal globules in a continuum polymer model with 10 bp monomers. The generated chromosome model is relaxed using Brownian dynamics and an energy function for modeling dsDNA as a twistable worm-like chain from (). In order to connect the chromosome model to a Martini-level representation, the model is transformed to a one-bead-per-base-pair resolution by spline interpolation. Rotation minimizing frames are then constructed along the chromosomal contour, providing a consistent reference to which the Martini DNA model can be backmapped (Wang et al., 2008). After adding an equilibrium twist along the frame’s tangent vector, Martini base pair templates matching the 543 kbp genome sequence are positioned along the chromosome following the local contour reference frame. By performing a short energy minimization the system is relaxed, resulting in a stable chromosome structure. The subsequent model consists of 543 kbps, which at a Martini resolution is equivalent to seven million beads. By implementing this backmapping procedure in Polyply, we are able to efficiently generate the coordinates for the chromosome in a force field agnostic manner. The overall chromosome building takes a matter of minutes, opening up the possibility of studying larger protein-DNA complexes like chromatin fibers and Escherichia coli chromosomes. The required topology files were generated from the sequence using the default Polyply methods.
Cytosol modeling
In order to model the cytosol, it is essential to have a complete picture of the bacterial proteome, including both protein structures and proteomics counts. The genome reduction leading to Syn3A limits the number of different proteins that have to be taken into account by only retaining 452 protein-coding genes. This minimal genome has been extensively characterized, and only 91 genes remain without an annotated function. A recent study by uses computational analyses to further elucidate the function of uncharacterized genes and work toward complete functional characterization of the proteome. By gaining a better understanding of the function of encoded proteins, we will be able to inform the spatial distribution of proteins in our whole-cell model.
From the 452 different proteins expressed by Syn3A, 281 are characterized as cytosolic proteins, 63 as trans-membrane proteins, 42 as peripheral membrane proteins, and the remaining 66 still have an unknown localization. As part of the computational gene characterization workflow, Bianchi et al. modeled the protein structures of the entire proteome using AlphaFold2 (). Martinize2 successfully converted all but one of the predicted protein structures (451) to a corresponding Martini model. Using Bentopy, the cytosolic protein models are packed into the intracellular volume alongside the chromosome and ribosomes. The number of copies of each protein is based on available proteomics data (; Thornburg et al., 2022); in total, around 60,000 proteins were distributed within a spherical volume with a diameter of 400 nm. Concerning the ribosomes, we used bacterial homologs that we had already generated previously (Uusitalo et al., 2017), placing 503 ribosomes in random orientations near the positions originally determined from the cryo-electron tomography map (). Single-stranded RNA fragments were not included at this stage.
The next major component of the cytosol are the small molecules that, together with enzymatic proteins, participate in the metabolic pathways. In the current model, we include only the metabolites for which Martini topologies were already available, primarily amino acids and nucleotide cofactors (Sousa et al., 2021), and which are present at high concentrations inside Syn3A. The metabolite models were automatically generated from the topology files using Polyply. Based on available metabolomic data (Thornburg et al., 2022), 1.7 million metabolites are distributed within the cytosol, approximately 55% of the metabolite count for the complete metabolome.
Constructing the envelope
Modeling the cell envelope of the Syn3A is a straightforward procedure since it is solely composed of a singular cytoplasmic membrane. Furthermore, experimental measurements indicate the absence of a cell capsule, drastically reducing the complexity of the cell boundary. The lipid membrane is constructed using TS2CG with a uniform lipid mixture across both membrane leaflets. It should be noted that since the minimal cell acquires membrane components through lipid synthesis from fatty acids and direct incorporation of lipids from its environment, the lipid composition of the cellular membrane heavily depends on the growth medium. We base our model on the lipidomics data presented by (Thornburg et al., 2022), indicating the presence of five main lipid types: cholesterol (59%), sphingomyelin (18%), cardiolipin (17%), phosphatidylcholines (4%), and phosphatidylglycerol (2%). In the absence of more detailed lipidomics data, all lipids are modeled with fully saturated palmitoyl tails. The total lipid count amounted to 1.3 million lipids.
Additionally, we randomly inserted membrane proteins into the cell membrane using TS2CG. From the available proteomics data, the number and types of membrane proteins are determined. While AlphaFold2 structure predictions can be used directly to model monomeric membrane proteins, experimental crystal structures are still required for the protein transport complexes. Martinize2 is again used to generate the Martini models for the membrane proteins. For simplicity, we selected five abundant protein complexes and distributed these uniformly over the membrane. In total, 2,200 protein complexes were embedded in the cell envelope, corresponding with the expected number of membrane proteins present on the surface of Syn3A.
Solvating and simulating the cell
Having modeled all the cell components, the final step in constructing a starting structure for subsequent simulation is defining the periodic simulation box and solvating the system. Considering the whole-cell model’s spherical shape, a logical choice for the periodic box is a rhombic dodecahedron. To solvate, a periodic water box is tiled across the cell model, removing the water beads that overlap with the model using a collision detection scheme. The system is neutralized by placing counter ions near the highly charged components in the cytosol, i.e., the chromosome and ribosomes; the overall negative charge is substantial, amounting to 3.2 million elementary charges. As part of the solvation procedure, we also replace an appropriate number of water beads with ion beads to establish an ion concentration of 135 mM NaCl across our system, mimicking the experimental buffer. Thus, we ended up with a system containing 447 million water beads (208 million inside, 239 outside of the cell), 8.5 million sodium, and 5.3 million chloride ions. Note that Martini CG water beads represent four real water molecules. The total bead count, including all biomolecules, adds up to 561 million beads. A snapshot of the full system is shown in Figure 2.
FIGURE 2
Having constructed a starting model for Syn3A, the current challenge is to perform an actual MD simulation. At the time being, this proved to be non-trivial. Gromacs (), the main MD engine to run Martini-based simulations, is having difficulties handling systems comprising hundreds of millions of particles, in particular featuring large molecules such as the genome spread over multiple domains. The Gromacs developer team is aware of this problem and is dedicated to solving it. Possible other software engines to consider are ddcMD (Zhang et al., 2020) and openMM (), both of which are supporting Martini and offer simulation speeds comparable to those of Gromacs.
Discussion
In the wake of a continuous rise in computing power, MD simulations have transitioned from studying idealized representations of biomolecular systems to modeling their full complexity. The culmination of this development would be simulations at the level of entire cells. As a proof of principle that we are ready to meet this challenge, we presented a model of the complete minimal cell JCVI-syn3A, constructed using the Martini ecosystem. The final simulation box comprises more than 560 million CG beads, representing over six billion atoms in the cell (Figure 2).
Before looking at the broader prospects of this endeavor, it is important to discuss a number of limitations of our approach. The current model uses the Martini 2 version of the force field since Martini models for nucleic acids, and other essential cellular components are still under development for the latest Martini 3 release. However, the methods described in this paper can be straightforwardly transferred to the latest version of Martini when validated models become available. With over 800 different bead types and a recalibrated interaction matrix, Martini 3 offers an improved framework for CG MD simulations (Souza et al., 2021). Nevertheless, inherent limitations of Martini, such as an inability to sample protein secondary conformational changes, remain. We do not anticipate that such changes are of primary importance in determining the cellular organization, but details of protein-protein and protein-lipid interactions might be affected. This problem could perhaps be resolved by using Go potentials (Poma et al., 2017; Souza et al., 2019), which are already integrated into Martinize2 and Bentopy.
Even though our in silico cell contains more than 500 unique CG molecules and thereby presumably qualifies as the most complex system simulated to date, it simplifies the composition of various cellular components of Syn3A. Firstly, limited by the availability of Martini models for the metabolites, only a small subset is currently included in the cytosol. Future iterations of our whole-cell model will include Martini models for the complete metabolome, which comprises about 188 different compounds, and are expected to benefit from the ongoing development of dedicated automatic topology builders (; Potter et al., 2021). Secondly, since AlphaFold2 was used to predict the protein structures of the whole proteome, only monomeric structures were initially available. Essential multimeric proteins like the ribosomes and membrane-embedded transport complexes are either left out or represented by homologous proteins with available experimental crystal structures. In the future, improved protein structure prediction algorithms will be used that also facilitate the modeling of multimeric protein structures. In addition, ongoing progress in the experimental characterization of the Syn3A proteome and lipidome, as well as the characterization of the spatial distributions of membrane proteins, will help further increase our model’s realism. A “living” list of the complete composition of our in silico cell can be found in our GitHub repository (marrink-lab, 2022).
Another issue is the fine-tuning of the amount of interior solvent (both water and ions), together with the lipid balance between the inner and outer leaflet. Previous works on large-scale membrane-enveloped systems (Pezeshkian et al., 2021; Vermaas et al., 2022) have shown that finding this balance is a non-trivial task. Unbalanced systems might experience strong osmotic pressures and membrane (curvature) stress, causing unwanted shape deformations all the way to membrane rupture. As a complicating factor, these effects may only appear after prolonged simulation times. Clearly, dedicated computational resources are required for the simulation of whole cells or cell organelles. The forthcoming generation of supercomputers and simulation software is becoming increasingly efficient, and billion-particle simulations have already been achieved (; ).
An important challenge is reaching timescales long enough to allow meaningful analysis of such large systems. Assuming dedicated computer time on current infrastructure, we anticipate that we can reach timescales of the order of 10–100 µs in the foreseeable future. Although this is typically considered a long enough simulation time for standard system sizes (e.g., a single membrane protein), it is clear that on the scale of an entire cell, we will not be able to equilibrate our system; the generated ensemble of configurations will remain dependent on our starting state. Equilibration will only happen locally, and multiple replicas will need to be generated to obtain statistically relevant data. Note that the 10–100 µs range offers a nice overlap with state-of-the-art experimental techniques. For example, advanced MINFLUX microscopy from the Hell lab enables the tracking of particles as small as 1–2 nm for 100 s of microseconds (; Schmidt et al., 2021). Besides, Lattice Microbes simulations of the Luthey-Schulten group (Roberts et al., 2013) use time steps of the order of microseconds, which allows for a potential feedback loop between these computational approaches.
Another major limitation is the fact that real cells operate out-of-equilibrium, driven by the import and export of nutrients and an intricate metabolic network of chemical reactions. In our approach, which is based on classical MD, we do not take this into account. We are therefore limited to studying non-reactive processes, i.e., those arising from the physical interactions among the constituents. The current composition of our cell is based on average concentrations of proteins and metabolites and thus reflects a steady-state. Coupling our classic approach with approaches taking into account reactivity, such as the aforementioned Lattice Microbes simulations or other metabolic network models (see below), in principle, could capture the non-equilibrium aspect of real cells.
Keeping these limitations in mind, simulations of the minimal cell with a molecular resolution will make it possible to study a wide range of new aspects. Modelling cellular processes and chemical transformations involves a hierarchy of interconnected scales that cannot be separated without causing artefacts. Behaviour emerging from the interaction of millions of different compounds is easily missed when systems are simplified. One might question to what extent one part of the cell affects another, given the limited timescales likely to be reached. If the various cellular subsystems act independently, one might better simulate those in isolation. To find out, one needs to simulate the complete system in addition to the smaller-scale subsystems. Our whole-cell simulation is only a first step, which will benefit from imminent improvements in high-performance computing to extend these simulations to longer timescales, up to the point where all parts of the cell may influence each other. Currently, the internal organization of the cytosol of Syn3A is primarily a black box. Our model will allow us to observe how proteins inside the cytosol interact with macromolecular structures such as ribosomes and chromosomes. Viewing the cytosol from this perspective, we can observe emerging heterogeneities and viscosity gradients, following in the footsteps of other realistic models of the cytoplasm of various cell types (McGuffee and Elcock, 2010; Yu et al., 2016; Oliveira Bortot et al., 2020). We can expect arising interaction patterns between proteins and metabolites, and probe the possible appearance of biomolecular condensates (; Rhine et al., 2020).
A simulation at the level of the entire cell allows us to characterize the extent to which the cell membrane affects (and is affected by) the cellular interior. If we consider a membrane zone with a thickness of 30 nm (∼20 nm of the membrane together with its embedded proteins, plus another 10 nm layer underneath), 40% of the total cell volume is part of this membrane zone. Our simulations will provide detailed insights into the nature and extent of depletion or crowding layers, and into the level of heterogeneity inside this membrane zone, providing information on the extent to which compounds are either enriched or depleted near the cell surface (Nawrocki et al., 2019). A full-cell membrane model might explain why the minimal cell grows on a diet of both saturated and unsaturated fatty acids, but not on a diet of just saturated ones as observed in lipidomics experiments from the Saenz lab (private communication). A related question is why the cell membrane contains such a high percentage of cholesterol (20%–60% dependent on growth medium); this is uncommon for bacterial membranes although generally Mycoplasma do contain some cholesterol for membrane stability.
Of special interest is the potential existence of dynamic highways, i.e., regions in the cell with greater mobility of the constituents, which may arise from crowding effects or liquid-liquid phase separation phenomena, or may be induced by proximity of the cell membrane. Such dynamic highways could be important in regulating transport in an otherwise glassy state of the cytoplasm. For regions of the cell showing particularly interesting behaviour, smaller systems can be extracted with the advanced TS2CG tool and simulated for extended timescales to increase the statistical relevance. Besides passively studying the cellular environment, holistic cell modeling poses the ideal computational sandbox in which we can introduce new components to the cellular environment. For instance, elucidating the non-specific interactions between the cytosol and drug candidates and showing how drug-receptor interactions affect the entire cell instead of just the receptor site.
Using a multiscale modeling approach, we could potentially explore cell dynamics at various stages in its life cycle. Compared to MD simulations, other low-resolution modeling approaches can more broadly explore timescales of several orders of magnitude longer. Integrating other computational models will make it possible to sprout MD simulations in interesting regimes observed with the lower-resolution models. The primary computational method we will focus on integrating into our framework is the whole-cell fully dynamical kinetic model developed by the Luthey-Schulten lab, which accounts for the metabolic pathways governing the cellular processes (Thornburg et al., 2022). By transferring structural information from the kinetic model into our high-resolution model, it will be possible to paint a more detailed picture of the cell’s internal organization and dynamics at specific points of the cell’s life cycle, including during cell fission.
Since most of the tools in the Martini ecosystem are force field agnostic, the workflow can also be applied to generate all-atom whole-cell models. Given the substantial increase in associated computational costs, it might be a wiser approach to only sample smaller subsystems at the all-atom level. These could be straightforwardly obtained from backmapping representative regions taken from the whole-cell CG model. A number of such backmapping tools, optimized for Martini, already exist (; Vickery and Stansfeld, 2021; ).
A final challenge lies in the analysis and interpretation of the complex high-dimensional massive data that will be generated. Clearly, it will be impossible to perform a comprehensive analysis on a whole-cell trajectory, and one needs to focus on specific research questions. However, the trajectories can nowadays be easily shared with the broader community via dedicated open-access repositories such as Zenodo (https://zenodo.org/), allowing others to perform whatever additional analysis they fancy. One can also envision the usage of data reduction schemes to efficiently analyse the whole-cell simulation. One possibility is storing only centers-of-mass movement of the non-aqueous components, which would facilitate the analysis of diffusional behavior, for instance. Another approach would be using a voxel-based method () to dynamically segment the whole-cell model into similarity regions, e.g., membrane periphery or chromosomal region. The system segmentation would allow for efficient quantitative comparison of the cytosolic properties within and between distinct regions of the cell. Furthermore, machine-learning can be invoked to extract interaction patterns and other emergent behavior that might be missed by standard analysis tools (Noé et al., 2020; Wang et al., 2020; ). We foresee that our data sets will generate novel ways of dealing with this unprecedented level of complexity.
In conclusion, we presented a roadmap toward whole-cell MD simulations, illustrated with the construction of the first MD model of an entire cell using our Martini ecosystem. The model represents a next level realized with the computational microscope, providing a complete picture of the cell and making it possible to relate molecular structures and interactions to cellular function directly. In the long term, our computational framework will enable us to study a wide variety of mesoscopic systems, possibly informing the design of fully synthetic cells (Olivi et al., 2021) and modeling cells with more complex internal structures.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: https://github.com/marrink-lab/.
Author contributions
SM and ZL-S conceived of the project. JS and FG constructed the model with help of PT, MK, BG, TB, and ZT. JS, FG, and SM wrote the manuscript with input from all authors.
Funding
SM acknowledges funding from the ERC with the Advanced grant 101053661 (“COMP-O-CELL”), and funding from NWO through the NWA grant “The limits to growth: The challenge to dissipate energy” and BaSyc (“Building a Synthetic Cell”) consortium. ZL-S acknowledges funding from NSF (MCB: 2221237 “Simulating a growing minimal cell: Integrating experiment and theory”; PHY: 1430124 “Center for the Physics of Living Cells”; PHY: 1505008 and 2014027 “Collaborative Research Network from Physics of Living Systems”).
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
AbrahamM. J.MurtolaT.SchulzR.PallS.SmithJ. C.HessB.et al (2015). ‘GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers’. SoftwareX1–2, 19–25. 10.1016/j.softx.2015.06.001
2
AlessandriR.BarnoudJ.GertsenA. S.PatmanidisI.de VriesA. H.SouzaP. C. T.et al (2022). Martini 3 coarse-grained force field: Small molecules. Adv. Theory Simulations5 (1), 2100391. 10.1002/adts.202100391
3
AndoT.BhamidimarriS. P.BrendingN.Colin-YorkH.CollinsonL.De JongeN.et al (2018). The 2018 correlative microscopy techniques roadmap. Correl. Microsc. Tech. roadmap’, J. Phys. D Appl. Phys.51 (44), 443001. 10.1088/1361-6463/aad055
4
BereauT.KremerK. (2015). Automated parametrization of the coarse-grained Martini force field for small organic molecules. J. Chem. Theory Comput.11 (6), 2783–2791. 10.1021/acs.jctc.5b00056
5
BhatN. G.BalajiS. (2020). Whole-cell modeling and simulation: A brief survey. New Gener. Comput.38 (1), 259–281. 10.1007/s00354-019-00066-y
6
BianchiD. M.PelletierJ. F.HutchisonC. A.GlassJ. I.Luthey-SchultenZ. (2022). ‘Toward the complete functional characterization of a minimal bacterial proteome’. J. Phys. Chem. B126 (36), 6820–6834. 10.1021/acs.jpcb.2c04188
7
Blender Online Community (2022). ‘Blender - a 3D modelling and rendering package’. Amsterdam: Blender Foundation. Stichting Blender Foundation. Available at: http://www.blender.org
8
BonvinA. M. J. J. (2021). 50 years of PDB: A catalyst in structural biology. Nat. Methods18 (5), 448–449. 10.1038/s41592-021-01138-y
9
BrackleyC. A.MorozovA. N.MarenduzzoD. (2014). ‘Models for twistable elastic polymers in Brownian dynamics, and their implementation for LAMMPS’. J. Chem. Phys.140 (13), 135103. 10.1063/1.4870088
10
BreuerM.EarnestT. M.MerrymanC.WiseK. S.SunL.LynottM. R.et al (2019). ‘Essential metabolism for a minimal cell’. eLife8, e36842. 10.7554/elife.36842
11
BruininksB. M. H.ThieA. S.SouzaP. C. T.WassenaarT. A.FarajiS.MarrinkS. J. (2021). ‘Sequential voxel-based leaflet segmentation of complex lipid morphologies’. J. Chem. Theory Comput.17 (12), 7873–7885. 10.1021/acs.jctc.1c00446
12
CastagnaJ.GuoX.SeatonM.O’CaisA. (2020). Towards extreme scale dissipative particle dynamics simulations using multiple GPGPUs. Comput. Phys. Commun.251, 107159. 10.1016/j.cpc.2020.107159
13
ChengY. (2018). ‘Single-particle cryo-EM—how did it get here and where will it go’. Science361 (6405), 876–880. 10.1126/science.aat4346
14
ChorevD. S.BakerL. A.WuD.Beilsten-EdmandsV.RouseS. L.Zeev-Ben-MordehaiT.et al (2018). ‘Protein assemblies ejected directly from native membranes yield complexes for mass spectrometry’. Science362 (6416), 829–834. 10.1126/science.aau0976
15
ChristieS.ShiX.SmithA. W. (2020). ‘Resolving membrane protein–protein interactions in live cells with pulsed interleaved excitation fluorescence cross-correlation spectroscopy’. Accounts Chem. Res.53 (4), 792–799. 10.1021/acs.accounts.9b00625
16
CorradiV.Mendez-VilluendasE.IngolfssonH. I.GuR. X.SiudaI.MeloM. N.et al (2018). ‘Lipid–Protein interactions are unique fingerprints for membrane proteins’. ACS Central Sci.4 (6), 709–717. 10.1021/acscentsci.8b00143
17
de JongD. H.SinghG.BennettW. F. D.ArnarezC.WassenaarT. A.SchaferL. V.et al (2013). ‘Improved parameters for the Martini coarse-grained protein force field’. J. Chem. Theory Comput.9 (1), 687–697. 10.1021/ct300646g
18
DommerA.CasalinoL.KearnsF.RosenfeldM.WauerN.AhnS. H.et al (2022). ‘#COVIDisAirborne: AI-enabled multiscale computational microscopy of delta SARS-CoV-2 in a respiratory aerosol’. Int. J. High Perform. Comput. Appl., 109434202211282. 10.1177/10943420221128233
19
DrorR. O.DirksR. M.GrossmanJ.XuH.ShawD. E. (2012). Biomolecular simulation: A computational microscope for molecular biology. Annu. Rev. Biophysics41 (1), 429–452. 10.1146/annurev-biophys-042910-155245
20
EastmanP.SwailsJ.ChoderaJ. D.McGibbonR. T.ZhaoY.BeauchampK. A.et al (2017). ‘OpenMM 7: Rapid development of high performance algorithms for molecular dynamics’. PLOS Comput. Biol.13 (7), e1005659. 10.1371/journal.pcbi.1005659
21
EilersY.TaH.GwoschK. C.BalzarottiF.HellS. W. (2018). ‘MINFLUX monitors rapid molecular jumps with superior spatiotemporal resolution’. Proc. Natl. Acad. Sci.115 (24), 6117–6122. 10.1073/pnas.1801672115
22
EllisR. J.MintonA. P. (2003). ‘Join the crowd’. Nature425 (6953), 27–28. 10.1038/425027a
23
GilbertB. R.ThornburgZ. R.LamV.RashidF. Z. M.GlassJ. I.VillaE.et al (2021). ‘Generating chromosome geometries in a minimal cell from cryo-electron tomograms and chromosome conformation capture maps’. Front. Mol. Biosci.8, 644133. 10.3389/fmolb.2021.644133
24
GrünewaldF.AlessandriR.KroonP. C.MonticelliL.SouzaP. C. T.MarrinkS. J. (2022). Polyply; a python suite for facilitating simulations of macromolecules and nanomaterials. Nat. Commun.13 (1), 68. 10.1038/s41467-021-27627-4
25
GrünewaldF.PuntM. H.JefferysE. E.VainikkaP. A.KonigM.VirtanenV.et al (2022). ‘Martini 3 coarse-grained force field for carbohydrates’. J. Chem. Theory Comput.18, 7555–7569. 10.1021/acs.jctc.2c00757
26
GuilhasB.WalterJ. C.RechJ.DavidG.WalliserN. O.PalmeriJ.et al (2020). ‘ATP-Driven separation of liquid phase condensates in bacteria’. Mol. Cell79 (2), 293–303.e4. 10.1016/j.molcel.2020.06.034
27
GuptaC.SarkarD.TielemanD. P.SingharoyA. (2022). ‘The ugly, bad, and good stories of large-scale biomolecular simulations’. Curr. Opin. Struct. Biol.73, 102338. 10.1016/j.sbi.2022.102338
28
HilpertC.BerangerL.SouzaP. C.VainikkaP. A.NietoV.MarrinkS. J.et al (2023). Facilitating CG simulations with MAD: The MArtini Database server. J. Chem. Inf. Model. 10.1021/acs.jcim.2c01375
29
HowardM. P.AndersonJ. A.NikoubashmanA.GlotzerS. C.PanagiotopoulosA. Z. (2016). ‘Efficient neighbor list calculation for molecular simulation of colloidal systems using graphics processing units’. Comput. Phys. Commun.203, 45–52. 10.1016/j.cpc.2016.02.003
30
HutchisonC. A.ChuangR. Y.NoskovV. N.Assad-GarciaN.DeerinckT. J.EllismanM. H.et al (2016). ‘Design and synthesis of a minimal bacterial genome’. Science351, aad6253. 10.1126/science.aad6253
31
JohnsonG. T.AutinL.Al-AlusiM.GoodsellD. S.SannerM. F.OlsonA. J. (2015). cellPACK: A virtual mesoscope to model and visualize structural systems biology. Nat. Methods12 (1), 85–91. 10.1038/nmeth.3204
32
JumperJ.EvansR.PritzelA.GreenT.FigurnovM.RonnebergerO.et al (2021). ‘Highly accurate protein structure prediction with AlphaFold’. Nature596 (7873), 583–589. 10.1038/s41586-021-03819-2
33
JungJ.NishimaW.DanielsM.BascomG.KobayashiC.AdedoyinA.et al (2019). ‘Scaling molecular dynamics beyond 100, 000 processor cores for large-scale biophysical simulations’. J. Comput. Chem.40 (21), 1919–1930. 10.1002/jcc.25840
34
KaptanS.VattulainenI. (2022). ‘Machine learning in the analysis of biomolecular simulations’. Adv. Phys. X7 (1), 2006080. 10.1080/23746149.2021.200608
35
KarrJ. R.SanghviJ.MacklinD.GutschowM.JacobsJ.BolivalB.et al (2012). ‘A whole-cell computational model predicts phenotype from genotype’. Cell150 (2), 389–401. 10.1016/j.cell.2012.05.044
36
KarrJ. R.TakahashiK.FunahashiA. (2015). ‘The principles of whole-cell modeling’. Curr. Opin. Microbiol.27, 18–24. 10.1016/j.mib.2015.06.004
37
KhalidS.RouseS. L. (2020). ‘Simulation of subcellular structures’. Curr. Opin. Struct. Biol.61, 167–172. 10.1016/j.sbi.2019.12.017
38
KroonP. C.GrünewaldF.BarnoudJ.van TilburgM.SouzaP. C. T.MarrinkS. J. (2022). ‘Martinize2 and vermouth: Unified framework for topology generation’. arXiv. 10.48550/arXiv.2212.01191
39
LeeE. H.HsinJ.SotomayorM.ComellasG.SchultenK. (2009). ‘Discovery through the computational microscope’. Structure17 (10), 1295–1306. 10.1016/j.str.2009.09.001
40
LópezC. A.RzepielaA. J.de VriesA. H.DijkhuizenL.HunenbergerP. H.MarrinkS. J. (2009). Martini coarse-grained force field: Extension to carbohydrates. J. Chem. Theory Comput.5 (12), 3195–3210. 10.1021/ct900313w
41
LópezC. A.ZhangX.AydinF.ShresthaR.VanQ. N.StanleyC. B.et al (2022). ‘Asynchronous reciprocal coupling of Martini 2.2 coarse-grained and CHARMM36 all-atom simulations in an automated multiscale framework’. J. Chem. Theory Comput.18 (8), 5025–5045. 10.1021/acs.jctc.2c00168
42
LorentJ. H.LeventalK. R.GanesanL.Rivera-LongsworthG.SezginE.DoktorovaM.et al (2020). ‘Plasma membranes are asymmetric in lipid unsaturation, packing and protein shape’. Nat. Chem. Biol.16 (6), 644–652. 10.1038/s41589-020-0529-6
43
LouisonK. A.DrydenI. L.LaughtonC. A. (2021). Glimps: A machine learning approach to resolution transformation for multiscale modeling. J. Chem. Theory Comput.17 (12), 7930–7937. 10.1021/acs.jctc.1c00735
44
Luthey-SchultenZ. (2021). ‘Integrating experiments, theory and simulations into whole-cell models’. Nat. Methods18 (5), 446–447. 10.1038/s41592-021-01150-2
45
Luthey-SchultenZ.ThornburgZ. R.GilbertB. R. (2022). ‘Integrating cellular and molecular structures and dynamics into whole-cell models’. Curr. Opin. Struct. Biol.75, 102392. 10.1016/j.sbi.2022.102392
46
MacklinD. N.Ahn-HorstT. A.ChoiH.RuggeroN. A.CarreraJ.MasonJ. C.et al (2020). Simultaneous cross-evaluation of heterogeneous E. coli datasets via mechanistic simulation. Science369, eaav3751. 10.1126/science.aav3751
47
MacklinD. N.RuggeroN. A.CovertM. W. (2014). ‘The future of whole-cell modeling’. Curr. Opin. Biotechnol.28, 111–115. 10.1016/j.copbio.2014.01.012
48
MaritanM.AutinL.KarrJ.CovertM. W.OlsonA. J.GoodsellD. S. (2022). ‘Building structural models of a whole Mycoplasma cell’. J. Mol. Biol.434 (2), 167351. 10.1016/j.jmb.2021.167351
49
MarrinkS. J.MonticelliL.MeloM. N.AlessandriR.TielemanD. P.SouzaP. C. T.et al (2022). Two decades of Martini: Better beads, broader scope. WIREs Comput. Mol. Sci., e1620. 10.1002/wcms.1620
50
MarrinkS. J.CorradiV.SouzaP. C.IngolfssonH. I.TielemanD. P.SansomM. S. (2019). ‘Computational modeling of realistic cell membranes’. Chem. Rev.119 (9), 6184–6226. 10.1021/acs.chemrev.8b00460
51
marrink-lab (2022). ‘Martini_Minimal_Cell’. Available at https://github.com/marrink-lab/Martini_Minimal_Cell.
52
McGuffeeS. R.ElcockA. H. (2010). Diffusion, crowding and protein stability in a dynamic molecular model of the bacterial cytoplasm. PLoS Comput. Biol.6 (3), e1000694. 10.1371/journal.pcbi.1000694
53
MondalJ.BrattonB. P.LiY.YethirajA.WeisshaarJ. (2011). ‘Entropy-Based mechanism of ribosome-nucleoid segregation in E. coli cells’. Biophysical J.100 (11), 2605–2613. 10.1016/j.bpj.2011.04.030
54
MosalagantiS.Obarska-KosinskaA.SiggelM.TaniguchiR.TuronovaB.ZimmerliC. E.et al (2022). ‘AI-based structure prediction empowers integrative structural analysis of human nuclear pores’. Science376, eabm9506. 10.1126/science.abm9506
55
NarasimhanS.FolkersG. E.BaldusM. (2020). When small becomes too big: Expanding the use of in-cell solid-state NMR spectroscopy. ChemPlusChem85 (4), 760–768. 10.1002/cplu.202000167
56
NawrockiG.ImW.SugitaY.FeigM. (2019). ‘Clustering and dynamics of crowded proteins near membranes and their influence on membrane bending’. Proc. Natl. Acad. Sci.116 (49), 24562–24567. 10.1073/pnas.1910771116
57
NoéF.TkatchenkoA.MullerK. R.ClementiC. (2020). ‘Machine learning for molecular simulation’. Annu. Rev. Phys. Chem.71 (1), 361–390. 10.1146/annurev-physchem-042018-052331
58
Oliveira BortotL.BashardaneshZ.van der SpoelD. (2020). Making soup: Preparing and validating models of the bacterial cytoplasm for molecular simulation. J. Chem. Inf. Model.60 (1), 322–331. 10.1021/acs.jcim.9b00971
59
OliviL.BergerM.CreyghtonR. N. P.De FranceschiN.DekkerC.MulderB. M.et al (2021). ‘Towards a synthetic cell cycle’. Nat. Commun.12 (1), 4531. 10.1038/s41467-021-24772-8
60
PezeshkianW.GrunewaldF.NarykovO.LuS.ArkhipovaV.SolodovnikovA.et al (2021). Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling. Biophysics. 10.1101/2021.09.15.459697
61
PezeshkianW.KonigM.WassenaarT. A.MarrinkS. J. (2020). ‘Backmapping triangulated surfaces to coarse-grained membrane models’. Nat. Commun.11 (1), 2296. 10.1038/s41467-020-16094-y
62
PomaA. B.CieplakM.TheodorakisP. E. (2017). ‘Combining the MARTINI and structure-based coarse-grained approaches for the molecular dynamics studies of conformational transitions in proteins’. J. Chem. Theory Comput.13 (3), 1366–1374. 10.1021/acs.jctc.6b00986
63
PotterT. D.BarrettE. L.MillerM. A. (2021). ‘Automated coarse-grained mapping algorithm for the Martini force field and benchmarks for membrane–water partitioning’. J. Chem. Theory Comput.17 (9), 5777–5791. 10.1021/acs.jctc.1c00322
64
ReadingE.HallZ.MartensC.HaghighiT.FindlayH.AhdashZ.et al (2017). ‘Interrogating membrane protein conformational dynamics within native lipid compositions’. Angew. Chem. Int. Ed.56 (49), 15654–15657. 10.1002/anie.201709657
65
RhineK.VidaurreV.MyongS. (2020). ‘RNA droplets’. Annu. Rev. Biophysics49 (1), 247–265. 10.1146/annurev-biophys-052118-115508
66
RobertsE.StoneJ. E.Luthey-SchultenZ. (2013). Lattice microbes: High-performance stochastic simulation method for the reaction-diffusion master equation. J. Comput. Chem.34 (3), 245–255. 10.1002/jcc.23130
67
SchmidtR.WeihsT.WurmC. A.JansenI.RehmanJ.SahlS. J.et al (2021). ‘MINFLUX nanometer-scale 3D imaging and microsecond-range tracking on a common fluorescence microscope’. Nat. Commun.12 (1), 1478. 10.1038/s41467-021-21652-z
68
SingharoyA.MaffeoC.Delgado-MagneroK. H.SwainsburyD. J.SenerM.KleinekathoferU.et al (2019). Atoms to phenotypes: Molecular design principles of cellular energy metabolism. Cell179, 1098–1111.e23. 10.1016/j.cell.2019.10.021
69
SousaF. M.LimaL. M. P.ArnarezC.PereiraM. M.MeloM. N. (2021). ‘Coarse-Grained parameterization of nucleotide cofactors and metabolites: Coarse-grained parameterization of nucleotide cofactors and metabolites: Protonation constants, partition coefficients, and model topologies. J. Chem. Inf. Model.61 (1), 335–346. 10.1021/acs.jcim.0c01077
70
SouzaP. C. T.AlessandriR.BarnoudJ.ThallmairS.FaustinoI.GrunewaldF.et al (2021). Martini 3: A general purpose force field for coarse-grained molecular dynamics. Nat. Methods18 (4), 382–388. 10.1038/s41592-021-01098-3
71
SouzaP. C. T.ThallmairS.MarrinkS. J.Mera-AdasmeR. (2019). ‘An allosteric pathway in copper, an allosteric pathway in copper, zinc superoxide dismutase unravels the molecular mechanism of the G93A amyotrophic lateral sclerosis-linked mutation. J. Phys. Chem. Lett.10 (24), 7740–7744. 10.1021/acs.jpclett.9b02868
72
ŠteflM.HerbstK.RubsamM.BendaA.KnopM. (2020). ‘Single-Color fluorescence lifetime cross-correlation spectroscopy in vivo’. Biophysical J.119 (7), 1359–1370. 10.1016/j.bpj.2020.06.039
73
ThornburgZ. R.BianchiD. M.BrierT. A.GilbertB. R.EarnestT. M.MeloM. C.et al (2022). ‘Fundamental behaviors emerge from simulations of a living minimal cell’. Cell185 (2), 345–360.e28. 10.1016/j.cell.2021.12.025
74
UusitaloJ. J.IngolfssonH. I.AkhshiP.TielemanD. P.MarrinkS. J. (2015). Martini coarse-grained force field: Extension to DNA. J. Chem. Theory Comput.11 (8), 3932–3945. 10.1021/acs.jctc.5b00286_FILE/CT5B00286_SI_001.PDF
75
UusitaloJ. J.IngolfssonH. I.MarrinkS. J.FaustinoI. (2017). Martini coarse-grained force field: Extension to RNA. Biophysical J.113 (2), 246–256. 10.1016/j.bpj.2017.05.043
76
VermaasJ. V.MayneC. G.ShinnE.TajkhorshidE. (2022). ‘Assembly and analysis of cell-scale membrane envelopes’. J. Chem. Inf. Model.62 (3), 602–617. 10.1021/acs.jcim.1c01050
77
VickeryO. N.StansfeldP. J. (2021). CG2AT2: An enhanced fragment-based approach for serial multi-scale molecular dynamics simulations. J. Chem. Theory Comput.17 (10), 6472–6482. 10.1021/acs.jctc.1c00295
78
WangW.JuttlerB.ZhengD.LiuY. (2008). ‘Computation of rotation minimizing frames’. ACM Trans. Graph. (TOG)18 (1), 1. 10.1145/1330511.1330513
79
WangY.Lamim RibeiroJ. M.TiwaryP. (2020). ‘Machine learning approaches for analyzing and enhancing molecular dynamics simulations’. Curr. Opin. Struct. Biol.61, 139–145. 10.1016/j.sbi.2019.12.016
80
WassenaarT. A.IngolfssonH. I.BockmannR. A.TielemanD. P.MarrinkS. J. (2015). Computational lipidomics with insane: A versatile tool for generating custom membranes for molecular simulations. J. Chem. Theory Comput.11 (5), 2144–2155. 10.1021/acs.jctc.5b00209
81
WietrzynskiW.SchafferM.TegunovD.AlbertS.KanazawaA.PlitzkoJ. M.et al (2020). ‘Charting the native architecture of Chlamydomonas thylakoid membranes with single-molecule precision’. eLife9, e53740. 10.7554/eLife.53740
82
YuI.MoriT.AndoT.HaradaR.JungJ.SugitaY.et al (2016). ‘Biomolecular interactions modulate macromolecular structure and dynamics in atomistic model of a bacterial cytoplasm’. eLife5, e19274. 10.7554/ELIFE.19274
83
ZhangX.SundramS.OppelstrupT.Kokkila-SchumacherS. I. L.CarpenterT. S.IngolfssonH. I.et al (2020). ddcMD: A fully GPU-accelerated molecular dynamics program for the Martini force field. J. Chem. Phys.153 (4), 045103. 10.1063/5.0014500
Summary
Keywords
JCVI-syn3A, minimal cell, Martini force field, integrative modeling, coarse grain, polyply
Citation
Stevens JA, Grünewald F, van Tilburg PAM, König M, Gilbert BR, Brier TA, Thornburg ZR, Luthey-Schulten Z and Marrink SJ (2023) Molecular dynamics simulation of an entire cell. Front. Chem. 11:1106495. doi: 10.3389/fchem.2023.1106495
Received
23 November 2022
Accepted
09 January 2023
Published
18 January 2023
Volume
11 - 2023
Edited by
Adolfo Poma, Polish Academy of Sciences, Poland
Reviewed by
Rosana Collepardo, University of Cambridge, United Kingdom
Roberto Covino, Frankfurt Institute for Advanced Studies, Germany
Updates
Copyright
© 2023 Stevens, Grünewald, van Tilburg, König, Gilbert, Brier, Thornburg, Luthey-Schulten and Marrink.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Siewert J. Marrink, s.j.marrink@rug.nl
This article was submitted to Theoretical and Computational Chemistry, a section of the journal Frontiers in Chemistry
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.