Aggregated Single-Walled Carbon Nanotubes Absorb and Deform Dopamine-Related Proteins Based on Molecular Dynamics Simulations

The interface hydrophilic-hydrophobic integration of fluorinated defective graphene towards biomedical applications

In biomedical fields, rational design of novel two-dimensional (2D) biomedical nanomaterials aims to precisely manipulate biomolecules, including efficient capture, structural-functional transformation, directional movement, and self-assembly. In this work, we innovatively proposed new graphene nanosheets and selected two representative proteins to explore their binding mechanisms, structural-functional transformation of proteins, and biological effects of the materials. Fluorinated defective graphene (FDG) exhibited highly efficient capture and structural-functional transformation for the receptor binding domain (RBD), and we observed its collapse phenomenon in 2D materials for the first time. For the main protease (Mpro), FDG achieved an optimal balance between efficient capture, immobilization, and structural disruption. Further studies showed that fluorination on oxygen-containing defect graphene significantly enhanced variances in water distribution, surface properties, and hydrogen bond networks on the material surface. This allowed amino acids to be confined to specific areas, achieving efficient capture and directional movement. Additionally, the adsorption behavior and interaction strength of peptides and deoxynucleotides on FDG further validated the possibility of self-assembly. In summary, we highlight FDG as an excellent biomedical material with hydrophilic-hydrophobic integration.

Unwinding DNA strands by single-walled carbon nanotubes: Molecular docking and MD simulation approach

Despite the growing research into the use of carbon nano-tubes (CNTs) in science and medicine, concerns about their potential toxicity remain insufficiently studied. This study utilizes molecular docking calculations combined by molecular dynamics simulations to investigate the dynamic intricacies of the interaction between single-walled carbon nanotubes (swCNTs) and double-stranded DNA (dsDNA). By examining the influence of swCNT characteristics such as length, radius, and chirality, our findings shed light on the complex interplay that shapes the binding affinity and stability of the dsDNA-swCNT complex. Molecular docking results identify a zigzag swCNT, with a radius of 0.16 Å and a length of 38 Å, as exhibiting the highest binding affinity with dsDNA (-23.9 kcal/mol). Comprehensive analyses, spanning docking results, binding energies, RMSD, radius of gyration, and potential of mean force (PMF) profiles, provide a detailed understanding of the denaturation dynamics. The PMF profiles reveal the thermodynamic feasibility of the DNA-CNT interaction, outlining distinct energy landscapes and barriers: when the selected swCNT binds within the dsDNA groove, the system becomes trapped at the first and second local energy minima, occurring at 1.48 nm and 1.00 nm, respectively. Intramolecular hydrogen bond calculations show a significant reduction, affirming the denaturing effect of swCNTs on DNA. Furthermore, the study reveals a significant reduction in the binding affinity of Ethidium Bromide (EB) to dsDNA following its interaction with swCNT, with a decrease in EB binding to dsDNA of approximately 13.2 %. This research offers valuable insights into the toxic effects of swCNTs on dsDNA, contributing to a rationalization of the cancerous potential of swCNTs.

Patterned graphene: An effective platform for adsorption, immobilization, and destruction of SARS-CoV-2 Mpro

To address the ongoing challenges posed by the SARS-CoV-2 and potentially stronger viruses in the future, the development of effective methods to fabricate patterned graphene (PG) and other precisely functional products has become a new research frontier. Herein, we modeled the "checkerboard" graphene (CG) and stripped graphene (SG) as representatives of PG, and studied their interaction mechanism with the target protein (Mpro) by molecular dynamics simulation. The calculation results on the binding strength and the root mean square deviation values of the active pocket revealed that PG is an effective platform for adsorption, immobilization, and destruction of Mpro. Specifically, CG is found to promote disruption of the active pocket for Mpro, but the presence of "checkerboard" oxidized regions inhibits the adsorption of Mpro. Meanwhile, the SG can effectively confine Mpro within the non-oxidized strips and enhances their binding strength, but doesn't play well on disrupting the active pocket. Our work not only elucidates the biological effects of PGs, but also provides guidance for their targeted and precise utilization in combating the SARS-CoV-2.

In Silico Investigation on the Molecular Behavior and Structural Stability of the Rosette Nanotubes as the Drug Vehicles for Paclitaxel, an Anti-Cancer Drug

Most anticancer drugs affect healthy cells in addition to cancer cells, causing severe side effects. Targeted delivery by nano-based drug delivery systems (NDDS) can reduce these severe side effects while maintaining therapeutic efficacy. This work introduced rosette nanotube (RNT) as a potential drug vehicle for paclitaxel (PTX) due to its self-assembling property, biocompatibility, amphiphilicity, and low toxicity. Molecular dynamics (MD) simulations aided with molecular mechanics Poisson Boltzmann surface area (MMPBSA) analysis are used here to investigate the molecular behavior and the loading energetics of each type of RNT (K1, xK1, and iEt-xK1) with PTX. Analysis showed that the most probable configuration of PTX is on either end of each RNT. The binding free energies (-117.74 to -69.29 kJ/mol) when PTX is closer to one end were stronger than when it is in the inner channel (-53.51 to -40.88 kJ/mol). The latter alludes to the encapsulation of the PTX by each RNT. Thus, loading is possible by encapsulation during the self-assembly process given the favorable estimated binding free energies. Based on the results, RNT has potential as a drug vehicle for PTX, which warrants further investigation.

Converting Nanotoxicity Data to Information Using Artificial Intelligence and Simulation

Decades of nanotoxicology research have generated extensive and diverse data sets. However, data is not equal to information. The question is how to extract critical information buried in vast data streams. Here we show that artificial intelligence (AI) and molecular simulation play key roles in transforming nanotoxicity data into critical information, i.e., constructing the quantitative nanostructure (physicochemical properties)-toxicity relationships, and elucidating the toxicity-related molecular mechanisms. For AI and molecular simulation to realize their full impacts in this mission, several obstacles must be overcome. These include the paucity of high-quality nanomaterials (NMs) and standardized nanotoxicity data, the lack of model-friendly databases, the scarcity of specific and universal nanodescriptors, and the inability to simulate NMs at realistic spatial and temporal scales. This review provides a comprehensive and representative, but not exhaustive, summary of the current capability gaps and tools required to fill these formidable gaps. Specifically, we discuss the applications of AI and molecular simulation, which can address the large-scale data challenge for nanotoxicology research. The need for model-friendly nanotoxicity databases, powerful nanodescriptors, new modeling approaches, molecular mechanism analysis, and design of the next-generation NMs are also critically discussed. Finally, we provide a perspective on future trends and challenges.

In-Situ Surface Modification of ITO Substrate via Bio-Inspired Mussel Chemistry for Organic Memory Devices

The development of organic memory devices, regarding factors such as structure construction, principle exploration, and material design, has become a powerful supplement to traditional silicon-based information storage. The in-situ growth of materials on substrate surfaces can achieve closer bonding between materials and electrodes. Bio-inspired by mussel chemistry, polydopamine (PDA) was self-assembled on a flexible substrate as a connecting layer, and 2-bromoiso-butyryl bromide (BiBB) was utilized as an initiator for the polymerization of an iridium complex via surface-initiated atom-transfer radical polymerization (SI-ATRP). A device with the structure of Al/PDA-PPy₃Ir/ITO was constructed after the deposition of aluminum. The device exhibited a nonvolatile rewritable memory characteristic with a turn-on voltage of -1.0 V and an ON/OFF current ratio of 6.3 × 10³. In addition, the memory performance of the Al/PDA-PPy₃Ir/ITO device remained stable at bending states due to the intrinsic flexibility of the active layer, which can be expanded into the establishment of flexible memory devices. Spectroscopy and electrochemical characterization suggested that the resistive memory properties of the device stemmed from charge transfer between PDA and iridium polymer in the active layer (PDA-PPy₃Ir) under an applied voltage.

Discovery of novel Glutaminase allosteric inhibitors through drug repurposing and comparative MMGB/PBSA and molecular dynamics simulation

GLS1 enzymes (Glutaminase C (GAC) and kidney-type Glutaminase (KGA)) are gaining prominence as a target for tumor treatment including lung, breast, kidney, prostate, and colorectal. To date, several medicinal chemistry studies are being conducted to develop new and effective inhibitors against GLS1 enzymes. Telaglenastat, a drug that targets the allosteric site of GLS1, has undergone clinical trials for the first time for the therapy of solid tumors and hematological malignancies. A comprehensive computational investigation is performed to get insights into the inhibition mechanism of the Telaglenastat. Some novel inhibitors are also proposed against GLS1 enzymes using the drug repurposing approach using 2D-fingerprinting virtual screening method against 2.4 million compounds, application of pharmacokinetics, Molecular Docking, and Molecular Dynamic (MD) Simulations. A TIP3P water box of 10 Å was defined to solvate both enzymes to improve MD simulation reliability. The dynamics results were validated further by the MMGB/PBSA binding free energy method, RDF, and AFD analysis. Results of these computational analysis revealed a stable binding affinity of Telaglenastat, as well as an FDA approved drug Astemizole (IC₅₀ ∼ 0.9 nM) and a novel para position oriented methoxy group containing Chembridge compound (Chem-64284604) that provides an effective inhibitory action against GAC and KGA.

The molecular origin of the electrostatic gating of single-molecule field-effect biosensors investigated by molecular dynamics simulations

Field-effect biosensors (bioFETs) offer a novel way to measure the kinetics of biomolecular events such as protein function and DNA hybridization at the single-molecule level on a wide range of time scales. These devices generate an electrical current whose fluctuations are correlated to the kinetics of the biomolecule under study. BioFETs are indeed highly sensitive to changes in the electrostatic potential (ESP) generated by the biomolecule. Here, using all-atom solvent explicit molecular dynamics simulations, we further investigate the molecular origin of the variation of this ESP for two prototypical cases of proteins or nucleic acids attached to a carbon nanotube bioFET: the function of the lysozyme protein and the hybridization of a 10-nt DNA sequence, as previously done experimentally. Our results show that the ESP changes significantly on the surface of the carbon nanotube as the state of these two biomolecules changes. More precisely, the ESP distributions calculated for these molecular states explain well the magnitude of the conductance fluctuations measured experimentally. The dependence of the ESP with salt concentration is found to agree with the reduced conductance fluctuations observed experimentally for the lysozyme, but to differ for the case of DNA, suggesting that other mechanisms might be at play in this case. Furthermore, we show that the carbon nanotube does not impact significantly the structural stability of the lysozyme, corroborating that the kinetic rates measured using bioFETs are similar to those measured by other techniques. For DNA, we find that the structural ensemble of the single-stranded DNA is significantly impacted by the presence of the nanotube, which, combined with the ESP analysis, suggests a stronger DNA-device interplay. Overall, our simulations strengthen the comprehension of the inner working of field-effect biosensors used for single-molecule kinetics measurements on proteins and nucleic acids.

The Inhibition of SARS-CoV-2 3CL Mpro by Graphene and Its Derivatives from Molecular Dynamics Simulations

At present, the most powerful new drugs for COVID-19 are antibody proteins. In addition, there are some star small molecule drugs. However, there are few studies on nanomaterials. Here, we study the intact graphene (IG), defective graphene (DG), and graphene oxide (GO) interacting with COVID-19 protein. We find that they show progressive inhibition of COVID-19 protein. By using molecular dynamics simulations, we study the interactions between SARS-CoV-2 3CL Mpro and graphene-related materials (GRMs): IG, DG, and GO. The results show that Mpro can be absorbed onto the surfaces of investigated materials. DG and GO interacted with Mpro more intensely, causing the decisive part of Mpro to become more flexible. Further analysis shows that compared to IG and GO, DG can inactivate Mpro and inhibit its expression effectively by destroying the active pocket of Mpro. Our work not only provides detailed and reliable theoretical guidance for the application of GRMs in treating with SARS-CoV-2 but also helps in developing new graphene-based anti-COVID-19 materials.

Tracking the interaction between single-wall carbon nanotube and SARS-Cov-2 spike glycoprotein: A molecular dynamics simulations study

COVID-19, a newly discovered type of coronavirus, is the cause of the pandemic infection that was first reported in Wuhan, China, in December 2019. One of the most critical problems in this regard is to identify innovative drugs that may reduce or manage this global health concern. Nanoparticles have shown a pivotal role in drug delivery systems in recent decades. The surface of nanoparticles could be covered by a layer composed of different biomolecules (e.g., proteins and macromolecules) following the incubation with a biological fluid. This protein-rich layer is called “Protein Corona.” In this study, an all-atom molecular dynamics simulation was used for investigating the monomeric B domain of the spike glycoprotein due to its role in the accessibility of the spike glycoprotein to single-wall carbon nanotubes (SWCNTs). The interaction energy values between the carbon nanotube and B domain of the viral spike glycoprotein were evaluated. The obtained results, based on Lennard-Jones potentials, demonstrated that SWCNTs had an affinity to the B domain of the S1 subunit in the spike glycoprotein. The adsorption of SWCNTs on the B domain surface led to a significant change in solvent-accessible surface, internal hydrogen bonds, and finally in the tertiary structure, which could provide a reasonable method to impede the interaction between the angiotensin-converting enzyme II and SARS-CoV-2 spike glycoprotein. A decrease in the mean square displacement of the B domain was shown after the adsorption of SWCNTs as a result of increasing the hydrophobic-hydrophilic properties of the B domain. The arrangement of SWCNTs on the B domain surface and their interaction using the 2-acetamido-2-deoxy-β-d-glucopyranose group (988, 991, and 992) demonstrated that a change in the affinity of the S1 subunit could be used as a barrier to viral replication. The analysis of the SWCNT-B domain complex indicated that the presence of SWCNTs is able to cause alterations in the S1 subunit of the spike protein, and these nanotubes could be employed for further in-vitro and in-vivo antiviral studies. Also, SWCNTs are able to be utilized in drug delivery systems.

How do proteins 'response' to common carbon nanomaterials?

Carbon nanomaterials are widely produced and applied in biological and environmental fields because of their outstanding physical and chemical properties, which pose a threat to the safety of living organisms and the ecological environment. Therefore, understanding how carbon nanomaterials and their derivatives work on organisms is becoming important. In recent years, more and more researchers have explored the damage of carbon nanomaterials to organisms at the molecular level. This review pays special emphasis on how proteins response to the main carbon nanomaterials (fullerene, carbon nanotubes, graphene and their derivatives). In addition, how to use the interaction between carbon nanomaterials and proteins to do some beneficial things for human and the development of safe nanomaterials is simply discussed. Finally, some suggestions have been made to lay a theoretical foundation for future research.

Ordered hydration layer mediated ice adsorption of a globular antifreeze protein: mechanistic insight

The ice binding surface of a type III AFP induces water ordering at lower temperature, which mediates its adsorption on the ice surface.

Translocation of Bioactive Molecules through Carbon Nanotubes Embedded in the Lipid Membrane

One of the major challenges of nanomedicine and gene therapy is the effective translocation of drugs and genes across cell membranes. In this study, we describe a systematic procedure that could be useful for efficient drug and gene delivery into the cell. Using fully atomistic molecular dynamics (MD) simulations, we show that molecules of various shapes, sizes, and chemistries can be spontaneously encapsulated in a single-walled carbon nanotube (SWCNT) embedded in a 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) lipid bilayer, as we have exemplified with dendrimers, asiRNA, ssDNA, and ubiquitin protein. We compute the free energy gain by the molecules upon their entry inside the SWCNT channel to quantify the stability of these molecules inside the channel as well as to understand the spontaneity of the process. The free energy profiles suggest that all molecules can enter the channel without facing any energy barrier but experience a strong energy barrier (≫k_BT) to translocate across the channel. We propose a theoretical model for the estimation of encapsulation and translocation times of the molecules. Whereas the model predicts the encapsulation time to be of the order of few nanoseconds, which match reasonably well with those obtained from the simulations, it predicts the translocation time to be astronomically large for each molecule considered in this study. This eliminates the possibility of passive diffusion of the molecules through the CNT-nanopore spanning across the membrane. To counter this, we put forward a mechanical method of ejecting the encapsulated molecules by pushing them with other free-floating SWCNTs of diameter smaller than the pore diameter. The feasibility of the proposed method is also demonstrated by performing MD simulations. The generic strategy described here should work for other molecules as well and hence could be potentially useful for drug- and gene-delivery applications.

Fullerene derivatives act as inhibitors of leukocyte common antigen based on molecular dynamics simulations

Fullerene-based molecules are being studied as potential inhibitors of protein tyrosine phosphatases due to their unique properties and low toxicity. However, the underlying molecular mechanism remains elusive. In this study, molecular dynamics (MD) simulations in conjunction with molecular docking calculations were utilized to investigate the binding effects of C₆₀, C₆₀(NH₂)₃₀, and C₆₀(OH)₃₀ on the enzymatic activity of CD45 (a receptor-like protein tyrosine phosphatase). Our results show that all the investigated molecules can be docked into the region between D1 and D2 domains of CD45, and stabilize the protein structure. The average number of residues that directly interact with the C₆₀(NH₂)₃₀ is two more than that of C₆₀(OH)₃₀, F819 and F820 (located in the loop connects α3 and β12), resulting in different effects of C₆₀(NH₂)₃₀ and C₆₀(OH)₃₀ on protein activity. Detailed MD simulation analyses show that transformation of the interaction network caused by C₆₀(NH₂)₃₀ is completely different from that of the control simulation due to the misfolding of α3. Furthermore, the movement of D1 active pocket and KNRY motif are most severely impaired by docking with C₆₀(NH₂)₃₀. Our simulation results illustrate that fullerene derivatives modified with amino groups exhibit conspicuous tumor inhibition to protein tyrosine phosphatases, and can act as effective inhibitors. Our results give insight into the inhibitory effects of fullerene-based molecules on protein tyrosine phosphatases and providing a theoretical basis for the design of effective inhibitors.

Fullerene-based molecules are being studied as potential inhibitors of protein tyrosine phosphatases due to their unique properties and low toxicity.