Water-Intake and Water-Molecule Paths to the Active Site of Secretory Phospholipase A2 Studied Using MD Simulations and the Tracking Tool AQUA-DUCT

Revisiting the reaction pathways for phospholipid hydrolysis catalyzed by phospholipase A2 with QM/MM methods

Secreted phospholipase A2 (sPLA2) is a Ca2+-dependent, widely distributed enzyme superfamily in almost all mammalian tissues and bacteria. It is also a critical component of the venom of nearly all snakes, as well as many invertebrate species. In non-venomous contexts, sPLA2 assumes significance in cellular signaling pathways by binding cell membranes and catalyzing ester bond hydrolysis at the sn-2 position of phospholipids. Elevated levels of GIIA sPLA2 have been detected in the synovial fluid of arthritis patients, where it exhibits a pro-inflammatory function. Consequently, identifying sPLA2 inhibitors holds promise for creating highly effective pharmaceutical treatments. Beyond arthritis, the similarities among GIIA sPLA2s offer an opportunity for developing treatments against snakebite envenoming, the deadliest neglected tropical disease. Despite decades of study, the details of PLA2 membrane-binding, substrate-binding, and reaction mechanism remain elusive, demanding a comprehensive understanding of the sPLA2 catalytic mechanism. This study explores two reaction mechanism hypotheses, involving one or two water molecules, and distinct roles for the Ca2+ cofactor. Our research focuses on the human synovial sPLA2 enzyme bound to lipid bilayers of varying phospholipid compositions, and employing adiabatic QM/MM and QM/MM MD umbrella sampling methods to energetically and geometrically characterize the structures found along both reaction pathways. Our studies demonstrate the higher frequency of productive conformations within the single-water pathway, also revealing a lower free energy barrier for hydrolyzing POPC. Furthermore, we observe that the TS of this concerted one-step reaction closely resembles transition state geometries observed in X-ray crystallography complexes featuring high-affinity transition state analogue inhibitors.

Candida antarctica lipase B performance in organic solvent at varying water activities studied by molecular dynamics simulations

Applications of lipases in low-water environments are found across a broad range of industries, including the pharmaceutical and oleochemical sectors. This includes condensation reactions in organic solvents where the enzyme activity has been found to depend strongly on both the solvent and the water activity (aw). Despite several experimental and computational studies, knowledge is largely empirical, and a general predictive approach is much needed. To close this gap, we chose native Candida antarctica lipase B (CALB) and two mutants thereof and used molecular dynamics (MD) simulations to gain a molecular understanding of the effect of aw on the specific activity of CALB in hexane. Based on the simulations, we propose four criteria to understand the performance of CALB in organic media, which is supported by enzyme kinetics experiments. First, the lipase must be stable in the organic solvent, which was the case for native CALB and the two mutants studied here. Secondly, water clusters that form and grow close to the active site must not block the path of substrate molecules into the active site. Thirdly, the lipase's lid must not cover the active site. Finally, mutations and changes in aw must not disrupt the geometry of the active site. We show that mutating specific residues close to the active site can hinder water cluster formation and growth, making the lipase resistant to changes in aw. Our computational screening criteria could potentially be used to screen in-silico designed variants, so only promising candidates could be pushed forward to characterisation.

Conformational Plasticity-Rigidity Axis of the Coagulation Factor VII Zymogen Elucidated by Atomistic Simulations of the N-Terminally Truncated Factor VIIa Protease Domain

The vast majority of coagulation factor VII (FVII), a trypsin-like protease, circulates as the inactive zymogen. Activated FVII (FVIIa) is formed upon proteolytic activation of FVII, where it remains in a zymogen-like state and it is fully activated only when bound to tissue factor (TF). The catalytic domains of trypsin-like proteases adopt strikingly similar structures in their fully active forms. However, the dynamics and structures of the available corresponding zymogens reveal remarkable conformational plasticity of the protease domain prior to activation in many cases. Exactly how ligands and cofactors modulate the conformational dynamics and function of these proteases is not entirely understood. Here, we employ atomistic simulations of FVIIa (and variants hereof, including a TF-independent variant and N-terminally truncated variants) to provide fundamental insights with atomistic resolution into the plasticity-rigidity interplay of the protease domain conformations that appears to govern the functional response to proteolytic and allosteric activation. We argue that these findings are relevant to the FVII zymogen, whose structure has remained elusive despite substantial efforts. Our results shed light on the nature of FVII and demonstrate how conformational dynamics has played a crucial role in the evolutionary adaptation of regulatory mechanisms that were not present in the ancestral trypsin. Exploiting this knowledge could lead to engineering of protease variants for use as next-generation hemostatic therapeutics.