A Coarse-Grained SPICA Makeover for Solvated and Bare Sodium and Chloride Ions.

Aqueous ionic solutions are pivotal in various scientific domains due to their natural prevalence and vital roles in biological and chemical processes. Molecular dynamics has emerged as an effective methodology for studying the dynamic behavior of these systems. While all-atomistic models have made significant strides in accurately representing and simulating these ions, the challenge persists in achieving precise models for coarse-grained (CG) simulations. Our study introduces two optimized models for sodium and chloride ions within the nonpolarizable surface property fitting coarse-grained force field (SPICA-FF) framework. The two models represent solvated ions, such as the original FF model, and unsolvated or bare ions. The nonbonded Lennard-Jones interactions were reparameterized to faithfully reproduce bulk properties, including density and surface tension, in sodium chloride solutions at varying concentrations. Notably, these optimized models replicate experimental surface tensions at high ionic strengths, a property not well-captured by the ions of the original model in the SPICA-FF. The optimized unsolvated model also proved successful in reproducing experimental osmotic pressure. Additionally, the newly reparameterized ion models capture hydrophobic interactions within sodium chloride solutions and show qualitative agreement when modeling structural changes in phospholipid bilayers, aligning with experimental observations. For aqueous solutions, these optimized models promise a more precise representation of the ion behavior.

Cryo-EM architecture of a near-native stretch-sensitive membrane microdomain.

Biological membranes are partitioned into functional zones termed membrane microdomains, which contain specific lipids and proteins1-3. The composition and organization of membrane microdomains remain controversial because few techniques are available that allow the visualization of lipids in situ without disrupting their native behaviour3,4. The yeast eisosome, composed of the BAR-domain proteins Pil1 and Lsp1 (hereafter, Pil1/Lsp1), scaffolds a membrane compartment that senses and responds to mechanical stress by flattening and releasing sequestered factors5-9. Here we isolated near-native eisosomes as helical tubules made up of a lattice of Pil1/Lsp1 bound to plasma membrane lipids, and solved their structures by helical reconstruction. Our structures reveal a striking organization of membrane lipids, and, using in vitro reconstitutions and molecular dynamics simulations, we confirmed the positioning of individual PI(4,5)P2, phosphatidylserine and sterol molecules sequestered beneath the Pil1/Lsp1 coat. Three-dimensional variability analysis of the native-source eisosomes revealed a dynamic stretching of the Pil1/Lsp1 lattice that affects the sequestration of these lipids. Collectively, our results support a mechanism in which stretching of the Pil1/Lsp1 lattice liberates lipids that would otherwise be anchored by the Pil1/Lsp1 coat, and thus provide mechanistic insight into how eisosome BAR-domain proteins create a mechanosensitive membrane microdomain.

In situ architecture of the ER-mitochondria encounter structure.

The endoplasmic reticulum and mitochondria are main hubs of eukaryotic membrane biogenesis that rely on lipid exchange via membrane contact sites1-3, but the underpinning mechanisms remain poorly understood. In yeast, tethering and lipid transfer between the two organelles is mediated by the endoplasmic reticulum-mitochondria encounter structure (ERMES), a four-subunit complex of unresolved stoichiometry and architecture4-6. Here we determined the molecular organization of ERMES within Saccharomyces cerevisiae cells using integrative structural biology by combining quantitative live imaging, cryo-correlative microscopy, subtomogram averaging and molecular modelling. We found that ERMES assembles into approximately 25 discrete bridge-like complexes distributed irregularly across a contact site. Each bridge consists of three synaptotagmin-like mitochondrial lipid binding protein domains oriented in a zig-zag arrangement. Our molecular model of ERMES reveals a pathway for lipids. These findings resolve the in situ supramolecular architecture of a major inter-organelle lipid transfer machinery and provide a basis for the mechanistic understanding of lipid fluxes in eukaryotic cells.

Development of a coarse-grained model for surface-functionalized gold nanoparticles: towards an accurate description of their aggregation behavior.

Understanding the dispersion stability and aggregation propensity of self-assembled monolayer gold NPs at a molecular level is crucial to guide their rational design and to inform about the optimal surface functionalization for specific applications. To reach this goal, in silico modeling via coarse-grained (CG) molecular dynamics (MD) simulations is a fundamental tool to complement the information acquired from experimental studies since CG modeling allows to get a deep knowledge of the molecular interactions that take place at the nanoscale in this kind of systems. Unfortunately, current CG models of monolayer-protected AuNPs present several drawbacks that limit their accuracy in certain scenarios. We here develop a CG model that is fully compatible and extends the SPICA/SDK (Shinoda-DeVane-Klein) force field. Our model allows reproducing the behavior of AuNPs functionalized with hydrophobic as well as charged and more hydrophilic ligands. This model improves upon results obtained with previously derived CG force fields and successfully describes NPs aggregation and self-assembly in aqueous solution.

Monovalent ion-mediated charge-charge interactions drive aggregation of surface-functionalized gold nanoparticles.

Monolayer-protected metal nanoparticles (NPs) are not only promising materials with a wide range of potential industrial and biological applications, but they are also a powerful tool to investigate the behaviour of matter at nanoscopic scales, including the stability of dispersions and colloidal systems. This stability is dependent on a delicate balance between attractive and repulsive interactions that occur in the solution, and it is described in quantitative terms by the classic Derjaguin-Landau-Vewey-Overbeek (DLVO) theory, that posits that aggregation between NPs is driven by van der Waals interactions and opposed by electrostatic interactions. To investigate the limits of this theory at the nanoscale, where the continuum assumptions required by the DLVO theory break down, here we investigate NP dimerization by computing the Potential of Mean Force (PMF) of this process using fully atomistic MD simulations. Serendipitously, we find that electrostatic interactions can lead to the formation of metastable NP dimers at physiological ion concentrations. These dimers are stabilized by complexes formed by negatively charged ligands belonging to distinct NPs that are bridged by positively charged monovalent ions present in solution. We validate our findings by collecting tomographic EM images of NPs in solution and by quantifying their radial distribution function, that shows a marked peak at interparticle distance comparable with that of MD simulations. Taken together, our results suggest that not only van der Waals interactions, but also electrostatic interactions mediated by monovalent ions at physiological concentrations, contribute to attraction between nano-sized charged objects at very short length scales.

Triglyceride lipolysis triggers liquid crystalline phases in lipid droplets and alters the LD proteome.

Lipid droplets (LDs) are reservoirs for triglycerides (TGs) and sterol-esters (SEs), but how these lipids are organized within LDs and influence their proteome remain unclear. Using in situ cryo-electron tomography, we show that glucose restriction triggers lipid phase transitions within LDs generating liquid crystalline lattices inside them. Mechanistically this requires TG lipolysis, which decreases the LD's TG:SE ratio, promoting SE transition to a liquid crystalline phase. Molecular dynamics simulations reveal TG depletion promotes spontaneous TG and SE demixing in LDs, additionally altering the lipid packing of the PL monolayer surface. Fluorescence imaging and proteomics further reveal that liquid crystalline phases are associated with selective remodeling of the LD proteome. Some canonical LD proteins, including Erg6, relocalize to the ER network, whereas others remain LD-associated. Model peptide LiveDrop also redistributes from LDs to the ER, suggesting liquid crystalline phases influence ER-LD interorganelle transport. Our data suggests glucose restriction drives TG mobilization, which alters the phase properties of LD lipids and selectively remodels the LD proteome.

Relative performance of Bayesian morphological clock and parsimony methods for phylogenetic reconstructions: Insights from the case of Myomiminae and Dryomyinae glirid rodents.

Extinct organisms provide crucial information about the origin and time of origination of extant groups. The importance of morphological phylogenetics for rigorously dating the tree of life is now widely recognized and has been revitalized by methodological developments such as the application of tip-dating Bayesian approaches. Traditionally, molecular clocks have been node calibrated. However, node calibrations are often unsatisfactory because they do not allow the fossil age to inform about phylogenetic hypothesis. The introduction of tip calibrations allow fossil species to be included alongside their living relatives, and the absence of molecular sequence data for these taxa to be remedied by supplementing the sequence alignments for living taxa with phenotype character matrices for both living and fossil taxa. Therefore, only phylogenetic analyses that take into account morphological characters can incorporate both fossil and extant species. Herein we present an unprecedented morphological dataset for a vast group of glirid rodents, to which different phylogenetic methodologies have been applied. We have compared the tree topologies resulting from traditional parsimony and Bayesian phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods seem to outperform parsimony with our dataset. The strict consensus tree recovered by tip dating invalidates the classic classification and allows dates to be proposed for the divergence and origin of the different clades.

Integrative Phylogenetics: Tools for Palaeontologists to Explore the Tree of Life.

The modern era of analytical and quantitative palaeobiology has only just begun, integrating methods such as morphological and molecular phylogenetics and divergence time estimation, as well as phenotypic and molecular rates of evolution. Calibrating the tree of life to geological time is at the nexus of many disparate disciplines, from palaeontology to molecular systematics and from geochronology to comparative genomics. Creating an evolutionary time scale of the major events that shaped biodiversity is key to all of these fields and draws from each of them. Different methodological approaches and data employed in various disciplines have traditionally made collaborative research efforts difficult among these disciplines. However, the development of new methods is bridging the historical gap between fields, providing a holistic perspective on organismal evolutionary history, integrating all of the available evidence from living and fossil species. Because phylogenies with only extant taxa do not contain enough information to either calibrate the tree of life or fully infer macroevolutionary dynamics, phylogenies should preferably include both extant and extinct taxa, which can only be achieved through the inclusion of phenotypic data. This integrative phylogenetic approach provides ample and novel opportunities for evolutionary biologists to benefit from palaeontological data to help establish an evolutionary time scale and to test core macroevolutionary hypotheses about the drivers of biological diversification across various dimensions of organisms.

Bayesian Morphological Clock versus Parsimony: An Insight into the Relationships and Dispersal Events of Postvacuum Cricetidae (Rodentia, Mammalia).

Establishing an evolutionary timeline is fundamental for tackling a great variety of topics in evolutionary biology, including the reconstruction of patterns of historical biogeography, coevolution, and diversification. However, the tree of life is pruned by extinction and molecular data cannot be gathered for extinct lineages. Until recently methodological challenges have prevented the application of tip-dating Bayesian approaches in morphology-based fossil-only data sets. Herein, we present a morphological data set for a group of cricetid rodents to which we apply an array of methods fairly new in paleontology that can be used by paleontologists for the analysis of entirely extinct clades. We compare the tree topologies obtained by traditional parsimony, time-calibrated, and noncalibrated Bayesian inference phylogenetic approaches and calculate stratigraphic congruence indices for each. Bayesian tip-dated clock methods outperform parsimony in the case of our data set, which includes highly homoplastic morphological characters. Regardless, all three topologies support the monophyly of Megacricetodontinae, Democricetodontinae, and Cricetodontinae. Dispersal and speciation events inferred through Bayesian Binary Markov chain Monte Carlo and biodiversity analyses provide evidence for a correlation between biogeographic events, climatic changes, and diversification in cricetids. [Bayesian tip-dating; Cricetidae; Miocene; morphological clock; paleobiodiversity; paleobiogeography; paleoecology; parsimony; STRAP.].

Recharging your fats: CHARMM36 parameters for neutral lipids triacylglycerol and diacylglycerol.

Neutral lipids (NLs) are an abundant class of cellular lipids. They are characterized by the total lack of charged chemical groups in their structure, and, as a consequence, they play a major role in intracellular lipid storage. NLs that carry a glycerol backbone, such as triacylglycerols (TGs) and diacylglycerols (DGs), are also involved in the biosynthetic pathway of cellular phospholipids, and they have recently been the subject of numerous structural investigations by means of atomistic molecular dynamics simulations. However, conflicting results on the physicochemical behavior of NLs were observed depending on the nature of the atomistic force field used. Here, we show that current phospholipid-derived CHARMM36 parameters for DGs and TGs cannot adequately reproduce interfacial properties of these NLs because of excessive hydrophilicity at the glycerol-ester region. By following a CHARMM36-consistent parameterization strategy, we develop improved parameters for both TGs and DGs that are compatible with both cutoff-based and particle mesh Ewald schemes for the treatment of Lennard-Jones interactions. We show that our improved parameters can reproduce interfacial properties of NLs and their behavior in more complex lipid assemblies. We discuss the implications of our findings in the context of intracellular lipid storage and NLs' cellular activity.

Investigating the structural properties of hydrophobic solvent-rich lipid bilayers.

In vitro reconstitutions of lipid membranes have proven to be an indispensable tool to rationalize their molecular complexity and to understand their role in countless cellular processes. However, amongst the various techniques used to reconstitute lipid bilayers in vitro, several approaches are not solvent-free, but rather contain residual hydrophobic solvents in between the two bilayer leaflets, generally as a consequence of the procedure used to generate the bilayer. To what extent the presence of these hydrophobic solvents modifies bilayer properties with respect to native, solvent-free, conditions remains an open question that has important implications for the appropriate interpretation of numerous experimental observations. Here, we thorouhgly characterize hydrophobic solvent-rich lipid bilayers using atomistic molecular dynamics simulations. Our data indicate that while the presence of hydrophobic solvents at high concentrations, such as hexadecane, has a significant effect on membrane thickness, their effects on surface properties, membrane order and lateral stress are quite moderate. Our results corroborate the validity of in vitro approaches as model systems for the investigations of biological membranes but raise a few cautionary aspects that must be considered when investigating specific membrane properties.

Protonation Equilibrium in the Active Site of the Photoactive Yellow Protein.

The role and existence of low-barrier hydrogen bonds (LBHBs) in enzymatic and protein activity has been largely debated. An interesting case is that of the photoactive yellow protein (PYP). In this protein, two short HBs adjacent to the chromophore, p-coumaric acid (pCA), have been identified by X-ray and neutron diffraction experiments. However, there is a lack of agreement on the chemical nature of these H-bond interactions. Additionally, no consensus has been reached on the presence of LBHBs in the active site of the protein, despite various experimental and theoretical studies having been carried out to investigate this issue. In this work, we perform a computational study that combines classical and density functional theory (DFT)-based quantum mechanical/molecular mechanical (QM/MM) simulations to shed light onto this controversy. Furthermore, we aim to deepen our understanding of the chemical nature and dynamics of the protons involved in the two short hydrogen bonds that, in the dark state of PYP, connect pCA with the two binding pocket residues (E46 and Y42). Our results support the existence of a strong LBHB between pCA and E46, with the H fully delocalized and shared between both the carboxylic oxygen of E46 and the phenolic oxygen of pCA. Additionally, our findings suggest that the pCA interaction with Y42 can be suitably described as a typical short ionic H-bond of moderate strength that is fully localized on the phenolic oxygen of Y42.

Pre-existing bilayer stresses modulate triglyceride accumulation in the ER versus lipid droplets.

Cells store energy in the form of neutral lipids (NLs) packaged into micrometer-sized organelles named lipid droplets (LDs). These structures emerge from the endoplasmic reticulum (ER) at sites marked by the protein seipin, but the mechanisms regulating their biogenesis remain poorly understood. Using a combination of molecular simulations, yeast genetics, and fluorescence microscopy, we show that interactions between lipids' acyl-chains modulate the propensity of NLs to be stored in LDs, in turn preventing or promoting their accumulation in the ER membrane. Our data suggest that diacylglycerol, which is enriched at sites of LD formation, promotes the packaging of NLs into LDs, together with ER-abundant lipids, such as phosphatidylethanolamine. On the opposite end, short and saturated acyl-chains antagonize fat storage in LDs and promote accumulation of NLs in the ER. Our results provide a new conceptual understanding of LD biogenesis in the context of ER homeostasis and function.

Accurate Estimation of Membrane Capacitance from Atomistic Molecular Dynamics Simulations of Zwitterionic Lipid Bilayers.

Lipid membranes are indispensable to life, and they regulate countless cellular processes. To investigate the properties of membranes under controlled conditions, numerous reconstitution methods have been developed over the last few decades. Several of these methods result in the formation of lipid bilayers containing residual hydrophobic molecules between the two monolayers. These contaminants might alter membrane properties, including bilayer thickness, that is usually inferred from measurements of membrane capacitance assuming a simple slab model. However, recent measurements on solvent-free bilayers raised significant questions on the reliability of this approach. To reconcile the observed discrepancies, we developed a protocol to predict membrane capacitance from the dielectric profile of lipid bilayers computed from molecular dynamics simulations. Our methodology shows excellent agreement against available data on solvent-free noncharged bilayers, and it confirms that the uniform slab model is a reliable approximation from which to infer membrane capacitance. We find that the effective electrical thickness contributing to membrane capacitance is different from the hydrophobic thickness inferred from X-ray scattering form factors. We apply our model to estimate the concentration of residual solvent in reconstituted systems, and we propose that our protocol could be used to infer membrane properties in the presence of hydrophobic solvents.

New insights on the role of ROS in the mechanisms of sonoporation-mediated gene delivery.

Reactive oxygen species (ROS) are hypothesized to play a role in the sonoporation mechanisms. Nevertheless, the acoustical phenomenon behind the ROS production as well as the exact mechanisms of ROS action involved in the increased cell membrane permeability are still not fully understood. Therefore, we investigated the key processes occurring at the molecular level in and around microbubbles subjected to ultrasound using computational chemistry methods. To confirm the molecular simulation predictions, we measured the ROS production by exposing SonoVue® microbubbles (MBs) to ultrasound using biological assays. To investigate the role of ROS in cell membrane permeabilization, cells were subjected to ultrasound in presence of MBs and plasmid encoding reporter gene, and the transfection level was assessed using flow cytometry. The molecular simulations showed that under sonoporation conditions, ROS can form inside the MBs. These radicals could easily diffuse through the MB shell toward the surrounding aqueous phase and participate in the permeabilization of nearby cell membranes. Experimental data confirmed that MBs favor spontaneous formation of a host of free radicals where HO was the main ROS species after US exposure. The presence of ROS scavengers/inhibitors during the sonoporation process decreased both the production of ROS and the subsequent transfection level without significant loss of cell viability. In conclusion, the exposure of MBs to ultrasound might be the origin of chemical effects, which play a role in the cell membrane permeabilization and in the in vitro gene delivery when generated in its proximity.

First levantine fossil murines shed new light on the earliest intercontinental dispersal of mice.

Recent extensive field prospecting conducted in the Upper Miocene of Lebanon resulted in the discovery of several new fossiliferous localities. One of these, situated in the Zahleh area (Bekaa Valley, central Lebanon) has yielded a particularly diverse vertebrate fauna. Micromammals constitute an important part of this assemblage because not only do they represent the first Neogene rodents and insectivores from Lebanon, but they are also the only ones from the early Late Miocene of the Arabian Peninsula and circumambient areas. Analyses of the murines from Zahleh reveal that they belong to a small-sized early Progonomys, which cannot be assigned to any of the species of the genus hitherto described. They are, thereby, shown to represent a new species: Progonomys manolo. Morphometric analyses of the outline of the first upper molars of this species suggest a generalist and omnivorous diet. This record sheds new light onto a major phenomenon in the evolutionary history of rodents, which is the earliest dispersal of mice. It suggests that the arrival of murines in Africa got under way through the Levant rather than via southern Europe and was monitored by the ecological requirements of Progonomys.

An Atomistic Look into Bio-inspired Nanoparticles and their Molecular Interactions with Cells.

Nanoparticles (NPs) have sizes that approach those of pathogens and they can interact with the membranes of eukaryotic cells in an analogous fashion. Typically, NPs are taken up by the cell via the plasma membrane by receptor-mediated processes and subsequently interact with various endomembranes. Unlike pathogens, however, NPs lack the remarkable specificity gained during the evolutionary process and their design and optimization remains an expensive and time-consuming undertaking, especially considering the limited information available on their molecular interactions with cells. In this context, molecular dynamics (MD) simulations have emered as a promising strategy to investigate the mechanistic details of the interaction of NPs with mammalian or viral membranes. In particular, MD simulations have been extensively used to study the uptake process of NPs into the cell, focusing on membrane vesiculation, endocytic routes, or passive permeation processes. While such work is certainly relevant for understanding NP-cell interactions, it remains very difficult to determine the correspondence between generic models and the actual NP. Here, we review how chemically-specific MD simulations can provide rational guidelines towards further bio-inspired NP optimization.

The Molecular Mechanism of the Catalase-like Activity in Horseradish Peroxidase.

Horseradish peroxidase (HRP) is one of the most relevant peroxidase enzymes, used extensively in immunochemistry and biocatalysis applications. Unlike the closely related catalase enzymes, it exhibits a low activity to disproportionate hydrogen peroxide (H2O2). The origin of this disparity remains unknown due to the lack of atomistic information on the catalase-like reaction in HRP. Using QM(DFT)/MM metadynamics simulations, we uncover the mechanism for reduction of the HRP Compound I intermediate by H2O2 at atomic detail. The reaction begins with a hydrogen atom transfer, forming a peroxyl radical and a Compound II-like species. Reorientation of the peroxyl radical in the active site, concomitant with the transfer of the second hydrogen atom, is the rate-limiting step, with a computed free energy barrier (18.7 kcal/mol, ∼ 6 kcal/mol higher than the one obtained for catalase) in good agreement with experiments. Our simulations reveal the crucial role played by the distal pocket residues in accommodating H2O2, enabling formation of a Compound II-like intermediate, similar to catalases. However, out of the two pathways for Compound II reduction found in catalases, only one is operative in HRP. Moreover, the hydrogen bond network in the distal side of HRP compensates less efficiently than in catalases for the energetic cost required to reorient the peroxyl radical at the rate-determining step. The distal Arg and a water molecule in the "wet" active site of HRP have a substantial impact on the reaction barrier, compared to the "dry" active site in catalase. Therefore, the lower catalase-like efficiency of heme peroxidases compared to catalases can be directly attributed to the different distal pocket architecture, providing hints to engineer peroxidases with a higher rate of H2O2 disproportionation.

Structural Determinants for the Binding of Morphinan Agonists to the µ-Opioid Receptor.

Atomistic descriptions of the µ-opioid receptor (µOR) noncovalently binding with two of its prototypical morphinan agonists, morphine (MOP) and hydromorphone (HMP), are investigated using molecular dynamics (MD) simulations. Subtle differences between the binding modes and hydration properties of MOP and HMP emerge from the calculations. Alchemical free energy perturbation calculations show qualitative agreement with in vitro experiments performed in this work: indeed, the binding free energy difference between MOP and HMP computed by forward and backward alchemical transformation is 1.2±1.1 and 0.8±0.8 kcal/mol, respectively, to be compared with 0.4±0.3 kcal/mol from experiment. Comparison with an MD simulation of µOR covalently bound with the antagonist ß-funaltrexamine hints to agonist-induced conformational changes associated with an early event of the receptor's activation: a shift of the transmembrane helix 6 relative to the transmembrane helix 3 and a consequent loss of the key R165-T279 interhelical hydrogen bond. This finding is consistent with a previous proposal suggesting that the R165-T279 hydrogen bond between these two helices indicates an inactive receptor conformation.

Keys to Lipid Selection in Fatty Acid Amide Hydrolase Catalysis: Structural Flexibility, Gating Residues and Multiple Binding Pockets.

The fatty acid amide hydrolase (FAAH) regulates the endocannabinoid system cleaving primarily the lipid messenger anandamide. FAAH has been well characterized over the years and, importantly, it represents a promising drug target to treat several diseases, including inflammatory-related diseases and cancer. But its enzymatic mechanism for lipid selection to specifically hydrolyze anandamide, rather than similar bioactive lipids, remains elusive. Here, we clarify this mechanism in FAAH, examining the role of the dynamic paddle, which is formed by the gating residues Phe432 and Trp531 at the boundary between two cavities that form the FAAH catalytic site (the "membrane-access" and the "acyl chain-binding" pockets). We integrate microsecond-long MD simulations of wild type and double mutant model systems (Phe432Ala and Trp531Ala) of FAAH, embedded in a realistic membrane/water environment, with mutagenesis and kinetic experiments. We comparatively analyze three fatty acid substrates with different hydrolysis rates (anandamide > oleamide > palmitoylethanolamide). Our findings identify FAAH's mechanism to selectively accommodate anandamide into a multi-pocket binding site, and to properly orient the substrate in pre-reactive conformations for efficient hydrolysis that is interceded by the dynamic paddle. Our findings therefore endorse a structural framework for a lipid selection mechanism mediated by structural flexibility and gating residues between multiple binding cavities, as found in FAAH. Based on the available structural data, this exquisite catalytic strategy for substrate specificity seems to be shared by other lipid-degrading enzymes with similar enzymatic architecture. The mechanistic insights for lipid selection might assist de-novo enzyme design or drug discovery efforts.