Engineering of 2D nanomaterials to trap and kill SARS-CoV-2: a new insight from multi-microsecond atomistic simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Khedri M
PROVIDER: S-EPMC8413075 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:The formation of membraneless organelles in cells commonly occurs via liquid-liquid phase separation (LLPS) and is in many cases driven by multivalent interactions between intrinsically disordered proteins (IDPs). Investigating the nature of these interactions, and their effect on dynamics within the condensed phase, is therefore of critical importance but very challenging for either simulation or experiment. Here, we study these interactions and their dynamics by pairing a novel multiscale simulation strategy with microsecond all-atom MD simulations of a condensed, IDP-rich phase. We simulate two IDPs this way, the low complexity domain of FUS and the N-terminal disordered domain of LAF-1, and find good agreement with experimental information about average density, water content, and residue-residue contacts. We go significantly beyond what is known from experiments by showing that ion partitioning within the condensed phase is largely driven by the charge distribution of the proteins and-in the cases considered-shows little evidence of preferential interactions of the ions with the proteins. Furthermore, we can probe the microscopic diffusive dynamics within the condensed phase, showing that water and ions are in dynamic equilibrium between dense and dilute phases, and their diffusion is reduced in the dense phase. Despite their high concentration in the condensate, the protein molecules also remain mobile, explaining the observed liquid-like properties of this phase. We finally show that IDP self-association is driven by a combination of nonspecific hydrophobic interactions as well as hydrogen bonds, salt bridges, and π-π and cation-π interactions. The simulation approach presented here allows the structural and dynamical properties of biomolecular condensates to be studied in microscopic detail and is generally applicable to single- and multicomponent systems of proteins and nucleic acids involved in LLPS.
| S-EPMC7879053 | biostudies-literature
Project description:The recent pandemic of COVID-19 has raised global health concerns. Preventing severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) activity in the body is a very promising method to overcome the COVID-19 pandemic. One of the prevention methods is constraining the binding process among the human cell receptor-ACE2 and coronavirus spike protein. In the research done, the effect of deformation of the spike protein structure, due to the covalent organic frameworks (COFs), in reducing the interactions of ACE2 and the spike protein by the computational method was investigated. In this regard, atomic analysis of the interactions of ACE2 and the spike protein is provided using a molecular dynamics simulation. First, we investigated the interactions of the three different COFs, including COF-78, DAAQ-TFP, and COF-OEt, with the spike protein by analyzing the bond energies, as well as structural changes of the spike protein. Then, intermolecular interactions of the deformed spike protein along with ACE2 were assessed to clarify the protein's fusion after the deformation. As indicated by the results, although all introduced COFs deformed the spike protein in an effective way, COF-78 showed the best performance in the prevention of spike protein-ACE2 interactions by changing the molecular structure of the protein. Indeed, the interaction analysis of the deformed spike protein by COF-78 with the ACE2 showed that their interactions had the lowest absolute value of energy, along with the least amount of hydrogen bonds, in which the compaction of the protein was lower compared to the other deformed proteins. Moreover, having a high contact area with an aqueous media as well as severe fluctuations during the simulation time confirmed the positive performance of COF-78. In the current study, we aimed to introduce novel materials and COVID-19 prevention methodology that can be used in face masks and for surface disinfection.
| S-EPMC8550915 | biostudies-literature
Project description:Cholesterol is widely known to alter the physical properties and permeability of membranes. Several prior works have implicated cell membrane cholesterol as a barrier to tissue oxygenation, yet a good deal remains to be explained with regard to the mechanism and magnitude of the effect. We use molecular dynamics simulations to provide atomic-resolution insight into the influence of cholesterol on oxygen diffusion across and within the membrane. Our simulations show strong overall agreement with published experimental data, reproducing the shapes of experimental oximetry curves with high accuracy. We calculate the upper-limit transmembrane oxygen permeability of a 1-palmitoyl,2-oleoylphosphatidylcholine phospholipid bilayer to be 52 ± 2 cm/s, close to the permeability of a water layer of the same thickness. With addition of cholesterol, the permeability decreases somewhat, reaching 40 ± 2 cm/s at the near-saturating level of 62.5 mol % cholesterol and 10 ± 2 cm/s in a 100% cholesterol mimic of the experimentally observed noncrystalline cholesterol bilayer domain. These reductions in permeability can only be biologically consequential in contexts where the diffusional path of oxygen is not water dominated. In our simulations, cholesterol reduces the overall solubility of oxygen within the membrane but enhances the oxygen transport parameter (solubility-diffusion product) near the membrane center. Given relatively low barriers to passing from membrane to membrane, our findings support hydrophobic channeling within membranes as a means of cellular and tissue-level oxygen transport. In such a membrane-dominated diffusional scheme, the influence of cholesterol on oxygen permeability is large enough to warrant further attention.
| S-EPMC5474842 | biostudies-other
Project description:Developing strategies for atomic-scale controlled synthesis of new two-dimensional (2D) functional materials will directly impact their applications. Here, using in situ aberration-corrected scanning transmission electron microscopy, we obtain direct insight into the homoepitaxial Frank-van der Merwe atomic layer growth mechanism of TiC single adlayers synthesized on surfaces of Ti3C2 MXene substrates with the substrate being the source material. Activated by thermal exposure and electron-beam irradiation, hexagonal TiC single adlayers form on defunctionalized surfaces of Ti3C2 MXene at temperatures above 500 °C, generating new 2D materials Ti4C3 and Ti5C4. The growth mechanism for a single TiC adlayer and the energies that govern atom migration and diffusion are elucidated by comprehensive density functional theory and force-bias Monte Carlo/molecular dynamics simulations. This work could lead to the development of bottom-up synthesis methods using substrates terminated with similar hexagonal-metal surfaces, for controllable synthesis of larger-scale and higher quality single-layer transition metal carbides.
| S-EPMC5995835 | biostudies-literature
Project description:Despite the development of massively parallel computing hardware including inexpensive graphics processing units (GPUs), it has remained infeasible to simulate the folding of atomistic proteins at room temperature using conventional molecular dynamics (MD) beyond the microsecond scale. Here, we report the folding of atomistic, implicitly solvated protein systems with folding times τ ranging from ∼10 μs to ∼100 ms using the weighted ensemble (WE) strategy in combination with GPU computing. Starting from an initial structure or set of structures, WE organizes an ensemble of GPU-accelerated MD trajectory segments via intermittent pruning and replication events to generate statistically unbiased estimates of rate constants for rare events such as folding; no biasing forces are used. Although the variance among atomistic WE folding runs is significant, multiple independent runs are used to reduce and quantify statistical uncertainty. Folding times are estimated directly from WE probability flux and from history-augmented Markov analysis of the WE data. Three systems were examined: NTL9 at low solvent viscosity (yielding τf = 0.8-9 μs), NTL9 at water-like viscosity (τf = 0.2-2 ms), and Protein G at low viscosity (τf = 3-200 ms). In all cases, the folding time, uncertainty, and ensemble properties could be estimated from WE simulation; for Protein G, this characterization required significantly less overall computing than would be required to observe a single folding event with conventional MD simulations. Our results suggest that the use and calibration of force fields and solvent models for precise estimation of kinetic quantities is becoming feasible.
| S-EPMC6660137 | biostudies-literature
Project description:Mixtures of long- and short-tail phosphatidylcholine lipids are known to self-assemble into a variety of aggregates combining flat bilayerlike and curved micellelike features, commonly called bicelles. Atomistic simulations of bilayer ribbons and perforated bilayers containing dimyristoylphosphatidylcholine (DMPC, di-C(14) tails) and dihexanoylphosphatidylcholine (DHPC, di-C(6) tails) have been carried out to investigate the partitioning of these components between flat and curved microenvironments and the stabilization of the bilayer edge by DHPC. To approach equilibrium partitioning of lipids on an achievable simulation timescale, configuration-bias Monte Carlo mutation moves were used to allow individual lipids to change tail length within a semigrand-canonical ensemble. Since acceptance probabilities for direct transitions between DMPC and DHPC were negligible, a third component with intermediate tail length (didecanoylphosphatidylcholine, di-C(10) tails) was included at a low concentration to serve as an intermediate for transitions between DMPC and DHPC. Strong enrichment of DHPC is seen at ribbon and pore edges, with an excess linear density of approximately 3 nm(-1). The simulation model yields estimates for the onset of edge stability with increasing bilayer DHPC content between 5% and 15% DHPC at 300 K and between 7% and 17% DHPC at 323 K, higher than experimental estimates. Local structure and composition at points of close contact between pores suggest a possible mechanism for effective attractions between pores, providing a rationalization for the tendency of bicelle mixtures to aggregate into perforated vesicles and perforated sheets.
| S-EPMC2884235 | biostudies-other
Project description:Here, molecular dynamics (MD) simulations were employed to explore the self-assembly of polymers and docetaxel (DTX) as an anticancer drug in the presence of nitrogen, phosphorous, and boron-nitrogen incorporated graphene and fullerene. The electrostatic potential and the Gibbs free energy of the self-assembled materials were used to optimize the atomic doping percentage of the N- and P-doped formulations at 10% and 50%, respectively. Poly lactic-glycolic acid (PLGA)- polyethylene glycol (PEG)-based polymeric nanoparticles were assembled in the presence of nanocarbons in the common (corresponding to the bulk environment) and interface of organic/aqueous solutions (corresponding to the microfluidic environment). Assessment of the modeling results (e.g., size, hydrophobicity, and energy) indicated that among the nanocarbons, the N-doped graphene nanosheet in the interface method created more stable polymeric nanoparticles (PNPs). Energy analysis demonstrated that doping with nanocarbons increased the electrostatic interaction energy in the self-assembly process. On the other hand, the fullerene-based nanocarbons promoted van der Waals intramolecular interactions in the PNPs. Next, the selected N-doped graphene nanosheet was utilized to prepare nanoparticles and explore the physicochemical properties of the nanosheets in the permeation of the resultant nanoparticles through cell-based lipid bilayer membranes. In agreement with the previous results, the N-graphene assisted PNP in the interface method and was translocated into and through the cell membrane with more stable interactions. In summary, the present MD simulation results demonstrated the success of 2D graphene dopants in the nucleation and growth of PLGA-based nanoparticles for improving anticancer drug delivery to cells, establishing new promising materials and a way to assess their performance that should be further studied.
| S-EPMC8564517 | biostudies-literature
Project description:The dynamics of polymer chains in the polymer/solid interphase region have been a point of debate in recent years. Its understanding is the first step towards the description and the prediction of the properties of a wide family of commercially used polymeric-based nanostructured materials. Here, we present a detailed investigation of the conformational and dynamical features of unentangled and mildly entangled cis-1,4-polybutadiene melts in the vicinity of amorphous silica surface via atomistic simulations. Accounting for the roughness of the surface, we analyze the properties of the polymer chains as a function of their distance from the silica slab, their conformations and the chain molecular weight. Unlike the case of perfectly flat and homogeneous surfaces, the monomeric translational motion parallel to the surface was affected by the presence of the silica slab up to distances comparable with the extension of the density fluctuations. In addition, the intramolecular dynamical heterogeneities in adsorbed chains were revealed by linking the conformations and the structure of the adsorbed chains with their dynamical properties. Strong dynamical heterogeneities within the adsorbed layer are found, with the chains possessing longer sequences of adsorbed segments ("trains") exhibiting slower dynamics than the adsorbed chains with short ones. Our results suggest that, apart from the density-dynamics correlation, the configurational entropy plays an important role in the dynamical response of the polymers confined between the silica slabs.
| S-EPMC8401719 | biostudies-literature
Project description:The Calcium Silicate Hydrate (C-S-H) nucleation is a crucial step during cement hydration and determines to a great extent the rheology, microstructure, and properties of the cement paste. Recent evidence indicates that the C-S-H nucleation involves at least two steps, yet the underlying atomic scale mechanism, the nature of the primary particles and their stability, or how they merge/aggregate to form larger structures is unknown. In this work, we use atomistic simulation methods, specifically DFT, evolutionary algorithms (EA), and Molecular Dynamics (MD), to investigate the structure and formation of C-S-H primary particles (PPs) from the ions in solution, and then discuss a possible formation pathway for the C-S-H nucleation. Our simulations indicate that even for small sizes the most stable clusters encode C-S-H structural motifs, and we identified a C4S4H2 cluster candidate to be the C-S-H basic building block. We suggest a formation path in which small clusters formed by silicate dimers merge into large elongated aggregates. Upon dehydration, the C-S-H basic building blocks can be formed within the aggregates, and eventually crystallize.
| S-EPMC10693585 | biostudies-literature
Project description:Molecular dynamics (MD) simulations of membranes are often hindered by the slow lateral diffusion of lipids and the limited time scale of MD. In order to study the dynamics of mixing and characterize the lateral distribution of lipids in converged mixtures, we report microsecond-long all-atom MD simulations performed on the special-purpose machine Anton. Two types of mixed bilayers, POPE:POPG (3:1) and POPC:cholesterol (2:1), as well as a pure POPC bilayer, were each simulated for up to 2 μs. These simulations show that POPE:POPG and POPC:cholesterol are each fully miscible at the simulated conditions, with the final states of the mixed bilayers similar to a random mixture. By simulating three POPE:POPG bilayers at different NaCl concentrations (0, 0.15, and 1 M), we also examined the effect of salt concentration on lipid mixing. While an increase in NaCl concentration is shown to affect the area per lipid, tail order, and lipid lateral diffusion, the final states of mixing remain unaltered, which is explained by the largely uniform increase in Na(+) ions around POPE and POPG. Direct measurement of water permeation reveals that the POPE:POPG bilayer with 1 M NaCl has reduced water permeability compared with those at zero or low salt concentration. Our calculations provide a benchmark to estimate the convergence time scale of all-atom MD simulations of lipid mixing. Additionally, equilibrated structures of POPE:POPG and POPC:cholesterol, which are frequently used to mimic bacterial and mammalian membranes, respectively, can be used as starting points of simulations involving these membranes.
| S-EPMC4196744 | biostudies-literature