polarizable force fields
Dotaz
Zobrazit nápovědu
Recent advances in polarizable force fields have revealed that major reparameterization is necessary when the polarization energy is treated explicitly. This study is focused on the torsional parameters, which are crucial for the accurate description of conformational equilibria in biomolecules. In particular, attention is paid to the influence of polarization on the (i) transferability of dihedral terms between molecules, (ii) transferability between different environments, and (iii) additivity of dihedral energies. To this end, three polarizable force fields based on the induced point dipole model designed for use in AMBER are tested, including two recent ff02 reparameterizations. Attention is paid to the contributions due to short range interactions (1-2, 1-3, and 1-4) within the four atoms defining the dihedral angle. The results show that when short range 1-2 and 1-3 polarization interactions are omitted, as for instance in ff02, the 1-4 polarization contribution is rather small and unlikely to improve the description of the torsional energy. Conversely, when screened 1-2 and 1-3 interactions are included, the polarization contribution is sizeable and shows potential to improve the transferability of parameters between different molecules and environments as well as the additivity of dihedral terms. However, to reproduce intramolecular polarization effects accurately, further fine-tuning of the short range damping of polarization is necessary.
The coarse-grained Martini force field is widely used in biomolecular simulations. Here we present the refined model, Martini 3 ( http://cgmartini.nl ), with an improved interaction balance, new bead types and expanded ability to include specific interactions representing, for example, hydrogen bonding and electronic polarizability. The updated model allows more accurate predictions of molecular packing and interactions in general, which is exemplified with a vast and diverse set of applications, ranging from oil/water partitioning and miscibility data to complex molecular systems, involving protein-protein and protein-lipid interactions and material science applications as ionic liquids and aedamers.
... The Four Forces of Nature 3 -- 1.2. Greek and Medieval Notions of Intermolecular Forces 3 -- 1.3. ... ... Intermolecular Force-Laws and Interaction -- Potentials: Long- and Short-Range Forces 9 -- 1.7. ... ... Solvation Forces, Structural Forces, and Hydration Forces 80 -- 4.7. ... ... “Body Forces” and “Surface Forces” 220 -- Problems and Discussion Topics 220 -- 12. ... ... Attractive Depletion Forces 593 -- 21.5. Attractive Hydrophobic Forces 595 -- 21.6. ...
Third edition xxx, 674 stran : 24 cm il. ;
- MeSH
- fyzikální chemie MeSH
- Publikační typ
- monografie MeSH
- Konspekt
- Fyzikální chemie
- NLK Obory
- chemie, klinická chemie
... electrical forces 4 -- 1.3 The electric field 5 -- 1.4 Representation of electric fields 8 -- 1.5 Gauss ... ... How conductors shape an electric field 29 -- 2.1 Electrical conductors 29 -- 2.2 The electric field within ... ... Static magnetic fields 100 -- 6.1 The magnetic field and Gauss’s law 100 -- 6.2 The interaction of a ... ... magnetic dipole and a magnetic field 104 -- 6.3 The magnetic field of a macroscopic current loop 106 ... ... force 158 -- 10.2 The motion of a free charged particle in a static magnetic field 159 x Contents -- ...
xii, 286 s.
Eigenmodes of the spherical and ellipsoidal dielectric electromagnetic resonator have been analysed. The sizes and shape of the resonators have been chosen to represent the shape of the interphase and dividing animal cell. Electromagnetic modes that have shape exactly suitable for positioning of the sufficiently large organelles in cell (centrosome, nucleus) have been identified. We analysed direction and magnitude of dielectrophoretic force exerted on large organelles by electric field of the modes. We found that the TM(1m1) mode in spherical resonator acts by centripetal force which drags the large organelles which have higher permittivity than the cytosol to the center of the cell. TM-kind of mode in the ellipsoidal resonator acts by force on large polarizable organelles in a direction that corresponds to the movement of the centrosomes (also nucleus) observed during the cell division, i.e. to the foci of the ellipsoidal cell. Minimal required force (10(-16) N), gradient of squared electric field and corresponding energy (10(-16) J) of the mode have been calculated to have biological significance within the periods on the order of time required for cell division. Minimal required energy of the mode, in order to have biological significance, can be lower in the case of resonance of organelle with the field of the cellular resonator mode. In case of sufficient energy in the biologically relevant mode, electromagnetic field of the mode will act as a positioning or steering mechanism for centrosome and nucleus in the cell, thus contribute to the spatial and dynamical self-organization in biological systems.
Fluorescence solvent relaxation experiments are based on the characterization of time-dependent shifts in the fluorescence emission of a chromophore, yielding polarity and viscosity information about the chromophore's immediate environment. A chromophore applied to a phospholipid bilayer at a well-defined location (with respect to the z-axis of the bilayer) allows monitoring of the hydration and mobility of the probed segment of the lipid molecules. Specifically, time-resolved fluorescence experiments, fluorescence quenching data and molecular dynamic (MD) simulations show that 6-lauroyl-2-dimethylaminonaphthalene (Laurdan) probes the hydration and mobility of the sn-1 acyl groups in a phosphatidylcholine bilayer. The time-dependent fluorescence shift (TDFS) of Laurdan provides information on headgroup compression and expansion induced by the addition of different amounts of cationic lipids to phosphatidylcholine bilayers. Those changes were predicted by previous MD simulations. Addition of truncated oxidized phospholipids leads to increased mobility and hydration at the sn-1 acyl level. This experimental finding can be explained by MD simulations, which indicate that the truncated chains of the oxidized lipid molecules are looping back into aqueous phase, hence creating voids below the glycerol level. Fluorescence solvent relaxation experiments are also useful in understanding salt effects on the structure and dynamics of lipid bilayers. For example, such experiments demonstrate that large anions increase hydration and mobility at the sn-1 acyl level of phosphatidylcholine bilayers, an observation which could not be explained by standard MD simulations. If polarizability is introduced into the applied force field, however, MD simulations show that big soft polarizable anions are able to interact with the hydrophilic/hydrophobic interface of the lipid bilayer, penetrating to the level probed by Laurdan, and that they expand and destabilize the bilayer making it more hydrated and mobile.
The lone-pair···π (lp···π) (deoxy)ribose···nucleobase stacking is a recurring interaction in Z-DNA and RNAs that is characterized by sub-van der Waals lp···π contacts (<3.0 Å). It is a part of the structural signature of CpG Z-step motifs in Z-DNA and r(UNCG) tetraloops that are known to behave poorly in molecular dynamics (MD) simulations. Although the exact origin of the MD simulation issues remains unclear, a significant part of the problem might be due to an imbalanced description of nonbonded interactions, including the characteristic lp···π stacking. To gain insights into the links between lp···π stacking and MD, we present an in-depth comparison between accurate large-basis-set double-hybrid Kohn-Sham density functional theory calculations DSD-BLYP-D3/ma-def2-QZVPP (DHDF-D3) and data obtained with the nonbonded potential of the AMBER force field (AFF) for NpN Z-steps (N = G, A, C, and U). Among other differences, we found that the AFF overestimates the DHDF-D3 lp···π distances by ∼0.1-0.2 Å, while the deviation between the DHDF-D3 and AFF descriptions sharply increases in the short-range region of the interaction. Based on atom-in-molecule polarizabilities and symmetry-adapted perturbation theory analysis, we inferred that the DHDF-D3 versus AFF differences partly originate in identical nucleobase carbon atom Lennard-Jones (LJ) parameters despite the presence/absence of connected electron-withdrawing groups that lead to different effective volumes or vdW radii. Thus, to precisely model the very short CpG lp···π contact distances, we recommend revision of the nucleobase atom LJ parameters. Additionally, we suggest that the large discrepancy between DHDF-D3 and AFF short-range repulsive part of the interaction energy potential may significantly contribute to the poor performances of MD simulations of nucleic acid systems containing Z-steps. Understanding where, and if possible why, the point-charge-type effective potentials reach their limits is vital for developing next-generation FFs and for addressing specific issues in contemporary MD simulations.
