Three distinct molecular surfaces in ephrin-A5 are essential for a functional interaction with EphA3

Eph receptor tyrosine kinases (Ephs) function as molecular relays that interact with cell surface-bound ephrin ligands to direct the position of migrating cells. Structural studies revealed that, through two distinct contact surfaces on opposite sites of each protein, Eph and ephrin binding domains assemble into symmetric, circular heterotetramers. However, Eph signal initiation requires the assembly of higher order oligomers, suggesting additional points of contact. By screening a random library of EphA3 binding-compromised ephrin-A5 mutants, we have now determined ephrin-A5 residues that are essential for the assembly of high affinity EphA3 signaling complexes. In addition to the two interfaces predicted from the crystal structure of the homologous EphB2 center dot ephrin-B2 complex, we identified a cluster of 10 residues on the ephrin-A5 E alpha-helix, the E-F loop, the underlying H beta-strand, as well as the nearby B - C loop, which define a distinct third surface required for oligomerization and activation of EphA3 signaling. Together with a corresponding third surface region identified recently outside of the minimal ephrin binding domain of EphA3, our findings provide experimental evidence for the essential contribution of three distinct protein-interaction interfaces to assemble functional EphA3 signaling complexes.

Molecular characterization of Osh6p, an oxysterol binding protein homolog in the yeast Saccharomyces cerevisiae

Oxysterol binding protein (OSBP) and its homologs have been shown to regulate lipid metabolism and vesicular transport. However, the exact molecular function of individual OSBP homologs remains uncharacterized. Here we demonstrate that the yeast OSBP homolog, Osh6p, bound phosphatidic acid and phosphoinositides via its N-terminal half containing the conserved OSBP-related domain (ORD). Using a green fluorescent protein fusion chimera, Osh6p was found to localize to the cytosol and patch-like or punctate structures in the vicinity of the plasma membrane. Further examination by domain mapping demonstrated that the N-terminal half was associated with FM4-64 positive membrane compartments; however, the C-terminal half containing a putative coiled-coil was localized to the nucleoplasm. Functional analysis showed that the deletion of OSH6 led to a significant increase in total cellular ergosterols, whereas OSH6 overexpression caused both a significant decrease in ergosterol levels and resistance to nystatin. Oleate incorporation into sterol esters was affected in OSH6 overexpressing cells. However, Lucifer yellow internalization, and FM4-64 uptake and transport were unaffected in both OSH6 deletion and overexpressing cells. Furthermore, osh6 Delta exhibited no defect in carboxypeptidase Y transport and maturation. Lastly, we demonstrated that both the conserved ORD and the putative coiled-coil motif were indispensable for the in vivo function of Osh6p. These data suggest that Osh6p plays a role primarily in regulating cellular sterol metabolism, possibly stero transport.

Molecular basis of intramolecular electron transfer in sulfite-oxidizing enzymes is revealed by high resolution structure of a heterodimeric complex of the catalytic molybdopterin subunit and a c-type cytochrome subunit

Sulfite-oxidizing molybdoenzymes convert the highly reactive and therefore toxic sulfite to sulfate and have been identified in insects, animals, plants, and bacteria. Although the well studied enzymes from higher animals serve to detoxify sulfite that arises from the catabolism of sulfur-containing amino acids, the bacterial enzymes have a central role in converting sulfite formed during dissimilatory oxidation of reduced sulfur compounds. Here we describe the structure of the Starkeya novella sulfite dehydrogenase, a heterodimeric complex of the catalytic molybdopterin subunit and a c-type cytochrome subunit, that reveals the molecular mechanism of intramolecular electron transfer in sulfite-oxidizing enzymes. The close approach of the two redox centers in the protein complex (Mo-Fe distance 16.6 angstrom) allows for rapid electron transfer via tunnelling or aided by the protein environment. The high resolution structure of the complex has allowed the identification of potential through-bond pathways for electron transfer including a direct link via Arg-55A and/or an aromatic-mediated pathway. A potential site of electron transfer to an external acceptor cytochrome c was also identified on the SorB subunit on the opposite side to the interaction with the catalytic SorA subunit.

Sociality and the rate of molecular evolution

The molecular clock does not tick at a uniform rate in all taxa but maybe influenced by species characteristics. Eusocial species (those with reproductive division of labor) have been predicted to have faster rates of molecular evolution than their nonsocial relatives because of greatly reduced effective population size; if most individuals in a population are nonreproductive and only one or few queens produce all the offspring, then eusocial animals could have much lower effective population sizes than their solitary relatives, which should increase the rate of substitution of nearly neutral mutations. An earlier study reported faster rates in eusocial honeybees and vespid wasps but failed to correct for phylogenetic nonindependence or to distinguish between potential causes of rate variation. Because sociality has evolved independently in many different lineages, it is possible to conduct a more wide-ranging study to test the generality of the relationship. We have conducted a comparative analysis of 25 phylogenetically independent pairs of social lineages and their nonsocial relatives, including bees, wasps, ants, termites, shrimps, and mole rats, using a range of available DNA sequences (mitochondrial and nuclear DNA coding for proteins and RNAs, and nontranslated sequences). By including a wide range of social taxa, we were able to test whether there is a general influence of sociality on rates of molecular evolution and to test specific predictions of the hypothesis: (1) that social species have faster rates because they have reduced effective population sizes; (2) that mitochondrial genes would show a greater effect of sociality than nuclear genes; and (3) that rates of molecular evolution should be correlated with the degree of sociality. We find no consistent pattern in rates of molecular evolution between social and nonsocial lineages and no evidence that mitochondrial genes show faster rates in social taxa. However, we show that the most highly eusocial Hymenoptera do have faster rates than their nonsocial relatives. We also find that social parasites (that utilize the workers from related species to produce their own offspring) have faster rates than their social relatives, which is consistent with an effect of lower effective population size on rate of molecular evolution. Our results illustrate the importance of allowing for phylogenetic nonindependence when conducting investigations of determinants of variation in rate of molecular evolution.

hnRNP A2, a potential ssDNA/RNA molecular adapter at the telomere

The heterogeneous nuclear ribonucleoprotein (hnRNP) A2 is a multi-tasking protein that acts in the cytoplasm and nucleus. We have explored the possibility that this protein is associated with telomeres and participates in their maintenance. Rat brain hnRNP A2 was shown to have two nucleic acid binding sites. In the presence of heparin one site binds single-stranded oligodeoxyribonucleotides irrespective of sequence but not the corresponding oligoribonucleotides. Both the hnRNP A2-binding cis-acting element for the cytoplasmic RNA trafficking element, A2RE, and the ssDNA telomere repeat match a consensus sequence for binding to a second sequence-specific site identified by mutational analysis. hnRNP A2 protected the telomeric repeat sequence, but not the complementary sequence, against DNase digestion: the glycine-rich domain was found to be necessary, but not sufficient, for protection. The N-terminal RRM (RNA recognition motif) and tandem RRMs of hnRNP A2 also bind the single-stranded, template-containing segment of telomerase RNA. hnRNP A2 colocalizes with telomeric chromatin in the subset of PML bodies that are a hallmark of ALT cells, reinforcing the evidence for hnRNPs having a role in telomere maintenance. Our results support a model in which hnRNP A2 acts as a molecular adapter between single-stranded telomeric repeats, or telomerase RNA, and another segment of ssDNA.

CandiVF - Candida albicans virulence factor database

Candida albicans is a pathogen commonly infecting patients who receive immunosuppressive drug therapy, long-term catheterization, or those who suffer from acquired immune deficiency syndrome (AIDS). The major factor accountable for pathogenicity of C. albicans is host immune status. Various virulence molecules, or factors, of are also responsible for the disease progression. Virulence proteins are published in public databases but they normally lack detailed functional annotations. We have developed CandiVF, a specialized database of C. albicans virulence factors (http://antigen.i2r.a-star.edu.sg/Templar/DB/CandiVF/) to facilitate efficient extraction and analysis of data aimed to assist research on immune responses, pathogenesis, prevention, and control of candidiasis. CandiVF contains a large number of annotated virulence proteins, including secretory, cell wall-associated, membrane, cytoplasmic, and nuclear proteins. This database has in-built bioinformatics tools including keyword and BLAST search, visualization of 3D-structures, HLA-DR epitope prediction, virulence descriptors, and virulence factors ontology.

Functional characterisation of an engineered multidomain human P450 E1 by molecular Lego

The human cytochrome P450s constitute an important family of monooxygenase enzymes that carry out essential roles in the metabolism of endogenous compounds and foreign chemicals. We present here results of a fusion between a human P450 enzyme and a bacterial reductase that for the first time is shown does not require the addition of lipids or detergents to achieve wild-type-like activities. The fusion enzyme, P450 2E1-BMR, contains the N-terminally modified residues 22-493 of the human P450 2E1 fused at the C-terminus to residues 473-1049 of the P450 BM3 reductase (BMR). The P450 2E1-BMR enzyme is active, self-sufficient and presents the typical marker activities of the native human P450 2E1: the hydroxylation of p-nitrophenol (K (M)=1.84 +/- 0.09 mM and k (cat) of 2.98 +/- 0.04 nmol of p-nitrocatechol formed per minute per nanomole of P450) and chlorzoxazone (K (M)=0.65 +/- 0.08 mM and k (cat) of 0.95 +/- 0.10 nmol of 6-hydroxychlorzoxazone formed per minute per nanomole of P450). A 3D model of human P450 2E1 was generated to rationalise the functional data and to allow an analysis of the surface potentials. The distribution of charges on the model of P450 2E1 compared with that of the FMN domain of BMR provides the ground for the understanding of the interaction between the fused domains. The results point the way to successfully engineer a variety of catalytically self-sufficient human P450 enzymes for drug metabolism studies in solution.

Mimicking the action of GroEL in molecular dynamics simulations: Application to the refinement of protein structures

Bacterial chaperonin, GroEL, together with its co-chaperonin, GroES, facilitates the folding of a variety of polypeptides. Experiments suggest that GroEL stimulates protein folding by multiple cycles of binding and release. Misfolded proteins first bind to an exposed hydrophobic surface on GroEL. GroES then encapsulates the substrate and triggers its release into the central cavity of the GroEL/ES complex for folding. In this work, we investigate the possibility to facilitate protein folding in molecular dynamics simulations by mimicking the effects of GroEL/ES namely, repeated binding and release, together with spatial confinement. During the binding stage, the (metastable) partially folded proteins are allowed to attach spontaneously to a hydrophobic surface within the simulation box. This destabilizes the structures, which are then transferred into a spatially confined cavity for folding. The approach has been tested by attempting to refine protein structural models generated using the ROSETTA procedure for ab initio structure prediction. Dramatic improvements in regard to the deviation of protein models from the corresponding experimental structures were observed. The results suggest that the primary effects of the GroEL/ES system can be mimicked in a simple coarse-grained manner and be used to facilitate protein folding in molecular dynamics simulations. Furthermore, the results Sur port the assumption that the spatial confinement in GroEL/ES assists the folding of encapsulated proteins.

Modelling the structure of latexin-carboxypeptidase A complex based on chemical cross-linking and molecular docking

We have determined the three-dimensional structure of the protein complex between latexin and carboxypeptidase A using a combination of chemical cross-linking, mass spectrometry and molecular docking. The locations of three intermolecular cross-links were identified using mass spectrometry and these constraints were used in combination with a speed-optimised docking algorithm allowing us to evaluate more than 3 x 10(11) possible conformations. While cross-links represent only limited structural constraints, the combination of only three experimental cross-links with very basic molecular docking was sufficient to determine the complex structure. The crystal structure of the complex between latexin and carboxypeptidase A4 determined recently allowed us to assess the success of this structure determination approach. Our structure was shown to be within 4 angstrom r.m.s. deviation of C alpha atoms of the crystal structure. The study demonstrates that cross-linking in combination with mass spectrometry can lead to efficient and accurate structural modelling of protein complexes.

Molecular crowding effects of linear polymers in protein solutions

Measurement of protein-polymer second virial coefficients (B-AP) by sedimentation equilibrium studies of carbonic anhydrase and cytochrome c in the presence of dextrans (T10-T80) has revealed an inverse dependence of B-AP upon dextran molecular mass that conforms well with the behaviour predicted for the excluded-volume interaction between a spherical protein solute A and a random-flight representation of the polymeric cosolute P. That model of the protein-polymer interaction is also shown to provide a reasonable description of published gel chromatographic and equilibrium dialysis data on the effect of polymer molecular mass on BAP for human serum albumin in the presence of polyethylene glycols, a contrary finding from analysis of albumin solubility measurements being rejected on theoretical grounds. Inverse dependence upon polymer chainlength is also the predicted excluded-volume effect on the strength of several types of macromolecular equilibria-protein isomerization, protein dimerization, and 1 : 1 complex formation between dissimilar protein reactants. It is therefore concluded that published experimental observations of the reverse dependence, preferential reaction enhancement within DNA replication complexes by larger polyethylene glycols, must reflect the consequences of cosolute chemical interactions that outweigh those of thermodynamic nonideality arising from excluded-volume effects. (c) 2005 Elsevier B.V. All rights reserved.

Molecular determinants of ginkgolide binding in the glycine receptor pore

Ginkgolides are potent blockers of the glycine receptor Cl- channel (GlyR) pore. We sought to identify their binding sites by comparing the effects of ginkgolides A, B and C and bilobalide on alpha 1, alpha 2, alpha 1 beta and alpha 2 beta GlyRs. Bilobalide sensitivity was drastically reduced by incorporation of the beta subunit. In contrast, the sensitivities to ginkgolides B and C were enhanced by beta subunit expression. However, ginkgolide A sensitivity was increased in the alpha 2 beta GlyR relative to the alpha 2 GlyR but not in the alpha 1 beta GlyR relative to the alpha 1 GlyR. We hypothesised that the subunit-specific differences were mediated by residue differences at the second transmembrane domain 2' and 6' pore-lining positions. The increased ginkgolide A sensitivity of the alpha 2 beta GlyR was transferred to the alpha 1 beta GlyR by the G2'A (alpha 1 to alpha 2 subunit) substitution. In addition, the alpha 1 subunit T6'F mutation abolished inhibition by all ginkgolides. As the ginkgolides share closely related structures, their molecular interactions with pore-lining residues were amenable to mutant cycle analysis. This identified an interaction between the variable R2 position of the ginkgolides and the 2' residues of both alpha 1 and beta subunits. These findings provide strong evidence for ginkgolides binding at the 2' pore-lining position.

Mutations in the regulatory subunit of yeast acetohydroxyacid synthase affect its activation by MgATP

Isoleucine, leucine and valine are synthesized via a common pathway in which the first reaction is catalysed by AHAS (acetohydroxyacid synthase; EC 2.2.1.6). This heterotetrameric enzyme is composed of a larger subunit that contains the catalytic machinery and a smaller subunit that plays a regulatory role. The RSU (regulatory subunit) enhances the activity of the CSU (catalytic sub unit) and mediates end-product inhibition by one or more of the branched-chain amino acids, usually valine. Fungal AHAS differs front that in other organisms in that the inhibition by valine is reversed by MgATP. The fungal AHAS RSU also differs from that in other organisms in that it contains a sequence insert. We suggest that this insert may form the MgATP-binding site and we have tested this hypothesis by mutating ten highly conserved amino acid residues of the yeast AHAS RSU. The modified subunits were tested for their ability to activate the yeast AHAS CSU, to confer sensitivity to valine inhibition and to mediate reversal of the inhibition by MgATP. All but one of the mutations resulted in substantial changes in the properties of the RSU. Unexpectedly, four of them gave a protein that required mgATP in order for strong stimulation of the CSU and valine inhibition to be observed. A model to explain this result is proposed. Five of the mutations abolished MgATP activation and are suggested to constitute the binding site for this modulator.

Molecular diversity in venom from the Australian Brown Snake, Pseudonaja textilis

Venom from the Australian elapid Pseudonaja textilis (Common or Eastern Brown snake), is the second most toxic snake venom known and is the most common cause of death from snake bite in Australia. This venom is known to contain a prothrombin activator complex, serine proteinase inhibitors, various phospholipase A(2)s, and pre-and postsynaptic neurotoxins. In this study, we performed a proteomic identification of the venom using two- dimensional gel electrophoresis, mass spectrometry, and de novo peptide sequencing. We identified most of the venom proteins including proteins previously not known to be present in the venom. In addition, we used immunoblotting and post-translational modification-specific enzyme stains and antibodies that reveal the complexity and regional diversity of the venom. Modifications observed include phosphorylation, gamma-carboxylation, and glycosylation. Glycoproteins were further characterized by enzymatic deglycosylation and by lectin binding specificity. The venom contains an abundance of glycoproteins with N-linked sugars that include glucose/mannose, N-acetylgalactosamine, N-acetylglucosamine, and sialic acids. Additionally there are multiple isoforms of mammalian coagulation factors that comprise a significant proportion of the venom. Indeed two of the identified proteins, a procoagulant and a plasmin inhibitor, are currently in development as human therapeutic agents.

Molecular dynamics characterization of n-octyl-beta-D-glucopyranoside micelle structure in aqueous solution

n-Octyl-beta-D-glueopyranoside (OG) is a non-ionic glycolipid, which is used widely in biotechnical and biochemical applications. All-atom molecular dynamics simulations from two different initial coordinates and velocities in explicit solvent have been performed to characterize the structural behaviour of an OG aggregate at equilibrium conditions. Geometric packing properties determined from the simulations and small angle neutron scattering experiment state that OG micelles are more likely to exist in a non-spherical shape, even at the concentration range near to the critical micelle concentration (0.025 M). Despite few large deviations in the principal moment of inertia ratios, the average micelle shape calculated from both simulations is a prolate ellipsoid. The deviations at these time scales are presumably the temporary shape change of a micelle. However, the size of the micelle and the accessible surface areas were constant during the simulations with the micelle surface being rough and partially elongated. Radial distribution functions computed for the hydroxyl oxygen atoms of an OG show sharper peaks at a minimum van der Waals contact distance than the acetal oxygen, ring oxygen, and anomeric carbon atoms. This result indicates that these atoms are pointed outwards at the hydrophilic/hydrophobic interface, form hydrogen bonds with the water molecules, and thus hydrate the micelle surface effectively. (c) 2005 Elsevier Inc. All rights reserved.

Molecular mechanisms for the variation of mitochondrial gene content and gene arrangement among chigger mites of the genus Leptotrombidium (Acari : Acariformes)

The gene content of a mitochondrial (mt) genome, i.e., 37 genes and a large noncoding region (LNR), is usually conserved in Metazoa. The arrangement of these genes and the LNR is generally conserved at low taxonomic levels but varies substantially at high levels. We report here a variation in mt gene content and gene arrangement among chigger mites of the genus Leptotrombidium. We found previously that the mt genome of Leptotrombidium pallidum has an extra gene for large-subunit rRNA (rrnL), a pseudo-gene for small-subunit rRNA (PrrnS), and three extra LNRs, additional to the 37 genes and an LNR typical of Metazoa. Further, the arrangement of mt genes of L. pallidum differs drastically from that of the hypothetical ancestor of the arthropods. To find to what extent the novel gene content and gene arrangement occurred in Leptotrombidium, we sequenced the entire or partial mt genomes of three other species, L. akamushi, L. deliense, and L. fletcheri. These three species share the arrangement of all genes with L. pallidum, except trnQ (for tRNA-glutamine). Unlike L. pallidum, however, these three species do not have extra rrnL or PrrnS and have only one extra LNR. By comparison between Leptotrombidium species and the ancestor of the arthropods, we propose that (1) the type of mt genome present in L. pallidum evolved from the type present in the other three Leptotrombidium species, and (2) three molecular mechanisms were involved in the evolution of mt gene content and gene arrangement in Leptotrombidium species.