Molecular dynamics simulations of lysozyme-lipid systems: probing the early steps of protein aggregation

Phospholipid-induced secondary structural changes of lysozyme polymorphic amyloid fibrils studied using vacuum-ultraviolet circular dichroism

The hallmark of amyloidosis, such as Alzheimer's disease and Parkinson's disease, is the deposition of amyloid fibrils in various internal organs. The onset of the disease is related to the strength of cytotoxicity caused by toxic amyloid species. Furthermore, amyloid fibrils show polymorphism, where some types of fibrils are cytotoxic while others are not. It is thus essential to understand the molecular mechanism of cytotoxicity, part of which is caused by the interaction between amyloid polymorphic fibrils and cell membranes. Here, using amyloid polymorphs of hen egg white lysozyme, which is associated with hereditary systemic amyloidosis, showing different levels of cytotoxicity and liposomes of DMPC and DMPG, changes in the secondary structure of the polymorphs and the structural state of phospholipid membranes caused by the interaction were investigated using vacuum-ultraviolet circular dichroism (VUVCD) and Laurdan fluorescence measurements, respectively. Analysis has shown that the more cytotoxic polymorph increases the antiparallel β-sheet content and causes more disorder in the membrane structure while the other less cytotoxic polymorph shows the opposite structural changes and causes less structural disorder in the membrane. These results suggest a close correlation between the structural properties of amyloid fibrils and the degree of structural disorder of phospholipid membranes, both of which are involved in the fundamental process leading to amyloid cytotoxicity.

Elucidating the Molecular Interactions between Lipids and Lysozyme: Evaporation Resistance and Bacterial Barriers for Dry Eye Disease

Dry eye disease (DED) is a chronic condition characterized by ocular dryness and inflammation. The tear film lipid layer (TFLL) is the outermost layer composed of lipids and proteins that protect the ocular surface. However, environmental contaminants can disrupt its structure, potentially leading to DED. Although the importance of tear proteins in the TFLL functionality has been clinically recognized, the molecular mechanisms underlying TFLL-protein interactions remain unclear. In this study, we investigated tear protein-lipid interactions and analyzed their role in the TFLL functionality. The results show that lysozyme (LYZ) increases the stability of the TFLL by reducing its surface tension and increasing its surface pressure, resulting in increased TFLL evaporation and bacterial invasion resistance, with improved wettability and lubrication performance. These findings highlight the critical role of LYZ in maintaining ocular health and provide potential avenues for investigating novel approaches to DED treatment and patient well-being.

Aggregation of Lysozyme in the Presence of a Mixed Bilayer of POPC and POPG

Understanding the molecular mechanisms by which amyloidogenic proteins interact with membranes is a challenging task. Amyloid accumulates from many human diseases have been observed to contain membrane lipids. In this work, coarse-grained molecular dynamics simulations have been used to inspect hen egg white lysozyme (HEWL) aggregation and membrane association in the presence of a pure POPC (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) bilayer and a POPC and POPG (1-palmitoyl-2-oleoyl-sn-glycero-3-phosphatidylglycerol) mixed bilayer. It was observed that, in both cases, two HEWLs formed aggregates. In the presence of a mixed bilayer, after aggregation, the aggregated system started to interact with the membrane. It has been found that one of the lysozymes which came closer to the mixed bilayer unfolded more. The process of the initial insertion of an aggregated system in the mixed bilayer has been analyzed. The structural rearrangements of the protein and lipids were analyzed as well along the course of the simulation. Although with a pure POPC bilayer, aggregation was observed, the aggregated system moved away from the membrane. We believe that our study will provide considerable insights into lysozyme aggregation in the presence of a membrane environment.

Molecular Dynamics Simulation to Uncover the Mechanisms of Protein Instability During Freezing

Freezing is a common process applied in the pharmaceutical industry to store and transport biotherapeutics. Herewith, multi-scale molecular dynamics simulations of Lactate dehydrogenase (LDH) protein in phosphate buffer with/without ice formation performed to uncover the still poorly understood mechanisms and molecular details of protein destabilization upon freezing. Both fast and slow ice growing conditions were simulated at 243 K from one or two-side of the simulation box, respectively. The rate of ice formation at all-atom simulations was crucial to LDH stability, as faster freezing rates resulted in enhanced structural stability maintained by a higher number of intramolecular hydrogen bonds, less flexible protein's residues, lower solvent accessibility and greater structural compactness. Further, protein aggregation investigated by coarse-grained simulations was verified to be initiated by extended protein structures and retained by electrostatic interactions of the salt bridges between charged residues and hydrogen bonds between polar residues of the protein. Lastly, the study of free energy of dissociation through steered molecular dynamics simulation revealed LDH was destabilized by the solvation of the hydrophobic core and the loss of hydrophobic interactions. For the first time, experimentally validated molecular simulations revealed the detailed mechanisms of LDH destabilization upon ice formation and cryoconcentration of solutes.

Study of Structural stability and formation mechanisms in DSPC and DPSM liposomes: A coarse-grained molecular dynamics simulation

Liposomes or biological vesicles can be created from cholesterol, phospholipid, and water. Their stability is affected by their phospholipid composition which can influence disease treatment and drug delivery efficacy. In this study, the effect of phospholipid type on the formation and stability of liposomes using coarse-grained molecular dynamics simulations is investigated. For this purpose, the simulation study of the DSPC (1,2-distearoyl-sn-glycero-3-phosphocholine) and DPSM (Egg sphingomyelin) lipids were considered. All simulations were carried out using the Gromacs software and Martini force field 2.2. Energy minimization (3000 steps) model, equilibrium at constant volume to adjust the temperature at 400 Kelvin and equilibrium at constant pressure to adjust the pressure, at atmospheric pressure (1 bar) have been validated. Microsecond simulations, as well as formation analysis including density, radial distribution function, and solvent accessible surface area, demonstrated spherical nanodisc structures for the DPSM and DSPC liposomes. The results revealed that due to the cylindrical geometric structure and small-size head group, the DSPC lipid maintained its perfectly spherical structure. However, the DPSM lipid showed a conical geometric structure with larger head group than other lipids, which allows the liposome to form a micelle structure. Although the DSPC and DPSM lipids used in the laboratory tests exhibit liposome and micelle behaviors, the simulation results revealed their nanodisc structures. Energy analysis including overall energy, Van der Waals interaction energy, and electrostatic interaction energy showed that DPSM liposome is more stable than DSPC liposome.

Computational investigation of a ternary model of SnoN-SMAD3-SMAD4 complex

The transforming growth factor β (TGFβ) plays an essential role in the regulation of cellular processes such as cell proliferation, migration, differentiation, and apoptosis by association with SMAD transcriptional factors that are regulated by the transcriptional regulator SnoN. The crystal structure of SnoN-SMAD4 reveals that SnoN can adopt two binding modes, the open and closed forms, at the interfaces of SMAD4 subunits. Accumulating evidence indicates that SnoN can interact with both SMAD3 and SMAD4 to form a ternary SnoN-SMAD3-SMAD4 complex in the TGFβ signaling pathway. However, how the interaction of SnoN with the SMAD3 and SMAD4 remains unclear. Here, molecular dynamics (MD) simulations and molecular modeling methods were performed to figure out this issue. The simulations reveal that SnoN^open exists in two, open and semi-closed, conformations. Molecular modeling and MD simulation studies suggest that the SnoN^closed form interferes with the SMAD3-SMAD4 protein; in contract, the SnoN^open can form a stable SnoN-SMAD3-SMAD4 complex. These mechanistic mechanisms may help elucidate the detailed engagement of SnoN with two SMAD3 and SMAD4 transcriptional factors in the regulation of TGFβ signaling pathway.

Characterization of heterogeneous intermediate ensembles on the guanidinium chloride-induced unfolding pathway of β-lactoglobulin

Folding pathway of β-LgA (β-lactoglobulin) evolves through the conformational α→β transition. The α→β transition is a molecular hallmark of various neurodegenerative diseases. Thus, β-LgA may serve as a good model for understanding molecular mechanism of protein aggregation involved in neurodegenerative diseases. Here, we studied the conformational dynamics of β-LgA in 6 M GdmCl at different temperatures using MD simulations. Structural order parameters such as RMSD, R_g, SASA, native contacts (Q), hydrophobic distal-matrix and free-energy landscape (FEL) were used to investigate the conformational transitions. Our results show that GdmCl destabilizes secondary and tertiary structure of β-LgA by weakening the hydrophobic interactions and hydrogen bond network. Multidimensional FEL shows the presence of different unfolding intermediates at 400 K. I1 is long-lived intermediate which has mostly intact native secondary structure, but loose tertiary structure. I2 is structurally compact intermediate formed after the partial loss of secondary structure. The transiently and infrequently buried evolution of W19 shows that intermediate conformational ensembles are structurally heterogeneous. We observed that the intermediate conformations are largely stabilized by non-native H-bonds. The outcome of this work provides the molecular details of intermediates trapped due to non-native interactions that may be regarded as pathogenic conformations involved in neurodegenerative diseases.Communicated by Ramaswamy H. Sarma.

Diffusional interaction behavior of NSAIDs in lipid bilayer membrane using molecular dynamics (MD) simulation: Aspirin and Ibuprofen

In this research, for the first time, molecular dynamics (MD) method was used to simulate aspirin and ibuprofen at various concentrations and in neutral and charged states. Effects of the concentration (dosage), charge state, and existence of an integral protein in the membrane on the diffusion rate of drug molecules into lipid bilayer membrane were investigated on 11 systems, for which the parameters indicating diffusion rate and those affecting the rate were evaluated. Considering the diffusion rate, a suitable score was assigned to each system, based on which, analysis of variance (ANOVA) was performed. By calculating the effect size of the indicative parameters and total scores, an optimum system with the highest diffusion rate was determined. Consequently, diffusion rate controlling parameters were obtained: the drug-water hydrogen bond in protein-free systems and protein-drug hydrogen bond in the systems containing protein.