Template-based modeling and ab-initio docking using CoDock in CAPRI

Integration of template-based modeling, global sampling and precise scoring is crucial for the development of molecular docking programs with improved accuracy. We combined template-based modeling and ab-initio docking protocol as hybrid docking strategy called CoDock for the docking and scoring experiments of the seventh CAPRI edition. For CAPRI rounds 38-45, we obtained acceptable or better models in the top 10 submissions for eight out of the 16 evaluated targets as predictors, nine out of the 16 targets as scorers. Especially, we submitted acceptable models for all of the evaluated protein-oligosaccharide targets. For the CASP13-CAPRI experiment (round 46), we obtained acceptable or better models in the top 5 submissions for 10 out of the 20 evaluated targets as predictors, 11 out of the 20 targets as scorers. The failed cases for our group were mainly the difficult targets and the protein-peptide systems in CAPRI and CASP13-CAPRI experiments. In summary, this CAPRI edition showed that our hybrid docking strategy can be efficiently adapted to the increasing variety of challenges in the field of molecular interactions.     

C1188D mutation abolishes specific recognition between MLL1-CXXC domain and CpG site by inducing conformational switch of flexible N-terminal

Mixed lineage leukemia protein (MLL1 protein) recognizes the CpG site via its CXXC domain and is frequently associated with leukemia. The specific recognition is abolished by C1188D mutation, which also prevents MLL-related leukemia. In this paper, multiple molecular dynamic (MD) simulations were performed to investigate the mechanism of recognition and influences of C1188D mutation. Started from fully dissociated DNA and MLL1-CXXC domain, remarkably, the center of mass (COM) of MLL1-CXXC domain quickly concentrates on the vicinity of the CpG site in all 53 short MD simulations. Extended simulations of the wild type showed that the native complex formed in 500 ns among 4 of 53 simulations. In contrast, the C1188D mutant COM distributed broadly around the DNA and the native complex was not observed in any of the extended simulations. Simulations on the apo MLL1-CXXC domain further suggest that the wild type protein remained predominantly in an open form that closely resembles its structure in the native complex whereas C1188D mutant formed predominantly compact structures in which the N- terminal bends to D1188. This conformational switch hinders the formation of encounter complex, thus abolishes the recognition. Our study also provides clues to the study mechanism of recognition, by the CXXC domain from proteins like DNA methyltransferase and ten-eleven translocation enzymes.     

Archiascomycetes: detection of a major new lineage within the Ascomycota

For phylogenetic analysis of the higher fungi, we sequenced the nuclear small subunit rRNA (18S rRNA) gene fromTaphrina populina, the type species of the genusTaphrina, andProtomyces lactucae-debilis. The molecular phylogeny inferred from these 2 sequences and 75 sequences from the DNA data bank divided the Ascomycota into three major lineages: the hemiascomycetes, the euascomycetes, and the archiascomycetes, newly described herein. The former two lineages are monophyletic, whereas the archiascomycetes, which originated first and are comprised ofTaphrina, Protomyces, Saitoella, Schizosaccharomyces, andPneumocystis, may not be monophyletic. Among the archiascomycetes, theTaphrina/Protomyces branch is monophyletic. Confirmation of the archiascomycetes as a monophyletic taxonomic class will require comparison of additional genetically defined characters.This work was supported in part by grants 05454030 from the Ministry of Education, Science, and Culture of Japan (to J. S.) and 4369 from the Japan Society for the Promotion of Science Fellowship Programs (to H. N.).     

Structural dynamics of the mitochondrial ADP/ATP carrier revealed by molecular dynamics simulation studies

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier has been recently crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). In the crystal structure, the six-transmembrane helix bundle that defines the nucleotide translocation pathway is closed on the matrix side due to sharp kinks in the odd-numbered helices. The closed conformation is further sealed by the loops protruding into the matrix that interact through an intricate network of charge-pairs. To gain insight into its structural dynamics we performed molecular dynamics (MD) simulation studies of the ADP/ATP carrier with and without its cocrystallized inhibitor. The two trajectories sampled a conformational space around two different configurations characterized by distinct salt-bridge networks with a significant shift from inter- to intrarepeat bonding on the matrix side in the absence of CATR. Analysis of the geometrical parameters defining the transmembrane helices showed that even-numbered helices can undergo a face rotation, whereas odd-numbered helices can undergo a change in the wobble angle with a conserved proline acting as molecular hinge. Our results provide new information on the dynamical properties of the ADP/ATP carrier and for the first time yield a detailed picture of a stable carrier conformation in absence of the inhibitor.     

Radiation damage in the bulk and at the surface

The UK environmental e-science initiative supports the development and modification of simulation tools used to study radiation damage effects. We discuss the development and modification to the DL_POLY molecular dynamics (MD) code. Using the newly developed tools, we study the effects of radiation damage related to the safe encapsulation of highly radioactive materials, including nuclear waste. We address the possible differences between the radiation damage in the bulk and at the surface of a material, and perform MD simulations of energetic events in zircon structure. We find that in the case of readily amorphizable material, the formation of a stable alternative covalent network reduces the possible effect of the surface on the damaged structure.     

厦门港纤小裸甲藻（Takayama pulchellum）的形态、生长及分子特征

在2003年厦门港发生的一次裸甲藻赤潮中分离得到了纤小裸甲藻（Takayama pulchellum）,并对其形态学进行了光镜和电镜的观察,该株藻细胞个体较小,长约20～23μm,宽约14～20μm.藻细胞外表卵圆形,细胞上壳有一条明显的S形顶沟,顶沟是特异性的反S形.与模式种相比,这一株裸甲藻细胞核更靠上部,大部分细胞横沟都没有偏移.纤小裸甲藻在盐度为28时比生长率最大,达到0.42,随着盐度下降,比生长率也跟着下降,当盐度下降到16时,生长率为0.盐度范围在8～16时,超过80%的藻类能存活48h;当盐度低于4时,生长率和存活率都为0.24～27℃是纤小裸甲藻的最适生长温度,比生长率超过0.50,当温度升到30℃时,生长率急剧下降到约0.35.毒素测定结果显示该株裸甲藻不含麻痹性贝毒和神经性贝毒.本株纤小裸甲藻大亚基D1～D2区序列长度为721bp,与基因库中该种的一株相似种同源性超过99%.对18株裸甲藻转录间隔区（ITS）序列建立的系统进化树显示,纤小裸甲藻和Karlodinium micrum的距离最近,通过ITS序列建立的系统发育树大致能够把Akashiwo属、Karenia属、Karlodinium属和Takayama属与Gymnodinium属区分开来.     

Response of SCP-2L domain of human MFE-2 to ligand removal: binding site closure and burial of peroxisomal targeting signal

In the study of the structure and function relationship of human MFE-2, we have investigated the dynamics of human MFE-2SCP-2L (hSCP-2L) and its response to ligand removal. A comparison was made with homologous rabbit SCP-2. Breathing and a closing motion are found, identifiable with an adjustment in size and a closing off of the binding pocket. Crucial residues for structural integrity have been identified. Particularly mobile areas of the protein are loop 1 that is connecting helices A and C in space, and helix D, next to the entrance of the pocket. In hSCP-2L, the binding pocket gets occupied by Phe93, which is making a tight hydrophobic contact with Trp36. In addition, it is found that the C-terminal peroxisomal targeting signal (PTS1) that is solvent exposed in the complexed structure becomes buried when no ligand is present. Moreover, an anti-correlation exists between burial of PTS1 and the size of the binding pocket. The results are in accordance with plant nsLTPs, where a similar accommodation of binding pocket size was found after ligand binding/removal. Furthermore, the calculations support the suggestion of a ligand-assisted targeting mechanism.     

Biochemical evidence accumulates across neurons to drive a network-level eruption

1. Download : Download high-res image (214KB)
2. Download : Download full-size image

     

An alternate sucrose binding mode in the E203Q Arabidopsis invertase mutant: an X-ray crystallography and docking study

In the present study, we report on the X-ray crystallographic structure of a GH32 invertase mutant, (i.e., the Arabidopsis thaliana cell-wall invertase 1-E203Q, AtcwINV1-mutant) in complex with sucrose. This structure was solved to reveal the features of sugar binding in the catalytic pocket. However, as demonstrated by the X-ray structure the sugar binding and the catalytic pocket arrangement is significantly altered as compared with what was expected based on previous X-ray structures on GH-J clan enzymes. We performed a series of docking and molecular dynamics simulations on various derivatives of AtcwINV1 to reveal the reasons behind this modified sugar binding. Our results demonstrate that the E203Q mutation introduced into the catalytic pocket triggers conformational changes that alter the wild type substrate binding. In addition, this study also reveals the putative productive sucrose binding modus in the wild type enzyme.     

New pyrazolylpyrazoline derivatives as dual acting antimalarial-antileishamanial agents: synthesis,biological evaluation and molecular modelling simulations

Promising inhibitory activities of the parasite multiplication were obtained upon evaluation of in vivo antimalarial activities of new pyrazolylpyrazoline derivatives against Plasmodium berghei infected mice. Further evaluation of 5b and 6a against chloroquine-resistant strain (RKL9) of P. falciparum showed higher potency than chloroquine. In vitro antileishmanial activity testing against Leishmania aethiopica promastigote and amastigote forms indicated that 5b, 6a and 7b possessed promising activity compared to miltefosine and amphotericin B deoxycholate. Moreover, antileishmanial activity reversal of the active compounds via folic and folinic acids showed comparable results to the positive control trimethoprim, indicating an antifolate mechanism via targeting leishmanial DHFR and PTR1. The compounds were non-toxic at 125, 250 and 500 mg/kg. In addition, docking of the most active compound against putative malarial target Pf-DHFR-TS and leishmanial PTR1 rationalised the observed activities. Molecular dynamics simulations confirmed a stable and high potential binding of 7a against leishmanial PTR1.