Understanding the Phosphorylation Mechanism by Using Quantum Chemical Calculations and Molecular Dynamics Simulations

Phosphorylation is one of the most frequent post-translational modifications on proteins. It regulates many cellular processes by modulation of phosphorylation on protein structure and dynamics. However, the mechanism of phosphorylation-induced conformational changes of proteins is still poorly understood. Here, we report a computational study of three representative groups of tyrosine in ADP-ribosylhydrolase 1, serine in BTG2, and serine in Sp100C by using six molecular dynamics (MD) simulations and quantum chemical calculations. Added phosphorylation was found to disrupt hydrogen bond, and increase new weak interactions (hydrogen bond and hydrophobic interaction) during MD simulations, leading to conformational changes. Quantum chemical calculations further indicate that the phosphorylation on tyrosine, threonine, and serine could decrease the optical band gap energy (<i>E</i>_gap), which can trigger electronic transitions to form or disrupt interactions easily. Our results provide an atomic and electronic description of how phosphorylation facilitates conformational and dynamic changes in proteins, which may be useful for studying protein function and protein design

Data_Sheet_4_Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study.ZIP

<p>Thermus thermophiles Argonaute (TtAgo) is a complex, which is consisted of 5′-phosphorylated guide DNA and a series of target DNA with catalytic activities at high temperatures. To understand why high temperatures are needed for the catalytic activities, three molecular dynamics simulations and binding free energy calculations at 310, 324, and 338K were performed for the TtAgo-DNA complex to explore the conformational changes between 16-mer guide DNA/15-mer target DNA and TtAgo at different temperatures. The simulation results indicate that a collapse of a small β-strand (residues 507–509) at 310 K caused Glu512 to move away from the catalytic residues Asp546 and Asp478, resulting in a decrease in catalytic activity, which was not observed in the simulations at 324 and 338 K. The nucleic acid binding channel became enlarged at 324 and 338K, thereby facilitating the DNA to slide in. Binding free energy calculations and hydrogen bond occupancy indicated that the interaction between TtAgo and the DNA was more stable at 324K and 338K than at 310 K. The DNA binding pocket residues Lys575 and Asn590 became less solvent accessible at 324 and 338K than at 310 K to influence hydrophilic interaction with DNA. Our simulation studies shed some light on the mechanism of TtAgo and explained why a high temperature was needed by TtAgo during gene editing of CRISPR.</p

Data_Sheet_3_Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study.ZIP

<p>Thermus thermophiles Argonaute (TtAgo) is a complex, which is consisted of 5′-phosphorylated guide DNA and a series of target DNA with catalytic activities at high temperatures. To understand why high temperatures are needed for the catalytic activities, three molecular dynamics simulations and binding free energy calculations at 310, 324, and 338K were performed for the TtAgo-DNA complex to explore the conformational changes between 16-mer guide DNA/15-mer target DNA and TtAgo at different temperatures. The simulation results indicate that a collapse of a small β-strand (residues 507–509) at 310 K caused Glu512 to move away from the catalytic residues Asp546 and Asp478, resulting in a decrease in catalytic activity, which was not observed in the simulations at 324 and 338 K. The nucleic acid binding channel became enlarged at 324 and 338K, thereby facilitating the DNA to slide in. Binding free energy calculations and hydrogen bond occupancy indicated that the interaction between TtAgo and the DNA was more stable at 324K and 338K than at 310 K. The DNA binding pocket residues Lys575 and Asn590 became less solvent accessible at 324 and 338K than at 310 K to influence hydrophilic interaction with DNA. Our simulation studies shed some light on the mechanism of TtAgo and explained why a high temperature was needed by TtAgo during gene editing of CRISPR.</p

Image_1_Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study.PDF

<p>Thermus thermophiles Argonaute (TtAgo) is a complex, which is consisted of 5′-phosphorylated guide DNA and a series of target DNA with catalytic activities at high temperatures. To understand why high temperatures are needed for the catalytic activities, three molecular dynamics simulations and binding free energy calculations at 310, 324, and 338K were performed for the TtAgo-DNA complex to explore the conformational changes between 16-mer guide DNA/15-mer target DNA and TtAgo at different temperatures. The simulation results indicate that a collapse of a small β-strand (residues 507–509) at 310 K caused Glu512 to move away from the catalytic residues Asp546 and Asp478, resulting in a decrease in catalytic activity, which was not observed in the simulations at 324 and 338 K. The nucleic acid binding channel became enlarged at 324 and 338K, thereby facilitating the DNA to slide in. Binding free energy calculations and hydrogen bond occupancy indicated that the interaction between TtAgo and the DNA was more stable at 324K and 338K than at 310 K. The DNA binding pocket residues Lys575 and Asn590 became less solvent accessible at 324 and 338K than at 310 K to influence hydrophilic interaction with DNA. Our simulation studies shed some light on the mechanism of TtAgo and explained why a high temperature was needed by TtAgo during gene editing of CRISPR.</p

Data_Sheet_2_Why Is a High Temperature Needed by Thermus thermophilus Argonaute During mRNA Silencing: A Theoretical Study.ZIP

<p>Thermus thermophiles Argonaute (TtAgo) is a complex, which is consisted of 5′-phosphorylated guide DNA and a series of target DNA with catalytic activities at high temperatures. To understand why high temperatures are needed for the catalytic activities, three molecular dynamics simulations and binding free energy calculations at 310, 324, and 338K were performed for the TtAgo-DNA complex to explore the conformational changes between 16-mer guide DNA/15-mer target DNA and TtAgo at different temperatures. The simulation results indicate that a collapse of a small β-strand (residues 507–509) at 310 K caused Glu512 to move away from the catalytic residues Asp546 and Asp478, resulting in a decrease in catalytic activity, which was not observed in the simulations at 324 and 338 K. The nucleic acid binding channel became enlarged at 324 and 338K, thereby facilitating the DNA to slide in. Binding free energy calculations and hydrogen bond occupancy indicated that the interaction between TtAgo and the DNA was more stable at 324K and 338K than at 310 K. The DNA binding pocket residues Lys575 and Asn590 became less solvent accessible at 324 and 338K than at 310 K to influence hydrophilic interaction with DNA. Our simulation studies shed some light on the mechanism of TtAgo and explained why a high temperature was needed by TtAgo during gene editing of CRISPR.</p

Exploring stereochemical specificity of phosphotriesterase by MM-PBSA and MM-GBSA calculation and steered molecular dynamics simulation

<p>Wild-type phosphotriesterase (PTE) prefers the <i>S</i>_P-enantiomers over the corresponding <i>R</i>_P-enantiomers by factors ranging from 10 to 90. To satisfy the binding modes of the PTE of <i>S</i>_P- and <i>R</i>_P-enantiomers, all-atom molecular dynamics simulations were carried out on two paraoxon <i>S</i>_P and <i>R</i>_P derivatives, namely, <i>S</i>_p-1 and <i>R</i>_p-1. Molecular mechanics Poisson–Boltzmann surface area and molecular mechanics generalized Born surface area (MM-PBSA and MM-GBSA) calculations indicated that His230 in <i>S</i>_p-1-PTE had a closer interaction with the substrate than that in <i>R</i>_p-1-PTE and that such interaction increased the catalytic efficiency of PTE for <i>S</i>_p-1. The steered molecular dynamics simulation indicated that, compared with <i>S</i>_p-1, <i>R</i>_p-1 in the unbinding (binding) may hinder some residue displacement, thus requiring more effort to escape the binding pocket of PTE. In addition, Trp131, Phe306, and Tyr309 are deemed important residues for the <i>S</i>_p-1 unbinding pathway via PTE, whereas Tyr309 alone is considered an important residue for the <i>R</i>_p-1 unbinding pathway. These results demonstrate the possibility of dramatically altering the stereoselectivity and overall reactivity of the native enzyme toward chiral substrates by modifying specific residues located within the active site of PTE.</p