Phosphorylation of an intrinsically disordered region of Ets1 shifts a multi-modal interaction ensemble to an auto-inhibitory state

The impact of IDR phosphorylation on the RNA binding profiles of proteins

Due to their capacity to mediate repetitive protein interactions, intrinsically disordered regions (IDRs) are crucial for the formation of various types of protein-RNA complexes. The functions of IDRs are strongly modulated by post-translational modifications (PTMs). Phosphorylation is the most common and well-studied modification of IDRs, which can alter homomeric or heteromeric interactions of proteins and impact their ability to phase separate. Moreover, phosphorylation can influence the RNA-binding properties of proteins, and recent studies demonstrated its selective impact on the global profiles of protein-RNA binding and regulation. These findings highlight the need for further integrative approaches to understand how signalling remodels protein-RNA networks in cells.

Effects of linker length on phase separation: lessons from the Rubisco-EPYC1 system of the algal pyrenoid

Biomolecular condensates are membraneless organelles formed via phase separation of macromolecules, typically consisting of bond-forming "stickers" connected by flexible "linkers". Linkers have diverse roles, such as occupying space and facilitating interactions. To understand how linker length relative to other lengths affects condensation, we focus on the pyrenoid, which enhances photosynthesis in green algae. Specifically, we apply coarse-grained simulations and analytical theory to the pyrenoid proteins of Chlamydomonas reinhardtii: the rigid holoenzyme Rubisco and its flexible partner EPYC1. Remarkably, halving EPYC1 linker lengths decreases critical concentrations by ten-fold. We attribute this difference to the molecular "fit" between EPYC1 and Rubisco. Varying Rubisco sticker locations reveals that the native sites yield the poorest fit, thus optimizing phase separation. Surprisingly, shorter linkers mediate a transition to a gas of rods as Rubisco stickers approach the poles. These findings illustrate how intrinsically disordered proteins affect phase separation through the interplay of molecular length scales.

Molecular Mechanisms of Functional Modulation of Transcriptional Coactivator PC4 via Phosphorylation on Its Intrinsically Disordered Region

To investigate the effects of phosphorylation on the function of the human positive cofactor 4 (PC4), an enhanced molecular dynamics (MD) simulation was performed. The simulation system consists of the N-terminal intrinsic disordered region (IDR) of PC4 and a complex that comprises the C-terminal acidic activation domain of a herpes simplex virion protein 16 (VP16ad) and a homodimer of the C-terminal structured core domain of PC4 (PC4ctd). An earlier report of an experimental study reported that the PC4-VP16ad interaction is modulated by incremental phosphorylation of the IDR. That report also proposed a dynamic model where phosphorylated serine residues of a segment (SEAC) in the IDR contact positively charged residues (lysin and arginine) of another segment (K-rich) in the IDR. This contact formation induced by the phosphorylation results in variation of PC4-VP16ad interaction. However, this contact formation has not yet been measured directly because it is transiently formed and because the SEAC and K-rich segments are unstructured with high flexibility. We performed two simulations to mimic the incremental phosphorylation: the IDR was not phosphorylated in one simulation and only partially phosphorylated in the other. Our simulation results indicate that the phosphorylation weakens the IDR-VP16ad contact considerably with the induction of a compact structure in the IDR. This structure was stabilized by electrostatic interactions between the phosphorylated serine residues of a segment and lysine or arginine residues of another segment in the IDR, but the conformational fluctuation of this compact structure was considerably large. Consequently, the present study supports the experimentally proposed dynamic model. Results of this study can be important for computational elucidation of the functional modulation of PC4.

ETS transcription factors: Multifaceted players from cancer progression to tumor immunity

The E26 transformation specific (ETS) family comprises 28 transcription factors, the majority of which are involved in tumor initiation and development. Serving as a group of functionally heterogeneous gene regulators, ETS factors possess a structurally conserved DNA-binding domain. As one of the most prominent families of transcription factors that control diverse cellular functions, ETS activation is modulated by multiple intracellular signaling pathways and post-translational modifications. Disturbances in ETS activity often lead to abnormal changes in oncogenicity, including cancer cell survival, growth, proliferation, metastasis, genetic instability, cell metabolism, and tumor immunity. This review systematically addresses the basics and advances in studying ETS factors, from their tumor relevance to clinical translational utility, with a particular focus on elucidating the role of ETS family in tumor immunity, aiming to decipher the vital role and clinical potential of regulation of ETS factors in the cancer field.

Non-Ewald methods for evaluating the electrostatic interactions of charge systems: similarity and difference

In molecular simulations, it is essential to properly calculate the electrostatic interactions of particles in the physical system of interest. Here we consider a method called the non-Ewald method, which does not rely on the standard Ewald method with periodic boundary conditions, but instead relies on the cutoff-based techniques. We focus on the physicochemical and mathematical conceptual aspects of the method in order to gain a deeper understanding of the simulation methodology. In particular, we take into account the reaction field (RF) method, the isotropic periodic sum (IPS) method, and the zero-multipole summation method (ZMM). These cutoff-based methods are based on different physical ideas and are completely distinguishable in their underlying concepts. The RF and IPS methods are "additive" methods that incorporate information outside the cutoff region, via dielectric medium and isotropic boundary condition, respectively. In contrast, the ZMM is a "subtraction" method that tries to remove the artificial effects, generated near the boundary, from the cutoff sphere. Nonetheless, we find physical and/or mathematical similarities between these methods. In particular, the modified RF method can be derived by the principle of neutralization utilized in the ZMM, and we also found a direct relationship between IPS and ZMM.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Zinc controls PML nuclear body formation through regulation of a paralog specific auto-inhibition in SUMO1

SUMO proteins are important regulators of many key cellular functions in part through their ability to form interactions with other proteins containing SUMO interacting motifs (SIMs). One characteristic feature of all SUMO proteins is the presence of a highly divergent intrinsically disordered region at their N-terminus. In this study, we examine the role of this N-terminal region of SUMO proteins in SUMO–SIM interactions required for the formation of nuclear bodies by the promyelocytic leukemia (PML) protein (PML-NBs). We demonstrate that the N-terminal region of SUMO1 functions in a paralog specific manner as an auto-inhibition domain by blocking its binding to the phosphorylated SIMs of PML and Daxx. Interestingly, we find that this auto-inhibition in SUMO1 is relieved by zinc, and structurally show that zinc stabilizes the complex between SUMO1 and a phospho-mimetic form of the SIM of PML. In addition, we demonstrate that increasing cellular zinc levels enhances PML-NB formation in senescent cells. Taken together, these results provide important insights into a paralog specific function of SUMO1, and suggest that zinc levels could play a crucial role in regulating SUMO1-SIM interactions required for PML-NB formation and function.

The auto-inhibition mechanism of transcription factor Ets-1 induced by phosphorylation on the intrinsically disordered region

As the most abundant post-translation modifications (PTMs), the phosphorylation usually occurred on the intrinsically disordered regions (IDRs). The regulation on the structures and interactions of IDRs induced by phosphorylation is critical to the function performing. The eukaryotic transcription factor 1 (Ets-1) is a member of transcription factor family, which participates in many important biological processes. The DNA-binding ability of Ets-1 is auto-inhibited by a disordered serine-rich region (SRR) on the Ets-1. The inhibition ability of SRR is greatly enhanced by the phosphorylation of the serine on the SRR. Nevertheless, the molecular mechanisms of the phosphorylation regulation on the structure and activity of Ets-1 are still unclear and under debates. By using both of the molecular simulations and biochemical experiments, we studied the molecule mechanism of phosphorylation regulation on the auto-inhibition of the Ets-1. The reasons of stabilization of Ets-1 core by phosphorylation on SRR region were elucidated. More important, the free energy landscapes (FEL) show that both of the steric hindrance and allosteric regulation are responsible for the DNA-binding inhibitory induced by phosphorylation, but the steric effects contribute greater than the allosteric regulation. The phosphorylation not only enhances the electrostatic interactions to facilitate the steric impedance, but also promotes the formation of hydrophobic residue clusters, which provide major driven force for the allosteric regulation. The structural basis of auto-inhibition of Ets-1 induced by the phosphorylation revealed in this study would great help the developing of inhibitor for the cancer therapy.

Extended ensemble simulations of a SARS-CoV-2 nsp1-5'-UTR complex

Nonstructural protein 1 (nsp1) of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is a 180-residue protein that blocks translation of host mRNAs in SARS-CoV-2-infected cells. Although it is known that SARS-CoV-2’s own RNA evades nsp1’s host translation shutoff, the molecular mechanism underlying the evasion was poorly understood. We performed an extended ensemble molecular dynamics simulation to investigate the mechanism of the viral RNA evasion. Simulation results suggested that the stem loop structure of the SARS-CoV-2 RNA 5’-untranslated region (SL1) binds to both nsp1’s N-terminal globular region and intrinsically disordered region. The consistency of the results was assessed by modeling nsp1-40S ribosome structure based on reported nsp1 experiments, including the X-ray crystallographic structure analysis, the cryo-EM electron density map, and cross-linking experiments. The SL1 binding region predicted from the simulation was open to the solvent, yet the ribosome could interact with SL1. Cluster analysis of the binding mode and detailed analysis of the binding poses suggest residues Arg124, Lys47, Arg43, and Asn126 may be involved in the SL1 recognition mechanism, consistent with the existing mutational analysis.

The pandemic of COVID-19 is still rampant all over the world as of 2021 June. SARS-CoV-2 (severe acute respiratory syndrome coronavirus 2), the causative pathogen of COVID-19, encodes a protein called nsp1 (nonstructural protein 1), which modulates and hijacks the ribosome of the infected host cells. With nsp1, infected human cells selectively translate SARS-CoV-2’s RNA, which increases the virus reproduction efficiency while evading the host immunity. Though it has been known that nsp1 recognizes characteristic stem-loop structure at 5’-end of SARS-CoV-2’s RNA (called SL1), the molecular mechanism underlying the recognition has been poorly understood. We investigated the mechanism of selective translation using the all-atom molecular dynamics simulation of nsp1-SL1 complex. Our simulation results suggest that the binding between nsp1 and SL1 is multi-modal. The results also imply that both the N-terminal globular part and the C-terminal flexible tail of nsp1 are involved in the binding. The residues involved in nsp1-SL1 binding coincides with the known mutant analyses of SARS-CoV-1 and SARS-CoV-2, as well as experimental evidence about nsp1-ribosome interactions.

Microsecond molecular dynamics suggest that a non-synonymous mutation, frequently observed in patients with mild symptoms in Tokyo, alters dynamics of the SARS-CoV-2 main protease

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which causes the coronavirus disease 2019 (COVID-19), spread rapidly around the globe. The main protease encoded by SARS-CoV-2 is essential for processing of the polyproteins translated from the viral RNA genome, making this protein a potential drug target. A recently reported mutation in the protease, P108S, may be responsible for milder symptoms observed in COVID-19 patients in Tokyo. Starting from a crystal structure of the SARS-CoV-2 main protease in the dimeric form, we performed triplicate 5.0-μs molecular dynamics simulations of the wild-type and P108S mutant. Our computational results suggest a link between the mutation P108S and dynamics of the catalytic sites in the main protease: The catalytic dyad become considerably inaccessible to substrates in the P108S mutant. Our results also demonstrate the potential of molecular dynamics simulations to complement experimental techniques and other computational methods, such as protein design calculations, which predict effects of mutations based on static crystal structures. Further studies are certainly necessary to quantitively understand the relationships between the P108S mutation and physical properties of the main protease, but the results of our study will immediately inform development of new protease inhibitors.

The diversity of molecular interactions involving intrinsically disordered proteins: A molecular modeling perspective

Intrinsically Disordered Proteins and Regions (IDPs/IDRs) are key components of a multitude of biological processes. Conformational malleability enables IDPs/IDRs to perform very specialized functions that cannot be accomplished by globular proteins. The functional role for most of these proteins is related to the recognition of other biomolecules to regulate biological processes or as a part of signaling pathways. Depending on the extent of disorder, the number of interacting sites and the type of partner, very different architectures for the resulting assemblies are possible. More recently, molecular condensates with liquid-like properties composed of multiple copies of IDPs and nucleic acids have been proven to regulate key processes in eukaryotic cells. The structural and kinetic details of disordered biomolecular complexes are difficult to unveil experimentally due to their inherent conformational heterogeneity. Computational approaches, alone or in combination with experimental data, have emerged as unavoidable tools to understand the functional mechanisms of this elusive type of assemblies. The level of description used, all-atom or coarse-grained, strongly depends on the size of the molecular systems and on the timescale of the investigated mechanism. In this mini-review, we describe the most relevant architectures found for molecular interactions involving IDPs/IDRs and the computational strategies applied for their investigation.

Fundamental Challenges and Outlook in Simulating Liquid-Liquid Phase Separation of Intrinsically Disordered Proteins

Intrinsically disordered proteins (IDPs) populate an ensemble of dynamic conformations, making their structural characterization by experiments challenging. Many IDPs undergo liquid-liquid phase separation into dense membraneless organelles with myriad cellular functions. Multivalent interactions in low-complexity IDPs promote the formation of these subcellular coacervates. While solution NMR, Förster resonance energy transfer (FRET), and small-angle X-ray scattering (SAXS) studies on IDPs have their own challenges, recent computational methods draw a rational trade-off to characterize the driving forces underlying phase separation. In this Perspective, we critically evaluate the scope of approximation-free field theoretic simulations, well-tempered ensemble methods, enhanced sampling techniques, coarse-grained force fields, and slab simulation approaches to offer an improved understanding of phase separation. A synergy between simulation length scale and model resolution would reduce the existing caveats and enable theories of polymer physics to elucidate finer details of liquid-liquid phase separation (LLPS). These computational advances offer promise for rigorous characterization of the IDP proteome and designing peptides with tunable material and self-assembly properties.

Molecular Dynamics Simulations in Drug Discovery and Pharmaceutical Development

Molecular dynamics (MD) simulations have become increasingly useful in the modern drug development process. In this review, we give a broad overview of the current application possibilities of MD in drug discovery and pharmaceutical development. Starting from the target validation step of the drug development process, we give several examples of how MD studies can give important insights into the dynamics and function of identified drug targets such as sirtuins, RAS proteins, or intrinsically disordered proteins. The role of MD in antibody design is also reviewed. In the lead discovery and lead optimization phases, MD facilitates the evaluation of the binding energetics and kinetics of the ligand-receptor interactions, therefore guiding the choice of the best candidate molecules for further development. The importance of considering the biological lipid bilayer environment in the MD simulations of membrane proteins is also discussed, using G-protein coupled receptors and ion channels as well as the drug-metabolizing cytochrome P450 enzymes as relevant examples. Lastly, we discuss the emerging role of MD simulations in facilitating the pharmaceutical formulation development of drugs and candidate drugs. Specifically, we look at how MD can be used in studying the crystalline and amorphous solids, the stability of amorphous drug or drug-polymer formulations, and drug solubility. Moreover, since nanoparticle drug formulations are of great interest in the field of drug delivery research, different applications of nano-particle simulations are also briefly summarized using multiple recent studies as examples. In the future, the role of MD simulations in facilitating the drug development process is likely to grow substantially with the increasing computer power and advancements in the development of force fields and enhanced MD methodologies.

Short Disordered Epitope of CRTAM Ig-Like V Domain as a Potential Target for Blocking Antibodies

Class-I Restricted T Cell-Associated Molecule (CRTAM) is a protein that is expressed after T cell activation. The interaction of CRTAM with its ligand, nectin-like 2 (Necl2), is required for the efficient production of IL-17, IL-22, and IFNγ by murine CD4 T cells, and it plays a role in optimal CD8 T and NK cell cytotoxicity. CRTAM promotes the pro-inflammatory cytokine profile; therefore, it may take part in the immunopathology of autoimmune diseases such as diabetes type 1 or colitis. Thus, antibodies that block the interaction between CRTAM and Necl2 would be useful for controlling the production of these inflammatory cytokines. In this work, using bioinformatics predictions, we identified three short disordered epitopes (sDE1-3) that are located in the Ig-like domains of murine CRTAM and are conserved in mammalian species. We performed a structural analysis by molecular dynamics simulations of sDE1 (QHPALKSSKY, Ig-like V), sDE2 (QRNGEKSVVK, Ig-like C1), and sDE3 (CSTERSKKPPPQI, Ig-like C1). sDE1, which is located within a loop of the contact interface of the heterotypic interaction with Nectl2, undergoes an order–disorder transition. On the contrary, even though sDE2 and sDE3 are flexible and also located within loops, they do not undergo order–disorder transitions. We evaluated the immunogenicity of sDE1 and sDE3 through the expression of these epitopes in chimeric L1 virus-like particles. We confirmed that sDE1 induces polyclonal antibodies that recognize the native folding of CRTAM expressed in activated murine CD4 T cells. In contrast, sDE3 induces polyclonal antibodies that recognize the recombinant protein hCRTAM-Fc, but not the native CRTAM. Thus, in this study, an exposed disordered epitope in the Ig-like V domain of CRTAM was identified as a potential site for therapeutic antibodies.

The human CRY1 tail controls circadian timing by regulating its association with CLOCK:BMAL1

Circadian rhythms are generated by interlocked transcription-translation feedback loops that establish cell-autonomous biological timing of ∼24 h. Mutations in core clock genes that alter their stability or affinity for one another lead to changes in circadian period. The human CRY1Δ11 mutant lengthens circadian period to cause delayed sleep phase disorder (DSPD), characterized by a very late onset of sleep. CRY1 is a repressor that binds to the transcription factor CLOCK:BMAL1 to inhibit its activity and close the core feedback loop. We previously showed how the PHR (photolyase homology region) domain of CRY1 interacts with distinct sites on CLOCK and BMAL1 to sequester the transactivation domain from coactivators. However, the Δ11 variant alters an intrinsically disordered tail in CRY1 downstream of the PHR. We show here that the CRY1 tail, and in particular the region encoded by exon 11, modulates the affinity of the PHR domain for CLOCK:BMAL1. The PHR-binding epitope in exon 11 is necessary and sufficient to disrupt the interaction between CRY1 and the subunit CLOCK. Moreover, PHR-tail interactions are conserved in the paralog CRY2 and reduced when either CRY is bound to the circadian corepressor PERIOD2. Discovery of this autoregulatory role for the mammalian CRY1 tail and conservation of PHR-tail interactions in both mammalian cryptochromes highlights functional conservation with plant and insect cryptochromes, which also utilize PHR-tail interactions to reversibly control their activity.

Molecular Interaction Mechanism of a 14-3-3 Protein with a Phosphorylated Peptide Elucidated by Enhanced Conformational Sampling

Enhanced conformational sampling, a genetic-algorithm-guided multidimensional virtual-system coupled molecular dynamics, can provide equilibrated conformational distributions of a receptor protein and a flexible ligand at room temperature. The distributions provide not only the most stable but also semistable complex structures and propose a ligand-receptor binding process. This method was applied to a system consisting of a receptor protein, 14-3-3ε, and a flexible peptide, phosphorylated myeloid leukemia factor 1 (pMLF1). The results present comprehensive binding pathways of pMLF1 to 14-3-3ε. We identified four thermodynamically stable clusters of MLF1 on the 14-3-3ε surface and free-energy barriers among some clusters. The most stable cluster includes two high-density spots connected by a narrow corridor. When pMLF1 passes the corridor, a salt-bridge relay (switching) related to the phosphorylated residue of pMLF1 occurs. Conformations in one high-density spot are similar to the experimentally determined complex structure. Three-dimensional distributions of residues in the intermolecular interface rationally explain the binding constant changes resulting from the alanine mutation experiment for the residues. We also performed a simulation of nonphosphorylated peptide and 14-3-3ε, which demonstrated that the complex structure was unstable, suggesting that phosphorylation of the peptide is crucially important for binding to 14-3-3ε.

myPresto/omegagene 2020: a molecular dynamics simulation engine for virtual-system coupled sampling

The molecular dynamics (MD) method is a promising approach for investigating the molecular mechanisms of microscopic phenomena. In particular, generalized ensemble MD methods can efficiently explore the conformational space with a rugged free-energy surface. However, the implementation and acquisition of technical knowledge for each generalized ensemble MD method are not straightforward for end-users. Here, we present a new version of the myPresto/omegagene software, which is an MD simulation engine tailored for a series of generalized ensemble methods, which are virtual-system coupled multicanonical MD (V-McMD), virtual-system coupled adaptive umbrella sampling (V-AUS), and virtual-system coupled canonical MD (VcMD). This program has been applied in several studies analyzing free-energy landscapes of a variety of molecular systems with all-atom simulations. The updated version provides new functionality for coarse-grained simulations powered by the hydrophobicity scale method. The software package includes a step-by-step tutorial document for enhanced conformational sampling of the poly-glutamine (poly-Q) oligomer expressed as a one-bead per residue model. The myPresto/omegagene software is freely available at the following URL: https://github.com/kotakasahara/omegagene under the Apache2 license.

Bioinformatic Analysis and Biophysical Characterization Reveal Structural Disorder in G0S2 Protein

G0S2 is a small protein of 103 residues in length that is involved in multiple cellular processes. To date, several reports have shown that G0S2 functions by making direct protein-protein interactions with key proteins. In lipolysis, G0S2 specifically interacts with adipose triglyceride lipase, inhibiting its activity and resulting in lipolysis being downregulated. In a similar way, G0S2 also participates in the regulation of apoptosis, cell proliferation, and oxidative phosphorylation; however, information regarding G0S2 structural and biophysical properties is limited. In this work, we conducted a comparative structural analysis of human and mouse G0S2 proteins. Bioinformatics suggests the presence of a disordered C-terminal region in human G0S2. Experimental characterization by size-exclusion chromatography and dynamic light scattering showed that human and mouse G0S2 have different hydrodynamic properties. In comparison to the mouse G0S2, which behaves similar to a globular protein, the human G0S2 shows an elongated conformation, most likely by displaying a disordered C-terminal region. Further analysis of hydrodynamic properties under denaturing conditions suggests the presence of a structural element in the mouse protein that undergoes an order to disorder transition at low urea concentration. Structural analysis by circular dichroism revealed that in native conditions, both proteins are mainly unstructured, showing the presence of beta sheet structures. Further analysis of CD data suggests that both proteins belong to the premolten globule family of intrinsically disordered proteins. We suggest that the intrinsic disorder observed in the G0S2 protein may facilitate its interaction with multiple partners in the regulation of cellular metabolism.

Multifunnel Energy Landscapes for Phosphorylated Translation Repressor 4E-BP2 and Its Mutants

Upon phosphorylation of specific sites, eukaryotic translation initiation factor 4E (eIF4E) binding protein 2 (4E-BP2) undergoes a fundamental structural transformation from a disordered state to a four-stranded β-sheet, leading to decreased binding affinity for its partner. This change reflects the significant effects of phosphate groups on the underlying energy landscapes of proteins. In this study, we combine high-temperature molecular dynamics simulations and discrete path sampling to construct energy landscapes for a doubly phosphorylated 4E-BP2_18-62 and two mutants (a single site mutant D33K and a double mutant Y54A/L59A). The potential and free energy landscapes for these three systems are multifunneled with the folded state and several alternative states lying close in energy, suggesting perhaps a multifunneled and multifunctional protein. Hydrogen bonds between phosphate groups and other residues not only stabilize these low-lying conformations to different extents but also play an important role in interstate transitions. From the energy landscape perspective, our results explain some interesting experimental observations, including the low stability of doubly phosphorylated 4E-BP2 and its moderate binding to eIF4E and the inability of phosphorylated Y54A/L59A to fold.

Quantitative proteomics indicate a strong correlation of mitotic phospho-/dephosphorylation with non-structured regions of substrates

Protein phosphorylation plays a critical role in the regulation and progression of mitosis. >40,000 phosphorylated residues and the associated kinases have been identified to date via proteomic analyses. Although some of these phosphosites are associated with regulation of either protein-protein interactions or the catalytic activity of the substrate protein, the roles of most mitotic phosphosites remain unclear. In this study, we examined structural properties of mitotic phosphosites and neighboring residues to understand the role of heavy phosphorylation in non-structured domains. Quantitative mass spectrometry analysis of mitosis-arrested and non-arrested HeLa cells revealed >4100 and > 2200 residues either significantly phosphorylated or dephosphorylated, respectively, at mitotic entry. The calculated disorder scores of amino acid sequences of neighboring individual phosphosites revealed that >70% of dephosphorylated phosphosites exist in disordered regions, whereas 50% of phosphorylated sites exist in non-structured domains. A clear inverse correlation was observed between probability of phosphorylation in non-structured domain and increment of phosphorylation in mitosis. These results indicate that at entry to mitosis, a significant number of phosphate groups are removed from non-structured domains and transferred to more-structured domains. Gene ontology term analysis revealed that mitosis-related proteins are heavily phosphorylated, whereas RNA-related proteins are both dephosphorylated and phosphorylated, suggesting that heavy phosphorylation/dephosphorylation in non-structured domains of RNA-binding proteins plays a role in dynamic rearrangement of RNA-containing organelles, as well as other intracellular environments.

Multidimensional virtual-system coupled canonical molecular dynamics to compute free-energy landscapes of peptide multimer assembly

An enhanced-sampling method termed multidimensional virtual-system coupled canonical molecular dynamics (mD-VcMD) method is developed. In many cases, generalized-ensemble methods realizing enhanced sampling, for example, adaptive umbrella sampling, apply an effective potential, which is derived from temporarily assumed canonical distribution as a function of one or more arbitrarily defined reaction coordinates. However, it is not straightforward to estimate the appropriate canonical distribution, especially for cases applying multiple reaction coordinates. The current method, mD-VcMD, does not rely on the form of the canonical distribution. Therefore, it is practically useful to explore a high-dimensional reaction-coordinate space. In this article, formulation of mD-VcMD and its evaluation with the simple molecular models consisting of three or four alanine peptides are presented. We confirmed that mD-VcMD efficiently searched 2D and 3D reaction-coordinate spaces defined as interpeptide distances. Direct comparisons with results of long-term canonical MD simulations revealed that mD-VcMD produces correct canonical ensembles. © 2019 Wiley Periodicals, Inc.

Robustness by intrinsically disordered C-termini and translational readthrough

During protein synthesis genetic instructions are passed from DNA via mRNA to the ribosome to assemble a protein chain. Occasionally, stop codons in the mRNA are bypassed and translation continues into the untranslated region (3′-UTR). This process, called translational readthrough (TR), yields a protein chain that becomes longer than would be predicted from the DNA sequence alone. Protein sequences vary in propensity for translational errors, which may yield evolutionary constraints by limiting evolutionary paths. Here we investigated TR in Saccharomyces cerevisiae by analysing ribosome profiling data. We clustered proteins as either prone or non-prone to TR, and conducted comparative analyses. We find that a relatively high frequency (5%) of genes undergo TR, including ribosomal subunit proteins. Our main finding is that proteins undergoing TR are highly expressed and have a higher proportion of intrinsically disordered C-termini. We suggest that highly expressed proteins may compensate for the deleterious effects of TR by having intrinsically disordered C-termini, which may provide conformational flexibility but without distorting native function. Moreover, we discuss whether minimizing deleterious effects of TR is also enabling exploration of the phenotypic landscape of protein isoforms.

Studies on Molecular Dynamics of Intrinsically Disordered Proteins and Their Fuzzy Complexes: A Mini-Review

The molecular dynamics (MD) method is a promising approach toward elucidating the molecular mechanisms of intrinsically disordered regions (IDRs) of proteins and their fuzzy complexes. This mini-review introduces recent studies that apply MD simulations to investigate the molecular recognition of IDRs. Firstly, methodological issues by which MD simulations treat IDRs, such as developing force fields, treating periodic boundary conditions, and enhanced sampling approaches, are discussed. Then, several examples of the applications of MD to investigate molecular interactions of IDRs in terms of the two kinds of complex formations; coupled-folding and binding and fuzzy complex. MD simulations provide insight into the molecular mechanisms of these binding processes by sampling conformational ensembles of flexible IDRs. In particular, we focused on all-atom explicit-solvent MD simulations except for studies of higher-order assembly of IDRs. Recent advances in MD methods, and computational power make it possible to dissect the molecular details of realistic molecular systems involving the dynamic behavior of IDRs.

Hydrogen bond donors and acceptors are generally depolarized in α-helices as revealed by a molecular tailoring approach

Hydrogen‐bond (H‐bond) interaction energies in α‐helices of short alanine peptides were systematically examined by precise density functional theory calculations, followed by a molecular tailoring approach. The contribution of each H‐bond interaction in α‐helices was estimated in detail from the entire conformation energies, and the results were compared with those in the minimal H‐bond models, in which only H‐bond donors and acceptors exist with the capping methyl groups. The former interaction energies were always significantly weaker than the latter energies, when the same geometries of the H‐bond donors and acceptors were applied. The chemical origin of this phenomenon was investigated by analyzing the differences among the electronic structures of the local peptide backbones of the α‐helices and those of the minimal H‐bond models. Consequently, we found that the reduced H‐bond energy originated from the depolarizations of both the H‐bond donor and acceptor groups, due to the repulsive interactions with the neighboring polar peptide groups in the α‐helix backbone. The classical force fields provide similar H‐bond energies to those in the minimal H‐bond models, which ignore the current depolarization effect, and thus they overestimate the actual H‐bond energies in α‐helices. © 2019 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.