Molecular dynamic studies on the impact of mutations on the structure,stability, and N‐terminal orientation of annexin A1: Implications for membrane aggregation

Multiple MD simulations were performed for the full‐length wild‐type A1, the full length A1 mutations S27E and S27A, as well as the N‐terminal peptide (AMVSEFLKQAWFIDNEEQEYIKTVKG S ²⁷ KGGPGSAVSPYPTFN) of wild‐type A1 and mutations S27E and S27A. The MD simulation trajectories of about 350 ns were generated and analyzed to examine the changes of core domain calcium binding affinity, core domain and N‐terminal domain structures, and N‐terminal domain orientation. Our results indicated that S27A and S27E mutations caused little changes on the calcium‐binding affinity of the core domain of A1. However, the S27A mutation made the N‐terminal domain of A1 less helical, and made the N‐terminal domain migrate faster toward the core domain; these impacts on A1 are beneficial to the membrane aggregation process. On the contrary, the S27E mutation made the N‐terminal domain of A1 more stable, and hindered the migration to the core domain; these changes on A1 are antagonistic for the membrane aggregation process. Our results using MD simulations provide an atomistic explanation for experimental observations that the S27E mutant showed a higher calcium concentration requirement and lower maximal extent of aggregation, while the wild‐type and two mutants S27E and S27A required identical calcium concentrations for liposome binding. Proteins 2014; 82:3327–3334. © 2014 Wiley Periodicals, Inc.     

Computational modeling and molecular dynamics simulations of mammalian cytoplasmic tyrosyl-tRNA synthetase and its complexes with substrates

Cytoplasmic tyrosyl-tRNA synthetase (TyrRS) is one of the key enzymes of protein biosynthesis. TyrRSs of pathogenic organisms have gained attention as potential targets for drug development. Identifying structural differences between various TyrRSs will facilitate the development of specific inhibitors for the TyrRSs of pathogenic organisms. However, there is a deficiency in structural data for mammalian cytoplasmic TyrRS in complexes with substrates. In this work, we constructed spatial structure of full-length Bos taurus TyrRS (BtTyrRS) and its complexes with substrates using the set of computational modeling techniques. Special attention was paid to BtTyrRS complexes with substrates [L-tyrosine, K⁺ and ATP:Mg²⁺] and intermediate products [tyrosyl-adenylate (Tyr-AMP), K⁺ and PP_i:Mg²⁺] with the different catalytic loop conformations. In order to analyze their dynamical properties, we performed 100 ns of molecular dynamics (MD) simulations. MD simulations revealed new structural data concerning the tyrosine activation reaction in mammalian TyrRS. Formation of strong interaction between Lys154 and γ-phosphate suggests the additional role of CP1 insertion as an important factor for ATP binding. The presence of a potassium-binding pocket within the active site of mammalian TyrRS compensates the absence of the second lysine in the KMSKS motif. Our data provide new details concerning a role of K⁺ ions at different stages of the first step of the tyrosylation reaction, including the coordination of substrates and involvement in the PP_i releasing. The results of this work suggest that differences between ATP-binding sites of mammalian and bacterial TyrRSs are meaningful and could be exploited in the drug design.     

Molecular dynamics simulations reveal how the reticulon-homology domain of the autophagy receptor RETREG1/FAM134B remodels membranes for efficient selective reticulophagy

ABSTRACT

The autophagy receptor for selective reticulophagy, RETREG1/FAM134B is essential for ER maintenance, and its dysfunction is associated with neuronal disorders, vascular dementia, or viral infections. The protein consists of the reticulon-homology domain (RHD) that is flanked at the N- and C-termini by an intrinsically disordered protein region (IDPR), where the C terminal IDPR carries the indispensable LC3-interacting region (LIR) motif for the interaction with LC3. The RHD of RETREG1 is presumed to play a role in membrane remodeling, but the absence of a known 3D structure of this domain so far prevented researchers from gaining mechanistic insights into how the RETREG1 RHD curves membranes, and thereby facilities reticulophagy. The recent study by Bhaskara et al., which is described in this editor’s corner article, used molecular dynamics (MD) simulations to create a structural model of the RETREG1 RHD. MD simulations along with in vitro liposome remodeling experiments reveal how the RHD domain acts on the ER membrane and, in concert with the C terminal IDPR, executes the function of RETREG1 in selective reticulophagy.

Abbreviations: ER, endoplasmic reticulum; IDPR, intrinsically disordered protein region; LIR, LC3-interacting region; MD, molecular dynamics; RHD, reticulon-homology domain; TM, transmembrane     

Crystal structure of the N‐terminal methyltransferase‐like domain of anamorsin

Anamorsin is a recently identified molecule that inhibits apoptosis during hematopoiesis. It contains an N‐terminal methyltransferase‐like domain and a C‐terminal Fe‐S cluster motif. Not much is known about the function of the protein. To better understand the function of anamorsin, we have solved the crystal structure of the N‐terminal domain at 1.8 Å resolution. Although the overall structure resembles a typical S‐adenosylmethionine (SAM) dependent methyltransferase fold, it lacks one α‐helix and one β‐strand. As a result, the N‐terminal domain as well as the full‐length anamorsin did not show S‐adenosyl‐l ‐methionine (AdoMet) dependent methyltransferase activity. Structural comparisons with known AdoMet dependent methyltransferases reveals subtle differences in the SAM binding pocket that preclude the N‐terminal domain from binding to AdoMet. The N‐terminal methyltransferase‐like domain of anamorsin probably functions as a structural scaffold to inhibit methyl transfers by out‐competing other AdoMet dependant methyltransferases or acts as bait for protein–protein interactions.Proteins 2014; 82:1066–1071. © 2013 Wiley Periodicals, Inc.     

The dimerization of PSGL‐1 is driven by the transmembrane domain via a leucine zipper motif

P‐selectin glycoprotein ligand‐1 (PSGL‐1) is a homodimeric mucin ligand that is important to mediate the earliest adhesive event during an inflammatory response by rapidly forming and dissociating the selectin‐ligand adhesive bonds. Recent research indicates that the noncovalent associations between the PSGL‐1 transmembrane domains (TMDs) can substitute for the C320‐dependent covalent bond to mediate the dimerization of PSGL‐1. In this article, we combined TOXCAT assays and molecular dynamics (MD) simulations to probe the mechanism of PSGL‐1 dimerization. The results of TOXCAT assays and Martini coarse‐grained molecular dynamics (CG MD) simulations demonstrated that PSGL‐1 TMDs strongly dimerized in a natural membrane and a leucine zipper motif was responsible for the noncovalent dimerization of PSGL‐1 TMD since mutations of the residues that occupied a or d positions in an (abcdefg)n leucine heptad repeat motif significantly reduced the dimer activity. Furthermore, we studied the effects of the disulfide bond on the PSGL‐1 dimer using MD simulations. The disulfide bond was critical to form the leucine zipper structure, by which the disulfide bond further improved the stability of the PSGL‐1 dimer. These findings provide insights to understand the transmembrane association of PSGL‐1 that is an important structural basis for PSGL‐1 preferentially binding to P‐selectin to achieve its biochemical and biophysical functions.     

Comparative structural dynamics of Tyrosyl-tRNA synthetase complexed with different substrates explored by molecular dynamics

Several molecular dynamics simulations of S. aureus Tyrosyl-tRNA synthetase (TyrRS) in its free form and complexed with Tyr, ATP, tyrosyl adenylate and inhibitor respectively have been carried out to investigate the ligand-linked conformational stability changes associated with its catalytic cycle. The results show that unliganded S. aureus TyrRS samples a more relaxed conformational space than substrate-bound TyrRS. There are three high flexibility regions encompassing residues 114–118, 128–133, and 226–238 respectively. The region which includes the KMSKS motif (KFGKS in S. aureus TyrRS) shows the highest difference in fluctuations. Hydrogen bond network formed by Tyr, ATP, tyrosyl adenylate and inhibitor with S. aureus TyrRS is discussed. Our simulations suggest the induced-fit conformational changes of the KMSKS loop as follows: the KMSKS loop of substrate-free S. aureus TyrRS adopts an open conformation. The tyrosine binds in the pocket with the KMSKS loop balancing between semi-open and open forms. The ATP binding induces the KMSKS loop to the open form. After the Tyr-AMP is formed, the first two residues of KMSKS loop twists in an anticlockwise direction and drives the loop in a conformation similar to the closed one, while those of the last three GKS residues adopt a conformation between semi-open and open conformation. This conformational change may probably be involved in the initial tRNA binding. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

In Silico Analysis Reveals Sequential Interactions and Protein Conformational Changes during the Binding of Chemokine CXCL-8 to Its Receptor CXCR1

Chemokine CXCL-8 plays a central role in human immune response by binding to and activate its cognate receptor CXCR1, a member of the G-protein coupled receptor (GPCR) family. The full-length structure of CXCR1 is modeled by combining the structures of previous NMR experiments with those from homology modeling. Molecular docking is performed to search favorable binding sites of monomeric and dimeric CXCL-8 with CXCR1 and a mutated form of it. The receptor-ligand complex is embedded into a lipid bilayer and used in multi ns molecular dynamics (MD) simulations. A multi-steps binding mode is proposed: (i) the N-loop of CXCL-8 initially binds to the N-terminal domain of receptor CXCR1 driven predominantly by electrostatic interactions; (ii) hydrophobic interactions allow the N-terminal Glu-Leu-Arg (ELR) motif of CXCL-8 to move closer to the extracellular loops of CXCR1; (iii) electrostatic interactions finally dominate the interaction between the N-terminal ELR motif of CXCL-8 and the EC-loops of CXCR1. Mutation of CXCR1 abrogates this mode of binding. The detailed binding process may help to facilitate the discovery of agonists and antagonists for rational drug design.     

Structure of dystrophia myotonica protein kinase

Dystrophia myotonica protein kinase (DMPK) is a serine/threonine kinase composed of a kinase domain and a coiled‐coil domain involved in the multimerization. The crystal structure of the kinase domain of DMPK bound to the inhibitor bisindolylmaleimide VIII (BIM‐8) revealed a dimeric enzyme associated by a conserved dimerization domain. The affinity of dimerisation suggested that the kinase domain alone is insufficient for dimerisation in vivo and that the coiled‐coil domains are required for stable dimer formation. The kinase domain is in an active conformation, with a fully‐ordered and correctly positioned αC helix, and catalytic residues in a conformation competent for catalysis. The conserved hydrophobic motif at the C‐terminal extension of the kinase domain is bound to the N‐terminal lobe of the kinase domain, despite being unphosphorylated. Differences in the arrangement of the C‐terminal extension compared to the closely related Rho‐associated kinases include an altered PXXP motif, a different conformation and binding arrangement for the turn motif, and a different location for the conserved NFD motif. The BIM‐8 inhibitor occupies the ATP site and has similar binding mode as observed in PDK1.     

NMR structural analysis of Sleeping Beauty transposase binding to DNA

The Sleeping Beauty (SB) transposon is the most widely used DNA transposon in genetic applications and is the only DNA transposon thus far in clinical trials for human gene therapy. In the absence of atomic level structural information, the development of SB transposon relied primarily on the biochemical and genetic homology data. While these studies were successful and have yielded hyperactive transposases, structural information is needed to gain a mechanistic understanding of transposase activity and guides to further improvement. We have initiated a structural study of SB transposase using Nuclear Magnetic Resonance (NMR) and Circular Dichroism (CD) spectroscopy to investigate the properties of the DNA‐binding domain of SB transposase in solution. We show that at physiologic salt concentrations, the SB DNA‐binding domain remains mostly unstructured but its N‐terminal PAI subdomain forms a compact, three‐helical structure with a helix‐turn‐helix motif at higher concentrations of NaCl. Furthermore, we show that the full‐length SB DNA‐binding domain associates differently with inner and outer binding sites of the transposon DNA. We also show that the PAI subdomain of SB DNA‐binding domain has a dominant role in transposase's attachment to the inverted terminal repeats of the transposon DNA. Overall, our data validate several earlier predictions and provide new insights on how SB transposase recognizes transposon DNA.     

Dissecting the multifunctional role of the N‐terminal domain of the Melon necrotic spot virus coat protein in RNA packaging,viral movement and interference with antiviral plant defence

The coat protein (CP) of Melon necrotic spot virus (MNSV) is structurally composed of three major domains. The middle S‐domain builds a robust protein shell around the viral genome, whereas the C‐terminal protruding domain, or P‐domain, is involved in the attachment of virions to the transmission vector. Here, we have shown that the N‐terminal domain, or R‐domain, and the arm region, which connects the R‐domain and S‐domain, are involved in different key steps of the viral cycle, such as cell‐to‐cell movement and the suppression of RNA silencing and pathogenesis through their RNA‐binding capabilities. Deletion mutants revealed that the CP RNA‐binding ability was abolished only after complete, but not partial, deletion of the R‐domain and the arm region. However, a comparison of the apparent dissociation constants for the CP RNA‐binding reaction of several partial deletion mutants showed that the arm region played a more relevant role than the R‐domain in in vitro RNA binding. Similar results were obtained in in vivo assays, although, in this case, full‐length CPs were required to encapsidate full‐length genomes. We also found that the R‐domain carboxyl portion and the arm region were essential for efficient cell‐to‐cell movement, for enhancement of Potato virus X pathogenicity, for suppression of systemic RNA silencing and for binding of small RNAs. Therefore, unlike other carmovirus CPs, the R‐domain and the arm region of MNSV CP have acquired, in addition to other essential functions such as genome binding and encapsidation functions, the ability to suppress RNA silencing by preventing systemic small RNA transport.     

Conformation dynamics of cyclic disulfides and selenosulfides in CXXC(U) (X = Gly,Ala) tetrapeptide redox motifs

Thioredoxin fold proteins often contain a Cys‐(Xxx)_n‐Cys(Sec) or CX_nC(U) motif, where the active cysteine (C) or selenocysteine (U) is bridged by X residues, which vary with protein function. The effect of the X residues on the conformation space of the oxidized disulfide and selenosulfide forms of the CXXC(U) motif has been investigated using molecular dynamics (MD) and density functional theory. Multi‐microsecond‐length MD simulations of the CGGC, CGAC, and CAGC cyclic peptides show that CGGC rings readily exchange between several conformations over the course of the simulation, but steric interactions with the methyl group of Ala limit the conformation space available to the cyclic peptide, especially for CGAC. The potential for the motif to be reduced, as measured by the energy of the lowest unoccupied molecular orbitals, is dependent upon the ring conformation. These results suggest that control of available conformations by the bridging residues and the protein tertiary structure may be important for defining the function of the CXXC motif. Theoretical ⁷⁷Se chemical shifts of the selenosulfide moiety are dependent upon the conformation and/or intramolecular Se···O interactions with the backbone carbonyl group of the C‐terminal U residue.     

Nuclear magnetic resonance‐based determination of dioxygen binding sites in protein cavities

The location and ligand accessibility of internal cavities in cysteine‐free wild‐type T4 lysozyme was investigated using O₂ gas‐pressure NMR spectroscopy and molecular dynamics (MD) simulation. Upon increasing the concentration of dissolved O₂ in solvent to 8.9 mM, O₂‐induced paramagnetic relaxation enhancements (PREs) to the backbone amide and side chain methyl protons were observed, specifically around two cavities in the C‐terminal domain. To determine the number of O₂ binding sites and their atomic coordinates from the 1/r⁶ distance dependence of the PREs, we established an analytical procedure using Akaike's Information Criterion, in combination with a grid‐search. Two O₂‐accessible sites were identified in internal cavities: One site was consistent with the xenon‐binding site in the protein in crystal, and the other site was established to be a novel ligand‐binding site. MD simulations performed at 10 and 100 mM O₂ revealed dioxygen ingress and egress as well as rotational and translational motions of O₂ in the cavities. It is therefore suggested that conformational fluctuations within the ground‐state ensemble transiently develop channels for O₂ association with the internal protein cavities.     

Human MICAL1: Activation by the small GTPase Rab8 and small‐angle X‐ray scattering studies on the oligomerization state of MICAL1 and its complex with Rab8

Human MICAL1 is a member of a recently discovered family of multidomain proteins that couple a FAD‐containing monooxygenase‐like domain to typical protein interaction domains. Growing evidence implicates the NADPH oxidase reaction catalyzed by the flavoprotein domain in generation of hydrogen peroxide as a second messenger in an increasing number of cell types and as a specific modulator of actin filaments stability. Several proteins of the Rab families of small GTPases are emerging as regulators of MICAL activity by binding to its C‐terminal helical domain presumably shifting the equilibrium from the free – auto‐inhibited – conformation to the active one. We here extend the characterization of the MICAL1–Rab8 interaction and show that indeed Rab8, in the active GTP‐bound state, stabilizes the active MICAL1 conformation causing a specific four‐fold increase of k_cat of the NADPH oxidase reaction. Kinetic data and small‐angle X‐ray scattering (SAXS) measurements support the formation of a 1:1 complex between full‐length MICAL1 and Rab8 with an apparent dissociation constant of approximately 8 μM. This finding supports the hypothesis that Rab8 is a physiological regulator of MICAL1 activity and shows how the protein region preceding the C‐terminal Rab‐binding domain may mask one of the Rab‐binding sites detected with the isolated C‐terminal fragment. SAXS‐based modeling allowed us to propose the first model of the free full‐length MICAL1, which is consistent with an auto‐inhibited conformation in which the C‐terminal region prevents catalysis by interfering with the conformational changes that are predicted to occur during the catalytic cycle.     

Functional domain analysis of the Mycoplasma pneumoniae co‐chaperone TopJ

Colonization of conducting airways of humans by the prokaryote Mycoplasma pneumoniae is mediated by a differentiated terminal organelle important in cytadherence, gliding motility and cell division. TopJ is a predicted J‐domain co‐chaperone also having domains unique to mycoplasma terminal organelle proteins and is essential for terminal organelle function, as well as stabilization of protein P24, which is required for normal initiation of terminal organelle formation. J‐domains activate the ATPase of DnaK chaperones, facilitating peptide binding and proper protein folding. We performed mutational analysis of the predicted J‐domain, central acidic and proline‐rich (APR) domain, and C‐terminal domain of TopJ and assessed the phenotypic consequences when introduced into an M. pneumoniae topJ mutant. A TopJ derivative with amino acid substitutions in the canonical J‐domain histidine–proline–aspartic acid motif restored P24 levels but not normal motility, morphology or cytadherence, consistent with a J‐domain co‐chaperone function. In contrast, TopJ derivatives having APR or C‐terminal domain deletions were less stable and failed to restore P24, but resulted in normal morphology, intermediate gliding motility and cytadherence levels exceeding that of wild‐type cells. Results from immunofluorescence microscopy suggest that both the APR and C‐terminal domains, but not the histidine–proline–aspartic acid motif, are critical for TopJ localization to the terminal organelle.     

The unstructured domain of colicin N kills Escherichia coli

Bacteria often produce toxins which kill competing bacteria. Colicins, produced by and toxic to Escherichia coli bacteria are three‐domain proteins so efficient that one molecule can kill a cell. The C‐terminal domain carries the lethal activity and the central domain is required for surface receptor binding. The N‐terminal domain, required for translocation across the outer membrane, is always intrinsically unstructured. It has always been assumed therefore that the C‐terminal cytotoxic domain is required for the bactericidal activity. Here we report the unexpected finding that in isolation, the 90‐residue unstructured N‐terminal domain of colicin N is cytotoxic. Furthermore it causes ion leakage from cells but, unlike known antimicrobial peptides (AMPs) with this property, shows no membrane binding behaviour. Finally, its activity remains strictly dependent upon the same receptor proteins (OmpF and TolA) used by full‐length colicin N. This mechanism of rapid membrane disruption, via receptor mediated binding of a soluble peptide, may reveal a new target for the development of highly specific antibacterials.     

Centromere localization and function of Mis18 requires Yippee‐like domain‐mediated oligomerization

Mis18 is a key regulator responsible for the centromere localization of the CENP‐A chaperone Scm3 in Schizosaccharomyces pombe and HJURP in humans, which establishes CENP‐A chromatin that defines centromeres. The molecular and structural determinants of Mis18 centromere targeting remain elusive. Here, by combining structural, biochemical, and yeast genetic studies, we show that the oligomerization of S. pombe Mis18, mediated via its conserved N‐terminal Yippee‐like domain, is crucial for its centromere localization and function. The crystal structure of the N‐terminal Yippee‐like domain reveals a fold containing a cradle‐shaped pocket that is implicated in protein/nucleic acid binding, which we show is required for Mis18 function. While the N‐terminal Yippee‐like domain forms a homodimer in vitro and in vivo, full‐length Mis18, including the C‐terminal α‐helical domain, forms a homotetramer in vitro. We also show that the Yippee‐like domains of human Mis18α/Mis18β interact to form a heterodimer, implying a conserved structural theme for Mis18 regulation.     

Highly differentiated motifs responsible for two cytokine activities of a split human tRNA synthetase.

While native human tyrosyl-tRNA synthetase (TyrRS) is inactive as a cell-signaling molecule, it can be split into two distinct cytokines. The enzyme is secreted under apoptotic conditions in culture where it is cleaved into an N-terminal fragment that harbors the catalytic site and into a C-domain fragment found only in the mammalian enzymes. The N-terminal fragment is an interleukin-8 (IL-8)-like cytokine, whereas the released C-domain is an endothelial-monocyte-activating polypeptide II (EMAP II)-like cytokine. Although the IL-8-like activity of the N-fragment depends on an ELR motif found in alpha-chemokines and conserved among mammalian TyrRSs, here we show that a similar (NYR) motif in the context of a lower eukaryote TyrRS does not confer the IL8-like activity. We also show that a heptapeptide from the C-domain has EMAP II-like chemotaxis activity for mononuclear phagocytes and polymorphonuclear leukocytes. Eukaryote proteins other than human TyrRS that have EMAP II-like domains have variants of the heptapeptide motif. Peptides based on these sequences are inactive as cytokines. Thus, the cytokine activities of split human TyrRS depend on highly differentiated motifs that are idiosyncratic to the mammalian system.     

Temperature‐accelerated molecular dynamics gives insights into globular conformations sampled in the free state of the AC catalytic domain

The catalytic domain of the adenyl cyclase (AC) toxin from Bordetella pertussis is activated by interaction with calmodulin (CaM), resulting in cAMP overproduction in the infected cell. In the X‐ray crystallographic structure of the complex between AC and the C terminal lobe of CaM, the toxin displays a markedly elongated shape. As for the structure of the isolated protein, experimental results support the hypothesis that more globular conformations are sampled, but information at atomic resolution is still lacking. Here, we use temperature‐accelerated molecular dynamics (TAMD) simulations to generate putative all‐atom models of globular conformations sampled by CaM‐free AC. As collective variables, we use centers of mass coordinates of groups of residues selected from the analysis of standard molecular dynamics (MD) simulations. Results show that TAMD allows extended conformational sampling and generates AC conformations that are more globular than in the complexed state. These structures are then refined via energy minimization and further unrestrained MD simulations to optimize inter‐domain packing interactions, thus resulting in the identification of a set of hydrogen bonds present in the globular conformations. Proteins 2014; 82:2483–2496. © 2014 Wiley Periodicals, Inc.     

Iron inhibits Escherichia coli topoisomerase I activity by targeting the first two zinc‐binding sites in the C‐terminal domain

Escherichia coli DNA topoisomerase I (TopA) contains a 67 kDa N‐terminal catalytic domain and a 30 kDa C‐terminal zinc‐binding region (ZD domain) which has three adjacent tetra‐cysteine zinc‐binding motifs. Previous studies have shown that E. coli TopA can bind both iron and zinc, and that iron binding in TopA results in failure to unwind the negatively supercoiled DNA. Here, we report that each E. coli TopA monomer binds one atom of iron via the first two zinc‐binding motifs in ZD domain and both the first and second zinc‐binding motifs are required for iron binding in TopA. The site‐directed mutagenesis studies further reveal that while the mutation of the third zinc‐binding motif has very little effect on TopA's activity, mutation of the first two zinc‐binding motifs in TopA greatly diminishes the topoisomerase activity in vitro and in vivo, indicating that the first two zinc‐binding motifs in TopA are crucial for its function. The DNA‐binding activity assay and intrinsic tryptophan fluorescence measurements show that iron binding in TopA may decrease the single‐stranded (ss) DNA‐binding activity of ZD domain and also change the protein structure of TopA, which subsequently modulate topoisomerase activity.