Investigation of glutathione as a natural antioxidant and multitarget inhibitor for Alzheimer’s disease: Insights from molecular simulations

Unveiling the dynamic and thermodynamic interactions of hydrocortisone with β-cyclodextrin and its methylated derivatives through insights from molecular dynamics simulations

Cyclodextrins (CDs) can enhance the stability and bioavailability of pharmaceutical compounds by encapsulating them within their cavities. This study utilized molecular dynamics simulations to investigate the interaction mechanisms between hydrocortisone (HC) and various methylated CD derivatives. The results reveal that the loading of HC into CD cavities follows different mechanisms depending on the degree and position of methylation. Loading into βCD and 6-MeβCD was more complete, with the hydroxyl groups of HC facing the primary hydroxyl rim (PHR) and the ketone side facing the secondary hydroxyl rim (SHR). In contrast, 2,3-D-MeβCD and 2,6-D-MeβCD showed a different loading mechanism, with the ketone side facing the PHR and the hydroxyl groups facing the SHR. The root mean square fluctuation (RMSF) analysis demonstrated that methylation increases the flexibility of CD heavy atoms, with 3-MeβCD and 2,3-D-MeβCD exhibiting the highest flexibility. However, upon inclusion of HC, 3-MeβCD, 2,3-D-MeβCD, 2-MeβCD, and 6-MeβCD showed a significant reduction in flexibility, suggesting a more rigid structure that effectively retains HC within their cavities. The radial distribution function revealed a significant reduction in the number of water molecules within the innermost layer of the methylated CD cavities, particularly in TMeβCD, indicating a decrease in polarity. The presence of HC led to the release of high-energy water molecules, creating more favorable conditions for HC loading. Conformational analysis showed that methylation caused a partial decrease in the area of the PHR, a significant decrease in the area of the middle rim, and a notable decrease in the area of the SHR. The loading of HC increased the area of the PHR in most derivatives, with the most pronounced increase observed in 2,6-D-MeβCD and 6-MeβCD. The analysis of interaction energies and binding free energies demonstrated that the binding of HC to methylated CD derivatives is thermodynamically more favorable than to βCD, with the strongest association observed for 6-MeβCD, 2-MeβCD, and 2,3-D-MeβCD.

Antioxidant properties of α-amino acids: a density functional theory viewpoint

The antioxidant properties of 21 proteinogenic amino acids (AAs) and 3,4-dioxophenylanine (DOPA) have been studied in implicit water using density functional theory (DFT). All the calculations have been performed according to three oxidation mechanisms: (1) hydrogen-atom transfer (HAT); (2) single electron transfer followed by proton transfer (SET-PT); and (3) sequential proton-loss electron transfer (SPLET). As a result, five AAs with the highest antioxidant capacity have been established: DOPA, selenocysteine (Sec), tyrosine (Tyr), cysteine (Cys), and tryptophan (Trp). Also, global reactivity in terms of hardness/softness has been evaluated, as well as Fukui indices of local reactivity. Trp has been determined as the most reactive molecule, whereas selenium atom of Sec has been established as the most reactive atom. All the findings are in agreement with the recent literature on both experimental and theoretical studies of amino acids antioxidant activity. However, to the best of my knowledge, the calculations for one electron redox reactions of zwitterionic amino acids in implicit water have been performed for the first time.

De novo antioxidant peptide design via machine learning and DFT studies

Antioxidant peptides (AOPs) are highly valued in food and pharmaceutical industries due to their significant role in human function. This study introduces a novel approach to identifying robust AOPs using a deep generative model based on sequence representation. Through filtration with a deep-learning classification model and subsequent clustering via the Butina cluster algorithm, twelve peptides (GP1-GP12) with potential antioxidant capacity were predicted. Density functional theory (DFT) calculations guided the selection of six peptides for synthesis and biological experiments. Molecular orbital representations revealed that the HOMO for these peptides is primarily localized on the indole segment, underscoring its pivotal role in antioxidant activity. All six synthesized peptides exhibited antioxidant activity in the DPPH assay, while the hydroxyl radical test showed suboptimal results. A hemolysis assay confirmed the non-hemolytic nature of the generated peptides. Additionally, an in silico investigation explored the potential inhibitory interaction between the peptides and the Keap1 protein. Analysis revealed that ligands GP3, GP4, and GP12 induced significant structural changes in proteins, affecting their stability and flexibility. These findings highlight the capability of machine learning approaches in generating novel antioxidant peptides.

Exploring the role of cyclodextrins as a cholesterol scavenger: a molecular dynamics investigation of conformational changes and thermodynamics

This study presents a comprehensive analysis of the cholesterol binding mechanism and conformational changes in cyclodextrin (CD) carriers, namely βCD, 2HPβCD, and MβCD. The results revealed that the binding of cholesterol to CDs was spontaneous and thermodynamically favorable, with van der Waals interactions playing a dominant role, while Coulombic interactions have a negligible contribution. The solubility of cholesterol/βCD and cholesterol/MβCD complexes was lower compared to cholesterol/2HPβCD complex due to stronger vdW and Coulombic repulsion between water and CDs. Hydrogen bonding was found to have a minor role in the binding process. The investigation of mechanisms and kinetics of binding demonstrated that cholesterol permeates into the CD cavities completely. Replicas consideration indicated that while the binding to 2HPβCD occurred perpendicularly and solely through positioning cholesterol's oxygen toward the primary hydroxyl rim (PHR), the mechanism of cholesterol binding to βCD and MβCD could take place with the orientation of oxygen towards both rims. Functionalization resulted in decreased cavity polarity, increased constriction tendency, and altered solubility and configuration of the carrier. Upon cholesterol binding, the CDs expanded, increasing the cavity volume in cholesterol-containing systems. The effects of cholesterol on the relative shape anisotropy (κ2) and asphericity parameter (b) in cyclodextrins were investigated. βCD exhibited a spherical structure regardless of cholesterol presence, while 2HPβCD and MβCD displayed more pronounced non-sphericity in the absence of cholesterol. Loading cholesterol transformed 2HPβCD and MβCD into more spherical shapes, with increased probabilities of higher κ2. MβCD showed a higher maximum peak of κ2 compared to 2HPβCD after cholesterol loading, while 2HPβCD maintained a significant maximum peak at 0.2 for b.

Elucidating anticancer drugs release from UiO-66 as a carrier through the computational approaches

The computational analysis of drug release from metal-organic frameworks (MOFs), specifically UiO-66, is the primary focus of this research. MOFs are recognized as nanocarriers due to their crystalline structure, porosity, and potential for added functionalities. The research examines the release patterns of three drugs: temozolomide, alendronate, and 5-fluorouracil, assessing various factors such as the drugs' distance from the UiO-66 centers, the interaction of drug functional groups with Zr metal ions, and the drug density throughout the nanocarrier. Findings reveal that 5-fluorouracil is located furthest from the UiO-66 center and exhibits the highest positive energy compared to the other drugs. Alendronate's density is observed to shift to the carrier surface, while 5-fluorouracil's density significantly decreases within the system. The drug density diminishes as the distance from the UiO-66 center of mass increases, suggesting a stronger positive interaction between the drugs and the nanocarrier. Moreover, Monte Carlo calculations were employed to load drugs onto the UiO-66 surface, leading to a substantial release of 5-fluorouracil from UiO-66. Quantum and Monte Carlo adsorption localization calculations were also conducted to gather data on the compounds' energy and geometry. This research underscores the potential of MOFs as nanocarriers for drug delivery and highlights the crucial role of temperature in regulating drug release from UiO-66. It provides insights into the complex dynamics of drug release and the factors influencing it, thereby emphasizing the promise of UiO-66 as a viable candidate for drug delivery. This work contributes to our understanding of UiO-66's role and sets the stage for improved performance optimization in the cancer treatment.

Insights from molecular dynamics and DFT calculations into the interaction of 1,4-benzodiazepines with 2-hydroxypropyl-βCD in a theoretical study

This study delves into the interaction between benzodiazepine (BZD) drugs and 2-hydroxypropyl-β-cyclodextrin (2HPβCD), a cyclodextrin (CD) known to improve drug delivery and enhance therapeutic outcomes. We find that the 2HPβCD's atoms become more rigid in the presence of chlordiazepoxide (CDP), clonazepam (CLZ), and diazepam (DZM), whereas they become more flexible in the presence of nordazepam (NDM) and nitrazepam (NZP). We also investigated the structure of 2HPβCD and found that loading these drugs increases both the area and volume of the 2HPβCD cavity, making it more suitable for drug delivery. Moreover, this research found that all drugs exhibited negative values for the binding free energy, indicating thermodynamic favorability and improved solubility. The binding free energy order of the BZDs was consistent in both molecular dynamics and Monte Carlo methods, with CDP and DZM having the highest affinity for binding. We also analyzed the contribution of different interaction energies in binding between the carrier and the drugs and found that Van der Waals energy is the primary component. Our results indicate that the number of hydrogen bonds between 2HPβCD/water slightly decreases in the presence of BZDs, but the hydrogen bond's quality remains constant.

Interactions between β-cyclodextrin as a carrier for anti-cancer drug delivery: a molecular dynamics simulation study

A series of molecular dynamics simulations were performed on 5-fluorouracil (5-Fu), Alendronate (Ald), and Temozolomide (TMZ) anticancer drugs in the presence and absence of β-cyclodextrin (βCD) as a carrier. Thermodynamic investigations showed that the van der Waals interaction energy was dominant in loading all drugs inside the βCD cavity. The sum of the interaction energies illustrated that the highest affinity was related to Ald (-136.5 kJ/mol), which in turn was due to the presence of bulky and charged atoms of phosphorus and oxygen, although TMZ (-115.92 kJ/mol) showed a very high affinity as well. At the same time, the hydrogen bond analysis also represented that Ald had the most hydrogen bond (1.97) with the highest half-life (3.13 ps) with βCD. Investigation of the root mean fluctuation (RMSF) indicated that all the drugs had a relatively rigid structure and maintain this rigidity during loading in the βCD cavity, and in the meantime, Ald was slightly more flexible than 5-Fu and TMZ. The area of ​the primary hydroxyl rim decreased in all drug-containing systems, which in turn was caused by the attractive interaction of drugs with oxygens in the primary hydroxyl rim. Especially for those drugs that were able to penetrate to the end of the primary hydroxyl rim of the βCD, that means TMZ and 5-Fu. Meanwhile, due to the lack of Ald penetration to the end of the primary hydroxyl rim, the area change in the Ald-containing system was less than in the two others.Communicated by Ramaswamy H. Sarma.

Machine-learning-assisted molecular design of phenylnaphthylamine-type antioxidants

In this study, a total of 302 molecular structures of phenylnaphthylamine antioxidants based on N-phenyl-1-naphthylamine and N-phenyl-2-naphthylamine skeletons with various substituents were modeled by exhaustive methods. Antioxidant parameters, including the hydrogen dissociation energy, solubility parameter, and binding energy, were calculated through molecular simulations. Then, a group decomposition scheme was determined to decompose 302 antioxidants. The antioxidant parameters and decomposition results constituted machine-learning data sets. Using an artificial neural network model, a correlation coefficient between the predicted and true values above 0.88 and an average relative error within 6% were achieved. Random forest models were used to analyze the factors affecting antioxidant activity from chemical and physical perspectives; the results showed that amino and alkyl groups were conducive to improving antioxidant performance. Moreover, substituent positions 1, 7, and 10 of N-phenyl-1-naphthylamine and 3, 7, and 10 of N-phenyl-2-naphthylamine were found to be the optimal positions for modifications to improve antioxidant activity. Two potentially efficient phenylnaphthylamine antioxidant structures were proposed and their antioxidant parameters were also calculated; the hydrogen dissociation energy and solubility parameter decreased by more than 9% and 7%, respectively, whereas the binding energy increased by more than 16% compared with the benchmark of N-phenyl-1-naphthylamine. These results indicate that molecular simulation and machine learning could provide alternative tools for the molecular design of new antioxidants.

In situ synthesis of chiral AuNCs with aggregation-induced emission using glutathione and ceria precursor nanosheets for glutathione biosensing

Au(i)–SG/Ce(CO3)2 NS conjugated nanoprobe was developed for selective GSH detection. The redox reaction between GSH and the NS could release Ce3+ ions to initiate the intense AIE of Au(i)–SG oligomers.