ReaxFF/AMBER-A Framework for Hybrid Reactive/Nonreactive Force Field Molecular Dynamics Simulations

Exploring hydrophobicity or hydrophilicity of borophene surface via reactive molecular dynamics simulation

Borophene, a novel two-dimensional material unveiled in 1998, has garnered significant interest among researchers due to its distinct mechanical and electrical characteristics. Efforts to experimentally synthesize borophene continue to captivate researchers' interest in recent years. Given the current lack of experimental studies on the interaction between water and the borophene surface, molecular dynamics simulation offers a valuable approach for predicting the substance's reactivity with water. Additionally, such simulations can assess the hydrophilicity and hydrophobicity of borophene, providing valuable insights into its properties. In our current research, we utilized reactive molecular dynamics simulation to investigate the wetting behavior of borophene. Our findings reveal that the borophene surface exhibits hydrophobic characteristics, demonstrating anisotropic wettability. Specifically, the water contact angle was calculated to be 149.11° along the zigzag direction and 148.4° along the armchair direction. The contour map of the interaction energy between a water molecule and the borophene surface revealed a notable energy barrier in the zigzag direction. This barrier contributes to the asymmetric spreading of the water droplet on the surface. Density profiles and radial pair distribution function (RDF) diagrams of the water droplet on the borophene surface further corroborated the hydrophobic nature of borophene by indicating a significant distance between the water droplet and the surface. Moreover, analysis of the number of hydrogen bonds demonstrated that borophene efficiently utilizes nearly all its capacity to form hydrogen bonds. Additionally, we compared the wettability of borophene with that of other two-dimensional materials, such as various graphene allotropes and phosphorene, which have been subjects of recent investigation.

ReaxFF molecular dynamics of graphene oxide/NaCl aqueous solution interfaces

In this work, the interaction of NaCl aqueous solution with graphene (G), graphene oxide (GO), and graphite oxide (GTO) is studied using the ReaxFF module of Amsterdam Modeling Suite (AMS) software. We consider four models using the NaCl aqueous solution, containing a graphene sheet (G), a single sheet of GO with epoxide and hydroxyl groups on its surface, 4 layers of GO to model GTO, and a bulk NaCl solution as a reference. The structural and dynamical properties of G, GO, and GTO were quantified by analyzing the functional groups, radial distribution functions, density profiles and diffusivities of water and ions. Due to the reactive force field, the systems underwent spontaneous modification of surface functional groups during the first 750 ps after which the structure stabilizes (the energy stabilizes in less than 400 ps). Pristine graphene in contact with the NaCl solution formed hydroxyl groups on the edges, i.e., converted to partially reduced graphene oxide. The epoxy groups (Oe) on the initial GO were rather unstable, leading to a reduction of their number, however, there was an increase in the number of hydroxyl groups (Oh), mainly at the edges. The interactions of NaCl with the carbon-based sheets are rather weak, including GO and GTO which are decorated with numerous functional groups. Diffusion coefficients of water agree with the available data, but discrepancies in Na+ and Cl- diffusivity compared to other references underscore the need for further development in the dynamic parameters of the reactive force field used. In essence, our research provides specific data previously unreported, laying a foundation for advancing water desalination system design. The study's novelty lies in its realistic approach to graphene/graphene oxide modification, comprehensive characterization, and the application of the reactive force field to explore the graphene oxide-NaCl aqueous interface, contributing to the development of a practical membrane system for water desalination.

Application of molecular dynamics simulation in self-assembled cancer nanomedicine

Self-assembled nanomedicine holds great potential in cancer theragnostic. The structures and dynamics of nanomedicine can be affected by a variety of non-covalent interactions, so it is essential to ensure the self-assembly process at atomic level. Molecular dynamics (MD) simulation is a key technology to link microcosm and macroscale. Along with the rapid development of computational power and simulation methods, scientists could simulate the specific process of intermolecular interactions. Thus, some experimental observations could be explained at microscopic level and the nanomedicine synthesis process would have traces to follow. This review not only outlines the concept, basic principle, and the parameter setting of MD simulation, but also highlights the recent progress in MD simulation for self-assembled cancer nanomedicine. In addition, the physicochemical parameters of self-assembly structure and interaction between various assembled molecules under MD simulation are also discussed. Therefore, this review will help advanced and novice researchers to quickly zoom in on fundamental information and gather some thought-provoking ideas to advance this subfield of self-assembled cancer nanomedicine.

Molecular Insight into the Pozzolanic Reaction of Metakaolin and Calcium Hydroxide

The reaction mechanism of the pozzolanic reaction of metakaolin (MK) from the atomic point of view has not yet been explored. To explain the process and mechanism of the pozzolanic reaction from the atomic point of view, molecular insight into the pozzolanic reaction of MK and calcium hydroxide (CH) was analyzed through the reaction molecular dynamics (MD) simulation. The results show that the pozzolanic reaction of MK and CH can be essentially regarded as the CH decomposition and penetration into MK. Also, the structure evolution after the pozzolanic reaction shows that the water molecules cannot penetrate the MK structure till the participation of Ca²⁺ and OH^- ions of CH. The Ca²⁺ and OH^- ions have strong interaction with MK and drill into the MK structure, followed by the destruction of a part of the MK structure and water penetration. The final structure of CH removed by MK can be regarded as the precursor of the CASH gel structure.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Effect of the R126C mutation on the structure and function of the glucose transporter GLUT1: A molecular dynamics simulation study

Glucose transporter 1 (GLUT1) is responsible for basal glucose uptake and is expressed in most tissues under normal conditions. GLUT1 mutations can cause early-onset absence epilepsy and myoclonus dystonia syndrome (MDS), with MDS potentially lethal. In this study, the effect of the R126C mutation, which is associated with MDS, on structural stability and substrate transport of GLUT1 was investigated. Various bioinformatics tools were used to predict the stability of GLUT1, revealing that the R126C mutation reduces the structural stability of GLUT1. Molecular dynamics (MD) simulations were used to further characterize the effect of the R126C mutation on GLUT1 structural stability. Based on the MD simulations, specific conformational changes and dominant motions of the GLUT1 mutant were characterized by Principal component analysis (PCA). The mutation disrupts hydrogen bonds between substrate-binding residues and glucose, thus likely reducing substrate affinity. The R126C mutation reduces the conformational stability of the protein, and fewer intramolecular hydrogen bonds were present in the mutated GLUT1 when compared with that of wild-type GLUT1. The mutation increased the free energy of glucose transport through GLUT1 significantly, especially at the mutation site, indicating that passage of glucose through the channel is hindered, and this mutant may even release cytoplasmic glucose. This study provides a detailed atomic-level explanation for the reduced structural stability and substrate transport capacity of a GLUT1 mutant. The results aid our understanding of the structure of GLUT1 and provide a framework for developing drugs to treat GLUT1-related diseases, such as MDS.

A Polarizable Cationic Dummy Metal Ion Model

A novel locally polarizable multisite model based on the original cation dummy atom (CDA) model is described for molecular dynamics simulations of ions in condensed phases. Polarization effects are introduced by the electronegativity equalization model (EEM) method where charges on the metal ion and its dummy atoms can fluctuate to respond to the environment. This model includes explicit polarization and ion-induced interactions and can be coupled with nonpolarizable or polarizable water models, making it more transferable to simpler force fields. This approach allows us to enhance the original fixed charge CDA model where the charge distribution cannot adapt to the local solvent structure. To illustrate the new CDApol model, we examined properties of the Zn2+, Al3+, and Zr4+ ions in aqueous solution. The polarizable model and Lennard-Jones parameters were refined for octahedrally coordinated Zn2+, Al3+, and Zr4+ CDAs to reproduce thermodynamic and geometrical properties. Using this locally polarizable model, we were able to obtain the experimental hydration free energy, ion-oxygen distance, and coordination number coupled with the standard 12-6 Lennard-Jones model. This model can be used in myriad additional applications where local polarization and charge transfer is important.

Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors

Bromodomain and extra-terminal domain (BET) subfamily is the most studied subfamily of bromodomain-containing proteins (BCPs) family which can modulate acetylation signal transduction and produce diverse physiological functions. Thus, the BET family can be treated as an alternative strategy for targeting androgen-receptor (AR)-driven cancers. In order to explore the effect of inhibitors binding to BRD4 (the most studied member of BET family), four 150 ns molecular dynamic simulations were performed (free BRD4, Cpd4-BRD4, Cpd9-BRD4 and Cpd19-BRD4). Docking studies showed that Cpd9 and Cpd19 were located at the active pocket, as well as Cpd4. Molecular dynamics (MD) simulations indicated that only Cpd19 binding to BRD4 can induce residue Trp81-Ala89 partly become α-helix during MD simulations. MM-GBSA calculations suggested that Cpd19 had the best binding effect with BRD4 followed by Cpd4 and Cpd9. Computational alanine scanning results indicated that mutations in Phe83 made the greatest effects in Cpd9-BRD4 and Cpd19-BRD4 complexes, showing that Phe83 may play crucial roles in Cpd9 and Cpd19 binding to BRD4. Our results can provide some useful clues for further BCPs family search.

Recent Advances for Improving the Accuracy, Transferability, and Efficiency of Reactive Force Fields

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.