This study aims to investigate the interactions appearing when the beta-2-glycoprotein-1 binds to a lipid bilayer. The inter- and intra-molecular forces acting between the two macromolecular systems have been investigated using a molecular dynamics simulation method. The importance of water bridges has also been addressed. Additionally, the viscoelastic response of the bilayer has been studied. In detail, the (saturated-chain) 1,2-dipalmitoyl-sn-glycero-3-phosphocholine (DPPC) and (unsaturated-chain) 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphoethanolamine (POPE) bilayers have been chosen to test their behavior near the protein. Both of the lipids have a polar head but different chemical structures and are similar to the main phospholipids present in the synovial fluid. This study is meaningful for further explaining the worsening friction properties in articular cartilage, as the inactivation of phospholipid bilayers by beta-2-glycoprotein-1 is believed to be a cause of the destruction of cartilage in most rheumatic diseases and osteoarthritis. It was found that the protein binds stronger to the DPPC bilayer than to the POPE, but in both cases, it has the potential to change the local bilayer stability. Nevertheless, the binding forces are placed within a small area (only a few lipids contribute to the binding, creating many interactions). However, together, they are not stronger than the covalent bonds between C-O, thus, potentially, it is possible to push the lipids into the bilayer but detaching the lipids' heads from the tail is not possible. Additionally, the protein causes water displacement from the vicinity of the bilayer, and this may be a contributor to the instability of the bilayer (disrupting the water bridges needed for the stabilization of the bilayer, especially in the case of DPPC where the heads are not so well stabilized by H-bonds as they are in POPE). Moreover, it was found that the diffusivity of lipids in the DPPC bilayer bound to the protein is significantly different from the diffusivity of the ones which are not in contact with the protein. The POPE bilayer is stiffer due to intramolecular interactions, which are stronger than in the DPPC; thus, the viscous to elastic effects in the POPE case are more significant than in the case of the DPPC. It is, therefore, harder to destabilize the POPE bilayer than the DPPC one.
Keywords: DPPC; POPE; beta-2-glycoprotein-1; hydrogen bonds; hydrophobic interactions; ionic interactions; molecular dynamics.