Competition between Intramolecular Hydrogen Bonds and Solvation in Phosphorylated Peptides: Simulations with Explicit and Implicit Solvent

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared with published findings. We simulated short (<20 residues) and long (>80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multiphosphorylated IDRs. Our results support and expand on previous observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

A curated rotamer library for common post-translational modifications of proteins

MOTIVATION

Sidechain rotamer libraries of the common amino acids of a protein are useful for folded protein structure determination and for generating ensembles of intrinsically disordered proteins (IDPs). However, much of protein function is modulated beyond the translated sequence through the introduction of post-translational modifications (PTMs).

RESULTS

In this work, we have provided a curated set of side chain rotamers for the most common PTMs derived from the RCSB PDB database, including phosphorylated, methylated, and acetylated sidechains. Our rotamer libraries improve upon existing methods such as SIDEpro, Rosetta, and AlphaFold3 in predicting the experimental structures for PTMs in folded proteins. In addition, we showcase our PTM libraries in full use by generating ensembles with the Monte Carlo Side Chain Entropy (MCSCE) for folded proteins, and combining MCSCE with the Local Disordered Region Sampling algorithms within IDPConformerGenerator for proteins with intrinsically disordered regions.

AVAILABILITY AND IMPLEMENTATION

The codes for dihedral angle computations and library creation are available at https://github.com/THGLab/ptm_sc.git.

Effects of Phosphorylation on Protein Backbone Dynamics and Conformational Preferences

Phosphorylations are the most common and extensively studied post-translational modification (PTM) of proteins in eukaryotes. They constitute a major regulatory mechanism, modulating protein function, protein-protein interactions, as well as subcellular localization. Phosphorylation sites are preferably located in intrinsically disordered regions and have been shown to trigger structural rearrangements and order-to-disorder transitions. They can therefore have a significant effect on protein backbone dynamics or conformation, but only sparse experimental data are available. To obtain a more general description of how and when phosphorylations have a significant effect on protein behavior, molecular dynamics (MD) currently provides the only suitable framework to study these effects at a large scale in atomistic detail. This study develops a systematic MD simulation framework to explore the influence of phosphorylations on the local backbone dynamics and conformational propensities of proteins. Through a series of glycine-backbone peptides, we studied the effects of amino acid residues including the three most common phosphorylations (Ser, Thr, and Tyr), on local backbone dynamics and conformational propensities. We further extended our study to investigate the interactions of all such residues between position i to positions i + 1, i + 2, i + 3, and i + 4 in such peptides. The final data set comprises structural ensembles for 3393 sequences with more than 1 μs of sampling for each ensemble. To validate the relevance of the results, the structural and conformational properties extracted from the MD simulations are compared to NMR data from the Biological Magnetic Resonance Data Bank. The systematic nature of this study enables the projection of the gained knowledge onto any phosphorylation site in the proteome and provides a general framework for the study of further PTMs. The full data set is publicly available, as a training and reference set.

Phosphorylation of disordered proteins tunes local and global intramolecular interactions

Protein post-translational modifications, such as phosphorylation, are important regulatory signals for diverse cellular functions. In particular, intrinsically disordered protein regions (IDRs) are subject to phosphorylation as a means to modulate their interactions and functions. Toward understanding the relationship between phosphorylation in IDRs and specific functional outcomes, we must consider how phosphorylation affects the IDR conformational ensemble. Various experimental techniques are suited to interrogate the features of IDR ensembles; molecular simulations can provide complementary insights and even illuminate ensemble features that may be experimentally inaccessible. Therefore, we sought to expand the tools available to study phosphorylated IDRs by all-atom Monte Carlo simulations. To this end, we implemented parameters for phosphoserine (pSer) and phosphothreonine (pThr) into the OPLS version of the continuum solvent model, ABSINTH, and assessed their performance in all-atom simulations compared to published findings. We simulated short (< 20 residues) and long (> 80 residues) phospho-IDRs that, collectively, survey both local and global phosphorylation-induced changes to the ensemble. Our simulations of four well-studied phospho-IDRs show near-quantitative agreement with published findings for these systems via metrics including changes to radius of gyration, transient helicity, and persistence length. We also leveraged the inherent advantage of sequence control in molecular simulations to explore the conformational effects of diverse combinations of phospho-sites in two multi-phosphorylated IDRs. Our results support and expand on prior observations that connect phosphorylation to changes in the IDR conformational ensemble. Herein, we describe phosphorylation as a means to alter sequence chemistry, net charge and charge patterning, and intramolecular interactions, which can collectively modulate the local and global IDR ensemble features.

An Inherent Difference between Serine and Threonine Phosphorylation: Phosphothreonine Strongly Prefers a Highly Ordered, Compact, Cyclic Conformation

Phosphorylation and dephosphorylation of proteins by kinases and phosphatases are central to cellular responses and function. The structural effects of serine and threonine phosphorylation were examined in peptides and in proteins, by circular dichroism, NMR spectroscopy, bioinformatics analysis of the PDB, small-molecule X-ray crystallography, and computational investigations. Phosphorylation of both serine and threonine residues induces substantial conformational restriction in their physiologically more important dianionic forms. Threonine exhibits a particularly strong disorder-to-order transition upon phosphorylation, with dianionic phosphothreonine preferentially adopting a cyclic conformation with restricted ϕ (ϕ ∼ -60°) stabilized by three noncovalent interactions: a strong intraresidue phosphate-amide hydrogen bond, an n → π* interaction between consecutive carbonyls, and an n → σ* interaction between the phosphate Oγ lone pair and the antibonding orbital of C-Hβ that restricts the χ2 side-chain conformation. Proline is unique among the canonical amino acids for its covalent cyclization on the backbone. Phosphothreonine can mimic proline's backbone cyclization via noncovalent interactions. The preferred torsions of dianionic phosphothreonine are ϕ,ψ = polyproline II helix > α-helix (ϕ ∼ -60°); χ1 = g-; χ2 ∼ +115° (eclipsed C-H/O-P bonds). This structural signature is observed in diverse proteins, including in the activation loops of protein kinases and in protein-protein interactions. In total, these results suggest a structural basis for the differential use and evolution of threonine versus serine phosphorylation sites in proteins, with serine phosphorylation typically inducing smaller, rheostat-like changes, versus threonine phosphorylation promoting larger, step function-like switches, in proteins.

In-Silico Selection of Aptamer Targeting SARS-CoV-2 Spike Protein

Aptamers are single-stranded, short DNA or RNA oligonucleotides that can specifically bind to various target molecules. To diagnose the infected cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) in time, numerous conventional methods are applied for viral detection via the amplification and quantification of DNA or antibodies specific to antigens on the virus. Herein, we generated a large number of mutated aptamer sequences, derived from a known sequence of receptor-binding domain (RBD)-1C aptamer, specific to the RBD of SARS-CoV-2 spike protein (S protein). Structural similarity, molecular docking, and molecular dynamics (MD) were utilized to screen aptamers and characterize the detailed interactions between the selected aptamers and the S protein. We identified two mutated aptamers, namely, RBD-1CM1 and RBD-1CM2, which presented better docking results against the S protein compared with the RBD-1C aptamer. Through the MD simulation, we further confirmed that the RBD-1CM1 aptamer can form the most stable complex with the S protein based on the number of hydrogen bonds formed between the two biomolecules. Based on the experimental data of quartz crystal microbalance (QCM), the RBD-1CM1 aptamer could produce larger signals in mass change and exhibit an improved binding affinity to the S protein. Therefore, the RBD-1CM1 aptamer, which was selected from 1431 mutants, was the best potential candidate for the detection of SARS-CoV-2. The RBD-1CM1 aptamer can be an alternative biological element for the development of SARS-CoV-2 diagnostic testing.

AMBER and CHARMM Force Fields Inconsistently Portray the Microscopic Details of Phosphorylation

Phosphorylation of serine, threonine, and tyrosine is one of the most frequently occurring and crucial post-translational modifications of proteins often associated with important structural and functional changes. We investigated the direct effect of phosphorylation on the intrinsic conformational preferences of amino acids as a potential trigger of larger structural events. We conducted a comparative study of force fields on terminally capped amino acids (dipeptides) as the simplest model for phosphorylation. Our bias-exchange metadynamics simulations revealed that all model dipeptides sampled a great heterogeneity of ensembles affected by introduction of mono- and dianionic phosphate groups. However, the detected changes in populations of backbone conformers and side-chain rotamers did not reveal a strong discriminatory shift in preferences, as could be anticipated for the bulky, charged phosphate group. Furthermore, the AMBER and CHARMM force fields provided inconsistent populations of individual conformers as well as net structural trends upon phosphorylation. Detailed analysis of ensembles revealed competition between hydration and formation of internal hydrogen bonds involving amide hydrogens and the phosphate group. The observed difference in hydration free energy and potential for hydrogen bonding in individual force fields could be attributed to the different partial atomic charges used in each force field and, hence, the different parametrization strategies. Nevertheless, conformational propensities and net structural changes upon phosphorylation are difficult to extract from experimental measurements, and existing experimental data provide limited guidance for force field assessment and further development.

Phosphorylation Increases Persistence Length and End-to-End Distance of a Segment of Tau Protein

Intrinsically disordered regions of proteins, which lack unique tertiary structure under physiological conditions, are enriched in phosphorylation sites and in significant local bias toward the polyproline II conformation. The overrepresented coincidence of this posttranslational regulatory signal and local conformational bias within unstructured regions raises a question: can phosphorylation serve to manipulate the conformational preferences of a disordered protein? In this study, we use time-resolved fluorescence resonance energy transfer and a, to our knowledge, novel data analysis method to directly measure the end-to-end distance distribution of a phosphorylatable peptide derived from the human microtubule associated protein tau. Our results show that phosphorylation at threonine or serine extends the end-to-end distance and increases the effective persistence length of the tested model peptides. Unexpectedly, the extension is independent of salt concentration, suggestive of a nonelectrostatic origin. The phosphorylation extension and stiffening effect provides a peptide-scale physical interpretation for the posttranslational regulation of the highly abundant protein-protein interactions found in disordered proteins, as well as a potential insight into the regulatory mechanism of the tau protein's microtubule binding activity.

Monte Carlo Sampling with Hierarchical Move Sets: POSH Monte Carlo

We present a new Monte Carlo method for sampling rugged energy landscapes that allows for efficient transitions across sparsely distributed local basins. The trial move consists of two steps. The first step is a large initial trial move, and the second step is a Monte Carlo trajectory generated using smaller trial moves. To maintain detailed balance, a reverse transition probability is estimated along a path that differs from the forward path. Since the forward and reverse transitions are different, we name the algorithm POSH (port out, starboard home) Monte Carlo. The process obeys detailed balance to the extent that the transition probabilities are correctly estimated. There is an optimal range of performance for a given energy landscape, which depends on how sparsely the low energy states of the system are distributed. For simple model systems, adequate precision is obtained over a large range of inner steps settings. Side chain sampling of residues in the binding region of progesterone antibody 1dba are studied, and show that significant improvement over a comparable standard protocol can be obtained using POSH sampling. To compare with experimental data, the phosphopeptide Ace-Gly-Ser-pSer-Ser-Nma is also studied, and the resulting NMR observables compare well with experiment. For the biomolecular systems studied, we show that POSH sampling generates precise distributions using the number of inner steps set up to 20.

QM/MM MD simulations of La(iii)–phosphopeptide complexes

Hybrid quantum mechanical/molecular mechanical simulations have been used to study the structural and dynamical properties of a La(iii)–phosphopeptide complex.

Role of electrostatic interactions in binding of peptides and intrinsically disordered proteins to their folded targets. 1. NMR and MD characterization of the complex between the c-Crk N-SH3 domain and the peptide Sos

Intrinsically disordered proteins (IDPs) often rely on electrostatic interactions to bind their structured targets. To obtain insight into the mechanism of formation of the electrostatic encounter complex, we investigated the binding of the peptide Sos (PPPVPPRRRR), which serves as a minimal model for an IDP, to the c-Crk N-terminal SH3 domain. Initially, we measured ¹⁵N relaxation rates at two magnetic field strengths and determined the binding shifts for the complex of Sos with wild-type SH3. We have also recorded a 3 μs molecular dynamics (MD) trajectory of this complex using the Amber ff99SB*-ILDN force field. The comparison of the experimental and simulated data shows that MD simulation consistently overestimates the strength of salt bridge interactions at the binding interface. The series of simulations using other advanced force fields also failed to produce any satisfactory results. To address this issue, we have devised an empirical correction to the Amber ff99SB*-ILDN force field whereby the Lennard-Jones equilibrium distance for the nitrogen-oxygen pair across the Arg-to-Asp and Arg-to-Glu salt bridges has been increased by 3%. Implementing this correction resulted in a good agreement between the simulations and the experiment. Adjusting the strength of salt bridge interactions removed a certain amount of strain contained in the original MD model, thus improving the binding of the hydrophobic N-terminal portion of the peptide. The arginine-rich C-terminal portion of the peptide, freed from the effect of the overstabilized salt bridges, was found to interconvert more rapidly between its multiple conformational states. The modified MD protocol has also been successfully used to simulate the entire binding process. In doing so, the peptide was initially placed high above the protein surface. It then arrived at the correct bound pose within ∼2 Å of the crystallographic coordinates. This simulation allowed us to analyze the details of the dynamic binding intermediate, i.e., the electrostatic encounter complex. However, an experimental characterization of this transient, weakly populated state remains out of reach. To overcome this problem, we designed the double mutant of c-Crk N-SH3 in which mutations Y186L and W169F abrogate tight Sos binding and shift the equilibrium toward the intermediate state resembling the electrostatic encounter complex. The results of the combined NMR and MD study of this engineered system will be reported in the next part of this paper.

OGlcNAcylation and phosphorylation have similar structural effects in α-helices: post-translational modifications as inducible start and stop signals in α-helices, with greater structural effects on threonine modification

OGlcNAcylation and phosphorylation are the major competing intracellular post-translational modifications of serine and threonine residues. The structural effects of both post-translational modifications on serine and threonine were examined within Baldwin model α-helical peptides (Ac-AKAAAAKAAAAKAAGY-NH₂ or Ac-YGAKAAAAKAAAAKAA-NH₂). At the N-terminus of an α-helix, both phosphorylation and OGlcNAcylation stabilized the α-helix relative to the free hydroxyls, with a larger induced structure for phosphorylation than for OGlcNAcylation, for the dianionic phosphate than for the monoanionic phosphate, and for modifications on threonine than for modifications on serine. Both phosphoserine and phosphothreonine resulted in peptides more α-helical than alanine at the N-terminus, with dianionic phosphothreonine the most α-helix-stabilizing residue here. In contrast, in the interior of the α-helix, both post-translational modifications were destabilizing with respect to the α-helix, with the greatest destabilization seen for threonine OGlcNAcylation at residue 5 and threonine phosphorylation at residue 10, with peptides containing either post-translational modification existing as random coils. At the C-terminus, both OGlcNAcylation and phosphorylation were destabilizing with respect to the α-helix, though the induced structural changes were less than in the interior of the α-helix. In general, the structural effects of modifications on threonine were greater than the effects on serine, because of both the lower α-helical propensity of Thr and the more defined induced structures upon modification of threonine than serine, suggesting threonine residues are particularly important loci for structural effects of post-translational modifications. The effects of serine and threonine post-translational modifications are analogous to the effects of proline on α-helices, with the effects of phosphothreonine being greater than those of proline throughout the α-helix. These results provide a basis for understanding the context-dependent structural effects of these competing protein post-translational modifications.

OGlcNAcylation and phosphorylation have opposing structural effects in tau: phosphothreonine induces particular conformational order

Phosphorylation and OGlcNAcylation are dynamic intracellular protein post-translational modifications that frequently are alternatively observed on the same serine and threonine residues. Phosphorylation and OGlcNAcylation commonly occur in natively disordered regions of proteins, and often have opposing functional effects. In the microtubule-associated protein tau, hyperphosphorylation is associated with protein misfolding and aggregation as the neurofibrillary tangles of Alzheimer's disease, whereas OGlcNAcylation stabilizes the soluble form of tau. A series of peptides derived from the proline-rich domain (residues 174-251) of tau was synthesized, with free Ser/Thr hydroxyls, phosphorylated Ser/Thr (pSer/pThr), OGlcNAcylated Ser/Thr, and diethylphosphorylated Ser/Thr. Phosphorylation and OGlcNAcylation were found by CD and NMR to have opposing structural effects on polyproline helix (PPII) formation, with phosphorylation favoring PPII, OGlcNAcylation opposing PPII, and the free hydroxyls intermediate in structure, and with phosphorylation structural effects greater than OGlcNAcylation. For tau196-209, phosphorylation and OGlcNAcylation had similar structural effects, opposing a nascent α-helix. Phosphomimic Glu exhibited PPII-favoring structural effects. Structural changes due to Thr phosphorylation were greater than those of Ser phosphorylation or Glu, with particular conformational restriction as the dianion, with mean (3)JαN = 3.5 Hz (pThr) versus 5.4 Hz (pSer), compared to 7.2, 6.8, and 6.2 Hz for Thr, Ser, and Glu, respectively, values that correlate with the backbone torsion angle ϕ. Dianionic phosphothreonine induced strong phosphothreonine amide protection and downfield amide chemical shifts (δmean = 9.63 ppm), consistent with formation of a stable phosphate-amide hydrogen bond. These data suggest potentially greater structural importance of threonine phosphorylation than serine phosphorylation due to larger induced structural effects.

The natively disordered loop of Bcl-2 undergoes phosphorylation-dependent conformational change and interacts with Pin1

Bcl-2 plays a central role in the regulation of apoptosis. Structural studies of Bcl-2 revealed the presence of a flexible and natively disordered loop that bridges the Bcl-2 homology motifs, BH3 and BH4. This loop is phosphorylated on multiple sites in response to a variety of external stimuli, including the microtubule-targeting drugs, paclitaxel and colchicine. Currently, the underlying molecular mechanism of Bcl-2 phosphorylation and its biological significance remain elusive. In this study, we investigated the molecular characteristics of this anti-apoptotic protein. To this end, we generated synthetic peptides derived from the Bcl-2 loop, and multiple Bcl-2 loop truncation mutants that include the phosphorylation sites. Our results demonstrate that S87 in the flexible loop of Bcl-2 is the primary phosphorylation site for JNK and ERK2, suggesting some sequence or structural specificity for the phosphorylation by these kinases. Our NMR studies and molecular dynamics simulation studies support indicate that phosphorylation of S87 induces a conformational change in the peptide. Finally, we show that the phosphorylated peptides of the Bcl-2 loop can bind Pin1, further substantiating the phosphorylation-mediated conformation change of Bcl-2.

Multisite phosphorylation disrupts arginine-glutamate salt bridge networks required for binding of cytoplasmic linker-associated protein 2 (CLASP2) to end-binding protein 1 (EB1)

A group of diverse proteins reversibly binds to growing microtubule plus ends through interactions with end-binding proteins (EBs). These +TIPs control microtubule dynamics and microtubule interactions with other intracellular structures. Here, we use cytoplasmic linker-associated protein 2 (CLASP2) binding to EB1 to determine how multisite phosphorylation regulates interactions with EB1. The central, intrinsically disordered region of vertebrate CLASP proteins contains two SXIP EB1 binding motifs that are required for EB1-mediated plus-end-tracking in vitro. In cells, both EB1 binding motifs can be functional, but most of the binding free energy results from nearby electrostatic interactions. By employing molecular dynamics simulations of the EB1 interaction with a minimal CLASP2 plus-end-tracking module, we find that conserved arginine residues in CLASP2 form extensive hydrogen-bond networks with glutamate residues predominantly in the unstructured, acidic C-terminal tail of EB1. Multisite phosphorylation of glycogen synthase kinase 3 (GSK3) sites near the EB1 binding motifs disrupts this electrostatic "molecular Velcro." Molecular dynamics simulations and (31)P NMR spectroscopy indicate that phosphorylated serines participate in intramolecular interactions with and sequester arginine residues required for EB1 binding. Multisite phosphorylation of these GSK3 motifs requires priming phosphorylation by interphase or mitotic cyclin-dependent kinases (CDKs), and we find that CDK- and GSK3-dependent phosphorylation completely disrupts CLASP2 microtubule plus-end-tracking in mitosis.

Phosphorylation alters backbone conformational preferences of serine and threonine peptides

Despite the notion that a control of protein function by phosphorylation works mainly by inducing its conformational changes, the phosphorylation effects on even small peptide conformation have not been fully understood yet. To study its possible effects on serine and threonine peptide conformations, we recently carried out pH- and temperature-dependent circular dichroism (CD) as well as (1)H NMR studies of the phosphorylated serine and threonine peptides and compared them with their unphosphorylated analogs. In the present article, by performing the self-consistent singular value decomposition analysis of the temperature-dependent CD spectra and by analyzing the (3)J(H(N),H(α)) coupling constants extracted from the NMR spectra, the populations of the polyproline II (PPII) and β-strand conformers of the phosphorylated Ser and Thr peptides are determined. As temperature is increased, the β-strand populations of both phosphorylated serine and threonine peptides increase. However, the dependences of PPII/β-strand population ratio on pH are different for these two cases. The phosphorylation of the serine peptide enhances the PPII propensity, whereas that of the threonine peptide has the opposite effect. This suggests that the serine and threonine phosphorylations can alter the backbone conformational propensity via direct but selective intramolecular hydrogen-bonding interactions with the peptide N--H groups. This clearly indicates that the phosphoryl group actively participates in modulating the peptide backbone conformations.

On the role of electrostatics in protein-protein interactions

The role of electrostatics in protein-protein interactions and binding is reviewed in this paper. A brief outline of the computational modeling, in the framework of continuum electrostatics, is presented and the basic electrostatic effects occurring upon the formation of the complex are discussed. The effect of the salt concentration and pH of the water phase on protein-protein binding free energy is demonstrated which indicates that the increase of the salt concentration tends to weaken the binding, an observation that is attributed to the optimization of the charge-charge interactions across the interface. It is pointed out that the pH-optimum (pH of optimal binding affinity) varies among the protein-protein complexes, and perhaps is a result of their adaptation to particular subcellular compartments. The similarities and differences between hetero- and homo-complexes are outlined and discussed with respect to the binding mode and charge complementarity.

Assessing protein loop flexibility by hierarchical Monte Carlo sampling

Loop flexibility is often crucial to protein biological function in solution. We report a new Monte Carlo method for generating conformational ensembles for protein loops and cyclic peptides. The approach incorporates the triaxial loop closure method which addresses the inverse kinematic problem for generating backbone move sets that do not break the loop. Sidechains are sampled together with the backbone in a hierarchical way, making it possible to make large moves that cross energy barriers. As an initial application, we apply the method to the flexible loop in triosephosphate isomerase that caps the active site, and demonstrate that the resulting loop ensembles agree well with key observations from previous structural studies. We also demonstrate, with 3 other test cases, the ability to distinguish relatively flexible and rigid loops within the same protein.

Structural basis for the complete loss of GSK3beta catalytic activity due to R96 mutation investigated by molecular dynamics study

Many Ser/Thr protein kinases, to be fully activated, are obligated to introduce a phospho-Ser/Thr in their activation loop. Presently, the similarity of activation loop between two crystal complexes, i.e. glycogen synthase kinase 3beta (GSK3beta)-AMPNP and GSK3beta-sulfate ion complex, indicates that the activation segment of GSK3beta is preformed requiring neither a phosphorylation event nor conformational changes. GSK3beta, when participated in glycogen synthesis and Wnt signaling pathways, possesses a unique feature with the preference of such substrate with a priming phosphate. Experimental mutagenesis proved that the residue arginine at amino acid 96 mutations to lysine (R96K) or alanine (R96A) selectively abolish activity on the substrates involved in glycogen synthesis signaling pathway. Based on two solved crystal structures, wild type (WT) and two mutants (R96K and R96A) GSK3beta-ATP-phospho-Serine (pSer) complexes were modeled. Molecular dynamics simulations and energy analysis were employed to investigate the effect of pSer involvement on the GSK3beta structure in WT, and the mechanisms of GSK3beta deactivation due to R96K and R96A mutations. The results indicate that the introduction of pSer to WT GSK3beta generates a slight lobe closure on GSK3beta without any remarkable changes, which may illuminate the experimental conclusion, whereas the conformations of GSK3beta and ATP undergo significant changes in two mutants. As to GSK3beta, the affected positions distribute over activation loop, alpha-helix, and glycine-rich loop. Based on coupling among the mentioned positions, the allosteric mechanisms for distorted ATP were proposed. Energy decomposition on the residues of activation loop identified the important residues Arg96 and Arg180 in anchoring the phosphate group.

Characterization of phosphorylated peptides using traveling wave-based and drift cell ion mobility mass spectrometry

Phosphorylation is one the most studied and important post translational modifications. Nano electrospray mass spectrometry coupled with traveling wave (T-Wave)-based ion mobility has been used to filter for phosphorylated peptides in tryptic protein digests. T-Wave parameters have been optimized to maximize the separation between phosphorylated and non-phosphorylated peptides. A method to calibrate the T-Wave device, to provide estimates of collision cross sections, is presented, and these estimates are in excellent agreement with values obtained on drift cell instrumentation. Phosphorylated peptides have smaller cross sections which enables their separation from non-phosphorylated peptides of the same m/z. Post-mobility fragmentation is used to obtain the primary sequence for peptides of interest. This approach is shown to have potential as an additional screen for phosphorylated peptides, where up to 40% of observed peptides can be eliminated from the study.

Elastic deformation and failure in protein filament bundles: Atomistic simulations and coarse-grained modeling

The synthetic peptide RAD16-II has shown promise in tissue engineering and drug delivery. It has been studied as a vehicle for cell delivery and controlled release of IGF-1 to repair infarcted cardiac tissue, and as a scaffold to promote capillary formation for an in vitro model of angiogenesis. The structure of RAD16-II is hierarchical, with monomers forming long beta-sheets that pair together to form filaments; filaments form bundles approximately 30-60 nm in diameter; branching networks of filament bundles form macroscopic gels. We investigate the mechanics of shearing between the two beta-sheets constituting one filament, and between cohered filaments of RAD16-II. This shear loading is found in filament bundle bending or in tensile loading of fibers composed of partial-length filaments. Molecular dynamics simulations show that time to failure is a stochastic function of applied shear stress, and that for a given loading time behavior is elastic for sufficiently small shear loads. We propose a coarse-grained model based on Langevin dynamics that matches molecular dynamics results and facilities extending simulations in space and time. The model treats a filament as an elastic string of particles, each having potential energy that is a periodic function of its position relative to the neighboring filament. With insight from these simulations, we discuss strategies for strengthening RAD16-II and similar materials.

Phosphorylation effect on the GSSS peptide conformation in water: Infrared, vibrational circular dichroism, and circular dichroism experiments and comparisons with molecular dynamics simulations

The phosphorylation effect on the small peptide conformation in water has not been clearly understood yet, despite the widely acknowledged notion that control of protein activity by phosphorylation works mainly by inducing conformational change. To elucidate the detailed mechanism, we performed infrared (IR) absorption and vibrational and electronic circular dichroism studies of both unphosphorylated and phosphorylated tetrapeptides, GSSS 1 and GSSpS 2. The solution structure of the tetrapeptide is found to be little dependent on the presence of the neutral or negatively charged phosphoryl group, and to be a mixture of extended structures including polyproline II (PII) and beta-sheet conformations. The additional band at 1598 cm(-1) in the amide I IR spectrum of the phosphorylated peptide GSSpS at neutral pD appears to be clear spectroscopic evidence for direct intramolecular hydrogen-bonding interaction between the side chain dianionic phosphoryl group and the backbone amide proton. On the basis of amide I IR band analyses, the authors found that the probability of finding the phosphoryl group strongly H bonded to the backbone proton in GSSpS is about 43% at pD 7.0 and 37 degrees C. Such a H-bonding interaction in GSSpS has the biological standard enthalpy and entropy of -15.1 kJ/mol and -51.2 J/K mol, respectively. Comparisons between the experimentally measured IR and VCD spectra and the numerically simulated ones suggested that the currently available force field parameters need to be properly modified. The results in this paper may shed light on an unknown mechanism of controlling the peptide conformation by phosphorylation.

Car−Parrinello Molecular Dynamics Study of Anharmonic Systems: A Mannich Base in Solution

A Car-Parrinello molecular dynamics study was performed for 4,5-dimethyl-2-(N,N-dimethylaminomethyl)phenol, a Mannich base, to investigate the vibrational properties in solution of its intramolecular hydrogen bond. The dynamic behavior of this hydrogen-bonded system was investigated using an explicit solvent model. Addition of a nonpolar solvent permitted inclusion of delicate environmental effects on the strongly anharmonic system which was studied from first principles. Molecular dynamics and a posteriori quantization of the O-H motion were applied to reproduce the vibrational features of the O-H stretching mode. Consistent application of Car-Parrinello dynamics based on the density functional theory with subsequent solution of the vibrational Schrödinger equation for the O-H stretching motion offers an effective method for strongly anharmonic systems, and this is supported by the comparison of the results with experimental spectra. As a further element of the intramolecular hydrogen bond study, the effects of deuteration were taken into account and a successful application of the O-H stretching mode quantization technique to the liquid phase is demonstrated. This provides a valuable computational methodology for investigations incorporating nuclear quantum effects in the liquid phase and enzyme active centers and can be used to investigate numerous systems that are not readily susceptible to experimental analysis.

Influence of Metal Cations on the Intramolecular Hydrogen-Bonding Network and pKa in Phosphorylated Compounds

The binding of the most common metal cations of cytoplasm (Li+, Na+, K+, Mg2+ and Ca2+) to a model molecule having an intramolecular hydrogen-bonding network, myo-inositol-2-monophosphate, was studied using first principles. A strong correlation between the conformation of metal inositol phosphate complexes with the type of metal cation, degree of deprotonation state, and the surrounding environment has been observed. On the basis of the hydrogen-bonding network analysis of the cation-phosphate complexes (Mn+-Ins(2)P1), the alkali cations show little effect on the conformational preference while the conformational preference for the binding of the alkaline earth cations is pH-dependent and solvent-dependent. For example, these calculations predict that Mg2+-Ins(2)P1(0) and Mg2+-Ins(2)P1(2-) favor the 1a/5e form while Mg2+-Ins(2)P1(1-) favors the 5a/1e conformation. The Ca2+-Ins(2)P1(2-) complex prefers the 1a/5e conformation in the gas phase and in a nonpolar protein environment, but inverts to the 5a/1e conformation upon entering the polar aqueous phase. The binding affinities of the cations and the pK(a) values for the cation-phosphate complexes are derived from thermodynamical analysis.

The use of a generalized born model for the analysis of protein conformational transitions: a comparative study with explicit solvent simulations for chemotaxis Y protein (CheY)

To investigate whether implicit solvent models are appropriate for mechanistic studies of conformational transition in proteins, a recently developed generalized Born model (GBSW) was applied to a small signaling protein, chemotaxis protein Y (CheY), with different combinations of the phosphorylation state and conformation of the system; the results were compared to explicit solvent simulations using a stochastic boundary condition. The subtle but distinct conformational transitions involved in CheY activation makes the system ideally suited for comparing implicit and explicit solvent models because these conformational transitions are potentially accessible in both types of simulations. The structural and dynamical properties analyzed include not only those localized to the active site region but also throughout the protein, such as sidechain methyl group order parameters, backbone hydrogen bonding lifetime and occupancy as well as principal components of the trajectories. Overall, many properties were well reproduced by the GBSW simulations when compared with the explicit solvent calculations, although a number of observations consistently point to the suggestion that the current parameterization of the GBSW model tends to overestimate hydrogen-bonding interactions involving both charged groups and (charge-neutral) backbone atoms. This deficiency led to overstabilization of certain secondary structural motifs and more importantly, qualitatively different behaviors for the active site groups (Thr 87, Ala 88, the beta4-alpha4 loop) in response to phosphorylation, when compared with explicit solvent simulations. The current study highlights the value of carrying out both explicit and implicit solvent simulations for complementary mechanistic insights in the analysis of conformational transition in biomolecules.

Phosphorylation-Induced Conformational Changes in Short Peptides Probed by Fluorescence Resonance Energy Transfer in the 10 Å Domain

Phosphorylation-induced conformational changes in short polypeptides were probed by a fluorescence resonance energy transfer (FRET) method by employing a short-distance FRET pair (R(0) approximately 10 A) based on tryptophan as natural donor and a 2,3-diazabicyclo[2.2.2]oct-2-ene-labeled asparagine (Dbo) as synthetic acceptor. Two substrates for kinases, LeuArgArgTrpSerLeuGly-Dbo (peptide I) and TrpLysArgThrLeuArgArg-Dbo (peptide II), were investigated, with serine and threonine, respectively, as phosphorylation sites. Steady-state and time-resolved fluorescence experiments in H(2)O revealed a decrease in FRET efficiency for peptide I and an increase for peptide II; this suggested that the effective distances between donor and acceptor increased and decreased, respectively. The same trends and similar absolute variations in effective donor-acceptor distances were observed in propylene glycol, a less polar and highly viscous solvent; this suggested that the variations are due to intrinsic structural preferences. Fitting of the time-resolved decay traces according to a distribution function model (Gaussian distribution) provided the mean donor-acceptor distances, which showed an increase upon phosphorylation for peptide I (from 9.7 to 10.5 A) and a decrease for peptide II (from 10.9 to 9.3 A) in H(2)O. The broadness (half-width) of the distributions, which provides a measure of the rigidity of the peptides, remained similar upon phosphorylation of peptide I (3.0 versus 3.1 A), but decreased for peptide II (from 3.1 to 0.73 A in H(2)O); this suggests a more compact, structured conformation upon phosphorylation of the latter peptide. The elongation of the peptide backbone (by ca. 0.7 A) for peptide I is attributed to an increase in steric demand upon phosphorylation, which favors an extended conformation. The contraction (by ca. 1.4 A) and structural rigidification of peptide II is attributed to attractive Coulombic interactions and hydrogen bonding between the phosphate group and the arginine residues.

Insights into unbinding mechanisms upon two mutations investigated by molecular dynamics study of GSK3β-axin complex: Role of packing hydrophobic residues

Glycogen synthase kinase 3beta (GSK 3beta) is a key component of several cellular processes including Wnt and insulin signalling pathways. The interaction of GSK3beta with scaffolding peptide axin is thought to be responsible for the effective phosphorylation of beta-catenin, the core effector of Wnt signaling, which has been linked with the occurrence of colon cancer and melanoma. It has been demonstrated that the binding of axin to GSK3beta is abolished by the single-point mutation of Val267 to Gly (V267G) in GSK3beta or Leu392 to Pro (L392P) in axin. Molecular dynamics (MD) simulations were performed on wild type (WT), V267G mutant and L392P one to elucidate the two unbinding mechanisms that occur through different pathways. Besides, rough energy and residue-based energy decomposition were calculated by MM_GBSA (molecular mechanical Generalized_Born surface area) approach to illuminate the instability of the two mutants. The MD simulations of the two mutants and WT reveal that the structure of GSK3beta remains unchanged, while axin moves away from the interfacial hydrophobic pockets in both two mutants. Axin exhibits positional shift in V267G mutant, whereas, losing the hydrogen bonds that are indispensable for stabilizing the helix structure of wild type axin, the helix of axin is distorted in L392P mutant. To conclude, both two mutants destroy the hydrophobic interaction that is essential to the stability of GSK3beta-axin complex.

Effects of i and i+3 residue identity on cis-trans isomerism of the aromatic(i+1)-prolyl(i+2) amide bond: implications for type VI beta-turn formation

Cis-trans isomerization of amide bonds plays critical roles in protein molecular recognition, protein folding, protein misfolding, and disease. Aromatic-proline sequences are particularly prone to exhibit cis amide bonds. The roles of residues adjacent to a tyrosine-proline residue pair on cis-trans isomerism were examined. A short series of peptides XYPZ was synthesized and cis-trans isomerism was analyzed. Based on these initial studies, a series of peptides XYPN, X = all 20 canonical amino acids, was synthesized and analyzed by NMR for i residue effects on cis-trans isomerization. The following effects were observed: (a) aromatic residues immediately preceding Tyr-Pro disfavor cis amide bonds, with K(trans/cis)= 5.7-8.0, W > Y > F; (b) proline residues preceding Tyr-Pro lead to multiple species, exhibiting cis-trans isomerization of either or both X-Pro amide bonds; and (c) other residues exhibit similar values of K(trans/cis) (= 2.9-4.2), with Thr and protonated His exhibiting the highest fraction cis. beta-Branched and short polar residues were somewhat more favorable in stabilizing the cis conformation. Phosphorylation of serine at the i position modestly increases the stability of the cis conformer. In addition, the effect of the i+3 residue was examined in a limited series of peptides TYPZ. NMR data indicated that aromatic residues, Pro, Asn, Ala, and Val at the i+3 residue all favor cis amide bonds, with aromatic residues and Asn favoring more compact phi at Tyr(cis) and Ala and Pro favoring more extended phi at Tyr(cis). D-Alanine at the i+3 position particularly disfavors cis amide bonds.

Conformational changes in protein loops and helices induced by post-translational phosphorylation

Post-translational phosphorylation is a ubiquitous mechanism for modulating protein activity and protein-protein interactions. In this work, we examine how phosphorylation can modulate the conformation of a protein by changing the energy landscape. We present a molecular mechanics method in which we phosphorylate proteins in silico and then predict how the conformation of the protein will change in response to phosphorylation. We apply this method to a test set comprised of proteins with both phosphorylated and non-phosphorylated crystal structures, and demonstrate that it is possible to predict localized phosphorylation-induced conformational changes, or the absence of conformational changes, with near-atomic accuracy in most cases. Examples of proteins used for testing our methods include kinases and prokaryotic response regulators. Through a detailed case study of cyclin-dependent kinase 2, we also illustrate how the computational methods can be used to provide new understanding of how phosphorylation drives conformational change, why substituting Glu or Asp for a phosphorylated amino acid does not always mimic the effects of phosphorylation, and how a phosphatase can “capture” a phosphorylated amino acid. This work illustrates how computational methods can be used to elucidate principles and mechanisms of post-translational phosphorylation, which can ultimately help to bridge the gap between the number of known sites of phosphorylation and the number of structures of phosphorylated proteins.

Many proteins are chemically modified after they are synthesized in the cell. These post-translational modifications can modulate the ability of a protein to perform chemical reactions and to interact with other proteins. At the cellular level, for example, these chemical modifications are critical for allowing the cell to respond to its environment and control its division. One of the most common mechanisms by which proteins can be modified is by phosphorylation—the addition of a phosphate group to an amino acid side chain of the protein. Thousands of proteins are known to be modified by phosphorylation, but only for a small minority of these do we have any detailed understanding of how the chemical modification regulates the function of the protein. The authors describe a computational method that can make testable predictions about the structural changes that occur in a protein induced by post-translational phosphorylation. Their results show that the method can produce structural models of the phosphorylated proteins with near-atomic accuracy, and provide insight into the energetics of conformational switches driven by phosphorylation. As such, the computational method complements experiments aimed at understanding the mechanisms of protein regulation by phosphorylation.