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

Distinct antithrombin activation modes for fondaparinux and natural heparins detected using millisecond hydrogen deuterium exchange and collision induced unfolding

Low molecular weight heparin and synthetic mimetics such as fondaparinux show different binding kinetics, protease specificity, and clinical effects. A combination of allosteric and template-mediated bridging mechanisms have been proposed to explain the differences in rate acceleration and specificity. The difficulty in working with heterogeneous heparin species has rendered a crystallographic interpretation of the differences in antithrombin activation between mimetics and natural heparin inaccessible. In this study, we examine the allosteric changes in antithrombin caused by binding fondaparinux, enoxaparin and depolymerized natural heparins using millisecond hydrogen deuterium exchange mass spectrometry (TRESI-HDX MS) and relate these conformational changes to complex stability in the gas phase using collision induced unfolding (CIU). This exploration reveals that in addition to the dynamic changes caused by fondaparinux, long chain heparins reduce structural flexibility proximal to Arg393, the cleavable residue in the reactive centre loop of the protein. These local changes in protein dynamics are associated with an increase in overall complex stability that increases with heparin chain length. Ultimately, these results shed light on the molecular mechanisms underlying differences in activity and specificity between heparin mimetics and natural heparins.

The Interaction of Factor Xa and IXa with Non-Activated Antithrombin in Michaelis Complex: Insights from Enhanced-Sampling Molecular Dynamics Simulations

The interaction between coagulation factors Xa and IXa and the activated state of their inhibitor, antithrombin (AT),have been investigated using X-ray diffraction studies. However, only mutagenesis data are available for non-activated AT. Our aim was to propose a model based on docking and advanced-sampling molecular dynamics simulations that can reveal the conformational behavior of the systems when AT is not binding a pentasaccharide. We built the initial structure for non-activated AT-FXa and AT-FIXa complexes using HADDOCK 2.4. The conformational behavior was studied using Gaussian accelerated molecular dynamics simulations. In addition to the docked complexes, two systems based on the X-ray structures were also simulated, with and without the ligand. The simulations revealed large variability in conformation for both factors. In the docking-based complex of AT-FIXa, conformations with stable Arg150-AT interactions can exist for longer time periods but the system also has a higher tendency for reaching states with very limited interaction with the "exosite" of AT. By comparing simulations with or without the pentasaccharide, we were able to gain insights into the effects of conformational activation on the Michaelis complexes. RMSF analysis and correlation calculations for the alpha-carbon atoms revealed important details of the allosteric mechanisms. Our simulations provide atomistic models for better understanding the conformational activation mechanism of AT against its target factors.

The mechanism of high affinity pentasaccharide binding to antithrombin, insights from Gaussian accelerated molecular dynamics simulations

The activity of antithrombin (AT), a serpin protease inhibitor, is enhanced by heparin and heparin analogs against its target proteases, mainly thrombin, factors Xa and IXa. Considerable amount of information is available on the multistep mechanism of the heparin pentasaccharide binding and conformational activation. However, much of the details were inferred from 'static' structures obtained by X-ray diffraction. Moreover, limited information is available for the early steps of binding mechanism other than kinetic studies with various ligands. To gain insights into these processes, we performed enhanced sampling molecular dynamics (MD) simulations using the Gaussian Accelerated Molecular Dynamics (GAMD) method, applied previously in drug binding studies. We were able to observe the binding of the pentasaccharide idraparinux to a 'non-activated' AT conformation in two separate trajectories with low root mean square deviation (RMSD) values compared to X-ray structures of the bound state. These trajectories along with further simulations of the AT-pentasaccharide complex provided insights into the mechanisms of multiple conformational transitions, including the expulsion of the hinge region, the extension of helix D and the conformational behavior of the reactive center loop (RCL). We could also confirm the high stability of helix P in non-activated AT conformations, such states might play an important role in heparin binding. 'Generalized correlation' matrices revealed possible paths of allosteric signal propagation to the binding sites for the target proteases, factors Xa and IXa. Enhanced MD simulations of ligand binding to AT may assist the design of new anticoagulant drugs.Communicated by Ramaswamy H. Sarma.

Singular value decomposition for the correlation of atomic fluctuations with arbitrary angle

Many proteins exhibit a critical property called allostery, which enables intra-molecular transmission of information between distal sites. Microscopically, allosteric response is closely related to correlated atomic fluctuations. Conventional correlation analysis correlates the atomic fluctuations at two sites by taking the dot product (DP) between the fluctuations, which accounts only for the parallel and antiparallel components. Here, we present a singular value decomposition (SVD) method that analyzes the correlation coefficient of fluctuation dynamics with an arbitrary angle between the correlated directions. In a model allosteric system, the second PDZ domain (PDZ2) in the human PTP1E protein, approximately one third of the strong correlations have near-perpendicular directions, which are underestimated in the conventional method. The discrimination becomes more prominent for residue pairs with larger separation. The results of the proposed SVD method are more consistent with the experimentally determined PDZ2 dynamics than those of conventional method. In addition, the SVD method improved the prediction accuracy of the allosteric sites in a dataset of 23 known allosteric monomer proteins. The proposed method may inspire extended investigation not only into allostery, but also into protein dynamics and drug design.

Conformational preludes to the latency transition in PAI-1 as determined by atomistic computer simulations and hydrogen/deuterium-exchange mass spectrometry

Both function and dysfunction of serine protease inhibitors (serpins) involve massive conformational change in their tertiary structure but the dynamics facilitating these events remain poorly understood. We have studied the dynamic preludes to conformational change in the serpin plasminogen activator inhibitor 1 (PAI-1). We report the first multi-microsecond atomistic molecular dynamics simulations of PAI-1 and compare the data with experimental hydrogen/deuterium-exchange data (HDXMS). The simulations reveal notable conformational flexibility of helices D, E and F and major fluctuations are observed in the W86-loop which occasionally leads to progressive detachment of β-strand 2 A from β-strand 3 A. An interesting correlation between C_α-RMSD values from simulations and experimental HDXMS data is observed. Helices D, E and F are known to be important for the overall stability of active PAI-1 as ligand binding in this region can accelerate or decelerate the conformational inactivation. Plasticity in this region may thus be mechanistically linked to the conformational change, possibly through facilitation of further unfolding of the hydrophobic core, as previously reported. This study provides a promising example of how computer simulations can help tether out mechanisms of serpin function and dysfunction at a spatial and temporal resolution that is far beyond the reach of any experiment.

Heparan sulfate and heparin interactions with proteins

Heparan sulfate (HS) polysaccharides are ubiquitous components of the cell surface and extracellular matrix of all multicellular animals, whereas heparin is present within mast cells and can be viewed as a more sulfated, tissue-specific, HS variant. HS and heparin regulate biological processes through interactions with a large repertoire of proteins. Owing to these interactions and diverse effects observed during in vitro, ex vivo and in vivo experiments, manifold biological/pharmacological activities have been attributed to them. The properties that have been thought to bestow protein binding and biological activity upon HS and heparin vary from high levels of sequence specificity to a dependence on charge. In contrast to these opposing opinions, we will argue that the evidence supports both a level of redundancy and a degree of selectivity in the structure-activity relationship. The relationship between this apparent redundancy, the multi-dentate nature of heparin and HS polysaccharide chains, their involvement in protein networks and the multiple binding sites on proteins, each possessing different properties, will also be considered. Finally, the role of cations in modulating HS/heparin activity will be reviewed and some of the implications for structure-activity relationships and regulation will be discussed.