Molecular Dynamics Simulation of Protein Biosurfactants

Insights into the emulsification mechanism of the surfactant-like protein oleosin

Oleosins are proteins with a unique central hydrophobic hairpin designed to stabilize lipid droplets (oleosomes) in plant seeds. For efficient droplet stabilization, the hydrophobic hairpin with a strong affinity for the apolar droplet core is flanked by hydrophilic arms on each side. This gives oleosins a unique surfactant-like shape making them a very interesting protein. In this study, we tested if isolated oleosins retain their ability to stabilize oil-in-water emulsions, and investigated the underlying stabilization mechanism. Due to their surfactant-like shape, oleosins when dispersed in aqueous buffers associated to micelle-like nanoparticles with a size of ∼33 nm. These micelles, in turn, clustered into larger aggregates of up to 20 µm. Micelle aggregation was more extensive when oleosins lacked charge. During emulsification, oleosin micelles and micelle aggregates dissociated and mostly individual oleosins adsorbed on the oil droplet interface. Oleosins prevented the coalescence of the oil droplets and if sufficiently charged, droplet flocculation as well.

Interdisciplinary Overview of Lipopeptide and Protein-Containing Biosurfactants

Biosurfactants are amphipathic molecules capable of lowering interfacial and superficial tensions. Produced by living organisms, these compounds act the same as chemical surfactants but with a series of improvements, the most notable being biodegradability. Biosurfactants have a wide diversity of categories. Within these, lipopeptides are some of the more abundant and widely known. Protein-containing biosurfactants are much less studied and could be an interesting and valuable alternative. The harsh temperature, pH, and salinity conditions that target organisms can sustain need to be understood for better implementation. Here, we will explore biotechnological applications via lipopeptide and protein-containing biosurfactants. Also, we discuss their natural role and the organisms that produce them, taking a glimpse into the possibilities of research via meta-omics and machine learning.

Nanoemulsions: Formulation, characterization, biological fate, and potential role against COVID-19 and other viral outbreaks

Viral diseases are emerging as global threats. Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), that causes coronavirus disease (COVID-19), has severe global impacts. Safety, dosage, and potency of vaccines recently approved for emergency use against SARS-CoV-2 need further evaluation. There is still no effective treatment against COVID-19; therefore, safe, and effective vaccines or therapeutics against SARS-CoV-2 are urgently needed. Oil-in-water nanoemulsions (O/W NEs) are emerging as sophisticated, protective, and therapeutic platforms. Encapsulation capacity, which offers better drug pharmacokinetics, coupled with the tunable surfaces present NEs as promising tools for pharmaceutical applications. The challenges facing drug discovery, and the advancements of NEs in drug delivery demonstrate the potential of NEs against evolving diseases, like COVID-19. Here we summarize current COVID-19 knowledge and discuss the composition, stability, preparation, characterization, and biological fate of O/W NEs. We also provide insights into NE structural-functional properties that may contribute to therapeutic or preventative solutions against COVID-19.

Molecular Dynamics Simulations of Nanostructures Formed by Hydrophobins and Oil in Seawater

Classical molecular dynamics simulations using the Martini coarse-grained force field were performed to study oil nanodroplets surrounded by fungal hydrophobin (HP) proteins in seawater. The class I EAS and the class II HFBII HPs were studied along with two model oils, namely, benzene and n-decane. Both HPs exhibit free energy minima at the oil-seawater interface, which is deeper in benzene compared to the n-decane systems. Larger constraint forces are required to keep both HPs within the n-decane phase compared to inside benzene, with HFBII being more affine to benzene compared to EAS. Smaller surface tensions are observed at benzene-seawater interfaces coated with HPs compared to their n-decane counterparts. In the latter the surface tension remains unchanged upon increases in the concentration of HPs, whereas in benzene systems adding more HPs lead to decreases in surface tension. EAS has a larger tendency to cluster together in the interface compared to HFBII, with both HPs having larger coordination numbers when surrounding benzene droplets compared to when they are around n-decane nanoblobs. The HP-oil nanostructures in seawater examined have radii of gyration ranging between 2 and 12 nm, where the n-decane structures are larger and have more irregular shapes compared to the benzene systems. The n-decane molecules within the nanostructures form a compact spherical core, with the HPs partially covering its surface and clustering together, conferring irregular shapes to the nanostructures. The EAS with n-decane structures are larger and have more irregular shapes compared to their HFBII counterparts. In contrast, in the HP-benzene structures both HPs tend to penetrate the oil part of the droplet. The HFBII-benzene structures having the larger oil/HP ratios examined tend to be more compact and spherical compared to their EAS counterparts; however, some of the HFBII-benzene systems that have smaller oil/HP ratios have a more elongated structure compared to their EAS counterparts. This simulation study provides insights into HP-oil nanostructures that are smaller than the oil droplets and gas bubbles recently studied in experiments and, thus, might be challenging to examine with experimental techniques.

Combining Experimental Isotherms, Minimalistic Simulations, and a Model to Understand and Predict Chemical Adsorption onto Montmorillonite Clays

An attractive approach to minimize human and animal exposures to toxic environmental contaminants is the use of safe and effective sorbent materials to sequester them. Montmorillonite clays have been shown to tightly bind diverse toxic chemicals. Due to their promise as sorbents to mitigate chemical exposures, it is important to understand their function and rapidly screen and predict optimal clay-chemical combinations for further testing. We derived adsorption free-energy values for a structurally and physicochemically diverse set of toxic chemicals using experimental adsorption isotherms performed in the current and previous studies. We studied the diverse set of chemicals using minimalistic MD simulations and showed that their interaction energies with calcium montmorillonite clays calculated using simulation snapshots in combination with their net charge and their corresponding solvent's dielectric constant can be used as inputs to a minimalistic model to predict adsorption free energies in agreement with experiments. Additionally, experiments and computations were used to reveal structural and physicochemical properties associated with chemicals that can be adsorbed to calcium montmorillonite clay. These properties include positively charged groups, phosphine groups, halide-rich moieties, hydrogen bond donor/acceptors, and large, rigid structures. The combined experimental and computational approaches used in this study highlight the importance and potential applicability of analogous methods to study and design novel advanced sorbent systems in the future, broadening their applicability for environmental contaminants.

Orientation and Conformation of Proteins at the Air-Water Interface Determined from Integrative Molecular Dynamics Simulations and Sum Frequency Generation Spectroscopy

Understanding the assembly of proteins at the air-water interface (AWI) informs the formation of protein films, emulsion properties, and protein aggregation. Determination of protein conformation and orientation at an interface is difficult to resolve with a single experimental or simulation technique alone. To date, the interfacial structure of even one of the most widely studied proteins, lysozyme, at the AWI remains unresolved. In this study, molecular dynamics (MD) simulations are used to determine if the protein adopts a side-on, head-on, or axial orientation at the AWI with two different forcefields, GROMOS-53a6 + SPC/E and a99SB-disp + TIP4P-D. Vibrational sum frequency generation (SFG) spectroscopy experiments and spectral SFG calculations validate consistency between the structure determined from MD and experiments. Overall, we show with strong agreement that lysozyme adopts an axial conformation at pH 7. Further, we provide molecular-level insight as to how pH influences the binding domains of lysozyme resulting in side-on adsorption near the isoelectric point of the lysozyme.

Dynamic heterogeneity in complex interfaces of soft interface-dominated materials

Complex interfaces stabilized by proteins, polymers or nanoparticles, have a much richer dynamics than those stabilized by simple surfactants. By subjecting fluid-fluid interfaces to step extension-compression deformations, we show that in general these complex interfaces have dynamic heterogeneity in their relaxation response that is well described by a Kohlrausch-Williams-Watts function, with stretch exponent β between 0.4-0.6 for extension, and 0.6-1.0 for compression. The difference in β between expansion and compression points to an asymmetry in the dynamics. Using atomic force microscopy and simulations we prove that the dynamic heterogeneity is intimately related to interfacial structural heterogeneity and show that the dominant mode for stretched exponential relaxation is momentum transfer between bulk and interface, a mechanism which has so far largely been ignored in experimental surface rheology. We describe how its rate constant can be determined using molecular dynamics simulations. These interfaces clearly behave like disordered viscoelastic solids and need to be described substantially different from the 2d homogeneous viscoelastic fluids typically formed by simple surfactants.