Molecular Dynamics Model Research Articles

Mechanical Properties of Small Quasi-Square Graphene Nanoflakes

The rise of straintronics—the possibility of fine-tuning the electronic properties of nanosystems by applying strain to them—has enhanced the interest in characterizing the mechanical properties of these systems when they are subjected to tensile (or compressive), shear and torsion strains. Four parameters are customarily used to describe the mechanical behavior of a macroscopic solid within the elastic regime: Young’s and shear moduli, the torsion constant and Poisson’s ratio. There are some relations among these quantities valid for elastic continuous isotropic systems that are being used for 2D nanocrystals without taking into account the non-continuous anisotropic nature of these systems. We present in this work computational results on the mechanical properties of six small quasi-square (aspect ratio between 0.9 and 1.1) graphene nanocrystals using the PM7 semiempirical method. We use the results obtained to test the validity of two relations derived for macroscopic homogeneous isotropic systems and sometimes applied to 2D systems. We show they are not suitable for these nanostructures and pinpoint the origin of some discrepancies in the elastic properties and effective thicknesses reported in the literature. In an attempt to recover one of these formulas, we introduce an effective torsional thickness for graphene analogous to the effective bending thickness found in the literature. Our results could be useful for fitting interatomic potentials in molecular mechanics or molecular dynamics models for finite carbon nanostructures, especially near their edges and for twisted systems.

Network pharmacology and experimental insights into STAT3 inhibition by novel isoxazole derivatives of piperic acid in triple negative breast cancer

STAT3 is a crucial member within a family of seven essential transcription factors. Elevated STAT3 levels have been identified in various cancer types, notably in breast cancer (BC). Consequently, inhibiting STAT3 is recognized as a promising and effective strategy for therapeutic intervention against breast cancer. We herein synthesize a library of isoxazole (PAIs) from piperic acid [2E, 4E)-5-(2H-1,3-Benzodioxol-5-yl) penta-2,4-dienoic acid] on treatment with propargyl bromide followed by oxime under prescribed reaction conditions. Piperic acid was obtained by hydrolysis of piperine extracted from Piper nigrum. First, we checked the binding potential of isoxazole derivatives with breast cancer target proteins by network pharmacology, molecular docking, molecular dynamic (MD) simulation and cytotoxicity analysis as potential anti-breast cancer (BC) agents. The multi-source databases were used to identify possible targets for isoxazole derivatives. A network of protein-protein interactions (PPIs) was generated by obtaining 877 target genes that overlapped gene symbols associated with isoxazole derivatives and BC. Molecular docking and MD modelling demonstrated a strong affinity between isoxazole derivatives and essential target genes. Further, the cell viability studies of isoxazole derivatives on the human breast carcinoma cell lines showed toxicity in all breast cancer cell lines. In summary, our study indicated that the isoxazole derivative showed the significant anticancer activity. The results highlight the prospective utility of isoxazole derivatives as new drug candidates for anticancer chemotherapy, suggesting route for the continued exploration and development of drugs suitable for clinical applications.

Investigation of Alkane-Assisted dodecylamine on flotation desilication of magnesite Tailings: Molecular dynamics simulations and model optimization

Investigating the reverse flotation desilication process of magnesite tailings (MT) holds substantial merit in the context of magnesium metal recovery from solid waste. However, scant attention has been devoted to exploring the influence of the carbon chain length of diverse alkanes on the reverse flocfloatation desilication of MT. This research systematically delved into the selective adsorption behavior exhibited by dodecylamine (DDA) and its compound collectors on the surfaces of quartz and magnesite. As the carbon chain length of alkanes increases, the efficacy of reverse flotation desilication of MT diminishes. Under optimized conditions, the concentration of MgO in the concentrate registers a 0.93% increment, accompanied by a 2.05% reduction in SiO2 content. The optimal MgO grade attains 42.68%. The augmenting carbon chain length of alkanes leads to a decline in the dispersibility of the compound collector and a reduction in the diffusion coefficient of the monolayer film at the alkane/water interface. The compound collector, comprising n-heptane and DDA, manifests a noteworthy synergistic promotion effect on the flotation desilication of MT. This effect reduces the multilayer adsorption of DDA, promotes more hydrogen bonding formation, and favors DDA adsorption on the quartz surface over magnesite. Consequently, issues related to froth entrapment and poor selectivity are alleviated. Longer alkane chains in the compound collectors result in a stronger van der Waals force effect between the carbon chains. This effect causes DDDA to slide more easily on the quartz surface, which is not conducive to enhancing the hydrophobicity of the surface. Furthermore, alkane-induced hydrophobic flocculation increases the floccules size of quartz, enhancing the likelihood of bubble-particle adhesion. This investigation elucidates the molecular-level mechanism underlying alkane hydrophobic flocculation with varying chain lengths, contributing a theoretical foundation and technical support for the large-scale MT management and the sustainable production of magnesium metal.

Micromechanical Dilution of PLA/PETG-Glass/Iron Nanocomposites: A More Efficient Molecular Dynamics Approach.

Polylactic acid (PLA) and poly(ethylene terephthalate glycol) (PETG) are popular thermoplastics used in additive manufacturing applications. The mechanical properties of PLA and PETG can be significantly improved by introducing fillers, such as glass and iron nanoparticles (NPs), into the polymer matrix. Molecular dynamics (MD) simulations with the reactive INTERFACE force field were used to predict the mechanical responses of neat PLA/PETG and PLA-glass/iron and PETG-glass/iron nanocomposites with relatively high loadings of glass/iron NPs. We found that the iron and glass NPs significantly increased the elastic moduli of the PLA matrix, while the PETG matrix exhibited modest increases in elastic moduli. This difference in reinforcement ability may be due to the slightly greater attraction between the glass/iron NP and PLA matrix. The NASA Multiscale Analysis Tool was used to predict the mechanical response across a range of volume percent glass/iron filler by using only the neat and highly loaded MD predictions as input. This provides a faster and more efficient approach than creating multiple MD models per volume percent per polymer/filler combination. To validate the micromechanics predictions, experimental samples incorporating hollow glass microspheres (MS) and carbonyl iron particles (CIP) into PLA/PETG were developed and tested for elastic modulus. The CIP produced a larger reinforcement in elastic modulus than the MS, with similar increases in elastic modulus between PLA/CIP and PETG/CIP at 7.77 vol % CIP. The micromechanics-based mechanical predictions compare excellently with the experimental values, validating the integrated micromechanical/MD simulation-based approach.

A Molecular Dynamics Analysis of the Thickness and Adhesion Characteristics of the Quasi-Liquid Layer at the Asphalt-Ice Interface.

The quasi-liquid layer (QLL), a microstructure located between ice and an adhering substrate, is critical in generating capillary pressure, which in turn influences ice adhesion behavior. This study employed molecular dynamics (MD) methods to obtain QLL thickness and utilized these measurements to estimate the adhesive strength between ice and asphalt. The research involved constructing an ice-QLL-asphalt MD model, encompassing four asphalt types and five temperature ranges from 250 K to 270 K. The QLL thickness was determined for various asphalts and temperatures using the tetrahedral order parameter gradient. Additionally, capillary pressure was calculated based on the QLL thickness and other geometric parameters obtained from the MD analysis. These findings were then compared with ice adhesion strength data acquired from pull-off tests. The results indicate that QLL thickness varies with different asphalt types and increases with temperature. At a constant temperature, the QLL thickness decreases in the order of the basal plane, primary prism plane, and secondary prism plane. Furthermore, the adhesion strength of the QLL diminishes as the temperature rises, attributed to the disruption of hydrogen bonds at lower temperatures. The greater the polarity of the asphalt's interface molecules, the stronger the adhesion strength and binding free energy. The MD simulations of the asphalt-ice interface offer insights into the atomic-scale adhesive properties of this interface, contributing to the enhancement in QLL property prediction and calibration at larger scales.

C9orf72 polyPR directly binds to various nuclear transport components

The disruption of nucleocytoplasmic transport (NCT) is an important mechanism in neurodegenerative diseases. In the case of C9orf72-ALS, trafficking of macromolecules through the nuclear pore complex (NPC) might get frustrated by the binding of C9orf72-translated arginine-containing dipeptide repeat proteins (R-DPRs) to the Kapβ family of nuclear transport receptors. Besides Kapβs, several other types of transport components have been linked to NCT impairments in R-DPR-expressed cells, but the molecular origin of these observations has not been clarified. Here, we adopt a coarse-grained molecular dynamics model at amino acid resolution to study the direct interaction between polyPR, the most toxic DPR, and various nuclear transport components to elucidate the binding mechanisms and provide a complete picture of potential polyPR-mediated NCT defects. We found polyPR to directly bind to several isoforms of the Impα family, CAS (the specific exporter of Impα) and RanGAP. We observe no binding between polyPR and Ran. Longer polyPRs at lower salt concentrations also make contact with RanGEF and NTF2. Analyzing the polyPR contact sites on the transport components reveals that polyPR potentially interferes with RanGTP/RanGDP binding, with nuclear localization signal (NLS)-containing cargoes (cargo-NLS) binding to Impα, with cargo-NLS release from Impα, and with Impα export from the nucleus. The abundance of polyPR-binding sites on multiple transport components combined with the inherent polyPR length dependence makes direct polyPR interference of NCT a potential mechanistic pathway of C9orf72 toxicity.

C9orf72 polyPR directly binds to various nuclear transport components.

The disruption of nucleocytoplasmic transport (NCT) is an important mechanism in neurodegenerative diseases. In the case of C9orf72-ALS, trafficking of macromolecules through the nuclear pore complex (NPC) might get frustrated by the binding of C9orf72-translated arginine-containing dipeptide repeat proteins (R-DPRs) to the Kapβ family of nuclear transport receptors. Besides Kapβs, several other types of transport components have been linked to NCT impairments in R-DPR-expressed cells, but the molecular origin of these observations has not been clarified. Here, we adopt a coarse-grained molecular dynamics model at amino acid resolution to study the direct interaction between polyPR, the most toxic DPR, and various nuclear transport components to elucidate the binding mechanisms and provide a complete picture of potential polyPR-mediated NCT defects. We found polyPR to directly bind to several isoforms of the Impα family, CAS (the specific exporter of Impα) and RanGAP. We observe no binding between polyPR and Ran. Longer polyPRs at lower salt concentrations also make contact with RanGEF and NTF2. Analyzing the polyPR contact sites on the transport components reveals that polyPR potentially interferes with RanGTP/RanGDP binding, with nuclear localization signal (NLS)-containing cargoes (cargo-NLS) binding to Impα, with cargo-NLS release from Impα, and with Impα export from the nucleus. The abundance of polyPR-binding sites on multiple transport components combined with the inherent polyPR length dependence makes direct polyPR interference of NCT a potential mechanistic pathway of C9orf72 toxicity.

Coupling strengthening mechanism of precipitate and solute atoms in edge dislocation motion of Al-Mg-Si alloy

To investigate the coupling strengthening mechanism of precipitate and Cu solute atoms in Al-Mg-Si alloy, the molecular dynamics model of edge dislocation motion interacted with precipitate and Cu solute atoms was established. The dislocation length is increased in the dislocation-precipitate interaction due to the bending and bowing of the edge dislocation around the precipitate. No coupling between precipitate and solute can be found in material strengthening in cases with lower solute concentrations. The exponential effect of the coupling effects between precipitate and solute on material strengthening is found when the solute concentration is increased above 0.3.at% at room temperature. The change of the edge dislocation caused by solute atoms around the precipitate is the key reason for the strengthening coupling effect. Furthermore, the effects of temperature variations on the coupling effect are more pronounced than the effects of shear strain rate variations. When the temperature is increased from 200K to 500K, the critical ratio of solute and precipitate strengthening is decreased from 1.5 to 0.4.

A reactive molecular dynamics model for uranium/hydrogen containing systems.

Uranium-based materials are valuable assets in the energy, medical, and military industries. However, understanding their sensitivity to hydrogen embrittlement is particularly challenging due to the toxicity of uranium and the computationally expensive nature of quantum-based methods generally required to study such processes. In this regard, we have developed a Chebyshev Interaction Model for Efficient Simulation (ChIMES) that can be employed to compute energies and forces of U and UH3 bulk structures with vacancies and hydrogen interstitials with accuracy similar to that of Density Functional Theory (DFT) while yielding linear scaling and orders of magnitude improvement in computational efficiency. We show that the bulk structural parameters, uranium and hydrogen vacancy formation energies, and diffusion barriers predicted by the ChIMES potential are in strong agreement with the reference DFT data. We then use ChIMES to conduct molecular dynamics simulations of the temperature-dependent diffusion of a hydrogen interstitial and determine the corresponding diffusion activation energy. Our model has particular significance in studies of actinides and other high-Z materials, where there is a strong need for computationally efficient methods to bridge length and time scales between experiments and quantum theory.

Improvement of the Interface Contact between Solid Polymer Electrolytes and Electrodes via Ultrasonic-Press Fusion

All solid-state lithium batteries have high energy density, good safety, and broad prospects. However, the poor compatibility of the electrolyte/electrode interface limits its development. To improve the interface contact performance of batteries, this study proposes an ultrasonic-press fusion method. The electrolyte/electrode structure was pressed by an ultrasonic sonotrode, and fused tightly via thermal and mechanical effects produced by ultrasonic-press fusion. The results showed that the charge transfer impedance at the battery interface was reduced by 80.3% in AC impedance testing, and the battery interface impedance performance was significantly improved. By observing the surface and interface morphology of the electrolyte, it was found that ultrasonic-press fusion improved the electrolyte interface contact. The electrolyte melted and produced microscale rheological behavior under the action of ultrasonic-press fusion, ultimately filling the capillary structure of anode in the battery, improving contact performance of the interface. Furthermore, a molecular dynamics model of the interface was constructed to calculate the contact performance. It was found that ultrasonic-press fusion generates thermal and mechanical effects. The thermal effect facilitates the electrolyte melts, and the mechanical effect causes the melt to vibrate, which coordinate and act on the electrolyte and cause the electrolyte to solder to the anode.

Mechanisms of Resistance to Tyrosine Kinase Inhibitors in ROS1 Fusion-Positive Nonsmall Cell Lung Cancer.

ROS1 fusion-positive (ROS1+) nonsmall cell lung cancer (NSCLC) patients are highly sensitive to tyrosine kinase inhibitor (TKI) treatments. However, acquired TKI resistance remains the major hurdle preventing patients from experiencing prolonged benefits. 107 advanced or metastatic ROS1+ NSCLC patients who progressed on crizotinib and lorlatinib were recruited. Tissue and plasma samples were collected at baseline (N = 50), postcrizotinib (N = 91), and postlorlatinib (N = 21), which were all subject to the 139-gene targeted next-generation DNA sequencing. Molecular dynamics modeling was performed to investigate the effects of ROS1 mutations on binding to different TKIs. In patients with postcrizotinib and postlorlatinib samples, an accumulation of on- and off-target resistance alterations after multiple TKI treatments was observed. ROS1 G2032R and MET amplification were the most common on-target and off-target alterations, respectively. Patients with CD74-ROS1 and SLC34A2-ROS1 had longer progression-free survival (PFS) (P < 0.001) and higher rates of resistance mutations (on-target, P = 0.001; off-target, P = 0.077) than other ROS1 fusion variants following crizotinib treatment. Ten distinct on-target resistance mutations were detected after TKI therapies, of which 4 were previously unreported (ROS1 L2010M, G1957A, D1988N, L1982V). Molecular dynamics simulations showed that all 4 mutations were refractory to crizotinib, while G1957A, D1988N, and L1982V were potentially sensitive to lorlatinib and entrectinib. This study provided a comprehensive portrait of TKI-resistance mechanisms in ROS1+ NSCLC patients. Using in silico simulations of TKI activity, novel secondary mutations that may confer TKI resistance were identified and may support clinical therapeutic decision-making.

Unraveling the Mechanical Property Decrease of Electrospun Spider Silk: A Molecular Dynamics Simulation Study.

This study investigated the impact of electric fields on Nephila clavipes spider silk using molecular dynamics modeling. Electric fields with varying amplitudes and directions were observed to disrupt the β sheet structure of spider silk and reduce its mechanical properties. However, a notable exception was observed when a 0.1 V/nm electric field was applied in the antiparallel direction, resulting in improvements in Young's modulus and ultimate tensile strength. The antiparallel direction was observed to be particularly sensitive to electric fields, causing disruptions in beta sheets and hydrogen bonds, which significantly influence the mechanical properties. This study demonstrates that spider silk maintains its structural integrity at 0.1 V/nm. Possibly, lowering the power levels of typical electrospinning machines can prevent secondary structural disruption. These findings provide valuable insights for enhancing silk fiber production and applications using natural silk proteins while shedding light on the impact of electric fields on other silk proteins. Finally, this study opens up possibilities for optimizing electrospinning processes to enhance performance in various silk electrospinning applications.

Dissipation and recovery in collagen fibrils under cyclic loading: A molecular dynamics study.

The hysteretic behavior exhibited by collagen fibrils, when subjected to cyclic loading, is known to result in both dissipation as well as accumulation of residual strain. On subsequent relaxation, partial recovery has also been reported. Cross-links have been considered to play a key role in overall mechanical properties. Here, we modify an existing coarse-grained molecular dynamics model for collagen fibril with initially cross-linked collagen molecules, which is known to reproduce the response to uniaxial strain, by incorporating reformation of cross-links to allow for possible recovery of the fibril. Using molecular dynamics simulations, we show that our model successfully replicates the key features observed in experimental data, including the movement of hysteresis loops, the time evolution of residual strains and energy dissipation, as well as the recovery observed during relaxation. We also show that the characteristic cycle number, describing the approach toward steady state, has a value similar to that in experiments. We also emphasize the vital role of the degree of cross-linking on the key features of the macroscopic response to cyclic loading.

Ion Conductivity in Salt-Doped Polymers: Combined Effects of Temperature and Salt Concentration.

We construct a coarse-grained molecular dynamics model based on poly(ethylene oxide) and lithium bis(trifluoromethane)sulfonimide salt to examine the combined effects of temperature and salt concentration on the transport properties. Salt doping notably slows the dynamics of polymer chains and reduces ion diffusivity, resulting in a glass transition temperature increase proportional to the salt concentration. The polymer diffusion is shown to be well represented by a modified Vogel-Fulcher-Tamman (M-VFT) equation that accounts for both the temperature and salt concentration dependence. Furthermore, we find that, at any temperature, the concentration dependence of the conductivity is well described by the product of its infinite dilution value and a correction factor accounting for the reduced segmental mobility with increasing salt concentration. These results highlight the important role of polymer segmental mobility in the salt concentration dependence of ion conductivity for temperatures near and above the glass transition.

Momentum correlation of light nuclei in Au + Au collisions at [formula omitted] GeV

Within the Ultra-relativistic Quantum Molecular Dynamics model (UrQMD) coupled with nucleon coalescence model and Mini-Spanning-Tree model, the yields of light nuclei have been stimulated in Au + Au collisions over an energy range of sNN=2.0∼7.7GeV and the momentum correlation functions of light nuclei pairs have been calculated by both the Lednický-Lyuboshitz and the Correlation After Burner methods. We compare the yields of light nuclei and their momentum correlation functions at midrapidity in this energy region with experimental data. It is found that there are differences between the results of the two models, and the coalescence method seems less valid at low collision energy. Furthermore, both the peak values of proton-proton correlation functions and the transition point of elliptic flows from out-of-plane to in-plane emission show a turning point around 3-4 GeV, which suggests that there is a relation between momentum correlation function and collective flow of particles.

In situ catalytic polymerization of LiNO3-containing PDOL electrolytes for high-energy quasi-solid-state lithium metal batteries

1,3-Dioxolane (DOL) can be induced to form polymer electrolytes by in situ ring opening polymerization at Lewis acid salts. However, the presence of ion–dipole interactions between NO3− and DOL suppresses this polymerization behavior. Herein, we report that the use of Sc(OTf)3 as an initiator can well disrupt this ion–dipole interaction. The in situ ring-opening polymerization of DOL is achieved to form a LiNO3-modified polyDOL-based electrolyte. New mechanism of Li migration in polymer electrolytes have been revealed by molecular dynamics modeling and constructed molecular interface models. It is found that the electrolyte interior is affected by the ion–dipole interactions of Li+ by NO3−, which made it easier for Li+ to be released from the polar ether-oxygen ligands on the polymer chain for rapid migration, thus exhibiting an ionic conductivity of 1.8 mS cm−1 and tLi+ of 0.78. In addition, NO3− preempts the binding sites around Li+, improves the coordination environment, and prioritizes the formation of a kinetically stable anion-derivative interface, which effectively mitigates the interfacial side reactions between the electrodes and the electrolytes. As a result, the assembled solid-state Li||LiFePO4 cell exhibits an impressive 152.3 mAh/g discharge capacity at 0.5C and maintains 80.3 % of the capacity after 450 cycles at 50 °C. This work not only opens a new avenue for designing high-performance gel polymer electrolytes for more metal-based batteries, but also provides valuable insights into understanding the ion migration mechanism.

Mechanistic QSAR modeling derived virtual screening, drug repurposing, ADMET and in-vitro evaluation to identify anticancer lead as lysine-specific demethylase 5a inhibitor

A lysine-specific demethylase is an enzyme that selectively eliminates methyl groups from lysine residues. KDM5A, also known as JARID1A or RBP2, belongs to the KDM5 Jumonji histone demethylase subfamily. To identify novel molecules that interact with the LSD5A receptor, we created a quantitative structure-activity relationship (QSAR) model. A group of 435 compounds was used in a study of the quantitative relationship between structure and activity to guess the IC50 values for blocking LASD5A. We used a genetic algorithm-multilinear regression-based quantitative structure-activity connection model to forecast the bioactivity (PIC50) of 1615 food and drug administration pharmaceuticals from the zinc database with the goal of repurposing clinically used medications. We used molecular docking, molecular dynamic simulation modelling, and molecular mechanics generalised surface area analysis to investigate the molecule’s binding mechanism. A genetic algorithm and multi-linear regression method were used to make six variable-based quantitative structure-activity relationship models that worked well (R 2 = 0.8521, Q 2 LOO = 0.8438, and Q 2 LMO = 0.8414). ZINC000000538621 was found to be a new hit against LSD5A after a quantitative structure-activity relationship-based virtual screening of 1615 zinc food and drug administration compounds. The docking analysis revealed that the hit molecule 11 in the KDM5A binding pocket adopted a conformation similar to the pdb-6bh1 ligand (docking score: −8.61 kcal/mol). The results from molecular docking and the quantitative structure-activity relationship were complementary and consistent. The most active lead molecule 11, which has shown encouraging results, has good absorption, distribution, metabolism, and excretion (ADME) properties, and its toxicity has been shown to be minimal. In addition, the MTT assay of ZINC000000538621 with MCF-7 cell lines backs up the in silico studies. We used molecular mechanics generalise borne surface area analysis and a 200-ns molecular dynamics simulation to find structural motifs for KDM5A enzyme interactions. Thus, our strategy will likely expand food and drug administration molecule repurposing research to find better anticancer drugs and therapies. Communicated by Ramaswamy H. Sarma

Unraveling Molecular Recognition of Glycan Ligands by Siglec-9 via NMR Spectroscopy and Molecular Dynamics Modeling.

Human sialic-acid-binding immunoglobulin-like lectin-9 (Siglec-9) is a glycoimmune checkpoint receptor expressed on several immune cells. Binding of Siglec-9 to sialic acid containing glycans (sialoglycans) is well documented to modulate its functions as an inhibitory receptor. Here, we first assigned the amino acid backbone of the Siglec-9 V-set domain (Siglec-9d1), using well-established triple resonance three-dimensional nuclear magnetic resonance (NMR) methods. Then, we combined solution NMR and molecular dynamic simulation methods to decipher the molecular details of the interaction of Siglec-9 with the natural ligands α2,3 and α2,6 sialyl lactosamines (SLN), sialyl Lewis X (sLeX), and 6-O sulfated sLeX and with two synthetically modified sialoglycans that bind with high affinity. As expected, Neu5Ac is accommodated between the F and G β-strands at the canonical sialic acid binding site. Addition of a heteroaromatic scaffold 9N-5-(2-methylthiazol-4-yl)thiophene sulfonamide (MTTS) at the C9 position of Neu5Ac generates new interactions with the hydrophobic residues located at the G-G' loop and the N-terminal region of Siglec-9. Similarly, the addition of the aromatic substituent (5-N-(1-benzhydryl-1H-1,2,3-triazol-4-yl)methyl (BTC)) at the C5 position of Neu5Ac stabilizes the conformation of the long and flexible B'-C loop present in Siglec-9. These results expose the underlying mechanism responsible for the enhanced affinity and specificity for Siglec-9 for these two modified sialoglycans and sheds light on the rational design of the next generation of modified sialoglycans targeting Siglec-9.

Influence of oxygen functional groups on the crack growth behavior of graphene with edge and middle crack: A molecular dynamics study

In the last decade, graphene has found a wide role and application in nanocomposites and nanodevices for its extraordinary properties. However, due to its low fracture toughness and brittleness, graphene is brittle, and crack growth occurs rapidly. As a result, crack growth behavior and control of its progression are important. In this study, molecular dynamics modeling was used to investigate the crack growth behavior and mechanical properties in functionalized graphene (with middle and edge crack) with different oxide groups (hydroxyl, epoxide, and carboxyl) at different percentages of graphene oxide (0–3 wt. %). The results show that functionalized graphene with oxide groups changes the properties of graphene, effectively reducing the crack growth in crack length (4.9 Å), crack opening displacement (1.4 Å), tensile strength for edge crack (74 to 54 Gpa), and middle crack (128 to 63 Gpa). Also, carboxyl and epoxide groups have a greater reduction in crack growth rate, opening displacement growth rate, and stress concentration. In addition, increase in the number of atoms in oxide groups around the initial crack (with radius 13 Å) reduces crack growth rate, opening displacement growth rate, and stress concentration. This study provides a broader insight into crack growth behavior in functionalized graphene with oxide groups.

Environmentally realistic dose of tire-derived metabolite 6PPD-Q exposure causes intestinal jejunum and ileum damage in mice via cannabinoid receptor-activated inflammation

N-(1,3-dimethylbutyl)-N′-phenyl-p-phenylenediamine quinone (6PPD-Q) is a quinone derivative of a common tire additive 6PPD, whose occurrence has been widely reported both in the environment and human bodies including in adults, pregnant women and children. Yet, knowledge on the potential intestinal toxicity of 6PPD-Q in mammals at environmentally relevant dose remain unknown. In this study, the effects of 6PPD-Q on the intestines of adult ICR mice were evaluated by orally administering environmentally relevant dose or lower levels of 6PPD-Q (0.1, 1, 10, and 100 μg/kg) for 21 days. We found that 6PPD-Q disrupted the integrity of the intestinal barrier, mostly in the jejunum and ileum, but not in the duodenum or colon, in a dose-dependent manner. Moreover, intestinal inflammation manifested with elevated levels of TNF-α, IL-1, and IL-6 mostly observed in doses at 10 and 100 μg/kg. Using reverse target screening technology combining molecular dynamic simulation modeling we identified key cannabinoid receptors including CNR2 activation to be potentially mediating the intestinal inflammation induced by 6PPD-Q. In summary, this study provides novel insights into the toxic effects of emerging contaminant 6PPD-Q on mammalian intestines and that the chemical may be a cannabinoid receptor agonist to modulate inflammation.