Next Article in Journal
A 1H-NMR-Based Metabonomic Study on the Anti-Depressive Effect of the Total Alkaloid of Corydalis Rhizoma
Previous Article in Journal
Synthesis and Physicochemical Characterization of the Impurities of Pemetrexed Disodium, an Anticancer Drug
Previous Article in Special Issue
QM/MM Calculations with deMon2k
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 20
Issue 6
10.3390/molecules200610032
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Yufei Yue
1,
Yuzhuo Chu
2 and
Hong Guo
1,3,*
1
Department of Biochemistry and Cellular and Molecular Biology, University of Tennessee, Knoxville, TN 37996, USA
2
School of Life Science and Biotechnology, Dalian University of Technology, Dalian 116024, China
3
UT/ORNL Center for Molecular Biophysics, Oak Ridge National Laboratory, Oak Ridge, TN 37830, USA
*
Author to whom correspondence should be addressed.
Molecules 2015, 20(6), 10032-10046; https://doi.org/10.3390/molecules200610032
Submission received: 3 March 2015 / Revised: 18 May 2015 / Accepted: 25 May 2015 / Published: 29 May 2015
(This article belongs to the Special Issue Multiscale Chemical Modeling Using Quantum Mechanics/Molecular Mechanics (QM/MM) Methods: Advances and Applications)
Browse Figures
Versions Notes

Abstract

:
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.

1. Introduction

Posttranslational methylation of histone proteins on their arginine residues is an epigenetic mark that plays a vital role in cell function and is related with cell disorders and diseases [1]. Protein arginine methyltransferases (PRMTs) catalyze the transfer of methyl group(s) from S-adenosyl-l-methionine (AdoMet) to the guanidine nitrogen of arginine residue, resulting in the reaction products containing methylarginine along with S-adenosylhomocysteine (SAH) [2]. In addition to the methylation of histone proteins, PRMTs could also modify a variety of other proteins [3]. Depending on types of PRMTs (I, II or III), the methylation products may contain asymmetric dimethylarginine (ADMA), symmetric dimethylarginine (SDMA) or monomethylarginine (MMA) (shown in Scheme 1) [1].
Molecules 20 10032 g009 550
Scheme 1. Methylation of arginine by different types of PRMTs. Type I PRMT can produce both monomethylarginine (MMA) and asymmetric dimethylarginine (ADMA). Type II PRMT can produce both MMA and symmetric dimethylarginine (SDMA). Type III PRMT can only produce MMA.
Scheme 1. Methylation of arginine by different types of PRMTs. Type I PRMT can produce both monomethylarginine (MMA) and asymmetric dimethylarginine (ADMA). Type II PRMT can produce both MMA and symmetric dimethylarginine (SDMA). Type III PRMT can only produce MMA.
Molecules 20 10032 g009
The Type II PRMT5 to be investigated in this work catalyzes the transfer of methyl groups to the two different ω-guanidino nitrogen atoms on the arginine residue of the target protein, producing the ω-NG, N′G symmetrically dimethylated arginine (SDMA) [4]. PRMT5 has a variety of substrates that include histones, transcription factors, splicesomal proteins and piRNA biogenesis related proteins, and this enzyme functions in both nucleus and cytoplasm [3]. SDMA may profoundly impact many biological processes including epigenetic control of gene expression [5], circadian rhythms [6,7], splicing regulation [8,9], germ cell development and pluripotency [10,11], and DNA damage response [12,13]. The Type II PRMT5, however, often share the common recognition sequence with the Type I PRMTs which add two methyl groups to the same ω-guanidino nitrogen atom (ADMA) [14]. Thus, the same target arginine may be either symmetrically or asymmetrically dimethylated. As such, SDMA and ADMA are isomeric protein posttranslational modifications with distinct and sometimes reversible biological effects [14]. One example is the methylation of arginine 3 on histone H4 (H4R3): symmetric dimethylation of H4R3 could repress gene expression [15,16], while asymmetric dimethylation of H4R3 is correlated with gene activation [17,18]. Recent studies have focused on understanding the enzymatic mechanisms that differentiate the two chemically isomeric but functionally antagonistic posttranslational modifications [14].
The crystal structures have been determined for several PRMTs, including PRMT1 [19], PRMT3 [20], PRMT6 [21], PRMT10 [22] of Type I, PRMT5 [23] of Type II and PRMT7 [24] of Type III. Computer simulations have been applied to investigate the catalytic mechanism of PRMT1 [25] and PRMT3 [26], and some important questions concerning the product specificity of Type I ADMA have been addressed. The computational approaches used in these earlier studies include molecular dynamics and free energy simulations (potential of mean force) with the hybrid quantum mechanical and molecular mechanical (QM/MM) potential that seem to be suitable to investigate the enzyme-catalyzed methyl transfer process and have also been widely used for some other methyltransferases [26,27,28,29,30,31,32].
In this study, the methylation reaction catalyzed by Type II PRMT5 is investigated by use of QM/MM MD and free energy simulations 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 involving the active site residues. Similar to the cases of Type I PRMT1 [25] and PRMT3 [26], a homologous Glu residue (Glu435) in PRMT5 seems to function as a general base catalyst during the catalysis, and the corresponding proton transfer is found to be somehow concerted with the methyl transfer process. However, unlike Type I PRMTs which energetically favors a single ω-guanidino nitrogen (Nη2) of arginine as the target for the both 1st and 2nd methylations [25,26], PRMT5 is found to have an energetic preference of targeting Nη1 for the first methylation and then targeting a different ω-guanidino nitrogen (Nη2) for the second methylation. The first and second methyl transfers are estimated to have free energy barriers of 19–20 and 18–19 kcal/mol, respectively, from the simulations. These results are consistent with the existing experimental data [23]. Our computational study provides a better understanding of the symmetric di-methylation mechanism that is different from the asymmetric di-methylation.

2. Results and Discussion

2.1. Comparison of the Active Site Structures of Type II PRMT5 and Type I PRMT3

The two invariant glutamate residues from the “double-E” loop found in all PRMTs [33] are structurally conserved in both Type I PRMT3 and Type II PRMT5 active sites (E435 and E444 in Figure 1A; E144 and E153 in Figure 1B). As is shown in Figure 1, the carboxylate side chain of E444 in PRMT5 (E153 in PRMT3) forms a stable salt bridge with the Nε and Nη1 atoms of the substrate arginine. The side chain of E435 in PRMT5 (E144 PRMT3) interacts with Nη1 and Nη2 of the substrate. These two Glu residues are required for the enzymatic activities; the mutation of either of them could greatly decrease the enzymatic activity [14].
Molecules 20 10032 g001 550
Figure 1. Comparison of the active sites of the crystal structures for Type-II PRMT5 (PDB ID: 4QGB) and Type-I PRMT3 (PDB ID:1F3L). (A) The active site of PRMT5; (B).The active site of PRMT3. The two glutamate residues, E435 (E144) and E444 (E153) in PRMT5 (PRMT3), are conserved in the both types of PRMTs. F327, K333, S578 and S439 in PRMT5 are the conserved residues among Type-II PRMTs. The corresponding residues in Type-I PRMT3 are M48, R54, H293 and Y148, respectively. The substrate arginine is labeled as Sub. Some distances are shown in blue with the unit of angstrom.
Figure 1. Comparison of the active sites of the crystal structures for Type-II PRMT5 (PDB ID: 4QGB) and Type-I PRMT3 (PDB ID:1F3L). (A) The active site of PRMT5; (B).The active site of PRMT3. The two glutamate residues, E435 (E144) and E444 (E153) in PRMT5 (PRMT3), are conserved in the both types of PRMTs. F327, K333, S578 and S439 in PRMT5 are the conserved residues among Type-II PRMTs. The corresponding residues in Type-I PRMT3 are M48, R54, H293 and Y148, respectively. The substrate arginine is labeled as Sub. Some distances are shown in blue with the unit of angstrom.
Molecules 20 10032 g001
Previous computational studies have identified E144 as the general base to accept proton from the arginine during the methyl transfer catalyzed by Type I PRMTs [25,26]. One interesting questions is whether the corresponding E435 from PRMT5 would play a similar role during the catalysis. Four other residues that are conserved in the active site of PRMT5 are F327, K333, S578 and S439; the corresponding residues in PRMT3 are M48, R54, Y148 and H293, respectively. Although the mutation of S439 and S578 of PRMT5 diminishes the enzymatic activity significantly [14], the exact role of these residues is not clear; both S439 and S578 seem to be far away from the methyl donating AdoMet (Figure 1A). Interestingly, the F327M mutant of PRMT5 could produce both ADMA and SDMA [14], indicating that F327 may occupy a key position in PRMT5 and its properties may be important in determining the product specificity. Moreover, K333 is in the vicinity of F327 and could form hydrogen bonds with both the carboxylate group of AdoMet and E435.
The previous computational studies have shown that for PRMT1 and PRMT3 the 1st and 2nd methyl transfers would be energetically more favorable with relatively lower barriers if Nη2 (see Figure 1B) is the methyl acceptor in each of the cases [25,26]. This is consistent with the fact that they are both Type I PRMTs and the methylation products contain asymmetric dimethylarginine (ADMA). In the active site of the PRMT5 structure (Figure 1A), Nη1 appears to be in a much better position for acting as the methyl acceptor compared to Nη2; this is in contrast with the case of PRMT3 where Nη2 seems to be in a better position for accepting the methyl group [26]. Indeed, the distance between Nη1 and CM is 3.2 Å in the active site of PRMT5 as compared to 4.9 Å between Nη2 and CM.

2.2. The First Methylation Catalyzed by PRMT5

The acceptor site for the first methylation process catalyzed by PRMT5 is studied by comparing the energetics of the methyl addition to Nη1 and Nη2, respectively, on the arginine via QM/MM MD and free energy simulations. As shown in Figure 2B, the free energy barrier for the methyl transfer to Nη1 is 20.4 kcal/mol, which is 9 kcal/mol lower than that for the methyl transfer to Nη2. This suggests that the 1st methyl group is likely to be transferred to the Nη1 atom, as the corresponding process is energetically more favorable. This is in contrast with the 1st methylation catalyzed by PRMT3 which favors Nη2 as the methyl acceptor [26].
Molecules 20 10032 g002 550
Figure 2. (A) Two possible sites for the 1st methylation. The substrate arginine is labeled as Sub; (B) Free energy profiles of the 1st methylation on Nη1 and Nη2. To Nη1: blue line with a free energy barrier of 20.4 kcal/mol; to Nη2: red dotted line with a free energy barrier of 29.4 kcal/mol. Error bars are shown on the free energy profile (within ± 0.12 kcal/mol).
Figure 2. (A) Two possible sites for the 1st methylation. The substrate arginine is labeled as Sub; (B) Free energy profiles of the 1st methylation on Nη1 and Nη2. To Nη1: blue line with a free energy barrier of 20.4 kcal/mol; to Nη2: red dotted line with a free energy barrier of 29.4 kcal/mol. Error bars are shown on the free energy profile (within ± 0.12 kcal/mol).
Molecules 20 10032 g002
The average active-site structures of the reactant complex and near transition state for the 1st methyl transfer are shown in Figure 3 and Figure 4 for the methyl transfer to Nη1 and Nη2, respectively. The interactions obtained from the simulations for the reactant complex generally resemble the interactions observed in the crystal structure. This seems to indicate that the conformation of the active site in the crystal structure is well maintained after initial 1.5 ns MD simulations. In the reactant complex for the methyl transfer to Nη1 (Figure 3A), Nη1 is well aligned with the methyl group of AdoMet, and this good alignment presumably leads to a relatively low barrier for the methyl transfer. The arginine on the substrate appears to be well stabilized through the salt bridge with E444 as well as hydrogen bonding interaction with E435. K333 may also help to adjust the orientation of E435 through a hydrogen bonding interaction. F327 seems to be involved in the π-cation interaction with the arginine in the reactant state. Moreover, as mentioned earlier the bulk size of this residue may prevent the formation of a good alignment between the transferable methyl group and Nη2 and interfere with the methyl transfer to Nη2. Indeed, the structures in Figure 1A shows that the distance between the transferable methyl group (CM) and Nη2 can be as much as 4.9 Å in the reactant complex. By contrast, the corresponding distance in PRMT3, which has M48 at the location rather than a Phe residue, is much smaller (e.g., 3.9 Å in Figure 1B). It is of interest to note from Figure 4B that the distance between F327 and the arginine increases significantly near the transition state during the methyl transfer to Nη2.
Molecules 20 10032 g003 550
Figure 3. The average active-site structures obtained from the MD and free energy simulations for the 1st methylation on Nη1. (A) Reactant complex; (B) Near the transition state. Some average distances in active sites are shown in angstrom. Note that E435 works as the proton acceptor based on the simulations.
Figure 3. The average active-site structures obtained from the MD and free energy simulations for the 1st methylation on Nη1. (A) Reactant complex; (B) Near the transition state. Some average distances in active sites are shown in angstrom. Note that E435 works as the proton acceptor based on the simulations.
Molecules 20 10032 g003
Molecules 20 10032 g004 550
Figure 4. The average active-site structure obtained from the MD and free energy simulations for the 1st methylation on Nη2. (A) Reactant complex; (B) Near the transition state.
Figure 4. The average active-site structure obtained from the MD and free energy simulations for the 1st methylation on Nη2. (A) Reactant complex; (B) Near the transition state.
Molecules 20 10032 g004
The 1D free energy simulations only used the reaction coordinate for the methyl transfer. Nevertheless, the proton transfer occurred near the transition state of methylation (Figure 3B and Figure 4B). To better understand the relationship between the methyl transfer and the proton transfer process, the 2D free energy simulations were performed (Figure 5). As is shown in Figure 5, the reaction path goes from point A (designated as the reactant complex) to point B (the transition state), and finally reaches point D as the product complex. Before reaching the transition state at point B, the proton has been basically transferred from Nη1 of the substrate arginine to the carboxyl oxygen of E435. And the free energy barrier is estimated to be 19–20 kcal/mol from the 2D free energy simulations, which is almost the same as that obtained from the 1D free energy simulations (20.4 kcal/mol). Moreover, point D in Figure 5 is about 6 kcal/mol lower than point C, suggesting that the deprotonated product is more stable in the active site than the protonated product and that the proton transfer is coupled with the methyl transfer during methylation.
Molecules 20 10032 g005 550
Figure 5. 2D free-energy contour map for the 1st methylation on Nη1. Rx: the reaction coordinate of the methyl transfer; Ry: the reaction coordinate for the proton transfer. Points A, B, and D designate the reactant, near transition state, and product complexes, respectively. Some contour lines are shown with respective energy values. The energy barrier is estimated to be around 19–20 kcal/mol.
Figure 5. 2D free-energy contour map for the 1st methylation on Nη1. Rx: the reaction coordinate of the methyl transfer; Ry: the reaction coordinate for the proton transfer. Points A, B, and D designate the reactant, near transition state, and product complexes, respectively. Some contour lines are shown with respective energy values. The energy barrier is estimated to be around 19–20 kcal/mol.
Molecules 20 10032 g005

2.3. The Second Methylation Catalyzed by PRMT5

With the first methyl group being transferred already, the 2nd methylation reaction can now be examined to determine the mechanism for the formation of the symmetrically dimethylated product. Previous experimental study on Caenorhabditis elegans PRMT5 (cPRMT5) has shown that the dimethylated product is generated through a distributive fashion [4] in which the peptide is released prior to rebinding to facilitate a second round of methylation. The distributive mechanism was further confirmed from the study of the human hPRMT5•MEP50 complex [34]. This is in contrast with the cases of some protein lysine methyltransferases for which the multiple rounds of methylation are believed to proceed processively without the release of the intermediates from the active sites. The kinetic data (Km and kcat) for cPRMT5, hPRMT5 and hPRMT5•MEP50 complex have also been obtained for a variety of un-methylated and mono-methylated substrates. For the hPRMT5•MEP50 complex for which the current investigation is based on, the differences in Km and kcat for the first and second methyl transfers are rather small and are beyond the accuracy of most of quantum mechanical approaches used in QM/MM studies. Although several experimental structures for PRMTs are available, the structure for the reactant complex of the second methyl transfer has not been determined. In our earlier investigations of protein lysine methyltransferases [27,28,29,30], we developed an approach based on the free energy simulations of the methyl transfer processes that allows us to determine whether or not the enzymes would be able to add methyl groups. This approach has been successfully applied to study the product specificity for a number of protein lysine methyltransferases, and its usefulness has also been confirmed on some other methyltransferases [26,31,32], including PRMT1 [26]. This approach is especially suitable for the current investigation because the Km values for the first and second methyl transfers are rather similar (see above), and the hypothetical non-existent methyl transfers can be presumably eliminated using this approach in the similar manner as the product specificity is determined [26,27,28,29,30]. In Table 1, the three different configurations for the mono-methyl arginine without destroying the salt bridge with E444 are given. The nitrogen atom already connecting to the methyl group is designated as Nη1 and that without the methyl group (i.e., NH2) designated as Nη2. As is demonstrated in Table 1, for each configuration there are two possible di-methylation products (SDMA and ADMA). The QM/MM free energy simulations have been performed for each of the six methylation processes. Table 1(I) shows that the free energy barriers were calculated to be 32.3 and 36 kcal/mol for the 2nd methyl transfer to the Nη1 and Nη2 atoms, respectively. Since either of these energy barriers seems to be too high, the 2nd methylation may not start from the configuration in Table 1(I). The similar argument can be made for the configuration in Table 1(III).
Table
Table 1. Free energy barriers for the 2nd methylation on either the mono-methylated nitrogen atom (designated as Nη1) or un-methylated nitrogen atom (designated as Nη2) of the mono-methylated arginine (MMA) with three possible configurations for the reactant complex without destroying the salt bridge involving E444. The only difference between Configurations II and III is whether the methyl group is pointing to the AdoMet side (III) or away from the AdoMet side (II). The conformation in (II) has the lowest energy barrier of 20.1 kcal/mol when the 2nd methylation occurs to Nη2 which will eventually produce the symmetrically dimethylated arginine (SDMA). Error bars are also determined, similar to those given in Figure 2 (within ± 0.2 kcal/mol).
Table 1. Free energy barriers for the 2nd methylation on either the mono-methylated nitrogen atom (designated as Nη1) or un-methylated nitrogen atom (designated as Nη2) of the mono-methylated arginine (MMA) with three possible configurations for the reactant complex without destroying the salt bridge involving E444. The only difference between Configurations II and III is whether the methyl group is pointing to the AdoMet side (III) or away from the AdoMet side (II). The conformation in (II) has the lowest energy barrier of 20.1 kcal/mol when the 2nd methylation occurs to Nη2 which will eventually produce the symmetrically dimethylated arginine (SDMA). Error bars are also determined, similar to those given in Figure 2 (within ± 0.2 kcal/mol).
Molecules 20 10032 i001
To Nη1To Nη2To Nη1To Nη2To Nη1To Nη2
32.3 kcal/mol36 kcal/mol31.3 kcal/mol20.1 kcal/mol43.9 kcal/mol32.5 kcal/mol
Table 1(II) shows that the second methyl transfer has the lowest free energy barrier (20.1 kcal/mol) if it is transferred to Nη2 (i.e., the SDMA product), while the free energy barrier for the second methyl transfer to Nη1 (i.e., the ADMA product) based on the configuration in Table 1(II) is 11.2 kcal/mol higher.
The average structures of the reactant complex and near transition state for the 2nd methyl transfer to Nη2 and Nη1 are exhibited in Figure 6 and Figure 7, respectively. The active-site structures of the reactant complex (Figure 6A) show that the Nη2 atom is well aligned with the methyl group of AdoMet. The relatively high barrier for the 2nd methyl transfer to Nη1 appears due in part to the steric hindrance of F327 (as in the case of the 1st methyl transfer), although other factors may be involved as well. One interesting observation for the second methyl transfer is that the proton has not been transferred to the general base E435 at the transition state (see Figure 6B and Figure 7B). This is in contrast with the case of the first methyl transfer (shown in Figure 3B and Figure 4B) where the proton transfer occurs before the methylation reaches the transition state. This observation is also reflected in the 2D free energy map involving both the 2nd methyl transfer and the proton transfer (Figure 8). As is shown in the reaction path illustrated in Figure 8, the reaction reaches the transition state at point B without the deprotonation of the monomethylarginine.
Molecules 20 10032 g006 550
Figure 6. Results from the simulations for the 2nd methylation on Nη2. (A) Active-site structure of the reactant complex; (B) Active-site structure near the transition state.
Figure 6. Results from the simulations for the 2nd methylation on Nη2. (A) Active-site structure of the reactant complex; (B) Active-site structure near the transition state.
Molecules 20 10032 g006
Molecules 20 10032 g007 550
Figure 7. Results from the MD simulations for the 2nd methylation on Nη1. (A) Active-site structure of the reactant complex; (B) Active-site structure near the transition state.
Figure 7. Results from the MD simulations for the 2nd methylation on Nη1. (A) Active-site structure of the reactant complex; (B) Active-site structure near the transition state.
Molecules 20 10032 g007
Molecules 20 10032 g008 550
Figure 8. 2D free-energy contour map for the 2nd methylation on Nη2. Points A, B, and D designate the reactant, near transition state, and product complex, respectively. Point C corresponds to the hypothetical product complex with proton not being transferred, which is only 2 kcal/mol higher than Point D in the free-energy contour map. Some contour lines are shown with respective energy value in kcal/mol. The energy bar is shown on the right. The energy barrier is estimated to be around 18–19 kcal/mol.
Figure 8. 2D free-energy contour map for the 2nd methylation on Nη2. Points A, B, and D designate the reactant, near transition state, and product complex, respectively. Point C corresponds to the hypothetical product complex with proton not being transferred, which is only 2 kcal/mol higher than Point D in the free-energy contour map. Some contour lines are shown with respective energy value in kcal/mol. The energy bar is shown on the right. The energy barrier is estimated to be around 18–19 kcal/mol.
Molecules 20 10032 g008
Another interesting observation from Figure 6B and Figure 7B is that E435 forms a salt bridge to both Nη1 and Nη2 near the transition state. This salt bridge may be relatively strong compared to that in the reactant complex and help to stabilize the transition state; it may also play a role for the delayed proton transfer mentioned earlier. The free energy barrier is estimated to be 18–19 kcal/mol by the 2D free energy simulation (Figure 8), consistent with 20.1 kcal/mol from the 1D free energy simulations (Table 1).

3. Experimental Section

The simulation coordinates for the reactant complexes of methyl transfers were based on the crystallographic complex (PDB ID: 4GQB, 2.06 Å) of PRMT5 [23] that contains AdoMet analog and the H4 peptide with the methylation targeting substrate, Arg3. The coordinates of PRMT3 for comparison is from the crystal structure (PDB ID:1F3L, 2.06 Å) and modified based on the previous work [26]. The QM/MM MD and free energy (potential of mean force, PMF)) simulations were applied for monitoring the methylation processes and determining the free energy profiles with the CHARMM program [35]. A water sphere based on a modified TIP3P water model [36] with radius(r) of 30 Å, centered at CZ of Arg3, was pre-equilibrated to the system. A stochastic boundary with a Poisson-Boltzmann charge-scaling scheme [37] was applied for the model. The reservoir region had r > 22 Å, and the buffer region had r equal to 20 Å ≤ r ≤ 22 Å. The reaction region had r ≤ 20 Å. The -CH2-CH2-S+(Me)-CH2- part of AdoMet, the side chain of substrate Arg3/monomethylated Arg3, the side chains of E435 and E444 were treated by QM and the rest of the system by MM. The resulting systems contained around 5500 atoms with about 800 water molecules. The all-hydrogen potential function (PARAM27) [38] was used for the MM region, the self-consistence charge density functional tight binding (SCC-DFTB) [39,40] method was used for the QM region. The link-atom approach [41] was applied to separate the QM and MM regions.
The initial structures for the entire resulting system were optimized by use of the steepest descent (SD) and adopted-basis Newton-Raphson (ABNR) methods. The systems were gradually heated from 50 to 310.15 K in 50 ps. A 1 fs time step was applied to integrate equation of motion. 1.5 ns QM/MM MD simulations were initially executed for each of the reactant complex. The reaction coordinate was defined as a linear combination of r(CM ... Sδ) and r(CM ··· Nη1/2), which is R = r(CM ··· Sδ) − r(CM ··· Nη1/2). The umbrella sampling method [42] in the CHARMM program along with the weighted histogram analysis method (WHAM) [43] was used to determine the free energy(PMF) change as a function of the reaction coordinate(s). The Monte Carlo Simulation bootstrapping integrated with WHAM [43] was applied 100 times to estimate the errors (shown in Figure 2 and indicated in the legend of Table 1) of the PMF profile. 20–22 simulation windows were saved for each methyl transfer process. And for each window 50 ps production runs were performed after 50 ps equilibration. The force constants of the harmonic biasing potentials used in the PMF simulations were 50–500 kcal·mol−1·Å−2. During the 50-ps production run of each window, structural coordinates were saved every 50 steps (1 fs/step), and consequently there are 1000 frames saved for each window. The distances shown in Figure 3, Figure 4, Figure 6 and Figure 7 are the means of corresponding distances with stand deviation (S.D.) collected from those 1000 frames. The 2D free energy (PMF) map were also determined with umbrella sampling method and two-dimensional-WHAM. The time step for simulations is 1 fs. The horizontal reaction coordinate for the 1st methylation on Nη1, Rx = r(CM ... Sδ) − r(CM … Nη1), was used to describe the methyl transfer process. And the vertical reaction coordinate, Ry = r(HH12 ... Nη1) − r(HH12 … OE1), was used to describe the proton transfer process. Similarly, for the 2nd methylations, the methyl transfer was explained by Rx = r(CM ... Sδ) − r(CM … Nη2) and the proton transfer by Ry = r(HH22 ... Nη2) − r(HH22 … OE1). Approximated 300–400 windows were used in the construction of 2D free energy map, with 50 ps production run following 50 ps equilibration for each window. The force constants for each window were in the range of 100–800 kcal·mol−1·Å−2.

4. Conclusions

In this study, the Type II PRMT5 was investigated by the QM/MM MD and free energy simulations. E435 was found to function as a general base catalyst for accepting a proton from the substrate arginine during the 1st and the 2nd methylation. This residue seems to play a similar role as E144 in Type I PRMT3. The simulations provide the detailed mechanism for the symmetric di-methylation by the enzyme and suggest the possible role for some other key residues as well during the catalysis. The kcat for the arginine methylation catalyzed by PRMT5 is measured to be between 10 to 50 h−1, indicating that the activation energy barrier is about 18 kcal/mol based on transition state theory [23]. Thus, our estimates of the free energy barriers for the methyl transfers (19–20 and 18–19 kcal/mol for the first and the second methylation, respectively) are reasonable compared with the experimental observations. Furthermore, the simulations suggest that the symmetric di-methylation by PRMT5 seems to be energetically favorable, in agreement with the fact that PRMT5 is Type II PRMT. The proposed mechanism here based on the simulation results is different from the asymmetric dimethylation [25,26].

Acknowledgments

We thank Martin Karplus for a gift of the CHARMM program. This work was supported by the National Science Foundation (Grant 0817940 to H.G.). We are also grateful for the computer resources (Newton) from University of Tennessee, Knoxville and from the Extreme Science and Engineering Discovery Environment (XSEDE), which is supported by National Science Foundation grant number ACI-1053575.

Author Contributions

YY and HG designed, performed research, analyzed the data and wrote the paper. CY contributed with the construction of initial PRMT3 structure. All authors read and approved the final manuscript.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bedford, M.T.; Clarke, S.G. Protein arginine methylation in mammals: Who, what, and why. Mol. Cell 2009, 33, 1–13. [Google Scholar] [CrossRef] [PubMed]
  2. Di Lorenzo, A.; Bedford, M.T. Histone arginine methylation. FEBS Lett. 2011, 585, 2024–2031. [Google Scholar] [CrossRef] [PubMed]
  3. Biggar, K.K.; Li, S.S. Non-histone protein methylation as a regulator of cellular signalling and function. Nat. Rev. Mol. Cell Biol. 2015, 16, 5–17. [Google Scholar] [CrossRef] [PubMed]
  4. Wang, M.; Xu, R.M.; Thompson, P.R. Substrate specificity, processivity, and kinetic mechanism of protein arginine methyltransferase 5. Biochemistry 2013, 52, 5430–5440. [Google Scholar] [CrossRef] [PubMed]
  5. Wysocka, J.; Allis, C.D.; Coonrod, S. Histone arginine methylation and its dynamic regulation. Front. Biosci. 2006, 11, 344–355. [Google Scholar] [CrossRef] [PubMed]
  6. Sanchez, S.E.; Petrillo, E.; Beckwith, E.J.; Zhang, X.; Rugnone, M.L.; Hernando, C.E.; Cuevas, J.C.; Godoy Herz, M.A.; Depetris-Chauvin, A.; Simpson, C.G.; et al. A methyl transferase links the circadian clock to the regulation of alternative splicing. Nature 2010, 468, 112–116. [Google Scholar] [CrossRef] [PubMed]
  7. Hong, S.; Song, H.R.; Lutz, K.; Kerstetter, R.A.; Michael, T.P.; McClung, C.R. Type II protein arginine methyltransferase 5 (prmt5) is required for circadian period determination in arabidopsis thaliana. Proc. Natl. Acad. Sci. USA 2010, 107, 21211–21216. [Google Scholar] [CrossRef] [PubMed]
  8. Meister, G.; Eggert, C.; Buhler, D.; Brahms, H.; Kambach, C.; Fischer, U. Methylation of sm proteins by a complex containing prmt5 and the putative u snrnp assembly factor picln. Curr. Biol. 2001, 11, 1990–1994. [Google Scholar] [CrossRef]
  9. Deng, X.; Gu, L.; Liu, C.; Lu, T.; Lu, F.; Lu, Z.; Cui, P.; Pei, Y.; Wang, B.; Hu, S.; et al. Arginine methylation mediated by the arabidopsis homolog of prmt5 is essential for proper pre-mrna splicing. Proc. Natl. Acad. Sci. USA 2010, 107, 19114–19119. [Google Scholar] [CrossRef] [PubMed]
  10. Nishida, K.M.; Okada, T.N.; Kawamura, T.; Mituyama, T.; Kawamura, Y.; Inagaki, S.; Huang, H.; Chen, D.; Kodama, T.; Siomi, H.; et al. Functional involvement of tudor and dprmt5 in the pirna processing pathway in drosophila germlines. EMBO J. 2009, 28, 3820–3831. [Google Scholar] [CrossRef] [PubMed]
  11. Gonsalvez, G.B.; Rajendra, T.K.; Tian, L.; Matera, A.G. The sm-protein methyltransferase, dart5, is essential for germ-cell specification and maintenance. Curr. Biol. 2006, 16, 1077–1089. [Google Scholar] [CrossRef] [PubMed]
  12. Jansson, M.; Durant, S.T.; Cho, E.C.; Sheahan, S.; Edelmann, M.; Kessler, B.; la Thangue, N.B. Arginine methylation regulates the p53 response. Nat. Cell Biol. 2008, 10, 1431–1439. [Google Scholar] [CrossRef] [PubMed]
  13. Yang, M.; Sun, J.; Sun, X.; Shen, Q.; Gao, Z.; Yang, C. Caenorhabditis elegans protein arginine methyltransferase prmt-5 negatively regulates DNA damage-induced apoptosis. PLoS Genet. 2009, 5, e1000514. [Google Scholar] [CrossRef] [PubMed]
  14. Sun, L.; Wang, M.; Lv, Z.; Yang, N.; Liu, Y.; Bao, S.; Gong, W.; Xu, R.M. Structural insights into protein arginine symmetric dimethylation by prmt5. Proc. Natl. Acad. Sci. USA 2011, 108, 20538–20543. [Google Scholar] [CrossRef] [PubMed]
  15. Wang, X.; Zhang, Y.; Ma, Q.; Zhang, Z.; Xue, Y.; Bao, S.; Chong, K. Skb1-mediated symmetric dimethylation of histone h4r3 controls flowering time in arabidopsis. EMBO J. 2007, 26, 1934–1941. [Google Scholar] [CrossRef] [PubMed]
  16. Zhao, Q.; Rank, G.; Tan, Y.T.; Li, H.; Moritz, R.L.; Simpson, R.J.; Cerruti, L.; Curtis, D.J.; Patel, D.J.; Allis, C.D.; et al. PRMT5-mediated methylation of histone H4R3 recruits DNMT3A, coupling histone and DNA methylation in gene silencing. Nat. Struct. Mol. Biol. 2009, 16, 304–311. [Google Scholar] [CrossRef] [PubMed]
  17. Wang, H.; Huang, Z.Q.; Xia, L.; Feng, Q.; Erdjument-Bromage, H.; Strahl, B.D.; Briggs, S.D.; Allis, C.D.; Wong, J.; Tempst, P.; et al. Methylation of histone H4 at arginine 3 facilitating transcriptional activation by nuclear hormone receptor. Science 2001, 293, 853–857. [Google Scholar] [CrossRef] [PubMed]
  18. Strahl, B.D.; Briggs, S.D.; Brame, C.J.; Caldwell, J.A.; Koh, S.S.; Ma, H.; Cook, R.G.; Shabanowitz, J.; Hunt, D.F.; Stallcup, M.R.; et al. Methylation of histone H4 at arginine 3 occurs in vivo and is mediated by the nuclear receptor coactivator prmt1. Curr. Biol. 2001, 11, 996–1000. [Google Scholar] [CrossRef]
  19. Zhang, X.; Cheng, X. Structure of the predominant protein arginine methyltransferase prmt1 and analysis of its binding to substrate peptides. Structure 2003, 11, 509–520. [Google Scholar] [CrossRef]
  20. Siarheyeva, A.; Senisterra, G.; Allali-Hassani, A.; Dong, A.; Dobrovetsky, E.; Wasney, G.A.; Chau, I.; Marcellus, R.; Hajian, T.; Liu, F.; et al. An allosteric inhibitor of protein arginine methyltransferase 3. Structure 2012, 20, 1425–1435. [Google Scholar] [CrossRef] [PubMed]
  21. Wang, C.; Zhu, Y.; Chen, J.; Li, X.; Peng, J.; Chen, J.; Zou, Y.; Zhang, Z.; Jin, H.; Yang, P.; et al. Crystal structure of arginine methyltransferase 6 from trypanosoma brucei. PLoS ONE 2014, 9, e87267. [Google Scholar] [CrossRef] [PubMed]
  22. Cheng, Y.; Frazier, M.; Lu, F.; Cao, X.; Redinbo, M.R. Crystal structure of the plant epigenetic protein arginine methyltransferase 10. J. Mol. Biol. 2011, 414, 106–122. [Google Scholar] [CrossRef] [PubMed]
  23. Antonysamy, S.; Bonday, Z.; Campbell, R.M.; Doyle, B.; Druzina, Z.; Gheyi, T.; Han, B.; Jungheim, L.N.; Qian, Y.; Rauch, C.; et al. Crystal structure of the human PRMT5:MEP50 complex. Proc. Natl. Acad. Sci. USA 2012, 109, 17960–17965. [Google Scholar] [CrossRef] [PubMed]
  24. Wang, C.; Zhu, Y.; Caceres, T.B.; Liu, L.; Peng, J.; Wang, J.; Chen, J.; Chen, X.; Zhang, Z.; Zuo, X.; et al. Structural determinants for the strict monomethylation activity by trypanosoma brucei protein arginine methyltransferase 7. Structure 2014, 22, 756–768. [Google Scholar] [CrossRef] [PubMed]
  25. Zhang, R.; Li, X.; Liang, Z.; Zhu, K.; Lu, J.; Kong, X.; Ouyang, S.; Li, L.; Zheng, Y.G.; Luo, C. Theoretical insights into catalytic mechanism of protein arginine methyltransferase 1. PLoS ONE 2013, 8, e72424. [Google Scholar] [CrossRef] [PubMed]
  26. Chu, Y.; Li, G.; Guo, H. QM/MM md and free energy simulations of the methylation reactions catalyzed by protein arginine methyltransferase prmt3. Can. J. Chem. 2013, 91, 605–612. [Google Scholar] [CrossRef]
  27. Guo, H.B.; Guo, H. Mechanism of histone methylation catalyzed by protein lysine methyltransferase set7/9 and origin of product specificity. Proc. Natl. Acad. Sci. USA 2007, 104, 8797–8802. [Google Scholar] [CrossRef] [PubMed]
  28. Xu, Q.; Chu, Y.Z.; Guo, H.B.; Smith, J.C.; Guo, H. Energy triplets for writing epigenetic marks: Insights from QM/MM free-energy simulations of protein lysine methyltransferases. Chemistry 2009, 15, 12596–12599. [Google Scholar] [CrossRef] [PubMed]
  29. Chu, Y.; Yao, J.; Guo, H. QM/MM MD and free energy simulations of G9a-like protein (glp) and its mutants: Understanding the factors that determine the product specificity. PLoS ONE 2012, 7, e37674. [Google Scholar] [CrossRef] [PubMed]
  30. Chu, Y.; Xu, Q.; Guo, H. Understanding energetic origins of product specificity of set8 from qm/mm free energy simulations: What causes the stop of methyl addition during histone lysine methylation? J. Chem. Theory Comput. 2010, 6, 1380–1389. [Google Scholar] [CrossRef]
  31. Yao, J.; Xu, Q.; Chen, F.; Guo, H. QM/MM free energy simulations of salicylic acid methyltransferase: Effects of stabilization of ts-like structures on substrate specificity. J. Phys. Chem. B 2010, 115, 389–396. [Google Scholar] [CrossRef] [PubMed]
  32. Yue, Y.; Guo, H. Quantum mechanical/molecular mechanical study of catalytic mechanism and role of key residues in methylation reactions catalyzed by dimethylxanthine methyltransferase in caffeine biosynthesis. J. Chem. Inf. Model. 2014, 54, 593–600. [Google Scholar] [CrossRef] [PubMed]
  33. Cheng, X.D.; Collins, R.E.; Zhang, X. Structural and sequence motifs of protein (histone) methylation enzymes. Annu. Rev. Biophys. Biomol. Struct. 2005, 34, 267–294. [Google Scholar] [CrossRef] [PubMed]
  34. Wang, M.; Fuhrmann, J.; Thompson, P.R. Protein arginine methyltransferase 5 catalyzes substrate dimethylation in a distributive fashion. Biochemistry 2014, 53, 7884–7892. [Google Scholar] [CrossRef] [PubMed]
  35. Brooks, B.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A program for macro-molecular energy, minimization, and dynamics calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar] [CrossRef]
  36. Jorgensen, W.; Chandrasekhar, J.; Madura, J.; Impey, R.; Klein, M. Comparison of simple potential functions for simulating liquid water. J. Chem. Phys. 1983, 79, 926–935. [Google Scholar] [CrossRef]
  37. Brooks, C.L.; Brunger, A.; Karplus, M. Active site dynamics in protein molecules: A stochastic boundary molecular-dynamics approach. Biopolymers 1985, 24, 843–865. [Google Scholar] [CrossRef] [PubMed]
  38. MacKerell, A.D.; Bashford, D.; Bellott, M.; Dunbrack, R.L.; Evanseck, J.D.; Field, M.J.; Fischer, S.; Gao, J.; Guo, H.; Ha, S.; et al. All-atom empirical potential for molecular modeling and dynamics studies of proteins. J. Phys. Chem. B 1998, 102, 3586–3616. [Google Scholar] [CrossRef] [PubMed]
  39. Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, T.; Suhai, S.; Seifert, G. Self-consistent-charge density functional tight-binding method for simulations of complex materials properties. Phys. Rev. B 1998, 58, 7260–7268. [Google Scholar] [CrossRef]
  40. Cui, Q.; Elstner, M.; Kaxiras, E.; Frauenheim, T.; Karplus, M. A QM/MM Implementation of the Self-Consistent Charge Density Functional Tight Binding (SCC-DFTB) Method. J. Phys. Chem. B 2001, 105, 569–585. [Google Scholar] [CrossRef]
  41. Field, M.J.; Bash, P.A.; Karplus, M. A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations. J. Comput. Chem. 1990, 11, 700–733. [Google Scholar] [CrossRef]
  42. Torrie, G.M.; Valleau, J.P. Monte-Carlo free energy estimates using non-Boltzmann sampling. Application to the subcritical Lennard Jones fluid. Chem. Phys. Lett. 1974, 28, 578–581. [Google Scholar] [CrossRef]
  43. Kumar, S.; Bouzida, D.; Swendsen, R.H.; Kollman, P.A.; Rosenberg, J.M. The weighted histogram analysis method for free energy calculations on biomolecules. I. The method. J. Comput. Chem. 1992, 13, 1011–1021. [Google Scholar] [CrossRef]
  • Sample Availability: Not apply.

Share and Cite

MDPI and ACS Style

Yue, Y.; Chu, Y.; Guo, H. Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations. Molecules 2015, 20, 10032-10046. https://doi.org/10.3390/molecules200610032

AMA Style

Yue Y, Chu Y, Guo H. Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations. Molecules. 2015; 20(6):10032-10046. https://doi.org/10.3390/molecules200610032

Chicago/Turabian Style

Yue, Yufei, Yuzhuo Chu, and Hong Guo. 2015. "Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations" Molecules 20, no. 6: 10032-10046. https://doi.org/10.3390/molecules200610032

APA Style

Yue, Y., Chu, Y., & Guo, H. (2015). Computational Study of Symmetric Methylation on Histone Arginine Catalyzed by Protein Arginine Methyltransferase PRMT5 through QM/MM MD and Free Energy Simulations. Molecules, 20(6), 10032-10046. https://doi.org/10.3390/molecules200610032

Article Metrics

Back to TopTop