Structural basis of the Tanford transition of bovine beta-lactoglobulin

Comparison and analysis of mechanism of β-lactoglobulin self-assembled gel carriers formed by different gelation methods

Based on the findings of our previous studies, a comprehensive comparative investigation of the quality and formation mechanism of gels obtained from protein self-assemblies induced by different methods is necessary. Self-assembled heat-induced gels had higher gel mechanical strength, and hydrophobic interactions played a greater role. Whether or not heat treatment was used to induce gel formation may play a more important role than the effect of divalent cations on gel formation. Hydrogen bonds played an important role in all gels formed using different gelation methods. Furthermore, Self-assembled cold-induced gels were considered to can load bioactive substances with different hydrophilicity properties due to the high water-holding capacity and the smooth, dense microstructure. Therefore, β-lactoglobulin fibrous and worm-like self-assembled cold-induced gels as a delivery material for hydrophilic bioactive substances (epigallocatechin gallate, vitamin B2) and amphiphilic bioactive substance (naringenin), with good encapsulation efficiency (91.92 %, 97.08 %, 96.72 %, 96.52 %, 98.94 %, 97.41 %, respectively) and slow-release performance.

Binding of Per- and Polyfluoroalkyl Substances to β-Lactoglobulin from Bovine Milk

Per- and polyfluoroalkyl substances (PFAS) are known for their high environmental persistence and potential toxicity. The presence of PFAS has been reported in many dairy products. However, the mechanisms underlying the accumulation of PFAS in these products remain unclear. Here, we used native mass spectrometry and molecular dynamics simulations to probe the interactions between 19 PFAS of environmental concern and two isoforms of the major bovine whey protein β-lactoglobulin (β-LG). We observed that six of these PFAS bound to both protein isoforms with low- to mid-micromolar dissociation constants. Based on quantitative, competitive binding experiments with endogenous ligands, PFAS can bind orthosterically and preferentially to β-LG's hydrophobic ligand-binding calyx. β-Cyclodextrin can also suppress binding of PFAS to β-LG owing to the ability of β-cyclodextrin to directly sequester PFAS from solution. This research sheds light on PFAS-β-LG binding, suggesting that such interactions could impact lipid-fatty acid transport in bovine mammary glands at high PFAS concentrations. Furthermore, our results highlight the potential use of β-cyclodextrin in mitigating PFAS binding, providing insights toward the development of strategies to reduce PFAS accumulation in dairy products and other biological systems.

Integrated spectroscopic and computational analyses unravel the molecular interaction of pesticide azinphos-methyl with bovine beta-lactoglobulin

Organophosphorus are typically hazardous chemicals used in the pharmaceutical, agricultural, and other industries. They pose a serious risk to human life and can be fatal upon direct exposure. Hence, studying the interaction between such compounds with proteins is crucial for environmental, health, and food safety. In this study, we investigated the interaction mechanism between azinphos-methyl (AZM) and β-lactoglobulin (BLG) at pH 7.4 using a combination of biophysical techniques. Intrinsic fluorescence investigations revealed that BLG fluorescence was quenched in the presence of increasing AZM concentrations. The quenching mechanism was identified as static, as evidenced by a decrease in the fluorescence quenching constant (1.25 × 104, 1.18 × 104, and 0.86 × 104 M-1) with an increase in temperatures. Thermodynamic calculations (ΔH > 0; ΔS > 0) affirmed the formation of a complex between AZM and BLG through hydrophobic interactions. The BLG's secondary structure was found to be increased due to AZM interaction. Ultraviolet -visible spectroscopy data showed alterations in BLG conformation in the presence of AZM. Molecular docking highlighted the significant role of hydrophobic interactions involving residues such as Val43, Ile56, Ile71, Val92, Phe105, and Met107 in the binding between BLG and AZM. A docking energy of -6.9 kcal mol-1, and binding affinity of 1.15 × 105 M-1 suggest spontaneous interaction between AZM and BLG with moderate to high affinity. These findings underscore the potential health risks associated with the entry of AZM into the food chain, emphasizing the need for further consideration of its impact on human health.

Reactivity and mechanism of the reactions of 4-methylbenzoquinone with amino acid residues in β-lactoglobulin: A kinetic and product investigation

Quinones, produced by the oxidation of phenolic compounds, covalently bind to nucleophilic groups on amino acids or proteins. In this study, the reactions of 4-methylbenzoquinone (4MBQ) with β-lactoglobulin (β-LG) and amino acids at neutral pH were investigated. LC-MS analysis revealed that Cys121 was likely the most modified residue in β-LG. Identification of reaction products by LC-MS/MS showed that Michael addition occurred in all reactions with amino acids tested. The formation of Schiff base and a di-adduct was found in His and Trp samples. Apparent second-order rate constants (k2) were determined at 25 °C and pH 7.0 by stopped-flow spectrophotometry. The rate of reactions decreased in the order: β-LG > His > Trp > Arg > Nα-acetyl His > Nα-acetyl Arg > Nα-acetyl Trp. The rate constants correlated with the pKa values of the amino acids, showing that the amount of unprotonated amine is the major factor determining the reactivity.

β-Lactoglobulin variants as potential carriers of pramoxine: Comprehensive structural and biophysical studies

β-Lactoglobulin (BLG) is a member of the lipocalin family. As other proteins from this group, BLG can be modified to bind specifically compounds of medical interests. The aim of this study was to evaluate the role of two mutations, L39Y and L58F, in the binding of topical anesthetic pramoxine (PRM) to β-lactoglobulin. Circular dichroism spectroscopy, isothermal titration calorimetry (ITC), and X-ray crystallography were used to understand the mechanisms of BLG-PRM interactions. Studies were performed for three new BLG mutants: L39Y, L58F, and L39Y/L58F. ITC measurements indicated a significant increase in the affinity to the PRM of variants L58F and L39Y. Measurements taken for the double mutant L39Y/L58F showed the additivity of two mutations leading to about 80-fold increase in the affinity to PRM in comparison to natural protein BLG from bovine milk. The determined crystal structures revealed that pramoxine is accommodated in the β-barrel interior of BLG mutants and stabilized by hydrophobic interactions. The observed additive effect of two mutations on drug binding opens the possibility for further designing of new BLG variants with high affinity to selected drugs.

Modeling protein-nucleic acid complexes with extremely large conformational changes using Flex-LZerD

Proteins and nucleic acids are key components in many processes in living cells, and interactions between proteins and nucleic acids are often crucial pathway components. In many cases, large flexibility of proteins as they interact with nucleic acids is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D atomic structures of such protein-nucleic acid complexes. When such structures are not yet experimentally determined, protein docking can be used to computationally generate useful structure models. However, such docking has long had the limitation that the consideration of flexibility is usually limited to small movements or to small structures. We previously developed a method of flexible protein docking which could model ordered proteins which undergo large-scale conformational changes, which we also showed was compatible with nucleic acids. Here, we elaborate on the ability of that pipeline, Flex-LZerD, to model specifically interactions between proteins and nucleic acids, and demonstrate that Flex-LZerD can model more interactions and types of conformational change than previously shown.

A model study on the site-specificity of (-)-epicatechin-induced reactions in β-lactoglobulin by high-resolution mass spectrometry in combination with bioinformatics

Polyphenol-protein reactions in model solutions of β-lactoglobulin (β-LG) incubated with (-)-epicatechin at 37 °C and 60 °C were monitored by microLC-timsTOF Pro-MS/MS combined with bioinformatics strategies. The addition of (-)-epicatechin to the model solutions resulted in changes in tryptic peptide profiles. Covalent bond formation between (-)-epicatechin o-quinones and β-LG was identified for the residues S27, S30, K60, C66, K69, and C160, with C160 being the predominant binding site. Furthermore, the incubation of β-LG with (-)-epicatechin significantly promoted oxidation, especially for the residues M7 and M24. The reaction of monomeric (-)-epicatechino-quinone at C160 was also identified in the milk chocolate sample. The adaptation of this study by extending the scope of the reaction products offers significant potential for comprehensive food profiling strategies.

The protein dynamics of bovine and caprine β-lactoglobulin differ as a function of pH

The properties of milk proteins differ between mammalian species. β-Lactoglobulin (βlg) proteins from caprine and bovine milk are sequentially and structurally highly similar, yet their physicochemical properties differ, particularly in response to pH. To resolve this conundrum, we compared the dynamics of both the monomeric and dimeric states for each homologue at pH 6.9 and 7.5 using hydrogen/deuterium exchange experiments. At pH 7.5, the rate of exchange is similar across both homologues, but at pH 6.9 the dimeric states of the bovine βlg B variant homologue have significantly more conformational flexibility compared with caprine βlg. Molecular dynamics simulations provide a mechanistic rationale for the experimental observations, revealing that variant-specific substitutions encode different conformational ensembles with different dynamic properties consistent with the hydrogen/deuterium exchange experiments. Understanding the dynamic differences across βlg homologues is essential to understand the different responses of these milks to processing, human digestion, and differences in immunogenicity.

Characterization of the interaction of metal-protoporphyrins photosensitizers with β- lactoglobulin

We investigated the interaction of a series of metal-protoporphyrins (PPIXs) with bovine β- lactoglobulin (BLG) using a combination of optical spectroscopy and computational simulations. Unlike other studies, the simulations were not merely used to rationalize the experimental data but were employed to refine the experimental data itself. The study was carried out at two pH values, 5 and 9, where BLG is known to have different conformation dictated by the so-called Tanford transition which occurs near pH 7.5. The transition is postulated to regulate access to the interior binding cavity of the protein, thus the pH variation was used as a parameter to investigate whether PPIXs access the central cavity of BLG. The results of our study show that indeed binding increases significantly at alkaline pH, however, the increased affinity is not due to the accessibility of the central cavity. Instead, binding appears to be determined by the tendency of PPIXs to form large inhomogeneous aggregates at acidic pH which hinders interactions with proteins. The binding site determined through a combination of experimental and computational methods is located at the interface between two BLG monomers where the long α-helix segment of the protein face each other. This region is rich in positively charged Lys residues that interact with the propionic acid chains of the protoporphyrins. Establishing the modality of binding between protoporphyrins and BLG would have important consequences for the use of BLG:PPIX complexes in applications such as artificial photoreceptors, artificial metallo-enzymes, delivery of photosensitizers for phototherapy and even solar energy conversion.

Gas-Phase Unfolding of Protein Complexes Distinguishes Conformational Isomers

Proteins can adopt different conformational states that are important for their biological function and, in some cases, can be responsible for their dysfunction. The essential roles that proteins play in biological systems make distinguishing the structural differences between these conformational states both fundamentally and practically important. Here, we demonstrate that collision-induced unfolding (CIU), in combination with ion mobility-mass spectrometry (IM-MS) measurements, distinguish subtly different conformational states for protein complexes. Using the open and closed states of the β-lactoglobulin (βLG) dimer as a model, we show that these two conformational isomers unfold during collisional activation to generate distinct states that are readily separated by IM-MS. Extensive molecular modeling of the CIU process reproduces the distinct unfolding intermediates and identifies the molecular details that explain why the two conformational states unfold in distinct ways. Strikingly, the open conformational state forms new electrostatic interactions upon collisional heating, while the closed state does not. These newly formed electrostatic interactions involve residues on the loop differentially positioned in the two βLG conformational isomers, highlighting that gas-phase unfolding pathways reflect aspects of solution structure. This combination of experiment and theory provides a path forward for distinguishing subtly different conformational isomers for protein complexes via gas-phase unfolding experiments. Our results also have implications for understanding how protein complexes dissociate in the gas phase, indicating that current models need to be refined to explain protein complex dissociation.

Computational Study of the pH-Dependent Ionic Environment around β-Lactoglobulin

Ions are involved in multiple biological processes and may exist bound to biomolecules or may be associated with their surface. Although the presence of ions in nucleic acids has traditionally gained more interest, ion-protein interactions, often with a marked dependency on pH, are beginning to gather attention. Here we present a detailed analysis on the binding and distribution of ions around β-lactoglobulin using a constant-pH MD (CpHMD) method, at a pH range 3-8, and compare it with the more traditional Poisson-Boltzmann (PB) model and the existing experimental data. Most analyses used ion concentration maps built around the protein, obtained from either the CpHMD simulations or PB calculations. The requirements of approximate charge neutrality and ionic strength equal to bulk, imposed on the MD box, imply that the absolute value of the ion excess should be half the protein charge, which is in agreement with experimental observation on other proteins ( Proc. Natl. Acad. Sci. U.S.A. 2021, 118, e2015879118) and lends support to this protocol. In addition, the protein total charge (including territorially bound ions) estimated with MD is in excellent agreement with electrophoretic measurements. Overall, the CpHMD simulations show good agreement with the nonlinear form of the PB (NLPB) model but not with its linear form, which involves a theoretical inconsistency in the calculation of the concentration maps. In several analyses, the observed pH-dependent trends for the counterions and co-ions are those generally expected, and the ion concentration maps correctly converge to the bulk ionic strength as one moves away from the protein. Despite the overall similarity, the CpHMD and NLPB approaches show some discrepancies when analyzed in more detail, which may be related to an apparent overestimation of counterion excess and underestimation of co-ion exclusion by the NLPB model, particularly at short distances from the protein.

Domain-Based Protein Docking with Extremely Large Conformational Changes

Proteins are key components in many processes in living cells, and physical interactions with other proteins and nucleic acids often form key parts of their functions. In many cases, large flexibility of proteins as they interact is key to their function. To understand the mechanisms of these processes, it is necessary to consider the 3D structures of such protein complexes. When such structures are not yet experimentally determined, protein docking has long been present to computationally generate useful structure models. However, protein docking has long had the limitation that the consideration of flexibility is usually limited to very small movements or very small structures. Methods have been developed which handle minor flexibility via normal mode or other structure sampling, but new methods are required to model ordered proteins which undergo large-scale conformational changes to elucidate their function at the molecular level. Here, we present Flex-LZerD, a framework for docking such complexes. Via partial assembly multidomain docking and an iterative normal mode analysis admitting curvilinear motions, we demonstrate the ability to model the assembly of a variety of protein-protein and protein-nucleic acid complexes.

Red emitting fluorogenic dye as an efficient turn-on probe for milk allergen

Development of simple, fast and non-destructive technique such as fluorescence based method for the quantification of milk allergens in various dairy products is a highly rewarding task. In this contribution, a red emitting fluorogenic dye, quinaldine red (QR) is reported for the detection and quantification of a milk allergen, beta lactoglobulin (β-LG) in milk and whey matrices, utilizing its high selectivity and sensitivity towards β-LG. Detail spectroscopic investigation reveals that binding of QR to the hydrophobic calyx site of β-LG protein substantially reduces the torsional agility and propensity of TICT state formation of QR, rendering the dye highly fluorescent in nature. This enables estimation of β-LG with LOD 52.1(±0.9) nM in buffer solution and 0.21(±0.01) μM in 5 % bovine milk matrix respectively. Additionally, high selectivity and sensitivity, excellent repeatability, quick response, and emission in the biologically favorable red spectral region make QR based fluorometric quantification of β-LG a highly attractive choice. Finally, the estimated β-LG concentrations in milk and whey matrices from fluorometric titration and densitometry methods are found to match excellently with each other, suggesting potential of QR as an efficient turn-on fluorescent probe for the quantification of β-LG (milk allergen) in various dairy products.

New ligand-binding sites identified in the crystal structures of β-lactoglobulin complexes with desipramine

The homodimeric β-lactoglobulin belongs to the lipocalin family of proteins that transport a wide range of hydrophobic molecules and can be modified by mutagenesis to develop specificity for novel groups of ligands. In this work, new lactoglobulin variants, FAF (I56F/L39A/M107F) and FAW (I56F/L39A/M107W), were produced and their interactions with the tricyclic drug desipramine (DSM) were studied using X-ray crystallography, calorimetry (ITC) and circular dichroism (CD). The ITC and CD data showed micromolar affinity of the mutants for DSM and interactions according to the classical one-site binding model. However, the crystal structures unambiguously showed that the FAF and FAW dimers are capable of binding DSM not only inside the β-barrel as expected, but also at the dimer interface and at the entrance to the binding pocket. The presented high-resolution crystal structures therefore provide important evidence of the existence of alternative ligand-binding sites in the β-lactoglobulin molecule. Analysis of the crystal structures highlighted the importance of shape complementarity for ligand recognition and selectivity. The binding sites identified in the crystal structures of the FAF-DSM and FAW-DSM complexes together with data from the existing literature are used to establish a systematic classification of the ligand-binding sites in the β-lactoglobulin molecule.

Interaction between pH-shifted ovalbumin and insoluble neohesperidin: Experimental and binding mechanism studies

In this study, ovalbumin (OVA) formed a complex with neohesperidin (NH) via a pH-shifting method. The NH-OVA complex self-assembled into NH-OVA nano-particles, which were then characterized and whose binding mechanism was evaluated by using multi-spectroscopic, thermodynamics, and molecular docking simulation methods. Fluorescence intensity decreased after OVA was complexed with NH. The binding constant of the OVA-NH complex was in the order of 6.32 × 10⁵ M^-1 suggesting that the complex is stable. Circular dichroism (CD) analysis showed that α -helix content increased, β-folding, β -turning, and irregular crimp content decreased after OVA and NH binding. Isothermal titration calorimetry results showed that hydrophobic interactions and hydrogen bonds made an important impact in the complex formation. The molecular docking results revealed that Van der Waals forces and hydrogen bonds contributed to the free binding energy of the complex. There were multiple possible surface binding sites between OVA with NH. The obtained results provide new insights into the interaction mechanism of OVA and NH, and as a vehicle for NH, the OVA has shown promising applications in functional foods.

Approach to Study pH-Dependent Protein Association Using Constant-pH Molecular Dynamics: Application to the Dimerization of β-Lactoglobulin

Protein-protein association is often mediated by electrostatic interactions and modulated by pH. However, experimental and computational studies have often overlooked the effect of association on the protonation state of the protein. In this work, we present a methodological approach based on constant-pH molecular dynamics (MD), which aims to provide a detailed description of a pH-dependent protein-protein association, and apply it to the dimerization of β-lactoglobulin (BLG). A selection of analyses is performed using the data generated by constant-pH MD simulations of monomeric and dimeric forms of bovine BLG, in the pH range 3-8. First, we estimate free energies of dimerization using a computationally inexpensive approach based on the Wyman-Tanford linkage theory, calculated in a new way through the use of thermodynamically based splines. The individual free energy contribution of each titratable site is also calculated, allowing for identification of relevant residues. Second, the correlations between the proton occupancies of pairs of sites are calculated (using the Pearson coefficient), and extensive networks of correlated sites are observed at acidic pH values, sometimes involving distant pairs. In general, strongly correlated sites are also slow proton exchangers and contribute significantly to the pH-dependency of the dimerization free energy. Third, we use ionic density as a fingerprint of protein charge distribution and observe electrostatic complementarity between the monomer faces that form the dimer interface, more markedly at the isoionic point (where maximum dimerization occurs) than at other pH values, which might contribute to guide the association. Finally, the pH-dependent dimerization modes are inspected using PCA, among other analyses, and two states are identified: a relaxed state at pH 4-8 (with the typical alignment of the crystallographic structure) and a compact state at pH 3-4 (with a tighter association and rotated alignment). This work shows that an approach based on constant-pH MD simulations can produce rich detailed pictures of pH-dependent protein associations, as illustrated for BLG dimerization.

A 200 nanoseconds all-atom simulation of the pH-dependent EF loop transition in bovine β-lactoglobulin. The role of the orientation of the E89 side chain

In silico molecular dynamics (MD) using crystallographic and NMR data was used to simulate the effects of the protonation state of E89 on the pH-dependent conformational rearrangement of the EF loop, also known as the Tanford transition, in a series of apo-β-lactoglobulin (BLG) structures. Compared to existing studies these simulations were carried out over a much longer time scale (200 ns where the stability of the transition can be evaluated) and used an explicit water model. We considered eight different entries from the Brookhaven Protein Data Bank (PDB) separated into two groups. We observed that fixing the protonation state of E89 prompts the transition of the EF loop only when its side chain is oriented under the loop and into the entrance of the interior cavity. The motion of the EF loop occurs mostly as a step-function and its timing varies greatly from ∼ 20 ns to ∼170 ns from the beginning of the simulation. Once the transition is completed, the protein appears to reach a stable conformation as in a true two-state transition. We also observed novel findings. When the transition occurs, the hydrogen bond between E89 and S116 is replaced with a salt bridge with Lys residues in the βC-CD loop-βD motif. This electrostatic interaction causes the distortion of this motif as well as the protrusion of the GH loop into the aperture of the cavity with the result of limiting the increase of its contour area despite the rotation of the EF loop.Communicated by Ramaswamy H. Sarma.

Simulation of pH-Dependent Conformational Transitions in Membrane Proteins: The CLC-ec1 Cl-/H+ Antiporter

Intracellular transport of chloride by members of the CLC transporter family involves a coupled exchange between a Cl⁻ anion and a proton (H⁺), which makes the transport function dependent on ambient pH. Transport activity peaks at pH 4.5 and stalls at neutral pH. However, a structure of the WT protein at acidic pH is not available, making it difficult to assess the global conformational rearrangements that support a pH-dependent gating mechanism. To enable modeling of the CLC-ec1 dimer at acidic pH, we have applied molecular dynamics simulations (MD) featuring a new force field modification scheme—termed an Equilibrium constant pH approach (ECpH). The ECpH method utilizes linear interpolation between the force field parameters of protonated and deprotonated states of titratable residues to achieve a representation of pH-dependence in a narrow range of physiological pH values. Simulations of the CLC-ec1 dimer at neutral and acidic pH comparing ECpH-MD to canonical MD, in which the pH-dependent protonation is represented by a binary scheme, substantiates the better agreement of the conformational changes and the final model with experimental data from NMR, cross-link and AFM studies, and reveals structural elements that support the gate-opening at pH 4.5, including the key glutamates Glu_in and Glu_ex.

Protein structure prediction using residue-resolved protection factors from hydrogen-deuterium exchange NMR

Hydrogen-deuterium exchange (HDX) measured by nuclear magnetic resonance (NMR) provides structural information for proteins relating to solvent accessibility and flexibility. While this structural information is beneficial, the data cannot be used exclusively to elucidate structures. However, the structural information provided by the HDX-NMR data can be supplemented by computational methods. In previous work, we developed an algorithm in Rosetta to predict structures using qualitative HDX-NMR data (categories of exchange rate). Here we expand on the effort, and utilize quantitative protection factors (PFs) from HDX-NMR for structure prediction. From observed correlations between PFs and solvent accessibility/flexibility measures, we present a scoring function to quantify the agreement with HDX data. Using a benchmark set of 10 proteins, an average improvement of 5.13 Å in root-mean-square deviation (RMSD) is observed for cases of inaccurate Rosetta predictions. Ultimately, seven out of 10 predictions are accurate without including HDX data, and nine out of 10 are accurate when using our PF-based HDX score.

Effect of Protein Genotypes on Physicochemical Properties and Protein Functionality of Bovine Milk: A Review

Milk protein comprises caseins (CNs) and whey proteins, each of which has different genetic variants. Several studies have reported the frequencies of these genetic variants and the effects of variants on milk physicochemical properties and functionality. For example, the C variant and the BC haplotype of αS1-casein (αS1-CN), β-casein (β-CN) B and A1 variants, and κ-casein (κ-CN) B variant, are favourable for rennet coagulation, as well as the B variant of β-lactoglobulin (β-lg). κ-CN is reported to be the only protein influencing acid gel formation, with the AA variant contributing to a firmer acid curd. For heat stability, κ-CN B variant improves the heat resistance of milk at natural pH, and the order of heat stability between phenotypes is BB > AB > AA. The A2 variant of β-CN is more efficient in emulsion formation, but the emulsion stability is lower than the A1 and B variants. Foaming properties of milk with β-lg variant B are better than A, but the differences between β-CN A1 and A2 variants are controversial. Genetic variants of milk proteins also influence milk yield, composition, quality and processability; thus, study of such relationships offers guidance for the selection of targeted genetic variants.

Conformational Changes of Whey and Pea Proteins upon Emulsification Approached by Front-Surface Fluorescence

Proteins are widely used to stabilize emulsions, and plant proteins have raised increasing interest for this purpose. The interfacial and emulsifying properties of proteins depend largely on their molecular properties. We used fluorescence spectroscopy to characterize the conformation of food proteins from different biological origins (dairy or pea) and transformation processes (commercial or lab-made isolates) in solution and at the oil-water interface. The fourth derivative of fluorescence spectra provided insights in the local environment of tryptophan (Trp) residues and thus in the protein structure. In emulsions, whey proteins adsorbed with their Trp-rich region at the oil-water interface. Proteins in the commercial pea isolate were present as soluble aggregates, and no changes in the local environment of the Trp residues were detected upon emulsification, suggesting that these structures adsorb without conformational changes. The lab-purified pea proteins were less aggregated and a Trp-free region of the vicilin adsorbed at the oil-water interface.

pH-Dependent complexation between β-lactoglobulin and lycopene: Multi-spectroscopy, molecular docking and dynamic simulation study

This study aims to investigate the effect of pH levels (pH 7.0 and pH 8.1) on binding ability of β-lactoglobulin (β-LG) with lycopene (LYC) and elucidate interaction mechanisms using multi-spectroscopy and molecular docking study. β-LG at pH 8.1 showed a stronger binding affinity to lycopene than that at pH 7.0 according to binding constant, binding number, energy transfer efficiency, and surface hydrophobicity. Lycopene bound to protein mainly by van der Waals force in the form of static quenching mode and preferred to interact with β-LG at the top of barrel for both pH levels. Molecular dynamic simulation revealed that β-LG/LYC complex at pH 8.1 was more stable than at pH 7.0. β-LG/LYC complexes formed at pH 8.1 showed significantly higher ABTS radical scavenging activity than samples at pH 7.0 (p < 0.05). Data obtained may contribute valuable information for preparing a whey protein-based delivery system for lycopene.

β-Lactoglobulin and Glycodelin: Two Sides of the Same Coin?

The two lipocalins, β-lactoglobulin (βLg) and glycodelin (Gd), are possibly the most closely related members of the large and widely distributed lipocalin family, yet their functions appear to be substantially different. Indeed, the function of β-lactoglobulin, a major component of ruminant milk, is still unclear although neonatal nutrition is clearly important. On the other hand, glycodelin has several specific functions in reproduction conferred through distinct, tissue specific glycosylation of the polypeptide backbone. It is also associated with some cancer outcomes. The glycodelin gene, PAEP, reflecting one of its names, progestagen-associated endometrial protein, is expressed in many though not all primates, but the name has now also been adopted for the β-lactoglobulin gene (HGNC, www.genenames.org). After a general overview of the two proteins in the context of the lipocalin family, this review considers the properties of each in the light of their physiological functional significance, supplementing earlier reviews to include studies from the past decade. While the biological function of glycodelin is reasonably well defined, that of β-lactoglobulin remains elusive.

Investigation of anomalous charge variant profile reveals discrete pH-dependent conformations and conformation-dependent charge states within the CDR3 loop of a therapeutic mAb

During the development of a therapeutic monoclonal antibody (mAb-1), the charge variant profile obtained by pH-gradient cation exchange chromatography (CEX) contained two main peaks, each of which exhibited a unique intrinsic fluorescence profile and demonstrated inter-convertibility upon reinjection of isolated peak fractions. Domain analysis of mAb-1 by CEX and liquid chromatography-mass spectrometry indicated that the antigen-binding fragment chromatographed as two separate peaks that had identical mass. Surface plasmon resonance binding analysis to antigen demonstrated comparable kinetics/affinity between these fractionated peaks and unfractionated starting material. Subsequent molecular modeling studies revealed that the relatively long and flexible complementarity-determining region 3 (CDR3) loop on the heavy chain could adopt two discrete pH-dependent conformations: an “open” conformation at neutral pH where the HC-CDR3 is largely solvent exposed, and a “closed” conformation at lower pH where the solvent exposure of a neighboring tryptophan in the light chain is reduced and two aspartic acid residues near the ends of the HC-CDR3 loop have atypical pKa values. The pH-dependent equilibrium between “open” and “closed” conformations of the HC-CDR3, and its proposed role in the anomalous charge variant profile of mAb-1, were supported by further CEX and hydrophobic interaction chromatography studies. This work is an example of how pH-dependent conformational changes and conformation-dependent changes to net charge can unexpectedly contribute to perceived instability and require thorough analytical, biophysical, and functional characterization during biopharmaceutical drug product development.

Change in conformational, digestive and immunological characteristics of bovine allergen β-lactoglobulin induced by metal ions in combination with heating

Aggregation of bovine β-lactoglobulin is affected easily by external factors. In this study, effects of metal ions combining with temperature on aggregation of β-lactoglobulin were explored. The conformational characteristics of aggregates were detected by environment scanning electron microscope, CD spectrum and free sulfhydryl group, respectively. Digestive and immunological characteristics were assessed by simulated digestion in vitro and ELISA respectively. The results showed that the morphology of β-lactoglobulin aggregates became more amorphous in Cu²⁺ and Mg²⁺ treated samples and more constricted in Zu²⁺-induced protein. Among them, Cu²⁺ altered the secondary structure of β-lactoglobulin aggregates and free sulfhydryl content most as well as that in gastric digestion. However, all ion-treated groups had similar digestive stability in intestinal digestion. Specially, Ca²⁺ and Mg²⁺ made the antigenicity and potential allergenicity of β-lactoglobulin aggregates decrease, which helps us understand the role of metal ions in immunological characteristics.

Binding Sites for Oligosaccharide Repeats from Lactic Acid Bacteria Exopolysaccharides on Bovine β-Lactoglobulin Identified by NMR Spectroscopy

Lactic acid bacterial exopolysaccharides (EPS) are used in the food industry to improve the stability and rheological properties of fermented dairy products. β-Lactoglobulin (BLG), the dominant whey protein in bovine milk, is well known to bind small molecules such as fatty acids, vitamins, and flavors, and to interact with neutral and anionic polysaccharides used in food and pharmaceuticals. While sparse data are available on the affinity of EPS-milk protein interactions, structural information on BLG-EPS complexes, including the EPS binding sites, is completely lacking. Here, binding sites on BLG variant A (BLGA), for oligosaccharides prepared by mild acid hydrolysis of two EPS produced by Streptococcus thermophilus LY03 and Lactobacillus delbrueckii ssp. bulgaricus CNRZ 1187, respectively, are identified by NMR spectroscopy and supplemented by isothermal titration calorimetry (ITC) and molecular docking of complexes. Evidence of two binding sites (site 1 and site 2) on the surface of BLGA is achieved for both oligosaccharides (LY03-OS and 1187-OS) through NMR chemical shift perturbations, revealing multivalency of BLGA for EPS. The affinities of LY03-OS and 1187-OS for BLGA gave K D values in the mM range obtained by both NMR (pH 2.65) and ITC (pH 4.0). Molecular docking suggested that the BLGA and EPS complexes depend on hydrogen bonds and hydrophobic interactions. The findings provide insights into how BLGA engages structurally different EPS-derived oligosaccharides, which may facilitate the design of BLG-EPS complexation, of relevance for formulation of dairy products and improve understanding of BLGA coacervation.

Energetic and structural effects of the Tanford transition on ligand recognition of bovine β-lactoglobulin

Bovine β-lactoglobulin, an abundant protein in whey, is a promising nanocarrier for peroral administration of drug-like hydrophobic molecules, a process that involves transit through the different acidic conditions of the human digestive tract. Among the several pH-induced conformational rearrangements that this lipocalin undergoes, the Tanford transition is particularly relevant. This transition, which occurs with a midpoint around neutral pH, involves a conformational change of the E-F loop that regulates accessibility to the primary binding site. The effect of this transition on the ligand binding properties of this protein has scarcely been explored. In this study, we carried out an energetic and structural characterization of β-lactoglobulin molecular recognition at pH values above and below the zone in which the Tanford transition occurs. The combined analysis of crystallographic, calorimetric, and molecular dynamics data sheds new light on the interplay between self-association, ligand binding, and the Tanford pre- and post-transition conformational states, revealing novel aspects underlying the molecular recognition mechanism of this enigmatic lipocalin.

Small angle X-ray scattering analysis of ligand-bound forms of tetrameric apolipoprotein-D

Human apolipoprotein-D (apoD) is a glycosylated lipocalin that plays a protective role in Alzheimer's disease due to its antioxidant function. Native apoD from human body fluids forms oligomers, predominantly a stable tetramer. As a lipocalin, apoD binds and transports small hydrophobic molecules such as progesterone, palmitic acid and sphingomyelin. Oligomerisation is a common trait in the lipocalin family and is affected by ligand binding in other lipocalins. The crystal structure of monomeric apoD shows no major changes upon progesterone binding. Here, we used small-angle X-ray scattering (SAXS) to investigate the influence of ligand binding and oxidation on apoD oligomerisation and conformation. As a solution-based technique, SAXS is well suited to detect changes in oligomeric state and conformation in response to ligand binding. Our results show no change in oligomeric state of apoD and no major conformational changes or subunit rearrangements in response to binding of ligands or protein oxidation. This highlights the highly stable structure of the native apoD tetramer under various physiologically relevant experimental conditions.

Characterization, spectroscopic and crystallographic analyses of β-lactoglobulin and docosahexaenoic acid nanocomplexes

In this work, we investigated the interaction of docosahexaenoic acid (DHA) with β-lactoglobulin (β-Lg) using spectroscopic and crystallographic methods. The fluorescence results showed that DHA formed complexes with β-Lg with a binding constant of 4.13 × 10⁴ M^-1. The secondary structure of β-Lg was not significantly (p > 0.05) changed after binding with DHA. Dynamic light scattering showed the particle size of β-Lg-DHA complexes was about 5 nm, the same as that of β-Lg alone. The turbidity of DHA in aqueous solution decreased after binding with β-Lg. The crystallographic results showed that DHA was bound at one site in the calyx of β-Lg and that the aliphatic chain was hidden inside the hydrophobic β-barrel while the carboxyl group was located at the calyx entrance. These findings indicate that β-Lg can act as an effective nanocarrier for DHA.

Bovine β-lactoglobulin peptides as novel carriers for flavonoids extracted with supercritical fluids from yellow onion skins

Our study describes in detail the binding mechanism between the main flavonoids that were extracted from onion skins by supercritical CO₂ and peptides from whey proteins, from the perspective of obtaining multifunctional ingredients, with health-promoting benefits. The supercritical CO₂ extract had 202.31 ± 11.56 mg quercetin equivalents/g DW as the major flavonoid and antioxidant activity of 404.93±1.39 mM Trolox/g DW. The experiments on thermolysin-derived peptides fluorescence quenching by flavonoids extract allowed estimating the binding parameters, in terms of binding constants, and the number of binding sites. The thermodynamic analysis indicated that the main forces involved in complex formation were hydrogen bonds and van der Waals interactions. Molecular docking tests indicated that peptide fluorescence quenching upon gradual addition of onion skin extract might be due to flavonoids binding by Val¹⁵ -Ser²¹ . All 7 to 14 amino acids long peptides appeared to have affinity toward quercetin-3,4'-O-diglucoside and quercetin-4'-O-monoglucoside. The study is important as a potential solution for reuse of valuable resources, underutilized, such as whey peptides and yellow onion skins flavonoids for efficient microencapsulation, as a holistic approach to deliver healthy and nutritious food. PRACTICAL APPLICATION: A growing interest was noticed in the last years in investigating the interactions between proteins and different biologically active compounds, such as to provide knowledge for efficient development of new food, pharmaceutical, and cosmetic products. Recent studies suggest that flavonoid-protein complexes may be designed to improve the functional performance of the flavonoids. The results obtained in this study bring certain benefits in terms of exploiting the bioactive potential of both flavonoids and bioactive peptides, for developing of formulas with improved functional properties.

Structural mechanism of the Tanford transition of bovine β-lactoglobulin through microsecond molecular dynamics simulations

X-ray analysis has provided structural data about a pH-driven conformational change in β-lactoglobulin (BLG) known as the Tanford transition, which occurs at around pH 7 and involves the EF loop, which acts as a lid closing the internal cavity of the protein below pH 7 and opening it above pH 7. NMR studies using wild-type BLG have encountered problems trying to explain the Tanford transition, however, they have provided important insight using a dimeric BLG mutant, revealing that the opening and closure of the EF loop consists of two types of motions in the microsecond and milliseconds timescales. This provides valuable information indicating that the dimeric state is a good model to study the Tanford transition, although the understanding of this structural change is still lacking at the atomic level. We performed microsecond molecular dynamics (MD) simulations starting from different conformations of BLG in the monomeric and dimeric state, with protonated and deprotonated E89, in order to explore the Tanford transition. Our results provide structural information for the transition from the closed to the open conformation in BLG and show it occurs in the dimeric state in the microsecond timescale, in line with the fast motion observed through NMR experiments. In addition, MD simulations coupled to MMGBSA approach indicated that the most populated conformer of BLG in the open state is able to bind ligands with similar affinity to that of BLG at neutral pH obtained through crystallographic experiments.Communicated by Ramaswamy H. Sarma.

On the utility of fluorescence-detection analytical ultracentrifugation in probing biomolecular interactions in complex solutions: a case study in milk

β-Lactoglobulin is the most abundant protein in the whey fraction of ruminant milks, yet is absent in human milk. It has been studied intensively due to its impact on the processing and allergenic properties of ruminant milk products. However, the physiological function of β-lactoglobulin remains unclear. Using the fluorescence-detection system within the analytical ultracentrifuge, we observed an interaction involving fluorescently labelled β-lactoglobulin in its native environment, i.e. cow and goat milk, for the first time. Co-elution experiments support that these β-lactoglobulin interactions occur naturally in milk and provide evidence that the interacting partners are immunoglobulins, while further sedimentation velocity experiments confirm that an interaction occurs between these molecules. The identification of these interactions, made possible through the use of fluorescence-detected analytical ultracentrifugation, provides possible clues to the long debated physiological function of this abundant milk protein.

Beta-lactoglobulin-nutrition allergen and nanotransporter of different nature ligands therapy with therapeutic action

β-lactoglobulin is one of the nutrition allergens present in the milk of many mammals, with the exception of human. This protein belongs to the family of lipocalins, consisting of nine antiparallel β-strands (β-A to β-I) and one α-helix. This structure allows it to serve as a nanotransporter of various nature ligands in a pH dependent manner, which allows us to confidently consider it as a reliable carrier of drugs directly into the intestine, bypassing the destructive acidic environment of the stomach. Based on the latest data, this review describes the currently known methods of reducing the allergenicity of beta-lactoglobulin, as well as the mechanisms and methods of forming complexes of this protein with ligands, which emphasizes its importance and versatility and explains the growing interest in studying its properties in recent decades, and also opens up prospects for its practical application in medicine and pharmaceuticals.

The Interaction of Bovine β-Lactoglobulin with Caffeic Acid: From Binding Mechanisms to Functional Complexes

In this study, the interaction of native and transglutaminase (Tgase) cross-linked β-lactoglobulin (β-LG) with caffeic acid (CA) was examined, aiming to obtain functional composites. Knowledge on the binding affinity and interaction mechanism was provided by performing fluorescence spectroscopy measurements, after heating the native and cross-linked protein at temperatures ranging from 25 to 95 °C. Regardless of the protein aggregation state, a static quenching mechanism of intrinsic fluorescence of β-LG by CA was established. The decrease of the Stern–Volmer constants with the temperature increase indicating the facile dissociation of the weakly bound complexes. The thermodynamic analysis suggested the existence of multiple contact types, such as Van der Waals’ force and hydrogen bonds, between β-LG and CA. Further molecular docking tests indicated the existence of various CA binding sites on the β-LG surface heat-treated at different temperatures. Anyway, regardless of the simulated temperature, the CA-β-LG assemblies appeared to be unstable. Compared to native protein, the CA-β-LG and CA-β-LG_Tgase complexes (ratio 1:1) exhibited significantly higher antioxidant activity and inhibitory effects on α-glucosidase, α-amylase, and pancreatic lipase, enzymes associated with metabolic syndrome. These findings might help the knowledge-based development of novel food ingredients with valuable biological properties.

Locating the binding sites of citric acid and gallic acid on milk β-lactoglobulin

β-Lactoglobulin (β-LG) is a member of lipocalin superfamily of transporters for small hydrophobic and hydrophilic molecules. We located the binding sites of citric acid and gallic acid on β-lactoglobulin (β-LG) in aqueous solution at physiological conditions, using spectroscopic methods, thermodynamic analysis and molecular modeling. Thermodynamic parameters ΔH⁰ -9.5 to -6.9 (kJ mol^-1), ΔS⁰ 23.9 to 13.6 (J mol^-1K^-1) and ΔG⁰ -14.5 to -13.6 (kJ mol^-1) showed that acid binds protein via ionic contacts with gallic acid forming stronger protein conjugates consistent with theoretical modeling. Different amino acids are involved in gallic acid and citric acid complexation. Protein conformation was altered with reduction of β-sheet from 58% (free protein) to 46-43% and a major increase in α-helix from 11% (free protein) to 29-23% and random coil structure in the acid-protein, indicating a partial protein destabilization. Communicated by Ramaswamy H. Sarma.

Influence of ultrasound pretreatment on the subsequent glycation of dietary proteins

The influence of ultrasound treatment on the subsequent glycation process of proteins is controversial. Glycation behaviors of bovine serum albumin (BSA), β-lactoglobulin (β-Lg) and β-casein (β-CN) after ultrasound pretreatment (UP) were compared by both evaluating glycation kinetics and analyzing structural changes of proteins. UP resulted in both unfolding and aggregation behavior in protein samples, which altered the accessibility of the Lys and Arg. Five cycles of UP up-regulated the glycation degree of BSA and β-Lg, possibly due to the unfolding behavior induced by UP, which exposed additional glycation sites. In contrast, 30 cycles of UP induced a dramatic increase (by 97.9 nm) in particle size of BSA, thus burying portions of glycation sites and suppressing the glycation process. Notably, UP had minimal influence on glycation kinetics of β-CN, due to its intrinsic disordered structure. Based on proteomics analysis, the preference of Lys and Arg during glycation was found to be changed by UP in BSA and β-Lg. Four, 3 and 3 unique carboxyethylated lysine residues were identified in glycated BSA after 0, 5 and 30 cycles of UP, respectively. This study suggests that the protein glycation can be affected by UP, depending on the ultrasonication duration and native structure of the protein.

Interactions of Whey Proteins with Metal Ions

Whey proteins tend to interact with metal ions, which have implications in different fields related to human life quality. There are two impacts of such interactions: they can provide opportunities for applications in food and nutraceuticals, but may lead to analytical challenges related to their study and outcomes for food processing, storage, and food interactions. Moreover, interactions of whey proteins with metal ions are complicated, requiring deep understanding, leading to consequences, such as metalloproteins, metallocomplexes, nanoparticles, or aggregates, creating a biologically active system. To understand the phenomena of metal–protein interactions, it is important to develop analytical approaches combined with studies of changes in the biological activity and to analyze the impact of such interactions on different fields. The aim of this review was to discuss chemistry of β-lactoglobulin, α-lactalbumin, and lactotransferrin, their interactions with different metal ions, analytical techniques used to study them and the implications for food and nutraceuticals.

Docking and Molecular Dynamics Predictions of Pesticide Binding to the Calyx of Bovine β-Lactoglobulin

Pesticides are used extensively in agriculture, and their residues in food must be monitored to prevent toxicity. The most abundant protein in cow’s milk, β-lactoglobulin (BLG), shows high affinity for diverse hydrophobic ligands in its central binding pocket, called the calyx. Several of the most frequently used pesticides are hydrophobic. To predict if BLG may be an unintended carrier for pesticides, we tested its ability to bind 555 pesticides and their isomers, for a total of 889 compounds, in a rigid docking screen. We focused on the analysis of 60 unique molecules belonging to the five pesticide classes defined by the World Health Organization, that docked into BLG’s calyx with ΔGs ranging from −8.2 to −12 kcal mol⁻¹, chosen by statistical criteria. These “potential ligands” were further analyzed using molecular dynamic simulations, and the binding energies were explored with Molecular Mechanics/Generalized Born/Surface Area (MMGBSA). Hydrophobic pyrethroid insecticides, like cypermethrin, were found to bind as deeply and tightly into the calyx as BLG’s natural ligand, palmitate; while polar compounds, like paraquat, were expelled. Our results suggest that BLG could be a carrier for pesticides, in particular for pyrethroid insecticides, allowing for their accumulation in cow’s milk beyond their solubility restrictions. This analysis opens possibilities for pesticide biosensor design based on BLG.

Structure-based design approach to rational site-directed mutagenesis of β-lactoglobulin

Recombinant proteins play an important role in medicine and have diverse applications in industrial biotechnology. Lactoglobulin has shown great potential for use in targeted drug delivery and body fluid detoxification because of its ability to bind a variety of molecules. In order to modify the biophysical properties of β-lactoglobulin, a series of single-site mutations were designed using a structure-based approach. A 3-dimensional structure alignment of homologous molecules led to the design of nine β-lactoglobulin variants with mutations introduced in the binding pocket region. Seven stable and correctly folded variants (L39Y, I56F, L58F, V92F, V92Y, F105L, M107L) were thoroughly characterized by fluorescence, circular dichroism, isothermal titration calorimetry, size-exclusion chromatography, and X-ray structural investigations. The effects of the amino acid substitutions were observed as slight rearrangements of the binding pocket geometry, but they also significantly influenced the global properties of the protein. Most of the mutations increased the thermal/chemical stability without altering the dimerization constant or pH-dependent conformational behavior. The crystal structures reveal that the I56F and F105L mutations reduced the depth of the binding pocket, which is advantageous since it can reduce the affinity to endogenous fatty acids. The F105L mutant created a unique binding mode for a fatty acid, supporting the idea that lactoglobulin can be altered to bind unique molecules. Selected variants possessing a unique combination of their individual properties can be used for further, more advanced mutagenesis, and the presented results support further research using β-lactoglobulin as a therapeutic delivery agent or a blood detoxifying molecule.