Uncovering Membrane-Bound Models of Coagulation Factors by Combined Experimental and Computational Approaches

Optimizing Coarse-Grained Models for Large-Scale Membrane Protein Simulation

Coarse-grained (CG) models have been developed for studying membrane proteins at physiologically relevant scales. Such methods, including popular CG lipid models, exhibit stability and efficiency at moderate scales, but they can become impractical or even unusable beyond a critical size due to various technical issues. Here, we report that these scale-dependent issues can arise from progressively slower relaxation dynamics and become confounded by unforeseen instabilities observed only at larger scales. To address these issues, we systemically optimized a 4-site solvent-free CG lipid model that is suitable for conducting micron-scale molecular dynamics simulations of membrane proteins under various membrane properties. We applied this lipid model to explore the long-range membrane deformation induced by a large mechanosensitive ion channel, PIEZO. We show that the optimized CG models are powerful in elucidating the structural and dynamic interplay between PIEZO and the membrane. Furthermore, we anticipate that our methodological insights can prove useful for resolving issues stemming from scale-dependent limitations of similar CG methodologies.

Membrane binding and lipid-protein interaction of the C2 domain from coagulation factor V

Anchoring of coagulation factors to anionic regions of the membrane involves the C2 domain as a key player. The rate of enzymatic reactions of the coagulation factors is increased by several orders of magnitude upon membrane binding. However, the precise mechanisms behind the rate acceleration remain unclear, primarily because of a lack of understanding of the conformational dynamics of the C2-containing factors and corresponding complexes. We elucidate the membrane-bound form of the C2 domain from human coagulation factor V (FV-C2) by characterizing its membrane binding the specific lipid-protein interactions. Employing all-atom molecular dynamics simulations and leveraging the highly mobile membrane-mimetic (HMMM) model, we observed spontaneous binding of FV-C2 to a phosphatidylserine (PS)-containing membrane within 2-25 ns across twelve independent simulations. FV-C2 interacted with the membrane through three loops (spikes 1-3), achieving a converged, stable orientation. Multiple HMMM trajectories of the spontaneous membrane binding provided extensive sampling and ample data to examine the membrane-induced effects on the conformational dynamics of C2 as well as specific lipid-protein interactions. Despite existing crystal structures representing presumed "open" and "closed" states of FV-C2, our results revealed a continuous distribution of structures between these states, with the most populated structures differing from both "open" and "closed" states observed in crystal environments. Lastly, we characterized a putative PS-specific binding site formed by K23, Q48, and S78 located in the groove enclosed by spikes 1-3 (PS-specificity pocket), suggesting a different orientation of a bound headgroup moiety compared to previous proposals based upon analysis of static crystal structures.

Factor VIII moiety of recombinant Factor VIII Fc fusion protein impacts Fc effector function and CD16+ NK cell activation

Recombinant Factor VIII-Fc fusion protein (rFVIIIFc) is an enhanced half-life therapeutic protein product used for the management of hemophilia A. Recent studies have demonstrated that rFVIIIFc interacts with Fc gamma receptors (FcγR) resulting in the activation or inhibition of various FcγR-expressing immune cells. We previously demonstrated that rFVIIIFc, unlike recombinant Factor IX-Fc (rFIXFc), activates natural killer (NK) cells via Fc-mediated interactions with FcγRIIIA (CD16). Additionally, we showed that rFVIIIFc activated CD16+ NK cells to lyse a FVIII-specific B cell clone. Here, we used human NK cell lines and primary NK cells enriched from peripheral blood leukocytes to study the role of the FVIII moiety in rFVIIIFc-mediated NK cell activation. Following overnight incubation of NK cells with rFVIIIFc, cellular activation was assessed by measuring secretion of the inflammatory cytokine IFNγ by ELISA or by cellular degranulation. We show that anti-FVIII, anti-Fc, and anti-CD16 all inhibited indicating that these molecules were involved in rFVIIIFc-mediated NK cell activation. To define which domains of FVIII were involved, we used antibodies that are FVIII domain-specific and demonstrated that blocking FVIII C1 or C2 domain-mediated membrane binding potently inhibited rFVIIIFc-mediated CD16+ NK cell activation, while targeting the FVIII heavy chain domains did not. We also show that rFVIIIFc binds CD16 with about five-fold higher affinity than rFIXFc. Based on our results we propose that FVIII light chain-mediated membrane binding results in tethering of the fusion protein to the cell surface, and this, together with increased binding affinity for CD16, allows for Fc-CD16 interactions to proceed, resulting in NK cellular activation. Our working model may explain our previous results where we observed that rFVIIIFc activated NK cells via CD16, whereas rFIXFc did not despite having identical IgG1 Fc domains.

Dynamic Nature of Staphylococcus aureus Type I Signal Peptidases

Molecular dynamics simulations are used to interrogate the dynamic nature of Staphylococcus aureus Type I signal peptidases, SpsA and SpsB, including the impact of the P29S mutation of SpsB. Fluctuations and plasticity- rigidity characteristics vary among the proteins, particularly in the extracellular domain. Intriguingly, the P29S mutation, which influences susceptibility to arylomycin antibiotics, affect the mechanically coupled motions in SpsB. The integrity of the active site is crucial for catalytic competency, and variations in sampled structural conformations among the proteins are consistent with diverse peptidase capabilities. We also explored the intricate interactions between the proteins and the model S. aureus membrane. It was observed that certain membrane-inserted residues in the loop around residue 50 (50s) and C-terminal loops, beyond the transmembrane domain, give rise to direct interactions with lipids in the bilayer membrane. Our findings are discussed in the context of functional knowledge about these signal peptidases, offering additional understanding of dynamic aspects relevant to some cellular processes with potential implications for drug targeting strategies.

Cholesterol and M2 Rendezvous in Budding and Scission of Influenza A Virus

The cholesterol of the host cell plasma membrane and viral M2 protein plays a crucial role in multiple stages of infection and replication of the influenza A virus. Cholesterol is required for the formation of heterogeneous membrane microdomains (or rafts) in the budozone of the host cell that serves as assembly sites for the viral components. The raft microstructures act as scaffolds for several proteins. Cholesterol may further contribute to the mechanical forces necessary for membrane scission in the last stage of budding and help to maintain the stability of the virus envelope. The M2 protein has been shown to cause membrane scission in model systems by promoting the formation of curved lipid bilayer structures that, in turn, can lead to membrane vesicles budding off or scission intermediates. Membrane remodeling by M2 is intimately linked with cholesterol as it affects local lipid composition, fluidity, and stability of the membrane. Thus, both cholesterol and M2 protein contribute to the efficient and proper release of newly formed influenza viruses from the virus-infected cells.

A systematic approach for evaluating the role of surface-exposed loops in trypsin-like serine proteases applied to the 170 loop in coagulation factor VIIa

Proteases play a major role in many vital physiological processes. Trypsin-like serine proteases (TLPs), in particular, are paramount in proteolytic cascade systems such as blood coagulation and complement activation. The structural topology of TLPs is highly conserved, with the trypsin fold comprising two β-barrels connected by a number of variable surface-exposed loops that provide a surprising capacity for functional diversity and substrate specificity. To expand our understanding of the roles these loops play in substrate and co-factor interactions, we employ a systematic methodology akin to the natural truncations and insertions observed through evolution of TLPs. The approach explores a larger deletion space than classical random or directed mutagenesis. Using FVIIa as a model system, deletions of 1–7 amino acids through the surface exposed 170 loop, a vital allosteric regulator, was introduced. All variants were extensively evaluated by established functional assays and computational loop modelling with Rosetta. The approach revealed detailed structural and functional insights recapitulation and expanding on the main findings in relation to 170 loop functions elucidated over several decades using more cumbersome crystallization and single deletion/mutation methodologies. The larger deletion space was key in capturing the most active variant, which unexpectedly had a six-amino acid truncation. This variant would have remained undiscovered if only 2–3 deletions were considered, supporting the usefulness of the methodology in general protease engineering approaches. Our findings shed further light on the complex role that surface-exposed loops play in TLP function and supports the important role of loop length in the regulation and fine-tunning of enzymatic function throughout evolution.

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.