Influence of Flexible Side-Chains on the Breathing Phase Transition of Pillared Layer MOFs: A Force Field Investigation

2020 ◽  
Author(s):  
Julian Keupp ◽  
Johannes P. Dürholt ◽  
Rochus Schmid
Keyword(s):  
First Principles ◽  
Good Accuracy ◽  
Field Investigation ◽  
Square Lattice ◽  
Side Chain ◽  
Second Step ◽  
Side Chains ◽  
Dynamics Simulations ◽  
Complex Phenomena ◽  
Closed Pore

The prototypical pillared layer MOFs, formed by a square lattice of paddle-<br>wheel units and connected by dinitrogen pillars, can undergo a breathing phase<br>transition by a “wine-rack” type motion of the square lattice. We studied this not<br>yet fully understood behavior using an accurate first principles parameterized force<br>field (MOF-FF) for larger nanocrystallites on the example of Zn 2 (bdc) 2 (dabco) [bdc:<br>benzenedicarboxylate, dabco: (1,4-diazabicyclo[2.2.2]octane)] and found clear indi-<br>cations for an interface between a closed and an open pore phase traveling through<br>the system during the phase transformation [Adv. Theory Simul. 2019, 2, 11]. In<br>conventional simulations in small supercells this mechanism is prevented by periodic<br>boundary conditions (PBC), enforcing a synchronous transformation of the entire<br>crystal. Here, we extend this investigation to pillared layer MOFs with flexible<br>side-chains, attached to the linker. Such functionalized (fu-)MOFs are experimen-<br>tally known to have different properties with the side-chains acting as fixed guest<br>molecules. First, in order to extend the parameterization for such flexible groups,<br>1a new parametrization strategy for MOF-FF had to be developed, using a multi-<br>structure force based fit method. The resulting parametrization for a library of<br>fu-MOFs is then validated with respect to a set of reference systems and shows very<br>good accuracy. In the second step, a series of fu-MOFs with increasing side-chain<br>length is studied with respect to the influence of the side-chains on the breathing<br>behavior. For small supercells in PBC a systematic trend of the closed pore volume<br>with the chain length is observed. However, for a nanocrystallite model a distinct<br>interface between a closed and an open pore phase is visible only for the short chain<br>length, whereas for longer chains the interface broadens and a nearly concerted trans-<br>formation is observed. Only by molecular dynamics simulations using accurate force<br>fields such complex phenomena can be studied on a molecular level.

Influence of Flexible Side-Chains on the Breathing Phase Transition of Pillared Layer MOFs: A Force Field Investigation

2020 ◽  
Author(s):  
Julian Keupp ◽  
Johannes P. Dürholt ◽  
Rochus Schmid
Keyword(s):  
First Principles ◽  
Good Accuracy ◽  
Field Investigation ◽  
Square Lattice ◽  
Side Chain ◽  
Second Step ◽  
Side Chains ◽  
Dynamics Simulations ◽  
Complex Phenomena ◽  
Closed Pore

The prototypical pillared layer MOFs, formed by a square lattice of paddle-<br>wheel units and connected by dinitrogen pillars, can undergo a breathing phase<br>transition by a “wine-rack” type motion of the square lattice. We studied this not<br>yet fully understood behavior using an accurate first principles parameterized force<br>field (MOF-FF) for larger nanocrystallites on the example of Zn 2 (bdc) 2 (dabco) [bdc:<br>benzenedicarboxylate, dabco: (1,4-diazabicyclo[2.2.2]octane)] and found clear indi-<br>cations for an interface between a closed and an open pore phase traveling through<br>the system during the phase transformation [Adv. Theory Simul. 2019, 2, 11]. In<br>conventional simulations in small supercells this mechanism is prevented by periodic<br>boundary conditions (PBC), enforcing a synchronous transformation of the entire<br>crystal. Here, we extend this investigation to pillared layer MOFs with flexible<br>side-chains, attached to the linker. Such functionalized (fu-)MOFs are experimen-<br>tally known to have different properties with the side-chains acting as fixed guest<br>molecules. First, in order to extend the parameterization for such flexible groups,<br>1a new parametrization strategy for MOF-FF had to be developed, using a multi-<br>structure force based fit method. The resulting parametrization for a library of<br>fu-MOFs is then validated with respect to a set of reference systems and shows very<br>good accuracy. In the second step, a series of fu-MOFs with increasing side-chain<br>length is studied with respect to the influence of the side-chains on the breathing<br>behavior. For small supercells in PBC a systematic trend of the closed pore volume<br>with the chain length is observed. However, for a nanocrystallite model a distinct<br>interface between a closed and an open pore phase is visible only for the short chain<br>length, whereas for longer chains the interface broadens and a nearly concerted trans-<br>formation is observed. Only by molecular dynamics simulations using accurate force<br>fields such complex phenomena can be studied on a molecular level.

scholarly journals The role of side chains in the interaction of new antitumor pyrimidoacridinetriones with DNA: molecular dynamics simulations.

2000 ◽  
Vol 47 (1) ◽  
pp. 47-57 ◽  
Author(s):  
J Mazerski ◽  
I Antonini ◽  
S Martelli
Keyword(s):  
Molecular Dynamics ◽  
Rational Design ◽  
Md Simulations ◽  
Ring System ◽  
Side Chain ◽  
Side Chains ◽  
Intercalation Complex ◽  
Highly Active ◽  
Simulation Techniques ◽  
Dynamics Simulations

Pyrimidoacridinetriones (PATs) are a new group of highly active antitumor compounds. It seems reasonable to assume that, like for some other acridine derivatives, intercalation into DNA is a necessary, however not a sufficient condition for antitumor activity of these compounds. Rational design of new compounds of this chemotype requires knowledge about the structure of the intercalation complex, as well as about interactions responsible for its stability. Computer simulation techniques such as molecular dynamics (MD) may provide valuable information about these problems. The results of MD simulations performed for three rationally selected PATs are presented in this paper. The compounds differ in the number and position of side chains. Each of the compounds was simulated in two systems: i) in water, and ii) in the intercalation complex with the dodecamer duplex d(GCGCGCGCGCGC)2. The orientation of the side chain in relation to the ring system is determined by the position of its attachment. Orientation of the ring system inside the intercalation cavity depends on the number and position of side chain(s). The conformations of the side chain(s) of all PATs studied in the intercalation complex were found to be very similar to those observed in water.

scholarly journals Side Chains and the Insufficient Lubrication of Water in Polyacrylamide Hydrogel—A New Insight

Polymers ◽  
2019 ◽  
Vol 11 (11) ◽  
pp. 1845 ◽  
Author(s):  
Lei ◽  
Zhou ◽  
Liu
Keyword(s):  
Mechanical Properties ◽  
Volume Fraction ◽  
Polymer Network ◽  
Polymer Chains ◽  
Side Chain ◽  
Side Chains ◽  
Polyacrylamide Hydrogels ◽  
Young’S Moduli ◽  
Property Changes ◽  
Dynamics Simulations

Existing theories cannot predict the mechanical property changes of polyacrylamide hydrogels with different water content because of the absence of side chains. In this study, polyacrylamide hydrogels are prepared and tested to investigate the side chain effect on their mechanical properties. First, the comparison between the effective chain density and total chain density provides proof of the large amount of side chains in the polymer network of PAAm hydrogel. We propose a practical chain density fraction to measure the side chain fraction. Then, the abnormal Young’s moduli-polymer volume fraction relationship reveals that side chains affect the mechanical properties of hydrogel through the insufficient lubrication of water. Water confined in narrow space within a molecular-level size can bear shear force to provide extra deformation resistance. A constitutive mode considering the effect of the insufficient lubrication of water is proposed. Combining this constitutive model with experimental results, we find that this insufficient lubrication of water exists even in equilibrium PAAm hydrogel. Molecular dynamics simulations reveal that this insufficient lubrication of water comes from the constraint of polymer chains. It also demonstrates that when there is insufficient lubrication, the rearrangement of water molecules leads to the persistent energy dissipation in the Mullins effect of PAAm hydrogel.

Modeling the Effects of O-sulfonation on the CID of Serine

2020 ◽  
Author(s):  
Kenneth Lucas ◽  
George Barnes
Keyword(s):  
Sulfo Group ◽  
High Energy ◽  
Empirical Method ◽  
Side Chain ◽  
Energy Structure ◽  
Direct Dynamics ◽  
Theoretical Mass ◽  
Intermolecular Proton Transfer ◽  
Dynamics Simulations ◽  
Semi Empirical

We present the results of direct dynamics simulations and DFT calculations aimed at elucidating the effect of \textit{O}-sulfonation on the collision induced dissociation for serine. Towards this end, direct dynamics simulations of both serine and sulfoserine were performed at multiple collision energies and theoretical mass spectra obtained. Comparisons to experimental results are favorable for both systems. Peaks related to the sulfo group are identified and the reaction dynamics explored. In particular, three significant peaks (m\z 106, 88, and 81) seen in the theoretical mass spectrum directly related to the sulfo group are analyzed as well as major peaks shared by both systems. Our analysis shows that the m\z 106 peaks result from intramolecular rearrangements, intermolecular proton transfer among complexes composed of initial fragmentation products, and at high energy side-chain fragmentation. The \mz 88 peak was found to contain multiple constitutional isomers, including a previously unconsidered, low energy structure. It was also seen that the RM1 semi empirical method was not able to obtain all of the major peaks seen in experiment for sulfoserine. In contrast, PM6 did obtain all major experimental peaks.

Real-Space X-Ray Pair Distribution Function Analysis and Molecular-Dynamics Modelling of Diflunisal Channel Hydrates

2020 ◽  
Author(s):  
Anuradha Pallipurath ◽  
Francesco Civati ◽  
Jonathan Skelton ◽  
Dean Keeble ◽  
Clare Crowley ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Distribution Function ◽  
Structural Dynamics ◽  
First Principles ◽  
Pair Distribution Function ◽  
Function Analysis ◽  
Real Space ◽  
X Ray ◽  
Dynamics Modelling ◽  
Dynamics Simulations

X-ray pair distribution function analysis is used with first-principles molecular dynamics simulations to study the co-operative H<sub>2</sub>O binding, structural dynamics and host-guest interactions in the channel hydrate of diflunisal.

Schistosomal Sulfotransferase Interaction with Oxamniquine Involves Hybrid Mechanism of Induced-fit and Conformational Selection

2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu
Keyword(s):  
Molecular Dynamics ◽  
Amino Acid ◽  
Binding Energy ◽  
Binding Site ◽  
Conformational Dynamics ◽  
Side Chain ◽  
Amino Acid Residues ◽  
Conformational Selection ◽  
Hybrid Mechanism ◽  
Dynamics Simulations

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.

scholarly journals Electrical Breakdown Mechanism of Transformer Oil with Water Impurity: Molecular Dynamics Simulations and First-Principles Calculations

Crystals ◽  
2021 ◽  
Vol 11 (2) ◽  
pp. 123
Author(s):  
Bin Cao ◽  
Ji-Wei Dong ◽  
Ming-He Chi
Keyword(s):  
Molecular Dynamics ◽  
Thermal Diffusion ◽  
First Principles ◽  
Electrical Breakdown ◽  
Transformer Oil ◽  
Water Molecules ◽  
Conducting Channel ◽  
Breakdown Mechanism ◽  
Water Impurity ◽  
Dynamics Simulations

Water impurity is the essential factor of reducing the insulation performance of transformer oil, which directly determines the operating safety and life of a transformer. Molecular dynamics simulations and first-principles electronic-structure calculations are employed to study the diffusion behavior of water molecules and the electrical breakdown mechanism of transformer oil containing water impurities. The molecular dynamics of an oil-water micro-system model demonstrates that the increase of aging acid concentration will exponentially expedite thermal diffusion of water molecules. Density of states (DOS) for a local region model of transformer oil containing water molecules indicates that water molecules can introduce unoccupied localized electron-states with energy levels close to the conduction band minimum of transformer oil, which makes water molecules capable of capturing electrons and transforming them into water ions during thermal diffusion. Subsequently, under a high electric field, water ions collide and impact on oil molecules to break the molecular chain of transformer oil, engendering carbonized components that introduce a conduction electronic-band in the band-gap of oil molecules as a manifestation of forming a conductive region in transformer oil. The conduction channel composed of carbonized components will be eventually formed, connecting two electrodes, with the carbonized components developing rapidly under the impact of water ions, based on which a large number of electron carriers will be produced similar to “avalanche” discharge, leading to an electrical breakdown of transformer oil insulation. The water impurity in oil, as the key factor for forming the carbonized conducting channel, initiates the electric breakdown process of transformer oil, which is dominated by thermal diffusion of water molecules. The increase of aging acid concentration will significantly promote the thermal diffusion of water impurities and the formation of an initial conducting channel, accounting for the degradation in dielectric strength of insulating oil containing water impurities after long-term operation of the transformer.

scholarly journals Effects of Amino Acid Side-Chain Length and Chemical Structure on Anionic Polyglutamic and Polyaspartic Acid Cellulose-Based Polyelectrolyte Brushes

Polymers ◽  
2021 ◽  
Vol 13 (11) ◽  
pp. 1789
Author(s):  
Dmitry Tolmachev ◽  
George Mamistvalov ◽  
Natalia Lukasheva ◽  
Sergey Larin ◽  
Mikko Karttunen
Keyword(s):  
Glutamic Acid ◽  
Chain Length ◽  
Chemical Structure ◽  
Side Chain ◽  
Side Chains ◽  
Side Chain Length ◽  
Polyelectrolyte Brushes ◽  
Grafting Density ◽  
The Difference ◽  
Cellulose Surface

We used atomistic molecular dynamics (MD) simulations to study polyelectrolyte brushes based on anionic α,L-glutamic acid and α,L-aspartic acid grafted on cellulose in the presence of divalent CaCl2 salt at different concentrations. The motivation is to search for ways to control properties such as sorption capacity and the structural response of the brush to multivalent salts. For this detailed understanding of the role of side-chain length, the chemical structure and their interplay are required. It was found that in the case of glutamic acid oligomers, the longer side chains facilitate attractive interactions with the cellulose surface, which forces the grafted chains to lie down on the surface. The additional methylene group in the side chain enables side-chain rotation, enhancing this effect. On the other hand, the shorter and more restricted side chains of aspartic acid oligomers prevent attractive interactions to a large degree and push the grafted chains away from the surface. The difference in side-chain length also leads to differences in other properties of the brush in divalent salt solutions. At a low grafting density, the longer side chains of glutamic acid allow the adsorbed cations to be spatially distributed inside the brush resulting in a charge inversion. With an increase in grafting density, the difference in the total charge of the aspartic and glutamine brushes disappears, but new structural features appear. The longer sides allow for ion bridging between the grafted chains and the cellulose surface without a significant change in main-chain conformation. This leads to the brush structure being less sensitive to changes in salt concentration.

Sign in / Sign up

Export Citation Format

Share Document