Protein pocket and ligand shape comparison and its application in virtual screening.

Understanding molecular recognition is one major requirement for drug discovery and design. Physicochemical and shape complementarity between two binding partners is the driving force during complex formation. In this study, the impact of shape within this process is analyzed. Protein binding pockets and co-crystallized ligands are represented by normalized principal moments of inertia ratios (NPRs). The corresponding descriptor space is triangular, with its corners occupied by spherical, discoid, and elongated shapes. An analysis of a selected set of sc-PDB complexes suggests that pockets and bound ligands avoid spherical shapes, which are, however, prevalent in small unoccupied pockets. Furthermore, a direct shape comparison confirms previous studies that on average only one third of a pocket is filled by its bound ligand, supplemented by a 50 % subpocket coverage. In this study, we found that shape complementary is expressed by low pairwise shape distances in NPR space, short distances between the centers-of-mass, and small deviations in the angle between the first principal ellipsoid axes. Furthermore, it is assessed how different binding pocket parameters are related to bioactivity and binding efficiency of the co-crystallized ligand. In addition, the performance of different shape and size parameters of pockets and ligands is evaluated in a virtual screening scenario performed on four representative targets.

A Brazilian family with hereditary inclusion body myopathy associated with Paget disease of bone and frontotemporal dementia

Inclusion body myopathy associated with Paget disease and frontotemporal dementia (IBMPFD) is a progressive and usually misdiagnosed autosomal dominant disorder. It is clinically characterized by a triad of features: proximal and distal myopathy, early onset Paget disease of bone (PDB), and frontotemporal dementia (FTD). It is caused by missense mutations in the valosin-containing protein (VCP) gene. We describe here the clinical and molecular findings of the first Brazilian family identified with IBMPFD. Progressive myopathy affecting the limb girdles was detected by clinical examination followed by muscle biopsy and creatine kinase measurement. PDB was suggested after anatomopathological bone examination and FTD was diagnosed by clinical, neuropsychological and language evaluations. Brain magnetic resonance revealed severe atrophy of the anterior temporal lobes, including the hippocampi. A R93C mutation in VCP was detected by direct sequencing screening in subject W (age 62) and in his mother. Four more individuals diagnosed with "dementia" were reported in this family. We also present a comprehensive genotype-phenotype correlation analysis of mutations in VCP in 182 patients from 29 families described in the literature and show that while IBM is a conspicuously penetrant symptom, PDB has a lower penetrance when associated with mutations in the AAAD1 domain and FTD has a lower penetrance when associated with mutations in the Junction (L1-D1) domain. Furthermore, the R93C mutation is likely to be associated with the penetrance of all the clinical symptoms of the triad.

Modeling of allergen proteins found in sea food products

Shellfish are a source of food allergens, and their consumption is the cause of severe allergic reactions in humans. Tropomyosins, a family of muscle proteins, have been identified as the major allergens in shellfish and mollusks species. Nevertheless, few experimentally determined three-dimensional structures are available in the Protein Data Base (PDB). In this study, 3D models of several homologous of tropomyosins present in marine shellfish and mollusk species (Chaf 1, Met e1, Hom a1, Per v1, and Pen a1) were constructed, validated, and their immunoglobulin E binding epitopes were identified using bioinformatics tools. All protein models for these allergens consisted of long alpha-helices. Chaf 1, Met e1, and Hom a1 had six conserved regions with sequence similarities to known epitopes, whereas Per v1 and Pen a1 contained only one. Lipophilic potentials of identified epitopes revealed a high propensity of hydrophobic amino acids in the immunoglobulin E binding site. This information could be useful to design tropomyosin-specific immunotherapy for sea food allergies.

Paragenetic history of the Ordovician Trenton Group carbonates, Southwestern Ontario

Geochemical examination of the rock matrix and cements from core material extracted from four oil wells within southwestern Ontario suggest various stages of diagenetic alteration and preservation of the Trenton Group carbonates. The geochemical compositions of Middle Ordovician (LMC) brachiopods reflect the physicochemical water conditions of the ambient depositional environment. The sediments appear to have been altered in the presence of mixed waters during burial in a relatively open diagenetic microenvironment. Conodont CAl determination suggests that the maturation levels of the Trenton Group carbonates are low and proceeded at temperatures of about 30 - 50°C within the shallow burial environment. The Trenton Group carbonates are characterized by two distinct stages of dolomitization which proceeded at elevated temperatures. Preexisting fracture patterns, and block faulting controlled the initial dolomitization of the precursor carbonate matrix. Dolomitization progressed In the presence of warm fluids (60 75°C) with physicochemical conditions characteristic of a progressively depleted basinal water. The matrix is mostly Idiotopic-S and Idiotopic-E dolomite, with Xenotopic-A dolomite dominating the matrix where fractures occur. The second stage of dolomitization involved hydrothermal basinal fluid(s) with temperatures of about 60 - 70°C. These are the postulated source for the saddle dolomite and blocky calcite cements occurring in pore space and fractures. Rock porosity was partly occluded by Idiotopic-E type dolomite. Late stage saddle dolomite, calcite, anhydrite, pyrite, marcasite and minor sphalerite and celestite cements effectively fill any remaining porosity within specific horizons. Based on cathode luminescence, precipitation of the different diagenetic phases probably proceeded in open diagenetic systems from chemically homogeneous fluids. Ultraviolet fluorescence of 11 the matrix and cements demonstrated that hydrocarbons were present during the earliest formation of saddle dolomite. Oxygen isotope values of -7.6 to -8.5 %0 (PDB), and carbon isotope values of - 0.5 and -3.0 %0 (PDB) from the latest stage dog-tooth calcite cement suggest that meteoric water was introduced into the system during their formation. This is estimated to have occurred at temperatures of about 25 - 40°C. Specific facies associations within the Trenton Group carbonates exhibit good hydrocarbon generating potential based on organic carbon preservation (1-3.5%). Thermal maturation and Lopatin burial-history evaluations suggest that hydrocarbons were generated within the Trenton Group carbonates some time after 300 Ma . Progressively depleted vanadium trends measured from hydrocarbon samples within southwestern Ontario suggests its potential use as a hydrocarbon migration indicator on local (within an oilfield) and on regional scales.

Calmodulin/KCa3.1 channel interactions as determinant to the KCa3.1 Ca2+ dependent gating : theoretical and experimental analyses

Differentes études ont montré que la sensibilité au Ca2+ du canal KCa3.1, un canal potassique indépendant du voltage, était conférée par la protéine calmoduline (CaM) liée de façon constitutive au canal. Cette liaison impliquerait la région C-lobe de la CaM et un domaine de $\ikca$ directement relié au segment transmembranaire S6 du canal. La CaM pourrait égalment se lier au canal de façon Ca2+ dépendante via une interaction entre un domaine de KCa3.1 du C-terminal (CaMBD2) et la région N-lobe de la CaM. Une étude fut entreprise afin de déterminer la nature des résidus responsables de la liaison entre le domaine CaMBD2 de KCa3.1 et la région N-lobe de la CaM et leur rôle dans le processus d'ouverture du canal par le Ca2+. Une structure 3D du complexe KCa3.1/CaM a d'abord été générée par modélisation par homologie avec le logiciel MODELLER en utilisant comme référence la structure cristalline du complexe SK2.2/CaM (PDB: 1G4Y). Le modèle ainsi obtenu de KCa3.1 plus CaM prévoit que le segment L361-S372 dans KCa3.1 devrait être responsable de la liaison dépendante du Ca2+ du canal avec la région N-lobe de la CaM via les résidus L361 et Q364 de KCa3.1 et E45, E47 et D50 de la CaM. Pour tester ce modèle, les résidus dans le segment L361-S372 ont été mutés en Cys et l'action du MTSET+ (chargé positivement) et MTSACE (neutre) a été mesurée sur l'activité du canal. Des enregistrements en patch clamp en configuration ``inside-out`` ont montré que la liaison du réactif chargé MTSET+ au le mutant Q364C entraîne une forte augmentation du courant, un effet non observé avec le MTSACE. De plus les mutations E45A et E47A dans la CaM, ont empêché l'augmentation du courant initié par MTSET+ sur le mutant Q364C. Une analyse en canal unitaire a confirmé que la liaison MTSET+ à Q364C cause une augmentation de la probabilité d'ouverture de KCa3.1 par une déstabilisation de l'état fermé du canal. Nous concluons que nos résultats sont compatibles avec la formation de liaisons ioniques entre les complexes chargés positivement Cys-MTSET+ à la position 364 de KCa3.1 et les résidus chargés négativement E45 et E47 dans la CaM. Ces données confirment qu'une stabilisation électrostatique des interactions CaM/KCa3.1 peut conduire à une augmentation de la probabilité d'ouverture du canal en conditions de concentrations saturantes de Ca2+.

Characterization and Pathogenicity of Vibrio cholerae and Vibrio vulnificus from Marine environments

The genus Vibrioof the family Vibrionaceae are Gram negative, oxidasepositive, rod- or curved- rodshaped facultative anaerobes, widespread in marine and estuarine environments. Vibrio species are opportunistic human pathogens responsible for diarrhoeal disease, gastroenteritis, septicaemia and wound infections and are also pathogens of aquatic organisms, causing infections to crustaceans, bivalves and fishes. In the present study, marine environmental samples like seafood and water and sediment samples from aquafarms and mangroves were screened for the presence of Vibrio species. Of the134 isolates obtained from the various samples, 45 were segregated to the genus Vibrio on the basis of phenotypic characterization.like Gram staining, oxidase test, MoF test and salinity tolerance. Partial 16S rDNA sequence analysis was utilized for species level identification of the isolates and the strains were identified as V. cholerae(N=21), V. vulnificus(N=18), V. parahaemolyticus(N=3), V. alginolyticus (N=2) and V. azureus (N=1). The genetic relatedness and variations among the 45 Vibrio isolates were elucidated based on 16S rDNA sequences. Phenotypic characterization of the isolates was based on their response to 12 biochemical tests namely Voges-Proskauers’s (VP test), arginine dihydrolase , tolerance to 3% NaCl test, ONPG test that detects β-galactosidase activity, and tests for utilization of citrate, ornithine, mannitol, arabinose, sucrose, glucose, salicin and cellobiose. The isolates exhibited diverse biochemical patterns, some specific for the species and others indicative of their environmental source.Antibiogram for the isolates was determined subsequent to testing their susceptibility to 12 antibiotics by the disc diffusion method. Varying degrees of resistance to gentamycin (2.22%), ampicillin(62.22%), nalidixic acid (4.44%), vancomycin (86.66), cefixime (17.77%), rifampicin (20%), tetracycline (42.22%) and chloramphenicol (2.22%) was exhibited. All the isolates were susceptible to streptomycin, co-trimoxazole, trimethoprim and azithromycin. Isolates from all the three marine environments exhibited multiple antibiotic resistance, with high MAR index value. The molecular typing methods such as ERIC PCR and BOX PCR revealed intraspecies relatedness and genetic heterogeneity within the environmental isolatesof V. cholerae and V. vulnificus. The 21 strains of V. choleraewere serogroupedas non O1/ non O139 by screening for the presence O1rfb and O139 rfb marker genes by PCR. The virulence/virulence associated genes namely ctxA, ctxB, ace, VPI, hlyA, ompU, rtxA, toxR, zot, nagst, tcpA, nin and nanwere screened in V. cholerae and V. vulnificusstrains.The V. vulnificusstrains were also screened for three species specific genes viz., cps, vvhand viu. In V. cholerae strains, the virulence associated genes like VPI, hlyA, rtxA, ompU and toxR were confirmed by PCR. All the isolates, except for strain BTOS6, harbored at least one or a combination of the tested genes and V. choleraestrain BTPR5 isolated from prawn hosted the highest number of virulence associated genes. Among the V. vulnificusstrains, only 3 virulence genes, VPI, toxR and cps, were confirmed out of the 16 tested and only 7 of the isolates had these genes in one or more combinations. Strain BTPS6 from aquafarm and strain BTVE4 from mangrove samples yielded positive amplification for the three genes. The toxRgene from 9 strains of V. choleraeand 3 strains of V. vulnificus were cloned and sequenced for phylogenetic analysis based on nucleotide and the amino acid sequences. Multiple sequence alignment of the nucleotide sequences and amino acid sequences of the environmental strains of V. choleraerevealed that the toxRgene in the environmental strains are 100% homologous to themselves and to the V. choleraetoxR gene sequence available in the Genbank database. The 3 strains of V. vulnificus displayed high nucleotide and amino acid sequence similarity among themselves and to the sequences of V. cholerae and V. harveyi obtained from the GenBank database, but exhibited only 72% homology to the sequences of its close relative V. vulnificus. Structure prediction of the ToxR protein of Vibrio cholerae strain BTMA5 was by PHYRE2 software. The deduced amino acid sequence showed maximum resemblance with the structure of DNA-binding domain of response regulator2 from Escherichia coli k-12 Template based homology modelling in PHYRE2 successfully modelled the predicted protein and its secondary structure based on protein data bank (PDB) template c3zq7A. The pathogenicity studies were performed using the nematode Caenorhabditiselegansas a model system. The assessment of pathogenicity of environmental strain of V. choleraewas conducted with E. coli strain OP50 as the food source in control plates, environmental V. cholerae strain BTOS6, negative for all tested virulence genes, to check for the suitability of Vibrio sp. as a food source for the nematode;V. cholerae Co 366 ElTor, a clinical pathogenic strain and V. cholerae strain BTPR5 from seafood (Prawn) and positive for the tested virulence genes like VPI, hlyA, ompU,rtxA and toxR. It was found that V. cholerae strain BTOS6 could serve as a food source in place of E. coli strain OP50 but behavioral aberrations like sluggish movement and lawn avoidance and morphological abnormalities like pharyngeal and intestinal distensions and bagging were exhibited by the worms fed on V. cholerae Co 366 ElTor strain and environmental BTPR5 indicating their pathogenicity to the nematode. Assessment of pathogenicity of the environmental strains of V. vulnificus was performed with V. vulnificus strain BTPS6 which tested positive for 3 virulence genes, namely, cps, toxRand VPI, and V. vulnificus strain BTMM7 that did not possess any of the tested virulence genes. A reduction was observed in the life span of worms fed on environmental strain of V. vulnificusBTMM7 rather than on the ordinary laboratory food source, E. coli OP50. Behavioral abnormalities like sluggish movement, lawn avoidance and bagging were also observed in the worms fed with strain BTPS6, but the pharynx and the intestine were intact. The presence of multi drug resistant environmental Vibrio strainsthat constitute a major reservoir of diverse virulence genes are to be dealt with caution as they play a decisive role in pathogenicity and horizontal gene transfer in the marine environments.

Protein Folding Simulations: Confinement, External Fields and Sequence Design

The present Thesis looks at the problem of protein folding using Monte Carlo and Langevin simulations, three topics in protein folding have been studied: 1) the effect of confining potential barriers, 2) the effect of a static external field and 3) the design of amino acid sequences which fold in a short time and which have a stable native state (global minimum). Regarding the first topic, we studied the confinement of a small protein of 16 amino acids known as 1NJ0 (PDB code) which has a beta-sheet structure as a native state. The confinement of proteins occurs frequently in the cell environment. Some molecules called Chaperones, present in the cytoplasm, capture the unfolded proteins in their interior and avoid the formation of aggregates and misfolded proteins. This mechanism of confinement mediated by Chaperones is not yet well understood. In the present work we considered two kinds of potential barriers which try to mimic the confinement induced by a Chaperon molecule. The first kind of potential was a purely repulsive barrier whose only effect is to create a cavity where the protein folds up correctly. The second kind of potential was a barrier which includes both attractive and repulsive effects. We performed Wang-Landau simulations to calculate the thermodynamical properties of 1NJ0. From the free energy landscape plot we found that 1NJ0 has two intermediate states in the bulk (without confinement) which are clearly separated from the native and the unfolded states. For the case of the purely repulsive barrier we found that the intermediate states get closer to each other in the free energy landscape plot and eventually they collapse into a single intermediate state. The unfolded state is more compact, compared to that in the bulk, as the size of the barrier decreases. For an attractive barrier modifications of the states (native, unfolded and intermediates) are observed depending on the degree of attraction between the protein and the walls of the barrier. The strength of the attraction is measured by the parameter $\epsilon$. A purely repulsive barrier is obtained for $\epsilon=0$ and a purely attractive barrier for $\epsilon=1$. The states are changed slightly for magnitudes of the attraction up to $\epsilon=0.4$. The disappearance of the intermediate states of 1NJ0 is already observed for $\epsilon =0.6$. A very high attractive barrier ($\epsilon \sim 1.0$) produces a completely denatured state. In the second topic of this Thesis we dealt with the interaction of a protein with an external electric field. We demonstrated by means of computer simulations, specifically by using the Wang-Landau algorithm, that the folded, unfolded, and intermediate states can be modified by means of a field. We have found that an external field can induce several modifications in the thermodynamics of these states: for relatively low magnitudes of the field ($<2.06 \times 10^8$ V/m) no major changes in the states are observed. However, for higher magnitudes than ($6.19 \times 10^8$ V/m) one observes the appearance of a new native state which exhibits a helix-like structure. In contrast, the original native state is a $\beta$-sheet structure. In the new native state all the dipoles in the backbone structure are aligned parallel to the field. The design of amino acid sequences constitutes the third topic of the present work. We have tested the Rate of Convergence criterion proposed by D. Gridnev and M. Garcia ({\it work unpublished}). We applied it to the study of off-lattice models. The Rate of Convergence criterion is used to decide if a certain sequence will fold up correctly within a relatively short time. Before the present work, the common way to decide if a certain sequence was a good/bad folder was by performing the whole dynamics until the sequence got its native state (if it existed), or by studying the curvature of the potential energy surface. There are some difficulties in the last two approaches. In the first approach, performing the complete dynamics for hundreds of sequences is a rather challenging task because of the CPU time needed. In the second approach, calculating the curvature of the potential energy surface is possible only for very smooth surfaces. The Rate of Convergence criterion seems to avoid the previous difficulties. With this criterion one does not need to perform the complete dynamics to find the good and bad sequences. Also, the criterion does not depend on the kind of force field used and therefore it can be used even for very rugged energy surfaces.

Factores perinatales para riesgo de sobrepeso y obesidad en niños de 4 a 5 años en Colombia 2010.

Introducción: La OMS revela que en 2010 alrededor de 43 millones de niños menores de 5 años presentan sobrepeso. En Colombia según la Encuesta Nacional de Situación Nutricional en Colombia en su versión 2005, mostraba una prevalencia general de sobrepeso de 3.1% niños de 0 a 4 años. Es una condición de salud de origen multifactorial en la que interviene factores genéticos, ambientales, maternos y perinatales. Objetivo: Establecer la asociación de riesgo entre el bajo peso al nacer y el desarrollo de sobrepeso y obesidad en niños de 4 a 5 años. Metodología: Se realizó un estudio observacional descriptivo retrospectivo de corte transversal con los datos nutricionales, maternos y perinatales de la Encuesta Nacional de Demografía en Salud del año 2010 en Colombia. Se analizó la asociación entre la variable independiente bajo peso al nacer con el desenlace sobrepeso y obesidad en menores de 4 a 5 años, usando como medida el IMC según la edad. Se realizaron análisis univariados, bivariados y de regresión logística con un modelo de riesgo según las variables que inciden en el desenlace y la variable independiente. Resultados: La muestra obtenida para el estudio fue de 2166 niños de 4 a 5 años de edad quienes cumplían los criterios de inclusión. La prevalencia de sobrepeso u obesidad en la primera infancia fue de 21.8% (472) y el bajo peso al nacer. Los resultados sugieren la asociación de bajo peso y sobrepeso u obesidad es de ORajustado= 0.560 (0.356 – 0.881). Conclusiones: Los resultados sugieren que existe una asociación como factor protector entre el bajo peso y el sobrepeso u obesidad en la primera infancia. Sin embargo, debido al comportamiento de las variables consideradas en la muestra no hay suficiente información para rechazar completamente la hipótesis nula.

The structure of echovirus type 12 bound to a two-domain fragment of its cellular attachment protein decay-accelerating factor (CD 55)

Echovirus type 12 (EV12), an enterovirus of the Picornaviridae family, uses the complement regulator, decay-accelerating factor (DAF, CD55) as a cellular receptor. We have calculated a three-dimensional reconstruction of EV12 bound to a fragment of DAF, consisting of short consensus repeat domains 3 and 4, from cryo-negative stain electron microscopy data (EMD #1057). This shows that, as for an earlier reconstruction of the related echovirus type 7 bound to DAF, attachment is not within the viral canyon but occurs close to the two-fold symmetry axes. Despite this general similarity, our reconstruction reveals a receptor interaction that is quite different from that observed for EV7. Fitting of the crystallographic co-ordinates for DAF34 and EV11 into the reconstruction shows a close agreement between the crystal structure of the receptor fragment and the density for the virus-bound receptor, allowing unambiguous positioning of the receptor with respect to the virion (PDB #1UPN). Our finding that the mode of virus-receptor interaction in EV12 is distinct from that seen for EV7 raises interesting questions regarding the evolution and biological significance of the DAF-binding phenotype in these viruses.

Both tumour-promoting and non-promoting phorbol esters inhibit 125I-EGF binding and stimulate the phosphorylation of a 80kd protein kinase C substrate protein in intact Swiss 3T3 cells

Sapintoxin A (SAP A) and 12-deoxyphorbol 13-phenylacetate (DOPP), are two biologically active but non-turnour-promoting phorbol esters that potently bind to and activate the phorbol ester receptor, protein kinase C (PKC). SAP A and DOPP cause a dose-dependent increase in the phosphorylation of an 80 kd (80K) substrate protein for PKC in Swiss 3T3 cells. A similar dose—response effect was seen with sapintoxin D (SAP D), the stage 2 promoting analogue of 12-O-tetradecanoylphorbol-13-acetate and the complete promoter phorbol 12,13-dibutyrate (PDB). The doses resulting in a half maximal phosphorylation of this protein (Ka were 20 nM (SAP A), 45 nM (DOPP), 23 nM (SAP D) and 37 nM (PDB). Both non-promoting and phorbol esters induced a dose-dependent inhibition of [125I]epidermal growth factor (EGF) binding to its receptor in Swiss 3T3 cells. The doses required for 50% inhibition of binding (Ki) were: 8 nM (SAP A), 16 nM (DOPP), 14 nM (SAP D) and 17 nM (PDB). The results clearly demonstrate that induction of phosphorylation of the Pu 80K phosphoprotein and inhibition of [125I]EGF binding in Swiss 3T3 cells following exposure to phorbol esters is independent of the tumour-promoting activity of these compounds. The fact that SAP A, DOPP, SAP D and PDB are mitogenic for a variety of cell types and that exposure to these compounds leads to 80K phosphorylation and inhibition of [125I]EGF binding, suggests that these early biological events may play a role in the mitogenic response induced by these compounds.

Molecular modelling studies of binding of DACD derivatives into G-Quadruplex DNA: comparison of force field and quantum polarized ligand docking methods

The DNA G-qadruplexes are one of the targets being actively explored for anti-cancer therapy by inhibiting them through small molecules. This computational study was conducted to predict the binding strengths and orientations of a set of novel dimethyl-amino-ethyl-acridine (DACA) analogues that are designed and synthesized in our laboratory, but did not diffract in Synchrotron light.Thecrystal structure of DNA G-Quadruplex(TGGGGT)4(PDB: 1O0K) was used as target for their binding properties in our studies.We used both the force field (FF) and QM/MM derived atomic charge schemes simultaneously for comparing the predictions of drug binding modes and their energetics. This study evaluates the comparative performance of fixed point charge based Glide XP docking and the quantum polarized ligand docking schemes. These results will provide insights on the effects of including or ignoring the drug-receptor interfacial polarization events in molecular docking simulations, which in turn, will aid the rational selection of computational methods at different levels of theory in future drug design programs. Plenty of molecular modelling tools and methods currently exist for modelling drug-receptor or protein-protein, or DNA-protein interactionssat different levels of complexities.Yet, the capasity of such tools to describevarious physico-chemical propertiesmore accuratelyis the next step ahead in currentresearch.Especially, the usage of most accurate methods in quantum mechanics(QM) is severely restricted by theirtedious nature. Though the usage of massively parallel super computing environments resulted in a tremendous improvement in molecular mechanics (MM) calculations like molecular dynamics,they are still capable of dealing with only a couple of tens to hundreds of atoms for QM methods. One such efficient strategy that utilizes thepowers of both MM and QM are the QM/MM hybrid methods. Lately, attempts have been directed towards the goal of deploying several different QM methods for betterment of force field based simulations, but with practical restrictions in place. One of such methods utilizes the inclusion of charge polarization events at the drug-receptor interface, that is not explicitly present in the MM FF.

Insights into the Specificity of Thioredoxin Reductase-Thioredoxin Interactions. A Structural and Functional Investigation of the Yeast Thioredoxin System

The enzymatic activity of thioredoxin reductase enzymes is endowed by at least two redox centers: a flavin and a dithiol/disulfide CXXC motif. The interaction between thioredoxin reductase and thioredoxin is generally species-specific, but the molecular aspects related to this phenomenon remain elusive. Here, we investigated the yeast cytosolic thioredoxin system, which is composed of NADPH, thioredoxin reductase (ScTrxR1), and thioredoxin 1 (ScTrx1) or thioredoxin 2 (ScTrx2). We showed that ScTrxR1 was able to efficiently reduce yeast thioredoxins (mitochondrial and cytosolic) but failed to reduce the human and Escherichia coli thioredoxin counterparts. To gain insights into this specificity, the crystallographic structure of oxidized ScTrxR1 was solved at 2.4 angstrom resolution. The protein topology of the redox centers indicated the necessity of a large structural rearrangement for FAD and thioredoxin reduction using NADPH. Therefore, we modeled a large structural rotation between the two ScTrxR1 domains (based on the previously described crystal structure, PDB code 1F6M). Employing diverse approaches including enzymatic assays, site-directed mutagenesis, amino acid sequence alignment, and structure comparisons, insights were obtained about the features involved in the species-specificity phenomenon, such as complementary electronic parameters between the surfaces of ScTrxR1 and yeast thioredoxin enzymes and loops and residues (such as Ser(72) in ScTrx2). Finally, structural comparisons and amino acid alignments led us to propose a new classification that includes a larger number of enzymes with thioredoxin reductase activity, neglected in the low/high molecular weight classification.

The CATH classification revisited-architectures reviewed and new ways to characterize structural divergence in superfamilies

The latest version of CATH (class, architecture, topology, homology) (version 3.2), released in July 2008 (http://www.cathdb.info), contains 1 14215 domains, 2178 Homologous superfamilies and 1110 fold groups. We have assigned 20 330 new domains, 87 new homologous superfamilies and 26 new folds since CATH release version 3.1. A total of 28 064 new domains have been assigned since our NAR 2007 database publication (CATH version 3.0). The CATH website has been completely redesigned and includes more comprehensive documentation. We have revisited the CATH architecture level as part of the development of a `Protein Chart` and present information on the population of each architecture. The CATHEDRAL structure comparison algorithm has been improved and used to characterize structural diversity in CATH superfamilies and structural overlaps between superfamilies. Although the majority of superfamilies in CATH are not structurally diverse and do not overlap significantly with other superfamilies, similar to 4% of superfamilies are very diverse and these are the superfamilies that are most highly populated in both the PDB and in the genomes. Information on the degree of structural diversity in each superfamily and structural overlaps between superfamilies can now be downloaded from the CATH website.

Fragment-guided approach to incorporating structural information into a CoMFA study: BACE-1 as an example

Alzheimer`s disease is an ultimately fatal neurodegenerative disease, and BACE-1 has become an attractive validated target for its therapy, with more than a hundred crystal structures deposited in the PDB. In the present study, we present a new methodology that integrates ligand-based methods with structural information derived from the receptor. 128 BACE-1 inhibitors recently disclosed by GlaxoSmithKline R&D were selected specifically because the crystal structures of 9 of these compounds complexed to BACE-1, as well as five closely related analogs, have been made available. A new fragment-guided approach was designed to incorporate this wealth of structural information into a CoMFA study, and the methodology was systematically compared to other popular approaches, such as docking, for generating a molecular alignment. The influence of the partial charges calculation method was also analyzed. Several consistent and predictive models are reported, including one with r (2) = 0.88, q (2) = 0.69 and r (pred) (2) = 0.72. The models obtained with the new methodology performed consistently better than those obtained by other methodologies, particularly in terms of external predictive power. The visual analyses of the contour maps in the context of the enzyme drew attention to a number of possible opportunities for the development of analogs with improved potency. These results suggest that 3D-QSAR studies may benefit from the additional structural information added by the presented methodology.

One, two or no record of late neoproterozoic glaciation in South-East Cameroon?

Severe climate changes culminating in at least three major glacial events have been recognized in the Neoproterozoic sedimentary record from many parts of the world Supportive to the global nature of these climatic shifts a considerable amount of data have been acquired from deposits exposed in Pan-African orogenic belts in southwestern and western Africa By comparison published data from the Pan-African belts in Central Africa are scarce We report here evidence of possibly two glacial events recorded in the Mintom Formation that is located on the margin of the Pan-African orogenic Yaounde belt in South-East Cameroon In the absence of reliable radiometric data only maximum and minimum age limits of 640 and 580 Ma respectively can at present be applied to the Mintom Formation The formation consists of two lithostratigraphic ensembles each subdivided in two members (i e in ascending stratigraphic order the Kol Metou Momibole and Atog Adjap Members) The basal ensemble exhibits a typical glacial to post-glacial succession It includes diamictites comprising cobbles and boulders in a massive argillaceous siltstone matrix and laminated siltstones followed by in sharp contact a 2 m-thick massive dolostone that yielded negative delta(13)C values (<-3 parts per thousand. V-PDB) similar to those reported for Marinoan cap carbonates elsewhere However uncertainty remains regarding the glacial influence on the siliciclastic facies because the diamictite is better explained as a mass-flow deposit and diagnostic features such as dropstones have not been seen in the overlying siltstones The Mintom Formation may thus provide an example of an unusual succession of non-glacial diamictite overlain by a truly glacial melt-related cap-carbonate We also report the recent discovery of ice-striated pavements on the structural surface cut in the Mintom Formation suggesting that glaciers developed after the latter had been deposited and deformed during the Pan-African orogeny Striations which consistently exhibit two principal orientations (N60 and N110) were identified in two different localities in the west of the study area on siltstones of the Kol Member and in the east on limestones of the Atog Adjap Member respectively N60-oriented striae indicate ice flow towards the WSW Assigning an age to these features remains problematical because they were not found associated with glaciogenic deposits Two hypotheses can equally be envisaged e either the striated surfaces are correlated (1) to the Gaskiers (or Neoproterozoic post-Gaskiers) glaciation and represent the youngest Ediacaran glacial event documented in the southern Yaounde belt or (2) to the Late Ordovician Hirnantian (Saharan) glaciation thereby providing new data about Hirnantian ice flows in Central Africa (C) 2010 Elsevier Ltd All rights reserved