1. |
- Lyubartsev, Alexander, 1962-, et al.
(författare)
-
Systematic coarse-graining of molecular models by the Newton inversion method
- 2010
-
Ingår i: Faraday discussions. - 1359-6640 .- 1364-5498. ; 144, s. 43-56
-
Tidskriftsartikel (refereegranskat)abstract
- Systematic construction of coarse-grained molecular models from detailed atomistic simulations, and even from ab initio simulations is discussed. Atomistic simulations are first performed to extract structural information about the system, which is then used to determine effective potentials for a coarse-grained model of the same system. The statistical-mechanical equations expressing the canonical properties in terms of potential parameters can be inverted and solved numerically according to the iterative Newton scheme. In our previous applications, known as the Inverse Monte Carlo, radial distribution functions were inverted to reconstruct pair potential, while in a more general approach the targets can be other canonical averages. We have considered several examples of coarse-graining; for the united atom water model we suggest an easy way to overcome the known problem of high pressure. Further, we have developed coarse-grained models for L- and D-prolines, dissolved here in an organic solvent (dimethylsulfoxide), keeping their enantiomeric properties from the corresponding all-atom proline model. Finally, we have revisited the previously developed coarse-grained lipid model based on an updated all-atomic force field. We use this model in large-scale meso-scale simulations demonstrating spontaneous formation of different structures, such as vesicles, micelles, and multi-lamellar structures, depending on thermodynamical conditions.
|
|
2. |
- Mirzoev, Alexander, 1984-, et al.
(författare)
-
Systematic implicit solvent coarse graining of DMPC lipids
-
Tidskriftsartikel (refereegranskat)abstract
- A 10-site coarse-grained implicit solvent model of DMPC phospholipid is developed with effective solvent-mediated potentials derived using the inverse Monte Carlo method. The potentials providing for the coarse-grained model the same site-site RDFs, bond and angle distributions as those computed in atomistic simulations were computed for four different lipid-water molar ratios.It was shown that there is a significant concentration dependence for both effective potentials and properties of lipid bilayers simulated using these potentials. Thus, effective potentials computed at low lipid concentration favour to more condensed and ordered structure of the bilayer with lower average area per lipid, while potentials obtained at higher lipid concentrations provide more fluid-like structure. Introduction of bending angle interaction into coarse-grained model makes effective potentials somewhat less concentration-dependent. The best agreement with the reference data and experiment was achieved using the set of potentials derived from atomistic simulations at 1:30 lipid:water molar ratio providing full hydration of DMPC lipids in bilayer. Despite theoretical limitations of the structure-based approach expressed in state point dependence of the effective potentials, all the resulting potentials provide a stable bilayer structure with correct partitioning of different lipid groups across the bilayer as well as acceptable values of the average lipid area and compressibility. Another important property demonstrated by the model is liquid/gel phase transition observed at lowering the temperature. In addition to bilayer simulations, the model has proven its robustness in modeling of self-aggregation of lipids from randomly dispersed solution to ordered bilayer structures, bicells and vesicles.
|
|
3. |
- Santosh, Mysore S., et al.
(författare)
-
Molecular Dynamics Investigation of Dipeptide - Transition Metal Salts in Aqueous Solutions
- 2010
-
Ingår i: Journal of Physical Chemistry B. - : American Chemical Society (ACS). - 1520-6106 .- 1520-5207. ; 114:49, s. 16632-16640
-
Tidskriftsartikel (refereegranskat)abstract
- Molecular dynamics (MD) simulations of glycylglycine dipeptide with transition metal ions (Mn2+, Fe2+, Co2+, Ni2+, Cu2+, and Zn2+) in aqueous solutions have been carried out to get an insight into the solvation structure, intermolecular interactions, and salt effects in these systems. The solvation structure and hydrogen bonding were described in terms of radial distribution function (RDF) and spatial distribution function (SDF). The dynamical properties of the solvation structure were also analyzed in terms of diffusion and residence times. The simulation results show the presence of a well-defined first hydration shell around the dipeptide, with water molecules forming hydrogen bonds to the polar groups of the dipeptide. This shell is, however, affected by the strong electric field of divalent metal ions, which at higher ion concentrations lead to the shift in the dipeptide−water RDFs. Higher salt concentrations lead also to increased residence times and slower diffusion rates. In general, smaller ions (Cu2+, Zn2+) demonstrate stronger binding to dipeptide than the larger ones (Fe2+, Mn2+). Simulations do not show any stronger association of peptide molecules indicating their dissolution in water. The above results may be of potential interest to future researchers on these molecular interactions.
|
|
4. |
- Sun, Tiedong, et al.
(författare)
-
A Bottom-Up Coarse-Grained Model for Nucleosome-Nucleosome Interactions with Explicit Ions
- 2022
-
Ingår i: Journal of Chemical Theory and Computation. - : American Chemical Society (ACS). - 1549-9618 .- 1549-9626. ; 18:6, s. 3948-3960
-
Tidskriftsartikel (refereegranskat)abstract
- The nucleosome core particle (NCP) is a large complex of 145–147 base pairs of DNA and eight histone proteins and is the basic building block of chromatin that forms the chromosomes. Here, we develop a coarse-grained (CG) model of the NCP derived through a systematic bottom-up approach based on underlying all-atom MD simulations to compute the necessary CG interactions. The model produces excellent agreement with known structural features of the NCP and gives a realistic description of the nucleosome–nucleosome attraction in the presence of multivalent cations (Mg(H2O)62+ or Co(NH3)63+) for systems comprising 20 NCPs. The results of the simulations reveal structural details of the NCP–NCP interactions unavailable from experimental approaches, and this model opens the prospect for the rigorous modeling of chromatin fibers.
|
|
5. |
- Agosta, Lorenzo, et al.
(författare)
-
Adsorption of Glycine on TiO2 in Water from On-the-fly Free-Energy Calculations and In Situ Electrochemical Impedance Spectroscopy
- 2024
-
Ingår i: Langmuir. - 0743-7463 .- 1520-5827. ; 40:23, s. 12009-12016
-
Tidskriftsartikel (refereegranskat)abstract
- We report here an experimental-computational study of hydrated TiO2 anatase nanoparticles interacting with glycine, where we obtain quantitative agreement of the measured adsorption free energies. Ab initio simulations are performed within the tight binding and density functional theory in combination with enhanced free-energy sampling techniques, which exploit the thermodynamic integration of the unbiased mean forces collected on-the-fly along the molecular dynamics trajectories. The experiments adopt a new and efficient setup for electrochemical impedance spectroscopy measurements based on portable screen-printed gold electrodes, which allows fast and in situ signal assessment. The measured adsorption free energy is −30 kJ/mol (both from experiment and calculation), with preferential interaction of the charged group which strongly adsorbs on the TiO2 bridging oxygens. This highlights the importance of the terminal amino groups in the adsorption mechanism of amino acids on hydrated metal oxides. The excellent agreement between computation and experiment for this amino acid opens the doors to the exploration of the interaction free energies for other moderately complex bionano systems.
|
|
6. |
- B. Brant Carvalho, Paulo H., 1990-, et al.
(författare)
-
Neutron scattering study of polyamorphic THF·17(H2O) : toward a generalized picture of amorphous states and structures derived from clathrate hydrates
- 2023
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry. - 1463-9076 .- 1463-9084. ; 25:21, s. 14981-14991
-
Tidskriftsartikel (refereegranskat)abstract
- From crystalline tetrahydrofuran clathrate hydrate, THF-CH (THF·17H2O, cubic structure II), three distinct polyamorphs can be derived. First, THF-CH undergoes pressure-induced amorphization when pressurized to 1.3 GPa in the temperature range 77-140 K to a form which, in analogy to pure ice, may be called high-density amorphous (HDA). Second, HDA can be converted to a densified form, VHDA, upon heat-cycling at 1.8 GPa to 180 K. Decompression of VHDA to atmospheric pressure below 130 K produces the third form, recovered amorphous (RA). Results from neutron scattering experiments and molecular dynamics simulations provide a generalized picture of the structure of amorphous THF hydrates with respect to crystalline THF-CH and liquid THF·17H2O solution (∼2.5 M). Although fully amorphous, HDA is heterogeneous with two length scales for water-water correlations (less dense local water structure) and guest-water correlations (denser THF hydration structure). The hydration structure of THF is influenced by guest-host hydrogen bonding. THF molecules maintain a quasiregular array, reminiscent of the crystalline state, and their hydration structure (out to 5 Å) constitutes ∼23H2O. The local water structure in HDA is reminiscent of pure HDA-ice featuring 5-coordinated H2O. In VHDA, the hydration structure of HDA is maintained but the local water structure is densified and resembles pure VHDA-ice with 6-coordinated H2O. The hydration structure of THF in RA constitutes ∼18 H2O molecules and the water structure corresponds to a strictly 4-coordinated network, as in the liquid. Both VHDA and RA can be considered as homogeneous.
|
|
7. |
- Barros Brant Carvalho, Paulo Henrique, et al.
(författare)
-
Neutron scattering study of polyamorphic THF ∙ (H2O)17 – toward a generalized picture of amorphous states and structures derived from clathrate hydrates
-
Annan publikation (övrigt vetenskapligt/konstnärligt)abstract
- From crystalline tetrahydrofuran clathrate hydrate, THF-CH (THF ∙ 17H2O, cubic structure II), three distinct polyamorphs can be derived. First, THF-CH undergoes pressure-induced amorphization when pressurized to 1.3 GPa in the temperature range 77–140 K to a form which, in analogy to pure ice, may be called high-density amorphous (HDA). Second, HDA can be converted to a densified form, very-HDA (VHDA), upon heat-cycling at 1.8 GPa to 180 K. Decompression of VHDA to atmospheric pressure below 130 K produces the third, recovered amorphous (RA) form. Results from a compilation of neutron scattering experiments and molecular dynamics simulations provide a generalized picture of the structure of amorphous THF hydrates with respect to crystalline THF-CH and liquid THF ∙ 17H2O solution (~2.5 M). The calculated density of (only in situ observable) HDA and VHDA at 2 GPa and 130 K is 1.287 and 1.328 g/cm3, respectively, whereas that of RA (at 1 atm) is 1.081 g/cm3. Although fully amorphous, HDA is heterogeneous with two length scales for water-water correlations (less dense local water structure) and guest-water correlations (denser THF hydration structure). The hydration structure of THF is influenced by guest-host hydrogen bonding. THF molecules maintain a quasiregular array, reminiscent of the crystalline state, and their hydration structure (out to 5 Å) constitutes ~23 H2O. The local water structure in HDA is reminiscent of pure HDA-ice, featuring 5-coordinated H2O. In VHDA, this structure is maintained but the local water structure is densified to resemble pure VHDA-ice with 6-coordinated H2O. The hydration structure of THF in RA constitutes ~18 H2O and the water structure corresponds to a strictly 4-coordinated network, as in the liquid. Both VHDA and RA can be considered as homogeneous, solid solutions of THF and water. The local water structure of water-rich (1:17) amorphous CHs resembles most that of the corresponding amorphous water ices when compared to guest-rich CHs, e.g., Ar ∙ ~6H2O. The proposed significance of different contributions of water local environments presents a simple view to justify neutron structure factor features.
|
|
8. |
- Dong, Kun, et al.
(författare)
-
Structural Evidence for the Ordered Crystallites of Ionic Liquid in Confined Carbon Nanotubes
- 2009
-
Ingår i: The Journal of Physical Chemistry C. - : American Chemical Society (ACS). - 1932-7447 .- 1932-7455. ; 113:23, s. 10013-10020
-
Tidskriftsartikel (refereegranskat)abstract
- Ionic liquids (ILs) are a class of new green materials that have attracted extensive attention in recent decades. Many novel properties not evident under normal conditions may appear when ionic liquids are confined to a nanometer scale. As was observed in the experiment, an anomalous phase behavior from liquid to high melting point perfect crystal occurred when 1-n-butyl-3-methylimidazolium hexafluorophosphate ([bmim][PF6]) ionic liquid was confined in a carbon nanotube. In this work, we performed molecular dynamics (MD) simulations for [bmim][PF6] ionic liquid and provided direct structural evidence that the ionic crystallizes in a carbon nanotube. The ordered ionic arrangement in both the radial and the axial directions can be observed inside the channels of the CNTs to induce the form of crystallites. The ionic stacking and distributing can be determined by the sizes of the CNTs. Hydrogen bonds remain the dominant interactions between cations and anions when the ionic liquid enters into the CNT from the bulk phase. The free energies as the thermal driven forces were calculated, and it is found that it is very difficult for a single anion to enter into the channel of the CNT spontaneously. A more favorable way is through an ion-pair in which a cation “pulls” an anion to enter into the channel of the CNT together. It is predicted that other ionic liquids that possess similar structures, even including the pyridinium-based ionic liquids, can show higher melting points when confined in CNTs.
|
|
9. |
- Ermilova, Inna, 1983-, et al.
(författare)
-
Modelling of interactions between Aβ(25-35) peptide and phospholipid bilayers : effects of cholesterol and lipid saturation
- 2020
-
Ingår i: RSC Advances. - : Royal Society of Chemistry (RSC). - 2046-2069. ; 10:7, s. 3902-3915
-
Tidskriftsartikel (refereegranskat)abstract
- Aggregation of amyloid beta (Aβ) peptides in neuronal membranes is a known promoter of Alzheimer’s disease. To gain insight into the molecular details of Aβ peptide aggregation and its effect on model neuronal membranes, we carried out molecular dynamics simulations of the Aβ(25–35) fragment of the amyloid precursor protein in phospholipid bilayers composed of either fully saturated or highly unsaturated lipids, in the presence or absence of cholesterol. It was found that the peptide does not penetrate through any of the considered membranes, but can reside in the headgroup region and upper part of the lipid tails showing a clear preference to a polyunsaturated cholesterol-free membrane. Due to the ordering and condensing effect upon addition of cholesterol, membranes become more rigid facilitating peptide aggregation on the surface. Except for the case of the cholesterol-free saturated lipid bilayer, the peptides have a small effect on the membrane structure and ordering. It was also found that the most “active” amino-acid for peptide–lipid and peptide–cholesterol interaction is methionine-35, followed by asparagine-27 and serine-26, which form hydrogen bonds between peptides and polar atoms of lipid headgroups. These amino acids are also primarily responsible for peptide aggregation. This work will be relevant for designing strategies to develop drugs to combat Alzheimer’s disease.
|
|
10. |
- Grote, Fredrik, 1995-, et al.
(författare)
-
Phase equilibrium, dynamics and rheology of phospholipid-ethanol mixtures : a combined molecular dynamics, NMR and viscometry study
- 2023
-
Ingår i: Physical Chemistry, Chemical Physics - PCCP. - : Royal Society of Chemistry (RSC). - 1463-9076 .- 1463-9084. ; 25:23, s. 15905-15915
-
Tidskriftsartikel (refereegranskat)abstract
- Binary mixtures of ethanol and phospholipids DOPC and DOPE have been investigated in a composition range relevant for topical drug delivery applications. This was done using a combined computer simulation and experimental approach where molecular dynamics simulations of ethanol-lipid mixtures with different compositions were performed. Several key properties including diffusion coefficients, longitudinal relaxation times, and shear viscosity were computed. In addition, diffusion coefficients, viscosities and NMR longitudinal relaxation times were measured experimentally for comparison and in order to validate the results from simulation. Diffusion coefficients and relaxation times obtained from simulations are in good agreement with results from NMR and computed viscosities are in reasonable agreement with viscometry experiments indicating that the simulations provide a realistic description of the ethanol-phospholipid mixtures. Structural changes in the simulated systems were investigated using an analysis based on radial distribution functions. This showed that the structure of ethanol-DOPC mixtures remains essentially unchanged in the investigated concentration range while ethanol-DOPE mixtures undergo structural rearrangements with the tendency for forming small aggregates on the 100 ns time scale consisting of less than 10 lipids. Although our simulations and experiments indicate that no larger aggregates form, they also show that DOPE has stronger aggregation tendency than DOPC. This highlights the importance of the character of the lipid headgroup for lipid aggregation in ethanol and gives new insights into phase equilibrium, dynamics and rheology that could be valuable for the development of advanced topical drug delivery formulations.
|
|