Unraveling Protein-Metabolite Interactions in Precision Nutrition: A Case Study of Blueberry-Derived Metabolites Using Advanced Computational Methods.

Metabolomics, the study of small-molecule metabolites within biological systems, has become a potent instrument for understanding cellular processes. Despite its profound insights into health, disease, and drug development, identifying the protein partners for metabolites, especially dietary phytochemicals, remains challenging. In the present study, we introduced an innovative in silico, structure-based target prediction approach to efficiently predict protein targets for metabolites. We analyzed 27 blood serum metabolites from nutrition intervention studies' blueberry-rich diets, known for their health benefits, yet with elusive mechanisms of action. Our findings reveal that blueberry-derived metabolites predominantly interact with Carbonic Anhydrase (CA) family proteins, which are crucial in acid-base regulation, respiration, fluid balance, bone metabolism, neurotransmission, and specific aspects of cellular metabolism. Molecular docking showed that these metabolites bind to a common pocket on CA proteins, with binding energies ranging from -5.0 kcal/mol to -9.0 kcal/mol. Further molecular dynamics (MD) simulations confirmed the stable binding of metabolites near the Zn binding site, consistent with known compound interactions. These results highlight the potential health benefits of blueberry metabolites through interaction with CA proteins.

A versatile delivery vehicle for cellular oxygen and fuels or metabolic sensor? A review and perspective on the functions of myoglobin.

A canonical view of the primary physiological function of myoglobin (Mb) is that it is an oxygen (O2) storage protein supporting mitochondrial oxidative phosphorylation, especially as the tissue O2 partial pressure (Po2) drops and Mb off-loads O2. Besides O2 storage/transport, recent findings support functions for Mb in lipid trafficking and sequestration, interacting with cellular glycolytic metabolites such as lactate (LAC) and pyruvate (PYR), and "ectopic" expression in some types of cancer cells and in brown adipose tissue (BAT). Data from Mb knockout (Mb-/-) mice and biochemical models suggest additional metabolic roles for Mb, especially regulation of nitric oxide (NO) pools, modulation of BAT bioenergetics, thermogenesis, and lipid storage phenotypes. From these and other findings in the literature over many decades, Mb's function is not confined to delivering O2 in support of oxidative phosphorylation but may serve as an O2 sensor that modulates intracellular Po2- and NO-responsive molecular signaling pathways. This paradigm reflects a fundamental change in how oxidative metabolism and cell regulation are viewed in Mb-expressing cells such as skeletal muscle, heart, brown adipocytes, and select cancer cells. Here, we review historic and emerging views related to the physiological roles for Mb and present working models illustrating the possible importance of interactions between Mb, gases, and small-molecule metabolites in regulation of cell signaling and bioenergetics.

Computational Analysis Reveals Unique Binding Patterns of Oxygenated and Deoxygenated Myoglobin to the Outer Mitochondrial Membrane.

Myoglobin (Mb) interaction with the outer mitochondrial membrane (OMM) promotes oxygen (O2) release. However, comprehensive molecular details on specific contact regions of the OMM with oxygenated (oxy-) and deoxygenated (deoxy-)Mb are missing. We used molecular dynamics (MD) simulations to explore the interaction of oxy- and deoxy-Mb with the membrane lipids of the OMM in two lipid compositions: (a) a typical whole membrane on average, and (b) specifically the cardiolipin-enriched cristae region (contact site). Unrestrained relaxations showed that on average, both the oxy- and deoxy-Mb established more stable contacts with the lipids typical of the cristae contact site, then with those of the average OMM. However, in steered detachment simulations, deoxy-Mb clung more tightly to the average OMM, and oxy-Mb strongly preferred the contact sites of the OMM. The MD simulation analysis further indicated that a non-specific binding, mediated by local electrostatic interactions, existed between charged or polar groups of Mb and the membrane, for stable interaction. To the best of our knowledge, this is the first computational study providing the molecular details of the direct Mb-mitochondria interaction that assisted in distinguishing the preferred localization of oxy- and deoxy-Mb on the OMM. Our findings support the existing experimental evidence on Mb-mitochondrial association and shed more insights on Mb-mediated O2 transport for cellular bioenergetics.

Exploring GPR109A Receptor Interaction with Hippuric Acid Using MD Simulations and CD Spectroscopy.

Previous research has indicated that various metabolites belonging to phenolic acids (PAs), produced by gut microflora through the breakdown of polyphenols, help in promoting bone development and protecting bone from degeneration. Results have also suggested that G-protein-coupled receptor 109A (GPR109A) functions as a receptor for those specific PAs such as hippuric acid (HA) and 3-(3-hydroxyphenyl) propionic acid (3-3-PPA). Indeed, HA has a molecular structural similarity with nicotinic acid (niacin) which has been shown previously to bind to GPR109A receptor and to mediate antilipolytic effects; however, the binding pocket and the structural nature of the interaction remain to be recognized. In the present study, we employed a computational strategy to elucidate the molecular structural determinants of HA binding to GPR109A and GPR109B homology models in understanding the regulation of osteoclastogenesis. Based on the docking and molecular dynamics simulation studies, HA binds to GPR109A similarly to niacin. Specifically, the transmembrane helices 3, 4 and 6 (TMH3, TMH4 and TMH6) and Extracellular loop 1 and 2 (ECL1 and ECL2) residues of GRP109A; R111 (TMH3), K166 (TMH4), ECL2 residues; S178 and S179, and R251 (TMH6), and residues of GPR109B; Y87, Y86, S91 (ECL1) and C177 (ECL2) contribute for HA binding. Simulations and Molecular Mechanics Poisson-Boltzmann solvent accessible area (MM-PBSA) calculations reveal that HA has higher affinity for GPR109A than for GPR109B. Additionally, in silico mutation analysis of key residues have disrupted the binding and HA exited out from the GPR109A protein. Furthermore, measurements of time-resolved circular dichroism spectra revealed that there are no major conformational changes in the protein secondary structure on HA binding. Taken together, our findings suggest a mechanism of interaction of HA with both GPR109A and GPR109B receptors.

Myoglobin-Pyruvate Interactions: Binding Thermodynamics, Structure-Function Relationships, and Impact on Oxygen Release Kinetics.

Myoglobin (Mb), besides its roles as an oxygen (O2) carrier/storage protein and nitric oxide NO scavenger/producer, may participate in lipid trafficking and metabolite binding. Our recent findings have shown that O2 is released from oxy-Mb upon interaction with lactate (LAC, anerobic glycolysis end-product). Since pyruvate (PYR) is structurally similar and metabolically related to LAC, we investigated the effects of PYR (aerobic glycolysis end-product) on Mb using isothermal titration calorimetry, circular dichroism, and O2-kinetic studies to evaluate PYR affinity toward Mb and to compare the effects of PYR and LAC on O2 release kinetics of oxy-Mb. Similar to LAC, PYR interacts with both oxy- and deoxy-Mb with a 1:1 stoichiometry. Time-resolved circular dichroism spectra revealed that there are no major conformational changes in the secondary structures of oxy- or deoxy-Mb during interactions with PYR or LAC. However, we found contrasting results with respect to binding affinities and substrate preference, where PYR has higher affinity toward deoxy-Mb when compared with LAC (which prefers oxy-Mb). Furthermore, PYR interaction with oxy-Mb releases a significantly lower amount of O2 than LAC. Taken together, our findings support the hypothesis that glycolytic end-products play a distinctive role in the Mb-rich tissues by serving as novel regulators of O2 availability, and/or by impacting other activities related to oxy-/deoxy-Mb toggling in resting vs. exercised or metabolically activated conditions.

Computational insights on molecular interactions of acifran with GPR109A and GPR109B.

Acifran is a well-known agonist of G-protein-coupled receptor protein, namely GPR109A. Acifran is primarily used in the treatment of dyslipidemia, myocardial infractions, and atherosclerosis in humans due to its lower vascular and metabolic side effects. However, experimental and computational studies on interaction sites of acifran with GPR proteins (GPR109A and GPR109B) are lacking. Our computational studies using docking and molecular dynamics simulation revealed that acifran binds distinctly to both GPR109A and GPR109B, but with lower affinity to the latter. The weak binding of acifran-GPR109B is mainly due to the presence of residues S91 and N94 in ECL1 and I178 amino acid in ECL2 region of GPR109B, whereas R111 and R251 residues in TMH3 and TMH6 are crucial for GPR109A-acifran complex stability. Additionally, molecular mechanics/Poisson-Boltzmann solvent accessible surface area (MM/PBSA) analysis revealed that both GPR109A- and GPR109B-acifran complexes are energetically stable with lower calculated binding free energy values for the latter. Energy-minimized structures of GPR109A-acifran and GPR109B-acifran complex.

Myoglobin Interaction with Lactate Rapidly Releases Oxygen: Studies on Binding Thermodynamics, Spectroscopy, and Oxygen Kinetics.

Myoglobin (Mb)-mediated oxygen (O2) delivery and dissolved O2 in the cytosol are two major sources that support oxidative phosphorylation. During intense exercise, lactate (LAC) production is elevated in skeletal muscles as a consequence of insufficient intracellular O2 supply. The latter results in diminished mitochondrial oxidative metabolism and an increased reliance on nonoxidative pathways to generate ATP. Whether or not metabolites from these pathways impact Mb-O2 associations remains to be established. In the present study, we employed isothermal titration calorimetry, O2 kinetic studies, and UV-Vis spectroscopy to evaluate the LAC affinity toward Mb (oxy- and deoxy-Mb) and the effect of LAC on O2 release from oxy-Mb in varying pH conditions (pH 6.0-7.0). Our results show that LAC avidly binds to both oxy- and deoxy-Mb (only at acidic pH for the latter). Similarly, in the presence of LAC, increased release of O2 from oxy-Mb was detected. This suggests that with LAC binding to Mb, the structural conformation of the protein (near the heme center) might be altered, which concomitantly triggers the release of O2. Taken together, these novel findings support a mechanism where LAC acts as a regulator of O2 management in Mb-rich tissues and/or influences the putative signaling roles for oxy- and deoxy-Mb, especially under conditions of LAC accumulation and lactic acidosis.

Structural and Energetic Insights Into the Interaction of Niacin With the GPR109A Receptor.

The transmembrane G-protein coupled receptor GPR109A has been previously shown to function as a receptor for niacin in mediating antilipolytic effects. Although administration of high doses of niacin has shown beneficial effects on lipid metabolism, however, it is often accompanied by disturbing side effects such as flushing, liver damage, glucose intolerance, or gastrointestinal problems. Thus, it is important to understand niacin-GPR109A interactions, which can be beneficial for the development of alternate drugs having antilipolytic effects with less or no side effects. To get into the structural insights on niacin binding to GPR109A, we have performed 100 nanoseconds long all-atom MD simulations of five niacin-GPR109A complexes (automatically docked pose 0, and randomly placed niacin in poses 1 to 4 in the receptor crevice) and analyzed using binding free energy calculations and H-bond analysis. Steered MD simulations were used to get an average force for niacin translocation between the bulk and the external crevice of the wild type and mutant (N86Y, W91 S, S178I, and triple mutant of all three residues) GPR109A receptors, as well as GPR109B (as a control that does not bind niacin). The H-bond analysis revealed that TMH3 residue R111 interacts with niacin in a total of 4 (poses 0 to 3) complexes, while residues C177, S178, and S179 contact niacin in complex pose 4, and all these complexes were energetically stable. According to steered MD simulations, all the GPR109A mutants and GPR109B required greater force than that of wild-type GPR109A to translocate in the external crevice, suggesting increased sterical obstacles. Thus, the residues N86 (at the junction of TMH2/ECL2), W91 (ECL2), R111 (TMH3), and ECL3 residues (C177, S178, S179) play an important role for optimal routing of niacin entry and to bind GPR109A.