Characterization of an Extensive Interface on Vitronectin for Binding to Plasminogen Activator Inhibitor-1: Adoption of Structure in an Intrinsically Disordered Region.

Small-angle neutron scattering (SANS) measurements were pursued to study human vitronectin, a protein found in tissues and the circulation that regulates cell adhesion/migration and proteolytic cascades that govern hemostasis and pericellular proteolysis. Many of these functions occur via interactions with its binding partner, plasminogen activator inhibitor-1 (PAI-1), the chief inhibitor of proteases that lyse and activate plasminogen. We focused on a region of vitronectin that remains uncharacterized from previous X-ray scattering, nuclear magnetic resonance, and computational modeling approaches and which we propose is involved in binding to PAI-1. This region, which bridges the N-terminal somatomedin B (SMB) domain with a large central ß-propeller domain of vitronectin, appears unstructured and has characteristics of an intrinsically disordered domain (IDD). The effect of osmolytes was evaluated using circular dichroism and SANS to explore the potential of the IDD to undergo a disorder-to-order transition. The results suggest that the IDD favors a more ordered structure under osmotic pressure; SANS shows a smaller radius of gyration (Rg) and a more compact fold of the IDD upon addition of osmolytes. To test whether PAI-1 binding is also coupled to folding within the IDD structure, a set of SANS experiments with contrast variation were performed on the complex of PAI-1 with a vitronectin fragment corresponding to the N-terminal 130 amino acids (denoted the SMB-IDD because it contains the SMB domain and IDD in linear sequence). Analysis of the SANS data using the Ensemble Optimization Method confirms that the SMB-IDD adopts a more compact configuration when bound to PAI-1. Calculated structures for the PAI-1:SMB-IDD complex suggest that the IDD provides an interaction surface outside of the primary PAI-1-binding site located within the SMB domain; this binding is proposed to lead to the assembly of higher-order structures of vitronectin and PAI-1 commonly found in tissues.

Resolving distinct molecular origins for copper effects on PAI-1.

Components of the fibrinolytic system are subjected to stringent control to maintain proper hemostasis. Central to this regulation is the serpin plasminogen activator inhibitor-1 (PAI-1), which is responsible for specific and rapid inhibition of fibrinolytic proteases. Active PAI-1 is inherently unstable and readily converts to a latent, inactive form. The binding of vitronectin and other ligands influences stability of active PAI-1. Our laboratory recently observed reciprocal effects on the stability of active PAI-1 in the presence of transition metals, such as copper, depending on the whether vitronectin was also present (Thompson et al. Protein Sci 20:353-365, 2011). To better understand the molecular basis for these copper effects on PAI-1, we have developed a gel-based copper sensitivity assay that can be used to assess the copper concentrations that accelerate the conversion of active PAI-1 to a latent form. The copper sensitivity of wild-type PAI-1 was compared with variants lacking N-terminal histidine residues hypothesized to be involved in copper binding. In these PAI-1 variants, we observed significant differences in copper sensitivity, and these data were corroborated by latency conversion kinetics and thermodynamics of copper binding by isothermal titration calorimetry. These studies identified a copper-binding site involving histidines at positions 2 and 3 that confers a remarkable stabilization of PAI-1 beyond what is observed with vitronectin alone. A second site, independent from the two histidines, binds metal and increases the rate of the latency conversion.

Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations.

Protein arginine methyltransferases (PRMTs) catalyze the transfer of the methyl group from S-adenosyl-l-methionine (AdoMet) to arginine residues. There are three types of PRMTs (I, II and III) that produce different methylation products, including asymmetric dimethylarginine (ADMA), symmetric dimethylarginine (SDMA) and monomethylarginine (MMA). Since these different methylations can lead to different biological consequences, understanding the origin of product specificity of PRMTs is of considerable interest. In this article, the quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed to study SDMA catalyzed by the Type II PRMT5 on the basis of experimental observation that the dimethylated product is generated through a distributive fashion. The simulations have identified some important interactions and proton transfers during the catalysis. Similar to the cases involving Type I PRMTs, a conserved Glu residue (Glu435) in PRMT5 is suggested to function as general base catalyst based on the result of the simulations. Moreover, our results show that PRMT5 has an energetic preference for the first methylation on Nη1 followed by the second methylation on a different ω-guanidino nitrogen of arginine (Nη2).The first and second methyl transfers are estimated to have free energy barriers of 19-20 and 18-19 kcal/mol respectively. The computer simulations suggest a distinctive catalytic mechanism of symmetric dimethylation that seems to be different from asymmetric dimethylation.

Understanding product specificity of protein lysine methyltransferases from QM/MM molecular dynamics and free energy simulations: the effects of mutation on SET7/9 beyond the Tyr/Phe switch.

The results of hybrid quantum mechanical/molecular mechanical (QM/MM) free energy (potential of mean force) simulations for methyl-transfer processes in SET7/9 and its Y245A mutant are compared to address the question concerning the change of the product specificity as well as catalytic efficiency due to the mutation. One of the key questions is whether or not the free energy profiles of methyl transfers may be used to predict the change of the product specificity as a result of the mutations for the residues that are not located at the Tyr/Phe switch position. The simulations show that while the wild-type SET7/9 is a monomethylase, the Y245→A mutation increases the ability of the enzyme to add more methyl groups on the target lysine (i.e., acting as a trimethylase). However, the first methyl-transfer process seems to become less efficient in the mutant compared to that in wild-type. All these results are consistent with experimental observations concerning the effects of the mutation on the product specificity and catalytic efficiency. Thus, the previous suggestion that the energetics of the methyl-transfer reactions may determine the product specificity, at least in some cases, is confirmed by the present work. Moreover, the dynamic information of the reactant complexes obtained from the QM/MM molecular dynamics simulations shows that the ability of the reactant complexes to form the reactive transition-state-like configurations may be used as an important indicator for the prediction of the product specificity of PKMTs, consistent with previous computational studies.

Dissecting molecular details and functional effects of the high-affinity copper binding site in plasminogen activator Inhibitor-1.

Plasminogen activator inhibitor-1 (PAI-1) is the primary inhibitor for plasminogen activators, tissue-type plasminogen activator (tPA) and urokinase-type plasminogen activator (uPA). As a unique member in the serine protease inhibitor (serpin) family, PAI-1 is metastable and converts to an inactive, latent structure with a half-life of 1-2 hr under physiological conditions. Unusual effects of metals on the rate of the latency conversion are incompletely understood. Previous work has identified two residues near the N-terminus, H2 and H3, which reside in a high-affinity copper-binding site in PAI-1 [Bucci JC, McClintock CS, Chu Y, Ware GL, McConnell KD, Emerson JP, Peterson CB (2017) J Biol Inorg Chem 22:1123-1,135]. In this study, neighboring residues, H10, E81, and H364, were tested as possible sites that participate in Cu(II) coordination at the high-affinity site. Kinetic methods, gel sensitivity assays, and isothermal titration calorimetry (ITC) revealed that E81 and H364 have different roles in coordinating metal and mediating the stability of PAI-1. H364 provides a third histidine in the metal-coordination sphere with H2 and H3. In contrast, E81 does not appear to be required for metal ligation along with histidines; contacts made by the side-chain carboxylate upon metal binding are perturbed and, in turn, influence dynamic fluctuations within the region encompassing helices D, E, and F and the W86 loop that are important in the pathway for the PAI-1 latency conversion. This investigation underscores a prominent role of protein dynamics, noncovalent bonding networks and ligand binding in controlling the stability of the active form of PAI-1.

Identification of a PAI-1-binding site within an intrinsically disordered region of vitronectin.

The serine protease inhibitor, plasminogen activator inhibitor Type-1 (PAI-1) is a metastable protein that undergoes an unusual transition to an inactive conformation with a short half-life of only 1-2 hr. Circulating PAI-1 is bound to a cofactor vitronectin, which stabilizes PAI-1 by slowing this latency conversion. A well-characterized PAI-1-binding site on vitronectin is located within the somatomedin B (SMB) domain, corresponding to the first 44 residues of the protein. Another PAI-1 recognition site has been identified with an engineered form of vitronectin lacking the SMB domain, yet retaining PAI-1 binding capacity (Schar, Blouse, Minor, Peterson. J Biol Chem. 2008;283:28487-28496). This additional binding site is hypothesized to lie within an intrinsically disordered domain (IDD) of vitronectin. To localize the putative binding site, we constructed a truncated form of vitronectin containing 71 amino acids from the N-terminus, including the SMB domain and an additional 24 amino acids from the IDD region. This portion of the IDD is rich in acidic amino acids, which are hypothesized to be complementary to several basic residues identified within an extensive vitronectin-binding site mapped on PAI-1 (Schar, Jensen, Christensen, Blouse, Andreasen, Peterson. J Biol Chem. 2008;283:10297-10309). Steady-state and stopped-flow fluorescence measurements demonstrate that the truncated form of vitronectin exhibits the same rapid biphasic association as full-length vitronectin and that the IDD hosts the elusive second PAI-1 binding site that lies external to the SMB domain of vitronectin.

QM/MM MD and Free Energy Simulation Study of Methyl Transfer Processes Catalyzed by PKMTs and PRMTs.

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here, we review the results of quantum mechanics/molecular mechanics molecular dynamics and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

QM/MM MD and free energy simulation study of methyl transfer processes catalyzed by PKMTs and PRMTs.

Methyl transfer processes catalyzed by protein lysine methyltransferases (PKMTs) and protein arginine methyltransferases (PRMTs) control important biological events including transcriptional regulation and cell signaling. One important property of these enzymes is that different PKMTs and PRMTs catalyze the formation of different methylated product (product specificity). These different methylation states lead to different biological outcomes. Here we review the results of quantum mechanics/molecular mechanics (QM/MM) molecular dynamics (MD) and free energy simulations that have been performed to study the reaction mechanism of PKMTs and PRMTs and the mechanism underlying the product specificity of the methyl transfer processes.

How Y357F, Y276F mutants affect the methylation activity of PRDM9: QM/MM MD and free energy simulations.

Histone methyltransferase PRDM9 catalyzes the methylation of H3K4me2 (histone 3 dimethylated lysine 4) to H3K4me3 (histone 3 trimethylated lysine 4) by transferring the methyl group from S-adenosyl methionine (AdoMet). PRDM9 is the major determinant of the meiotic recombination hotspot and the enrichment of H3K4me3 at the hotspot defines the initiation site of meiotic recombination. In PRDM9, two conserved tyrosine residues Tyr357 and Tyr276 surrounding the amino group of the substrate lysine may influence the methylation activity through hydrogen bond interactions with AdoMet or the substrate lysine. In this study, quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations were performed to reveal the methylation processes catalyzed by wild type PRDM9, its Y357F, and Y276F mutants, respectively. The different roles of Tyr357 and Tyr276 in the methylation activity of PRDM9 were also investigated and compared. The calculated free energy barriers of the methyl transfers suggest that the Y276F mutation decreases the catalytic activity of the methyl transfer, while the Y357F mutation does not change the catalytic activity of the methyl transfer. The reactant complex conformations generated in the QM/MM MD simulations show that the reactive configuration can be formed in the Y357F mutant but not in the Y276F mutant.

QM/MM MD and free energy simulations of G9a-like protein (GLP) and its mutants: understanding the factors that determine the product specificity.

Certain lysine residues on histone tails could be methylated by protein lysine methyltransferases (PKMTs) using S-adenosyl-L-methionine (AdoMet) as the methyl donor. Since the methylation states of the target lysines play a fundamental role in the regulation of chromatin structure and gene expression, it is important to study the property of PKMTs that allows a specific number of methyl groups (one, two or three) to be added (termed as product specificity). It has been shown that the product specificity of PKMTs may be controlled in part by the existence of specific residues at the active site. One of the best examples is a Phe/Tyr switch found in many PKMTs. Here quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations are performed on wild type G9a-like protein (GLP) and its F1209Y and Y1124F mutants for understanding the energetic origin of the product specificity and the reasons for the change of product specificity as a result of single-residue mutations at the Phe/Tyr switch as well as other positions. The free energy barriers of the methyl transfer processes calculated from our simulations are consistent with experimental data, supporting the suggestion that the relative free energy barriers may determine, at least in part, the product specificity of PKMTs. The changes of the free energy barriers as a result of the mutations are also discussed based on the structural information obtained from the simulations. The results suggest that the space and active-site interactions around the ε-amino group of the target lysine available for methyl addition appear to among the key structural factors in controlling the product specificity and activity of PKMTs.

The effects of nitrogen-heme-iron coordination on substrate affinities for cytochrome P450 2E1.

A descriptor based computational model was developed for cytochrome P450 2E1 (CYP2E1) based on inhibition constants determined for inhibition of chlorzoxazone, or 4-nitrophenol, metabolism. An empirical descriptor for type II binding was developed and tested for a series of CYP2E1 inhibitors. Inhibition constants where measured for 51 different compounds. A fast 2-dimensional predictive model was developed based on 40 compounds, and tested on 8 compounds of diverse structure. The trained model (n=40) had an r(2) value of 0.76 and an RMSE of 0.48. The correlation between the predicted and actual pK(i) values of the test set of compounds not included in the model gives an r(2) value of 0.78. The features that described binding include heme coordination (type II binding), molecular volume, octanol/water partition coefficient, solvent accessible surface area, and the sum of the atomic polarizabilities. The heme coordination parameter assigns an integer between 0 and 6 depending on structure, and is a new descriptor, based on simple quantum chemical calculations with correction for steric effects. The type II binding parameter was found to be important in obtaining a good correlation between predicted and experimental inhibition constants increasing the r(2) value from 0.38 to 0.77.