Revised AMBER parameters for bioorganic phosphates

On the Power and Challenges of Atomistic Molecular Dynamics to Investigate RNA Molecules

RNA molecules play a vital role in biological processes within the cell, with significant implications for science and medicine. Notably, the biological functions exerted by specific RNA molecules are often linked to the RNA conformational ensemble. However, the experimental characterization of such three-dimensional RNA structures is challenged by the structural heterogeneity of RNA and by its multiple dynamic interactions with binding partners such as small molecules, proteins, and metal ions. Consequently, our current understanding of the structure-function relationship of RNA molecules is still limited. In this context, we highlight molecular dynamics (MD) simulations as a powerful tool to complement experimental efforts on RNAs. Despite the recognized limitations of current force fields for RNA MD simulations, examining the dynamics of selected RNAs has provided valuable functional insights into their structures.

nP-Collabs: Investigating Counterion-Mediated Bridges in the Multiply Phosphorylated Tau-R2 Repeat

Tau is an intrinsically disordered (IDP) microtubule-associated protein (MAP) that plays a key part in microtubule assembly and organization. The function of tau can be regulated by multiple phosphorylation sites. These post-translational modifications are known to decrease the binding affinity of tau for microtubules, and abnormal tau phosphorylation patterns are involved in Alzheimer's disease. Using all-atom molecular dynamics simulations, we compared the conformational landscapes explored by the tau R2 repeat domain (which comprises a strong tubulin binding site) in its native state and with multiple phosphorylations on the S285, S289, and S293 residues, with four different standard force field (FF)/water model combinations. We find that the different parameters used for the phosphate groups (which can be more or less flexible) in these FFs and the specific interactions between bulk cations and water lead to the formation of a specific type of counterion bridge, termed nP-collab (for nphosphate collaboration, with n being an integer), where counterions form stable structures binding with two or three phosphate groups simultaneously. The resulting effect of nP-collabs on the tau-R2 conformational space differs when using sodium or potassium cations and is likely to impact the peptide overall dynamics and how this MAP interacts with tubulins. We also investigated the effect of phosphoresidue spacing and ionic concentration by modeling polyalanine peptides containing two phosphoserines located one-six residues apart. Three new metrics specifically tailored for IDPs (proteic Menger curvature, local curvature, and local flexibility) were introduced, which allow us to fully characterize the impact of nP-collabs on the dynamics of disordered peptides at the residue level.

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.

Assessing RNA atomistic force fields via energy landscape explorations in implicit solvent

Predicting the structure and dynamics of RNA molecules still proves challenging because of the relative scarcity of experimental RNA structures on which to train models and the very sensitive nature of RNA towards its environment. In the last decade, several atomistic force fields specifically designed for RNA have been proposed and are commonly used for simulations. However, it is not necessarily clear which force field is the most suitable for a given RNA molecule. In this contribution, we propose the use of the computational energy landscape framework to explore the energy landscape of RNA systems as it can bring complementary information to the more standard approaches of enhanced sampling simulations based on molecular dynamics. We apply the EL framework to the study of a small RNA pseudoknot, the Aquifex aeolicus tmRNA pseudoknot PK1, and we compare the results of five different RNA force fields currently available in the AMBER simulation software, in implicit solvent. With this computational approach, we can not only compare the predicted 'native' states for the different force fields, but the method enables us to study metastable states as well. As a result, our comparison not only looks at structural features of low energy folded structures, but provides insight into folding pathways and higher energy excited states, opening to the possibility of assessing the validity of force fields also based on kinetics and experiments providing information on metastable and unfolded states.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-024-01202-9.

Binding of Phosphate Species to Ca2+ and Mg2+ in Aqueous Solution

Phosphate derivatives and their interaction with metal cations are involved in many important biological phenomena, so an accurate characterization of the phosphate-metal interaction is necessary to properly understand the role of phosphate-metal contacts in mediating biological function. Herein, we improved the standard 12-6 Lennard-Jones (LJ) potential via the usage of the 12-6-4 LJ model, which incorporates ion-induced dipole interactions. Via parameter scanning, we fine-tuned the 12-6-4 LJ polarizability values to obtain accurate absolute binding free energies for the phosphate anions H2PO4-, HPO42-, PO43- coordinating with Ca2+ and Mg2+. First, we modified the phosphate 12-6-4 LJ parameters to reproduce the solvation free energies of the series of phosphate anions using the thermodynamic integration (TI) method. Then, using the potential mean force (PMF) method, the polarizability of the metal-phosphate interaction was obtained. We show that the free energy profiles of phosphate ions coordinated to Ca2+ and Mg2+ generally show similar trends at longer metal-phosphate distances, while the absolute binding energy values increased with deprotonation. The resulting parameters demonstrate the flexibility of the 12-6-4 LJ-type nonbonded model and its usefulness in accurately describing cation-anion interactions.

The effect of pseudoknot base pairing on cotranscriptional structural switching of the fluoride riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of Escherichia coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37°C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65°C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

Atomistic simulations of RNA duplex thermal denaturation: Sequence- and forcefield-dependence

Double-stranded RNA is the end-product of template-based replication, and is also the functional state of some biological RNAs. Similarly to proteins and DNA, they can be denatured by temperature, with important physiological and technological implications. Here, we use an in silico strategy to probe the thermal denaturation of RNA duplexes. Following previous results that were obtained on a few different duplexes, and which nuanced the canonical 2-state picture of nucleic acid denaturation, we here specifically address three different aspects that greatly improve our description of the temperature-induced dsRNA separation. First, we investigate the effect of the spatial distribution of weak and strong base-pairs among the duplex sequence. We show that the deviations from the two-state dehybridization mechanism are more pronounced when a strong core is flanked with weak extremities, while duplexes with a weak core but strong extremities exhibit a two-state behavior, which can be explained by the key role played by base fraying. This was later verified by generating artificial hairpin or circular states containing one or two locked duplex extremities, which results in an important reinforcement of the entire HB structure of the duplex and higher melting temperatures. Finally, we demonstrate that our results are little sensitive to the employed combination of RNA and water forcefields. The trends in thermal stability among the different sequences as well as the observed unfolding mechanisms (and the deviations from a two-state scenario) remain the same regardless of the employed atomistic models. However, our study points to possible limitations of recent reparametrizations of the Amber RNA forcefield, which sometimes results in duplexes that readily denature under ambient conditions, in contradiction with available experimental results.

4'-C-Acetamidomethyl-2'-O-methoxyethyl Nucleic Acid Modifications Improve Thermal Stability, Nuclease Resistance, Potency, and hAgo2 Binding of Small Interfering RNAs

In this study, we designed the 4'-C-acetamidomethyl-2'-O-methoxyethyl (4'-C-ACM-2'-O-MOE) uridine and thymidine modifications, aiming to test them into small interfering RNAs. Thermal melting studies revealed that incorporating a single 4'-C-ACM-2'-O-MOE modification in the DNA duplex reduced thermal stability. In contrast, an increase in thermal stability was observed when the modification was introduced in DNA:RNA hybrid and in siRNAs. Thermal destabilization in DNA duplex was attributed to unfavorable entropy, which was mainly compensated by the enthalpy factor to some extent. A single 4'-C-ACM-2'-O-MOE thymidine modification at the penultimate position of the 3'-end of dT20 oligonucleotides in the presence of 3'-specific exonucleases, snake venom phosphodiesterase (SVPD), demonstrated significant stability as compared to monomer modifications including 2'-O-Me, 2'-O-MOE, and 2'-F. In gene silencing studies, we found that the 4'-C-ACM-2'-O-MOE uridine or thymidine modifications at the 3'-overhang in the passenger strand in combination with two 2'-F modifications exhibited superior RNAi activity. The results suggest that the dual modification is well tolerated at the 3'-end of the passenger strand, which reflects better siRNA stability and silencing activity. Interestingly, 4'-C-ACM-2'-O-MOE-modified siRNAs showed considerable gene silencing even after 96 h posttransfection; it showed that our modification could induce prolonged gene silencing due to improved metabolic stability. Molecular modeling studies revealed that the introduction of the 4'-C-ACM-2'-O-MOE modification at the 3'-end of the siRNA guide strand helps to anchor the strand within the PAZ domain of the hAgo2 protein. The overall results indicate that the 4'-C-ACM-2'-O-MOE uridine and thymidine modifications are promising modifications to improve the stability, potency, and hAgo2 binding of siRNAs.

Molecular Driving Forces in the Self-Association of Silaffin Peptide R5 from MD Simulations

The 19-residue silaffin-R5 peptide has been widely studied for its ability to precipitate uniform SiO2 particles through mild temperature and pH pathways, in the absence of any organic solvents. There is consensus that post-translational modification (PTM) of side chains has a large impact on the biomineralization process. Thus, it is imperative to understand the precise mechanisms that dictate the formation of SiO2 from R5 peptide, including the effects of PTM on peptide aggregation and peptide-surface adsorption. In this work, we use molecular dynamics (MD) simulations to study the aggregation of R5 dimer with multiple PTMs, with the presence of different ions in solution. Since this system has strong interactions with deep metastable states, we use parallel bias metadynamics with partitioned families to efficiently sample the different states of the system. We find that peptide aggregation is a prerequisite for biomineralization. We observe that the electrostatic interactions are essential in the R5 dimer aggregation; for wild type R5 that only has positively charged residues, phosphate ions HPO4 2- in the solution form a bridge between two peptides and are essential for peptide aggregation.

N3-Methyluridine and 2'-O-Alkyl/2'-Fluoro-N3-methyluridine functionalized nucleic acids improve nuclease resistance while maintaining duplex geometry

Herein, we report the synthesis of 2'-O-alkyl/2'-fluoro-N3-methyluridine (2'-O-alkyl/2'-F-m3U) phosphoramidites and their incorporation in DNA and RNA oligonucleotides. The duplex binding affinity and base discrimination studies showed that all 2'-O-alkyl/2'-F-m3U modifications significantly decreased the thermal stability and base-pairing discrimination ability. Serum stability study of dT20 with 2'-O-alkyl-m3U modification exhibited excellent nuclease resistance when incubated with 3'-exonucleases (SVPD) or 5'-exonucleases (PDE-II) as compared to m3U, 2'-F, 2'-OMe modified oligonucleotides. MD simulation studies with RNA tetradecamer duplexes illustrated that the m3U and 2'-O-methyl-m3U modifications reduce the duplex stabilities by disrupting the Watson-Crick hydrogen bonding and base-stacking interactions. Further molecular modelling investigations demonstrated that the 2'-O-propyl-m3U modification exhibits steric interactions with amino acid residues in the active site of 3'- and 5'-exonuclease, leading to enhanced stability. These combined data indicate that the 2'-modified-m3U nucleotides can be used as a promising tool to enhance the stability, silencing efficiency, and drug-like properties of antisense/siRNA-based therapeutics.

van der Waals Parameter Scanning with Amber Nucleic Acid Force Fields: Revisiting Means to Better Capture the RNA/DNA Structure through MD

Molecular dynamics simulations can be used in combination with experimental techniques to uncover the intricacies of biomolecular structure, dynamics, and the resulting interactions. However, many noncanonical nucleic acid structures have proven to be challenging to replicate in accurate agreement with experimental data, often attributed to known force field deficiencies. A common force field criticism is the handling of van der Waals (vdW) parameters, which have not been updated since the regular use of Ewald's methods became routine. This work dives into the effects of minute vdW radii shifts on RNA tetranucleotide, B-DNA, and Z-DNA model systems described by commonly used Amber force fields. Using multidimensional replica exchange molecular dynamics (M-REMD), the GACC RNA tetranucleotide demonstrated changes in the structural distribution between the NMR minor and anomalous structure populations based on the O2' vdW radii scanning. However, no significant change in the NMR Major conformation population was observed. There were minimal changes in the B-DNA structure but there were more substantial improvements in Z-DNA structural descriptions, specifically with the Tumuc1 force field. This occurred with both LJbb vdW radii adjustments and incorporation of the CUFIX nonbonded parameter modifications. Though the limited vdW modifications tested did not provide a universal fix to the challenge of simulating the various known nucleic acid structures, they do provide direction and a greater understanding for future force field development efforts.

The Effect of Pseudoknot Base Pairing on Cotranscriptional Structural Switching of the Fluoride Riboswitch

A central question in biology is how RNA sequence changes influence dynamic conformational changes during cotranscriptional folding. Here we investigated this question through the study of transcriptional fluoride riboswitches, non-coding RNAs that sense the fluoride anion through the coordinated folding and rearrangement of a pseudoknotted aptamer domain and a downstream intrinsic terminator expression platform. Using a combination of E. coli RNA polymerase in vitro transcription and cellular gene expression assays, we characterized the function of mesophilic and thermophilic fluoride riboswitch variants. We showed that only variants containing the mesophilic pseudoknot function at 37 °C. We next systematically varied the pseudoknot sequence and found that a single wobble base pair is critical for function. Characterizing thermophilic variants at 65 °C through Thermus aquaticus RNA polymerase in vitro transcription showed the importance of this wobble pair for function even at elevated temperatures. Finally, we performed all-atom molecular dynamics simulations which supported the experimental findings, visualized the RNA structure switching process, and provided insight into the important role of magnesium ions. Together these studies provide deeper insights into the role of riboswitch sequence in influencing folding and function that will be important for understanding of RNA-based gene regulation and for synthetic biology applications.

RNAs undergo phase transitions with lower critical solution temperatures

Co-phase separation of RNAs and RNA-binding proteins drives the biogenesis of ribonucleoprotein granules. RNAs can also undergo phase transitions in the absence of proteins. However, the physicochemical driving forces of protein-free, RNA-driven phase transitions remain unclear. Here we report that various types of RNA undergo phase separation with system-specific lower critical solution temperatures. This entropically driven phase separation is an intrinsic feature of the phosphate backbone that requires Mg2+ ions and is modulated by RNA bases. RNA-only condensates can additionally undergo enthalpically favourable percolation transitions within dense phases. This is enabled by a combination of Mg2+-dependent bridging interactions between phosphate groups and RNA-specific base stacking and base pairing. Phase separation coupled to percolation can cause dynamic arrest of RNAs within condensates and suppress the catalytic activity of an RNase P ribozyme. Our work highlights the need to incorporate RNA-driven phase transitions into models for ribonucleoprotein granule biogenesis.

Exploring the impact of mutation and post-translational modification on α-Synuclein: Insights from molecular dynamics simulations with and without copper

We report molecular dynamics simulations of two modifications to α-Synuclein, namely A53T mutation and phosphorylation at Ser129, which have been observed in Parkinson's disease patients. Both modifications are close to known metal binding sites, so as well as each modified peptide we also study Cu(II) bound to N-terminal and C-terminal residues. We show that A53T is predicted to cause increased β-sheet content of the peptide, with a persistent β-hairpin between residues 35-55 particularly notable. Phosphorylation has less effect on secondary structure but is predicted to significantly increase the size of the peptide, especially when bound to Cu(II), which is ascribed to reduced interaction of C-terminal sequence with central non-amyloid component. In addition, estimate of binding free energy to Cu(II) indicates A53T has little effect on metal-ion affinity, whereas phosphorylation markedly enhances the strength of binding. We suggest that the predicted changes in spatial extent and secondary structure of α-Synuclein may have implications for aggregation into Lewy bodies.

Simple Adjustment of Intranucleotide Base-Phosphate Interaction in the OL3 AMBER Force Field Improves RNA Simulations

Molecular dynamics (MD) simulations represent an established tool to study RNA molecules. The outcome of MD studies depends, however, on the quality of the force field (ff). Here we suggest a correction for the widely used AMBER OL3 ff by adding a simple adjustment of the nonbonded parameters. The reparameterization of the Lennard-Jones potential for the -H8···O5'- and -H6···O5'- atom pairs addresses an intranucleotide steric clash occurring in the type 0 base-phosphate interaction (0BPh). The nonbonded fix (NBfix) modification of 0BPh interactions (NBfix0BPh modification) was tuned via a reweighting approach and subsequently tested using an extensive set of standard and enhanced sampling simulations of both unstructured and folded RNA motifs. The modification corrects minor but visible intranucleotide clash for the anti nucleobase conformation. We observed that structural ensembles of small RNA benchmark motifs simulated with the NBfix0BPh modification provide better agreement with experiments. No side effects of the modification were observed in standard simulations of larger structured RNA motifs. We suggest that the combination of OL3 RNA ff and NBfix0BPh modification is a viable option to improve RNA MD simulations.

DNA groove preference shift upon phosphorylation of a protamine-like cationic peptide

Protamines, arginine-rich DNA-binding proteins, are responsible for chromatin compaction in sperm cells, but their DNA groove preference, major or minor, is not clearly identified. We herein study the DNA groove preference of a short protamine-like cationic peptide before and after phosphorylation, using all-atom molecular dynamics and umbrella sampling simulations. According to various thermodynamic and structural analyses, a peptide in its non-phosphorylated native state prefers the minor groove over the major groove, but phosphorylation of the peptide bound to the minor groove not only reduces its binding affinity but also brings a serious deformation of the minor groove, eliminating the minor-groove preference. As protamines are heavily phosphorylated before binding to DNA, we expect that the structurally disordered phosphorylated protamines would prefer major grooves to enter into DNA during spermatogenesis.

Machine Learning Subtle Conformational Change due to Phosphorylation in Intrinsically Disordered Proteins

Phosphorylation of intrinsically disordered proteins/regions (IDPs/IDRs) has a profound effect in biological functions such as cell signaling, protein folding or unfolding, and long-range allosteric effects. However, here we focus on two IDPs, namely 83-residue IDR transcription factor Ash1 and 92-residue long N-terminal region of CDK inhibitor Sic1 protein, found in Saccharomyces cerevisiae, for which experimental measurements of average conformational properties, namely, radius of gyration and structure factor, indicate negligible changes upon phosphorylation. Here, we show that a judicious dissection of conformational ensemble via combination of unsupervised machine learning and extensive molecular dynamics (MD) trajectories can highlight key differences and similarities among the phosphorylated and wild-type IDP. In particular, we develop Markov state model (MSM) using the latent-space dimensions of an autoencoder, trained using multi-microsecond long MD simulation trajectories. Examination of structural changes among the states, prior to and upon phosphorylation, captured several similarities and differences in their backbone contact maps, secondary structure, and torsion angles. Hydrogen bonding analysis revealed that phosphorylation not only increases the number of hydrogen bonds but also switches the pattern of hydrogen bonding between the backbone and side chain atoms with the phosphorylated residues. We also observe that although phosphorylation introduces salt bridges, there is a loss of the cation-π interaction. Phosphorylation also improved the probability for long-range hydrophobic contacts and also enhanced interaction with water molecules and improved the local structure of water as evident from the geometric order parameters. The observations on these machine-learnt states gave important insights, as it would otherwise be difficult to determine experimentally which is important, if we were to understand the role of phosphorylation of IDPs in their biological functions.

Modulation of Structure and Dynamics of Cardiac Troponin by Phosphorylation and Mutations Revealed by Molecular Dynamics Simulations

Adrenaline acts on β1 receptors in the heart muscle to enhance contractility, increase the heart rate, and increase the rate of relaxation (lusitropy) via activation of the cyclic AMP-dependent protein kinase, PKA. Phosphorylation of serines 22 and 23 in the N-terminal peptide of cardiac troponin I is responsible for lusitropy. Mutations associated with cardiomyopathy suppress the phosphorylation-dependent change. Key parts of troponin responsible for this modulatory system are disordered and cannot be resolved by conventional structural approaches. We performed all-atom molecular dynamics simulations (5 × 1.5 μs runs) of the troponin core (419 amino acids) in the presence of Ca2+ in the bisphosphorylated and unphosphorylated states for both wild-type troponin and the troponin C (cTnC) G159D mutant. PKA phosphorylation affects troponin dynamics. There is significant rigidification of the structure involving rearrangement of the cTnI(1-33)-cTnC interaction and changes in the distribution of the cTnC helix A/B angle, troponin I (cTnI) switch peptide (149-164) docking, and the angle between the regulatory head and ITC arm domains. The familial dilated cardiomyopathy cTnC G159D mutation whose Ca2+ sensitivity is not modulated by cTnI phosphorylation exhibits a structure inherently more rigid than the wild type, with phosphorylation reversing the direction of all metrics relative to the wild type.

Mono-phosphorylation at Ser4 of barrier-to-autointegration factor (Banf1) significantly reduces its DNA binding capability by inducing critical changes in its local conformation and DNA binding surface

Barrier-to-autointegration factor (Banf1) is a small DNA-bridging protein. The binding status of Banf1 to DNA is regulated by its N-terminal phosphorylation and dephosphorylation, which plays a critical role in cell proliferation. Banf1 can be phosphorylated at Ser4 into mono-phosphorylated Banf1, which is further phosphorylated at Thr3 to form di-phosphorylated Banf1. It was observed decades ago that mono-phosphorylated Banf1 cannot bind to DNA. However, the underlying molecular- and atomic-level mechanisms remain unclear. A clear understanding of these mechanisms will aid in interfering with the cell proliferation process for better global health. Herein, we explored the detailed atomic bases of unphosphorylated Banf1-DNA binding and how mono- and di-phosphorylation of Banf1 impair these atomic bases to eliminate its DNA-binding capability, followed by exploring the DNA-binding capability of mono- and di-phosphorylation Banf1, using comprehensive and systematic molecular modelling and molecular dynamics simulations. This work presented in detail the residue-level binding energies, hydrogen bonds and water bridges between Banf1 and DNA, some of which have not been reported. Moreover, we revealed that mono-phosphorylation of Banf1 causes its N-terminal secondary structure changes, which in turn induce significant changes in Banf1's DNA binding surface, thus eliminating its DNA-binding capability. At the atomic level, we also uncovered the alterations in interactions due to the induction of mono-phosphorylation that result in the N-terminal secondary structure changes of Banf1. Additionally, our modelling showed that phosphorylated Banf1 with their dominant N-terminal secondary structures bind to DNA with a significantly lower affinity and the docked binding pose are not stable in MD simulations. These findings help future studies in predicting effect of mutations in Banf1 on its DNA-binding capability and open a novel avenue for the development of therapeutics such as cancer drugs, targeting cell proliferation by inducing conformational changes in Banf1's N-terminal domain.

Histone H3 serine-57 is a CHK1 substrate whose phosphorylation affects DNA repair

Histone post-translational modifications promote a chromatin environment that controls transcription, DNA replication and repair, but surprisingly few phosphorylations have been documented. We report the discovery of histone H3 serine-57 phosphorylation (H3S57ph) and show that it is implicated in different DNA repair pathways from fungi to vertebrates. We identified CHK1 as a major human H3S57 kinase, and disrupting or constitutively mimicking H3S57ph had opposing effects on rate of recovery from replication stress, 53BP1 chromatin binding, and dependency on RAD52. In fission yeast, mutation of all H3 alleles to S57A abrogated DNA repair by both non-homologous end-joining and homologous recombination, while cells with phospho-mimicking S57D alleles were partly compromised for both repair pathways, presented aberrant Rad52 foci and were strongly sensitised to replication stress. Mechanistically, H3S57ph loosens DNA-histone contacts, increasing nucleosome mobility, and interacts with H3K56. Our results suggest that dynamic phosphorylation of H3S57 is required for DNA repair and recovery from replication stress, opening avenues for investigating the role of this modification in other DNA-related processes.

Molecular-dynamics simulations of macromolecular diffraction, part I: Preparation of protein crystal simulations

Molecular-dynamics (MD) simulations of protein crystals enable the prediction of structural and dynamical features of both the protein and the solvent components of macromolecular crystals, which can be validated against diffraction data from X-ray crystallographic experiments. The simulations have been useful for studying and predicting both Bragg and diffuse scattering in protein crystallography; however, the preparation is not yet automated and includes choices and tradeoffs that can impact the results. Here we examine some of the intricacies and consequences of the choices involved in setting up MD simulations of protein crystals for the study of diffraction data, and provide a recipe for preparing the simulations, packaged in an accompanying Jupyter notebook. This article and the accompanying notebook are intended to serve as practical resources for researchers wishing to put these models to work.

Modulating Secondary Structure Motifs Through Photo-Labile Peptide Staples

Peptide-protein interactions (PPIs) are facilitated by the well-defined three-dimensional structure of bioactive peptides, interesting compounds for the development of new therapeutic agents. Their secondary structure and thus their propensity to engage in PPIs can be influenced by the introduction of peptide staples on the side chains. In particular, light-controlled staples based on azobenzene photoswitches and their structural influence on helical peptides have been studied extensively. In contrast, photolabile staples bearing photocages as a structural key motif, have mainly been used to block supramolecular interactions. Their influence on the secondary structure of the target peptide is under-investigated. Thus, in this study we use a combination of spectroscopic techniques and in silico simulations to systematically study a series of helical peptides with varying length of the photo-labile staple to obtain a detailed insight into the structure-property relationship in such photoresponsive biomolecules.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Toward a Molecular Mechanism of Complementary RNA Duplexes Denaturation

RNA duplexes are relatively rare but play very important biological roles. As an end-product of template-based RNA replication, they also have key implications for hypothetical primitive forms of life. Unless they are specifically separated by enzymes, these duplexes denature upon a temperature increase. However, mechanistic and kinetic aspects of RNA (and DNA) duplex thermal denaturation remain unclear at the microscopic level. We propose an in silico strategy that probes the thermal denaturation of RNA duplexes and allows for an extensive conformational space exploration along a wide temperature range with atomistic precision. We show that this approach first accounts for the strong sequence and length dependence of the duplexes melting temperature, reproducing the trends seen in the experiments and predicted by nearest-neighbor models. The simulations are then instrumental at providing a molecular picture of the temperature-induced strand separation. The textbook canonical "all-or-nothing" two-state model, very much inspired by the protein folding mechanism, can be nuanced. We demonstrate that a temperature increase leads to significantly distorted but stable structures with extensive base-fraying at the extremities, and that the fully formed duplexes typically do not form around melting. The duplex separation therefore appears as much more gradual than commonly thought.

Deconstructing allostery by computational assessment of the binding determinants of allosteric PTP1B modulators

Fragment-based drug discovery is an established methodology for finding hit molecules that can be elaborated into lead compounds. However it is currently challenging to predict whether fragment hits that do not bind to an orthosteric site could be elaborated into allosteric modulators, as in these cases binding does not necessarily translate into a functional effect. We propose a workflow using Markov State Models (MSMs) with steered molecular dynamics (sMD) to assess the allosteric potential of known binders. sMD simulations are employed to sample protein conformational space inaccessible to routine equilibrium MD timescales. Protein conformations sampled by sMD provide starting points for seeded MD simulations, which are combined into MSMs. The methodology is demonstrated on a dataset of protein tyrosine phosphatase 1B ligands. Experimentally confirmed allosteric inhibitors are correctly classified as inhibitors, whereas the deconstructed analogues show reduced inhibitory activity. Analysis of the MSMs provide insights into preferred protein-ligand arrangements that correlate with functional outcomes. The present methodology may find applications for progressing fragments towards lead molecules in FBDD campaigns.

Computational Characterization of the Binding Properties of the HIV1-Neutralizing Antibody PG16 and Design of PG16-Derived CDRH3 Peptides

PG16 is a broadly neutralizing antibody that binds to the gp120 subunit of the HIV-1 Env protein. The major interaction site is formed by the unusually long complementarity determining region (CDR) H3. The CDRH3 residue Tyr100H is known to represent a tyrosine sulfation site; however, this modification is not present in the experimental complex structure of PG16 with full-length HIV-1 Env. To investigate the role of sulfation for this complex, we modeled the sulfation of Tyr100H and compared the dynamics and energetics of the modified and unmodified complex by molecular dynamics simulations at the atomic level. Our results show that sulfation does not affect the overall conformation of CDRH3, but still enhances gp120 interactions both at the site of modification and for the neighboring residues. This stabilization affects not only protein-protein contacts, but also the interactions between PG16 and the gp120 glycan shield. Furthermore, we also investigated whether PG16-CDRH3 is a suitable template for the development of peptide mimetics. For a peptide spanning residues 93-105 of PG16, we obtained an experimental EC₅₀ value of 3nm for the binding of gp120 to the peptide. This affinity can be enhanced by almost one order of magnitude by artificial disulfide bonding between residues 99 and 100F. In contrast, any truncation results in significantly lower affinity, suggesting that the entire peptide segment is involved in gp120 recognition. Given their high affinity, it should be possible to further optimize the PG16-derived peptides as potential inhibitors of HIV invasion.

Atomistic Picture of Opening-Closing Dynamics of DNA Holliday Junction Obtained by Molecular Simulations

Holliday junction (HJ) is a noncanonical four-way DNA structure with a prominent role in DNA repair, recombination, and DNA nanotechnology. By rearranging its four arms, HJ can adopt either closed or open state. With enzymes typically recognizing only a single state, acquiring detailed knowledge of the rearrangement process is an important step toward fully understanding the biological function of HJs. Here, we carried out standard all-atom molecular dynamics (MD) simulations of the spontaneous opening-closing transitions, which revealed complex conformational transitions of HJs with an involvement of previously unconsidered "half-closed" intermediates. Detailed free-energy landscapes of the transitions were obtained by sophisticated enhanced sampling simulations. Because the force field overstabilizes the closed conformation of HJs, we developed a system-specific modification which for the first time allows the observation of spontaneous opening-closing HJ transitions in unbiased MD simulations and opens the possibilities for more accurate HJ computational studies of biological processes and nanomaterials.

Sensitivity of the RNA Structure to Ion Conditions as Probed by Molecular Dynamics Simulations of Common Canonical RNA Duplexes

RNA molecules play a key role in countless biochemical processes. RNA interactions, which are of highly diverse nature, are determined by the fact that RNA is a highly negatively charged polyelectrolyte, which leads to intimate interactions with an ion atmosphere. Although RNA molecules are formally single-stranded, canonical (Watson-Crick) duplexes are key components of folded RNAs. A double-stranded (ds) RNA is also important for the design of RNA-based nanostructures and assemblies. Despite the fact that the description of canonical dsRNA is considered the least problematic part of RNA modeling, the imperfect shape and flexibility of dsRNA can lead to imbalances in the simulations of larger RNAs and RNA-containing assemblies. We present a comprehensive set of molecular dynamics (MD) simulations of four canonical A-RNA duplexes. Our focus was directed toward the characterization of the influence of varying ion concentrations and of the size of the solvation box. We compared several water models and four RNA force fields. The simulations showed that the A-RNA shape was most sensitive to the RNA force field, with some force fields leading to a reduced inclination of the A-RNA duplexes. The ions and water models played a minor role. The effect of the box size was negligible, and even boxes with a small fraction of the bulk solvent outside the RNA hydration sphere were sufficient for the simulation of the dsRNA.

Structural investigation of CDCA3-Cdh1 protein-protein interactions using in vitro studies and molecular dynamics simulation

The anaphase-promoting complex/cyclosome (APC/C) ubiquitin ligase and its cofactor, Cdh1, regulate the expression of several cell-cycle proteins and their functions during mitosis. Levels of the protein cell division cycle-associated protein 3 (CDCA3), which is functionally required for mitotic entry, are regulated by APC/C^Cdh1 . CDCA3 is an intrinsically disordered protein and contains both C-terminal KEN box and D-box recognition motifs, enabling binding to Cdh1. Our previous findings demonstrate that CDCA3 has a phosphorylation-dependent non-canonical ABBA-like motif within the linker region bridging these two recognition motifs and is required for efficient binding to Cdh1. Here, we sought to identify and further characterize additional residues that participate within this ABBA-like motif using detailed in vitro experiments and in silico modeling studies. We identified the role of H-bonds, hydrophobic and ionic interactions across the CDCA3 ABBA-like motif in the linker region between KEN and D-box motifs. This linker region adopts a well-defined structure when bound to Cdh1 in the presence of phosphorylation. Upon alanine mutation, the structure of this region is lost, leading to higher flexibility, and alteration in affinities due to binding to alternate sites on Cdh1. Our findings identify roles for the anchoring residues in the non-canonical ABBA-like motif to promote binding to the APC/C^Cdh1 and regulation of CDCA3 protein levels.

Structural predictions of protein-DNA binding: MELD-DNA

Structural, regulatory and enzymatic proteins interact with DNA to maintain a healthy and functional genome. Yet, our structural understanding of how proteins interact with DNA is limited. We present MELD-DNA, a novel computational approach to predict the structures of protein-DNA complexes. The method combines molecular dynamics simulations with general knowledge or experimental information through Bayesian inference. The physical model is sensitive to sequence-dependent properties and conformational changes required for binding, while information accelerates sampling of bound conformations. MELD-DNA can: (i) sample multiple binding modes; (ii) identify the preferred binding mode from the ensembles; and (iii) provide qualitative binding preferences between DNA sequences. We first assess performance on a dataset of 15 protein-DNA complexes and compare it with state-of-the-art methodologies. Furthermore, for three selected complexes, we show sequence dependence effects of binding in MELD predictions. We expect that the results presented herein, together with the freely available software, will impact structural biology (by complementing DNA structural databases) and molecular recognition (by bringing new insights into aspects governing protein-DNA interactions).

Virtual Screening of Hepatitis B Virus Pre-Genomic RNA as a Novel Therapeutic Target

The global burden imposed by hepatitis B virus (HBV) infection necessitates the discovery and design of novel antiviral drugs to complement existing treatments. One attractive and underexploited therapeutic target is ε, an ~85-nucleotide (nt) cis-acting regulatory stem-loop RNA located at the 3'- and 5'-ends of the pre-genomic RNA (pgRNA). Binding of the 5'-end ε to the viral polymerase protein (P) triggers two early events in HBV replication: pgRNA and P packaging and reverse transcription. Our recent solution nuclear magnetic resonance spectroscopy structure of ε permits structure-informed drug discovery efforts that are currently lacking for P. Here, we employ a virtual screen against ε using a Food and Drug Administration (FDA)-approved compound library, followed by in vitro binding assays. This approach revealed that the anti-hepatitis C virus drug Daclatasvir is a selective ε-targeting ligand. Additional molecular dynamics simulations demonstrated that Daclatasvir targets ε at its flexible 6-nt priming loop (PL) bulge and modulates its dynamics. Given the functional importance of the PL, our work supports the notion that targeting ε dynamics may be an effective anti-HBV therapeutic strategy.

Atomistic Insights into Organization of RNA-Loaded Lipid Nanoparticles

RNA-based therapies have shown promise in a wide range of applications, from cancer therapy, treatment of inherited diseases to vaccination. Encapsulation of RNA into ionizable lipid (IL) containing lipid nanoparticles (LNPs) has enabled its safe and targeted delivery. We present here the simulations of the self-assembly process of pH-sensitive RNA-carrying LNPs and their internal morphology. At low pH, the simulations confirm a lipid core encapsulating RNA in the hexagonal phase. Our all-atom and coarse-grained simulations show that an RNA molecule inside an LNP is protected from interactions with ions by being enveloped in the charged ILs. At neutral pH, representing the environment after LNP administration into human tissues, LNPs expelled most of the encapsulated RNA and water and formed separate bulk IL-rich and ordered the helper-lipid-rich phase. Helper lipids arranged themselves to be in contact with RNA or water. The presented models provide atomistic understanding of the LNP structure and open a way to investigate them in silico, varying the LNP composition or interacting with other biostructures aiming at increasing the efficiency of RNA-based medicine.

Phosphorylation of the Hsp90 Co-Chaperone Hop Changes its Conformational Dynamics and Biological Function

The molecular chaperones Hsp90 and Hsp70 and their regulatory co-chaperone Hop play a key role at the crossroads of the folding pathways of numerous client proteins by forming fine-tuned multiprotein complexes. Alterations of the biomolecules involved may functionally impact the chaperone machinery: here, we integrate simulations and experiments to unveil how Hop conformational fitness and interactions can be controlled by the perturbation of just one residue. Specifically, we unveil how mechanisms mediated by Hop residue Y354 control Hop open and closed states, which affect binding of Hsp70/Hsp90. Phosphorylation or mutation of Hop-Y354 are shown to favor structural ensembles that are indeed not optimal for stable interactions with Hsp90 and Hsp70. This disfavors cellular accumulation of the stringent Hsp90 clients glucocorticoid receptor and the viral tyrosine kinase v-Src, with detrimental effects on v-Src activity. Our results show how the post-translational modification of a specific residue in Hop provides a regulation mechanism for the larger chaperone complex of which it is part. In this framework, the effects of one single alteration are amplified at the cellular level through the perturbation of protein-interaction networks.

Uridylation of the histone mRNA stem-loop weakens binding interactions with SLBP while maintaining interactions with 3'hExo

Histone mRNA degradation is controlled by the unique 3' stem-loop of histone mRNA and the stem-loop binding protein (SLBP). As part of this process, the 3' stem-loop is trimmed by the histone-specific 3' exonuclease (3'hExo) and uridylated by the terminal uridylyl transferase 7 (TUT7), creating partially degraded intermediates with short uridylations. The role of these uridylations in degradation is not fully understood. Our work examines changes in the stability of the ternary complex created by trimming and uridylation of the stem-loop to better understand the role of this process in the histone mRNA life cycle. In this study, we used fluorescence polarization and electrophoretic mobility shift assays to demonstrate that both SLBP and 3'hExo can bind to uridylated and partially degraded stem-loop intermediates, although with lower affinity. We further characterized this complex by performing 1-µs molecular dynamics simulations using the AMBER force field and Nanoscale Molecular Dynamics (NAMD). These simulations show that while uridylation helps maintain the overall shape of the stem-loop, the combination of uridylation and dephosphorylation of the TPNK motif in SLBP disrupts key RNA-protein interactions. They also demonstrate that uridylation allows 3'hExo to maintain contact with the stem-loop after partial degradation and plays a role in disrupting key base pairs in partially degraded histone mRNA intermediates. Together, these experiments and simulations suggest that trimming by 3'hExo, uridylation, and SLBP dephosphorylation weakens both RNA-protein interactions and the stem-loop itself. Our results further elucidate the role of uridylation and SLBP dephosphorylation in the early stages of histone mRNA degradation.

Polarizable Molecular Block Model: Toward the Development of an Induced Dipole Force Field for DNA

For flexible and highly ionized macromolecules such as DNA, it is important to correctly evaluate the intramolecular polarization in an induced dipole force field. In a proposed polarizable molecular block (PMB) model, a large molecule is divided into several molecular blocks. The atomic charges of the blocks are optimized by using the respective electrostatic potentials (ESPs) on the molecular surface. By using the capped hydrogen removal operation, the total charge of the blocks is controlled exactly to have an integer charge. The atomic polarizabilities of the blocks are optimized by using the respective polarized one-electron potentials that are the differences between ESPs with and without an external test charge. Induced dipole-charge interactions between the blocks are all included, but those interactions within the blocks are strictly excluded. All dipole-dipole interactions are included, but the damping functions are applied to the close dipole-dipole pairs. Several types of damping (simple scaling, exponential, linear, and Gaussian) are evaluated. The validity of the PMB model was verified by using trinucleotide duplexes which have A-, B-, and Z-DNA forms. The reference energies of trinucleotide duplexes including counterions (GGT3Na-ACC3Na, GAC3Na-GTC3Na, and GCG3Na-CGC3Na) are calculated using ωB97XD/aug-cc-pVDZ. All damping types reproduced well the reference interaction energies, dipole moments, and ESPs. Among them, the simple scaling with strong attenuation to 1-2 atomic pairs showed the highest stability against the polarization catastrophe. This study shows that it is possible to develop a high-quality polarizable force field by treating the intramolecular polarization on a block-by-block basis.

Effect of phosphorylation of protamine-like cationic peptide on the binding affinity to DNA

Protamines are more arginine-rich and more basic than histones and are responsible for providing a highly compacted shape to the sperm heads in the testis. Phosphorylation and dephosphorylation are two events that occur in the late phase of spermatogenesis before the maturation of sperms. In this work, we have studied the effect of phosphorylation of protamine-like cationic peptides using all-atom molecular dynamics simulations. Through thermodynamic analyses, we found that phosphorylation reduces the binding efficiency of such cationic peptides on DNA duplexes. Peptide phosphorylation leads to a less efficient DNA condensation, due to a competition between DNA-peptide and peptide-peptide interactions. We hypothesize that the decrease of peptide bonds between DNA together with peptide self-assembly might allow an optimal re-organization of chromatin and an efficient condensation through subsequent peptide dephosphorylation. Based on the globular and compact conformations of phosphorylated peptides mediated by arginine-phosphoserine H-bonding, we furthermore postulate that phosphorylated protamines could more easily intrude into chromatin and participate to histone release through disruption of histone-histone and histone-DNA binding during spermatogenesis.

Design, Synthesis, and Characterization of a Novel 2'-5'-Linked Amikacin-Binding Aptamer: An Experimental and MD Simulation Study

To extend the approach of using RNA aptamers as transient protective groups for the synthesis of novel small-molecule drug derivatives from the existing aminoglycosides, we incorporated 2'-5' phosphodiester backbone modification in a structurally known neomycin RNA aptamer and studied the binding of a series of aminoglycosides using isothermal calorimetry (ITC) and molecular dynamics (MD) simulation. Experimental characterization of amikacin, a commercially available and widely used aminoglycoside for treating bacterial infections, shows that the aptamer A1 with a 2'-5' linkage between G15 and U16 exhibits a sevenfold increase in binding affinity with a lower binding energy compared to the native aptamer. Molecular dynamics (MD) simulation studies rationalize that this noncanonical linkage generates a narrower binding pocket by creating a superspiral RNA helical structure, which improves the ligand's fit in the binding pocket. These results provide new insights into applying 2'-5' linkages to diversify functional RNA aptamers as noncovalent protective groups in the synthesis of aminoglycoside derivatives, which can be further extended to other current drug molecules and complex natural compounds to make new pools of drug candidates more efficiently.

Dynamics and Allosteric Information Pathways of Unphosphorylated c-Cbl

Human c-Cbl is a RING-type ligase and plays a central role in the protein degradation cascade. To elucidate its conformational changes related to substrate binding, we performed molecular dynamics simulations of different variants/states of c-Cbl for a cumulative time of 68 μs. Our simulations demonstrate that before the substrate binds, the RING domain samples a broad set of conformational states at a biologically relevant salt concentration, including the closed, partially open, and fully open states, whereas substrate binding leads to a restricted conformational sampling. Phe378 and the C-terminal region play an essential role in stabilizing the partially open state. To visualize the allosteric signal transmission pathways from the substrate-binding site to the 40 Å apart RING domain and identify the critical residues for allostery, we have created a subgraph from the optimal and suboptimal paths. Redundant paths are seen in the SH2 domain where the substrate binds, while the major bottlenecks are found at the junction between the SH2 domain and the linker helix region as well as that between the SH2 domain and the 4H bundle. These bottlenecks separate the paths into two overall routes. The nodes/residues at the bottlenecks on the subgraph are considered allosteric hot spots. This subgraph approach provides a general tool for network visualization and determination of critical residues for allostery. The structurally and allosterically critical residues identified in our work are testable and would provide valuable insights into the emerging strategies for drug discovery, such as targeted protein degradation.

Decoding the Identification Mechanism of an SAM-III Riboswitch on Ligands through Multiple Independent Gaussian-Accelerated Molecular Dynamics Simulations

S-Adenosyl-l-methionine (SAM)-responsive riboswitches play a central role in the regulation of bacterial gene expression at the level of transcription attenuation or translation inhibition. In this study, multiple independent Gaussian-accelerated molecular dynamics simulations were performed to decipher the identification mechanisms of SAM-III (S_MK) on ligands SAM, SAH, and EEM. The results reveal that ligand binding highly affects the structural flexibility, internal dynamics, and conformational changes of SAM-III. The dynamic analysis shows that helices P3 and P4 as well as two junctions J23 and J24 of SAM-III are highly susceptible to ligand binding. Analyses of free energy landscapes suggest that ligand binding induces different free energy profiles of SAM-III, which leads to the difference in identification sites of SAM-III on ligands. The information on ligand-nucleotide interactions not only uncovers that the π-π, cation-π, and hydrogen bonding interactions drive identification of SAM-III on the three ligands but also reveals that different electrostatic properties of SAM, SAH, and EEM alter the active sites of SAM-III. Meanwhile, the results also verify that the adenine group of SAM, SAH, and EEM is well recognized by conserved nucleotides G7, A29, U37, A38, and G48. We expect that this study can provide useful information for understanding the applications of SAM-III in chemical, synthetic RNA biology, and biomedical fields.

Multiscale Modeling of Phosphate···π Contacts in RNA U-Turns Exposes Differences between Quantum-Chemical and AMBER Force Field Descriptions

Phosphate···π, also called anion···π, contacts occur between nucleobases and anionic phosphate oxygens (OP2) in r(GNRA) and r(UNNN) U-turn motifs (N = A,G,C,U; R = A,G). These contacts were investigated using state-of-the-art quantum-chemical methods (QM) to characterize their physicochemical properties and to serve as a reference to evaluate AMBER force field (AFF) performance. We found that phosphate···π interaction energies calculated with the AFF for dimethyl phosphate···nucleobase model systems are less stabilizing in comparison with double-hybrid DFT and that minimum contact distances are larger for all nucleobases. These distance stretches are also observed in large-scale AFF vs QM/MM computations and classical molecular dynamics (MD) simulations on several r(gcGNRAgc) tetraloop hairpins when compared to experimental data extracted from X-ray/cryo-EM structures (res. ≤ 2.5 Å) using the WebFR3D bioinformatic tool. MD simulations further revealed shifted OP2/nucleobase positions. We propose that discrepancies between the QM and AFF result from a combination of missing polarization in the AFF combined with too large AFF Lennard-Jones (LJ) radii of nucleobase carbon atoms in addition to an exaggerated short-range repulsion of the r^-12 LJ repulsive term. We compared these results with earlier data gathered on lone pair···π contacts in CpG Z-steps occurring in r(UNCG) tetraloops. In both instances, charge transfer calculations do not support any significant n → π* donation effects. We also investigated thiophosphate···π contacts that showed reduced stabilizing interaction energies when compared to phosphate···π contacts. Thus, we challenge suggestions that the experimentally observed enhanced thermodynamic stability of phosphorothioated r(GNRA) tetraloops can be explained by larger London dispersion.

Allosteric rescue of catalytically impaired ATP phosphoribosyltransferase variants links protein dynamics to active-site electrostatic preorganisation

ATP phosphoribosyltransferase catalyses the first step of histidine biosynthesis and is controlled via a complex allosteric mechanism where the regulatory protein HisZ enhances catalysis by the catalytic protein HisG_S while mediating allosteric inhibition by histidine. Activation by HisZ was proposed to position HisG_S Arg56 to stabilise departure of the pyrophosphate leaving group. Here we report active-site mutants of HisG_S with impaired reaction chemistry which can be allosterically restored by HisZ despite the HisZ:HisG_S interface lying ~20 Å away from the active site. MD simulations indicate HisZ binding constrains the dynamics of HisG_S to favour a preorganised active site where both Arg56 and Arg32 are poised to stabilise leaving-group departure in WT-HisG_S. In the Arg56Ala-HisG_S mutant, HisZ modulates Arg32 dynamics so that it can partially compensate for the absence of Arg56. These results illustrate how remote protein-protein interactions translate into catalytic resilience by restoring damaged electrostatic preorganisation at the active site.

Spontaneous binding of single-stranded RNAs to RRM proteins visualized by unbiased atomistic simulations with a rescaled RNA force field

Recognition of single-stranded RNA (ssRNA) by RNA recognition motif (RRM) domains is an important class of protein-RNA interactions. Many such complexes were characterized using nuclear magnetic resonance (NMR) and/or X-ray crystallography techniques, revealing ensemble-averaged pictures of the bound states. However, it is becoming widely accepted that better understanding of protein-RNA interactions would be obtained from ensemble descriptions. Indeed, earlier molecular dynamics simulations of bound states indicated visible dynamics at the RNA-RRM interfaces. Here, we report the first atomistic simulation study of spontaneous binding of short RNA sequences to RRM domains of HuR and SRSF1 proteins. Using a millisecond-scale aggregate ensemble of unbiased simulations, we were able to observe a few dozen binding events. HuR RRM3 utilizes a pre-binding state to navigate the RNA sequence to its partially disordered bound state and then to dynamically scan its different binding registers. SRSF1 RRM2 binding is more straightforward but still multiple-pathway. The present study necessitated development of a goal-specific force field modification, scaling down the intramolecular van der Waals interactions of the RNA which also improves description of the RNA-RRM bound state. Our study opens up a new avenue for large-scale atomistic investigations of binding landscapes of protein-RNA complexes, and future perspectives of such research are discussed.

Structural dynamics of the N‐terminal SH2 domain of PI3K in its free and CD28 ‐bound states

Protein conformational changes with fluctuations are fundamental aspects of protein-protein interactions (PPIs); understanding these motions is required for the rational design of PPI-regulating compounds. Src homology 2 (SH2) domains are commonly found in adapter proteins involved in signal transduction and specifically bind to consensus motifs of proteins containing phosphorylated tyrosine (pY). Here, we analysed the interaction between the N-terminal SH2 domain (nSH2) of the regulatory subunit in phosphoinositide 3-kinase (PI3K) and the cytoplasmic region of the T-cell co-receptor, CD28, using NMR and molecular dynamics (MD) simulations. First, we assigned the backbone signals of nSH2 on ¹ H-¹⁵ N heteronuclear single quantum coherence spectra in the absence or presence of the CD28 phosphopeptide, SDpYMNMTPRRPG. Chemical shift perturbation experiments revealed allosteric changes at the BC loop and the C-terminal region of nSH2 upon CD28 binding. NMR relaxation experiments showed a conformational exchange associated with CD28 binding in these regions. The conformational stabilisation of the C-terminal region correlated with the regulation of PI3K catalytic function. Further, using ¹⁹ F- and ³¹ P-labelled CD28 phosphopeptide, we analysed the structural dynamics of CD28 and demonstrated that the aromatic ring of the pY residue fluctuated between multiple conformations upon nSH2 binding. Our MD simulations largely explained the NMR results and the structural dynamics of nSH2 and CD28 in both bound and unbound states. Notably, in addition to its major conformation, we detected a minor conformation of nSH2 in the CD28 bound state that may explain the allosteric conformational change in the BC loop.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg²⁺ -free state characterized by shortened β2/β3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Molecular Dynamics Simulations of Protein RNA Complexes by Using an Advanced Electrostatic Model

Protein-RNA interactions are integral to the biological functions of RNA. It is well recognized that molecular dynamics (MD) simulations of protein-RNA complexes are more challenging than those of each component. The difficulty arises from the strong electrostatic interactions and the delicate balance between various types of physical forces at the interface. Previously, MD simulations of protein-RNA complexes have predominantly employed fixed-charge force fields. Although force field modifications have been developed to address problems identified in the simulations, some protein-RNA structures are still hard to reproduce by simulations. Here, we present MD simulations of two representative protein-RNA complexes using the AMOEBA polarizable force field. The van der Waals parameters were refined to reproduce accurate quantum-mechanical data of base-base and base-amino acid interactions. It was found that the refined parameters produced a more stable hydrogen-bond network in the interface. One of the complexes remained stable during the short simulations, whereas it could quickly break down in previous simulations using fixed-charge force fields. There was reversible breaking and formation of hydrogen bonds that are observed in the crystal structure, which may indicate the difference in solution and crystal structures. While further improvement and validation of the force fields are still needed, this work demonstrates that polarizable force fields are promising for the study of protein-RNA complexes.

Substrate spectrum of PPM1D in the cellular response to DNA double-strand breaks

PPM1D is a p53-regulated protein phosphatase that modulates the DNA damage response (DDR) and is frequently altered in cancer. Here, we employed chemical inhibition of PPM1D and quantitative mass spectrometry-based phosphoproteomics to identify the substrates of PPM1D upon induction of DNA double-strand breaks (DSBs) by etoposide. We identified 73 putative PPM1D substrates that are involved in DNA repair, regulation of transcription, and RNA processing. One-third of DSB-induced S/TQ phosphorylation sites are dephosphorylated by PPM1D, demonstrating that PPM1D only partially counteracts ATM/ATR/DNA-PK signaling. PPM1D-targeted phosphorylation sites are found in a specific amino acid sequence motif that is characterized by glutamic acid residues, high intrinsic disorder, and poor evolutionary conservation. We identified a functionally uncharacterized protein Kanadaptin as ATM and PPM1D substrate upon DSB induction. We propose that PPM1D plays a role during the response to DSBs by regulating the phosphorylation of DNA- and RNA-binding proteins in intrinsically disordered regions.

•

MS-based phosphoproteomic profiling of PPM1D substrates in U2OS and HCT116 cells

•

PPM1D counteracts ATM in the cellular response to DNA double-strand breaks

•

PPM1D target sites localize to glutamic acid-rich regions with high intrinsic disorder

•

Kanadaptin is a putative DNA damage response factor regulated by ATM and PPM1D

Molecular Insights into the Regulation of 3-Phosphoinositide-Dependent Protein Kinase 1: Modeling the Interaction between the Kinase and the Pleckstrin Homology Domains

The 3-phosphoinositide-dependent protein kinase 1 (PDK1) K465E mutant kinase can still activate protein kinase B (PKB) at the membrane in a phosphatidylinositol-3,4,5-trisphosphate (PIP₃, PtdIns(3,4,5)P₃) independent manner. To understand this new PDK1 regulatory mechanism, docking and molecular dynamics calculations were performed for the first time to simulate the wild-type kinase domain-pleckstrin homology (PH) domain complex with PH-in and PH-out conformations. These simulations were then compared to the PH-in model of the KD-PH(mutant K465E) PDK1 complex. Additionally, three KD-PH complexes were simulated, including a substrate analogue bound to a hydrophobic pocket (denominated the PIF-pocket) substrate-docking site. We find that only the PH-out conformation, with the PH domain well-oriented to interact with the cellular membrane, is active for wild-type PDK1. In contrast, the active conformation of the PDK1 K465E mutant is PH-in, being ATP-stable at the active site while the PIF-pocket is more accessible to the peptide substrate. We corroborate that both the docking-site binding and the catalytic activity are in fact enhanced in knock-in mouse samples expressing the PDK1 K465E protein, enabling the phosphorylation of PKB in the absence of PIP₃ binding.

Structure of a 14-3-3ε:FOXO3apS253 Phosphopeptide Complex Reveals 14-3-3 Isoform-Specific Binding of Forkhead Box Class O Transcription Factor (FOXO) Phosphoproteins

The transcriptional activity of Forkhead Box O3 (FOXO3a) is inactivated by AKT-mediated phosphorylation on Serine 253 (S253), which enables FOXO3a binding to 14-3-3. Phosphorylated FOXO3a binding to 14-3-3 facilitates the nuclear exclusion of FOXO3a, causing cancer cell proliferation. The FOXO3a/14-3-3 interaction has, therefore, emerged as an important therapeutic target. Here, we report a comprehensive analysis using fluorescence polarization, isothermal titration calorimetry, small-angle X-ray scattering, X-ray crystallography, and molecular dynamics simulations to gain molecular-level insights into the interaction of FOXO3a^pS253 phosphopeptide with 14-3-3ε. A high-resolution structure of the fluorophore-labeled FOXO3a^pS253:14-3-3ε complex revealed a distinct mode of interaction compared to other 14-3-3 phosphopeptide complexes. FOXO3a^pS253 phosphopeptide showed significant structural difference in the positions of the -3 and -4 Arg residues relative to pSer, compared to that of a similar phosphopeptide, FOXO1^pS256 bound to 14-3-3σ. Moreover, molecular dynamics studies show that the significant structural changes and molecular interactions noticed in the crystal structure of FOXO3a^pS253:14-3-3ε are preserved over the course of the simulation. Thus, this study reveals structural differences between the binding to 14-3-3 isoforms of FOXO1^pS256 versus FOXO3a^pS253, providing a framework for the rational design of isoform-specific FOXO/14-3-3 protein-protein interaction inhibitors for therapy.

Conformational heterogeneity of UCAAUC RNA oligonucleotide from molecular dynamics simulations, SAXS, and NMR experiments

We describe the conformational ensemble of the single-stranded r(UCAAUC) oligonucleotide obtained using extensive molecular dynamics (MD) simulations and Rosetta's FARFAR2 algorithm. The conformations observed in MD consist of A-form-like structures and variations thereof. These structures are not present in the pool generated using FARFAR2. By comparing with available nuclear magnetic resonance (NMR) measurements, we show that the presence of both A-form-like and other extended conformations is necessary to quantitatively explain experimental data. To further validate our results, we measure solution X-ray scattering (SAXS) data on the RNA hexamer and find that simulations result in more compact structures than observed from these experiments. The integration of simulations with NMR via a maximum entropy approach shows that small modifications to the MD ensemble lead to an improved description of the conformational ensemble. Nevertheless, we identify persisting discrepancies in matching experimental SAXS data.