Effects of Linking Topology on the Shear Response of Connected Ring Polymers: Catenanes and Bonded Rings Flow Differently

We perform computer simulations of mechanically linked (poly[2]catenanes, PC) and chemically bonded (bonded rings, BR) pairs of self-avoiding ring polymers in steady shear. We find that BRs develop a novel motif, termed gradient tumbling, rotating around the gradient axis. For the PCs the rings are stretched and display another new pattern, termed slip tumbling. The dynamics of BRs is continuous and oscillatory, whereas that of PCs is intermittent between slip-tumbling attempts. Our findings demonstrate the interplay between topology and hydrodynamics in dilute solutions of connected polymers.

In silico study of the impact of oxidation on pyruvate transmission across the hVDAC1 protein channel

The overexpression of voltage dependent anion channels (VDACs), particularly VDAC1, in cancer cells compared to normal cells, plays a crucial role in cancer cell metabolism, apoptosis regulation, and energy homeostasis. In this study, we used molecular dynamics (MD) simulations to investigate the effect of a low level of VDAC1 oxidation (induced e.g., by cold atmospheric plasma (CAP)) on the pyruvate (Pyr) uptake by VDAC1. Inhibiting Pyr uptake through VDAC1 can suppress cancer cell proliferation. Our primary target was to study the translocation of Pyr across the native and oxidized forms of hVDAC1, the human VDAC1. Specifically, we employed MD simulations to analyze the hVDAC1 structure by modifying certain cysteine residues to cysteic acids and methionine residues to methionine sulfoxides, which allowed us to investigate the effect of oxidation. Our results showed that the free energy barrier for Pyr translocation through the native and oxidized channel was approximately 4.3 ± 0.7 kJ mol-1 and 10.8 ± 1.8 kJ mol-1, respectively. An increase in barrier results in a decrease in rate of Pyr permeation through the oxidized channel. Thus, our results indicate that low levels of CAP oxidation reduce Pyr translocation, resulting in decreased cancer cell proliferation. Therefore, low levels of oxidation are likely sufficient to treat cancer cells given the inhibition of Pyr uptake.

Electrochromic Sensor Augmented with Machine Learning for Enzyme-Free Analysis of Antioxidants

Our study presents an electrochromic sensor that operates without the need for enzymes or multiple oxidant reagents. This sensor is augmented with machine learning algorithms, enabling the identification, classification, and prediction of six different antioxidants with high accuracy. We utilized polyaniline (PANI), Prussian blue (PB), and copper-Prussian blue analogues (Cu-PBA) at their respective oxidation states as electrochromic materials (ECMs). By designing three readout channels with these materials, we were able to achieve visual detection of antioxidants without relying on traditional "lock and key" specific interactions. Our sensing approach is based on the direct electrochemical reactions between oxidized electrochromic materials (ECMsox) as electron acceptors and various antioxidants, which act as electron donors. This interaction generates unique fingerprint patterns by switching the ECMsox to reduced electrochromic materials (ECMsred), causing their colors to change. Through the application of density functional theory (DFT), we demonstrated the molecular-level basis for the distinct multicolor patterns. Additionally, machine learning algorithms were employed to correlate the optical patterns with RGB data, enabling complex data analysis and the prediction of unknown samples. To demonstrate the practical applications of our design, we successfully used the EC sensor to diagnose antioxidants in serum samples, indicating its potential for the on-site monitoring of antioxidant-related diseases. This advancement holds promise for various applications, including the real-time monitoring of antioxidant levels in biological samples, the early diagnosis of antioxidant-related diseases, and personalized medicine. Furthermore, the success of our electrochromic sensor design highlights the potential for exploring similar strategies in the development of sensors for diverse analytes, showcasing the versatility and adaptability of this approach.

Dynamics of fluid bilayer vesicles: Soft meshes and robust curvature energy discretization

Continuum models like the Helfrich Hamiltonian are widely used to describe fluid bilayer vesicles. Here we study the molecular dynamics compatible dynamics of the vertices of two-dimensional meshes representing the bilayer, whose in-plane motion is only weakly constrained. We show (i) that Jülicher's discretization of the curvature energy offers vastly superior robustness for soft meshes compared to the commonly employed expression by Gommper and Kroll and (ii) that for sufficiently soft meshes, the typical behavior of fluid bilayer vesicles can emerge even if the mesh connectivity remains fixed throughout the simulations. In particular, soft meshes can accommodate large shape transformations, and the model can generate the typical ℓ^{-4} signal for the amplitude of surface undulation modes of nearly spherical vesicles all the way up to the longest wavelength modes. Furthermore, we compare results for Newtonian, Langevin, and Brownian dynamics simulations of the mesh vertices to demonstrate that the internal friction of the membrane model is negligible, making it suitable for studying the internal dynamics of vesicles via coupling to hydrodynamic solvers or particle-based solvent models.

Counterintuitive properties of evolutionary measures: A stochastic process study in cyclic population structures with periodic environments

Local environmental interactions are a major factor in determining the success of a new mutant in structured populations. Spatial variations in the concentration of genotype-specific resources change the fitness of competing strategies locally and thus can drastically change the outcome of evolutionary processes in unintuitive ways. The question is how such local environmental variations in network population structures change the condition for selection and fixation probability of an advantageous (or deleterious) mutant. We consider linear graph structures and focus on the case where resources have a spatial periodic pattern. This is the simplest model with two parameters, length scale and fitness scales, representing heterogeneity. We calculate fixation probability and fixation times for a constant population birth-death process as fitness heterogeneity and period vary. Fixation probability is affected by not only the level of fitness heterogeneity but also spatial scale of resources variations set by period of distribution T. We identify conditions for which a previously a deleterious mutant (in a uniform environment) becomes beneficial as fitness heterogeneity is increased. We observe cases where the fixation probability of both mutant and resident types are more than their neutral value, 1/N, simultaneously. This coincides with exponential increase in time to fixation which points to potential coexistence of resident and mutant types. Finally, we discuss the effect of the 'fitness shift' where the fitness function of two types has a phase difference. We observe significant increases (or decreases) in the fixation probability of the mutant as a result of such phase shift.

Programmable Transport of C60 by Straining Graphene Substrate

Controlling the maneuverability of nanocars and molecular machines on the surface is essential for the targeted transportation of materials and energy at the nanoscale. Here, we evaluate the motion of fullerene, as the most popular candidate for use as a nanocar wheel, on the graphene nanoribbons with strain gradients based on molecular dynamics (MD), and theoretical approaches. The strain of the examined substrates linearly decreases by 20%, 16%, 12%, 8%, 4%, and 2%. MD calculations were performed with the open source LAMMPS solver. The essential physics of the interactions is captured by Lennard-Jones and Tersoff potentials. The motion of C60 on the graphene nanoribbon is simulated in canonical ensemble, which is implanted by using a Nose-Hoover thermostat. Since the potential energy of C60 is lower on the unstrained end of nanoribbons, this region is energetically more favorable for the molecule. As the strain gradient of the surface increases, the trajectories of the motion and the C60 velocity indicate more directed movements along the gradient of strain on the substrate. Based on the theoretical relations, it was shown that the driving force and diffusion coefficient of the C60 motion respectively find linear and quadratic growth with the increase of strain gradient, which is confirmed by MD simulations. To understand the effect of temperature, at each strain gradient of substrate, the simulations are repeated at the temperatures of 100, 200, 300, and 400 K. The large ratio of longitudinal speed to the transverse speed of fullerene at 100 and 200 K refers to the rectilinear motion of molecule at low temperatures. Using successive strain gradients on the graphene in perpendicular directions, we steered the motion of C60 to the desired target locations. The programmable transportation of nanomaterials on the surface has a significant role in different processes at the nanoscale, such as bottom-up assembly.

Collective movement and thermal stability of fullerene clusters on the graphene layer

Understanding the motion characteristics of fullerene clusters on the graphene surface is critical for designing surface manipulation systems. Toward this purpose, using the molecular dynamics method, we evaluated six clusters of fullerenes including 1, 2, 3, 5, 10, and 25 molecules on the graphene surface, in the temperature range of 25 to 500 K. First, the surface motion of clusters is studied at 200 K and lower temperatures, in which fullerenes remain as a single group. The trajectories of the motion as well as the diffusion coefficients indicate the reduction of surface mobility as a response to the increase of the fullerene number. The clusters show normal diffusion at the temperature of 25 K, while they follow the super-diffusion regime at higher temperatures. The separation of fullerenes occurs at 300 K and higher temperatures. Due to the increase of vdW attraction with the increase of the fullerene number, the separation of fullerenes in larger clusters occurs at higher temperatures. The thermal energy at 500 K is sufficient to divide the large C60 clusters into smaller clusters. This energy level is related to the saturation of the interaction energy experienced by individual fullerenes, which can be estimated from the potential energy analysis. The results of simulations reveal that the separation occurs at the edge of clusters. Moreover, we studied the thermal stability of multilayer fullerene clusters on graphene. The simulation results indicate the tendency of multilayer clusters to locate on the surface, which implies the wetting property of C60s on the graphene layer.

A modified Jarzynski free-energy estimator to eliminate non-conservative forces and its application in nanoparticle-membrane interactions

Computational methods to understand interactions in bio-complex systems are however limited to time-scales typically much shorter than in Nature. For example, on the nanoscale level, interactions between nanoparticles (NPs)/molecules/peptides and membranes are central in complex biomolecular processes such as membrane-coated NPs or cellular uptake. This can be remedied by the application of e.g. Jarzynski's equality where thermodynamic properties are extracted from non-equilibrium simulations. Although, the out of equilibrium work leads to non-conservative forces. We here propose a correction Pair Forces method, that removes these forces. Our proposed method is based on the calculation of pulling forces in backward and forward directions for the Jarzynski free-energy estimator using steered molecular dynamics simulation. Our results show that this leads to much improvement for NP-membrane translocation free energies. Although here we have demonstrated the application of the method in molecular dynamics simulation, it could be applied for experimental approaches.

Conformation- and phosphorylation-dependent electron tunnelling across self-assembled monolayers of tau peptides

We report on charge transport across self-assembled monolayers (SAMs) of short tau peptides by probing the electron tunneling rates and quantum mechanical simulation. We measured the electron tunneling rates across SAMs of carboxyl-terminated linker molecules (C₆H₁₂O₂S) and short cis-tau (CT) and trans-tau (TT) peptides, supported on template-stripped gold (Au^TS) bottom electrode, with Eutectic Gallium-Indium (EGaIn)(EGaIn) top electrode. Measurements of the current density across thousands of Au^TS/linker/tau//Ga₂O₃/EGaIn single-molecule junctions show that the tunneling current across CT peptide is one order of magnitude lower than that of TT peptide. Quantum mechanical simulation demonstrated a wider energy bandgap of the CT peptide, as compared to the TT peptide, which causes a reduction in its electron tunneling current. Our findings also revealed the critical role of phosphorylation in altering the charge transport characteristics of short peptides; more specifically, we found that the presence of phosphate groups can reduce the energy band gap in tau peptides and alter their electrical properties. Our results suggest that conformational and phosphorylation of short peptides (e.g., tau) can significantly change their charge transport characteristics and energy levels.

Effect of Cysteine Oxidation in SARS-CoV-2 Receptor-Binding Domain on Its Interaction with Two Cell Receptors: Insights from Atomistic Simulations

Binding of the SARS-CoV-2 S-glycoprotein to cell receptors is vital for the entry of the virus into cells and subsequent infection. ACE2 is the main cell receptor for SARS-CoV-2, which can attach to the C-terminal receptor-binding domain (RBD) of the SARS-CoV-2 S-glycoprotein. The GRP78 receptor plays an anchoring role, which attaches to the RBD and increases the chance of other RBDs binding to ACE2. Although high levels of reactive oxygen and nitrogen species (RONS) are produced during viral infections, it is not clear how they affect the RBD structure and its binding to ACE2 and GRP78. In this research, we apply molecular dynamics simulations to study the effect of oxidation of the highly reactive cysteine (Cys) amino acids of the RBD on its binding to ACE2 and GRP78. The interaction energy of both ACE2 and GRP78 with the whole RBD, as well as with the RBD main regions, is compared in both the native and oxidized RBDs. Our results show that the interaction energy between the oxidized RBD and ACE2 is strengthened by 155 kJ/mol, increasing the binding of the RBD to ACE2 after oxidation. In addition, the interaction energy between the RBD and GRP78 is slightly increased by 8 kJ/mol after oxidation, but this difference is not significant. Overall, these findings highlight the role of RONS in the binding of the SARS-CoV-2 S-glycoprotein to host cell receptors and suggest an alternative mechanism by which RONS could modulate the entrance of viral particles into the cells.

Nanocar Swarm Movement on Graphene Surface

Investigation of the nanomachines swarm motion is useful in the design of molecular transportation systems as well as in understanding the assembly process on the surface. Here, we evaluate the...

Polarization Switching Method for Effective Free Energy Calculation of Membrane Translocation on the Nano-scale

Free-energy calculations are crucial for investigating biomolecular interactions on the Nano-scale level. However, in theoretical studies, the neglect of electronic polarization can jeopardize their accuracy and correct predictive capabilities, specifically...

Monte Carlo simulation of a lattice model for the dynamics of randomly branching double-folded ring polymers

Supercoiled DNA, crumpled interphase chromosomes, and topologically constrained ring polymers often adopt treelike, double-folded, randomly branching configurations. Here we study an elastic lattice model for tightly double-folded ring polymers, which allows for the spontaneous creation and deletion of side branches coupled to a diffusive mass transport, which is local both in space and on the connectivity graph of the tree. We use Monte Carlo simulations to study systems falling into three different universality classes: ideal double-folded rings without excluded volume interactions, self-avoiding double-folded rings, and double-folded rings in the melt state. The observed static properties are in good agreement with exact results, simulations, and predictions of Flory theory for randomly branching polymers. For example, in the melt state rings adopt compact configurations and exhibit territorial behavior. In particular, we show that the emergent dynamics is in excellent agreement with a recent scaling theory and illustrate the qualitative differences with the familiar reptation dynamics of linear chains.

S494 O-glycosylation site on the SARS-CoV-2 RBD affects the virus affinity to ACE2 and its infectivity; a molecular dynamics study

SARS-CoV-2 is a strain of Coronavirus family that caused the ongoing pandemic of COVID-19. Several studies showed that the glycosylation of virus spike (S) protein and the Angiotensin-Converting Enzyme 2 (ACE2) receptor on the host cell is critical for the virus infectivity. Molecular Dynamics (MD) simulations were used to explore the role of a novel mutated O-glycosylation site (D494S) on the Receptor Binding Domain (RBD) of S protein. This site was suggested as a key mediator of virus-host interaction. By exploring the dynamics of three O-glycosylated models and the control systems of unglcosylated S4944 and S494D complexes, it was shown that the decoration of S494 with elongated O-glycans results in stabilized interactions on the direct RBD-ACE2. Calculation of the distances between RBD and two major H1, H2 helices of ACE2 and the interacting pairs of amino acids in the interface showed that the elongated O-glycan maintains these interactions by forming several polar contacts with the neighbouring residues while it would not interfere in the direct binding interface. Relative binding free energy of RBD-ACE2 is also more favorable in the O-glycosylated models with longer glycans. The increase of RBD binding affinity to ACE2 depends on the size of attached O-glycan. By increasing the size of O-glycan, the RBD-ACE2 binding affinity will increase. Hence, this crucial factor must be taken into account for any further inhibitory approaches towards RBD-ACE2 interaction.

Impact of temporal correlations on high risk outbreaks of independent and cooperative SIR dynamics

We first propose a quantitative approach to detect high risk outbreaks of independent and coinfective SIR dynamics on three empirical networks: a school, a conference and a hospital contact network. This measurement is based on the k-means clustering method and identifies proper samples for calculating the mean outbreak size and the outbreak probability. Then we systematically study the impact of different temporal correlations on high risk outbreaks over the original and differently shuffled counterparts of each network. We observe that, on the one hand, in the coinfection process, randomization of the sequence of the events increases the mean outbreak size of high-risk cases. On the other hand, these correlations do not have a consistent effect on the independent infection dynamics, and can either decrease or increase this mean. Randomization of the daily pattern correlations has no strong impact on the size of the outbreak in either the coinfection or the independent spreading cases. We also observe that an increase in the mean outbreak size does not always coincide with an increase in the outbreak probability; therefore, we argue that merely considering the mean outbreak size of all realizations may lead us into falsely estimating the outbreak risks. Our results suggest that some sort of contact randomization in the organizational level in schools, events or hospitals might help to suppress the spreading dynamics while the risk of an outbreak is high.

Locomotion of the C60-based nanomachines on graphene surfaces

We provide a comprehensive computational characterization of surface motion of two types of nanomachines with four C₆₀ "wheels": a flexible chassis Nanocar and a rigid chassis Nanotruck. We study the nanocars' lateral and rotational diffusion as well as the wheels' rolling motion on two kinds of graphene substrates-flexible single-layer graphene which may form surface ripples and an ideally flat graphene monolayer. We find that the graphene surface ripples facilitate the translational diffusion of Nanocar and Nanotruck, but have little effect on their surface rotation or the rolling of their wheels. The latter two types of motion are strongly affected by the structure of the nanomachines instead. Surface diffusion of both nanomachines occurs preferentially via a sliding mechanism whereas the rolling of the "wheels" contributes little. The axial rotation of all "wheels" is uncorrelated.

Thermal conductivity of the cell membrane in the presence of cholesterol and amyloid precursor protein

The cell membrane is responsible for the transportation of heat between inside and outside the cell. Whether the thermal properties of the cell membrane are affected by the cholesterol concentration or the membrane proteins has not been investigated so far. Although the experimental measurement of the membrane thermal conductivity was not available until very recently, computational methods have been widely used for this purpose. In this study, we carry out molecular dynamics simulations to investigate the relation between the concentration of cholesterol and the thermal conductivity of a model membrane. Our results suggest an increase in the membrane thermal conductivity upon increasing the concentration of cholesterol in the membrane. Moreover, we find that the asymmetric distribution of cholesterol in the two membrane leaflets decreases thermal conductivity. We also find a rectification effect when heat flows in opposite directions through a model membrane decorated with the amyloid precursor protein. The results of this study apply to the advancement of selective treatment methods, as well as the development of new materials such as biological rectifiers.

Glycan-mediated functional assembly of IL-1RI: structural insights into completion of the current description for immune response

Interleukin 1 Receptor type I (IL-1RI) is a multi-domain transmembrane receptor that triggers the inflammatory response. Understanding its detailed mechanism of action is crucial for treating immune disorders. IL-1RI is activated upon formation of its functional assembly that occurs by binding of the IL-1 cytokine and the accessory protein (Il-1RAcP) to it. X-ray crystallography, small-Angle X-ray Scattering and molecular dynamics simulation studies showed that IL-1RI adopts two types of 'compact' and 'extended' conformational states in its dynamical pattern. Furthermore, glycosylation has shown to play a critical role in its activation process. Here, classical and accelerated atomistic molecular dynamics were carried out to examine the role of full glycosylation of IL-1RI and IL-1RAcP in arrangement of the functional assembly. Simulations showed that the 'compact' and 'extended' IL-1RI form two types of 'cytokine-inaccessible-non-signaling' and 'cytokine-accessible-signaling' assemblies with the IL-1RacP, respectively that are both abiding in the presence of glycans. Suggesting that the cytokine binding to IL-1RI is not required for the formation of IL-1RI-IL-1RAcP complex and the 'compact' complex could act as a down-regulatory mechanism. The 'extended' complex is maintained by formation of several persistent hydrogen bonds between the IL-1RI-IL-1RAcP inter-connected glycans. Taken together, it was shown that full glycosylation regulates formation of the IL-1RI functional assembly and play critical role in cytokine biding and triggering the IL-1RI involved downstream pathways in the cell.Communicated by Ramaswamy H. Sarma.

GAG positioning on IL-1RI; A mechanism regulated by dual effect of glycosylation

IL-1RI is the signaling receptor for the IL-1 family of cytokines that are involved in establishment of the innate and acquired immune systems. Glycosylated extracellular (EC) domain of the IL-1RI binds to agonist such as IL-1β or antagonist ligands and the accessory protein to form the functional signaling complex. Dynamics and ligand binding of the IL-1RI is influenced by presence of the glycosaminoglycans (GAGs) of the EC matrix. Here a combination of molecular dockings and molecular dynamics simulations of the unglycosylated, partially N-glycosylated and fully N-glycosylated IL-1RI EC domain in the apo, GAG-bound and IL-1β-bound states were carried out to explain the co-occurring dynamical effect of receptor's glycosylation and GAGs. It was shown that the IL-1RI adopts two types of "extended" and "locked" conformations in its dynamical pattern, and glycosylation maintains the receptor in the latter form. Maintaining the receptor in the locked conformation disfavors IL-1β binding by burying its two binding site on the IL-1RI EC domain. Glycosylation disfavors GAG binding to the extended IL-1RI EC domain by sterically limiting the GAGs degrees of freedom in targeting its binding site, while it favors GAG binding to the locked IL-1RI by favorable packing interactions.

Effects of topological constraints on linked ring polymers in solvents of varying quality

We investigate the effects of topological constraints in catenanes composed of interlinked ring polymers on their size in a good solvent as well as on the location of their θ-point when the solvent quality is worsened. We mainly focus on poly[n]catenanes consisting of n ring polymers each of length m interlocked in a linear fashion. Using molecular dynamics simulations, we study the scaling of the poly[n]catenane's radius of gyration in a good solvent, assuming in general that R_g∼m^μn^ν and we find that μ = 0.65 ± 0.02 and ν = 0.60 ± 0.01 for the range of n and m considered. These findings are further rationalized with the help of a mean-field Flory-like theory yielding the values of μ = 16/25 and ν = 3/5, consistent with the numerical results. We show that individual rings within catenanes feature a surplus swelling due to the presence of N_L topological links. Furthermore, we consider poly[n]catenanes in solvents of varying quality and we demonstrate that the presence of topological links leads to an increase of its θ-temperature in comparison to isolated linear and ring chains with the following ordering: T > T > T. Finally, we show that the presence of links similarly raises the θ-temperature of a single linked ring in comparison to an unlinked one, bringing its θ-temperature close to the one of a poly[n]catenane.

Ruddlesden-Popper Phases of Methylammonium-Based Two-Dimensional Perovskites with 5-Ammonium Valeric Acid AVA2MA n-1Pb nI3 n+1 with n = 1, 2, and 3

5-Ammonium valeric acid (AVA) is a frequently used additive in the preparation of lead halide perovskites. However, its microscopic role as passivating, cross-linking, or templating agent is far from clear. In this work, we provide density functional theory-based structural models for the Ruddlesden-Popper (RP) phases of AVA₂(CH₃NH₃) _n-1Pb _nI_{3 n+1} for n = 1, 2, and 3 and validate with experimental data on polycrystalline samples for n = 1. The structural and electronic properties of the AVA-based RP phases are compared to the ones of other linker families. In contrast to aromatic and aliphatic spacers without additional functional groups, the RP phases of AVA are characterized by the formation of a regular and stable H-bonding network between the carbonyl head groups of adjacent AVA molecules in opposite layers. Because of these additional interactions, the penetration depth of the organic layer into the perovskite sheet is reduced with direct consequences for its crystalline phase. The possibility of forming strong interlinker hydrogen bonds may lead to an enhanced thermal stability.

Mechanistic Understanding of the Interactions between Nano-Objects with Different Surface Properties and α-Synuclein

Aggregation of the natively unfolded protein α-synuclein (α-syn) is key to the development of Parkinson's disease (PD). Some nanoparticles (NPs) can inhibit this process and in turn be used for treatment of PD. Using simulation strategies, we show here that α-syn self-assembly is electrostatically driven. Dimerization by head-to-head monomer contact is triggered by dipole-dipole interactions and subsequently stabilized by van der Waals interactions and hydrogen bonds. Therefore, we hypothesized that charged nano-objects could interfere with this process and thus prevent α-syn fibrillation. In our simulations, positively and negatively charged graphene sheets or superparamagnetic iron oxide NPs first interacted with α-syn's N/C terminally charged residues and then with hydrophobic residues in the non-amyloid-β component (61-95) region. In the experimental setup, we demonstrated that the charged nano-objects have the capacity not only to strongly inhibit α-syn fibrillation (both nucleation and elongation) but also to disaggregate the mature fibrils. Through the α-syn fibrillation process, the charged nano-objects induced the formation of off-pathway oligomers.

Molecular interaction of fibrinogen with zeolite nanoparticles

Fibrinogen is one of the key proteins that participate in the protein corona composition of many types of nanoparticles (NPs), and its conformational changes are crucial for activation of immune systems. Recently, we demonstrated that the fibrinogen highly contributed in the protein corona composition at the surface of zeolite nanoparticles. Therefore, understanding the interaction of fibrinogen with zeolite nanoparticles in more details could shed light of their safe applications in medicine. Thus, we probed the molecular interactions between fibrinogen and zeolite nanoparticles using both experimental and simulation approaches. The results indicated that fibrinogen has a strong and thermodynamically favorable interaction with zeolite nanoparticles in a non-cooperative manner. Additionally, fibrinogen experienced a substantial conformational change in the presence of zeolite nanoparticles through a concentration-dependent manner. Simulation results showed that both E- and D-domain of fibrinogen are bound to the EMT zeolite NPs via strong electrostatic interactions, and undergo structural changes leading to exposing normally buried sequences. D-domain has more contribution in this interaction and the C-terminus of γ chain (γ377-394), located in D-domain, showed the highest level of exposure compared to other sequences/residues.

Molecular Dynamics Simulations of Orai Reveal How the Third Transmembrane Segment Contributes to Hydration and Ca2+ Selectivity in Calcium Release-Activated Calcium Channels

Calcium release-activated calcium (CRAC) channels open upon depletion of Ca²⁺ from the endoplasmic reticulum, and when open, they are permeable to a selective flux of calcium ions. The atomic structure of Orai, the pore domain of CRAC channels, from Drosophila melanogaster has revealed many details about conduction and selectivity in this family of ion channels. However, it is still unclear how residues on the third transmembrane helix can affect the conduction properties of the channel. Here, molecular dynamics and Brownian dynamics simulations were employed to analyze how a conserved glutamate residue on the third transmembrane helix (E262) contributes to selectivity. The comparison between the wild-type and mutated channels revealed a severe impact of the mutation on the hydration pattern of the pore domain and on the dynamics of residues K270, and Brownian dynamics simulations proved that the altered configuration of residues K270 in the mutated channel impairs selectivity to Ca²⁺ over Na⁺. The crevices of water molecules, revealed by molecular dynamics simulations, are perfectly located to contribute to the dynamics of the hydrophobic gate and the basic gate, suggesting a possible role in channel opening and in selectivity function.

Development of a Virtual Cell Model to Predict Cell Response to Substrate Topography

Cells can sense and respond to changes in the topographical, chemical, and mechanical information in their environment. Engineered substrates are increasingly being developed that exploit these physical attributes to direct cell responses (most notably mesenchymal stem cells) and therefore control cell behavior toward desired applications. However, there are very few methods available for robust and accurate modeling that can predict cell behavior prior to experimental evaluations, and this typically means that many cell test iterations are needed to identify best material features. Here, we developed a unifying computational framework to create a multicomponent cell model, called the "virtual cell model" that has the capability to predict changes in whole cell and cell nucleus characteristics (in terms of shape, direction, and even chromatin conformation) on a range of cell substrates. Modeling data were correlated with cell culture experimental outcomes in order to confirm the applicability of the virtual cell model and demonstrating the ability to reflect the qualitative behavior of mesenchymal stem cells. This may provide a reliable, efficient, and fast high-throughput approach for the development of optimized substrates for a broad range of cellular applications including stem cell differentiation.

Topology of polymer chains under nanoscale confinement

Spatial confinement limits the conformational space accessible to biomolecules but the implications for bimolecular topology are not yet known. Folded linear biopolymers can be seen as molecular circuits formed by intramolecular contacts. The pairwise arrangement of intra-chain contacts can be categorized as parallel, series or cross, and has been identified as a topological property. Using molecular dynamics simulations, we determine the contact order distributions and topological circuits of short semi-flexible linear and ring polymer chains with a persistence length of l_p under a spherical confinement of radius R_c. At low values of l_p/R_c, the entropy of the linear chain leads to the formation of independent contacts along the chain and accordingly, increases the fraction of series topology with respect to other topologies. However, at high l_p/R_c, the fraction of cross and parallel topologies are enhanced in the chain topological circuits with cross becoming predominant. At an intermediate confining regime, we identify a critical value of l_p/R_c, at which all topological states have equal probability. Confinement thus equalizes the probability of more complex cross and parallel topologies to the level of the more simple, non-cooperative series topology. Moreover, our topology analysis reveals distinct behaviours for ring- and linear polymers under weak confinement; however, we find no difference between ring- and linear polymers under strong confinement. Under weak confinement, ring polymers adopt parallel and series topologies with equal likelihood, while linear polymers show a higher tendency for series arrangement. The radial distribution analysis of the topology reveals a non-uniform effect of confinement on the topology of polymer chains, thereby imposing more pronounced effects on the core region than on the confinement surface. Additionally, our results reveal that over a wide range of confining radii, loops arranged in parallel and cross topologies have nearly the same contact orders. Such degeneracy implies that the kinetics and transition rates between the topological states cannot be solely explained by contact order. We expect these findings to be of general importance in understanding chaperone assisted protein folding, chromosome architecture, and the evolution of molecular folds.

Gaussian theory for spatially distributed self-propelled particles

Obtaining a reduced description with particle and momentum flux densities outgoing from the microscopic equations of motion of the particles requires approximations. The usual method, we refer to as truncation method, is to zero Fourier modes of the orientation distribution starting from a given number. Here we propose another method to derive continuum equations for interacting self-propelled particles. The derivation is based on a Gaussian approximation (GA) of the distribution of the direction of particles. First, by means of simulation of the microscopic model, we justify that the distribution of individual directions fits well to a wrapped Gaussian distribution. Second, we numerically integrate the continuum equations derived in the GA in order to compare with results of simulations. We obtain that the global polarization in the GA exhibits a hysteresis in dependence on the noise intensity. It shows qualitatively the same behavior as we find in particles simulations. Moreover, both global polarizations agree perfectly for low noise intensities. The spatiotemporal structures of the GA are also in agreement with simulations. We conclude that the GA shows qualitative agreement for a wide range of noise intensities. In particular, for low noise intensities the agreement with simulations is better as other approximations, making the GA to an acceptable candidates of describing spatially distributed self-propelled particles.

Stiffer double-stranded DNA in two-dimensional confinement due to bending anisotropy

Using analytical approach and Monte Carlo (MC) simulations, we study the elastic behavior of the intrinsically twisted elastic ribbons with bending anisotropy, such as double-stranded DNA (dsDNA), in two-dimensional (2D) confinement. We show that, due to the bending anisotropy, the persistence length of dsDNA in 2D conformations is always greater than three-dimensional (3D) conformations. This result is in consistence with the measured values for DNA persistence length in 2D and 3D in equal biological conditions. We also show that in two dimensions, an anisotropic, intrinsically twisted polymer exhibits an implicit twist-bend coupling, which leads to the transient curvature increasing with a half helical turn periodicity along the bent polymer.

Topology of internally constrained polymer chains

An interacting partner can provide external control over folding rates and realized topologies.

Active Brownian particles and run-and-tumble particles separate inside a maze

A diverse range of natural and artificial self-propelled particles are known and are used nowadays. Among them, active Brownian particles (ABPs) and run-and-tumble particles (RTPs) are two important classes. We numerically study non-interacting ABPs and RTPs strongly confined to different maze geometries in two dimensions. We demonstrate that by means of geometrical confinement alone, ABPs are separable from RTPs. By investigating Matryoshka-like mazes with nested shells, we show that a circular maze has the best filtration efficiency. Results on the mean first-passage time reveal that ABPs escape faster from the center of the maze, while RTPs reach the center from the rim more easily. According to our simulations and a rate theory, which we developed, ABPs in steady state accumulate in the outermost region of the Matryoshka-like mazes, while RTPs occupy all locations within the maze with nearly equal probability. These results suggest a novel technique for separating different types of self-propelled particles by designing appropriate confining geometries without using chemical or biological agents.

Extreme bendability of DNA double helix due to bending asymmetry

Experimental data of the DNA cyclization (J-factor) at short length scales exceed the theoretical expectation based on the wormlike chain (WLC) model by several orders of magnitude. Here, we propose that asymmetric bending rigidity of the double helix in the groove direction can be responsible for extreme bendability of DNA at short length scales and it also facilitates DNA loop formation at these lengths. To account for the bending asymmetry, we consider the asymmetric elastic rod (AER) model which has been introduced and parametrized in an earlier study [B. Eslami-Mossallam and M. R. Ejtehadi, Phys. Rev. E 80, 011919 (2009)]. Exploiting a coarse grained representation of the DNA molecule at base pair (bp) level and using the Monte Carlo simulation method in combination with the umbrella sampling technique, we calculate the loop formation probability of DNA in the AER model. We show that the DNA molecule has a larger J-factor compared to the WLC model which is in excellent agreement with recent experimental data.

Vortex with fourfold defect lines in a simple model of self-propelled particles

We study the formation of a vortex with fourfold symmetry in a minimal model of self-propelled particles, confined inside a squared box, using computer simulations and also theoretical analysis. In addition to the vortex pattern, we observe five other regimes in the system: a homogeneous gaseous phase, band structures, moving clumps, moving clusters, and vibrating rings. All six regimes emerge from controlling the strength of noise and from the contribution of repulsion and alignment interactions. We study the shape of the vortex and its symmetry in detail. The pattern shows exponential defect lines where incoming and outgoing flows of particles collide. We show that alignment and repulsion interactions between particles are necessary to form such patterns. We derive hydrodynamical equations with an introduction of the "small deviation" technique to describe the vortex phase. The method is applicable to other systems as well. Finally, we compare the theory with the results of both computer simulations and an experiment using Quincke rotors. A good agreement between the three is observed.

Cooperation within von Willebrand factors enhances adsorption mechanism

von Willebrand factor (VWF) is a naturally collapsed protein that participates in primary haemostasis and coagulation events. The clotting process is triggered by the adsorption and conformational changes of the plasma VWFs localized to the collagen fibres found near the site of injury. We develop coarse-grained models to simulate the adsorption dynamics of VWF flowing near the adhesive collagen fibres at different shear rates and investigate the effect of factors such as interaction and cooperativity of VWFs on the success of adsorption events. The adsorption probability of a flowing VWF confined to the receptor field is enhanced when it encounters an adhered VWF in proximity to the collagen receptors. This enhancement is observed within a wide range of shear rates and is mostly controlled by the attractive van der Waals interactions rather than the hydrodynamic interactions among VWF monomers. The cooperativity between the VWFs acts as an effective mechanism for enhancing VWF adsorption to the collagen fibres. Additionally, this implies that the adsorption of such molecules is nonlinearly dependent on the density of flowing VWFs. These findings are important for studies of primary haemostasis as well as general adsorption dynamics processes in polymer physics.

Diffusion and self-assembly of C60 molecules on monolayer graphyne sheets

The motion of a fullerene (C₆₀) on 5 different types of graphyne is studied by all-atom molecular dynamics simulations and compared with former studies on the motion of C₆₀ on graphene. The motion shows a diffusive behavior which consists of either a continuous motion or discrete movements between trapping sites depending on the type of the graphyne sheet. For graphyne-4 and graphyne-5, fullerenes could detach from the surface of the graphyne sheet at room temperature which was not reported for similar cases on graphene sheets. Collective motion of a group of fullerenes interacting with a graphyne studied and it is shown that fullerenes exhibit stable assemblies. Depending on the type of graphyne, these assemblies can have either single or double layers. The mobility of the assembled structures is also dependent on the type of the graphyne sheet. The observed properties of the motion suggests novel applications for the complexes of fullerene and monolayer graphynes.

Understanding the nanoparticle-protein corona complexes using computational and experimental methods

Nanoparticles (NP) have capability to adsorb proteins from biological fluids and form protein layer, which is called protein corona. As the cell sees corona coated NPs, the protein corona can dictate biological response to NPs. The composition of protein corona is varied by physicochemical properties of NPs including size, shape, surface chemistry. Processing of protein adsorption is dynamic phenomena; to that end, a protein may desorb or leave a surface vacancy that is rapidly filled by another protein and cause changes in the corona composition mainly by the Vroman effect. In this review, we discuss the interaction between NP and proteins and the available techniques for identification of NP-bound proteins. Also we review current developed computational methods for understanding the NP-protein complex interactions.

Low-Reynolds-number predator

To generalize simple bead-linker model of swimmers to higher dimensions and to demonstrate the chemotaxis ability of such swimmers, here we introduce a low-Reynolds predator, using a two-dimensional triangular bead-spring model. Two-state linkers as mechanochemical enzymes expand as a result of interaction with particular activator substances in the environment, causing the whole body to translate and rotate. The concentration of the chemical stimulator controls expansion versus the contraction rate of each arm and so affects the ability of the body for diffusive movements; also the variation of activator substance's concentration in the environment breaks the symmetry of linkers' preferred state, resulting in the drift of the random walker along the gradient of the density of activators. External food or danger sources may attract or repel the body by producing or consuming the chemical activators of the organism's enzymes, inducing chemotaxis behavior. Generalization of the model to three dimensions is straightforward.

Equilibrium state of a cylindrical particle with flat ends in nematic liquid crystals

A continuum theory is employed to numerically study the equilibrium orientation and defect structures of a circular cylindrical particle with flat ends under a homeotropic anchoring condition in a uniform nematic medium. Different aspect ratios of this colloidal geometry from thin discotic to long rodlike shapes and several colloidal length scales ranging from mesoscale to nanoscale are investigated. We show that the equilibrium state of this colloidal geometry is sensitive to the two geometrical parameters: aspect ratio and length scale of the particle. For a large enough mesoscopic particle, there is a specific asymptotic equilibrium angle associated to each aspect ratio. Upon reducing the particle size to nanoscale, the equilibrium angle follows a descending or ascending trend in such a way that the equilibrium angle of a particle with the aspect ratio bigger than 1:1 (a discotic particle) goes to a parallel alignment with respect to the far-field nematic, whereas the equilibrium angle for a particle with the aspect ratio 1:1 and smaller (a rodlike particle) tends toward a perpendicular alignment to the uniform nematic direction. The discrepancy between the equilibrium angles of the mesoscopic and nanoscopic particles originates from the significant differences between their defect structures. The possible defect structures related to mesoscopic and nanoscopic colloidal particles of this geometry are also introduced.

Particle selection through topographic templates in nematic colloids

Liquid crystal colloids have been proposed as suitable candidates for responsive photonic crystals. Large scale growth of such colloidal systems is, however, a challenge and recently template-assisted assembly has been proposed to guide the growth of colloidal crystals, with controlled symmetries, in nematic liquid crystals. Known for their long-range anisotropic interactions, these colloidal systems are stabilized typically at the center of the cells due to strong particle-wall repulsion from the confining substrates. This behaviour is dramatically changed in the presence of topographic patterning. Here we propose the use of topographic modulation of surfaces to select and localize particles in nematic colloids. By considering convex and concave deformations of one of the confining surfaces we show that the colloid-flat surface repulsion may be enhanced or switched into an attraction. In particular, we find that when the colloidal particles have the same anchoring conditions as the patterned surfaces, they are strongly attracted to concave dimples, while if they exhibit different anchoring conditions they are pinned at the top of convex protrusions. Although dominated by elastic interactions the first mechanism is reminiscent of the depletion induced attraction or of the key-lock mechanism, while the second is specific to liquid crystal colloids. These long-ranged, highly tunable, surface-colloid interactions contribute to the development of template-assisted assembly of large colloidal crystals, with well defined symmetries, as required for applications.

Confinement dynamics of a semiflexible chain inside nano-spheres

We study the conformations of a semiflexible chain, confined in nano-scaled spherical cavities, under two distinct processes of confinement. Radial contraction and packaging are employed as two confining procedures. The former method is performed by gradually decreasing the diameter of a spherical shell which envelopes a confined chain. The latter procedure is carried out by injecting the chain inside a spherical shell through a hole on the shell surface. The chain is modeled with a rigid body molecular dynamics simulation and its parameters are adjusted to DNA base-pair elasticity. Directional order parameter is employed to analyze and compare the confined chain and the conformations of the chain for two different sizes of the spheres are studied in both procedures. It is shown that for the confined chains in the sphere sizes of our study, they appear in spiral or tennis-ball structures, and the tennis-ball structure is more likely to be observed in more compact confinements. Our results also show that the dynamical procedure of confinement and the rate of the confinement are influential parameters of the structure of the chain inside spherical cavities.

Confined nematic liquid crystal between two spherical boundaries with planar anchoring

Nematic shells of liquid crystals have been provided in microscales. Defect structures in the shells are very essential in the electro-optical applications of such colloidal objects. We have numerically minimized the free energy of symmetric and asymmetric spherical shells of the nematic liquid crystal. Considering degenerate planar anchoring on the surfaces and isotropic nematic elasticity, a variety of defect structures are observed by controlling or varying the thicknesses of the shell and its degree of asymmetry. In symmetric shells, our calculations show that boojums (bipolar) defects appear in thick shells and tetrahedral (baseball) defects in thin shells. In asymmetric shells, while we are in the bipolar regime, the boojums defects transform to trigonal configurations. Free energy landscape shows that in this regime the inner droplet is not stable in the center and it is trapped in an off-center minimum energy position. For the case of thin shells, there are two degenerate director textures with similar tetrahedral configuration of the disclination lines. The levels are split in asymmetric shells. The stability of the inner droplet in the center position depends on director texture. It is stable for one texture and unstable for the other one. For an unstable pattern there is no minimum energy position for the inner droplet and it moves until it touches the outer boundary.

Rigid-body molecular dynamics of DNA inside a nucleosome

The majority of eukaryotic DNA, about three quarter, is wrapped around histone proteins forming so-called nucleosomes. To study nucleosomal DNA we introduce a coarse-grained molecular dynamics model based on sequence-dependent harmonic rigid base pair step parameters of DNA and nucleosomal binding sites. Mixed parametrization based on all-atom molecular dynamics and crystallographic data of protein-DNA structures is used for the base pair step parameters. The binding site parameters are adjusted by experimental B-factor values of the nucleosome crystal structure. The model is then used to determine the energy cost for placing a twist defect into the nucleosomal DNA which allows us to use Kramers theory to calculate nucleosome sliding caused by such defects. It is shown that the twist defect scenario together with the sequence-dependent elasticity of DNA can explain the slow time scales observed for nucleosome mobility along DNA. With this method we also show how the twist defect mechanism leads to a higher mobility of DNA in the presence of sin mutations near the dyad axis. Finally, by performing simulations on 5s rDNA, 601, and telomeric base pair sequences, it is demonstrated that the current model is a powerful tool to predict nucleosome positioning.

Definition of the persistence length in the coarse-grained models of DNA elasticity

By considering the detailed structure of DNA in the base pair level, two possible definitions of the persistence length are compared. One definition is related to the orientation of the terminal base pairs, and the other is based on the vectors which connect two adjacent base pairs at each end of the molecule. It is shown that although these definitions approach each other for long DNA molecules, they are dramatically different on short length scales. We show analytically that the difference mostly comes from the shear flexibility of the molecule and can be used to measure the shear modulus of DNA.

Gating and conduction of nano-channel forming proteins: a computational approach

Monitoring conformational changes in ion channels is essential to understand their gating mechanism. Here, we explore the structural dynamics of four outer membrane proteins with different structures and functions in the slowest nonzero modes of vibration. Normal mode analysis was performed on the modified elastic network model of channel in the membrane. According to our results, when membrane proteins were analyzed in the dominant mode, the composed pores, TolC and α-hemolysin showed large motions at the intramembrane β-barrel region while, in other porins, OmpA and OmpF, largest motions observed in the region of external flexible loops. A criterion based on equipartition theorem was used to measure the possible amplitude of vibration in channel forming proteins. The current approach complements theoretical and experimental techniques including HOLE, Molecular Dynamics (MD), and voltage clamp used to address the channel's structure and dynamics and provides the means to conduct a theoretical simultaneous study of the structure and function of the channel. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:3.