Indentification and characterization of a novel CA^(2+)-binding protein in avian erythrocytes

A novel Ca^(2+)-binding protein with Mr of 23 K (designated p23) has been identified in avian erythrocytes and thrombocytes. p23 localizes to the marginal bands (MBs), centrosomes and discrete sites around the nuclear membrane in mature avian erythrocytes. p23 appears to bind Ca^(2+) directly and its interaction with subcellular organelles seems to be modulated by intracellular [Ca^(2+)]. However, its unique protein sequence lacks any known Ca^(2+)-binding motif. Developmental analysis reveals that p23 association to its target structures occurs only at very late stages of bone marrow definitive erythropoeisis. In primitive erythroid cells, p23 distributes diffusely in the cytoplasm and lacks any distinct localization. It is postulated that p23 association to subcellular structures may be induced in part by decreased intracellular [Ca^(2+)]. In vitro and in vivo experiments indicate that p23 does not appear to act as a classical microtubule-associated protein (MAP) but p23 homologues appear to be expressed in MB-containing cells of a variety of species from different vertebrate classes. It has been hypothesized that p23 may play a regulatory role in MB stabilization in a Ca^(2+)-dependent manner.

Binucleated (bnbn) turkey erythrocytes were found to express a truncated p23 variant (designated p21) with identical subcellular localization as p23 except immunostaining reveals the presence of multi-centrosomes in bnbn cells. The p21 sequence has a 62 amino acid deletion at the C-terminus and must therefore have an additional ~40 amino acids at the N-terminus. In addition, p21 seems to have lost the ability to bind Ca^(2+) and its supramolecular interactions are not modulated by intracellular [Ca^(2+)]. These apparent differences between p23 and p21 raised the possibility that the p23/p21 allelism could be the Bn/bn genotype. However, genetic analysis suggested that p23/p21 allelism had no absolute correlation with the Bn/bn genotype.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Investigating the Role of the GET3-GET4/GET5 Interaction During Tail-Anchor Protein Targeting

The proper targeting of membrane proteins is essential to the viability of all cells. Tail-anchored (TA) proteins, defined as having a single transmembrane helix at their C-terminus, are post-translationally targeted to the endoplasmic reticulum (ER) membrane by the GET pathway (Guided Entry of TA proteins). In the yeast pathway, the handover of TA substrates is mediated by the heterotetrameric Get4/Get5 (Get4/5) complex, which tethers the co-chaperone Sgt2 to the central targeting factor, the Get3 ATPase. Although binding of Get4/5 to Get3 is critical for efficient TA targeting, the mechanisms by which Get4 regulates Get3 are unknown. To understand the molecular basis of Get4 function, we used a combination of structural biology, biochemistry, and cell biology. Get4/5 binds across the Get3 dimer interface, in an orientation only compatible with a closed Get3, providing insight into the role of nucleotide in complex formation. Additionally, this structure reveals two functionally distinct binding interfaces for anchoring and ATPase regulation, and loss of the regulatory interface leads to strong defects in vitro and in vivo. Additional crystal structures of the Get3-Get4/5 complex give rise to an alternate conformation, which represents an initial binding interaction mediated by electrostatics that facilitates the rate of subsequent inhibited complex formation. This interface is supported by an in-depth kinetic analysis of the Get3-Get4/5 interaction confirming the two-step complex formation. These results allow us to generate a refined model for Get4/5 function in TA targeting.

Insights into mechanism kinematics for protein motion simulation

Background: The high demanding computational requirements necessary to carry out protein motion simulations make it difficult to obtain information related to protein motion. On the one hand, molecular dynamics simulation requires huge computational resources to achieve satisfactory motion simulations. On the other hand, less accurate procedures such as interpolation methods, do not generate realistic morphs from the kinematic point of view. Analyzing a protein's movement is very similar to serial robots; thus, it is possible to treat the protein chain as a serial mechanism composed of rotational degrees of freedom. Recently, based on this hypothesis, new methodologies have arisen, based on mechanism and robot kinematics, to simulate protein motion. Probabilistic roadmap method, which discretizes the protein configurational space against a scoring function, or the kinetostatic compliance method that minimizes the torques that appear in bonds, aim to simulate protein motion with a reduced computational cost. Results: In this paper a new viewpoint for protein motion simulation, based on mechanism kinematics is presented. The paper describes a set of methodologies, combining different techniques such as structure normalization normalization processes, simulation algorithms and secondary structure detection procedures. The combination of all these procedures allows to obtain kinematic morphs of proteins achieving a very good computational cost-error rate, while maintaining the biological meaning of the obtained structures and the kinematic viability of the obtained motion. Conclusions: The procedure presented in this paper, implements different modules to perform the simulation of the conformational change suffered by a protein when exerting its function. The combination of a main simulation procedure assisted by a secondary structure process, and a side chain orientation strategy, allows to obtain a fast and reliable simulations of protein motion.

Evolution of structures and specific recognitions in aminoacyl-tRNA synthetases

The aminoacyl-tRNA synthetases (AARS) are very important during the protein biosynthesis, which can make the gene sequence be accurately translated into the protein sequence by the specific recognition between AARS and tRNA/amino acids. However, the recog

Different protein tyrosine phosphatase superfamilies resulting from different gene reading frames

Protein tyrosine phosphatases (PTPs) are comprised of two superfamilies, the phosphatase I superfamily containing a single low-molecular-weight PTP (lmwPTP) family and the phosphatase II superfamily including both the higher-molecular-weight PTP (hmwPTP) and the dual-specificity phosphatase (DSP) families. The phosphatase I and H superfamilies are often considered to be the result of convergent evolution. The PTP sequence and structure analyses indicate that lmwPTPs, hmwPTPs, and DSPs share similar structures, functions, and a common signature motif, although they have low sequence identities and a different order of active sites in sequence or a circular permutation. The results of this work suggest that lmwPTPs and hmwPTPs/DSPs are remotely related in evolution. The earliest ancestral gene of PTPs could be from a short fragment containing about 90similar to120 nucleotides or 30similar to40 residues; however, a probable full PTP ancestral gene contained one transcript unit with two lmwPTP genes. All three PTP families may have resulted from a common ancestral gene by a series of duplications, fusions, and circular permutations. The circular permutation in PTPs is caused by a reading frame difference, which is similar to that in DNA methyltransferases. Nevertheless, the evolutionary mechanism of circular permutation in PTP genes seems to be more complicated than that in DNA methyltransferase genes. Both mechanisms in PTPs and DNA methyltransferases can be used to explain how some protein families and superfamilies came to be formed by circular permutations during molecular evolution.

Origin and evolution of new exons in the rodent zinc finger protein 39 gene

The origin of new structures and functions is an important process in evolution. In the past decades, we have obtained some preliminary knowledge of the origin and evolution of new genes. However, as the basic unit of genes, the origin and evolution of exons remain unclear. Because young exons retain the footprints of origination, they can be good materials for studying origin and evolution of new exons. In this paper, we report two young exons in a zinc finger protein gene of rodents. Since they are unique sequences in mouse and rat genome and no homologous sequences were found in the orthologous genes of human and pig, the young exons might originate after the divergence of primates and rodents through exonization of intronic sequences. Strong positive selection was detected in the new exons between mouse and rat, suggesting that these exons have undergone significant functional divergence after the separation of the two species. On the other hand, population genetics data of mouse demonstrate that the new exons have been subject to functional constraint, indicating an important function of the new exons in mouse. Functional analyses suggest that these new exons encode a nuclear localization signal peptide, which may mediate new ways of nuclear protein transport. To our knowledge, this is the first example of the origin and evolution of young exons.

The morphology of decorated amyloid fibers is controlled by the conformation and position of the displayed protein.

Self-assembled structures capable of mediating electron transfer are an attractive scientific and technological goal. Therefore, systematic variants of SH3-Cytochrome b(562) fusion proteins were designed to make amyloid fibers displaying heme-b(562) electron transfer complexes. TEM and AFM data show that fiber morphology responds systematically to placement of b(562) within the fusion proteins. UV-vis spectroscopy shows that, for the fusion proteins under test, only half the fiber-borne b(562) binds heme with high affinity. Cofactor binding also improves the AFM imaging properties and changes the fiber morphology through changes in cytochrome conformation. Systematic observations and measurements of fiber geometry suggest that longitudinal registry of subfilaments within the fiber, mediated by the interaction and conformation of the displayed proteins and their interaction with surfaces, gives rise to the observed morphologies, including defects and kinks. Of most interest is the role of small molecule modulation of fiber structure and mechanical stability. A minimum complexity model is proposed to capture and explain the fiber morphology in the light of these results. Understanding the complex interplay between these factors will enable a fiber design that supports longitudinal electron transfer.

Stability and molecular modeling of triplex DNA inhibiting DNA binding protein binding to the core promoter of hepatitis B virus.

Two three-dimensional structure models of the 21nt oligodeoxyribonucleotides, CPI (G3TG-2TGT2G5TG2TGT) and CP3 (TGTG2TGST2GTG2TG3), were constructed by InsightII (MSI) software in IRIS Indigo2 (SGI) workstation using the crystal structure of TAT tripler formation as the template. The initial structures subsequently were minimized by molecular mechanics. The final structures were believed as the dominant conformation. The results showed that the energy of CP1 is lower than that of CP3, and the former is more stable than the latter. Moreover, the results further proved that the 21nt oligodeoxyribo-nucleotide CP1 stably combines with the core promoter (Cp) fragment of hepatitis B virus (HBV) to form a tripler DNA, and CP1 specifically inhibits a specific cellular factor (DNA binding protein) binding to Cp fragment. These results indicated that specific repression of gene transcription of HBV DNA might be possible by tripler-formation DNA.

Local Frustration Determines Molecular and Macroscopic Helix Structures

The effect of displaying cytochromes from an amyloid fibre is modelled as perturbation of -strands in a bilayer of helical -sheets, thereby explaining the spiral morphology of decorated amyloid and the dynamic response of morphology to cytochrome conformation. The morphology of the modelled fibre, which consists of minimal energy assemblies of rigid building blocks containing two anisotropic interacting units, depends primarily on the rigid constraints between units rather than the soft interactions between them. The framework is a discrete version of the bilayered frustration principle that drives morphology in Bauhinia seedpods. We show that self-assembly of frustrated long range structures can occur if the building blocks themselves are internally frustrated, e.g. amyloid morphology is governed by the conformation of the misfolded protein nucleating the fibre. Our model supports the idea that any peptide sequence can form amyloid if bilayers can form first, albeit stabilised by additional material such as chaperones or cytochromes. Analysis of experimentally derived amyloid structures supports our conclusions and suggests a range of frustration effects, which natural amyloid fibres may exploit. From this viewpoint, amyloid appears as a molecular example of a more general universal bilayered frustration principle, which may have profound implications for materials design using fibrous systems. Our model provides quantitative guidance for such applications. The relevance to longer length scales was proved by designing the morphology of a series of macroscopic magnetic stacks. Finally, this work leads to the idea of mixing controlled morphologically defined species to generate higher-order assembly and complex functional behaviour. The systematic kinking of decorated fibres and the nested frustration of the Bauhinia seed pod are two outstanding examples.

Synthesis, characterization, and self-assembly of protein lysozyme monolayer-stabilized gold nanoparticles

Lysozyme monolayer-protected gold nanoparticles (Au NPs) which are hydrophilic and biocompatible and show excellent colloidal stability at low temperature, ca. 4 degrees C, were synthesized in aqueous medium by chemical reduction of HAuCl4 with NaBH4 in the presence of a familiar small enzyme, lysozyme. UV-vis spectra, transmission electron microscopy (TEM), atomic force microscopy, and X-ray photoelectron spectroscopy characterization of the as prepared nanoparticles revealed the formation of well-dispersed An NPs of ca. 2 nm diameter. Moreover, the color change of the An NP solution as well as UV-vis spectroscopy and TEM measurements have also demonstrated the occurrence of Ostwald ripening of the nanoparticles at low temperature. Further characterization with Fourier transform infrared spectroscopy (FTIR) and dynamic light scattering indicated the formation of a monolayer of lysozyme molecules on the particle surface. FTIR data also indicated the intactness of the protein molecules coated on An NPs. All the characterization results showed that the monodisperse An NPs are well-coated directly with lysozyme. Driven by the dipole-dipole attraction, the protein-stabilized Au NPs self-assembled into network structures and nanowires upon aging under ambient temperature.

Single-molecule dynamics reveals cooperative binding-folding in protein recognition

The study of associations between two biomolecules is the key to understanding molecular function and recognition. Molecular function is often thought to be determined by underlying structures. Here, combining a single-molecule study of protein binding with an energy-landscape-inspired microscopic model, we found strong evidence that biomolecular recognition is determined by flexibilities in addition to structures. Our model is based on coarse-grained molecular dynamics on the residue level with the energy function biased toward the native binding structure ( the Go model). With our model, the underlying free-energy landscape of the binding can be explored. There are two distinct conformational states at the free-energy minimum, one with partial folding of CBD itself and significant interface binding of CBD to Cdc42, and the other with native folding of CBD itself and native interface binding of CBD to Cdc42. This shows that the binding process proceeds with a significant interface binding of CBD with Cdc42 first, without a complete folding of CBD itself, and that binding and folding are then coupled to reach the native binding state.

Protein structural characterization by scanning tunneling microscopy with electrochemistry

The structure characterization of proteins or enzymes by STM on electrochemically prepared HOPG surface studied in this laboratory is reviewed. The serial structures of Hb were observed. The differences between the denaturation and inactivation of HRP were investigated by in situ and ex situ STM. The structural variation of Hb in an organic solvent was imaged while protein denaturation was easily observed in a polar solvent.

Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.

Genome-wide survey of putative serine/threonine protein kinases in cyanobacteria

Background: Serine/threonine kinases (STKs) have been found in an increasing number of prokaryotes, showing important roles in signal transduction that supplement the well known role of two-component system. Cyanobacteria are photoautotrophic prokaryotes able to grow in a wide range of ecological environments, and their signal transduction systems are important in adaptation to the environment. Sequence information from several cyanobacterial genomes offers a unique opportunity to conduct a comprehensive comparative analysis of this kinase family. In this study, we extracted information regarding Ser/Thr kinases from 21 species of sequenced cyanobacteria and investigated their diversity, conservation, domain structure, and evolution. Results: 286 putative STK homologues were identified. STKs are absent in four Prochlorococcus strains and one marine Synechococcus strain and abundant in filamentous nitrogen-fixing cyanobacteria. Motifs and invariant amino acids typical in eukaryotic STKs were conserved well in these proteins, and six more cyanobacteria- or bacteria-specific conserved residues were found. These STK proteins were classified into three major families according to their domain structures. Fourteen types and a total of 131 additional domains were identified, some of which are reported to participate in the recognition of signals or substrates. Cyanobacterial STKs show rather complicated phylogenetic relationships that correspond poorly with phylogenies based on 16S rRNA and those based on additional domains. Conclusion: The number of STK genes in different cyanobacteria is the result of the genome size, ecophysiology, and physiological properties of the organism. Similar conserved motifs and amino acids indicate that cyanobacterial STKs make use of a similar catalytic mechanism as eukaryotic STKs. Gene gain-and-loss is significant during STK evolution, along with domain shuffling and insertion. This study has established an overall framework of sequence-structure-function interactions for the STK gene family, which may facilitate further studies of the role of STKs in various organisms.