Understanding the Origins of Product Specificity of Protein Methyltransferases from QM/MM MD and Free Energy Simulations

Protein lysine methyltransferases (PKMTs) catalyze the methylation of certain lysine residues on histone tails using S-adenosyl-L-methionine (AdoMet) as the methyl donor. Regulation of chromatin structure and gene expression through histone lysine methylation depends on the degree of methylation. Therefore, it is of importance to understand the features of PKMTs that control how many methyl groups would be added to the target lysine (product specificity). In my dissertation, I have applied quantum mechanical/molecular mechanical (QM/MM) molecular dynamics (MD) and free energy simulations to understand the origin of the product specificities of SET8, GLP and their mutants. The relative free energy barriers of the methyl transfers calculated from our simulations are consistent with experimental data on product specificity, suggesting that they may govern the product specificity of PKMTs energetically. Our simulation results on the effects of the Tyr/Phe switch mutations in SET8 and GLP and the Y1124F mutation of GLP suggest that the space around the target lysine may be of importance. Protein arginine methyltransferases (PRMTs) catalyze the methylation of peptidyl arginines using AdoMet as the methyl donor. The product specificity of PRMTs often refers to their ability to form either asymmetric dimethylarginine (ADMA) (type I PRMTs) or symmetric dimethylarginine (SDMA) (type II PRMTs), and different product formation could lead to different biological consequences. The QM/MM MD and free energy simulations have been performed on one member of type I PRMTs, PRMT3. The simulation results suggest that the interactions of the active site residues make the same nitrogen atom of the guanidine to be the target of both 1st and 2nd methylations, leading to the formation of the asymmetric product. It is also proposed that Glu326 may function as a general base during the methyl transfer. The results of the MD simulations on the reactant complexes demonstrate that the ability of the reactant complexes to form the reactive configurations for the methyl transfer could be an important factor for the product specificity

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

Copyright © 2019 American Chemical Society. 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

© 2017, The Author(s). 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

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

Multiplexed Quantum Dot Labeling of Activated c-Met Signaling in Castration-Resistant Human Prostate Cancer

The potential application of multiplexed quantum dot labeling (MQDL) for cancer detection and prognosis and monitoring therapeutic responses has attracted the interests of bioengineers, pathologists and cancer biologists. Many published studies claim that MQDL is effective for cancer biomarker detection and useful in cancer diagnosis and prognosis, these studies have not been standardized against quantitative biochemical and molecular determinations. In the present study, we used a molecularly characterized human prostate cancer cell model exhibiting activated c-Met signaling with epithelial to mesenchymal transition (EMT) and lethal metastatic progression to bone and soft tissues as the gold standard, and compared the c-Met cell signaling network in this model, in clinical human prostate cancer tissue specimens and in a castration-resistant human prostate cancer xenograft model. We observed c-Met signaling network activation, manifested by increased phosphorylated c-Met in all three. The downstream survival signaling network was mediated by NF-κB and Mcl-1 and EMT was driven by receptor activator of NF-κB ligand (RANKL), at the single cell level in clinical prostate cancer specimens and the xenograft model. Results were confirmed by real-time RT-PCR and western blots in a human prostate cancer cell model. MQDL is a powerful tool for assessing biomarker expression and it offers molecular insights into cancer progression at both the cell and tissue level with high degree of sensitivity

The evolutionarily conserved long non‐coding RNA <i>LINC00261</i> drives neuroendocrine prostate cancer proliferation and metastasis <i>via</i> distinct nuclear and cytoplasmic mechanisms

Metastatic neuroendocrine prostate cancer (NEPC) is a highly aggressive disease, whose incidence is rising. Long noncoding RNAs (lncRNAs) represent a large family of disease- and tissue-specific transcripts, most of which are still functionally uncharacterized. Thus, we set out to identify the highly conserved lncRNAs that play a central role in NEPC pathogenesis. To this end, we performed transcriptomic analyses of donor-matched patient-derived xenograft models (PDXs) with immunohistologic features of prostate adenocarcinoma (AR+/PSA+) or NEPC (AR-/SYN+/CHGA+ ) and through differential expression analyses identified lncRNAs that were upregulated upon neuroendocrine transdifferentiation. These genes were prioritized for functional assessment based on the level of conservation in vertebrates. Here, LINC00261 emerged as the top gene with over 3229-fold upregulation in NEPC. Consistently, LINC00261 expression was significantly upregulated in NEPC specimens in multiple patient cohorts. Knockdown of LINC00261 in PC-3 cells dramatically attenuated its proliferative and metastatic abilities, which are explained by parallel downregulation of CBX2 and FOXA2 through distinct molecular mechanisms. In the cell cytoplasm, LINC00261 binds to and sequesters miR-8485 from targeting the CBX2 mRNA, while inside the nucleus, LINC00261 functions as a transcriptional scaffold to induce SMAD-driven expression of the FOXA2 gene. For the first time, these results demonstrate hyperactivation of the LINC00261-CBX2-FOXA2 axes in NEPC to drive proliferation and metastasis, and that LINC00261 may be utilized as a therapeutic target and a biomarker for this incurable disease

The non-coding transcriptome as a dynamic regulator of cancer metastasis.

Since the discovery of microRNAs, non-coding RNAs (NC-RNAs) have increasingly attracted the attention of cancer investigators. Two classes of NC-RNAs are emerging as putative metastasis-related genes: long non-coding RNAs (lncRNAs) and small nucleolar RNAs (snoRNAs). LncRNAs orchestrate metastatic progression through several mechanisms, including the interaction with epigenetic effectors, splicing control and generation of microRNA-like molecules. In contrast, snoRNAs have been long considered "housekeeping" genes with no relevant function in cancer. However, recent evidence challenges this assumption, indicating that some snoRNAs are deregulated in cancer cells and may play a specific role in metastasis. Interestingly, snoRNAs and lncRNAs share several mechanisms of action, and might synergize with protein-coding genes to generate a specific cellular phenotype. This evidence suggests that the current paradigm of metastatic progression is incomplete. We propose that NC-RNAs are organized in complex interactive networks which orchestrate cellular phenotypic plasticity. Since plasticity is critical for cancer cell metastasis, we suggest that a molecular interactome composed by both NC-RNAs and proteins orchestrates cancer metastasis. Interestingly, expression of lncRNAs and snoRNAs can be detected in biological fluids, making them potentially useful biomarkers. NC-RNA expression profiles in human neoplasms have been associated with patients' prognosis. SnoRNA and lncRNA silencing in pre-clinical models leads to cancer cell death and/or metastasis prevention, suggesting they can be investigated as novel therapeutic targets. Based on the literature to date, we critically discuss how the NC-RNA interactome can be explored and manipulated to generate more effective diagnostic, prognostic, and therapeutic strategies for metastatic neoplasms

Association Analysis of Driver Gene–Related Genetic Variants Identified Novel Lung Cancer Susceptibility Loci with 20,871 Lung Cancer Cases and 15,971 Controls

BACKGROUND:A substantial proportion of cancer driver genes (CDG) are also cancer predisposition genes. However, the associations between genetic variants in lung CDGs and the susceptibility to lung cancer have rarely been investigated. METHODS:We selected expression-related single-nucleotide polymorphisms (eSNP) and nonsynonymous variants of lung CDGs, and tested their associations with lung cancer risk in two large-scale genome-wide association studies (20,871 cases and 15,971 controls of European descent). Conditional and joint association analysis was performed to identify independent risk variants. The associations of independent risk variants with somatic alterations in lung CDGs or recurrently altered pathways were investigated using data from The Cancer Genome Atlas (TCGA) project. RESULTS:We identified seven independent SNPs in five lung CDGs that were consistently associated with lung cancer risk in discovery (P < 0.001) and validation (P < 0.05) stages. Among these loci, rs78062588 in TPM3 (1q21.3) was a new lung cancer susceptibility locus (OR = 0.86, P = 1.65 × 10-6). Subgroup analysis by histologic types further identified nine lung CDGs. Analysis of somatic alterations found that in lung adenocarcinomas, rs78062588[C] allele (TPM3 in 1q21.3) was associated with elevated somatic copy number of TPM3 (OR = 1.16, P = 0.02). In lung adenocarcinomas, rs1611182 (HLA-A in 6p22.1) was associated with truncation mutations of the transcriptional misregulation in cancer pathway (OR = 0.66, P = 1.76 × 10-3). CONCLUSIONS:Genetic variants can regulate functions of lung CDGs and influence lung cancer susceptibility. IMPACT:Our findings might help unravel biological mechanisms underlying lung cancer susceptibility

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