Competitive and Cooperative Interactions Mediate RNA Transfer from Herpesvirus Saimiri ORF57 to the Mammalian Export Adaptor ALYREF

The essential herpesvirus adaptor protein HVS ORF57, which has homologs in all other herpesviruses, promotes viral mRNA export by utilizing the cellular mRNA export machinery. ORF57 protein specifically recognizes viral mRNA transcripts, and binds to proteins of the cellular transcription-export (TREX) complex, in particular ALYREF. This interaction introduces viral mRNA to the NXF1 pathway, subsequently directing it to the nuclear pore for export to the cytoplasm. Here we have used a range of techniques to reveal the sites for direct contact between RNA and ORF57 in the absence and presence of ALYREF. A binding site within ORF57 was characterized which recognizes specific viral mRNA motifs. When ALYREF is present, part of this ORF57 RNA binding site, composed of an a-helix, binds preferentially to ALYREF. This competitively displaces viral RNA from the a-helix, but contact with RNA is still maintained by a flanking region. At the same time, the flexible N-terminal domain of ALYREF comes into contact with the viral RNA, which becomes engaged in an extensive network of synergistic interactions with both ALYREF and ORF57. Transfer of RNA to ALYREF in the ternary complex, and involvement of individual ORF57 residues in RNA recognition, were confirmed by UV cross-linking and mutagenesis. The atomic-resolution structure of the ORF57-ALYREF interface was determined, which noticeably differed from the homologous ICP27-ALYREF structure. Together, the data provides the first site-specific description of how viral mRNA is locked by a herpes viral adaptor protein in complex with cellular ALYREF, giving herpesvirus access to the cellular mRNA export machinery. The NMR strategy used may be more generally applicable to the study of fuzzy protein-protein-RNA complexes which involve flexible polypeptide regions

In Silico Elucidation of the Recognition Dynamics of Ubiquitin

Elucidation of the mechanism of biomacromolecular recognition events has been a topic of intense interest over the past century. The inherent dynamic nature of both protein and ligand molecules along with the continuous reshaping of the energy landscape during the binding process renders it difficult to characterize this process at atomic detail. Here, we investigate the recognition dynamics of ubiquitin via microsecond all-atom molecular dynamics simulation providing both thermodynamic and kinetic information. The high-level of consistency found with respect to experimental NMR data lends support to the accuracy of the in silico representation of the conformational substates and their interconversions of free ubiquitin. Using an energy-based reweighting approach, the statistical distribution of conformational states of ubiquitin is monitored as a function of the distance between ubiquitin and its binding partner Hrs-UIM. It is found that extensive and dense sampling of conformational space afforded by the µs MD trajectory is essential for the elucidation of the binding mechanism as is Boltzmann sampling, overcoming inherent limitations of sparsely sampled empirical ensembles. The results reveal a population redistribution mechanism that takes effect when the ligand is at intermediate range of 1–2 nm from ubiquitin. This mechanism, which may be depicted as a superposition of the conformational selection and induced fit mechanisms, also applies to other binding partners of ubiquitin, such as the GGA3 GAT domain

Unique Structure and Dynamics of the EphA5 Ligand Binding Domain Mediate Its Binding Specificity as Revealed by X-ray Crystallography, NMR and MD Simulations

10.1371/journal.pone.0074040PLoS ONE89-POLN

Atomic-Level Characterization of the Activation Mechanism of SERCA by Calcium

We have performed molecular dynamics (MD) simulations to elucidate, in atomic detail, the mechanism by which the sarcoplasmic reticulum Ca2+-ATPase (SERCA) is activated by Ca2+. Crystal structures suggest that activation of SERCA occurs when the cytoplasmic head-piece, in an open (E1) conformation stabilized by Ca2+, undergoes a large-scale open-to-closed (E1 to E2) transition that is induced by ATP binding. However, spectroscopic measurements in solution suggest that these structural states (E1 and E2) are not tightly coupled to biochemical states (defined by bound ligands); the closed E2 state predominates even in the absence of ATP, in both the presence and absence of Ca2+. How is this loose coupling consistent with the high efficiency of energy transduction in the Ca2+-ATPase? To provide insight into this question, we performed long (500 ns) all-atom MD simulations starting from the open crystal structure, including a lipid bilayer and water. In both the presence and absence of Ca2+, we observed a large-scale open-to-closed conformational transition within 400 ns, supporting the weak coupling between structural and biochemical states. However, upon closer inspection, it is clear that Ca2+ is necessary and sufficient for SERCA to reach the precise geometrical arrangement necessary for activation of ATP hydrolysis. Contrary to suggestions from crystal structures, but in agreement with solution spectroscopy, the presence of ATP is not required for this activating transition. Principal component analysis showed that Ca2+ reshapes the free energy landscape of SERCA to create a path between the open conformation and the activated closed conformation. Thus the malleability of the free energy landscape is essential for SERCA efficiency, ensuring that ATP hydrolysis is tightly coupled to Ca2+ transport. These results demonstrate the importance of real-time dynamics in the formation of catalytically competent conformations of SERCA, with broad implications for understanding enzymatic catalysis in atomic detail

Perturbation-Response Scanning Reveals Ligand Entry-Exit Mechanisms of Ferric Binding Protein

We study apo and holo forms of the bacterial ferric binding protein (FBP) which exhibits the so-called ferric transport dilemma: it uptakes iron from the host with remarkable affinity, yet releases it with ease in the cytoplasm for subsequent use. The observations fit the “conformational selection” model whereby the existence of a weakly populated, higher energy conformation that is stabilized in the presence of the ligand is proposed. We introduce a new tool that we term perturbation-response scanning (PRS) for the analysis of remote control strategies utilized. The approach relies on the systematic use of computational perturbation/response techniques based on linear response theory, by sequentially applying directed forces on single-residues along the chain and recording the resulting relative changes in the residue coordinates. We further obtain closed-form expressions for the magnitude and the directionality of the response. Using PRS, we study the ligand release mechanisms of FBP and support the findings by molecular dynamics simulations. We find that the residue-by-residue displacements between the apo and the holo forms, as determined from the X-ray structures, are faithfully reproduced by perturbations applied on the majority of the residues of the apo form. However, once the stabilizing ligand (Fe) is integrated to the system in holo FBP, perturbing only a few select residues successfully reproduces the experimental displacements. Thus, iron uptake by FBP is a favored process in the fluctuating environment of the protein, whereas iron release is controlled by mechanisms including chelation and allostery. The directional analysis that we implement in the PRS methodology implicates the latter mechanism by leading to a few distant, charged, and exposed loop residues. Upon perturbing these, irrespective of the direction of the operating forces, we find that the cap residues involved in iron release are made to operate coherently, facilitating release of the ion

Mathematical treatment of adiabatic fast passage pulses for the computation of nuclear spin relaxation rates in proteins with conformational exchange

Although originally designed for broadband inversion and decoupling in NMR spectroscopy, recent methodological developments have introduced adiabatic fast passage (AFP) pulses into the field of protein dynamics. AFP pulses employ a frequency sweep, and have not only superior inversion properties with respect to offset effects, but they are also easily implemented into a pulse sequence. As magnetization is dragged from the +z to the −z direction, Larmor precession is impeded since magnetization becomes spin-locked, which is a potentially useful feature for the investigation of microsecond to millisecond dynamics. A major drawback of these pulses as theoretical prediction is concerned, however, results from their time-dependent offset: simulations of spin density matrices under the influence of a time-dependent Hamiltonian with non-commuting elements are costly in terms of computational time, rendering data analysis impracticable. In this paper we suggest several ways to reduce the computational time without compromising accuracy with respect to effects such as cross-correlated relaxation and modulation of the chemical shift

PPARγ population shift produces disease-related changes in molecular networks associated with metabolic syndrome

Peroxisome proliferator-activated receptor gamma (PPARγ) is a key regulator of adipocyte differentiation and has an important role in metabolic syndrome. Phosphorylation of the receptor's ligand-binding domain at serine 273 has been shown to change the expression of a large number of genes implicated in obesity. The difference in gene expression seen when comparing wild-type phosphorylated with mutant non-phosphorylated PPARγ may have important consequences for the cellular molecular network, the state of which can be shifted from the healthy to a stable diseased state. We found that a group of differentially expressed genes are involved in bi-stable switches and form a core network, the state of which changes with disease progression. These findings support the idea that bi-stable switches may be a mechanism for locking the core gene network into a diseased state and for efficiently propagating perturbations to more distant regions of the network. A structural analysis of the PPARγ–RXRα dimer complex supports the hypothesis of a major structural change between the two states, and this may represent an important mechanism leading to the differential expression observed in the core network

Computer-Based Screening of Functional Conformers of Proteins

A long-standing goal in biology is to establish the link between function, structure, and dynamics of proteins. Considering that protein function at the molecular level is understood by the ability of proteins to bind to other molecules, the limited structural data of proteins in association with other bio-molecules represents a major hurdle to understanding protein function at the structural level. Recent reports show that protein function can be linked to protein structure and dynamics through network centrality analysis, suggesting that the structures of proteins bound to natural ligands may be inferred computationally. In the present work, a new method is described to discriminate protein conformations relevant to the specific recognition of a ligand. The method relies on a scoring system that matches critical residues with central residues in different structures of a given protein. Central residues are the most traversed residues with the same frequency in networks derived from protein structures. We tested our method in a set of 24 different proteins and more than 260,000 structures of these in the absence of a ligand or bound to it. To illustrate the usefulness of our method in the study of the structure/dynamics/function relationship of proteins, we analyzed mutants of the yeast TATA-binding protein with impaired DNA binding. Our results indicate that critical residues for an interaction are preferentially found as central residues of protein structures in complex with a ligand. Thus, our scoring system effectively distinguishes protein conformations relevant to the function of interest

Histidine Hydrogen-Deuterium Exchange Mass Spectrometry for Probing the Microenvironment of Histidine Residues in Dihydrofolate Reductase

Histidine Hydrogen-Deuterium Exchange Mass Spectrometry (His-HDX-MS) determines the HDX rates at the imidazole C(2)-hydrogen of histidine residues. This method provides not only the HDX rates but also the pK(a) values of histidine imidazole rings. His-HDX-MS was used to probe the microenvironment of histidine residues of E. coli dihydrofolate reductase (DHFR), an enzyme proposed to undergo multiple conformational changes during catalysis.Using His-HDX-MS, the pK(a) values and the half-lives (t(1/2)) of HDX reactions of five histidine residues of apo-DHFR, DHFR in complex with methotrexate (DHFR-MTX), DHFR in complex with MTX and NADPH (DHFR-MTX-NADPH), and DHFR in complex with folate and NADP+ (DHFR-folate-NADP+) were determined. The results showed that the two parameters (pK(a) and t(1/2)) are sensitive to the changes of the microenvironment around the histidine residues. Although four of the five histidine residues are located far from the active site, ligand binding affected their pK(a), t(1/2) or both. This is consistent with previous observations of ligand binding-induced distal conformational changes on DHFR. Most of the observed pK(a) and t(1/2) changes could be rationalized using the X-ray structures of apo-DHFR, DHFR-MTX-NADPH, and DHFR-folate-NADP+. The availability of the neutron diffraction structure of DHFR-MTX enabled us to compare the protonation states of histidine imidazole rings.Our results demonstrate the usefulness of His-HDX-MS in probing the microenvironments of histidine residues within proteins

The Mechanism of Ubiquitination in the Cullin-RING E3 Ligase Machinery: Conformational Control of Substrate Orientation

In cullin-RING E3 ubiquitin ligases, substrate binding proteins, such as VHL-box, SOCS-box or the F-box proteins, recruit substrates for ubiquitination, accurately positioning and orienting the substrates for ubiquitin transfer. Yet, how the E3 machinery precisely positions the substrate is unknown. Here, we simulated nine substrate binding proteins: Skp2, Fbw7, β-TrCP1, Cdc4, Fbs1, TIR1, pVHL, SOCS2, and SOCS4, in the unbound form and bound to Skp1, ASK1 or Elongin C. All nine proteins have two domains: one binds to the substrate; the other to E3 ligase modules Skp1/ASK1/Elongin C. We discovered that in all cases the flexible inter-domain linker serves as a hinge, rotating the substrate binding domain, optimally and accurately positioning it for ubiquitin transfer. We observed a conserved proline in the linker of all nine proteins. In all cases, the prolines pucker substantially and the pucker is associated with the backbone rotation toward the E2/ubiquitin. We further observed that the linker flexibility could be regulated allosterically by binding events associated with either domain. We conclude that the flexible linker in the substrate binding proteins orients the substrate for the ubiquitin transfer. Our findings provide a mechanism for ubiquitination and polyubiquitination, illustrating that these processes are under conformational control