Computational methods for studying serpin conformational change and structural plasticity

Reactive centre loop dynamics and serpin specificity

Serine proteinase inhibitors (serpins), typically fold to a metastable native state and undergo a major conformational change in order to inhibit target proteases. However, conformational lability of the native serpin fold renders them susceptible to misfolding and aggregation, and underlies misfolding diseases such as α₁-antitrypsin deficiency. Serpin specificity towards its protease target is dictated by its flexible and solvent exposed reactive centre loop (RCL), which forms the initial interaction with the target protease during inhibition. Previous studies have attempted to alter the specificity by mutating the RCL to that of a target serpin, but the rules governing specificity are not understood well enough yet to enable specificity to be engineered at will. In this paper, we use conserpin, a synthetic, thermostable serpin, as a model protein with which to investigate the determinants of serpin specificity by engineering its RCL. Replacing the RCL sequence with that from α1-antitrypsin fails to restore specificity against trypsin or human neutrophil elastase. Structural determination of the RCL-engineered conserpin and molecular dynamics simulations indicate that, although the RCL sequence may partially dictate specificity, local electrostatics and RCL dynamics may dictate the rate of insertion during protease inhibition, and thus whether it behaves as an inhibitor or a substrate. Engineering serpin specificity is therefore substantially more complex than solely manipulating the RCL sequence, and will require a more thorough understanding of how conformational dynamics achieves the delicate balance between stability, folding and function required by the exquisite serpin mechanism of action.

How pyridoxal 5'-phosphate differentially regulates human cytosolic and mitochondrial serine hydroxymethyltransferase oligomeric state

Adaptive metabolic reprogramming gives cancer cells a proliferative advantage. Tumour cells extensively use glycolysis to sustain anabolism and produce serine, which not only refuels the one-carbon units necessary for the synthesis of nucleotide precursors and for DNA methylation, but also affects the cellular redox homeostasis. Given its central role in serine metabolism, serine hydroxymethyltransferase (SHMT), a pyridoxal 5'-phosphate (PLP)-dependent enzyme, is an attractive target for tumour chemotherapy. In humans, the cytosolic isoform (SHMT1) and the mitochondrial isoform (SHMT2) have distinct cellular roles, but high sequence identity and comparable catalytic properties, which may complicate development of successful therapeutic strategies. Here, we investigated how binding of the cofactor PLP controls the oligomeric state of the human isoforms. The fact that eukaryotic SHMTs are tetrameric proteins while bacterial SHMTs function as dimers may suggest that the quaternary assembly in eukaryotes provides an advantage to fine-tune SHMT function and differentially regulate intertwined metabolic fluxes, and may provide a tool to address the specificity problem. We determined the crystal structure of SHMT2, and compared it to the apo-enzyme structure, showing that PLP binding triggers a disorder-to-order transition accompanied by a large rigid-body movement of the two cofactor-binding domains. Moreover, we demonstrated that SHMT1 exists in solution as a tetramer, both in the absence and presence of PLP, while SHMT2 undergoes a dimer-to-tetramer transition upon PLP binding. These findings indicate an unexpected structural difference between the two human SHMT isoforms, which opens new perspectives for understanding their differing behaviours, roles or regulation mechanisms in response to PLP availability in vivo.

Dynamic properties of the native free antithrombin from molecular dynamics simulations: computational evidence for solvent- exposed Arg393 side chain

While antithrombin (AT) has small basal inhibitory activity, it reaches its full inhibitory potential against activated blood coagulation factors, FXa, FIXa, and FIIa (thrombin), via an allosteric and/or template (bridging) mechanism by the action of heparin, heparan sulfate, or heparin-mimetic pentasaccharides (PS). From the numerous X-ray structures available for different conformational states of AT, only indirect and incomplete conclusions can be drawn on the inherently dynamic properties of AT. As a typical example, the basal inhibitory activity of AT cannot be interpreted on the basis of "non-activated" free antithrombin X-ray structures since the Arg393 side chain, playing crucial role in antithrombin-proteinase interaction, is not exposed. In order to reveal the intrinsic dynamic properties and the reason of basal inhibitory activity of antithrombin, 2 μs molecular dynamics simulations were carried out on its native free-forms. It was shown from the simulation trajectories that the reactive center loop which is functioning as "bait" for proteases, even without any biasing potential can populate conformational state in which the Arg393 side chain is solvent exposed. It is revealed from the trajectory analysis that the peptide sequences correspond to the helix D extension, and new helix P formation can be featured with especially large root-mean-square fluctuations. Mutual information analyses of the trajectory showed remarkable (generalized) correlation between those regions of antithrombin which changed their conformations as the consequence of AT-PS complex formation. This suggests that allosteric information propagation pathways are present even in the non-activated native form of AT.

Cofactor-dependent conformational heterogeneity of GAD65 and its role in autoimmunity and neurotransmitter homeostasis

The human neuroendocrine enzyme glutamate decarboxylase (GAD) catalyses the synthesis of the inhibitory neurotransmitter gamma-aminobutyric acid (GABA) using pyridoxal 5'-phosphate as a cofactor. GAD exists as two isoforms named according to their respective molecular weights: GAD65 and GAD67. Although cytosolic GAD67 is typically saturated with the cofactor (holoGAD67) and constitutively active to produce basal levels of GABA, the membrane-associated GAD65 exists mainly as the inactive apo form. GAD65, but not GAD67, is a prevalent autoantigen, with autoantibodies to GAD65 being detected at high frequency in patients with autoimmune (type 1) diabetes and certain other autoimmune disorders. The significance of GAD65 autoinactivation into the apo form for regulation of neurotransmitter levels and autoantibody reactivity is not understood. We have used computational and experimental approaches to decipher the nature of the holo → apo conversion in GAD65 and thus, its mechanism of autoinactivation. Molecular dynamics simulations of GAD65 reveal coupling between the C-terminal domain, catalytic loop, and pyridoxal 5'-phosphate-binding domain that drives structural rearrangement, dimer opening, and autoinactivation, consistent with limited proteolysis fragmentation patterns. Together with small-angle X-ray scattering and fluorescence spectroscopy data, our findings are consistent with apoGAD65 existing as an ensemble of conformations. Antibody-binding kinetics suggest a mechanism of mutually induced conformational changes, implicating the flexibility of apoGAD65 in its autoantigenicity. Although conformational diversity may provide a mechanism for cofactor-controlled regulation of neurotransmitter biosynthesis, it may also come at a cost of insufficient development of immune self-tolerance that favors the production of GAD65 autoantibodies.