RESUMO
Increased electricity usage over the past several decades has accelerated the need for efficient high-voltage power transmission with reliable insulating materials. Cross-linked polyethylene (XLPE) prepared via dicumyl peroxide (DCP) cross-linking has emerged as the insulator of choice for modern power cables. Although DCP cross-linking generates the desired XLPE product in high yield, other by-products are also produced. One such by-product, acetophenone, is particularly intriguing due to its aromaticity and positive electron affinity. In this work, constrained density functional theory (C-DFT) was utilized to develop an e-ReaxFF force field suitable for describing the acetophenone radical anion. Initial parameters were taken from the 2021 Akbarian e-ReaxFF force field, which was developed to describe XLPE chemistry. Then, C-DFT geometry optimizations were performed wherein an excess electron was constrained to each atom of acetophenone. The resulting C-DFT energy values for the various electronic positions were added to the e-ReaxFF training set. Next, an analogous set of structures was energy-minimized using e-ReaxFF, and equilibrium mixture compositions for the two methods were compared at multiple temperatures. Iterative fitting against C-DFT energy data was performed until satisfactory agreement was achieved. To test force field performance, molecular dynamics simulations were performed in e-ReaxFF and the resulting electronic distributions were qualitatively compared to unconstrained-DFT spin density data. By expanding our e-ReaxFF force field for XLPE, namely, adding the capability to describe acetophenone and its interactions with an excess electron, we move one step closer to a comprehensive molecular understanding of XLPE chemistry in a high-voltage power cable.
RESUMO
Reactive force fields provide an affordable model for simulating chemical reactions at a fraction of the cost of quantum mechanical approaches. However, classically accounting for chemical reactivity often comes at the expense of accuracy and transferability, while computational cost is still large relative to nonreactive force fields. In this Perspective, we summarize recent efforts for improving the performance of reactive force fields in these three areas with a focus on the ReaxFF theoretical model. To improve accuracy, we describe recent reformulations of charge equilibration schemes to overcome unphysical long-range charge transfer, new ReaxFF models that account for explicit electrons, and corrections for energy conservation issues of the ReaxFF model. To enhance transferability we also highlight new advances to include explicit treatment of electrons in the ReaxFF and hybrid nonreactive/reactive simulations that make it possible to model charge transfer, redox chemistry, and large systems such as reverse micelles within the framework of a reactive force field. To address the computational cost, we review recent work in extended Lagrangian schemes and matrix preconditioners for accelerating the charge equilibration method component of ReaxFF and improvements in its software performance in LAMMPS.
RESUMO
Protein-nucleic acid interactions are central to a variety of biological processes, many of which involve large-scale conformational changes that lead to bending of the nucleic acid helix. Here, we focus on the nonsequence-specific protein TRBP, whose double-stranded RNA-binding domains (dsRBDs) interact with the A-form geometry of double-stranded RNA (dsRNA). Crystal structures of dsRBD-dsRNA interactions suggest that the dsRNA helix must bend in such a way that its major groove expands to conform to the dsRBD's binding surface. We show through isothermal titration calorimetry experiments that dsRBD2 of TRBP binds dsRNA with a temperature-independent observed binding affinity (KD â¼500 nM). Furthermore, a near-zero observed heat capacity change (ΔCp = 70 ± 40 cal·mol(-1)·K(-1)) suggests that large-scale conformational changes do not occur upon binding. This result is bolstered by molecular-dynamics simulations in which dsRBD-dsRNA interactions generate only modest bending of the RNA along its helical axis. Overall, these results suggest that this particular dsRBD-dsRNA interaction produces little to no change in the A-form geometry of dsRNA in solution. These results further support our previous hypothesis, based on extensive gel-shift assays, that TRBP preferentially binds to sites of nearly ideal A-form structure while being excluded from sites of local deformation in the RNA helical structure. The implications of this mechanism for efficient micro-RNA processing will be discussed.
