Post-stroke depression: mechanisms, translation and therapy

The interaction between depression and stroke is highly complex. Post-stroke depression (PSD) is among the most frequent neuropsychiatric consequences of stroke. Depression also negatively impacts stroke outcome with increased morbidity, mortality and poorer functional recovery. Antidepressants such as the commonly prescribed selective serotonin reuptake inhibitors improve stroke outcome, an effect that may extend far beyond depression, e.g., to motor recovery. The main biological theory of PSD is the amine hypothesis. Conceivably, ischaemic lesions interrupt the projections ascending from midbrain and brainstem, leading to a decreased bioavailability of the biogenic amines--serotonin (5HT), dopamine (DA) and norepinephrine (NE). Acetylcholine would also be involved. So far, preclinical and translational research on PSD is largely lacking. The implementation and characterization of suitable animal models is clearly a major prerequisite for deeper insights into the biological basis of post-stroke mood disturbances. Equally importantly, experimental models may also pave the way for the discovery of novel therapeutic targets. If we cannot prevent stroke, we shall try to limit its long-term consequences. This review therefore presents animal models of PSD and summarizes potential underlying mechanisms including genomic signatures, neurotransmitter and neurotrophin signalling, hippocampal neurogenesis, cellular plasticity in the ischaemic lesion, secondary degenerative changes, activation of the hypothalamo-pituitary-adrenal (HPA) axis and neuroinflammation. As stroke is a disease of the elderly, great clinical benefit may especially accrue from deciphering and targeting basic mechanisms underlying PSD in aged animals.     

Threading with chemostructural restrictions method for predicting fold and functionally significant residues: application to dipeptidylpeptidase IV (DPP-IV)

We present a new method for more accurate modeling of protein structure, called threading with chemostructural restrictions. This method addresses those cases in which a target sequence has only remote homologues of known structure for which sequence comparison methods cannot provide accurate alignments. Although remote homologues cannot provide an accurate model for the whole chain, they can be used in constructing practically useful models for the most conserved-and often the most interesting-part of the structure. For many proteins of interest, one can suggest certain chemostructural patterns for the native structure based on the available information on the structural superfamily of the protein, the type of activity, the sequence location of the functionally significant residues, and other factors. We use such patterns to restrict (1) a number of possible templates, and (2) a number of allowed chain conformations on a template. The latter restrictions are imposed in the form of additional template potentials (including terms acting as sequence anchors) that act on certain residues. This approach is tested on remote homologues of alpha/beta-hydrolases that have significant structural similarity in the positions of their catalytic triads. The study shows that, in spite of significant deviations between the model and the native structures, the surroundings of the catalytic triad (positions of C(alpha) atoms of 20-30 nearby residues) can be reproduced with accuracy of 2-3 A. We then apply the approach to predict the structure of dipeptidylpeptidase IV (DPP-IV). Using experimentally available data identifying the catalytic triad residues of DPP-IV (David et al., J Biol Chem 1993;268:17247-17252); we predict a model structure of the catalytic domain of DPP-IV based on the 3D fold of prolyl oligopeptidase (Fulop et al., Cell 1998;94:161-170) and use this structure for modeling the interaction of DPP-IV with inhibitor.     

Calculation of protein ionization equilibria with conformational sampling: pK(a) of a model leucine zipper, GCN4 and barnase.

The use of conformational ensembles provided by nuclear magnetic resonance (NMR) experiments or generated by molecular dynamics (MD) simulations has been regarded as a useful approach to account for protein motions in the context of pK(a) calculations, yet the idea has been tested occasionally. This is the first report of systematic comparison of pK(a) estimates computed from long multiple MD simulations and NMR ensembles. As model systems, a synthetic leucine zipper, the naturally occurring coiled coil GCN4, and barnase were used. A variety of conformational averaging and titration curve-averaging techniques, or combination thereof, was adopted and/or modified to investigate the effect of extensive global conformational sampling on the accuracy of pK(a) calculations. Clustering of coordinates is proposed as an approach to reduce the vast diversity of MD ensembles to a few structures representative of the average electrostatic properties of the system in solution. Remarkable improvement of the accuracy of pK(a) predictions was achieved by the use of multiple MD simulations. By using multiple trajectories the absolute error in pK(a) predictions for the model leucine zipper was reduced to as low as approximately 0.25 pK(a) units. The validity, advantages, and limitations of explicit conformational sampling by MD, compared with the use of an average structure and a high internal protein dielectric value as means to improve the accuracy of pK(a) calculations, are discussed.     

Identification of a novel class in the alpha/beta hydrolase fold superfamily: the N-myc differentiation-related proteins

The alpha/beta hydrolases constitute a large protein superfamily that mainly consists of enzymes that catalyze a diverse range of reactions. These proteins exhibit the alpha/beta hydrolase fold, the essential features of which have recently been delineated: the presence of at least five parallel beta-strands, a catalytic triad in a specific order (nucleophile-acid-histidine), and a nucleophilic elbow. Because of the difficulties experimentally in identifying protein structures, we have used a Bayesian computational algorithm (PROBE) to identify the members of this superfamily based on distant sequence relationships. We found that the presence of five sequence motifs, which contain residues important for substrate binding and stabilization of the fold, are required for membership in this superfamily. The superfamily consists of at least 909 members, including the N-myc downstream regulated proteins, which are believed to be involved in cell differentiation. Unlike most of the other superfamily members, the N-myc downstream regulated proteins have never been proposed to possess the alpha/beta hydrolase fold and do not appear to be hydrolases.     

First report of the epiphytic benthic dinoflagellates Coolia canariensis and Coolia malayensis in the waters off Jeju Island, Korea: morphology and rDNA sequences

Coolia spp. are epiphytic and benthic dinoflagellates. Herein, we report for the first time, the occurrence of Coolia canariensis and Coolia malayensis in Korean waters. The morphology of the Korean strains of C. canariensis and C. malayensis isolated from the waters off Jeju Island, Korea was similar to that of the original Canary lslands strains and Malaysian strains, respectively. We found several pores and a line of small knobs on the pore plate, and perforations within the large pores of both C. canariensis and C. malayensis. The plates of the Korean strains of C. canariensis and C. malayensis were arranged in a Kofoidian series of Po, 3′, 7′′, 6c, 6s, 5′′′, and 2′′′′, and Po, 3′, 7′′, 7c, 6–7s, 5′′′, and 2′′′′, respectively. When properly aligned, the large subunit (LSU) rDNA sequence of the Korean strain of C. canariensis was identical to that of the Biscayan strains, but it was 2–3% different from the Canary lslands strain VGO0775 and the Australian strain. In addition, the sequences of small subunit (SSU) and/or LSU rDNA from the two Korean strains of C. malayensis were < 1% different from the Malaysian strains of C. malayensis and the Florida strain CCMP1345 and New Zealand strain CAWD39 (“Coolia monotis”). In phylogenetic trees based on LSU rDNA sequences, the Korean strains of C. malayensis belonged to a clade including the Malaysian strains and these two strains. Therefore, based on genealogical analyses, we suggest that the Korean strain of C. canariensis is closely related to two Atlantic strains and the Australian strain, whereas the Korean strains of C. malayensis are related to the Malaysian strains of C. malayensis and the Florida and New Zealand strains.     

Protein--protein docking of electron transfer complexes: cytochrome c oxidase and cytochrome c

Electron transferring protein complexes form only transiently and the crystal structures of electron transfer protein--protein complexes involving cytochrome c could so far be determined only for the pairs of yeast cytochrome c peroxidase (CcP) with iso-1-cytochrome c (iso-1-cyt c) and with horse heart cytochrome c (cyt c). This article presents models from computational docking for complexes of cytochrome c oxidase (COX) from Paracoccus denitrificans with horse heart cytochrome c, and with its physiological counterpart cytochrome c552 (c552). Initial docking is performed with the FTDOCK program, which permits an exhaustive search of translational and rotational space. A filtering procedure is then applied to reduce the number of complexes to a manageable number. In a final step of structural and energetic refinement, the complexes are optimized by rigid-body energy minimization with the molecular mechanics package CHARMM. This methodology was first tested on the CcP:iso-1-cyt c complex, in which the complex with the lowest CHARMM energy has an RMSD from the crystal structure of only 1.8 A (C(alpha) carbon atoms). Notably, the crystal conformation has an even lower energy. The same procedure was then applied to COX:cyt c and COX:c552. The lowest-energy COX:cyt c complex is very similar to a docking model previously described for the complex of bovine cytochrome c oxidase with horse heart cytochrome c. For the COX:c552 complex, cytochrome c552 is found in two different orientations, depending on whether it is docked against COX from a two-subunit or from a four-subunit crystal structure, respectively. Both conformations are discussed critically in the light of the available experimental data.     

Prediction of protein secondary structures from conformational biases

We use LINUS (the "Local Independently Nucleated Units of Structure"), a procedure developed by Srinivasan and Rose, to provide a physical interpretation of and predict the secondary structures of proteins. The secondary structure type at a given site is identified by the largest conformational bias during short simulations. We examine the rate of successful prediction as a function of temperature and the interaction window. At high temperatures, there is a large propensity for the establishment of beta-strands whereas alpha-helices appear only when the temperature is lower than a certain threshold value. It is found that there exists an optimal temperature at which the correct secondary structures are predicted most accurately. We find that this temperature is close to the peak temperature of the specific heat. Changing the interaction window or carrying out longer simulations approaching equilibrium lead to little change in the optimal success rate. Our findings are in accord with the observation by Srinivasan and Rose that the secondary structures are mainly determined by local interactions and appear in the early stage of folding.     

Establishment and characterization of Roberts syndrome and SC phocomelia model medaka (Oryzias latipes)

Roberts syndrome and SC phocomelia (RBS/SC) are genetic autosomal recessive syndromes caused by establishment of cohesion 1 homolog 2 ( ESCO 2) mutation. RBS/SC appear to have a variety of clinical features, even with the same mutation of the ESCO2 gene. Here, we established and genetically characterized a medaka model of RBS/SC by reverse genetics. The RBS/SC model was screened from a mutant medaka library produced by the Targeting Induced Local Lesions in Genomes method. The medaka mutant carrying the homozygous mutation at R80S in the conserved region of ESCO2 exhibited clinical variety (i.e. developmental arrest with craniofacial and chromosomal abnormalities and embryonic lethality) as characterized in RBS/SC. Moreover, widespread apoptosis and downregulation of some gene expression, including notch1a, were detected in the R80S mutant. The R80S mutant is the animal model for RBS/SC and a valuable resource that provides the opportunity to extend knowledge of ESCO2. Downregulation of some gene expression in the R80S mutant is an important clue explaining non-correlation between genotype and phenotype in RBS/SC.     

Combining mutations in HIV-1 protease to understand mechanisms of resistance

HIV-1 develops resistance to protease inhibitors predominantly by selecting mutations in the protease gene. Studies of resistant mutants of HIV-1 protease with single amino acid substitutions have shown a range of independent effects on specificity, inhibition, and stability. Four double mutants, K45I/L90M, K45I/V82S, D30N/V82S, and N88D/L90M were selected for analysis on the basis of observations of increased or decreased stability or enzymatic activity for the respective single mutants. The double mutants were assayed for catalysis, inhibition, and stability. Crystal structures were analyzed for the double mutants at resolutions of 2.2-1.2 A to determine the associated molecular changes. Sequence-dependent changes in protease-inhibitor interactions were observed in the crystal structures. Mutations D30N, K45I, and V82S showed altered interactions with inhibitor residues at P2/P2', P3/P3'/P4/P4', and P1/P1', respectively. One of the conformations of Met90 in K45I/L90M has an unfavorably close contact with the carbonyl oxygen of Asp25, as observed previously in the L90M single mutant. The observed catalytic efficiency and inhibition for the double mutants depended on the specific substrate or inhibitor. In particular, large variation in cleavage of p6(pol)-PR substrate was observed, which is likely to result in defects in the maturation of the protease from the Gag-Pol precursor and hence viral replication. Three of the double mutants showed values for stability that were intermediate between the values observed for the respective single mutants. D30N/V82S mutant showed lower stability than either of the two individual mutations, which is possibly due to concerted changes in the central P2-P2' and S2-S2' sites. The complex effects of combining mutations are discussed.     

Hybridization and population structure of the Culex pipiens complex in the islands of Macaronesia

The Culex pipiens complex includes two widespread mosquito vector species, Cx. pipiens and Cx. quinquefasciatus. The distribution of these species varies in latitude, with the former being present in temperate regions and the latter in tropical and subtropical regions. However, their distribution range overlaps in certain areas and interspecific hybridization has been documented. Genetic introgression between these species may have epidemiological repercussions for West Nile virus (WNV) transmission. Bayesian clustering analysis based on multilocus genotypes of 12 microsatellites was used to determine levels of hybridization between these two species in Macaronesian islands, the only contact zone described in West Africa. The distribution of the two species reflects both the islands' biogeography and historical aspects of human colonization. Madeira Island displayed a homogenous population of Cx. pipiens, whereas Cape Verde showed a more intriguing scenario with extensive hybridization. In the islands of Brava and Santiago, only Cx. quinquefasciatus was found, while in Fogo and Maio high hybrid rates (～40%) between the two species were detected. Within the admixed populations, second-generation hybrids (～50%) were identified suggesting a lack of isolation mechanisms. The observed levels of hybridization may locally potentiate the transmission to humans of zoonotic arboviruses such as WNV.     

High-resolution electron microscopical study of cyst walls of Entamoeba spp

Knowledge of the fine structural organization, molecular composition and permeability properties of the cell surface of intestinal protozoan cysts is important to understand the biologic basis of their resistance. Recent studies on the biology of the cyst walls of Entamoeba histolytica and Entamoeba invadens have considerably advanced knowledge on the cellular processes involved in the transport and surface deposition of the main cyst wall components. Using transmission electron microscopy, cytochemistry, scanning electron microscopy and freeze-fracture techniques, we have obtained new information. In mature cysts the permeability of Entamoeba cysts is limited to small molecules not by the cyst wall, but by the plasma membrane, as demonstrated with the use of ruthenium red as an electron-dense tracer. Cell walls of E. histolytica cysts are made up of five to seven layers of unordered fibrils 7-8 nm thick. Alcian blue stains a regular mesh of fibrils approximately 4 nm thick, running perpendicularly to the cyst wall. In addition, abundant ionogenic groups are seen in cyst walls treated with cationized ferritin. In the mature cysts of E. histolytica and E. invadens small cytoplasmic vesicles with granular material were in close contact with the plasma membrane, suggesting a process of fusion and deposition of granular material to the cell wall. The plasma membrane of mature cysts is devoid of intramembrane particles when analyzed with the freeze-fracture technique. When viewed with scanning electron microscopy the surface of E. histolytica cysts clearly differs from that of Entamoeba coli and E. invadens.     

Molecular dynamics study of a hyperthermophilic and a mesophilic rubredoxin

In recent years, increased interest in the origin of protein thermal stability has gained attention both for its possible role in understanding the forces governing the folding of a protein and for the design of new highly stable engineered biocatalysts. To study the origin of thermostability, we have performed molecular dynamics simulations of two rubredoxins, from the mesophile Clostridium pasteurianum and from the hyperthermophile Pyrococcus furiosus. The simulations were carried out at two temperatures, 300 and 373 K, for each molecule. The length of the simulations was within the range of 6-7.2 ns. The rubredoxin from the hyperthermophilic organism was more flexible than its mesophilic counterpart at both temperatures; however, the overall flexibility of both molecules at their optimal growth temperature was the same, despite 59% sequence homology. The conformational space sampled by both molecules was larger at 300 K than at 373 K. The essential dynamics analysis showed that the principal overall motions of the two molecules are significantly different. On the contrary, each molecule showed similar directions of motion at both temperatures.     

The sterol-binding protein Kes1/Osh4p is a regulator of polarized exocytosis

Oxysterol-binding protein (OSBP)-related protein Kes1/ Osh4p is implicated in nonvesicular sterol transfer between membranes in Saccharomyces cerevisiae. However, we found that Osh4p associated with exocytic vesicles that move from the mother cell into the bud, where Osh4p facilitated vesicle docking by the exocyst tethering complex at sites of polarized growth on the plasma membrane. Osh4p formed complexes with the small GTPases Cdc42p, Rho1p and Sec4p, and the exocyst complex subunit Sec6p, which was also required for Osh4p association with vesicles. Although Osh4p directly affected polarized exocytosis, its role in sterol trafficking was less clear. Contrary to what is predicted for a sterol-transfer protein, inhibition of sterol binding by the Osh4p Y97F mutation did not cause its inactivation. Rather, OSH4(Y97F) is a gain-of-function mutation that causes dominant lethality. We propose that in response to sterol binding and release Osh4p promotes efficient exocytosis through the co-ordinate regulation of Sac1p, a phosphoinositide 4-phosphate (PI4P) phosphatase, and the exocyst complex. These results support a model in which Osh4p acts as a sterol-dependent regulator of polarized vesicle transport, as opposed to being a sterol-transfer protein.