A molecular dynamics simulation study on the conformational stability of amylose-linoleic acid complex in water

Combined role of stearic acid and maleic anhydride in the development of thermoplastic starch-based materials with ultrahigh ductility and durability

The diverse properties reported for starch-based materials indicate their potential for use in the preparation of biodegradable flexible actuators. However, their natural brittleness and lack of durability after modification limit their practical application. Therefore, we propose a strategy for preparing flexible starch-based composites. The results of macro/micro property characterizations and molecular dynamics simulations indicated that using starch, maleic anhydride, and stearic acid (SA), the mobility of the starch chains was enhanced and retrogradation was inhibited through the synergistic effects induced by chain breaking, complex formation with SA, and esterification of the starch molecules. In addition, the elongation at break of the modified starch (MS) reached 2070 %, and considerable ductility (>1000 %) as well as well-complexed structure were maintained after six months. Furthermore, the MS was able to undergo self-healing after fracture or a temperature-controlled stiffness transition. Moreover, it underwent complete degradation in soil within 30 d. Finally, an actuator was prepared by doping the MS with nano-Fe3O4 particles to realize a dual magnetic and optical response. Dynamic monitoring was also achieved based on the electrical signal, thereby demonstrating the broad application scope of this material in the development of biodegradable flexible actuators.

Investigation of the effect of ultrasonication on starch-fatty acid complexes and the stabilization mechanism

The complexation of physically modified starch with fatty acids is favorable for the production of resistant starch. However, there is a lack of information on the effect of ultrasonication (UC) on the structure and properties of starch complexes and the molecular mechanism of the stabilization. Here, the multi-scale structure and in vitro digestive properties of starch-fatty acid complexes before and after UC were investigated, and the stabilization mechanisms of starch and fatty acids were explored. The results showed that the physicochemical properties and multi-scale structure of the starch-fatty acid complexes significantly changed with the type of fatty acids. The solubility and swelling power of the starch-fatty acid complexes were significantly decreased after UC (P < 0.05), which facilitated the binding of starch with fatty acids. The XRD results revealed that after the addition of fatty acids, the starch-fatty acid complexes showed typical V-shaped complexes. In addition, the starch-fatty acid complexes showed a significant increase in complexing index, improved short-range ordering and enhanced thermal stability. However, the differences in the structure and properties of the fatty acids themselves resulted in no significant improvement in the multi-scale structure of maize starch-palmitic acid by UC. In terms of digestibility, especially the complexes after UC were more compact in structure, which increased the difficulty of enzymatic digestion and thus slowed down the digestion process. DFT calculations and combined with FT-IR analysis showed that non-covalent interactions such as hydrogen bonding and hydrophobic interactions were the main driving force for the formation of the complexes, with binding energies (lauric acid, myristic acid and palmitic acid) of -30.50, -22.14 and -14.10 kcal/mol, respectively. Molecular dynamics simulations further confirmed the molecular mechanism of inclusion complex formation and stabilization. This study is important for the regulation of starchy foods by controlling processing conditions, and provides important information on the role of fatty acids in the regulation of starch complexes and the binding mechanism.

A new insight into the key matrix components for aftertaste in Ampelopsis grossedentata (vine tea) infusion: From the intensity and duration of taste profiles using non-targeted metabolomics and molecular simulation

The aftertaste with a prolonged duration in ampelopsis grossedentata infusion (AGTI) is easily perceived, however, its formation mechanism is unclear. Therefore, aftertaste-A and richness were confirmed as the characteristic aftertaste of AGTI through sensory evaluation and electronic tongue. Moreover, 5-KETE, theobromine, etc., metabolites were identified as the differential components between AGTI and green tea infusion. Among them, p-coumaroyl quinic acid, xanthine etc., and proline, dihydromyricetin, etc., components contributed more to the formation of aftertaste-A and richness, respectively. Further, the bonding between characteristic metabolites for aftertaste in AGTI with their receptors were shown to be more stable using molecular docking, compared to metabolites related to typical taste profiles. The aftertaste in AGTI was more easily perceived by saltiness components or in NaCl system by molecular simulation. This study offers novel insight into the interaction mechanism of aftertaste in tea infusion and will contribute to further study on aftertaste for other foods.

Water/oil interfacial behaviors of soy hull polysaccharide revealed by molecular dynamics simulation and particle tracking microrheology

The soy hull polysaccharide (SHP) exhibits excellent interfacial activity and holds potential as an emulsifier for emulsions. To reveal the behavior of SHP at the water/oil (W/O) interface in situ, molecular dynamics (MD) simulations and particle tracking microrheology were used in this study. The results of MD reveal that SHP molecular spontaneously move toward the interface and rhamnogalacturonan-I initiates this movement, while its galacturonic acids on it act as anchors to immobilize the SHP molecules at the W/O interface. Microrheology results suggest that SHP forms microgels at the W/O interface, with the lattices of the microgels continually undergoing dynamic changes. At low concentrations of SHP and short interfacial formation time, the network of the microgels is weak and dominated by viscous properties. However, when SHP reaches 0.75 % and the interfacial formation time is about 60 min, the microgels show perfect elasticity, which is beneficial for stabilizing emulsions.

Application of the molecular dynamics simulation GROMACS in food science

Food comprises proteins, lipids, sugars and various other molecules that constitute a multicomponent biological system. It is challenging to investigate microscopic changes in food systems solely by performing conventional experiments. Molecular dynamics (MD) simulation serves as a crucial bridge in addressing this research gap. The Groningen Machine for Chemical Simulations (GROMACS) is an open-source, high-performing molecular dynamics simulation software that plays a significant role in food science research owing to its high flexibility and powerful functionality; it has been used to explore the molecular conformations and the mechanisms of interaction between food molecules at the microcosmic level and to analyze their properties and functions. This review presents the workflow of the GROMACS software and emphasizes the recent developments and achievements in its applications in food science research, thus providing important theoretical guidance and technical support for obtaining an in-depth understanding of the properties and functions of food.

Computational modelling of extrusion process temperatures on the interactions between black soldier fly larvae protein and corn flour starch

Insects such as the black soldier fly (BSF) are recently being studied as food sources to address concerns about how to meet the food demand of the growing world population, as conventional production lines for meat proteins are currently unsustainable sources. Studies have been conducted evaluating the use of insect proteins to produce extruded foods such as expanded snacks and meat analogues. However, this field of study is still quite new and not much has been studied beyond digestibility and growth performance. The purpose of this work was to evaluate the compatibility of protein extracted from BSF flour with corn flour starch within an extruded balanced shrimp feed model through molecular dynamics simulations, for which cohesive energy density and solubility parameter (δ) of both components were determined. The calculations' results for the protein molecule systems yielded an average δ of 14.961 MPa0.5, while the δ for starch was calculated to be 23.166 MPa0.5. The range of difference between both δ (10 > δ > 7) suggests that the interaction of the BSF protein with corn starch is of a semi-miscible nature. These results suggest that it is possible to obtain a stable starch-protein mixture through the extrusion process.

Study on multiscale structures and digestibility of cassava starch and medium-chain fatty acids complexes using molecular simulation techniques

Effect of complexation of three medium-chain fatty acids (octanoic, decylic and lauric acid, OA, DA and LA, respectively) on structural characteristics, physicochemical properties and digestion behaviors of cassava starch (CS) was investigated. Current study indicated that LA was more easily to combine with CS (complex index 88.9%), followed by DA (80.9%), which was also consistent with their corresponding complexed lipids content. Following the investigation of morphology, short-range ordered structure, helical structure, crystalline/amorphous region and fractal dimension of the various complexes, all cassava starch-fatty acids complexes (CS-FAs) were characterized with a flaked morphology rather than a round morphology in native starch (control CS). X-ray diffraction demonstrated that all CS-FAs had a V-type crystalline structure, and nuclear magnetic resonance spectroscopy confirmed that the complexes made from different fatty acids displayed similar V6 or V7 type polymorphs. Interestingly, small-angle X-ray scattering analysis revealed that α value became greater following increased carbon chain length of fatty acids, indicating the formation of a more ordered fractal structure in the aggregates. Changes in rheological parameters G' and G'' indicated that starch complexed with fatty acids was more likely to form a gel network, but difference among three CS-FAs complexes was significant, which might be contributed to their corresponding hydrophobicity and hydrophilicity raised from individual fatty acids. Importantly, digestion indicated that CS-LA complexes had the lowest hydrolysis degree, followed by the greatest RS content, indicating the importance of chain length of fatty acids for manipulating the fine structure and functionality of the complexes.

Analysis of the relationship between starch molecular conformation and enzymatic hydrolysis efficiency

Resistant starch (RS) is important in controlling diabetes. The primary objective of this study is to examine the impact of molecular conformation on the enzymatic hydrolysis efficiency of starch by α-amylase. And the interactions between starch molecules with different conformations and α-amylase were analysed by using molecule dynamics simulation and molecular docking. It was found, the natural conformational starch molecule was hydrolysed from the middle of the starch chain by α-amylase, producing polysaccharides. The bent PS-conformational starch molecules with multiple O2-O3 intramolecular hydrogen bonds produced by high-pressure was hydrolysed from the head of the starch chain to produce glucose, which is not conducive to RS formation. The stretched H-conformation without intramolecular hydrogen bonds produced by heat treatment was not hydrolysed by α-amylase. However, it occupied the active groove and formed strong interactions with α-amylase, which prevented other starch molecules from binding to α-amylase, thus reducing hydrolysis efficiency. Moreover, the total interaction energies between the three starch molecules and α-amylase were approximately 78 kJ/mol. And several hydrogen bonds were formed between the starch molecules and α-amylase, which provides evidence for the continuous sliding hydrolysis hypothesis of α-amylase. Moreover, these results provide an important reference for elucidating the mechanism of RS formation.

Evaluating solubility, stability, and inclusion complexation of oxyresveratrol with various β-cyclodextrin derivatives using advanced computational techniques and experimental validation

Oxyresveratrol (OXY), a natural stilbenoid in mulberry fruits, is known for its diverse pharmacological properties. However, its clinical use is hindered by low water solubility and limited bioavailability. In the present study, the inclusion complexes of OXY with β-cyclodextrin (βCD) and its three analogs, dimethyl-β-cyclodextrin (DMβCD), hydroxypropyl-β-cyclodextrin (HPβCD) and sulfobutylether-β-cyclodextrin (SBEβCD), were investigated using in silico and in vitro studies. Molecular docking revealed two binding orientations of OXY, namely, 4',6'-dihydroxyphenyl (A-form) and 5,7-benzenediol ring (B-form). Molecular Dynamics simulations suggested the formation of inclusion complexes with βCDs through two distinct orientations, with OXY/SBEβCD exhibiting maximum atom contacts and the lowest solvent-exposed area in the hydrophobic cavity. These results corresponded well with the highest binding affinity observed in OXY/SBEβCD when assessed using the MM/GBSA method. Beyond traditional simulation methods, Ligand-binding Parallel Cascade Selection Molecular Dynamics method was employed to investigate how the drug enters and accommodates within the hydrophobic cavity. The in silico results aligned with stability constants: SBEβCD (2060 M-1), HPβCD (1860 M-1), DMβCD (1700 M-1), and βCD (1420 M-1). All complexes exhibited a 1:1 binding mode (AL type), with SBEβCD enhancing OXY solubility (25-fold). SEM micrographs, DSC thermograms, FT-IR and 1H NMR spectra confirm the inclusion complex formation, revealing novel surface morphologies, distinctive thermal behaviors, and new peaks. Notably, the inhibitory impact on the proliferation of breast cancer cell lines, MCF-7, exhibited by inclusion complexes particularly OXY/DMβCD, OXY/HPβCD, and OXY/SBEβCD were markedly superior compared to that of OXY alone.

Investigation of the formation mechanism of the pepper starch-piperine complex

In this study, a new carrier for loading piperine was prepared using pepper starch, and its interaction mechanism was investigated. The porous pepper starch-piperine complex (PPS-PIP) showed higher loading efficiency (76.15 %) compared to the porous corn starch-piperine complex (PCS-PIP (52.34 %)). This may be ascribed to the hemispherical shell structure of porous pepper starch (PPS) compared to the porous structure of porous corn starch (PCS) based on the SEM result. PPS-PIP had smaller particle size (10.53 μm), higher relative crystallinity (38.95 %), and better thermal stability (87.45 °C) than PCS-PIP (17.37 μm, 32.17 %, 74.35 °C). Fourier transform infrared spectroscopy (FTIR) results implied that piperine not only forms a complex with amylose but may also be physically present in porous starch. This was demonstrated by the short-range order and X-ray type. Molecular dynamics simulations confirmed that hydrogen bonding is the primary interaction between amylose and piperine. Besides the formation of the amylose-piperine complex, some of the piperine is also present in physical form.

Insights into formation and stability mechanism of V7-type short amylose-resveratrol complex using molecular dynamics simulation and molecular docking

Pre-formed V-type amylose as a kind of wall material has been reported to carry polyphenols, while the interaction mechanism between V-type amylose and polyphenol is still elusive. In this work, the formation and stability mechanism of a V7-type short amylose-resveratrol complex was investigated via isothermal titration calorimetry, molecular dynamics, and molecular docking. The results presented that two stoichiometric ratios of resveratrol to short amylose were calculated to 0.120 and 0.800, and the corresponding main driving force was hydrogen bonding and hydrophobic interaction, respectively. The folding and unfolding conformation of V7-type short amylose chains appeared alternately during the simulation. Resveratrol tended to be bound in the short amylose helix between 40 ns and 80 ns to form a more stable complex. Hydrogen bonds between resveratrol molecule and O6 at the 22nd glucose molecule/O2 at the 24th glucose molecules and hydrophobic interaction between resveratrol molecule and glucose molecules (19th, 20th, 21st and 23rd) could be found.

Multiscale structure changes and mechanism of polyphenol-amylose complexes modulated by polyphenolic structures

This study was designed to investigate the effect of polyphenolic structure on the interaction strength and process between polyphenols (gallic acid (GA), epigallocatechin gallate (EGCG) and tannic acid (TA)) and amylose (AM). The results of Fourier transform infrared spectroscopy, isothermal titration calorimetry, X-ray photoelectron spectroscopy and molecular dynamic simulation (MD) suggested that the interactions between the three polyphenols and AM were noncovalent, spontaneous, low-energy and driven by enthalpy, which would be enhanced with increasing amounts of pyrogallol groups in the polyphenols. The results of turbidity, particle size and appearance of the complex solution showed that the interaction process between polyphenols and AM could be divided into three steps and would be advanced by increasing the number of pyrogallol groups in the polyphenols. At the same time, MD was intuitively employed to exhibit the interaction process between amylose and polyphenols, and it revealed that the interaction induced the aggregation of amylose and that the agglomeration degree of amylose increased with increasing number of pyrogallol groups at polyphenols. Last, the SEM and TGA results showed that TA/AM complexes had the tightest structure and the highest thermal stability (TA/AM˃EGCG/AM˃GA/AM), which could be attributed to TA having five pyrogallol groups.

Gallic acid forms V-amylose complex structure with starch through hydrophobic interaction

This study aimed to investigate the role of amylose and amylopectin in the formation of starch-polyphenol complex and elucidate the interaction mechanisms. Gallic acid (GA) was used to complex with maize starch with various amylose contents. Results showed GA formed V-type crystals with normal maize starch (NMS) and high amylose maize starch (HAMS), while higher relative crystallinity was exhibited in HAMS-GA complexes than NMS counterparts. Molecular structure analysis revealed more amylose in GA-starch complexes than in treated starch counterparts without GA, and this was more apparent in HAMS than NMS, implying amylose is preferred to complex with GA than amylopectin. FTIR detected higher R1047/1022 value in starch-GA complexes than their starch counterparts without GA, suggesting increased short-range ordered structrure of complexes. Typical signatures of hydrophobic interactions were further revealed by isothermal titration calorimetry, indicating the complexation of GA to starch is mainly through hydrophobic bonds. More binding sites were observed for HAMS (72.50) than NMS (11.33), which proves the preferences of amylose to bind with GA. Molecular dynamics simulated the complexation of GA to amylose, and confirmed hydrophobic bond is the main interaction force. These findings would provide guidance for precise design and utilization of starch-polyphenol complexes in functional foods.

Optimizing biodegradable plastics: Molecular dynamics insights into starch plasticization with glycerol and oleic acid

Petroleum-based plastics dominate everyday life, necessitating the exploration of natural polymers as alternatives. Starch, abundant and biodegradable, is a promising raw material. However, understanding the molecular mechanisms underlying starch plasticization has proven challenging. To address this, we employ molecular dynamics simulations, focusing on amylose as a model. Our comprehensive evaluation revealed that chain size affects solubility, temperature influenced diffusivity and elastic properties, and oleic acid expressed potential as an alternative plasticizer. Furthermore, blending glycerol or oleic acid with water suggested the enhancement amylose's elasticity. These findings contribute to the design of sustainable and improved biodegradable plastics.

Characterization, in vitro digestibility and release properties of starch-linoleic acid-sodium alginate composite film

This study aimed to improve the complexing degree, digestibility and controlled release properties of the potato starch (PS)-linoleic acid (LA) complexes by encapsulating PS-LA complexes to sodium alginate (AG) beads. The results revealed that AG had a positive effect on the complexing index, R1047/1022 values, relative crystallinity, enthalpy and morphological structure of PS-LA-AG films, especially for PS-LA-AG film with the PS-LA: AG of 5:1. The in vitro digestion and hydrolysis kinetic analysis indicated that AG addition reduced the digestibility of PS-LA-AG films to a higher slowly digestible starch content and resistant starch content and a lower equilibrium hydrolysis percentage and kinetic constant. Furthermore, in vivo release study of PS-LA-AG films indicated a restrained release in simulated gastrointestinal conditions. Consequently, the results indicated that AG addition significantly improved the inclusion efficiency for the complex formation between PS and LA, which was beneficial for the design of resistant films to entrap and control release of unsaturated fatty.

Biological regulation on iodine using nano-starch for preventing thyroid dysfunction

The growing incidence of thyroid disease triggered by excess iodine uptake poses a severe health threat throughout the world. Extracellular interference therapies impede iodine transport across the sodium-iodide symporter (NIS) membrane protein and thus prevent excessive iodine uptake by thyroid cells, which may lessen the occurrence of disease. Herein, we for the first time utilized nano-starch particles (St NPs) to regulate iodine transport across the NIS protein of thyroid cells by using extracellular interference therapy. By precisely encapsulating iodine within the cavity of a glucan α-helix via hydrogen bonding, extracellular St NPs prevented excess iodine uptake by thyroid cells in vitro and in vivo; this down-regulated the expression of NIS protein (0.06-fold) and autophagy protein LC3B-II (0.35-fold). We also found that St NPs regulated the metabolic pathway of iodine in zebrafish. We believe this proposed strategy offers a novel insight into controlling iodine uptake by the thyroid and indicates a new direction for preventing iodine-induced thyroid disease.

Revealing the mechanism underlying the effects of γ-aminobutyric acid-dioscorin interactions on dioscorin structure and emulsifying properties by molecular dynamic simulations

Many studies have shown that γ-aminobutyric acid (GABA) exhibits various beneficial biological activities, including gut-modulating, neuro-stimulating, and cardio-protecting activities. Naturally, GABA exists in small amounts in yam, which is primarily synthesized by the decarboxylation of L-glutamic acid in the presence of glutamate decarboxylase. Dioscorin, the major tuber storage protein of yam, has been shown to have good solubility and emulsifying activity. However, how GABA interacts with dioscorin and affects their properties has yet to be clarified. In this research, the physicochemical and emulsifying properties of GABA-fortified dioscorin, which was dried by spray drying and freeze drying, were studied. As results, the freeze-dried (FD) dioscorin produced more stable emulsions, while the spray-dried (SD) dioscorin adsorbed more rapidly to oil/water (O/W) interface. The fluorescence spectroscopy, ultraviolet spectroscopy and circular dichroism spectroscopy showed that GABA changed the structure of dioscorin, by exposing its hydrophobic groups. The addition of GABA significantly promoted the adsorption of dioscorin to the O/W interface and prevented droplets coalescence. The results of molecular dynamics simulation (MD) showed that GABA destroyed the H-bond network between dioscorin and water, increased surface hydrophobicity and finally improved the emulsifying properties of dioscorin.

More simple, efficient and accurate food research promoted by intermolecular interaction approaches: A review

The investigation of intermolecular interactions has become increasingly important in many studies, mainly by combining different analytical approaches to reveal the molecular mechanisms behind specific experimental phenomena. From spectroscopic analysis to sophisticated molecular simulation techniques like molecular docking, molecular dynamics (MD) simulation, and quantum chemical calculations (QCC), the mechanisms of intermolecular interactions are gradually being characterized more clearly and accurately, leading to revolutionary advances. This article aims to review the progression in the main techniques involving intermolecular interactions in food research and the corresponding experimental results. Finally, we discuss the significant impact that cutting-edge molecular simulation technologies may have on the future of conducting deeper exploration. Applications of molecular simulation technology may revolutionize the food research, making it possible to design new future foods with precise nutrition and desired properties.

Determining changes in crystallinity of rice starch after heat-moisture treatment using terahertz spectroscopy

To investigate the potential of Fourier-transform terahertz (FT-THz) spectroscopy to follow crystalline structure changes in rice starch after heat-moisture treatment (HMT), we measured the crystallinity by X-ray diffraction (XRD) spectra and found its correlation with THz spectra. According to A-type crystal structure and Vh-type crystalline structure of amylose-lipid complex (ALC) in rice starch, crystallinity is divided into A-type and Vh-type. The intensity of second derivative spectra peak at 9.0 THz was highly correlated with both A-type and Vh-type crystallinity. Additionally, other three peaks at 10.5 THz, 12.2 THz, and 13.1 THz were also sensitive to Vh-type crystalline structure. These results indicate that after HMT, the crystallinity of ALC (Vh-type) and A-type starch can be quantified using THz peaks.

The Mechanism Underlying the Amylose-Zein Complexation Process and the Stability of the Molecular Conformation of Amylose-Zein Complexes in Water Based on Molecular Dynamics Simulation

The aim of this study was to employ molecular dynamics simulations to elucidate the mechanism involved in amylose–zein complexation and the stability of the molecular conformation of amylose–zein complexes in water at the atomic and molecular levels. The average root mean square deviation and radius of gyration were lower for amylose–zein complexes (1.11 nm and 1 nm, respectively) than for amylose (2.13 nm and 1.2 nm, respectively), suggesting a significantly higher conformational stability for amylose–zein complexes than for amylose in water. The results of radial distribution function, solvent-accessible surface area, and intramolecular and intermolecular hydrogen bonds revealed that the amylose–zein interaction inhibited water permeation into the amylose cavity, leading to enhanced conformational stabilities of the V-type helical structure of amylose and the amylose–zein complexes. Furthermore, the amylose in amylose–zein complexes displayed the thermodynamically stable 4C1 conformation. These findings can provide theoretical guidance in terms of the application of protein on starch processing aiming to improve the physicochemical and functional properties of starch (such as swelling capacity, pasting properties, and digestibility) for developing novel low-digestibility starch–protein products.

Simultaneous improvement of the physical and biological properties of starch films by incorporating steviol glycoside-based solid dispersion

Bioactive compounds are frequently incorporated into polysaccharides (e.g., starch) to form active biodegradable films for food packaging, but some of them are water insoluble (e.g., curcumin, CUR) that will make the films have undesirable performance. Herein, CUR was successfully solubilized into the aqueous starch film solution by steviol glycoside (STE, a natural sweetener)-based solid dispersion. The mechanisms of solubilization and film formation were explored through molecular dynamic simulation and various characterization methods. The results showed that the amorphous state of CUR combined with micellar encapsulation of STE achieved the solubilization of CUR. STE and starch chains cooperated to form the film via hydrogen bonding, while CUR was uniformly and densely distributed within the film in the form of needle-like microcrystals. The as-prepared film exhibited high flexibility, great moisture barrier, and excellent UV barrier (UV transmittance: ∼0 %). Compared with the film containing CUR alone, the as-prepared film possessed higher release efficiency, antibacterial activity, and pH response sensitivity due to the assistance of STE. Hence, the introduction of STE-based solid dispersion can simultaneously improve the biological and physical properties of starch films, which provides a green, nontoxic, and facile strategy for the perfect integration of hydrophobic bioactive compounds and polysaccharide-based films.

Molecular mechanism of epicatechin gallate binding with carboxymethyl β-glucan and its effect on antibacterial activity

The non-covalent binding between flavanols and polysaccharides has impacts on their bioactivities, but the binding mechanism is less understood. This work aimed to unveil the non-covalent interactions between epicatechin gallate (ECG) and anionic carboxymethyl Poria cocos polysaccharide (CMP_N) at the structural and molecular level based on the synergistic antibacterial effect between them. The results suggested that there was hydrogen bonding, hydrophobic and electrostatic interaction between ECG and CMP_N, which was also supported by the results of molecular dynamics simulations. The resulting changes in physicochemical properties enhanced the antibacterial activity of the ECG-CMP_N mixture. More specifically, through two-dimensional Fourier transform infrared correlation spectrum (2D-FT-IR) and nuclear magnetic resonance spectroscopy (NMR) analysis, COO^- in CMP_N carboxymethyl and CO in ECG galloyl had the highest response priority and binding strength in the interaction, allowing us to conclude the critical functional groups that affect the non-covalent interactions of polysaccharide and flavanols and their bioactivities.

Pseudocatalytic Hydrogels with Intrinsic Antibacterial and Photothermal Activities for Local Treatment of Subcutaneous Abscesses and Breast Tumors

The intimate relationship between bacteria and tumors has triggered a lot of activities in the development and design of bioactive materials to concurrently respond to antitumor and antibacterial demands. Herein, a pseudocatalytic hydrogel (AM-I@Agar) with intrinsic antibacterial and photothermal activities, synthesized by incorporating prefabricated amylose-iodine nanoparticles into low-melting-point agarose hydrogel, is explored as a bioactive agent for local treatment of subcutaneous abscesses and breast tumors. The AM-I@Agar hydrogel depicts the ability of pseudocatalytic O₂ generation from H₂ O₂ to alleviate hypoxia. Meanwhile, the AM-I@Agar hydrogel exhibits temperature self-regulation features, beneficial for avoiding thermal injury during photothermal therapy owing to thermochromic properties. Upon local injection into a subcutaneous abscess, methicillin-resistant Staphylococcus aureus is effectively eliminated by the AM-I@Agar hydrogel, and complete skin recovery is achieved in 8 d, demonstrating much better antibacterial effects compared with penicillin, a small-molecule antibiotic. AM-I/5-FU@Agar hydrogel, obtained after loading 5-fluorouracil (5-FU), significantly inhibits tumors in both normal 4T1 tumor-bearing mice and MRSA-infected 4T1 tumor-bearing mice models via a synergistic photothermal-chemo effect, and shows treatment efficiency superior to that achieved with photothermal therapy or 5-FU alone. This work provides a concept for the design and development of bioactive agents for potential management of bacteria-associated cancer.

A review on application of molecular simulation technology in food molecules interaction

Molecular simulation is a new technology to analyze the interaction between molecules. This review mainly summarizes the application of molecular simulation technology in the food industry. This technology has been employed to assess structural changes of biomolecules, the interaction between components, and the mechanism of physical and chemical property alterations. These conclusions provide a deeper understanding of the molecular interaction mechanism in foods, break through the limitations of scientific experiments and avoid blind and time-consuming scientific research. In this paper, the advantages and development trends of molecular simulation technology in the food research field are described. This methodology can be used to contribute to further studies of the mechanism of molecular interactions in food, confirm experimental results and provide new ideas for research in the field of food sciences.

Influence of different kinds of fatty acids on the behavior, structure and digestibility of high amylose maize starch-fatty acid complexes

BACKGROUND

The formation of starch-lipid complexes is of interest to food processing and human nutrition. Fatty acid (FA) structure is important for the formation and structure of starch-FA complexes. However, there is limited research regarding the complexing behavior between amylose and different kinds of FAs, as well as the relationship between fine structures and digestibility of the formed complexes. This study aimed to investigate the behavior, fine structure, and digestibility of complexes formed between high amylose maize starch (HMS) and FA having various chain lengths and unsaturation degrees.

RESULTS

Complexes containing different FA structures showed V_6III -type crystals. Complexes containing 18-carbon unsaturated FAs displayed significantly higher complexing index (P < 0.05) than other complexes. Complexes containing 12-carbon FAs and 18-carbon FAs with one unsaturation degree showed a higher degree of structural order and resistant starch (RS) content than other complexes. The 12-carbon FAs exhibited a higher binding degree with helical cavity of amylose than other FAs. Additionally, 10-carbon and 18-carbon saturated FAs tended to combine with HMS outside amylose helices more than other FAs. Laser confocal micro-Raman imaging revealed that the physically embedded 10-carbon and 18-carbon saturated FAs showed heterogeneous distribution in complexes, and that the complexed 18-carbon FAs with one unsaturation degree exhibited homogeneous distribution.

CONCLUSION

The behavior, structural order and digestibility of complexes could be regulated by FA structure. The 12-carbon FAs and 18-carbon FAs with one unsaturation degree were more suitable for the production of HMS-FA complexes with higher structural order and RS content than other FAs. © 2022 Society of Chemical Industry.

Formation and molecular dynamics simulation of inclusion complex of large-ring cyclodextrin and 4-terpineol

In the present study, the formation and structure of the inclusion compound of large-ring cyclodextrin and 4-terpineol were obtained through different experiments and molecular dynamics (MD) simulation. The analysis of FTIR, ¹ H-NMR, and thermodynamic results confirmed the formation of clathrates. Analysis of molecular structure (root-mean-square deviation and radius of gyration), solubility, and interaction energy (Coul, H bond) based on MD simulations further clarified the nature of the clathrate and the conformational changes caused by guest molecules as well as inclusion complexes process trends. The inclusion complex reportedly has a new crystal structure with improved thermal stability. PRACTICAL APPLICATION: This is the first work to demonstrate the complex formation between 4-terpineol and large-ring cyclodextrin by molecular dynamics simulation. Molecular dynamics simulation confirmed the formation of inclusion complexes theoretically. Conformational changes of the molecules and the formation of complexes with improved thermal stability were observed. Complexing with large-ring cyclodextrin can be used as an effective means to encapsulate the aroma/flavor compounds.

Mechanisms Underlying the Formation of Amylose- Lauric Acid-β-Lactoglobulin Complexes: Experimental and Molecular Dynamics Studies

The aim of the present study was to reveal the mechanisms underlying the formation of ternary complexes with a model system of amylose (AM), lauric acid (LA), and β-lactoglobulin (βLG) using experimental studies and molecular dynamics (MD) simulations. Experimental analyses showed that hydrophobic interactions and hydrogen bonds contributed more than electrostatic forces to the formation of the AM-LA-βLG complex. MD simulations indicated that interactions between AM and βLG through electrostatic forces and hydrogen bonds, and to a less extent van der Waals forces, and interactions between AM and LA through van der Waals forces, were mostly responsible for complex formation. The combination of experimental results and MD simulations has provided new mechanistic insights and led us to conclude that hydrophobic interactions, van der Waals forces between AM and LA, and van der Waals forces and hydrogen bonds between AM and βLG were the main driving forces for the formation of the AM-LA-βLG complex.

Exploration of the process and mechanism of magnesium chloride induced starch gelatinization

As a new starch gelatinization method, salt induced gelatinization can not only reduce energy consumption but also impart special physicochemical properties to starch gel. In this study, the process and mechanism of MgCl₂ induced starch gelatinization were explored. The results showed that, potato starch could be gelatinized after a treatment of 4 mol/L MgCl₂ for 3 h. The gelatinization started with the slight damage of outer shells, then the internal molecules leached out through the cracks or holes to form gel, finally the outer shells disintegrated. During the gelatinization process, the viscosity and granule size gradually increased after 0.5 h, while the original crystallinity disappeared rapidly in 0.5 h. Besides, MgCl₂ significantly increased the electrostatic interaction, then made starch molecules closer to each other and become denser, which may have close relationship with the appearance of the cracks and the disappearance of crystallization. Moreover, MgCl₂ enhanced the hydration and increased the binding free energy of starch molecules, then promoted starch gelatinization and accelerated the destruction of starch structure, which may be the critical factors of the starch gelatinization induced by MgCl₂. The results will provide reference for the research and application of salt induced gelatinization.

KOH/thiourea aqueous solution: A potential solvent for studying the dissolution mechanism and chain conformation of corn starch

Non-derivatizing, high-efficiency and low-toxicity solvents are important for studying the dissolution behavior and potential applications of starch. In this study, we investigated the starch dissolution mechanism and molecular conformation in KOH/thiourea aqueous solutions and compared these with KOH/urea and KOH aqueous solutions. Solubility analysis revealed that the KOH/thiourea solution demonstrates a better ability to dissolve corn starch than KOH/urea and KOH solutions. Rheological behavior and dynamic and static light scattering indicated that starch is stable in KOH/thiourea solution and exists as a regular star structure. Fourier transform infrared spectroscopy, ¹³C NMR, and molecular dynamics simulations indicated that hydrated K⁺ and OH^- destroy the strong starch hydrogen bond interactions; thiourea hydrate self-assembles into a shell surrounding the starch-KOH complex through interaction with KOH, whereas there is no direct strong interaction between urea and KOH. Therefore, adding thiourea to a KOH solution can promote dissolution and prevent self-aggregation of the starch chain.

Evidence for amylose inclusion complexes with multiple acyl chain lipids using solid-state NMR and theoretical approaches

Amylose is known to form inclusion complexes in the presence of hydrophobic guests. Among lipids, only single-chain fatty acids have been reported as possible guests with the surrounding amylose in a well-defined V-helix conformation. Using experimental ¹³C solid-state NMR, we studied the formation of inclusion complexes between amylose and a variety of multiple-chains lipids of increasing complexity. Molecular dynamics simulations and calculations of ¹³C isotropic chemical shifts using the Density Functional Theory approach were performed to support the interpretation of experimental spectra. We provide unambiguous evidences that amylose forms inclusion complexes with lipids bearing multiple acyl chains. Amylose conformations around these lipids are characterized by {ϕ,ψ} anomeric bond dihedral angles near {115°,105°}. In the ¹³C NMR spectra, this translates into C₁ and C₄ chemical shifts of 102.5 ppm and 81.1 ppm, regardless of the helical conformation of the amylose surrounding the acyl chains.

Engineered Stable 5-Hydroxymethylfurfural Oxidase (HMFO) from 8BxHMFO Variant of Methylovorus sp. MP688 through B-Factor Analysis

What is known as Furan-2,5-dicarboxylic acid (FDCA) is an attractive compound since it has similar properties to terephthalic acid. Further, 5-hydroxymethylfurfural oxidase (HMFO) is an enzyme, which could convert HMF to FDCA directly. Most wild types of HMFO have low activity on the oxidation of HMF to FDCA. The variant of 8BxHFMO from Methylovorus sp. MP688 was the only reported enzyme that was able to perform FDCA production. However, the stabilization of 8BxHMFO is still not that satisfactory, and further improvement is necessary for the industrial application of the enzyme. In this work, stability-enhanced HMFO from 8BxHFMO was engineered through employing B-factor analysis. The mutation libraries were created based on the NNK degeneracy of residues with the top ten highest B-factor value, and two of the effective mutants were screened out through the high throughput selection with the horseradish peroxidase (HRP)-Tyr assay. The mutants Q319K and N44G show a significantly increased yield of FDCA in the reaction temperature range of 30 to 40 °C. The mutant Q319K shows the best performance at 35 °C with a FDCA yield of 98% (the original 8BxHMFO was only 85%), and a half-life exceeding 72 h. Moreover, molecular dynamic simulation indicates that more hydrogen bonds are formed in the mutants, which improves the stability of the protein structure. The method could enhance the design of more stable biocatalysts; and provides potential for the further optimization and utilization of HMFO in biotechnological processes.

Complexation of 26-Mer Amylose with Egg Yolk Lipids with Different Numbers of Tails Using a Molecular Dynamics Simulation

A molecular dynamics simulation of mixtures of 26-mer amylose with three different egg yolk lipids, namely, cholesterol, triglyceride and 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC), demonstrated the formation of a stable complex. The 26-mer amylose fluctuated between a coiled and an extended helical conformation. The complex was a V-type amylose complex, with the hydrophobic tail of the lipids being inside the hydrophobic helical cavity of the amylose. The number of glucose units per turn was six for the two helical regions of the amylose-POPC complex and the palmitoyl tail region of the amylose-triglyceride complex. This value was eight for the cholesterol and the two-tail helical region in the amylose-triglyceride complex. Two tails of the POPC were in two different hydrophobic helical regions of the 26-mer amylose, whereas the palmitoyl tail of the triglyceride lay in one hydrophobic helical region and the linoleoyl and oleoyl tails both lay in another helical region, and the cross-sectional area of the latter was larger than the former to accommodate the two tails. The radii of the gyration of the complex were lower for all three cases compared to that of one single amylose. In addition, the stability of the complexes was ranked in the following order: POPC < cholesterol < triglyceride, with their average binding energy being −97.83, −134.09, and −198.35 kJ/mol, respectively.

Molecular simulation, characteristics and mechanism of thermal-responsive acetylated amylose V-type helical complexes

To explore the thermal-responsive characteristics of acetylated amylose-guest V-type helical complexes (AAGHCs) and their potential use as thermal-responsive drug carriers, different types of AAGHCs were built, in which acetylated amylose was used as a host, and iodine, propofol, or hexane was utilized as the guest molecule. Their thermal-responsive characteristics were investigated through molecular dynamic (MD) simulation and corresponding experiments. MD simulation showed that the thermal-responsive helix-unfolding and guest-release behavior in AAGHCs, and the complete unfolding of AAGHC could be divided into brewing, triggering and collapsing periods. Energy analysis revealed that the Lana-Jones potential is an important binding energy that bridges host and guest molecules and enhances the stability of the helix. The various types or number of guests showed different binding energies. The stronger the binding energy, higher is the temperature required to trigger the unfolding of the helix and the releasing of guests. FT-IR and X-ray diffraction analyses confirmed the structures of AAGHCs. The change in hydrated size, and UV-VIS absorption of AAGHCs at high temperatures both confirmed the thermal-responsiveness of AAGHCs. The fluorescence fluctuation of loaded 7-hydroxycoumarin reflected the same thermal-responsive process and mechanism as MD simulation. This study provides meaningful theoretical guidance for the design of thermal-responsive drug carriers based on acetylated amylose-guest V-type helical complexes.

Using geometric criteria to study helix-like structures produced in molecular dynamics simulations of single amylose chains in water

Amylose is a linear polymer chain of α-d-glucose units connected through α(1 → 4) glycosidic bonds. Experimental studies show that in non-polar solvents, single amylose chains form helical structures containing precise H-bond patterns. However, both experimental and computational studies indicate that these perfectly H-bonded helices are not stable in pure water. Nevertheless, amylose chains are observed to form helix-like structures in molecular dynamics (MD) simulations that exhibit imperfect H-bond patterns. In this paper, we study the structure of amylose chains in water using MD simulations to identify and characterize these “imperfect” helical structures. To this end we devise geometry-based criteria to define imperfect helical structures in amylose chains. Using this approach, the propensity of amylose chains to form these structures is quantified as a function of chain length and solvent temperature. This analysis also uncovers both short and long time helix-breaking mechanisms such as band-flips and kinks in the chain. This geometric approach to defining imperfect helices thus allows us to give new insight into the secondary structure of single amylose chains in spite of imperfect H-bond patterns.

We introduce a geometrical approach to capture and study helix-like structures in MD simulations of single amylose chains in water.

Multi-scale stabilization mechanism of pickering emulsion gels based on dihydromyricetin/high-amylose corn starch composite particles

For Pickering emulsifying effect, starch must be subjected to the pretreatments of acid hydrolysis, esterification, which are complicated and eco-unfriendly. In this study, a practical and green strategyto fabricate Pickering emulsion gels with dihydromyricetin (DMY)/high-amylose corn starch (HCS) composite particles was introduced for the first time. The DMY content in composite particles and the amount of addition of composite particles had obvious synergistic effect on the formation and properties of emulsion gels. The obtained emulsion gels were not sensitive to ionic strength, which could be attributed to emulsifying capacity and viscosity effect of composite particles. The spectral analysis confirmed the presence of DMY/amylose host-guest supramolecules. The molecular simulation of the supramolecular complexes in the oil-water system indicated that these complexes could spontaneously aggregate and anchor to the oil-water interface, reducing the interfacial tension. Based on experimental and theoretical results, the multi-scale relationship of "molecular interaction-particle characteristics-gel properties" was established.

Physicochemical properties of heat-moisture treated, stearic acid complexed starch: The effect of complexation time and temperature

Starch modification has been extensively studied to alter its physicochemical properties based on human needs. Lowering the digestion rate of starch is one of the interests in food science research, since when it is nutritionally improved, it can reduce the risk of human chronic diseases. In this study, heat-moisture treatment (HMT) followed by inclusion complexation with stearic acid at various temperatures and times was applied to improve the functional properties of starch. Thermal analysis suggested the formation of type I and type II complexes after complexation at 90 °C, indicated by a endothermal peak at 107 and 122 °C, respectively, while native starch after complexation only resulted in type I complexes. The formation of crystalline complexes was also confirmed by XRD showing peaks at 2θ = 13.1° and 20.1°. Furthermore, the modified starch displayed a higher pasting temperature, considerably less swelling and significantly lower viscosity behavior. This implied that the starch granules were thermally and mechanically more stable. The granular appearance of the modified starch was confirmed with light microscopy that presented more intact granules and less ruptured granules, even after heating to 90 °C. This study offers a way to upgrade the nutritional properties of starch.

Characterizing moisture uptake and plasticization effects of water on amorphous amylose starch models using molecular dynamics methods

Dynamics and thermophysical properties of amorphous starch were explored using molecular dynamics (MD) simulations. Using the OPLS3e force field, simulations of short amylose chains in water were performed to determine force field accuracy. Using well-tempered metadynamics, a free energy map of the two glycosidic angles of an amylose molecule was constructed and compared with other modern force fields. Good agreement of torsional sampling for both solvated and amorphous amylose starch models was observed. Using combined grand canonical Monte Carlo (GCMC)/MD simulations, a moisture sorption isotherm curve is predicted along with temperature dependence. Concentration-dependent activation energies for water transport agree quantitatively with previous experiments. Finally, the plasticization effect of moisture content on amorphous starch was investigated. Predicted glass transition temperature (T_g) depression as a function of moisture content is in line with experimental trends. Further, our calculations provide a value for the dry T_g for amorphous starch, a value which no experimental value is available.

Encapsulation of caffeine into starch matrices: Bitterness evaluation and suppression mechanism

In this study, caffeine (CA) was encapsulated into food-grade starch matrices, including swelled starch (SS), porous starch (PS), and V-type starch (VS). The bitterness of the microcapsules and suppression mechanisms were investigated using an electronic tongue, molecular dynamics (MD) simulation and the in vitro release kinetics of CA. All the CA-loaded microcapsules showed a lower bitterness intensity than the control. The MD results proved that the weak interactions between starch and CA resulted in a moderate CA release rate for SS-CA microcapsules. The PS-CA microcapsule presented the longest CA release, up to 40 min, whereas the VS-CA microcapsule completely released CA in 9 min. The CA release rate was found to be related to the microcapsule structure and rehydration properties. A low CA bitterness intensity could be attributed to a delay in the CA release rate and resistance to erosion of the microcapsules. The results of this work are valuable for improving starch-based microcapsules (oral-targeted drug-delivery systems) by suppressing the bitterness of alkaloid compounds.

Diversity of the Hydrogen Bond Network and Its Impact on NMR Parameters of Amylose B Polymorph: A Study Using Molecular Dynamics and DFT Calculations Within Periodic Boundary Conditions

Classical molecular dynamics simulations have been combined with quantum (DFT) calculations of ¹³C NMR parameters in order to relate the experimental spectrum of the double-helix form of the amylose B-polymorph in highly crystalline conditions not only to its 3D structure but also to the arrangement of atoms in the crystal lattice. Structures obtained from these simulations or from geometry optimization procedures at the DFT level have shown the presence of hydrogen bond networks between sugars of the same helix or between residues of the two chains of the double helix. ¹³C NMR parameter calculations have revealed the impact of such a network on the chemical shifts of carbon atoms. In addition, DFT calculations using periodic boundary conditions were compulsory to highlight the presence of two types of sugar within the crystal sample. It allows us to confirm, theoretically, the experimental hypothesis that the existence of two distinct sugar types in the NMR spectrum is a consequence of crystal packing.