Molecular dynamics simulation of a phospholipase A2-substrate complex

Characterization of a novel acidic phospholipase A2 isolated from the venom of Bothrops mattogrossensis: From purification to structural modeling

Phospholipases A2 (PLA2s) are highly prevalent in Bothrops snake venom and play a crucial role in inflammatory responses and immune cell activation during envenomation. Despite their significance, the specific role of PLA2s from Bothrops mattogrossensis venom (BmV) in inflammation is not fully understood. This study sought to isolate and characterize a novel acidic PLA2 from BmV, designated BmPLA2-A, and to evaluate its effects on human umbilical vein endothelial cells (HUVECs), with a specific focus on cytotoxicity, adhesion, and detachment. BmPLA2-A was isolated through a multi-step chromatographic procedure, involving cation exchange (CM-Sepharose), hydrophobic interaction (n-butyl-Sepharose-HP), and reversed-phase (C-18) chromatography. SDS-PAGE analysis revealed a single protein band of approximately 15 kDa. The primary structure of BmPLA2-A was determined by LC-MS/MS, while its tertiary structure was modeled using AlphaFold. Enzymatic activity was verified with the synthetic substrate 4N3OBA. Molecular dynamics simulations were conducted to further investigate the catalytic mechanism of BmPLA2-A at the molecular level. In vitro assays on HUVECs revealed that BmPLA2-A neither induce cytokine release (IL-6, IL-8, IL-1β, TNF) nor affected cell viability, adhesion, or detachment. The characteristics of BmPLA2-A are consistent with those of acidic Asp-49 PLA2 enzymes, highlighting its potential involvement in the cytotoxic and inflammatory effects of the venom.

The interaction of phospholipase A2 with a phospholipid bilayer: coarse-grained molecular dynamics simulations

A number of membrane-active enzymes act in a complex environment formed by the interface between a lipid bilayer and bulk water. Although x-ray diffraction studies yield structures of isolated enzyme molecules, a detailed characterization of their interactions with the interface requires a measure of how deeply such a membrane-associated protein penetrates into a lipid bilayer. Here, we apply coarse-grained (CG) molecular dynamics (MD) simulations to probe the interaction of porcine pancreatic phospholipase A2 (PLA2) with a lipid bilayer containing palmitoyl-oleoyl-phosphatidyl choline and palmitoyl-oleoyl-phosphatidyl glycerol molecules. We also used a configuration from a CG-MD trajectory to initiate two atomistic (AT) MD simulations. The results of the CG and AT simulations are evaluated by comparison with available experimental data. The membrane-binding surface of PLA2 consists of a patch of hydrophobic residues surrounded by polar and basic residues. We show this proposed footprint interacts preferentially with the anionic headgroups of the palmitoyl-oleoyl-phosphatidyl glycerol molecules. Thus, both electrostatic and hydrophobic interactions determine the location of PLA2 relative to the bilayer. From a general perspective, this study demonstrates that CG-MD simulations may be used to reveal the orientation and location of a membrane-surface-bound protein relative to a lipid bilayer, which may subsequently be refined by AT-MD simulations to probe more detailed interactions.

Modeling of enzyme-substrate complexes for the metalloproteases MMP-3, ADAM-9 and ADAM-10

The matrix metalloproteases (MMPs) and the ADAMs (A Disintegrin And Metalloprotease domain) are proteolytic enzyme families containing a catalytic zinc ion, that are implicated in a variety of normal and pathological processes involving tissue remodeling and cancer. Synthetic MMP inhibitors have been designed for applications in pathological situations. However, a greater understanding of substrate binding and the catalytic mechanism is required so that more effective and selective inhibitors may be developed for both experimental and clinical purposes. By modeling a natural substrate spanning P4-P4' in complex with the catalytic domains, we aim to compare substrate-specificities between Stromelysin-1 (MMP-3), ADAM-9 and ADAM-10, with the aid of molecular dynamics simulations. Our results show that the substrate retains a favourable antiparallel beta-sheet conformation on the P-side in addition to the well-known orientation of the P'-region of the scissile bond, and that the primary substrate selectivity is dominated by the sidechains in the S1' pocket and the S2/S3 region. ADAM-9 has a hydrophobic residue as the central determinant in the S1' pocket, while ADAM-10 has an amphiphilic residue, which suggests a different primary specificity. The S2/S3 pocket is largely hydrophobic in all three enzymes. Inspired by our molecular dynamics calculations and supported by a large body of literature, we propose a novel, hypothetical, catalytic mechanism where the Zn-ion polarizes the oxygens from the catalytic glutamate to form a nucleophile, leading to a tetrahedral oxyanion anhydride transition state.

A molecular dynamics investigation of mono and dimeric states of the outer membrane enzyme OMPLA

OMPLA is a phospholipase found in the outer membranes of many Gram-negative bacteria. Enzyme activation requires calcium-induced dimerisation plus bilayer perturbation. As the conformation of OMPLA in the different crystal forms (monomer versus dimer; with/without bound Ca(2+)) is remarkably similar we have used multi-nanosecond molecular dynamics (MD) simulations to probe possible differences in conformational dynamics that may be related to enzyme activation. Simulations of calcium-free monomeric OMPLA, of the Ca(2+)-bound dimer, and of the Ca(2+)-bound dimer with a substrate analogue covalently linked to the active site serine have been performed, all with the protein embedded in a phospholipid (POPC) bilayer. All simulations were stable, but differences in the dynamic behaviour of the protein between the various states were observed. In particular, the stability of the active site and the hydrophobic substrate-binding cleft varied. Dimeric OMPLA is less flexible than monomeric OMPLA, especially around the active site. In the absence of bound substrate analogue, the hydrophobic substrate-binding cleft of dimeric OMPLA collapses. A model is proposed whereby the increased stability of the active site in dimeric OMPLA is a consequence of the local ordering of water around the nearby calcium ion. The observed collapse of the substrate-binding cleft may explain the experimentally observed occurrence of multiple dimer conformations of OMPLA, one of which is fully active while the other shows significantly reduced activity.

Changes in a phospholipid bilayer induced by the hydrolysis of a phospholipase A2 enzyme: a molecular dynamics simulation study

Phospholipase A2 (PLA2) enzymes are important in numerous physiological processes. Their function at lipid-water interfaces is also used as a biophysical model for protein-membrane interactions. These enzymes catalyze the hydrolysis of the sn-2 bonds of various phospholipids and the hydrolysis products are known to increase the activity of the enzymes. Here, we have applied molecular dynamics (MD) simulations to study the membrane properties in three compositionally different systems that relate to PLA2 enzyme action. One-nanosecond simulations were performed for a 1-palmitoyl-2-linoleoyl-sn-glycero-3-phosphatidylcholine (PLPC) bilayer and for two of its PLA2-hydrolyzed versions, i.e., bilayers consisting of lysophospholipids and of either free charged linoleate or free uncharged linoleic acid molecules. The results revealed loosening of the structure in the hydrolyzed bilayer due to increased mobility of the molecules in the direction normal to the bilayer. This loss of integrity due to the hydrolysis products is in accord with observations that not only the presence of hydrolysis products, but also a variety of other perturbations of the membrane may activate PLA2. Additionally, changes were observed in other structural parameters and in the electrostatic potential across the membrane-water interface. These changes are discussed in relation to the simulation methodology and the experimental observations of PLA2-hydrolyzed membranes.

Catalytic significance of the specificity of divalent cations as KS* and kcat* cofactors for secreted phospholipase A2

Calcium is required for the substrate binding and for the chemical step of the interfacial catalytic turnover cycle of pancreatic phospholipase A2 (PLA2), but not for the binding of the enzyme to the interface. The role of calcium and other divalent cations (C) is analyzed for the effect on the substrate binding and kcat* for the chemical step. The cofactor role of 3d-cations(II) (C) for the hydrolysis of dimyristoylphosphatidylmethanol (DMPM) vesicles is characterized as an equilibrium dissociation constant for the interfacial binary (E*C) and ternary (E*CL) complexes of PLA2 and substrate mimics (L). Of the cations(II) that promote the binding of a mimic to the enzyme at the interface (E*), only a subgroup supports the chemical step. For example, Cd, Zn, and Cu form ternary E*CL complexes with kcat* of <1 s-1, compared to the rate of >100 s-1 with Ca, Fe, Mn, Co, and Ni. Oxygen exchange from H218O to the products of hydrolysis of DMPM incorporates one 18O in myristate. Incorporation of the first and second 18O occurs during the incubation of both the products of hydrolysis in H218O with PLA2 and Ca, but not with Zn. The cation-dependent changes in the UV difference spectrum, associated with the formation of E*C and E*CL, suggest that the changes are mainly due to catalytic His-48, and possibly Tyr-52 and Tyr-73, and are different with Ca as opposed to Zn. These results and simulations suggest considerable plasticity in the calcium binding and catalytic site environment. It is proposed that the higher ground state stability of the E*CS complex with the inhibitory cations increases the effective activation energy. For the chemical step, calcium coordinated with a nucleophilic water and the ester carbonyl oxygen facilitates the near-attack geometry in the E*CaS, and the His-48.Asp-99 pair acts as a proton acceptor. As a prelude to establishing the catalytic mechanism, factors controlling the energetically demanding transition state are also discussed.

X-ray crystal structure and molecular dynamics simulation of bovine pancreas phospholipase A2-n-dodecylphosphorylcholine complex

The crystal structure of n-dodecylphosphorylcholine (n-C12PC)-bovine pancreas phospholipase A2 (PLA2) complex provided the following structural characteristics: (1) the dodecyl chain of n-C12PC was located at the PLA2 N-terminal helical region by hydrophobic interactions, which corresponds to the binding pocket of 2-acyl fatty acid chain (beta-chain) of the substrate phospholipid, (2) the region from Lys-53 to Lys-56 creates a choline-receiving pocket of n-C12PC and (3) the N-terminal group of Ala-1 shifts significantly toward the Tyr-52 OH group by the binding of the n-C12PC inhibitor. Since the accuracy of the X-ray analysis (R = 0.275 at 2.3 A resolution) was insufficient to establish these important X-ray insights, the complex structure was further investigated through the molecular dynamics (MD) simulation, assuming a system in aqueous solution at 310K. The MD simulation covering 176 ps showed that the structural characteristics observed by X-ray analysis are intrinsic and also stable in the dynamic state. Furthermore, the MD simulation made clear that the PLA2 binding pocket is large enough to permit the conformational fluctuation of the n-C12PC hydrocarbon chain.