Molecular dynamics study of the mechanical properties of drug loaded model systems: A comparison of a polymersome with a bilayer.

In this work we implement a new methodology to study structural and mechanical properties of systems having spherical and planar symmetries throughout Molecular Dynamics simulations. This methodology is applied here to a drug delivery system based in polymersomes, as an example. The chosen model drug was the local anesthetic prilocaine due to previous parameterization within the used coarse grain scheme. In our approach, mass density profiles (MDPs) are used to obtain key structural parameters of the systems, and pressure profiles are used to estimate the curvature elastic parameters. The calculation of pressure profiles and radial MPDs required the development of specific methods, which were implemented in an in-house built version of the GROMACS 2018 code. The methodology presented in this work is applied to characterize poly(ethylene oxide)-poly(butadiene) polymersomes and bilayers loaded with the model drug prilocaine. Our results show that structural properties of the polymersome membrane could be obtained from bilayer simulations, with significantly lower computational cost compared to whole polymersome simulations, but the bilayer simulations are insufficient to get insights on their mechanical aspects, since the elastic parameters are canceled out for the complete bilayer (as consequence of the symmetry). The simulations of entire polymersomes, although more complex, offer a complementary approach to get insights on the mechanical behavior of the systems.

OpenEP: an open-source simulator for electroporation-based tumor treatments.

Electroporation (EP), the increase of cell membrane permeability due to the application of electric pulses, is a universal phenomenon with a broad range of applications. In medicine, some of the foremost EP-based tumor treatments are electrochemotherapy (ECT), irreversible electroporation, and gene electrotransfer (GET). The electroporation phenomenon is explained as the formation of cell membrane pores when a transmembrane cell voltage reaches a threshold value. Predicting the outcome of an EP-based tumor treatment consists of finding the electric field distribution with an electric threshold value covering the tumor (electroporated tissue). Threshold and electroporated tissue are also a function of the number of pulses, constituting a complex phenomenon requiring mathematical modeling. We present OpenEP, an open-source specific purpose simulator for EP-based tumor treatments, modeling among other variables, threshold, and electroporated tissue variations in time. Distributed under a free/libre user license, OpenEP allows the customization of tissue type; electrode geometry and material; pulse type, intensity, length, and frequency. OpenEP facilitates the prediction of an optimal EP-based protocol, such as ECT or GET, defined as the critical pulse dosage yielding maximum electroporated tissue with minimal damage. OpenEP displays a highly efficient shared memory implementation by taking advantage of parallel resources; this permits a rapid prediction of optimal EP-based treatment efficiency by pulse number tuning.

Fluorescence correlation spectroscopy reveals the dynamics of kinesins interacting with organelles during microtubule-dependent transport in cells.

Microtubule-dependent motors usually work together to transport organelles through the crowded intracellular milieu. Thus, transport performance depends on how motors organize on the cargo. Unfortunately, the lack of methodologies capable of measuring this organization in cells determines that many aspects of the collective action of motors remain elusive. Here, we combined fluorescence fluctuations and single particle tracking techniques to address how kinesins organize on rod-like mitochondria moving along microtubules in cells. This methodology simultaneously provides mitochondria trajectories and EGFP-tagged kinesin-1 intensity at different mitochondrial positions with millisecond resolution. We show that kinesin exchange at the mitochondrion surface is within ~100 ms and depends on the organelle speed. During anterograde transport, the mitochondrial leading tip presents slower motor exchange in comparison to the rear tip. In contrast, retrograde mitochondria show similar exchange rates of kinesins at both tips. Numerical simulations provide theoretical support to these results and evidence that motors do not share the load equally during intracellular transport.

Faithfulness-boost effect: Loyal teammate selection correlates with skill acquisition improvement in online games.

The problem of skill acquisition is ubiquitous and fundamental to life. Most tasks in modern society involve the cooperation with other subjects. Notwithstanding its fundamental importance, teammate selection is commonly overlooked when studying learning. We exploit the virtually infinite repository of human behavior available in Internet to study a relevant topic in anthropological science: how grouping strategies may affect learning. We analyze the impact of team play strategies in skill acquisition using a turn-based game where players can participate individually or in teams. We unveil a subtle but strong effect in skill acquisition based on the way teams are formed and maintained during time. "Faithfulness-boost effect" provides a skill boost during the first games that would only be acquired after thousands of games. The tendency to play games in teams is associated with a long-run skill improvement while playing loyally with the same teammate significantly accelerates short-run skill acquisition.

Mechanical properties of drug loaded diblock copolymer bilayers: A molecular dynamics study.

In this work, we present results of coarse-grained simulations to study the encapsulation of prilocaine (PLC), both neutral and protonated, on copolymer bilayers through molecular dynamics simulations. Using a previously validated membrane model, we have simulated loaded bilayers at different drug concentrations and at low (protonated PLC) and high (neutral PLC) pH levels. We have characterized key structural parameters of the loaded bilayers in order to understand the effects of encapsulation of PLC on the bilayer structure and mechanical properties. Neutral PLC was encapsulated in the hydrophobic region leading to a thickness increase, while the protonated species partitioned between the water phase and the poly(ethylene oxide)-poly(butadiene) (PBD) interface, relaxing the PBD region and leading to a decrease in the thickness. The tangential pressures of the studied systems were calculated, and their components were decomposed in order to gain insights on their compensation. In all cases, it is observed that the loading of the membrane does not significantly decrease the stability of the bilayer, indicating that the system could be used for drug delivery.

Mapping the dynamical organization of the cell nucleus through fluorescence correlation spectroscopy.

The hierarchical organization of the cell nucleus into specialized open reservoirs and the nucleoplasm overcrowding impose restrictions to the mobility of biomolecules and their interactions with nuclear targets. These properties determine that many nuclear functions such as transcription, replication, splicing or DNA repair are regulated by complex, dynamical processes that do not follow simple rules. Advanced fluorescence microscopy tools and, in particular, fluorescence correlation spectroscopy (FCS) provide complementary and exquisite information on the dynamics of fluorescent labeled molecules moving through the nuclear space and are helping us to comprehend the complexity of the nuclear structure. Here, we describe how FCS methods can be applied to reveal the dynamical organization of the nucleus in live cells. Specifically, we provide instructions for the preparation of cellular samples with fluorescent tagged proteins and detail how FCS can be easily instrumented in commercial confocal microscopes. In addition, we describe general rules to set the parameters for one and two-color experiments and the required controls for these experiments. Finally, we review the statistical analysis of the FCS data and summarize the use of numerical simulations as a complementary approach that helps us to understand the complex matrix of molecular interactions network within the nucleus.

Diffusion of single dye molecules in hydrated TiO2 mesoporous films.

Mesoporous oxide films are attractive frameworks in technological areas such as catalysis, sensing, environmental protection, and photovoltaics. Herein, we used fluorescence correlation spectroscopy to explore how the pore dimensions of hydrated TiO2 mesoporous calcined films modulate the molecular diffusion. Rhodamine B molecules in mesoporous films follow a Fickian process 2-3 orders slower compared to the probe in water. The mobility increases with the pore and neck radii reaching an approximately constant value for a neck radius >2.8 nm. However, the pore size does not control the dye diffusion at low ionic strength emphasizing the relevance of the probe interactions with the pore walls on dye mobility. In conclusion, our results show that the thermal conditioning of TiO2 mesoporous films provides an exceptional tool for controlling the pore and neck radii on the nanometer scale and has a major impact on molecular diffusion within the mesoporous network.

Diblock copolymer bilayers as model for polymersomes: A coarse grain approach.

This paper presents a new model for polymersomes developed using a poly(ethylene oxide)-poly(butadiene) diblock copolymer bilayer. The model is based on a coarse-grained approach using the MARTINI force field. Since no MARTINI parameters exist for poly(butadiene), we have refined these parameters using quantum mechanical calculations and molecular dynamics simulations. The model has been validated using extensive molecular dynamics simulations in systems with several hundred polymer units and reaching up to 6 µs. These simulations show that the copolymer coarse grain model self-assemble into bilayers and that NPT and NPNγT ensemble runs reproduce key structural and mechanical experimental properties for different copolymer length chains with a similar hydrophilic weight fraction.

Minimum free-energy paths for the self-organization of polymer brushes.

A methodology to calculate minimum free-energy paths based on the combination of a molecular theory and the improved string method is introduced and applied to study the self-organization of polymer brushes under poor solvent conditions. Polymer brushes in a poor solvent cannot undergo macroscopic phase separation due to the physical constraint imposed by the grafting points; therefore, they microphase separate forming aggregates. Under some conditions, the theory predicts that the homogeneous brush and the aggregates can exist as two different minima of the free energy. The theoretical methodology introduced in this work allows us to predict the minimum free-energy path connecting these two minima as well as the morphology of the system along the path. It is shown that the transition between the homogeneous brush and the aggregates may involve a free-energy barrier or be barrierless depending on the relative stability of the two morphologies and the chain length and grafting density of the polymer. In the case where a free-energy barrier exists, one of the morphologies is a metastable structure and, therefore, the properties of the brush as the quality of the solvent is cycled are expected to display hysteresis. The theory is also applied to study the adhesion/deadhesion transition between two opposing surfaces modified by identical polymer brushes and it is shown that this process may also require surpassing a free-energy barrier.

Long-Lived Binding of Sox2 to DNA Predicts Cell Fate in the Four-Cell Mouse Embryo.

Transcription factor (TF) binding to DNA is fundamental for gene regulation. However, it remains unknown how the dynamics of TF-DNA interactions change during cell-fate determination in vivo. Here, we use photo-activatable FCS to quantify TF-DNA binding in single cells of developing mouse embryos. In blastocysts, the TFs Oct4 and Sox2, which control pluripotency, bind DNA more stably in pluripotent than in extraembryonic cells. By contrast, in the four-cell embryo, Sox2 engages in more long-lived interactions than does Oct4. Sox2 long-lived binding varies between blastomeres and is regulated by H3R26 methylation. Live-cell tracking demonstrates that those blastomeres with more long-lived binding contribute more pluripotent progeny, and reducing H3R26 methylation decreases long-lived binding, Sox2 target expression, and pluripotent cell numbers. Therefore, Sox2-DNA binding predicts mammalian cell fate as early as the four-cell stage. More generally, we reveal the dynamic repartitioning of TFs between DNA sites driven by physiological epigenetic changes. VIDEO ABSTRACT.

Exploring the dynamics of cell processes through simulations of fluorescence microscopy experiments.

Fluorescence correlation spectroscopy (FCS) methods are powerful tools for unveiling the dynamical organization of cells. For simple cases, such as molecules passively moving in a homogeneous media, FCS analysis yields analytical functions that can be fitted to the experimental data to recover the phenomenological rate parameters. Unfortunately, many dynamical processes in cells do not follow these simple models, and in many instances it is not possible to obtain an analytical function through a theoretical analysis of a more complex model. In such cases, experimental analysis can be combined with Monte Carlo simulations to aid in interpretation of the data. In response to this need, we developed a method called FERNET (Fluorescence Emission Recipes and Numerical routines Toolkit) based on Monte Carlo simulations and the MCell-Blender platform, which was designed to treat the reaction-diffusion problem under realistic scenarios. This method enables us to set complex geometries of the simulation space, distribute molecules among different compartments, and define interspecies reactions with selected kinetic constants, diffusion coefficients, and species brightness. We apply this method to simulate single- and multiple-point FCS, photon-counting histogram analysis, raster image correlation spectroscopy, and two-color fluorescence cross-correlation spectroscopy. We believe that this new program could be very useful for predicting and understanding the output of fluorescence microscopy experiments.

GPU Accelerated Implementation of Density Functional Theory for Hybrid QM/MM Simulations.

The hybrid simulation tools (QM/MM) evolved into a fundamental methodology for studying chemical reactivity in complex environments. This paper presents an implementation of electronic structure calculations based on density functional theory. This development is optimized for performing hybrid molecular dynamics simulations by making use of graphic processors (GPU) for the most computationally demanding parts (exchange-correlation terms). The proposed implementation is able to take advantage of modern GPUs achieving acceleration in relevant portions between 20 to 30 times faster than the CPU version. The presented code was extensively tested, both in terms of numerical quality and performance over systems of different size and composition.