Study of bovine β-casein at water/lipid interface by molecular modeling

The role of hydrophobic interactions in the molten globule state of globular protein modulated by surfactants

In order to highlight the role of hydrophobic interactions in the molten globule (MG) state of globular protein modulated by surfactants, the interactions of bovine α-lactalbumin (α-LA) with alkyl trimethylammonium bromides (CnTAB, n = 10, 12, 14, and 16) have been studied by experimental and theoretical techniques. Isothermal titration calorimetry (ITC) showed that the enthalpy changes (ΔH) and area of the enthalpogram increased with increasing the chain length of CnTAB. The result of fluorescence, circular dichroism (CD) and 1H nuclear magnetic resonance (NMR) spectrum suggested that C10TAB and C12TAB unfolded α-LA partially, C14TAB reconstructed protein with a native-like secondary structure content, and C16TAB induced an MG state α-LA. The SAXS results confirmed that the tertiary structure of α-LA was disrupted by C16TAB forming an MG state complex with a micelle-like structure even at the surfactants concentrations below CMC. As indicated by MD results, the β-domain and unstructured region(s) were involved in the MG state α-LA modulated by CnTAB. This work not only provides molecular insights into the role of hydrophobic interactions in the MG state of a globular protein but also helps understand the mechanism of preparing α-LA based biomacromolecule modulated by hydrophobic interactions.

A theoretical and experimental investigation of the effect of sodium dodecyl sulfate on the structural and conformational properties of bovine β-casein

A predicted three-dimensional structure of bovine β-casein was constructed using homology modeling with the aid of MODELLER and I-TASSER programs, with the validity and reliability of the models evaluated according to stereochemical qualities and small angle X-ray scattering. By comparing the results obtained from the two models using the CRYSOL program, an optimal model of the β-casein structure derived from I-TASSER was selected and used in subsequent molecular dynamics (MD) analysis. 300 ns MD simulations of β-casein in water and in the presence of different SDS concentrations at 300 K were performed. The results of the MD simulations indicated that SDS molecules played a dual role in modifying the conformation of β-casein at 300 K. Concentrations of SDS below its CMC (1 mM), at which only the monomer form of SDS was present, induced β-casein to lose its secondary structure by converting helices into random coils; however the conformation of the complex was still comparable with that of native β-casein. In the presence of 10 mM SDS (above its CMC), the helical content of β-casein was increased along with reduced random coils, and the structural rearrangement led to a more compact conformation. The latter change is likely related to the hydrophobic interactions that dominate the binding of the C-terminal region, along with the anchoring of sulfate groups of SDS on the positively charged N-terminal portion via electrostatic attraction. Hydrogen bonding supplemented the SDS-induced stabilization of β-casein. A correlated "necklace and bead" model, in which the micelles nucleate on the protein hydrophobic sites, was proposed for the structure of β-casein-SDS complexes.

The self-assembly, aggregation and phase transitions of food protein systems in one, two and three dimensions

The aggregation of proteins is of fundamental relevance in a number of daily phenomena, as important and diverse as blood coagulation, medical diseases, or cooking an egg in the kitchen. Colloidal food systems, in particular, are examples that have great significance for protein aggregation, not only for their importance and implications, which touches on everyday life, but also because they allow the limits of the colloidal science analogy to be tested in a much broader window of conditions, such as pH, ionic strength, concentration and temperature. Thus, studying the aggregation and self-assembly of proteins in foods challenges our understanding of these complex systems from both the molecular and statistical physics perspectives. Last but not least, food offers a unique playground to study the aggregation of proteins in three, two and one dimensions, that is to say, in the bulk, at air/water and oil/water interfaces and in protein fibrillation phenomena. In this review we will tackle this very ambitious task in order to discuss the current understanding of protein aggregation in the framework of foods, which is possibly one of the broadest contexts, yet is of tremendous daily relevance.

Low-resolution 1H spin-spin relaxation of n-decane/water emulsions stabilized by beta-casein

A low-resolution 1H NMR relaxometry study on the dynamics of an n-decane/water emulsion stabilized by beta-casein is presented. Spin-spin (transverse) relaxation time constants (T2) were used to assess relative mobilities of emulsion components, by a selective deuteration procedure. Data analysis allowed the emulsion investigated to be described by a heterogeneous collection of dynamically distinct populations. A major population of n-decane molecules presented an average mobility that very nearly approached that of pure solvent, which is compatible with its occurrence in the emulsion continuous microphase. beta-Casein molecules displayed a prevalent population with significantly decreased mobility as compared to the free protein in solution, which is in accordance with the protein location at the oil/water interface. Also, a major H2O population with significantly lower average T2 as compared to the pure liquid was detected and has been assigned to interfacial water.

The effect of ethanol and nitric oxide on the N-nitrosodimethylamine formation in HepG2 cells overexpressing CYP2E1

The influence of lipopolysaccharide (LPS) and the nitric oxide synthase (iNOS) inhibitor--N-nitro-L-arginine methyl ester (L-NAME)--on the formation of N-nitrosodimethylamine (NDMA) by HepG2 cells, engineered to overexpress CYP2E1, was assessed and compared with data from empty vector-transfected cells. HepG2 cells produced significant amounts of NDMA but its levels in the culture media of cells overexpressing CYP2E1 was significantly lower than in empty-vector transfected cells. LPS increased the formation of NDMA, the expression of the iNOS and the production of the nitric oxide (NO). On the other hand, L-NAME significantly decreased NDMA levels. The results above indicate that the synthesis of NDMA by HepG2 cells depends on NO production. Furthermore, ethanol did not affect iNOS expression but decreased NDMA levels in CYP2E1-transfected cells below the detection limit. It is probably caused by the increased N-nitrosodimethylamine metabolism. In conclusion, HepG2 cells' ability to synthesize NO with simultaneous CYP2E1 activation may lead to an increase of carcinogenic products of the NDMA metabolism.

Brownian dynamics simulation of adsorbed layers of interacting particles subjected to large extensional deformation

We present Brownian dynamics simulations of the compression and expansion of monolayers adsorbed at a planar interface. The surface-active species are modelled as monodisperse spherical particles that can form particle-particle elastic bonds. The objective is to model the large compression and expansion of viscoelastic protein films investigated in Langmuir trough experiments. We determine the stress-strain response of the system and the associated microstructural changes induced by the large deformation of the interface as a function of particle adsorption energy, and bond breakability and stiffness. We also study the effect of the velocity of compression and the type of compression (uniaxial or homogeneous) on the mechanism of collapse of the adsorbed films. Furthermore, we present simulations on complex mixed systems containing both bond-forming particles (modelling protein) and nonbond-forming particles (modelling surfactant). We find that the preferential desorption of one type of particle or the other, upon compression, is sensitive to the extent of bond breakability of the bond-forming species.

Model system for testing the efficacy of antioxidants in muscle foods

The objective of this research was to study the effect of the antioxidants, delta-tocopherol, butylated hydroxyanisole (BHA), butylated hydroxytoluene (BHT), tertiary butylhydroquinone (TBHQ), and propyl gallate in a model system of lean muscle and canola oil and to compare the effects with those in minced herring. Two carrier solvents with different dielectric constants (epsilon), ethanol (epsilon = 24) and oil (epsilon= 2), were used. Oxidation was measured using thiobarbituric acid reactive substances (TBARS) and sensory analysis. In both the lean muscle-canola oil model system and in herring muscle, the hydrophilic antioxidants, propyl gallate and TBHQ, were more effective in providing oxidative stability than the lipophilic antioxidants, delta-tocopherol and BHT. The oxidative stability of a cod muscle-canola oil system in the presence of propyl gallate, and delta-tocopherol was not affected by the dielectric constant of the carrier solvent, while BHA was more effective as an antioxidant when added in the polar solvent ethanol.

Competitive adsorption of proteins and low-molecular-weight surfactants: computer simulation and microscopic imaging

Proteins and low-molecular-weight (LMW) surfactants are used in the food industry as emulsifying (and foaming) ingredients and as stabilizers. These attributes are related to their ability to adsorb at fluid-fluid (and gas-fluid) interfaces lowering the interfacial (and surface) tension of liquids. Hence, the study of the properties of adsorbed layers of these molecules can be expected to lead to a better understanding of their effect on food products. Direct proof of the validity of mesoscopic models of systems of proteins and LMW surfactants can only be achieved by quantitative theoretical predictions being tested against both macroscopic and mesoscopic experiments. Computer simulation constitutes one of the few available tools to predict mathematically the behaviour of models of realistic complexity. Furthermore, experimental techniques such as atomic force microscopy (AFM) now allow high resolution imaging of these systems, providing the mesoscopic scale measurements to compare with the simulations. In this review, we bring together a number of related findings that have been generated at this mesoscopic level over the past few years. A useful simple model consisting of spherical particles interacting via bonded and unbonded forces is described, and the derived computer simulation results are compared against those from the imaging experiments. Special attention is paid to the adsorption of binary mixtures of proteins, mixtures of LMW surfactants, and also protein+surfactant mixed systems. We believe that further development of these mathematically well-defined physical models is necessary in order to achieve a proper understanding of the key physico-chemical processes involved.

Dependence of the Interfacial Behavior of β-Casein on Phosphoserine Residues

The role of the phosphoserine residues on the dynamical and structural properties of beta-casein was studied by molecular dynamics of the protein in water/lipid interfacial regions. The initial protein structure adopted in the modeling was that proposed for bovine beta-casein A2, where the five phosphoserine residues, originally present in its primary structure, were partially or totally substituted by serine residues. The simulations revealed a dependence of the interfacial behavior of beta-casein on the phosphorylation grade. When only partially dephosphorylated, the protein showed a similar behavior as that observed for the original beta-casein reported in previous work. During dynamics, the protein migrated from the aqueous environment towards the lipid medium, and remained attached to the interface separating both media. Quite different was the dynamics of the totally dephosphorylated beta-casein, that did not perceive the interface and immersed incessantly into lipid medium. The results suggest that the phosphoserine residues appear to be, in fact, intrinsically related to the mechanisms of beta-casein emulsion stabilization.

Nonspecific interaction forces at water-membrane interface by forced molecular dynamics simulations

Nonspecific interactions are the main driving forces for the behavior of molecules with great affinity for biologic membranes. To investigate not only the molecular details of these interactions but to estimate their magnitude as well, the theoretical method of Forced Molecular Dynamics Simulations, based on the Atomic Force Spectroscopy experimental technique, was applied. In this approach, an additional one-dimensional elastic force, representing the cantilever probe, was incorporated to the force field of a Molecular Dynamics computational program. This force represents a spring fixed on one end to a selected atom of the molecule; the other end of the spring is displaced at constant velocity to pull the molecule out of the membrane. The force experimented by the molecule due to the spring, is proportional to the spring elongation relative to its equilibrium position. This value is registered during the entire simulation, and its maximum value will determine the molecule-membrane interaction force. Nonexplicit medium simulations were carried out. Polar and apolar media were considered according to their polarizability degree and a specific dielectric constant value was assigned. In this approach, the membrane was considered as the apolar region limited by two flat surfaces with a polar aqueous medium. The potential energy discontinuity at the interfaces was smoothed by considering the polarization-induced effects using the image method. The results of this methodology are presented using a small system, a single Alanine amino acid model, which enables extended simulations in a microsecond time scale. The confinement of this amino acid at the interface reduces its degrees of freedom and forces it to adopt one of the six defined conformations. A correlation between these stable structures at the water-membrane interface and the interaction force value was determined.

Protein micelles from lipoxygenase 3

Heat-induced conformational changes in lipoxygenase 3 were characterized by differential scanning calorimetry. The positions of the observed transitions were sensitive to the composition of the buffer. In particular, lipoxygenase 3 heated in carbonate buffer at pH 8.0 formed large soluble aggregates. Variable-temperature circular dichroism revealed that the formation of the aggregates was not accompanied by the unfolding of the C-terminal domain, which is composed primarily of alpha-helix. The aggregates were investigated using size exclusion chromatography, native polyacrylamide gel electrophoresis, dynamic light scattering, and electron microscopy. The data were consistent with the formation of roughly spherical particles with an average hydrodynamic radius of 26 nm and an approximate composite molecular weight of 10,000,000 Da. To account for the formation of soluble aggregates from lipoxygenase 3, we propose that hydrophobic amino acid residues are exposed by unfolding of the N-terminal beta-barrel domain of the protein resulting in the formation of protein micelles with a hydrophilic surface composed of the C-terminal domains.

Molecular modeling and dynamics of the sodium channel inactivation gate

The intracellular linker L(III-IV) of voltage-gated sodium channels is known to be involved in their mechanism of inactivation. Its primary sequence is well conserved in sodium channels from different tissues and species. However, the role of charged residues in this region, first thought to play an important role in inactivation, has not been well identified, whereas the IFM triad (I1488-M1490) has been characterized as the crucial element for inactivation. In this work, we constructed theoretical models and performed molecular dynamics simulations, exploring the role of L(III-IV)-charged residues in the presence of a polar/nonpolar planar interface represented by a dielectric discontinuity. From structural predictions, two alpha-helical segments are proposed. Moreover, from dynamics simulations, a time-conserved motif is detected and shown to play a relevant role in guiding the inactivation particle toward its receptor site.

Folding interpenetration in a gliadin model: the role of the characteristic octapeptide motif

A model of a rheologically relevant protein, omega-gliadin, is proposed and studied in this work by means of molecular dynamics techniques. The model is based on an octapeptide repeat motif that is experimentally described as characteristic of that protein and as constituting it almost entirely. The initial molecular structure consisted of 20 such repeats. It was optimized and the dynamics developed along 980 ps, at dielectric constant epsilon = 80. Remarkable structural features were observed for the model built, such as an elongated twisted tubular overall structure with a peculiar interpenetrating folding pattern, of a very regular character, organized strand formation, topologically segregated sites on the outer surface with an alternate hydrophilic/hydrophobic character and a hydrophilic inner cavity. Dynamics produced significantly more relaxed structures, but was not able to change the main geometric features presented by the original structure. Preliminary attempts of correlating some structural/dynamic aspects observed for the model with features of gliadin rheological behavior are presented.