Application of molecular dynamic simulation to study food proteins: A review

Exploring high hydrostatic pressure effects on anthocyanin binding to serum albumin and food-derived transferrins

This study investigated the effects of high hydrostatic pressure (HHP) on the binding interactions of cyanindin-3-O-glucoside (C3G) to bovine serum albumin, human serum albumin (HSA), bovine lactoferrin, and ovotransferrin. Fluorescence quenching revealed that HHP reduced C3G-binding affinity to HSA, while having a largely unaffected role for the other proteins. Notably, pretreating HSA at 500 MPa significantly increased its dissociation constant with C3G from 24.7 to 34.3 μM. Spectroscopic techniques suggested that HSA underwent relatively pronounced tertiary structural alterations after HHP treatments. The C3G-HSA binding mechanisms under pressure were further analyzed through molecular dynamics simulation. The localized structural changes in HSA under pressure might weaken its interaction with C3G, particularly polar interactions such as hydrogen bonds and electrostatic forces, consequently leading to a decreased binding affinity. Overall, the importance of pressure-induced structural alterations in proteins influencing their binding with anthocyanins was highlighted, contributing to optimizing HHP processing for anthocyanin-based products.

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.

Integrating serum pharmacochemistry, network pharmacology and untargeted metabolomics strategies to reveal the material basis and mechanism of action of Feining keli in the treatment of chronic bronchitis

ETHNOPHARMACOLOGICAL RELEVANCE

Feining keli (FNKL) is herbal preparation mainly made from Senecio cannabifolius Less., In recent years, more and more studies have found that FNKL has excellent therapeutic effects on chronic bronchitis (CB). Nevertheless, its pharmacodynamic material basis and mechanism of action are still unknown.

AIM OF THE STUDY

This study aimed to explore the pharmacodynamic material basis and mechanism of action of FNKL in treating CB.

MATERIALS AND METHODS

The CB rat model was induced using nasal drops of lipopolysaccharide (LPS) in combination with smoking. Various assessments including behavioral and body mass examination, lung index measurement, enzyme linked immunosorbent assay (ELISA), as well as histological analyses using hematoxylin and eosin (H&E) and Masson staining were conducted to validate the reliability of the CB model. The serum components of FNKL in CB rats were identified using ultra-high-performance liquid chromatography Orbitrap Exploris mass spectrometer (UHPLC-OE-MS). Network pharmacology was used to predict the network of action of the active ingredients in FNKL based on these serum components. Signaling pathways were enriched and analyzed, and molecular docking was conducted for key targets. Molecular dynamics simulations were performed using GROMACS software. The mechanism was confirmed through a series of experiments including Western blot (WB), immunofluorescence (IF), and reverse transcription (RT)-PCR. Additionally, untargeted metabolomics was employed to identify biomarkers and relevant metabolic pathways associated with the treatment of CB with FNKL.

RESULTS

In CB rats, FNKL improved body mass, lung index, and pathological damage of lung tissues. It also decreased interleukin (IL)-6, tumor necrosis factor-alpha (TNF-α), malonaldehyde (MDA) levels, and percentage of lung collagen fiber area. Furthermore, FNKL increased IL-10 and superoxide dismutase (SOD) levels, which helped alleviate bronchial inflammation in the lungs. A total of 70 FNKL chemical components were identified in CB rat serum. Through network pharmacology analysis, 5 targets, such as PI3K, AKT, NF-κB, HIF-1α, and MYD88, were identified as key targets of FNKL in the treatment of CB. Additionally, the key signaling pathways identified were PI3K/AKT pathway、NF-κB/MyD88 pathway、HIF-1α pathway. WB, IF, and RT-PCR experiments were conducted to confirm the findings. Molecular docking studies demonstrated successful docking of 16 potential active components with 5 key targets. Additionally, molecular dynamics simulations indicated the stability of quercetin-3-galactoside and HIF-1α. Metabolomics analysis revealed that FNKL primarily regulated pathways related to alpha-linolenic acid metabolism, primary bile acid biosynthesis, bile secretion, arachidonic acid metabolism, neuroactive ligand-receptor interaction, and folate biosynthesis. Furthermore, the expression levels of traumatic acid, traumatin, alpha linolenic acid, cholic acid, 2-arachidonoylglycerol, deoxycholic acid, 7,8-dihydroneopterin, and other metabolites were found to be regulated.

CONCLUSION

FNKL exhibits positive therapeutic effects on CB, with quercetin-3-galactoside identified as a key active component. The mechanism of FNKL's therapeutic action on CB involves reducing inflammatory response, oxidative stress, and regulating metabolism, and its molecular mechanism was better elucidated in a holistic manner. This study serves as a reference for understanding the pharmacodynamic material basis and mechanism of action of FNKL in treating CB, and provides avenues for exploring the effects of compounded herbal medicines on CB.

Computational Chemistry Strategies to Investigate the Antioxidant Activity of Flavonoids-An Overview

Despite several decades of research, the beneficial effect of flavonoids on health is still enigmatic. Here, we focus on the antioxidant effect of flavonoids, which is elementary to their biological activity. A relatively new strategy for obtaining a more accurate understanding of this effect is to leverage computational chemistry. This review systematically presents various computational chemistry indicators employed over the past five years to investigate the antioxidant activity of flavonoids. We categorize these strategies into five aspects: electronic structure analysis, thermodynamic analysis, kinetic analysis, interaction analysis, and bioavailability analysis. The principles, characteristics, and limitations of these methods are discussed, along with current trends.

Improving the Quality of Wheat Flour Bread by a Thermophilic Xylanase with Ultra Activity and Stability Reconstructed by Ancestral Sequence and Computational-Aided Analysis

Xylanase is an essential component used to hydrolyze the xylan in wheat flour to enhance the quality of bread. Presently, cold-activated xylanase is popularly utilized to aid in the development of dough. In this study, ancestral sequence reconstruction and molecular docking of xylanase and wheat xylan were used to enhance the activity and stability of a thermophilic xylanase. The results indicated that the ancestral enzyme TmxN3 exhibited significantly improved activity and thermal stability. The Vmax increased by 2.7 times, and the catalytic efficiency (Kcat/Km) increased by 1.7 times in comparison to TmxB. After being incubated at 100 °C for 120 min, it still retained 87.3% of its activity, and the half-life in 100 °C was 330 min, while the wild type xylanase was only 55 min. This resulted in an improved shelf life of bread, while adding TmxN3 considerably enhanced its quality with excellent volume and reduced hardness, chewiness, and gumminess. The results showed that the hardness was reduced by 55.2%, the chewiness was reduced by 40.11%, and the gumminess was reduced by 53.52%. To facilitate its industrial application, we further optimized the production conditions in a 5L bioreactor, and the xylanase activity reached 1.52 × 106 U/mL culture.

Interaction of the extracellular protease from Staphylococcus xylosus with meat proteins elucidated via spectroscopic and molecular docking

This study aimed to assess the effect of an external protease secreted by Staphylococcus (S.) xylosus on the hydrolysis and flavor properties of meat protein. The results indicated that the protease significantly increased the solubility of myofibrillar proteins (MPs) and sarcoplasmic proteins (SPs) in water (P < 0.05), and altered their surface hydrophobicity and secondary structure. The results of micromorphological and free amino acids analyses suggested that the protease degraded the large and insoluble meat protein aggregates into small molecular proteins with uniform distribution and amino acids, especially glycine, glutamic acid, leucine, and cysteine. Moreover, the protease-catalyzed hydrolysis promoted the formation of some volatile compounds in the MPs and SPs. Additionally, molecular docking analysis suggested that hydrogen bond and hydrophobic interaction promoted the formation of a S. xylosus protease/meat protein complex. These results provided a basis for the future application of S. xylosus protease in meat products.

Pinobanksin from peony seed husk: A flavonoid with the potential to inhibit the proliferation of SH-SY5Y

Pinobanksin, as one of the flavonoids, has powerful biological activities but has been under-recognized. In this study, we optimized the extraction method of phragmites from peony seed shells by using organic solvent extraction. The yield of PSMS was 10.54 ± 0.13% under the conditions of ethanol volume fraction 70%, extraction temperature 70°C, material-liquid ratio 1:25 g/mL, and extraction time 60 min; the optimized PSMS could be effectively separated in S-8 macroporous resin coupled with C18. The relative content of PSMS was increased from 0.42% in PSMS to 92.53% after C18 purification; the antioxidant activity test revealed that pinobanksin could exert antioxidant ability by binding catalase (CAT) enzyme. Second, it was found that pinobanksin could effectively inhibit the proliferation of SH-SY5Y cells, mainly by binding to BCL2-associated X (BAX), B-cell lymphoma-2 (BCL-2), and cyclin-dependent Kinase 4/6 (CDK4/6) to produce more hydrogen bonds to inhibit their activities. This study confirms the medicinal potential of pinobanksin and provides the basis for the proper understanding of pinobanksin and the development of related products.

Integrating network pharmacology and animal experimental validation to investigate the action mechanism of oleanolic acid in obesity

BACKGROUND

Obesity, a condition associated with the development of widespread cardiovascular disease, metabolic disorders, and other health complications, has emerged as a significant global health issue. Oleanolic acid (OA), a pentacyclic triterpenoid compound that is widely distributed in various natural plants, has demonstrated potential anti-inflammatory and anti-atherosclerotic properties. However, the mechanism by which OA fights obesity has not been well studied.

METHOD

Network pharmacology was utilized to search for potential targets and pathways of OA against obesity. Molecular docking and molecular dynamics simulations were utilized to validate the interaction of OA with core targets, and an animal model of obesity induced by high-fat eating was then employed to confirm the most central of these targets.

RESULTS

The network pharmacology study thoroughly examined 42 important OA targets for the treatment of obesity. The key biological processes (BP), cellular components (CC), and molecular functions (MF) of OA for anti-obesity were identified using GO enrichment analysis, including intracellular receptor signaling, intracellular steroid hormone receptor signaling, chromatin, nucleoplasm, receptor complex, endoplasmic reticulum membrane, and RNA polymerase II transcription Factor Activity. The KEGG/DAVID database enrichment study found that metabolic pathways, PPAR signaling pathways, cancer pathways/PPAR signaling pathways, insulin resistance, and ovarian steroidogenesis all play essential roles in the treatment of obesity and OA. The protein-protein interaction (PPI) network was used to screen nine main targets: PPARG, PPARA, MAPK3, NR3C1, PTGS2, CYP19A1, CNR1, HSD11B1, and AGTR1. Using molecular docking technology, the possible binding mechanism and degree of binding between OA and each important target were validated, demonstrating that OA has a good binding potential with each target. The molecular dynamics simulation's Root Mean Square Deviation (RMSD), and Radius of Gyration (Rg) further demonstrated that OA has strong binding stability with each target. Additional animal studies confirmed the significance of the core target PPARG and the core pathway PPAR signaling pathway in OA anti-obesity.

CONCLUSION

Overall, our study utilized a multifaceted approach to investigate the value and mechanisms of OA in treating obesity, thereby providing a novel foundation for the identification and development of natural drug treatments.

Research advances of molecular docking and molecular dynamic simulation in recognizing interaction between muscle proteins and exogenous additives

During storage and processing, muscle proteins, e.g. myosin and myoglobin, will inevitably undergo degeneration, which is thus accompanied by quality deterioration of muscle foods. Some exogenous additives have been widely used to interact with muscle proteins to stabilize the quality of muscle foods. Molecular docking and molecular dynamics simulation (MDS) are regarded as promising tools for recognizing dynamic molecular information at atomic level. Molecular docking and MDS can explore chemical bonds, specific binding sites, spatial structure changes, and binding energy between additives and muscle proteins. Development and workflow of molecular docking and MDS are systematically summarized in this review. Roles of molecular simulations are, for the first time, comprehensively discussed in recognizing the interaction details between muscle proteins and exogenous additives aimed for stabilizing color, texture, flavor, and other properties of muscle foods. Finally, research directions of molecular docking and MDS for improving the qualities of muscle foods are discussed.

Computational study to investigate Proteus mirabilis proteomes for multi-epitope vaccine construct design

Proteus mirabilis is a gram-negative bacterium particularly known for its unique swarming ability. The swarming gives the bacteria ability to enhance adherence to the catheter surface and epithelium cells of the urethra to cause catheter associated urinary tract infections. P. mirabilis has evolved resistant to antibiotics. Additionally, there is an approved vaccine against P. mirabilis, thus demanding for identification of new vaccine targets. This gram-negative bacterium consists of 19,502 core proteins, out of which 19,063 are redundant proteins and remaining 439 are non-redundant proteins. The non-redundant proteins have 21 proteins present on the cell surface out of which 11 proteins are virulent. Antigenicity analysis predicted only 2 proteins as antigenic (fimbrial biogenesis outer membrane usher protein and ligand-gated channel protein). Four and seven B-cells epitopes were predicted from the former and later proteins, respectively. The predicted B-cells epitopes were used for T- cells epitopes prediction. The predicted epitopes were linked to each other through GPGPG linkers and joined with cholera toxin beta subunit adjuvant. A multi-epitopes vaccine construct consisting of 226 residues was docked with MHC-I, MHC-II and TLR-4. The best docked complex in each case has binding energy of -714.6, -744.6 and -829.5 kcal/mol, respectively. Moreover, the docking results were validated through molecular dynamics simulation and binding free energies estimation. The net energy of -137.2 kcal/mol was calculated for vaccine-MHC-I complex, -133.39 kcal/mol for vaccine-MHC-II and -158.68 kcal/mol for vaccine-TLR-4 complex. The designed vaccine construct could provoke immune responses against targeted pathogen and may be used in experimental testing.Communicated by Ramaswamy H. Sarma.

Computational insight into stability-enhanced systems of anthocyanin with protein/peptide

Anthocyanins, which belong to the flavonoid group, are commonly found in the organs of plants native to South and Central America. However, these pigments are unstable under conditions of varying pH, heat, etc., which limits their potential applications. One method for preserving the stability of anthocyanins is through encapsulation using proteins or peptides. Nevertheless, the complex and diverse structure of these molecules, as well as the limitation of experimental technologies, have hindered a comprehensive understanding of the encapsulation processes and the mechanisms by which stability is enhanced. To address these challenges, computational methods, such as molecular docking and molecular dynamics simulation have been used to study the binding affinity and dynamics of interactions between proteins/peptides and anthocyanins. This review summarizes the mechanisms of interaction between these systems, based on computational approaches, and highlights the role of proteins and peptides in the stability enhancement of anthocyanins. It also discusses the current limitations of these methods and suggests possible solutions.

A cheminformatics-biophysics correlate to identify promising lead molecules against matrix metalloproteinase-2 (MMP-2) enzyme: A promising anti-cancer target

Matrix metalloproteinase-2 (MMP-2) is an endopeptidase enzyme that is devoted to extracellular matrix proteins degradation. The enzyme is warranted as promising drugs target for different light threating diseases such as arthritis, cancer and fibrosis. Herein, in this study, three drug molecules: CMNPD8322, CMNPD8320, and CMNPD8318 were filtered as high affinity binding compounds with binding energy score of -9.75 kcal/mol, -9.11 kcal/mol, -9.05 kcal/mol, respectively. The control binding energy score was -9.01 kcal/mol. The compounds docked deeply inside the pocket interacting with S1 pocket residues. The docked complexes dynamics in real time at cellular environment was then done to decipher the stable binding conformation and intermolecular interactions network. The compounds complexes achieved very stable dynamics with root mean square deviation (RMSD) with mean value of around 2-3 Å compared to control complex that showed higher fluctuations of 5 Å. The simulation trajectories frames based binding free energy demonstrated all the compounds-MMP-2 complexes reported highly stable energy, particularly the van der Waals energy dominate the overall net energy. Similarly, the complexes revalidation of WaterSwap based energies also disclosed the complexes highly stable in term docked conformation. Also, the compounds illustrated the compounds favorable pharmacokinetics and were non-toxic and non-mutagenic. Thus, the compounds might be used thorough experimental assays to confirm compounds selective biological potency against MMP-2 enzyme.

Application of molecular dynamics simulation for exploring the roles of plant biomolecules in promoting environmental health

Understanding the dynamic changes of plant biomolecules is vital for exploring their mechanisms in the environment. Molecular dynamics (MD) simulation has been widely used to study structural evolution and corresponding properties of plant biomolecules at the microscopic scale. Here, this review (i) outlines structural properties of plant biomolecules, and the crucial role of MD simulation in advancing studies of the biomolecules; (ii) describes the development of MD simulation in plant biomolecules, determinants of simulation, and analysis parameters; (iii) introduces the applications of MD simulation in plant biomolecules, including the response of the biomolecules to multiple stresses, their roles in corrosive environments, and their contributions in improving environmental health; (iv) reviews techniques integrated with MD simulation, such as molecular biology, quantum mechanics, molecular docking, and machine learning modeling, which bridge gaps in MD simulation. Finally, we make suggestions on determination of force field types, investigation of plant biomolecule mechanisms, and use of MD simulation in combination with other techniques. This review provides comprehensive summaries of the mechanisms of plant biomolecules in the environment revealed by MD simulation and validates it as an applicable tool for bridging gaps between macroscopic and microscopic behavior, providing insights into the wide application of MD simulation in plant biomolecules.

Allergenicity of peanut allergens and its dependence on the structure

Food allergies are a global food safety problem. Peanut allergies are common due, in part, to their popular utilization in the food industry. Peanut allergy is typically an immunoglobulin E-mediated reaction, and peanuts contain 17 allergens belonging to different families in peanut. In this review, we first introduce the mechanisms and management of peanut allergy, followed by the basic structures of associated allergens. Subsequently, we summarize methods of epitope localization for peanut allergens. These methods can be instrumental in speeding up the discovery of allergenicity-dependent structures. Many attempts have been made to decrease the allergenicity of peanuts. The structures of hypoallergens, which are manufactured during processing, were analyzed to strengthen the desensitization process and allergen immunotherapy. The identification of conformational epitopes is the bottleneck in both peanut and food allergies. Further, the identification and modification of such epitopes will lead to improved strategies for managing and preventing peanut allergy. Combining traditional wet chemistry research with structure simulation studies will help in the epitopes' localization.

Usnic acid as potential inhibitors of BCL2 and P13K protein through network pharmacology-based analysis, molecular docking and molecular dynamic simulation

Usnic acid (UA) lately piqued the interest of researchers for its extraordinary biological characteristics, including anticancer activity. Here, the mechanism was clarified through network pharmacology,molecular docking and molecular dynamic simulation. Sixteen proteins were selected through network pharmacology study as they are probable to interact with UA. Out of these proteins, 13 were filtered from PPI network analysis based on their significance of interactions (p < 0.05). KEGG pathway analysis has also aided us in determining the three most significant protein targets for UA, which are BCL2, PI3KCA and PI3KCG. Therefore molecular docking and molecular dynamic (MD) simulations throughout 100 ns were performed for usnic acid onto the three proteins mentioned. However, UA's docking score in all proteins is lower than their co-crystalised ligand, especially for BCL2 (-36.5158 kcal/mol) and PI3KCA (-44.5995 kcal/mol) proteins. The only exception is PI3KCG which has comparable results with the co-crystallised ligand with (-41.9351 kcal/mol). Furthermore, MD simulation has also revealed that usnic acid does not stay fit in the protein throughout the simulation trajectory for PI3KCA protein evident from RMSF and RMSD plots. Nevertheless, it still poses good ability in inhibiting BCL2 and PI3KCG protein in MD simulation. In the end, usnic acid has exhibited good potential in the inhibition of PI3KCG proteins, rather than the other proteins mentioned. Thus further study on structural modification of usnic acid could enhance the ability of usnic acid in the inhibition of PI3KCG as anti-colorectal and anti-small cell lung cancer drug candidate.Communicated by Ramaswamy H. Sarma.

Current Trends and Changes in Use of Membrane Molecular Dynamics Simulations within Academia and the Pharmaceutical Industry

There has been an almost exponential increase in the use of molecular dynamics simulations in basic research and industry over the last 5 years, with almost a doubling in the number of publications each year. Many of these are focused on neurological membranes, and biological membranes in general, applied to the medical industry. A smaller portion have utilized membrane simulations to answer more basic questions related to the function of specific proteins, chemicals or biological processes. This review covers some newer studies, alongside studies from the last two decades, to determine changes in the field. Some of these are basic, while others are more profound, such as multi-component embedded membrane machinery. It is clear that many facets of the discipline remain the same, while the focus on and uses of the technology are broadening in scope and utilization as a general research tool. Analysis of recent literature provides an overview of the current methodologies, covers some of the recent trends or advances and tries to make predictions of the overall path membrane molecular dynamics will follow in the coming years. In general, the overview presented is geared towards the general scientific community, who may wish to introduce the use of these methodologies in light of these changes, making molecular dynamic simulations more feasible for general scientific or medical research.

Insights into glycosidic bond specificity of an engineered selective α-L-rhamnosidase N12-Rha via activity assays and molecular modelling

αL-rhamnosidase (EC 3.2.1.40) has been widely used in food processing and pharmaceutical preparation. The recombinant α-L-rhamnosidase N12-Rha from Aspergillus niger JMU-TS528 had significantly higher catalytic activity on α-1,6 glycosidic bond than α-1,2 glycosidic bond, and had no activity on α-1,3 glycosidic bond. The activities of hydrolyzed hesperidin and naringin were 7240 U/mL and 945 U/mL, respectively, which are 10.63 times that of native α-L-rhamnosidase. The activity could maintain more than 80% at pH 3–6 and 40–60℃. Quantum chemistry calculations showed that charge difference of the C-O atoms of the α-1,2, α-1,3 and α-1,6 bonds indicated that α-1,6 bond is most easily broken and α-1,3 bond is the most stable. Molecular dynamics simulations revealed that the key residue Trp359 that may affect substrate specificity and the main catalytic sites of N12-Rha are located in the (α/α)6-barrel domain.

Effects of High Pressure Processing and Thermal Treatment on the Interaction between α-Lactalbumin and Pelargonium-3-Glucoside

In this study, high pressure processing (HPP) and thermal treatment were comparatively evaluated by examining their impacts on the binding behavior and interaction between α-lactalbumin (α-La) and pelargonium-3-glucoside (P3G) under pH values of 6.0, 7.4, and 8.0. The methods of circular dichroism spectroscopy, fluorescence quenching, dynamic light scattering, and molecular simulation were used to characterize the effects of processing-induced changes in protein structure, size distribution, binding site conformation, and residue charges on their binding characteristics between them. The results indicated that the thermal treatments significantly increased the quenching constants of the complex at pH 7.4/8.0 and 60/80 °C, as well as the accessible fraction of protein at pH 8.0/80 °C. Both HPP and thermal treatments increased the random coil content and showed limited effects on the α-helix and β-sheet contents of α-La and caused the aggregation of the complex to varying degrees. Molecular dynamic simulation and docking analyses revealed that the binding site of the complex did not change under different processing conditions, but the solvent-accessible surface area varied under different conditions.

Designing of a Novel Multi-Antigenic Epitope-Based Vaccine against E. hormaechei: An Intergraded Reverse Vaccinology and Immunoinformatics Approach

Enterobacter hormaechei is involved in multiple hospital-associated infections and is resistant to beta-lactam and tetracycline antibiotics. Due to emerging antibiotics resistance in E. hormaechei and lack of licensed vaccine availability, efforts are required to overcome the antibiotics crisis. In the current research study, a multi-epitope-based vaccine against E. hormaechei was designed using reverse vaccinology and immunoinformatic approaches. A total number of 50 strains were analyzed from which the core proteome was extracted. One extracellular (curlin minor subunit CsgB) and two periplasmic membrane proteins (flagellar basal-body rod protein (FlgF) and flagellar basal body P-ring protein (FlgI) were prioritized for B and T-cell epitope prediction. Only three filtered TPGKMDYTS, GADMTPGKM and RLSAESQAT epitopes were used when designing the vaccine construct. The epitopes were linked via GPGPG linkers and EAAAK linker-linked cholera toxin B-subunit adjuvant was used to enhance the immune stimulation efficacy of the vaccine. Docking studies of the vaccine construct with immune cell receptors revealed better interactions, vital for generating proper immune reactions. Docked complexes of vaccine with MHC-I, MHC-II and Tool-like receptor 4 (TLR-4) reported the lowest binding energy of −594.1 kcal/mol, −706.7 kcal/mol, −787.2 kcal/mol, respectively, and were further subjected to molecular dynamic simulations. Net binding free energy calculations also confirmed that the designed vaccine has a strong binding affinity for immune receptors and thus could be a good vaccine candidate for future experimental investigations.

Molecular Interactions Associated with Coagulation of Organic Pollutants by 2S Albumin of Plant Proteins: A Computational Approach

The removal of organic pollutants is a major challenge in wastewater treatment technologies. Coagulation by plant proteins is a promising technique for this purpose. The use of these proteins has been experimentally investigated and reported in the literature. However, the determination of the molecular interactions of these species is experimentally challenging and the computational approach offers a suitable alternative in gathering useful information for this system. The present study used a molecular dynamic simulation approach to predict the potentials of using Moringa oleifera (MO), Arachis hypogaea, Bertholletia excelsa, Brassica napus, and Helianthus annuus plant proteins for the coagulation of organic pollutants and the possible mechanisms of coagulation of these proteins. The results showed that the physicochemical and structural properties of the proteins are linked to their performance. Maximum coagulation of organic molecules to the proteins is between 50-100%. Among five proteins studied for coagulation, Brassica napus and Helianthus annuus performed better than the well-known MO protein. The amino acid residues interacting with the organic molecules play a significant role in the coagulation and this is peculiar with each plant protein. Hydrogen bond and π-interactions dominate throughout the protein-pollutants molecular interactions. The reusability of the proteins after coagulation derived from their structural quality analysis along with the complexes looks promising and most of them are better than that of the MO. The results showed that the seed proteins studied have good prediction potentials to be used for the coagulation of organic pollutants from the environment, as well as the insights into their molecular activities for bioremediation.

Computational Methods for the Interaction between Cyclodextrins and Natural Compounds: Technology, Benefits, Limitations, and Trends

Cyclodextrins (CDs) have a hollow structure with a hydrophobic interior and hydrophilic exterior. Forming inclusion complexes with CDs will maximize the bioavailability of natural compounds and enable active components to be processed into functional foods, medicines, additives, and so forth. However, experimental methods cannot explain CD-guest binding at the atomic level. Different models have been recently developed to simulate the interaction between CDs and guests to study the binding conformation and analyze noncovalent forces. This review paper summarizes modeling methods of CD-natural compound complexes. The methods include quantitative structure-activity relationships, molecular docking, molecular dynamics simulations, and quantum-chemical calculations. The applications of these methods to enhance the solubility and bioactivities of guest molecules, assist material transportation, and promote compound extraction are also discussed. The purpose of this review is to explore interaction mechanisms of CDs and guests and to help expand new applications of CDs.

Functional Performance of Plant Proteins

Increasingly, consumers are moving towards a more plant-based diet. However, some consumers are avoiding common plant proteins such as soy and gluten due to their potential allergenicity. Therefore, alternative protein sources are being explored as functional ingredients in foods, including pea, chickpea, and other legume proteins. The factors affecting the functional performance of plant proteins are outlined, including cultivars, genotypes, extraction and drying methods, protein level, and preparation methods (commercial versus laboratory). Current methods to characterize protein functionality are highlighted, including water and oil holding capacity, protein solubility, emulsifying, foaming, and gelling properties. We propose a series of analytical tests to better predict plant protein performance in foods. Representative applications are discussed to demonstrate how the functional attributes of plant proteins affect the physicochemical properties of plant-based foods. Increasing the protein content of plant protein ingredients enhances their water and oil holding capacity and foaming stability. Industrially produced plant proteins often have lower solubility and worse functionality than laboratory-produced ones due to protein denaturation and aggregation during commercial isolation processes. To better predict the functional performance of plant proteins, it would be useful to use computer modeling approaches, such as quantitative structural activity relationships (QSAR).

Heat-induced conversion of multiscale molecular structure of natural food nutrients: A review

Thermal process is the most important way of treating foods. Heat energy inputted into the natural food system induces the depolymerization of multi-scale structures of matrix, and causes the intramolecular and intermolecular interactions of different nutrients. It attacks and breaks the original polymeric molecule structures and the functional properties of macronutrients such as carbohydrates, proteins and lipids. Micronutrients such as vitamins and other novel functional ingredients are also thermally converted. The heat-induced conversions of nutrients are slightly or totally with discrepancy in simple-, simulated- and real-food systems, respectively. Thus, this review aims to extensively summarize the heat-induced structural characteristics, thermal conversion pathways and pyrolysis mechanism of nutrients both in simple and complex food matrices. The structural change of each nutrient and its thermal reaction kinetics depend on the molecule structure and polymeric characteristic of the unit substances in the system.

An Overview of Molecular Dynamics Simulation for Food Products and Processes

Molecular dynamics (MD) simulation is a particularly useful technique in food processing. Normally, food processing techniques can be optimized to favor the creation of higher-quality, safer, more functional, and more nutritionally valuable food products. Modeling food processes through the application of MD simulations, namely, the Groningen Machine for Chemical Simulations (GROMACS) software package, is helpful in achieving a better understanding of the structural changes occurring at the molecular level to the biomolecules present in food products during processing. MD simulations can be applied to define the optimal processing conditions required for a given food product to achieve a desired function or state. This review presents the development history of MD simulations, provides an in-depth explanation of the concept and mechanisms employed through the running of a GROMACS simulation, and outlines certain recent applications of GROMACS MD simulations in the food industry for the modeling of proteins in food products, including peanuts, hazelnuts, cow’s milk, soybeans, egg whites, PSE chicken breast, and kiwifruit.

In silico analysis of energy interactions between nociceptin/orfanin FQ receptor and two antagonists with potential antidepressive action

This study addresses the binding energies of NOPR-ligand complexes and presents the main amino acid residues involved in the interaction between these complexes.

Encapsulation of an anticancer drug Isatin inside a host nano-vehicle SWCNT: a molecular dynamics simulation

The use of carbon nanotubes as anticancer drug delivery cargo systems is a promising modality as they are able to perforate cellular membranes and transport the carried therapeutic molecules into the cellular components. Our work describes the encapsulation process of a common anticancer drug, Isatin (1H-indole-2,3-dione) as a guest molecule, in a capped single-walled carbon nanotube (SWCNT) host with chirality of (10,10). The encapsulation process was modelled, considering an aqueous solution, by a molecular dynamics (MD) simulation under a canonical NVT ensemble. The interactions between the atoms of Isatin were obtained from the DREIDING force filed. The storage capacity of the capped SWCNT host was evaluated to quantify its capacity to host multiple Isatin molecules. Our results show that the Isatin can be readily trapped inside the volume cavity of the capped SWCNT and it remained stable, as featured by a reduction in the van der Waals forces between Isatin guest and the SWCNT host (at approximately - 30 kcal mol^-1) at the end of the MD simulation (15 ns). Moreover, the free energy of encapsulation was found to be - 34 kcal mol^-1 suggesting that the Isatin insertion procedure into the SWCNT occurred spontaneously. As calculated, a capped SWCNT (10,10) with a length of 30 Å, was able to host eleven (11) molecules of Isatin, that all remained steadily encapsulated inside the SWCNT volume cavity, showing a potential for the use of carbon nanotubes as drug delivery cargo systems.

Performance-Based Screening of Porous Materials for Carbon Capture

Computational screening methods have changed the way new materials and processes are discovered and designed. For adsorption-based gas separations and carbon capture, recent efforts have been directed toward the development of multiscale and performance-based screening workflows where we can go from the atomistic structure of an adsorbent to its equilibrium and transport properties at different scales, and eventually to its separation performance at the process level. The objective of this work is to review the current status of this new approach, discuss its potential and impact on the field of materials screening, and highlight the challenges that limit its application. We compile and introduce all the elements required for the development, implementation, and operation of multiscale workflows, hence providing a useful practical guide and a comprehensive source of reference to the scientific communities who work in this area. Our review includes information about available materials databases, state-of-the-art molecular simulation and process modeling tools, and a complete catalogue of data and parameters that are required at each stage of the multiscale screening. We thoroughly discuss the challenges associated with data availability, consistency of the models, and reproducibility of the data and, finally, propose new directions for the future of the field.

Theoretical Encapsulation of Fluorouracil (5-FU) Anti-Cancer Chemotherapy Drug into Carbon Nanotubes (CNT) and Boron Nitride Nanotubes (BNNT)

INTRODUCTION

Chemotherapy with anti-cancer drugs is considered the most common approach for killing cancer cells in the human body. However, some barriers such as toxicity and side effects would limit its usage. In this regard, nano-based drug delivery systems have emerged as cost-effective and efficient for sustained and targeted drug delivery. Nanotubes such as carbon nanotubes (CNT) and boron nitride nanotubes (BNNT) are promising nanocarriers that provide the cargo with a large inner volume for encapsulation. However, understanding the insertion process of the anti-cancer drugs into the nanotubes and demonstrating drug-nanotube interactions starts with theoretical analysis.

METHODS

First, interactions parameters of the atoms of 5-FU were quantified from the DREIDING force field. Second, the storage capacity of BNNT (8,8) was simulated to count the number of drugs 5-FU encapsulated inside the cavity of the nanotubes. In terms of the encapsulation process of the one drug 5-FU into nanotubes, it was clarified that the drug 5-FU was more rapidly adsorbed into the cavity of the BNNT compared with the CNT due to the higher van der Waals (vdW) interaction energy between the drug and the BNNT.

RESULTS

The obtained values of free energy confirmed that the encapsulation process of the drug inside the CNT and BNNT occurred spontaneously with the free energies of -14 and -25 kcal·mol-1, respectively.

DISCUSSION

However, the lower value of the free energy in the system containing the BNNT unraveled more stability of the encapsulated drug inside the cavity of the BNNT comparing the system having CNT. The encapsulation of Fluorouracil (5-FU) anti-cancer chemotherapy drug (commercial name: Adrucil®) into CNT (8,8) and BNNT (8,8) with the length of 20 Å in an aqueous solution was discussed herein applying molecular dynamics (MD) simulation.

α-Helical Antimicrobial Peptide Encapsulation and Release from Boron Nitride Nanotubes: A Computational Study

Introduction

Antimicrobial peptides are potential therapeutics as anti-bacteria, anti-viruses, anti-fungi, or anticancers. However, they suffer from a short half-life and drug resistance which limit their long-term clinical usage.

Methods

Herein, we captured the encapsulation of antimicrobial peptide HA-FD-13 into boron nitride nanotube (BNNT) (20,20) and its release due to subsequent insertion of BNNT (14,14) with molecular dynamics simulation.

Results

The peptide-BNNT (20,20) van der Waals (vdW) interaction energy decreased to −270 kcal·mol⁻¹ at the end of the simulation (15 ns). However, during the period of 0.2–1.8 ns, when half of the peptide was inside the nanotube, the encapsulation was paused due to an energy barrier in the vicinity of BNNT and subsequently the external intervention, such that the self-adjustment of the peptide allowed full insertion. The free energy of the encapsulation process was −200.12 kcal·mol⁻¹, suggesting that the insertion procedure occurred spontaneously.

Discussion

Once the BNNT (14,14) entered into the BNNT (20,20), the peptide was completely released after 83.8 ps. This revealed that the vdW interaction between the BNNT (14,14) and BNNT (20,20) was stronger than between BNNT (20,20) and the peptide; therefore, the BNNT (14,14) could act as a piston pushing the peptide outside the BNNT (20,20). Moreover, the sudden drop in the vdW energy between nanotubes to the value of the −1300 Kcal·mol⁻¹ confirmed the self-insertion of the BNNT (14,14) into the BNNT (20,20) and correspondingly the release of the peptide.

Molecular dynamics simulation for mechanism revelation of the safety and nutrition of lipids and derivatives in food: State of the art

Molecular dynamics (MD) simulation has proved to be a powerful tool in the study of proteins, nucleic acids, lipids, and carbohydrates et al. in fields of health, nutrition, and food science. In particular, MD simulation has been employed in the investigation of various lipid systems such as triglycerides, phospholipid membranes, etc. Due to the continuous updating of computing resources and the development of new MD simulation methods and force field parameters, the simulation's time and size scale of lipids system has increased by several orders of magnitude. However, MD simulation cannot be used for systems invovle chemical reactions. These greatly limit its further application in the field of lipid research. This paper reviews the progress and development of MD simulation, especially for the application of MD simulation in different lipid systems. In this paper, MD simulation and its general workflow was briefly introduced firstly. Subsequently, the application of MD simulation in various lipid systems was reviewed in-depth. Finally, the limitation and future prospects of MD simulation in lipid research were also discussed. This review provided new insights into the investigation of MD simulation, and a novel thought for lipid study. We believe that MD simulation will exhibit more and more great advantages in the investigation of lipids in the future due to the development of novlel methods.

Boron Nitride Nanotube as an Antimicrobial Peptide Carrier: A Theoretical Insight

INTRODUCTION

Nanotube-based drug delivery systems have received considerable attention because of their large internal volume to encapsulate the drug and the ability to penetrate tissues, cells, and bacteria. In this regard, understanding the interaction between the drug and the nanotube to evaluate the encapsulation behavior of the drug in the nanotube is of crucial importance.

METHODS

In this work, the encapsulation process of the cationic antimicrobial peptide named cRW3 in the biocompatible boron nitride nanotube (BNNT) was investigated under the Canonical ensemble (NVT) by molecular dynamics (MD) simulation.

RESULTS

The peptide was absorbed into the BNNT by van der Waals (vdW) interaction between cRW3 and the BNNT, in which the vdW interaction decreased during the simulation process and reached the value of -142.7 kcal·mol-1 at 4 ns.

DISCUSSION

The increase in the potential mean force profile of the encapsulated peptide during the pulling process of cRW3 out of the nanotube showed that its insertion into the BNNT occurred spontaneously and that the inserted peptide had the desired stability. The energy barrier at the entrance of the BNNT caused a pause of 0.45 ns when half of the peptide was inside the BNNT during the encapsulation process. Therefore, during this period, the peptide experienced the weakest movement and the smallest conformational changes.

Structural and functional changes of the protein β-lactoglobulin under thermal and electrical processing conditions

In the present study we have tried to explore the effect of static external electric field of strength 3.0 V/nm on the conformational changes adopted by the protein β-lactoglobulin. We have chosen different temperatures viz. 300 K, 400 K and 450 K to evaluate the temperature dependent effect of electric field. We have observed that combined effect of high temperature and static external electric field show significant changes on the structural conformation of the protein which in turn may affect the functional properties of the protein. Calculations of root mean square deviations reveal that both helical and β-sheet regions of the protein are noticeably affected at high temperature. We have used solvent accessible surface area (SASA) and dipole moment values to explain that there is changes in hydrophobicity of the protein surface due to presence of external electric field. The study reveals that electric field in combination with high temperature can be used to alter the conformation of the protein and the effect of external electric field is more pronounced at high temperature than that of low temperature. The study provides a better understanding of the conformational changes adopted by the protein under the stress of external electric field and high temperature and provide guidance to choose optimum conditions for processing without loss of nutritional properties.

Sensorial Perception of Astringency: Oral Mechanisms and Current Analysis Methods

Understanding consumers' food choices and the psychological processes involved in their preferences is crucial to promote more mindful eating regulation and guide food design. Fortifying foods minimizing the oral dryness, rough, and puckering associated with many functional ingredients has been attracting interest in understanding oral astringency over the years. A variety of studies have explored the sensorial mechanisms and the food properties determining astringency perception. The present review provides a deeper understanding of astringency, a general view of the oral mechanisms involved, and the exciting variety of the latest methods used to direct and indirectly quantify and simulate the astringency perception and the specific mechanisms involved.

Functional Properties of Extracted Protein from Edible Insect Larvae and Their Interaction with Transglutaminase

Global concern about food supply shortage has increased interest on novel food sources. Among them, edible insects have been studied as a potential major food source. This study aimed to improve the functional properties of protein solutions extracted from Protaetia brevitarsis (PB) by use of transglutaminase (TG) as a cross-linking agent. After various incubation times (10, 20, 30, 60, and 90 min) with TG, the protein solutions were assessed with regard to their amino acid composition, protein nutritional quality, pH, color (yellowness), molecular weight distribution, thermal stability, foam ability (capacity and stability), and emulsion ability (capacity and stability). Incubation with TG changed the amino acid composition of the proteins and shifted the molecular weight distribution towards higher values, while improving the rest of the aforementioned properties. Since the incubation time for 90 min decreased the protein functionality, the optimum incubation time for cross-linking PB-derived protein with TG is 60 min. The application of TG to edible insect proteins ultimately increases its functionality and allows for the development of novel insect processing technology.

Simulations of Temperature and Pressure Unfolding in Soy Allergen Gly m 4 Using Molecular Modeling

This study aims at understanding the changes in the structural properties of Gly m 4 soy allergen as a result of the influence of the temperature and pressure deviations. The primary emphasis was placed on analyzing the surface properties of suspected linear and conformational epitopes present in the Gly m 4 allergen. All three epitopes of Gly m 4 were studied, and the results showed that the molecule has significant structural changes in terms of solvent-accessible surface area (SASA) and radius of gyration, which showed that the increased pressures resulted in compaction. However, at lower temperatures and higher pressures (300 K and 6 kbar), swelling in the molecule was observed with a significant increase in the surface area. The study also tracked the changes in surface areas of individual residues that are part of the selected epitopes. Residues, such as D-27 and T-51, were found to have significant changes in their SASA as a result of temperature and pressure deviations.