Model building by comparison at CASP3: using expert knowledge and computer automation

Rotamer Dynamics: Analysis of Rotamers in Molecular Dynamics Simulations of Proteins

Given by χ torsional angles, rotamers describe the side-chain conformations of amino acid residues in a protein based on the rotational isomers (hence the word rotamer). Constructed rotamer libraries, based on either protein crystal structures or dynamics studies, are the tools for classifying rotamers (torsional angles) in a way that reflect their frequency in nature. Rotamer libraries are routinely used in structure modeling and evaluation. In this perspective article, we would like to encourage researchers to apply rotamer analyses beyond their traditional use. Molecular dynamics (MD) of proteins highlight the in silico behavior of molecules in solution and thus can identify favorable side-chain conformations. In this article, we used simple computational tools to study rotamer dynamics (RD) in MD simulations. First, we isolated each frame in the MD trajectories in separate Protein Data Bank files via the cpptraj module in AMBER. Then, we extracted torsional angles via the Bio3D module in R language. The classification of torsional angles was also done in R according to the penultimate rotamer library. RD analysis is useful for various applications such as protein folding, study of rotamer-rotamer relationship in protein-protein interaction, real-time correlation between secondary structures and rotamers, study of flexibility of side chains in binding site for molecular docking preparations, use of RD as guide in functional analysis and study of structural changes caused by mutations, providing parameters for improving coarse-grained MD accuracy and speed, and many others. Major challenges facing RD to emerge as a new scientific field involve the validation of results via easy, inexpensive wet-lab methods. This realm is yet to be explored.

Theoretical and experimental study of the antifreeze protein AFP752, trehalose and dimethyl sulfoxide cryoprotection mechanism: correlation with cryopreserved cell viability

In this work the physico-chemical properties of selected cryoprotectants (antifreeze protein TrxA-AFP752, trehalose and dimethyl sulfoxide) were correlated with their impact on the constitution of ice and influence on frozen/thawed cell viability. The freezing processes and states of investigated materials solutions were described and explained from a fundamental point of view using ab-initio modelling (molecular dynamics, DFT), Raman spectroscopy, Differential Scanning Calorimetry and X-Ray Diffraction. For the first time, in this work we correlated the microscopic view (modelling) with the description of the frozen solution states and put these results in the context of human skin fibroblast viability after freezing and thawing. DMSO and AFP had different impacts on their solution's freezing process but in both cases the ice crystallinity size was considerably reduced. DMSO and AFP treatment in different ways improved the viability of frozen/thawed cells.

Discovery of catalases in members of the Chlamydiales order

Catalase is an important virulence factor for survival in macrophages and other phagocytic cells. In Chlamydiaceae, no catalase had been described so far. With the sequencing and annotation of the full genomes of Chlamydia-related bacteria, the presence of different catalase-encoding genes has been documented. However, their distribution in the Chlamydiales order and the functionality of these catalases remain unknown. Phylogeny of chlamydial catalases was inferred using MrBayes, maximum likelihood, and maximum parsimony algorithms, allowing the description of three clade 3 and two clade 2 catalases. Only monofunctional catalases were found (no catalase-peroxidase or Mn-catalase). All presented a conserved catalytic domain and tertiary structure. Enzymatic activity of cloned chlamydial catalases was assessed by measuring hydrogen peroxide degradation. The catalases are enzymatically active with different efficiencies. The catalase of Parachlamydia acanthamoebae is the least efficient of all (its catalytic activity was 2 logs lower than that of Pseudomonas aeruginosa). Based on the phylogenetic analysis, we hypothesize that an ancestral class 2 catalase probably was present in the common ancestor of all current Chlamydiales but was retained only in Criblamydia sequanensis and Neochlamydia hartmannellae. The catalases of class 3, present in Estrella lausannensis and Parachlamydia acanthamoebae, probably were acquired by lateral gene transfer from Rhizobiales, whereas for Waddlia chondrophila they likely originated from Legionellales or Actinomycetales. The acquisition of catalases on several occasions in the Chlamydiales suggests the importance of this enzyme for the bacteria in their host environment.

α(2)-µ-Globulin fragment (a2-f) from kidneys of male rats

The structure of α₂-µ-Globulin fragment (A2-f) is not known.α₂-µ-Globulin fragment (A2-f) is a 15.5 kDa protein that binds equimolar amount of fatty acids in male rat kidneys. The expression of this protein has been shown to change in response to druginduced and genetic hypertension which suggests that it plays an important role in renal fatty acid metabolism under pathological conditions as well as normal conditions. A2-f has sequence homology with amino acid 28-178 of α₂-µ-Globulin (A2U) that is synthesized pre-dominantly in the male rat liver and is present in the urine. It is believed that unusual structural features permit A2-f to be targeted to the proximal tubule cell; to escape lysosomal degradation in liver and to enter the cytosol of proximal tubule cells of the kidneys. Homology modeling has been employed to determine the structural elements of this protein and they have been compared with the published structure of A2U. Results suggest differences between the structure of A2-f and its precursor protein A2U.

Defining and predicting structurally conserved regions in protein superfamilies

MOTIVATION

The structures of homologous proteins are generally better conserved than their sequences. This phenomenon is demonstrated by the prevalence of structurally conserved regions (SCRs) even in highly divergent protein families. Defining SCRs requires the comparison of two or more homologous structures and is affected by their availability and divergence, and our ability to deduce structurally equivalent positions among them. In the absence of multiple homologous structures, it is necessary to predict SCRs of a protein using information from only a set of homologous sequences and (if available) a single structure. Accurate SCR predictions can benefit homology modelling and sequence alignment.

RESULTS

Using pairwise DaliLite alignments among a set of homologous structures, we devised a simple measure of structural conservation, termed structural conservation index (SCI). SCI was used to distinguish SCRs from non-SCRs. A database of SCRs was compiled from 386 SCOP superfamilies containing 6489 protein domains. Artificial neural networks were then trained to predict SCRs with various features deduced from a single structure and homologous sequences. Assessment of the predictions via a 5-fold cross-validation method revealed that predictions based on features derived from a single structure perform similarly to ones based on homologous sequences, while combining sequence and structural features was optimal in terms of accuracy (0.755) and Matthews correlation coefficient (0.476). These results suggest that even without information from multiple structures, it is still possible to effectively predict SCRs for a protein. Finally, inspection of the structures with the worst predictions pinpoints difficulties in SCR definitions.

AVAILABILITY

The SCR database and the prediction server can be found at http://prodata.swmed.edu/SCR.

CONTACT

91huangi@gmail.com or grishin@chop.swmed.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Online.

Eukaryotic translation elongation factor 1A (eEF1A) domain I from S. cerevisiae is required but not sufficient for inter-species complementation

Ethanolamine phosphoglycerol (EPG) is a protein modification attached exclusively to eukaryotic elongation factor 1A (eEF1A). In mammals and plants, EPG is linked to conserved glutamate residues located in eEF1A domains II and III, whereas in the unicellular eukaryote Trypanosoma brucei, only domain III is modified by a single EPG. A biosynthetic precursor of EPG and structural requirements for EPG attachment to T. brucei eEF1A have been reported, but nothing is known about the EPG modifying enzyme(s). By expressing human eEF1A in T. brucei, we now show that EPG attachment to eEF1A is evolutionarily conserved between T. brucei and Homo sapiens. In contrast, S. cerevisiae eEF1A, which has been shown to lack EPG is not modified in T. brucei. Furthermore, we show that eEF1A cannot functionally complement across species when using T. brucei and S. cerevisiae as model organisms. However, functional complementation in yeast can be obtained using eEF1A chimera containing domains II or III from other species. In contrast, yeast domain I is strictly required for functional complementation in S. cerevisiae.

Identification and bioinformatics characterization of translation initiation complex eIF4F components and poly(A)-binding protein from Plasmodium falciparum

Protein synthesis in eukaryotes initiates with binding of the multisubunit translation initiation complex eIF4F. This complex contains eIF4E, eIF4A and eIF4G. eIF4E directly interacts with the cap structure, eIF4A is an RNA helicase and eIF4G acts as a scaffold for the complex. eIF4G contains the binding sites for both the subunits i.e., eIF4A and eIF4E and it also interacts with poly(A)-binding protein (PABP). In present study we have identified and characterized the main components of the eIF4F complex i.e., eIF4E, eIF4A and eIF4G and PABP from Plasmodium falciparum. Molecular modeling of PfeIF4E, PfeIF4G and PfPABP confirms that they contain all the characteristic conserved structural features. We have annotated some of the genes of P. falciparum and as a result these studies demonstrate that the components of translation initiation complex are highly conserved. Therefore these studies will contribute to understand the basic biology and components of translation complex in P. falciparum.

In silico analysis for transcription factors with Zn(II)(2)C(6) binuclear cluster DNA-binding domains in Candida albicans

A total of 6047 open reading frames in the Candida albicans genome were screened for Zn(II)₂C₆-type zinc cluster proteins (or binuclear cluster proteins) involved in DNA recognition. These fungal proteins are transcription regulators of genes involved in a wide range of cellular processes, including metabolism of different compounds such as sugars or amino acids, as well as multi-drug resistance, control of meiosis, cell wall architecture, etc. The selection criteria used in the sequence analysis were the presence of the CysX₂CysX₆CysX_5-16CysX₂CysX_6-8Cys motif and a putative nuclear localization signal. Using this approach, 70 putative Zn(II)₂C₆ transcription factors have been found in the genome of C. albicans.

Evolutionary dynamics of the mS952 intron: a novel mitochondrial group II intron encoding a LAGLIDADG homing endonuclease gene

Examination of the mitochondrial small subunit ribosomal RNA (rns) gene of five species of the fungal genus Leptographium revealed that the gene has been invaded at least once at position 952 by a group II intron encoding a LAGLIDADG homing endonuclease gene. Phylogenetic analyses of the intron and homing endonuclease sequences indicated that each element in Leptographium species forms a single clade and is closely related to the group II intron/homing endonuclease gene composite element previously reported at position 952 of the mitochondrial rns gene of Cordyceps species and of Cryphonectria parasitica. The results of an intron survey of the mt rns gene of Leptographium species superimposed onto the phylogenetic analysis of the host organisms suggest that the composite element was transmitted vertically in Leptographium lundbergii. However, its stochastic distribution among strains of L. wingfieldii, L. terebrantis, and L. truncatum suggests that it has been horizontally transmitted by lateral gene transfer among these species, although the random presence of the intron may reflect multiple random loss events. A model is proposed describing the initial invasion of the group II intron in the rns gene of L. lundbergii by a LAGLIDADG homing endonuclease gene and subsequent evolution of this gene to recognize a novel DNA target site, which may now promote the mobility of the intron and homing endonuclease gene as a composite element.

Human Nek6 is a monomeric mostly globular kinase with an unfolded short N-terminal domain

BACKGROUND

The NIMA-related kinases (Neks) are widespread among eukaryotes. In mammalians they represent an evolutionarily conserved family of 11 serine/threonine kinases, with 40-45% amino acid sequence identity to the Aspergillus nidulans mitotic regulator NIMA within their catalytic domains. Neks have cell cycle-related functions and were recently described as related to pathologies, particularly cancer, consisting in potential chemotherapeutic targets. Human Nek6, -7 and -9 are involved in the control of mitotic spindle formation, acting together in a mitotic kinase cascade, but their mechanism of regulation remain elusive.

RESULTS

In this study we performed a biophysical and structural characterization of human Nek6 with the aim of obtaining its low resolution and homology models. SAXS experiments showed that hNek6 is a monomer of a mostly globular, though slightly elongated shape. Comparative molecular modeling together with disorder prediction analysis also revealed a flexible disordered N-terminal domain for hNek6, which we found to be important to mediate interactions with diverse partners. SEC-MALS experiments showed that hNek6 conformation is dependent on its activation/phosphorylation status, a higher phosphorylation degree corresponding to a bigger Stokes radius. Circular dichroism spectroscopy confirmed our in silico predictions of secondary structure content and thermal stability shift assays revealed a slightly higher stability of wild-type hNek6 compared to the activation loop mutant hNek6(S206A).

CONCLUSIONS

Our data present the first low resolution 3D structure of hNek6 protein in solution. SAXS, comparative modeling and SEC-MALS analysis revealed that hNek6 is a monomeric kinase of slightly elongated shape and a short unfolded N-terminal domain.

The EpsE flagellar clutch is bifunctional and synergizes with EPS biosynthesis to promote Bacillus subtilis biofilm formation

Many bacteria inhibit motility concomitant with the synthesis of an extracellular polysaccharide matrix and the formation of biofilm aggregates. In Bacillus subtilis biofilms, motility is inhibited by EpsE, which acts as a clutch on the flagella rotor to inhibit motility, and which is encoded within the 15 gene eps operon required for EPS production. EpsE shows sequence similarity to the glycosyltransferase family of enzymes, and we demonstrate that the conserved active site motif is required for EPS biosynthesis. We also screen for residues specifically required for either clutch or enzymatic activity and demonstrate that the two functions are genetically separable. Finally, we show that, whereas EPS synthesis activity is dominant for biofilm formation, both functions of EpsE synergize to stabilize cell aggregates and relieve selective pressure to abolish motility by genetic mutation. Thus, the transition from motility to biofilm formation may be governed by a single bifunctional enzyme.

Bacteria form persistent and antibiotic-resistant cell aggregates known as biofilms. Biofilms can form in environmental settings on plant and animal tissues, in industrial settings on pipes and the hulls of ships, and in clinical settings on catheters and medical devices. Biofilms are characterized by two features: the cells within the aggregates are non-motile, and they produce an extracellular polysaccharide (EPS) matrix. We have found a bifunctional enzyme EpsE that contributes to both features of biofilm formation in Bacillus subtilis. EpsE interacts with the flagella rotor to inhibit motility and also cooperates with other enzymes to synthesize the EPS matrix. Thus, the transition from motility to biofilm formation may be governed by a single bifunctional protein. In the past decade, research on biofilms has been focused on biofilm eradication. Understanding how cells transition into the biofilm state may provide additional approaches of preventing the formation of a biofilm in the first place.

Mutational dynamics of murine angiogenin duplicates

Background

Angiogenin (Ang) is a protein involved in angiogenesis by inducing the formation of blood vessels. The biomedical importance of this protein has come from findings linking mutations in Ang to cancer progression and neurodegenerative diseases. These findings highlight the evolutionary constrain on Ang amino acid sequence. However, previous studies comparing human Angiogenin with homologs from other phylogenetically related organisms have led to the conclusion that Ang presents a striking variability. Whether this variability has an adaptive value per se remains elusive. Understanding why many functional Ang paralogs have been preserved in mouse and rat and identifying functional divergence mutations at these copies may explain the relationship between mutations and function. In spite of the importance of testing this hypothesis from the evolutionarily and biomedical perspectives, this remains yet unaccomplished. Here we test the main mutational dynamics driving the evolution and function of Ang paralogs in mammals.

Results

We analysed the phylogenetic asymmetries between the different Ang gene copies in mouse and rat in the context of vertebrate Ang phylogeny. This analysis shows strong evidence in support of accelerated evolution in some Ang murine copies (mAng). This acceleration is not due to non-functionalisation because constraints on amino acid replacements remain strong. We identify many of the amino acid sites involved in signal localization and nucleotide binding by Ang to have evolved under diversifying selection. Compensatory effects of many of the mutations at these paralogs and their key structural location in or nearby important functional regions support a possible functional shift (functional divergence) in many Ang copies. Similarities between 3D-structural models for mAng copies suggest that their divergence is mainly functional.

Conclusions

We identify the main evolutionary dynamics shaping the variability of Angiogenin in vertebrates and highlight the plasticity of this protein after gene duplication. Our results suggest functional divergence among mAng paralogs. This puts forward mAng as a good system candidate for testing functional plasticity of such an important protein while stresses caution when using mouse as a model to infer the consequences of mutations in the single Ang copy of humans.

Identification and structural characterization of FYVE domain-containing proteins of Arabidopsis thaliana

BACKGROUND

FYVE domains have emerged as membrane-targeting domains highly specific for phosphatidylinositol 3-phosphate (PtdIns(3)P). They are predominantly found in proteins involved in various trafficking pathways. Although FYVE domains may function as individual modules, dimers or in partnership with other proteins, structurally, all FYVE domains share a fold comprising two small characteristic double-stranded beta-sheets, and a C-terminal alpha-helix, which houses eight conserved Zn2+ ion-binding cysteines. To date, the structural, biochemical, and biophysical mechanisms for subcellular targeting of FYVE domains for proteins from various model organisms have been worked out but plant FYVE domains remain noticeably under-investigated.

RESULTS

We carried out an extensive examination of all Arabidopsis FYVE domains, including their identification, classification, molecular modeling and biophysical characterization using computational approaches. Our classification of fifteen Arabidopsis FYVE proteins at the outset reveals unique domain architectures for FYVE containing proteins, which are not paralleled in other organisms. Detailed sequence analysis and biophysical characterization of the structural models are used to predict membrane interaction mechanisms previously described for other FYVE domains and their subtle variations as well as novel mechanisms that seem to be specific to plants.

CONCLUSIONS

Our study contributes to the understanding of the molecular basis of FYVE-based membrane targeting in plants on a genomic scale. The results show that FYVE domain containing proteins in plants have evolved to incorporate significant differences from those in other organisms implying that they play a unique role in plant signaling pathways and/or play similar/parallel roles in signaling to other organisms but use different protein players/signaling mechanisms.

E93R substitution of Escherichia coli FtsZ induces bundling of protofilaments, reduces GTPase activity, and impairs bacterial cytokinesis

Recently, we found that divalent calcium has no detectable effect on the assembly of Mycobacterium tuberculosis FtsZ (MtbFtsZ), whereas it strongly promoted the assembly of Escherichia coli FtsZ (EcFtsZ). While looking for potential calcium binding residues in EcFtsZ, we found a mutation (E93R) that strongly promoted the assembly of EcFtsZ. The mutation increased the stability and bundling of the FtsZ protofilaments and produced a dominating effect on the assembly of the wild type FtsZ (WT-FtsZ). Although E93R-FtsZ was found to bind to GTP similarly to the WT-FtsZ, it displayed lower GTPase activity than the WT-FtsZ. E93R-FtsZ complemented for its wild type counterpart as observed by a complementation test using JKD7-1/pKD3 cells. However, the bacterial cells became elongated upon overexpression of the mutant allele. We modeled the structure of E93R-FtsZ using the structures of MtbFtsZ/Methanococcus jannaschi FtsZ (MjFtsZ) dimers as templates. The MtbFtsZ-based structure suggests that the Arg(93)-Glu(138) salt bridge provides the additional stability, whereas the effect of mutation appears to be indirect (allosteric) if the EcFtsZ dimer is similar to that of MjFtsZ. The data presented in this study suggest that an increase in the stability of the FtsZ protofilaments is detrimental for the bacterial cytokinesis.

Cardiomyopathy-causing deletion K210 in cardiac troponin T alters phosphorylation propensity of sarcomeric proteins

Ca(2+) desensitization of myofilaments is indicated as a primary mechanism for the pathogenesis of familial dilated cardiomyopathy (DCM) associated with the deletion of lysine 210 (DeltaK210) in cardiac troponin T (cTnT). DeltaK210 knock-in mice closely recapitulate the clinical phenotypes documented in patients with this mutation. Considerable evidence supports the proposition that phosphorylation of cardiac sarcomeric proteins is a key modulator of function and may exacerbate the effect of the deletion. In this study we investigate the impact of K210 deletion on phosphorylation propensity of sarcomeric proteins. Analysis of cardiac myofibrils isolated from DeltaK210 hearts identified a decrease in phosphorylation of cTnI (46%), cTnT (30%) and MyBP-C (32%) compared with wild-type controls. Interestingly, immunoblot analyses with phospho-specific antibodies show augmented phosphorylation of cTnT-Thr(203) (28%) and decreased phosphorylation of cTnI-Ser(23/24) (41%) in mutant myocardium. In vitro kinase assays indicate that DeltaK210 increases phosphorylation propensity of cTnT-Thr(203) three-fold, without changing cTnI-Ser(23/24) phosphorylation. Molecular modeling of cTnT-DeltaK210 structure reveals changes in the electrostatic environment of cTnT helix (residues 203-224) that lead to a more basic environment around Thr(203), which may explain the enhanced PKC-dependent phosphorylation. In addition, yeast two-hybrid assays indicate that cTnT-DeltaK210 binds stronger to cTnI compared with cTnT-wt. Collectively, our observations suggest that cardiomyopathy-causing DeltaK210 has far-reaching effects influencing cTnI-cTnT binding and posttranslational modifications of key sarcomeric proteins.

A second actin-like MamK protein in Magnetospirillum magneticum AMB-1 encoded outside the genomic magnetosome island

Magnetotactic bacteria are able to swim navigating along geomagnetic field lines. They synthesize ferromagnetic nanocrystals that are embedded in cytoplasmic membrane invaginations forming magnetosomes. Regularly aligned in the cytoplasm along cytoskeleton filaments, the magnetosome chain effectively forms a compass needle bestowing on bacteria their magnetotactic behaviour. A large genomic island, conserved among magnetotactic bacteria, contains the genes potentially involved in magnetosome formation. One of the genes, mamK has been described as encoding a prokaryotic actin-like protein which when it polymerizes forms in the cytoplasm filamentous structures that provide the scaffold for magnetosome alignment. Here, we have identified a series of genes highly similar to the mam genes in the genome of Magnetospirillum magneticum AMB-1. The newly annotated genes are clustered in a genomic islet distinct and distant from the known magnetosome genomic island and most probably acquired by lateral gene transfer rather than duplication. We focused on a mamK-like gene whose product shares 54.5% identity with the actin-like MamK. Filament bundles of polymerized MamK-like protein were observed in vitro with electron microscopy and in vivo in E. coli cells expressing MamK-like-Venus fusions by fluorescence microscopy. In addition, we demonstrate that mamK-like is transcribed in AMB-1 wild-type and DeltamamK mutant cells and that the actin-like filamentous structures observed in the DeltamamK strain are probably MamK-like polymers. Thus MamK-like is a new member of the prokaryotic actin-like family. This is the first evidence of a functional mam gene encoded outside the magnetosome genomic island.

Heterotypic humoral and cellular immune responses following Norwalk virus infection

Norovirus immunity is poorly understood as the limited data available on protection after infection are often contradictory. In contrast to the more prominent GII noroviruses, GI norovirus infections are less frequent in outbreaks. The GI noroviruses display very complex patterns of heterotypic immune responses following infection, and many individuals are highly susceptible to reinfection. To study the immune responses and mechanisms of GI.1 persistence, we built structural models and recombinant virus-like particles (VLPs) of five GI strains: GI.1-1968, GI.1-2001, GI.2-1999, GI.3-1999, and GI.4-2000. Structural models of four GI genotype capsid P domain dimers suggested that intragenotype structural variation is limited, that the GI binding pocket is mostly preserved between genotypes, and that a conserved, surface-exposed epitope may allow for highly cross-reactive immune responses. GI VLPs bound to histo-blood group antigens (HBGAs) including fucose, Lewis, and A antigens. Volunteers infected with GI.1-1968 (n = 10) had significant increases between prechallenge and convalescent reactive IgG for all five GI VLPs measured by enzyme immunoassay. Potential cross-neutralization of GI VLPs was demonstrated by convalescent-phase serum cross-blockade of GI VLP-HBGA interaction. Although group responses were significant for all GI VLPs, each individual volunteer demonstrated a unique VLP blockade pattern. Further, peripheral blood mononuclear cells (PBMCs) were stimulated with each of the VLPs, and secretion of gamma interferon (IFN-gamma) was measured. As seen with blockade responses, IFN-gamma secretion responses differed by individual. Sixty percent responded to at least one GI VLP, with only two volunteers responding to GI.1 VLP. Importantly, four of five individuals with sufficient PBMCs for cross-reactivity studies responded more robustly to other GI VLPs. These data suggest that preexposure history and deceptive imprinting may complicate PBMC and B-cell immune responses in some GI.1-1968-challenged individuals and highlight a potential complication in the design of efficacious norovirus vaccines.

Evidence for an activation domain at the amino terminus of simian immunodeficiency virus Vpx

Vpx and Vpr are related lentiviral accessory proteins that enhance virus replication in macrophages and dendritic cells. Both proteins are packaged into virions and mediate their effects in the target cell through an interaction with an E3 ubiquitin ligase that contains DCAF1 and DDB1. When introduced into primary macrophages and dendritic cells in viruslike particles, Vpx can enhance the efficiency of a subsequent infection. Here, we confirm the ability of Vpx to enhance simian immunodeficiency virus (SIV) and human immunodeficiency virus type 1 (HIV-1) infection of macrophages up to 100-fold by using single-cycle reporter viruses and by pretreatment of the cells with Vpx-containing viruslike particles. Vpx was also active in differentiated THP-1 cells but not in other cell lines. Induction of an antiviral state in macrophages with type I interferon significantly magnified the effect of Vpx on HIV-1 infection, suggesting that Vpx helps the virus to overcome an inducible intracellular restriction. Quantitative PCR quantitation of SIV and HIV-1 reverse transcripts in newly infected macrophages showed that the block was at an early step in reverse transcription. In spite of its structural similarity, Vpr was inactive. This difference allowed us to map the functional domains of Vpx with a panel of Vpr/Vpx chimeras. Analysis of the chimeras demonstrated that the amino-terminal domain of Vpx is important for the enhancement of infection. Fine mapping of the region indicated that amino acids at positions 9, 12, and 15 to 17 were required. Although the mutants failed to enhance infection, they retained their ability to interact with DCAF1. These findings suggest that the Vpx amino terminus contains an activation domain that serves as the binding site for a cellular restriction factor.

Arginine, a key residue for the enhancing ability of an antifreeze protein of the beetle Dendroides canadensis

Antifreeze proteins (AFPs) can produce a difference between the nonequilibrium freezing point and the melting point, termed thermal hysteresis (TH). The TH activity of an antifreeze protein (AFP) depends on the specific AFP and its concentration as well as the presence of cosolutes including low molecular mass solutes and/or proteins. We recently identified series of carboxylates and polyols as efficient enhancers for an AFP from the beetle Dendroides canadensis. In this study, we chemically modified DAFP-1 using the arginine-specific reagent 1,2-cyclohexanedione. We demonstrated that 1,2-cyclohexanedione specifically modifies one arginine residue and the modified DAFP-1 loses its enhancing ability completely or partially in the presence of previously identified enhancers. The stronger the enhancement ability of the enhancer on the native DAFP-1, the stronger the enhancement effect of the enhancer on the modified DAFP-1. The weaker enhancers (e.g., glycerol) completely lose their enhancement effect on the modified DAFP-1 due to their inability to compete with 1,2-cyclohexanedione for the arginine residue. Regeneration of the arginine residue using hydroxylamine fully restored the enhancing ability of DAFP-1. These studies indicated that an arginine residue is critical for the enhancing ability of DAFP-1 and the guanidinium group of the arginine residue is important for its interaction with the enhancers, where the general mechanism of arginine-ligand interaction is borne. This work may initiate a complete mechanistic study of the enhancement effect in AFPs.

Kinetic and functional analysis of L-threonine kinase, the PduX enzyme of Salmonella enterica

The PduX enzyme of Salmonella enterica is an l-threonine kinase used for the de novo synthesis of coenzyme B(12) and the assimilation of cobyric acid. PduX with an N-terminal histidine tag (His(8)-PduX) was produced in Escherichiacoli and purified. The recombinant enzyme was soluble and active. Kinetic analysis indicated a steady-state Ordered Bi Bi complex mechanism in which ATP is the first substrate to bind. Based on a multiple sequence alignment of PduX homologues and other GHMP (galactokinase, homoserine kinase, mevalonate kinase, and phosphomevalonate kinase) family members, 14 PduX variants having changes at 10 conserved serine/threonine and aspartate/glutamate sites were constructed by site-directed mutagenesis. Each variant was produced in E. coli and purified. Comparison of the circular dichroism spectra and kinetic properties of the PduX variants with those of the wild-type enzyme indicated that Glu-24 and Asp-135 are needed for proper folding, Ser-99 and Glu-132 are used for ATP binding, and Ser-253 and Ser-255 are critical to l-threonine binding whereas Ser-100 is essential to catalysis, but its precise role is uncertain. The studies reported here are the first to investigate the kinetic and catalytic mechanisms of l-threonine kinase from any organism.

The BARD1 C-terminal domain structure and interactions with polyadenylation factor CstF-50

The BARD1 N-terminal RING domain binds BRCA1 while the BARD1 C-terminal ankyrin and tandem BRCT repeat domains bind CstF-50 to modulate mRNA processing and RNAP II stability in response to DNA damage. Here we characterize the BARD1 structural biochemistry responsible for CstF-50 binding. The crystal structure of the BARD1 BRCT domain uncovers a degenerate phosphopeptide binding pocket lacking the key arginine required for phosphopeptide interactions in other BRCT proteins. Small angle X-ray scattering together with limited proteolysis results indicates that ankyrin and BRCT domains are linked by a flexible tether and do not adopt a fixed orientation relative to one another. Protein pull-down experiments utilizing a series of purified BARD1 deletion mutants indicate that interactions between the CstF-50 WD-40 domain and BARD1 involve the ankyrin-BRCT linker but do not require ankyrin or BRCT domains. The structural plasticity imparted by the ANK-BRCT linker helps to explain the regulated assembly of different protein BARD1 complexes with distinct functions in DNA damage signaling including BARD1-dependent induction of apoptosis plus p53 stabilization and interactions. BARD1 architecture and plasticity imparted by the ANK-BRCT linker are suitable to allow the BARD1 C-terminus to act as a hub with multiple binding sites to integrate diverse DNA damage signals directly to RNA polymerase.

Functional characterization of a human aquaporin 0 mutation that leads to a congenital dominant lens cataract

The aquaporin (AQP) transmembrane proteins facilitate the movement of water across the plasma membrane. In the lens, AQP0 is expressed in fiber cells and AQP1 in the epithelium. Recently, two individuals were identified with congenital polymorphic autosomal dominant cataract, due to a single nucleotide base deletion mutation in the lens AQP0. The deletion modified the reading frame resulting in the addition of a premature stop codon. In the present study, we examined the water permeability properties, trafficking and dominant negative effects as well as cytotoxicity due to the mutant AQP0 (Delta213-AQP0) protein. The membrane water permeability (P(w)) of Delta213-AQP0 expressing oocytes (14+/-1 microm/s) was significantly lower than those expressing WT-AQP0 (25+/-3 microm/s). P(w) of water injected control oocytes was 13+/-2 microm/s. Co-expression of WT-AQP0 with Delta213-AQP0 significantly lowered the P(w) (18+/-3 microm/s) compared to WT-AQP0. With or without the EGFP tag, WT-AQP0 protein localized in the plasma membranes of oocytes and cultured cells whereas Delta213-AQP0 was retained in the ER. Forster Resonance Energy Transfer (FRET) showed that WT-AQP0 partly localized with the co-expressed Delta213-AQP0. Co-localization studies suggest that the mutant AQP0 gained its dominant function by trapping the WT-AQP0 in the ER through hetero-oligomerization. Incubating the cells with chemical chaperones, namely, TMAO and DMSO, did not correct the folding/trafficking defects. Cell death in the Delta213-AQP0 expressing cells was due to necrosis caused by the accumulation of Delta213-AQP0 protein in the ER in cytotoxic proportions. The data show that replacement of the distal end of the 6th TM domain and the C-terminal domain of AQP0 due to the deletion mutation resulted in the impairment of cell membrane P(w), localization of the mutant protein in the ER without trafficking to the plasma membrane, and cytotoxicity due to the accumulation of the mutant protein. Cataracts in patients with this mutation might have resulted from the above mentioned consequences.

A novel mutation in murine hepatitis virus nsp5, the viral 3C-like proteinase, causes temperature-sensitive defects in viral growth and protein processing

Sequencing and reversion analysis of murine hepatitis virus (MHV) temperature-sensitive (ts) viruses has identified putative ts mutations in the replicase nonstructural proteins (nsp's) of these coronaviruses. In this study, reverse transcriptase PCR sequencing of the RNA genome of an isolate of the MHV ts virus Alb ts6, referred to as Alb/ts/nsp5/V148A, identified a putative ts mutation in nsp5 (T10651C, Val148Ala), the viral 3C-like proteinase (3CLpro). The introduction of the T10651C mutation into the infectious MHV clone resulted in the recovery of a mutant virus, the nsp5/V148A virus, that demonstrated reduced growth and nsp5 proteinase activity identical to that of Alb/ts/nsp5/V148A at the nonpermissive temperature. Sequence analysis of 40 degrees C revertants of Alb/ts/nsp5/V148A identified primary reversion to Ala148Val in nsp5, as well as two independent second-site mutations resulting in Ser133Asn and His134Tyr substitutions in nsp5. The introduction of the Ser133Asn or His134Tyr substitution into the cloned nsp5/V148A mutant virus background resulted in the recovery of viruses with increased growth fitness and the partial restoration of nsp5 activity at the nonpermissive temperature. Modeling of the nsp5 structure of Alb/ts/nsp5/V148A predicted that the Val148Ala mutation alters residue 148 interactions with residues of the substrate binding S1 subsite of the nsp5 active-site cavity. This study identifies novel residues in nsp5 that may be important for regulating substrate specificity and nsp5 proteinase activity.

Mechanisms of GII.4 norovirus persistence in human populations

BACKGROUND

Noroviruses are the leading cause of viral acute gastroenteritis in humans, noted for causing epidemic outbreaks in communities, the military, cruise ships, hospitals, and assisted living communities. The evolutionary mechanisms governing the persistence and emergence of new norovirus strains in human populations are unknown. Primarily organized by sequence homology into two major human genogroups defined by multiple genoclusters, the majority of norovirus outbreaks are caused by viruses from the GII.4 genocluster, which was first recognized as the major epidemic strain in the mid-1990s. Previous studies by our laboratory and others indicate that some noroviruses readily infect individuals who carry a gene encoding a functional alpha-1,2-fucosyltransferase (FUT2) and are designated "secretor-positive" to indicate that they express ABH histo-blood group antigens (HBGAs), a highly heterogeneous group of related carbohydrates on mucosal surfaces. Individuals with defects in the FUT2 gene are termed secretor-negative, do not express the appropriate HBGA necessary for docking, and are resistant to Norwalk infection. These data argue that FUT2 and other genes encoding enzymes that regulate processing of the HBGA carbohydrates function as susceptibility alleles. However, secretor-negative individuals can be infected with other norovirus strains, and reinfection with the GII.4 strains is common in human populations. In this article, we analyze molecular mechanisms governing GII.4 epidemiology, susceptibility, and persistence in human populations.

METHODS AND FINDINGS

Phylogenetic analyses of the GII.4 capsid sequences suggested an epochal evolution over the last 20 y with periods of stasis followed by rapid evolution of novel epidemic strains. The epidemic strains show a linear relationship in time, whereby serial replacements emerge from the previous cluster. Five major evolutionary clusters were identified, and representative ORF2 capsid genes for each cluster were expressed as virus-like particles (VLPs). Using salivary and carbohydrate-binding assays, we showed that GII.4 VLP-carbohydrate ligand binding patterns have changed over time and include carbohydrates regulated by the human FUT2 and FUT3 pathways, suggesting that strain sensitivity to human susceptibility alleles will vary. Variation in surface-exposed residues and in residues that surround the fucose ligand interaction domain suggests that antigenic drift may promote GII.4 persistence in human populations. Evidence supporting antigenic drift was obtained by measuring the antigenic relatedness of GII.4 VLPs using murine and human sera and demonstrating strain-specific serologic and carbohydrate-binding blockade responses. These data suggest that the GII.4 noroviruses persist by altering their HBGA carbohydrate-binding targets over time, which not only allows for escape from highly penetrant host susceptibility alleles, but simultaneously allows for immune-driven selection in the receptor-binding region to facilitate escape from protective herd immunity.

CONCLUSIONS

Our data suggest that the surface-exposed carbohydrate ligand binding domain in the norovirus capsid is under heavy immune selection and likely evolves by antigenic drift in the face of human herd immunity. Variation in the capsid carbohydrate-binding domain is tolerated because of the large repertoire of similar, yet distinct HBGA carbohydrate receptors available on mucosal surfaces that could interface with the remodeled architecture of the capsid ligand-binding pocket. The continuing evolution of new replacement strains suggests that, as with influenza viruses, vaccines could be targeted that protect against norovirus infections, and that continued epidemiologic surveillance and reformulations of norovirus vaccines will be essential in the control of future outbreaks.

Preservation of protein clefts in comparative models

BACKGROUND

Comparative, or homology, modelling of protein structures is the most widely used prediction method when the target protein has homologues of known structure. Given that the quality of a model may vary greatly, several studies have been devoted to identifying the factors that influence modelling results. These studies usually consider the protein as a whole, and only a few provide a separate discussion of the behaviour of biologically relevant features of the protein. Given the value of the latter for many applications, here we extended previous work by analysing the preservation of native protein clefts in homology models. We chose to examine clefts because of their role in protein function/structure, as they are usually the locus of protein-protein interactions, host the enzymes' active site, or, in the case of protein domains, can also be the locus of domain-domain interactions that lead to the structure of the whole protein.

RESULTS

We studied how the largest cleft of a protein varies in comparative models. To this end, we analysed a set of 53507 homology models that cover the whole sequence identity range, with a special emphasis on medium and low similarities. More precisely we examined how cleft quality - measured using six complementary parameters related to both global shape and local atomic environment, depends on the sequence identity between target and template proteins. In addition to this general analysis, we also explored the impact of a number of factors on cleft quality, and found that the relationship between quality and sequence identity varies depending on cleft rank amongst the set of protein clefts (when ordered according to size), and number of aligned residues.

CONCLUSION

We have examined cleft quality in homology models at a range of seq.id. levels. Our results provide a detailed view of how quality is affected by distinct parameters and thus may help the user of comparative modelling to determine the final quality and applicability of his/her cleft models. In addition, the large variability in model quality that we observed within each sequence bin, with good models present even at low sequence identities (between 20% and 30%), indicates that properly developed identification methods could be used to recover good cleft models in this sequence range.

In silico identification of functional divergence between the multiple groEL gene paralogs in Chlamydiae

Background

Heat-shock proteins are specialized molecules performing different and essential roles in the cell including protein degradation, folding and trafficking. GroEL is a 60 Kda heat-shock protein ubiquitous in bacteria and has been regarded as an important molecule implicated in chronic inflammatory processes caused by Chlamydiae infections. GroEL in Chlamydiae became duplicated at the origin of the Chlamydiae lineage presenting three distinct molecular chaperones, namely the original protein GroEL1 (Ct110), and its paralogous proteins GroEL2 (Ct604) and GroEL3 (Ct755). These chaperones present differential and independent expressions during the different stages of Chlamydiae infections and have been suggested to present differential physiological and regulatory roles.

Results

In this comprehensive in silico study we show that GroEL protein paralogs have diverged functionally after the different gene duplication events and that this divergence has occurred mainly between GroEL3 and GroEL1. GroEL2 presents an intermediate functional divergence pattern from GroEL1. Our results point to the different protein-protein interaction patterns between GroEL paralogs and known GroEL protein clients supporting their functional divergence after groEL gene duplication. Analysis of selective constraints identifies periods of adaptive evolution after gene duplication that led to the fixation of amino acid replacements in GroEL protein domains involved in the interaction with GroEL protein clients.

Conclusion

We demonstrate that GroEL protein copies in Chlamydiae species have diverged functionally after the gene duplication events. We also show that functional divergence has occurred in important functional regions of these GroEL proteins and that very probably have affected the ancestral GroEL regulatory role and protein-protein interaction patterns with GroEL client proteins. Most of the amino acid replacements that have affected interaction with protein clients and that were responsible for the functional divergence between GroEL paralogs were fixed by adaptive evolution after the groEL gene duplication events.

Breakage-reunion domain of Streptococcus pneumoniae topoisomerase IV: crystal structure of a gram-positive quinolone target

The 2.7 A crystal structure of the 55-kDa N-terminal breakage-reunion domain of topoisomerase (topo) IV subunit A (ParC) from Streptococcus pneumoniae, the first for the quinolone targets from a gram-positive bacterium, has been solved and reveals a 'closed' dimer similar in fold to Escherichia coli DNA gyrase subunit A (GyrA), but distinct from the 'open' gate structure of Escherichia coli ParC. Unlike GyrA whose DNA binding groove is largely positively charged, the DNA binding site of ParC exhibits a distinct pattern of alternating positively and negatively charged regions coincident with the predicted positions of the grooves and phosphate backbone of DNA. Based on the ParC structure, a new induced-fit model for sequence-specific recognition of the gate (G) segment by ParC has been proposed. These features may account for the unique DNA recognition and quinolone targeting properties of pneumococcal type II topoisomerases compared to their gram-negative counterparts.

In silico pathway reconstruction: Iron-sulfur cluster biogenesis in Saccharomyces cerevisiae

Background

Current advances in genomics, proteomics and other areas of molecular biology make the identification and reconstruction of novel pathways an emerging area of great interest. One such class of pathways is involved in the biogenesis of Iron-Sulfur Clusters (ISC).

Results

Our goal is the development of a new approach based on the use and combination of mathematical, theoretical and computational methods to identify the topology of a target network. In this approach, mathematical models play a central role for the evaluation of the alternative network structures that arise from literature data-mining, phylogenetic profiling, structural methods, and human curation. As a test case, we reconstruct the topology of the reaction and regulatory network for the mitochondrial ISC biogenesis pathway in S. cerevisiae. Predictions regarding how proteins act in ISC biogenesis are validated by comparison with published experimental results. For example, the predicted role of Arh1 and Yah1 and some of the interactions we predict for Grx5 both matches experimental evidence. A putative role for frataxin in directly regulating mitochondrial iron import is discarded from our analysis, which agrees with also published experimental results. Additionally, we propose a number of experiments for testing other predictions and further improve the identification of the network structure.

Conclusion

We propose and apply an iterative in silico procedure for predictive reconstruction of the network topology of metabolic pathways. The procedure combines structural bioinformatics tools and mathematical modeling techniques that allow the reconstruction of biochemical networks. Using the Iron Sulfur cluster biogenesis in S. cerevisiae as a test case we indicate how this procedure can be used to analyze and validate the network model against experimental results. Critical evaluation of the obtained results through this procedure allows devising new wet lab experiments to confirm its predictions or provide alternative explanations for further improving the models.

Characterization of the Moraxella catarrhalis opa-like protein, OlpA, reveals a phylogenetically conserved family of outer membrane proteins

Moraxella catarrhalis is a human-restricted pathogen that can cause respiratory tract infections. In this study, we identified a previously uncharacterized 24-kDa outer membrane protein with a high degree of similarity to Neisseria Opa protein adhesins, with a predicted beta-barrel structure consisting of eight antiparallel beta-sheets with four surface-exposed loops. In striking contrast to the antigenically variable Opa proteins, the M. catarrhalis Opa-like protein (OlpA) is highly conserved and constitutively expressed, with 25 of 27 strains corresponding to a single variant. Protease treatment of intact bacteria and isolation of outer membrane vesicles confirm that the protein is surface exposed yet does not bind host cellular receptors recognized by neisserial Opa proteins. Genome-based analyses indicate that OlpA and Opa derive from a conserved family of proteins shared by a broad array of gram-negative bacteria.

Structural biology and bioinformatics in drug design: opportunities and challenges for target identification and lead discovery

Impressive progress in genome sequencing, protein expression and high-throughput crystallography and NMR has radically transformed the opportunities to use protein three-dimensional structures to accelerate drug discovery, but the quantity and complexity of the data have ensured a central place for informatics. Structural biology and bioinformatics have assisted in lead optimization and target identification where they have well established roles; they can now contribute to lead discovery, exploiting high-throughput methods of structure determination that provide powerful approaches to screening of fragment binding.

SCER and SCAN: two novel high-prevalence antigens in the Scianna blood group system

BACKGROUND

More than 20 years ago, two probands were described whose red blood cells (RBCs) typed Sc:1,-2,3. Their serum samples contained alloantibodies reactive with all RBCs tested except those of the Sc:-1,-2,-3 phenotype. Cloning of the Scianna gene allowed us to determine the molecular bases of these samples.

STUDY DESIGN AND METHODS

In a collaborative effort, the two probands' samples and also two Sc:-1,-2,-3 samples were obtained from frozen storage. All 11 SC (ERMAP) exons and their flanking regions were sequenced.

RESULTS

The two probands with antibodies to Scianna-related antigens were homozygous, respectively, for an ERMAP(R81Q) allele caused by a G to A substitution at nucleotide 242 in the ERMAP gene and for an ERMAP(H26Y,G35S) allele, in which the G35S substitution was caused by a G to A substitution at nucleotide 103. Two patients with the Sc:-1,-2,-3 phenotype both carried ERMAP(R332X) alleles caused by a C to T substitution at nucleotide 994 that differed at one nucleotide position in the noncoding region of exon 11. In eight samples carrying orphan low-prevalence antigens, no ERMAP variants were detected that could be implicated in Scianna antigen expression.

CONCLUSION

SCER and SCAN expanded the Scianna blood group system to seven antigens, have been assigned the ISBT numbers 013.006 (Sc6) and 013.007 (Sc7), and were associated with ERMAP(R81Q) and ERMAP(G35S) proteins, respectively. ERMAP(R332X) is a second molecular basis for the Sc(null) phenotype. The eight low-prevalence antigens By, To(a), Pt(a), Re(a), Je(a), Li(a), SARA, and Sk(a) do not belong to the Scianna blood group system.

Minimalist explicit solvation models for surface loops in proteins

We have performed molecular dynamics simulations of protein surface loops solvated by explicit water, where a prime focus of the study is the small numbers (e.g., ~100) of explicit water molecules employed. The models include only part of the protein (typically 500 - 1000 atoms), and the water molecules are restricted to a region surrounding the loop. In this study, the number of water molecules (N(w)) is systematically varied, and convergence with large N(w) is monitored to reveal N(w)(min), the minimum number required for the loop to exhibit realistic (fully hydrated) behavior. We have also studied protein surface coverage, as well as diffusion and residence times for water molecules as a function of N(w). A number of other modeling parameters are also tested. These include the number of environmental protein atoms explicitly considered in the model, as well as two ways to constrain the water molecules to the vicinity of the loop (where we find one of these methods to perform better when N(w) is small). The results (for RMSD and its fluctuations for four loops) are further compared to much larger, fully solvated systems (using ~10,000 water molecules under periodic boundary conditions and Ewald electrostatics), and to results for the GBSA implicit solvation model. We find that the loop backbone can stabilize with a surprisingly small number of water molecules (as low as 5 molecules per amino acid residue). The side chains of the loop require somewhat larger N(w), where the atomic fluctuations become too small if N(w) is further reduced. Thus, in general, we find adequate hydration to occur at roughly 12 water molecules per residue. This is an important result, because at this hydration level, computational times are comparable to those required for GBSA. Therefore these "minimalist explicit models" can provide a viable and potentially more accurate alternative. The importance of protein loop modeling is discussed in the context of these, and other, loop models, along with other challenges including the relevance of appropriate free energy simulation methodology for assessment of conformational stability.

Identification of the rctA gene, which is required for repression of conjugative transfer of rhizobial symbiotic megaplasmids

An analysis of the conjugative transfer of pRetCFN42d, the symbiotic plasmid (pSym) of Rhizobium etli, has revealed a novel gene, rctA, as an essential element of a regulatory system for silencing the conjugative transfer of R. etli pSym by repressing the transcription of conjugal transfer genes in standard laboratory media. The rctA gene product lacks sequence conservation with other proteins of known function but may belong to the winged-helix DNA-binding subfamily of transcriptional regulators. Similar to that of many transcriptional repressors, rctA transcription seems to be positively autoregulated. rctA expression is greatly reduced upon overexpression of another gene, rctB, previously identified as a putative activator of R. etli pSym conjugal transfer. Thus, rctB seems to counteract the repressive action of rctA. rctA homologs are present in at least three other bacterial genomes within the order Rhizobiales, where they are invariably located adjacent to and divergently transcribed from putative virB-like operons. We show that similar to that of R. etli pSym, conjugative transfer of the 1.35-Mb symbiotic megaplasmid A of Sinorhizobium meliloti is also subjected to the inhibitory action of rctA. Our data provide strong evidence that the R. etli and S. meliloti pSym plasmids are indeed self-conjugative plasmids and that this property would only be expressed under optimal, as yet unknown conditions that entail inactivation of the rctA function. The rctA gene seems to represent novel but probably widespread regulatory systems controlling the transfer of conjugative elements within the order Rhizobiales.

The modeled structure of the RNA dependent RNA polymerase of GBV-C virus suggests a role for motif E in Flaviviridae RNA polymerases

Background

The Flaviviridae virus family includes major human and animal pathogens. The RNA dependent RNA polymerase (RdRp) plays a central role in the replication process, and thus is a validated target for antiviral drugs. Despite the increasing structural and enzymatic characterization of viral RdRps, detailed molecular replication mechanisms remain unclear. The hepatitis C virus (HCV) is a major human pathogen difficult to study in cultured cells. The bovine viral diarrhea virus (BVDV) is often used as a surrogate model to screen antiviral drugs against HCV. The structure of BVDV RdRp has been recently published. It presents several differences relative to HCV RdRp. These differences raise questions about the relevance of BVDV as a surrogate model, and cast novel interest on the "GB" virus C (GBV-C). Indeed, GBV-C is genetically closer to HCV than BVDV, and can lead to productive infection of cultured cells. There is no structural data for the GBV-C RdRp yet.

Results

We show in this study that the GBV-C RdRp is closest to the HCV RdRp. We report a 3D model of the GBV-C RdRp, developed using sequence-to-structure threading and comparative modeling based on the atomic coordinates of the HCV RdRp structure. Analysis of the predicted structural features in the phylogenetic context of the RNA polymerase family allows rationalizing most of the experimental data available. Both available structures and our model are explored to examine the catalytic cleft, allosteric and substrate binding sites.

Conclusion

Computational methods were used to infer evolutionary relationships and to predict the structure of a viral RNA polymerase. Docking a GTP molecule into the structure allows defining a GTP binding pocket in the GBV-C RdRp, such as that of BVDV. The resulting model suggests a new proposition for the mechanism of RNA synthesis, and may prove useful to design new experiments to implement our knowledge on the initiation mechanism of RNA polymerases.

Mutations in the Saccharomyces cerevisiae LSM1 gene that affect mRNA decapping and 3' end protection

The decapping of eukaryotic mRNAs is a key step in their degradation. The heteroheptameric Lsm1p-7p complex is a general activator of decapping and also functions in protecting the 3' ends of deadenylated mRNAs from a 3'-trimming reaction. Lsm1p is the unique member of the Lsm1p-7p complex, distinguishing that complex from the functionally different Lsm2p-8p complex. To understand the function of Lsm1p, we constructed a series of deletion and point mutations of the LSM1 gene and examined their effects on phenotype. These studies revealed the following: (i) Mutations affecting the predicted RNA-binding and inter-subunit interaction residues of Lsm1p led to impairment of mRNA decay, suggesting that the integrity of the Lsm1p-7p complex and the ability of the Lsm1p-7p complex to interact with mRNA are important for mRNA decay function; (ii) mutations affecting the predicted RNA contact residues did not affect the localization of the Lsm1p-7p complex to the P-bodies; (iii) mRNA 3'-end protection could be indicative of the binding of the Lsm1p-7p complex to the mRNA prior to activation of decapping, since all the mutants defective in mRNA 3' end protection were also blocked in mRNA decay; and (iv) in addition to the Sm domain, the C-terminal domain of Lsm1p is also important for mRNA decay function.

Recognition of phosphorylated-Smad2-containing complexes by a novel Smad interaction motif

Transforming growth factor beta (TGF-beta) superfamily members signal via complexes of activated Smads, comprising phosphorylated receptor-regulated Smads, such as Smad2 and Smad3, and Smad4. These complexes are recruited to DNA by specific transcription factors. The forkhead/winged-helix transcription factors, XFast-1/XFoxH1a and XFast-3/XFoxH1b, bind an activated Smad heterotrimer comprising two Smad2s and one Smad4. Here we identify a novel Smad2 interaction motif, the Fast/FoxH1 motif (FM), present in all known Fast/FoxH1 family members, N-terminal to the common Smad interaction motif (SIM). The FM is necessary and sufficient to bind active Smad2/Smad4 complexes. The FM differs from the SIM since it discriminates between Smad2 and Smad3, and moreover only binds phosphorylated Smad2 in the context of activated Smad complexes. It is the first Smad interaction motif with this property. Site-directed mutagenesis indicates that the binding site for the FM on a Smad2/Smad4 heterotrimer is a hydrophobic pocket that incorporates the Smad/Smad interface. We demonstrate that the presence of an FM and SIM in the Fast/FoxH1 proteins allows them to compete efficiently for activated Smad2/Smad4 complexes with transcription factors such as Mixer that only contain a SIM. This establishes a hierarchy of Smad-interacting transcription factors, determined by their affinity for active Smad complexes.

Significance of conservative asparagine residues in the thermal hysteresis activity of carrot antifreeze protein

The approximately 24-amino-acid leucine-rich tandem repeat motif (PXXXXXLXXLXXLXLSXNXLXGXI) of carrot antifreeze protein comprises most of the processed protein and should contribute at least partly to the ice-binding site. Structural predictions using publicly available online sources indicated that the theoretical three-dimensional model of this plant protein includes a 10-loop beta-helix containing the approximately 24-amino-acid tandem repeat. This theoretical model indicated that conservative asparagine residues create putative ice-binding sites with surface complementarity to the 1010 prism plane of ice. We used site-specific mutagenesis to test the importance of these residues, and observed a distinct loss of thermal hysteresis activity when conservative asparagines were replaced with valine or glutamine, whereas a large increase in thermal hysteresis was observed when phenylalanine or threonine residues were replaced with asparagine, putatively resulting in the formation of an ice-binding site. These results confirmed that the ice-binding site of carrot antifreeze protein consists of conservative asparagine residues in each beta-loop. We also found that its thermal hysteresis activity is directly correlated with the length of its asparagine-rich binding site, and hence with the size of its ice-binding face.

Molecular modeling of the membrane targeting of phospholipase C pleckstrin homology domains

Phospholipases C (PLCs) reversibly associate with membranes to hydrolyze phosphatidylinositol-4, 5-bisphosphate (PI[4,5]P(2)) and comprise four main classes: beta, gamma, delta, and epsilon. Most eukaryotic PLCs contain a single, N-terminal pleckstrin homology (PH) domain, which is thought to play an important role in membrane targeting. The structure of a single PLC PH domain, that from PLCdelta1, has been determined; this PH domain binds PI(4,5)P(2) with high affinity and stereospecificity and has served as a paradigm for PH domain functionality. However, experimental studies demonstrate that PH domains from different PLC classes exhibit diverse modes of membrane interaction, reflecting the dissimilarity in their amino acid sequences. To elucidate the structural basis for their differential membrane-binding specificities, we modeled the three-dimensional structures of all mammalian PLC PH domains by using bioinformatic tools and calculated their biophysical properties by using continuum electrostatic approaches. Our computational analysis accounts for a large body of experimental data, provides predictions for those PH domains with unknown functions, and indicates functional roles for regions other than the canonical lipid-binding site identified in the PLCdelta1-PH structure. In particular, our calculations predict that (1). members from each of the four PLC classes exhibit strikingly different electrostatic profiles than those ordinarily observed for PH domains in general, (2). nonspecific electrostatic interactions contribute to the membrane localization of PLCdelta-, PLCgamma-, and PLCbeta-PH domains, and (3). phosphorylation regulates the interaction of PLCbeta-PH with its effectors through electrostatic repulsion. Our molecular models for PH domains from all of the PLC classes clearly demonstrate how a common structural fold can serve as a scaffold for a wide range of surface features and biophysical properties that support distinctive functional roles.

EVA: Evaluation of protein structure prediction servers

EVA (http://cubic.bioc.columbia.edu/eva/) is a web server for evaluation of the accuracy of automated protein structure prediction methods. The evaluation is updated automatically each week, to cope with the large number of existing prediction servers and the constant changes in the prediction methods. EVA currently assesses servers for secondary structure prediction, contact prediction, comparative protein structure modelling and threading/fold recognition. Every day, sequences of newly available protein structures in the Protein Data Bank (PDB) are sent to the servers and their predictions are collected. The predictions are then compared to the experimental structures once a week; the results are published on the EVA web pages. Over time, EVA has accumulated prediction results for a large number of proteins, ranging from hundreds to thousands, depending on the prediction method. This large sample assures that methods are compared reliably. As a result, EVA provides useful information to developers as well as users of prediction methods.

Reversal of DNA alkylation damage by two human dioxygenases

The Escherichia coli AlkB protein protects against the cytotoxicity of methylating agents by repair of the DNA lesions 1-methyladenine and 3-methylcytosine, which are generated in single-stranded stretches of DNA. AlkB is an alpha-ketoglutarate- and Fe(II)-dependent dioxygenase that oxidizes the relevant methyl groups and releases them as formaldehyde. Here, we identify two human AlkB homologs, ABH2 and ABH3, by sequence and fold similarity, functional assays, and complementation of the E. coli alkB mutant phenotype. The levels of their mRNAs do not appear to correlate with cell proliferation but tissue distributions are different. Both enzymes remove 1-methyladenine and 3-methylcytosine from methylated polynucleotides in an alpha-ketoglutarate-dependent reaction, and act by direct damage reversal with the regeneration of the unsubstituted bases. AlkB, ABH2, and ABH3 can also repair 1-ethyladenine residues in DNA with the release of acetaldehyde.

Structural characterization of the human proteome

This paper reports an analysis of the encoded proteins (the proteome) of the genomes of human, fly, worm, yeast, and representatives of bacteria and archaea in terms of the three-dimensional structures of their globular domains together with a general sequence-based study. We show that 39% of the human proteome can be assigned to known structures. We estimate that for 77% of the proteome, there is some functional annotation, but only 26% of the proteome can be assigned to standard sequence motifs that characterize function. Of the human protein sequences, 13% are transmembrane proteins, but only 3% of the residues in the proteome form membrane-spanning regions. There are substantial differences in the composition of globular domains of transmembrane proteins between the proteomes we have analyzed. Commonly occurring structural superfamilies are identified within the proteome. The frequencies of these superfamilies enable us to estimate that 98% of the human proteome evolved by domain duplication, with four of the 10 most duplicated superfamilies specific for multicellular organisms. The zinc-finger superfamily is massively duplicated in human compared to fly and worm, and occurrence of domains in repeats is more common in metazoa than in single cellular organisms. Structural superfamilies over- and underrepresented in human disease genes have been identified. Data and results can be downloaded and analyzed via web-based applications at http://www.sbg.bio.ic.ac.uk.

Mechanochemical ATPases and transcriptional activation

Transcriptional activator proteins that act upon the sigma54-containing form of the bacterial RNA polymerase belong to the extensive AAA+ superfamily of ATPases, members of which are found in all three kingdoms of life and function in diverse cellular processes, often via chaperone-like activities. Formation and collapse of the transition state of ATP for hydrolysis appears to engender the interaction of the activator proteins with sigma54 and leads to the protein structural transitions needed for RNA polymerase to isomerize and engage with the DNA template strand. The common oligomeric structures of AAA+ proteins and the creation of the active site for ATP hydrolysis between protomers suggest that the critical changes in protomer structure required for productive interactions with sigma54-holoenzyme occur as a consequence of sensing the state of the gamma-phosphate of ATP. Depending upon the form of nucleotide bound, different functional states of the activator are created that have distinct substrate and chaperone-like binding activities. In particular, interprotomer ATP interactions rely upon the use of an arginine finger, a situation reminiscent of GTPase-activating proteins.

Different Smad2 partners bind a common hydrophobic pocket in Smad2 via a defined proline-rich motif

Transforming growth factor-beta (TGF-beta)/activin-induced Smad2/Smad4 complexes are recruited to different promoter elements by transcription factors, such as Fast-1 or the Mix family proteins Mixer and Milk, through a direct interaction between Smad2 and a common Smad interaction motif (SIM) in the transcription factors. Here we identify residues in the SIM critical for Mixer-Smad2 interaction and confirm their functional importance by demonstrating that only Xenopus and zebrafish Mix family members containing a SIM with all the correct critical residues can bind Smad2 and mediate TGF-beta-induced transcriptional activation in vivo. We identify significant sequence similarity between the SIM and the Smad-binding domain (SBD) of the membrane-associated protein SARA (Smad anchor for receptor activation). Molecular modelling, supported by mutational analyses of Smad2 and the SIM and the demonstration that the SARA SBD competes directly with the SIM for binding to Smad2, indicates that the SIM binds Smad2 in the same hydrophobic pocket as does the proline-rich rigid coil region of the SARA SBD. Thus, different Smad2 partners, whether cytoplasmic or nuclear, interact with the same binding pocket in Smad2 through a common proline-rich motif.

A cofactor-dependent phosphoglycerate mutase homolog from Bacillus stearothermophilus is actually a broad specificity phosphatase

The distribution of phosphoglycerate mutase (PGM) activity in bacteria is complex, with some organisms possessing both a cofactor-dependent and a cofactor-independent PGM and others having only one of these enzymes. Although Bacillus species contain only a cofactor-independent PGM, genes homologous to those encoding cofactor-dependent PGMs have been detected in this group of bacteria, but in at least one case the encoded protein lacks significant PGM activity. Here we apply sequence analysis, molecular modeling, and enzymatic assays to the cofactor-dependent PGM homologs from B. stearothermophilus and B. subtilis, and show that these enzymes are phosphatases with broad substrate specificity. Homologs from other gram-positive bacteria are also likely to possess phosphatase activity. These studies clearly show that the exploration of genomic sequences through three-dimensional modeling is capable of producing useful predictions regarding function. However, significant methodological improvements will be needed before such analysis can be carried out automatically.

CODA: a combined algorithm for predicting the structurally variable regions of protein models

CODA, an algorithm for predicting the variable regions in proteins, combines FREAD a knowledge based approach, and PETRA, which constructs the region ab initio. FREAD selects from a database of protein structure fragments with environmentally constrained substitution tables and other rule-based filters. FREAD was parameterized and tested on over 3000 loops. The average root mean square deviation ranged from 0.78 A for three residue loops to 3.5 A for eight residue loops on a nonhomologous test set. CODA clusters the predictions from the two independent programs and makes a consensus prediction that must pass a set of rule-based filters. CODA was parameterized and tested on two unrelated separate sets of structures that were nonhomologous to one another and those found in the FREAD database. The average root mean square deviation in the test set ranged from 0.76 A for three residue loops to 3.09 A for eight residue loops. CODA shows a general improvement in loop prediction over PETRA and FREAD individually. The improvement is far more marked for lengths six and upward, probably as the predictive power of PETRA becomes more important. CODA was further tested on several model structures to determine its applicability to the modeling situation. A web server of CODA is available at http://www-cryst.bioc.cam.ac.uk/~charlotte/Coda/search_coda.html.

Construction of a full three-dimensional model of the transpeptidase domain of Streptococcus pneumoniae PBP2x starting from its Calpha-atom coordinates

A new method is described for generating all-atom protein structures from Calpha-atom information. The method, which combines both local structural trace alignments and comparative side chain modeling with ab initio side chain modeling, makes use of both the virtual-bond and the dipole-path methods. Provided that 3D structures of structurally and functionally related proteins exist, the method presented here is highly suitable for generating all-atom coordinates of partly solved, low-resolution crystal structures. Particularly the active site region can be modeled accurately with this procedure, which enables investigation of the binding modes of different classes of ligands with molecular dynamics simulations. The method is applied to the trace of Streptococcus pneumoniae, in order to construct an all-atom structure of the transpeptidase domain. Since after generation of full coordinates of the transpeptidase domain the structure had been solved to 2.4 A resolution, new X-ray coordinates for the worst modeled loop (residues T370 to M386; 17 out of a total number of 351 residues constituting the transpeptidase domain) were incorporated, as kindly provided by Dr. Dideberg. The structure was relaxed with molecular dynamics simulations and simulated annealing methods. The RMS deviation between the 144 aligned Calpha-atoms and the corresponding ones in the originally solved 3.5 A resolution crystal structure was 0.98. The 351 Calpha-atoms of the whole transpeptidase domain of the final model showed an RMS deviation of 1.58. The Ramachandran plot showed that 79.3% of the residues are in the most favored regions, with only 1.0% occurring in disallowed regions. The model presented here can be used to investigate the three-dimensional influences of mutations around the active site of PBP2x.