Per|Mut: Spatially Resolved Hydration Entropies from Atomistic Simulations.
Ontology highlight
ABSTRACT:
SUBMITTER: Heinz LP
PROVIDER: S-EPMC8047778 | biostudies-literature |
REPOSITORIES: biostudies-literature
Similar Datasets
Project description:Predicting efficacy and optimal drug delivery strategies for small molecule and biological therapeutics is challenging due to the complex interactions between diverse cell types in different tissues that determine disease outcome. Here we present a new methodology to simulate inflammatory disease manifestation and test potential intervention strategies in silico using agent-based computational models. Simulations created using this methodology have explicit spatial and temporal representations, and capture the heterogeneous and stochastic cellular behaviours that lead to emergence of pathology or disease resolution. To demonstrate this methodology we have simulated the prototypic murine T cell-mediated autoimmune disease experimental autoimmune encephalomyelitis, a mouse model of multiple sclerosis. In the simulation immune cell dynamics, neuronal damage and tissue specific pathology emerge, closely resembling behaviour found in the murine model. Using the calibrated simulation we have analysed how changes in the timing and efficacy of T cell receptor signalling inhibition leads to either disease exacerbation or resolution. The technology described is a powerful new method to understand cellular behaviours in complex inflammatory disease, permits rational design of drug interventional strategies and has provided new insights into the role of TCR signalling in autoimmune disease progression.
| S-EPMC3828403 | biostudies-literature
Project description:Spatially resolved transcriptomics (SRT)-specific computational methods are often developed, tested, validated, and evaluated in silico using simulated data. Unfortunately, existing simulated SRT data are often poorly documented, hard to reproduce, or unrealistic. Single-cell simulators are not directly applicable for SRT simulation as they cannot incorporate spatial information. We present SRTsim, an SRT-specific simulator for scalable, reproducible, and realistic SRT simulations. SRTsim not only maintains various expression characteristics of SRT data but also preserves spatial patterns. We illustrate the benefits of SRTsim in benchmarking methods for spatial clustering, spatial expression pattern detection, and cell-cell communication identification.
| S-EPMC9983268 | biostudies-literature
Project description:In this molecular dynamics simulation study, we analyze intermolecular vibrations in the hydration shell of a solvated enyzme, the membrane type 1-matrix metalloproteinase, with high spatial resolution. Our approach allows us to characterize vibrational signatures of the local hydrogen bond network, the translational mobility of water molecules, as well as the molecular entropy, in specific local environments. Our study demonstrates the heterogeneity of water properties within the hydration shell of a complex biomolecule. We define a classification scheme based on the vibrational density of states that allows us to distinguish separate classes of hydration water species and facilitates the description of hydration water properties at distinct hydration sites. The results demonstrate that no single characteristic of the protein surface is sufficient to determine the properties of nearby water. The protein surface geometry, quantified here by the number of protein atoms in the vicinity of a hydration water molecule, as well as the chemical nature of a solvated protein functional group, influences dynamic and thermodynamic properties of solvating water molecules.
| S-EPMC5607456 | 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:Using molecular dynamics free energy simulations with TIP3P explicit solvent, we compute the hydration free energies of 504 neutral small organic molecules and compare them to experiments. We find, first, good general agreement between the simulations and the experiments, with an RMS error of 1.24 kcal/mol over the whole set (i.e., about 2 kT) and a correlation coefficient of 0.89. Second, we use an automated procedure to identify systematic errors for some classes of compounds, and suggest some improvements to the force field. We find that alkyne hydration free energies are particularly poorly predicted due to problems with a Lennard-Jones well depth, and find that an alternate choice for this well depth largely rectifies the situation. Third, we study the non-polar component of hydration free energies - that is, the part that is not due to electrostatics. While we find that repulsive and attractive components of the non-polar part both scale roughly with surface area (or volume) of the solute, the total non-polar free energy does not scale with the solute surface area or volume, because it is a small difference between large components and is dominated by the deviations from the trend. While the methods used here are not new, this is a more extensive test than previous explicit solvent studies, and the size of the test set allows identification of systematic problems with force field parameters for particular classes of compounds. We believe that the computed free energies and components will be valuable to others in future development of force fields and solvation models.
| S-EPMC2701304 | biostudies-literature
Project description:We explore the spatial and temporal spread of the novel SARS-CoV-2 virus under containment measures in three European countries based on fits to data of the early outbreak. Using data from Spain and Italy, we estimate an age dependent infection fatality ratio for SARS-CoV-2, as well as risks of hospitalization and intensive care admission. We use them in a model that simulates the dynamics of the virus using an age structured, spatially detailed agent based approach, that explicitly incorporates governmental interventions and changes in mobility and contact patterns occurred during the COVID-19 outbreak in each country. Our simulations reproduce several of the features of its spatio-temporal spread in the three countries studied. They show that containment measures combined with high density are responsible for the containment of cases within densely populated areas, and that spread to less densely populated areas occurred during the late stages of the first wave. The capability to reproduce observed features of the spatio-temporal dynamics of SARS-CoV-2 makes this model a potential candidate for forecasting the dynamics of SARS-CoV-2 in other settings, and we recommend its application in low and lower-middle income countries which remain understudied.
| S-EPMC8323901 | 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:The celebrated Miller experiments reported on the spontaneous formation of amino acids from a mixture of simple molecules reacting under an electric discharge, giving birth to the research field of prebiotic chemistry. However, the chemical reactions involved in those experiments have never been studied at the atomic level. Here we report on, to our knowledge, the first ab initio computer simulations of Miller-like experiments in the condensed phase. Our study, based on the recent method of treatment of aqueous systems under electric fields and on metadynamics analysis of chemical reactions, shows that glycine spontaneously forms from mixtures of simple molecules once an electric field is switched on and identifies formic acid and formamide as key intermediate products of the early steps of the Miller reactions, and the crucible of formation of complex biological molecules.
| S-EPMC4183268 | biostudies-literature
Project description:In our recent efforts to map protein surfaces using mixed-solvent molecular dynamics (MixMD) (Ghanakota, P.; Carlson, H. A. Moving Beyond Active-Site Detection: MixMD Applied to Allosteric Systems. J. Phys. Chem. B 2016, 120, 8685-8695), we were able to successfully capture active sites and allosteric sites within the top-four most occupied hotspots. In this study, we describe our approach for estimating the thermodynamic profile of the binding sites identified by MixMD. First, we establish a framework for calculating free energies from MixMD simulations, and we compare our approach to alternative methods. Second, we present a means to obtain a relative ranking of the binding sites by their configurational entropy. The theoretical maximum and minimum free energy and entropy values achievable under such a framework along with the limitations of the techniques are discussed. Using this approach, the free energy and relative entropy ranking of the top-four MixMD binding sites were computed and analyzed across our allosteric protein targets: Abl Kinase, Androgen Receptor, Pdk1 Kinase, Farnesyl Pyrophosphate Synthase, Chk1 Kinase, Glucokinase, and Protein Tyrosine Phosphatase 1B.
| S-EPMC6549510 | biostudies-literature
Project description:Coarse-graining of fully atomistic molecular dynamics simulations is a long-standing goal in order to allow the description of processes occurring on biologically relevant timescales. For example, the prediction of pathways, rates and rate-limiting steps in protein-ligand unbinding is crucial for modern drug discovery. To achieve the enhanced sampling, we perform dissipation-corrected targeted molecular dynamics simulations, which yield free energy and friction profiles of molecular processes under consideration. Subsequently, we use these fields to perform temperature-boosted Langevin simulations which account for the desired kinetics occurring on multisecond timescales and beyond. Adopting the dissociation of solvated sodium chloride, trypsin-benzamidine and Hsp90-inhibitor protein-ligand complexes as test problems, we reproduce rates from molecular dynamics simulation and experiments within a factor of 2-20, and dissociation constants within a factor of 1-4. Analysis of friction profiles reveals that binding and unbinding dynamics are mediated by changes of the surrounding hydration shells in all investigated systems.
| S-EPMC7286908 | biostudies-literature