A mutation in VPS35, encoding a subunit of the retromer complex, causes late-onset Parkinson disease.

To identify rare causal variants in late-onset Parkinson disease (PD), we investigated an Austrian family with 16 affected individuals by exome sequencing. We found a missense mutation, c.1858G>A (p.Asp620Asn), in the VPS35 gene in all seven affected family members who are alive. By screening additional PD cases, we saw the same variant cosegregating with the disease in an autosomal-dominant mode with high but incomplete penetrance in two further families with five and ten affected members, respectively. The mean age of onset in the affected individuals was 53 years. Genotyping showed that the shared haplotype extends across 65 kilobases around VPS35. Screening the entire VPS35 coding sequence in an additional 860 cases and 1014 controls revealed six further nonsynonymous missense variants. Three were only present in cases, two were only present in controls, and one was present in cases and controls. The familial mutation p.Asp620Asn and a further variant, c.1570C>T (p.Arg524Trp), detected in a sporadic PD case were predicted to be damaging by sequence-based and molecular-dynamics analyses. VPS35 is a component of the retromer complex and mediates retrograde transport between endosomes and the trans-Golgi network, and it has recently been found to be involved in Alzheimer disease.

A leucine-rich repeat assembly approach for homology modeling of the human TLR5-10 and mouse TLR11-13 ectodomains.

So far, 13 groups of mammalian Toll-like receptors (TLRs) have been identified. Most TLRs have been shown to recognize pathogen-associated molecular patterns from a wide range of invading agents and initiate both innate and adaptive immune responses. The TLR ectodomains are composed of varying numbers and types of leucine-rich repeats (LRRs). As the crystal structures are currently missing for most TLR ligand-binding ectodomains, homology modeling enables first predictions of their three-dimensional structures on the basis of the determined crystal structures of TLR ectodomains. However, the quality of the predicted models that are generated from full-length templates can be limited due to low sequence identity between the target and templates. To obtain better templates for modeling, we have developed an LRR template assembly approach. Individual LRR templates that are locally optimal for the target sequence are assembled into multiple templates. This method was validated through the comparison of a predicted model with the crystal structure of mouse TLR3. With this method, we also constructed ectodomain models of human TLR5, TLR6, TLR7, TLR8, TLR9, and TLR10 and mouse TLR11, TLR12, and TLR13 that can be used as first passes for a computational simulation of ligand docking or to design mutation experiments. This template assembly approach can be extended to other repetitive proteins.

TollML: a database of toll-like receptor structural motifs.

Toll-like receptors (TLRs) play a key role in the innate immune system. TLRs recognize pathogen-associated molecular patterns and initiate an intracellular kinase cascade to induce an immediate defensive response. During recent years TLRs have become the focus of tremendous research interest. A central repository for the growing amount of relevant TLR sequence information has been created. Nevertheless, structural motifs of most sequenced TLR proteins, such as leucine-rich repeats (LRRs), are poorly annotated in the established databases. A database that organizes the structural motifs of TLRs could be useful for developing pattern recognition programs, structural modeling and understanding functional mechanisms of TLRs. We describe TollML, a database that integrates all of the TLR sequencing data from the NCBI protein database. Entries were first divided into TLR families (TLR1-23) and then semi-automatically subdivided into three levels of structural motif categories: (1) signal peptide (SP), ectodomain (ECD), transmembrane domain (TD) and Toll/IL-1 receptor (TIR) domain of each TLR; (2) LRRs of each ECD; (3) highly conserved segment (HCS), variable segment (VS) and insertions of each LRR. These categories can be searched quickly using an easy-to-use web interface and dynamically displayed by graphics. Additionally, all entries have hyperlinks to various sources including NCBI, Swiss-Prot, PDB, LRRML and PubMed in order to provide broad external information for users. The TollML database is available at http://tollml.lrz.de.

Unraveling the molecular mechanisms of nitrogenase conformational protection against oxygen in diazotrophic bacteria.

BACKGROUND: G. diazotrophicus and A. vinelandii are aerobic nitrogen-fixing bacteria. Although oxygen is essential for the survival of these organisms, it irreversibly inhibits nitrogenase, the complex responsible for nitrogen fixation. Both microorganisms deal with this paradox through compensatory mechanisms. In A. vinelandii a conformational protection mechanism occurs through the interaction between the nitrogenase complex and the FeSII protein. Previous studies suggested the existence of a similar system in G. diazotrophicus, but the putative protein involved was not yet described. This study intends to identify the protein coding gene in the recently sequenced genome of G. diazotrophicus and also provide detailed structural information of nitrogenase conformational protection in both organisms. RESULTS: Genomic analysis of G. diazotrophicus sequences revealed a protein coding ORF (Gdia0615) enclosing a conserved "fer2" domain, typical of the ferredoxin family and found in A. vinelandii FeSII. Comparative models of both FeSII and Gdia0615 disclosed a conserved beta-grasp fold. Cysteine residues that coordinate the 2[Fe-S] cluster are in conserved positions towards the metallocluster. Analysis of solvent accessible residues and electrostatic surfaces unveiled an hydrophobic dimerization interface. Dimers assembled by molecular docking presented a stable behaviour and a proper accommodation of regions possibly involved in binding of FeSII to nitrogenase throughout molecular dynamics simulations in aqueous solution. Molecular modeling of the nitrogenase complex of G. diazotrophicus was performed and models were compared to the crystal structure of A. vinelandii nitrogenase. Docking experiments of FeSII and Gdia0615 with its corresponding nitrogenase complex pointed out in both systems a putative binding site presenting shape and charge complementarities at the Fe-protein/MoFe-protein complex interface. CONCLUSIONS: The identification of the putative FeSII coding gene in G. diazotrophicus genome represents a large step towards the understanding of the conformational protection mechanism of nitrogenase against oxygen. In addition, this is the first study regarding the structural complementarities of FeSII-nitrogenase interactions in diazotrophic bacteria. The combination of bioinformatic tools for genome analysis, comparative protein modeling, docking calculations and molecular dynamics provided a powerful strategy for the elucidation of molecular mechanisms and structural features of FeSII-nitrogenase interaction.

Inhibition of Toll-like receptors TLR4 and 7 signaling pathways by SIGIRR: a computational approach.

Toll-like receptors (TLRs) belong to the Toll-like receptor/interleukin-1 receptor (TLR/IL-1R) superfamily which is defined by a common cytoplasmic Toll/interleukin-1 receptor (TIR) domain. TLRs recognize pathogen-associated molecular patterns and initiate an intracellular kinase cascade to trigger an immediate defensive response. SIGIRR (single immunoglobulin interleukin-1 receptor-related molecule), another member of the TLR/IL-1R superfamily, acts as a negative regulator of MyD88-dependent TLR signaling. It attenuates the recruitment of MyD88 adaptors to the receptors with its intracellular TIR domain. Thus, SIGIRR is a highly important molecule for the therapy of autoimmune diseases caused by TLRs. So far, the structural mechanism of interactions between SIGIRR, TLRs and adaptor molecules is unclear. To develop a working hypothesis for this interaction, we constructed three-dimensional models for the TIR domains of TLR4, TLR7, MyD88 and SIGIRR based on computational modeling. Through protein-protein docking analysis, we developed models of essential complexes involved in the TLR4 and 7 signaling and the SIGIRR inhibiting processes. We suggest that SIGIRR may exert its inhibitory effect through blocking the molecular interface of TLR4, TLR7 and the MyD88 adaptor mainly via its BB-loop region.

First-principles simulation of photoreactions in biological systems.

First-principles simulations start to be applicable to the photochemistry and photophysics in biological systems. In this review the prerequisites for investigating such excited state phenomena in large systems are outlined. Generally, a quantum mechanical description of the electronic structure is combined with molecular dynamics simulations, which allows to describe the motion of the atoms in the field produced by the quantum-mechanical potential. Like this, bonds can be formed and broken, that is, chemical reactions can be simulated. The review focuses on applications of first-principles molecular dynamics to photoactive proteins.

Homology modeling of human Toll-like receptors TLR7, 8, and 9 ligand-binding domains.

Toll-like receptors (TLRs) play a key role in the innate immune system. The TLR7, 8, and 9 compose a family of intracellularly localized TLRs that signal in response to pathogen-derived nucleic acids. So far, there are no crystallographic structures for TLR7, 8, and 9. For this reason, their ligand-binding mechanisms are poorly understood. To enable first predictions of the receptor-ligand interaction sites, we developed three-dimensional structures for the leucine-rich repeat ectodomains of human TLR7, 8, and 9 based on homology modeling. To achieve a high sequence similarity between targets and templates, structural segments from all known TLR ectodomain structures (human TLR1/2/3/4 and mouse TLR3/4) were used as candidate templates for the modeling. The resulting models support previously reported essential ligand-binding residues. They also provide a basis to identify three potential receptor dimerization mechanisms. Additionally, potential ligand-binding residues are identified using combined procedures. We suggest further investigations of these residues through mutation experiments. Our modeling approach can be extended to other members of the TLR family or other repetitive proteins.

LRRML: a conformational database and an XML description of leucine-rich repeats (LRRs).

BACKGROUND: Leucine-rich repeats (LRRs) are present in more than 6000 proteins. They are found in organisms ranging from viruses to eukaryotes and play an important role in protein-ligand interactions. To date, more than one hundred crystal structures of LRR containing proteins have been determined. This knowledge has increased our ability to use the crystal structures as templates to model LRR proteins with unknown structures. Since the individual three-dimensional LRR structures are not directly available from the established databases and since there are only a few detailed annotations for them, a conformational LRR database useful for homology modeling of LRR proteins is desirable. DESCRIPTION: We developed LRRML, a conformational database and an extensible markup language (XML) description of LRRs. The release 0.2 contains 1261 individual LRR structures, which were identified from 112 PDB structures and annotated manually. An XML structure was defined to exchange and store the LRRs. LRRML provides a source for homology modeling and structural analysis of LRR proteins. In order to demonstrate the capabilities of the database we modeled the mouse Toll-like receptor 3 (TLR3) by multiple templates homology modeling and compared the result with the crystal structure. CONCLUSION: LRRML is an information source for investigators involved in both theoretical and applied research on LRR proteins. It is available at http://zeus.krist.geo.uni-muenchen.de/~lrrml.

Gene characterization and predicted protein structure of the mitochondrial chaperonin HSP10 of Trypanosoma cruzi.

Heat shock protein (HSP) 10 is a member of the highly conserved group of molecular chaperons, which are necessary for efficient folding of many proteins in normal and stress conditions and have been implicated in several human diseases. We have characterized the HSP10 genes of Trypanosoma cruzi, the causative agent of Chagas' disease. After sequence analysis of clones obtained from the T. cruzi Genome Initiative, we show that the T. cruzi HSP10 coding region is 300 bp long, encoding a polypeptide of 100 amino acids with highest sequence identity (83%) to HSP10 of Trypanosoma brucei and lowest (28%) to HSP10 of Escherichia coli. The T. cruzi HSP10 genes are arranged in 3 tandemly repeated copies, which give rise to a major mRNA of 1.0 kb that remains unaltered during heat shock; a smaller mRNA species is induced at 37 degrees C by alternate polyadenylation. Finally, the presence of a conserved 5-amino acid residue deletion in trypanosomatid HSP10s led us to generate a molecular model of the T. cruzi HSP10 structure. The oligomeric assembly of this model shows some peculiar characteristics that may have functional significance.

Mutational analysis and molecular modeling of the binding of Staphylococcus aureus enterotoxin C2 to a murine T cell receptor Vbeta10 chain.

We investigated the role of the beta chain HV4 region in the binding of a Vbeta10 T cell receptor to superantigen S. aureus enterotoxin C2 (SEC2)/MHC class II complexes. Residues 6971 of the Cw3/1.1 TCR Vbeta10 chain, derived from an H-2K(d)-restricted cytotoxic clone, were individually changed to alanine using site-directed mutagenesis, and mutated TCR beta chains were transfected along with the wild-type TCR alpha chain into a TCR alpha(-)beta(-) T hybridoma. SEC2/MHC recognition was measured by IL-2 production. Alanine substitutions in the HV4beta region, either did not affect (Ser69 and Lys71), or increased the recognition of SEC2/HLA-DR1 complex (Asp70), arguing against a general and direct role for the HV4beta region in superantigenrecognition. A theoretical-computational model of the SEC2/TCR Vbeta10 chain complex was constructed and predicted the presence of a unique salt bridge between Vbeta Asp30 and SEC2 Lys103. A perfect correlation was found between the likely presence of this salt bridge and the capacity of the HV4beta and previously obtained CDR1beta alanine mutants to induce an equal or greater response than the wild-type TCR.