Database resources of the National Center for Biotechnology Information

In addition to maintaining the GenBank® nucleic acid sequence database, the National Center for Biotechnology Information (NCBI) provides data analysis and retrieval resources that operate on the data in GenBank and a variety of other biological data made available through NCBI’s Web site. NCBI data retrieval resources include Entrez, PubMed, LocusLink and the Taxonomy Browser. Data analysis resources include BLAST, Electronic PCR, OrfFinder, RefSeq, UniGene, HomoloGene, Database of Single Nucleotide Polymorphisms (dbSNP), Human Genome Sequencing, Human MapViewer, GeneMap’99, Human–Mouse Homology Map, Cancer Chromosome Aberration Project (CCAP), Entrez Genomes, Clusters of Orthologous Groups (COGs) database, Retroviral Genotyping Tools, Cancer Genome Anatomy Project (CGAP), SAGEmap, Gene Expression Omnibus (GEO), Online Mendelian Inheri­tance in Man (OMIM), the Molecular Modeling Database (MMDB) and the Conserved Domain Database (CDD). Augmenting many of the Web applications are custom implementations of the BLAST program optimized to search specialized data sets. All of the resources can be accessed through the NCBI home page at: http://www.ncbi.nlm.nih.gov.

Selenoprotein oxidoreductase with specificity for thioredoxin and glutathione systems

Thioredoxin (Trx) and glutathione (GSH) systems are considered to be two major redox systems in animal cells. They are reduced by NADPH via Trx reductase (TR) or oxidized GSH (GSSG) reductase and further supply electrons for deoxyribonucleotide synthesis, antioxidant defense, and redox regulation of signal transduction, transcription, cell growth, and apoptosis. We cloned and characterized a pyridine nucleotide disulfide oxidoreductase, Trx and GSSG reductase (TGR), that exhibits specificity for both redox systems. This enzyme contains a selenocysteine residue encoded by the TGA codon. TGR can reduce Trx, GSSG, and a GSH-linked disulfide in in vitro assays. This unusual substrate specificity is achieved by an evolutionary conserved fusion of the TR and glutaredoxin domains. These observations, together with the biochemical probing and molecular modeling of the TGR structure, suggest a mechanism whereby the C-terminal selenotetrapeptide serves a role of a protein-linked GSSG and shuttles electrons from the disulfide center within the TR domain to either the glutaredoxin domain or Trx.

Hybrid characteristics of oligonucleotides consisting of isonucleoside 2′,5′-anhydro-3′-deoxy-3′-(thymin-1-yl)-d-mannitol with different linkage modes

Oligonucleotides consisting of the isonucleoside repeating unit 2′,5′-anhydro-3′-deoxy-3′-(thymin-1-yl)-d-mannitol (4) were synthesized with the monomeric unit 4 incorporated into oligonucleotides as 1′→4′ linkage 4a (oligomer I) or 6′→4′ linkage 4b (oligomer II). The hybrid properties of the two oligonucleotides I and II with their complementary strands were investigated by thermal denaturation and CD spectra. Oligonucleotide I (4a) formed a stable duplex with d(A)14 with a slightly reduced Tm value of 36.6°C, relative to 38.2°C for the control duplex d(T)14/d(A)14, but oligomer II (4b) failed to hybridize with a DNA complementary single strand. The spectrum of the duplex oligomer I/d(A)14 showed a positive CD band at 217 nm and a negative CD band at 248 nm attributable to a B-like conformation. Molecular modeling showed that in the case of oligomer I the C6′ hydroxy group of each unit could be located in the groove area when hybridized to the DNA single strand, which might contribute additional hydrogen bonding to the stability of duplex formation.

Redesigning secondary structure to invert coenzyme specificity in isopropylmalate dehydrogenase.

Rational engineering of enzymes involves introducing key amino acids guided by a knowledge of protein structure to effect a desirable change in function. To date, all successful attempts to change specificity have been limited to substituting individual amino acids within a protein fold. However, the infant field of protein engineering will only reach maturity when changes in function can be generated by rationally engineering secondary structures. Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to systematically invert the coenzyme specificity of Thermus thermophilus isopropylmalate dehydrogenase from a 100-fold preference for NAD to a 1000-fold preference for NADP. The engineered mutant, which is twice as active as wild type, contains four amino acid substitutions and an alpha-helix and loop that replaces the original beta-turn. These results demonstrate that rational engineering of secondary structures to produce enzymes with novel properties is feasible.

Structural analysis of p185c-neu and epidermal growth factor receptor tyrosine kinases: oligomerization of kinase domains.

The epidermal growth factor receptor (EGFR) and p185c-neu proteins associate as dimers to create an efficient signaling assembly. Overexpression of these receptors together enhances their intrinsic kinase activity and concomitantly results in oncogenic cellular transformation. The ectodomain is able to stabilize the dimer, whereas the kinase domain mediates biological activity. Here we analyze potential interactions of the cytoplasmic kinase domains of the EGFR and p185c-neu tyrosine kinases by homology molecular modeling. This analysis indicates that kinase domains can associate as dimers and, based on intermolecular interaction calculations, that heterodimer formation is favored over homodimers. The study also predicts that the self-autophosphorylation sites located within the kinase domains are not likely to interfere with tyrosine kinase activity, but may regulate the selection of substrates, thereby modulating signal transduction. In addition, the models suggest that the kinase domains of EGFR and p185c-neu can undergo higher order aggregation such as the formation of tetramers. Formation of tetrameric complexes may explain some of the experimentally observed features of their ligand affinity and hetero-receptor internalization.

Zn2+ interaction with Alzheimer amyloid beta protein calcium channels.

The Alzheimer disease 40-residue amyloid beta protein (AbetaP[1-40]) forms cation-selective channels across acidic phospholipid bilayer membranes with spontaneous transitions over a wide range of conductances ranging from 40 to 4000 pS. Zn2+ has been reported to bind to AbetaP[1-40] with high affinity, and it has been implicated in the formation of amyloid plaques. We now report the functional consequences of such Zn2+ binding for the AbetaP[1-40] channel. Provided the AbetaP[1-40] channel is expressed in the low conductance (<400 pS) mode, Zn2+ blocks the open channel in a dose- dependent manner. For AbetaP[1-40] channels in the giant conductance mode (>400 pS), Zn2+ doses in the millimolar range were required to exert substantial blockade. The Zn2+ chelator o-phenanthroline reverses the blockade. We also found that Zn2+ modulates AbetaP[1-40] channel gating and conductance only from one side of the channel. These data are consistent with predictions of our recent molecular modeling studies on AbetaP[1-40] channels indicating asymmetric Zn(2+)-AbetaP[1-40] interactions at the entrance to the pore.

A highly active decarboxylating dehydrogenase with rationally inverted coenzyme specificity.

The isocitrate dehydrogenase of Escherichia coli, which lacks the Rossmann fold common to other dehydrogenases, displays a 7000-fold preference for NADP over NAD (calculated as the ratio of kcat/Km). Guided by x-ray crystal structures and molecular modeling, site-directed mutagenesis has been used to introduce six substitutions in the adenosine binding pocket that systematically shift coenzyme preference toward NAD. The engineered enzyme displays an 850-fold preference for NAD over NADP, which exceeds the 140-fold preference displayed by a homologous NAD-dependent enzyme. Of the six mutations introduced, only one is identical in all related NAD-dependent enzyme sequences--strict adherence to homology as a criterion for replacing these amino acids impairs function. Two additional mutations at remote sites improve performance further, resulting in a final mutant enzyme with kinetic characteristics and coenzyme preference comparable to naturally occurring homologous NAD-dependent enzymes.

A general model of invariant chain association with class II major histocompatibility complex proteins.

The binding of invariant chain to major histocompatibility complex (MHC) proteins is an important step in processing of MHC class II proteins and in antigen presentation. The question of how invariant chain can bind to all MHC class II proteins is central to understanding these processes. We have employed molecular modeling to predict the structure of class II-associated invariant chain peptide (CLIP)-MHC protein complexes and to ask whether the predicted mode of association could be general across all MHC class II proteins. CLIP fits identically into the MHC class II alleles HLA-DR3, I-Ak, I-Au, and I-Ad, with a consistent pattern of hydrogen bonds, contacts, and hydrophobic burial and without bad contacts. Our model predicts the burial of CLIP residues Met-91 and Met-99 in the deep P1 and P9 anchor pockets and other detailed interactions, which we have compared with available data. The predicted pattern of I-A allele-specific effects on CLIP binding is very similar to that observed experimentally by alanine-scanning mutations of CLIP. Together, these results indicate that CLIP may bind in a single, general way across products of MHC class II alleles.

Síntese de potencial inibidor de acetilcolinesterase para tratamento da Doença de Alzheimer

A doença de Alzheimer (DA) é a forma mais comum de demência, representando cerca de 80% dos casos. A DA é caracterizada por um processo de declínio progressivo e irreversível das funções cognitivas e da memória, que se estende para a desorganização do comportamento. Atualmente, 46,8 milhões de pessoas em todo o mundo foram diagnosticadas com demência. Embora vários fatores tenham sido implicados na DA, sua etiologia ainda não é completamente conhecida. Do ponto de vista neuropatológico, é observado no cérebro de indivíduos com DA atrofia cortical difusa, presença de grande número de placas senis, emaranhados neurofibrilares, processo inflamatório e perda neuronal. A progressão dos sintomas está associada a mudanças estruturais nas sinapses colinérgicas em certas regiões do cérebro, que consequentemente, apresentam neurotransmissão colinérgica reduzida. Os vários eventos patológicos interligados contribuem para o avanço da doença e direcionam diversas pesquisas na busca por tratamentos multialvos com base no processo multifatorial de DA. Assim o presente trabalho descreve a síntese de derivados híbridos dual binding site de donepezila-tacrina (fármacos inibidores de acetilcolinesterase), com potencial para agir em dois alvos terapêuticos pela (i) inibição da acetilcolinesterase em ambos os sítios ativo e periférico, como demonstrado pelos estudos de modelagem molecular, e (ii) na agregação do peptídeo A? neurotóxico induzido pela acetilcolinesterase, na tentativa de interromper a progressão da doença. A estratégia sintética envolveu a condensação da 5,6-dimetóxiindanona com a unidade 4-piperidinil carbaldeído, a qual forneceu o intermediário 5,6- dimetóxindan-1-ona-4-piperidinil-metileno-1-[(4-cloroquinolin-2-il)metil], seguido de redução da dupla ligação, gerada na reação de condensação anterior, e substituição do átomo de cloro-quinolina por amino para obtenção do produto final, ou manutenção da função olefina, seguido de substituição do átomo de cloro-quinolina por azido ou amino, gerando cinco híbridos estruturalmente correlacionados. Os híbridos foram testados em ensaio de inibição de acetilcolinesterase e butirilcolinesterase pelo método de Ellman, e o híbrido insaturado, contendo a função amino-quinolina foi o mais ativo da série com IC50 na faixa de nanomolar (0,014 ?M). Futuramente, os intermediários da reação e produto final serão submetidos ao ensaio de inibição da agregação do peptídeo A? neurotóxico pelo método da tioflavina T. Neste trabalho, também são descritos os testes de predição in vitro para permeação pela barreira hematoencefálica, bem como sua absorção intestinal, pelo método PAMPA.

Structural basis and prediction of substrate specificity in protein serine/threonine kinases

The large number of protein kinases makes it impractical to determine their specificities and substrates experimentally. Using the available crystal structures, molecular modeling, and sequence analyses of kinases and substrates, we developed a set of rules governing the binding of a heptapeptide substrate motif (surrounding the phosphorylation site) to the kinase and implemented these rules in a web-interfaced program for automated prediction of optimal substrate peptides, taking only the amino acid sequence of a protein kinase as input. We show the utility of the method by analyzing yeast cell cycle control and DNA damage checkpoint pathways. Our method is the only available predictive method generally applicable for identifying possible substrate proteins for protein serine/threonine kinases and helps in silico construction of signaling pathways. The accuracy of prediction is comparable to the accuracy of data from systematic large-scale experimental approaches.

Active site mutations and substrate inhibition in human sulfotransferase 1A1 and 1A3

Human SULT1A1 is primarily responsible for sulfonation of xenobiotics, including the activation of promutagens, and it has been implicated in several forms of cancer. Human SULT1A3 has been shown to be the major sulfotransferase that sulfonates dopamine. These two enzymes shares 93% amino acid sequence identity and have distinct but overlapping substrate preferences. The resolution of the crystal structures of these two enzymes has enabled us to elucidate the mechanisms controlling their substrate preferences and inhibition. The presence of two p-nitrophenol (pNP) molecules in the crystal structure of SULT1A1 was postulated to explain cooperativity at low and inhibition at high substrate concentrations, respectively. In SULT1A1, substrate inhibition occurs with pNP as the substrate but not with dopamine. For SULT1A3, substrate inhibition is found for dopamine but not with pNP. We investigated how substrate inhibition occurs in these two enzymes using molecular modeling, site-directed mutagenesis, and kinetic analysis. The results show that residue Phe-247 of SULT1A1, which interacts with both p-nitrophenol molecules in the active site, is important for substrate inhibition. Mutation of phenylalanine to leucine at this position in SULT1A1 results in substrate inhibition by dopamine. We also propose, based on modeling and kinetic studies, that substrate inhibition by dopamine in SULT1A3 is caused by binding of two dopamine molecules in the active site. © 2004 by The American Society for Biochemistry and Molecular Biology, Inc.

Human sulfotransferases and their role in chemical metabolism

Sulfonation is an important reaction in the metabolism of numerous xenobiotics, drugs, and endogenous compounds. A supergene family of enzymes called sulfotransferases (SULTs) catalyze this reaction. In most cases, the addition of a sulfonate moiety to a compound increases its water solubility and decreases its biological activity. However, many of these enzymes are also capable of bioactivating procarcinogens to reactive electrophiles. In humans three SULT families, SULT1, SULT2, and SULT4, have been identified that contain at least thirteen distinct members. SULTs have a wide tissue distribution and act as a major detoxification enzyme system in adult and the developing human fetus. Nine crystal structures of human cytosolic SULTs have now been determined, and together with site-directed mutagenesis experiments and molecular modeling, we are now beginning to understand the factors that govern distinct but overlapping substrate specificities. These studies have also provided insight into the enzyme kinetics and inhibition characteristics of these enzymes. The regulation of human SULTs remains as one of the least explored areas of research in the field, though there have been some recent advances on the molecular transcription mechanism controlling the individual SULT promoters. Interindividual variation in sulfonation capacity may be important in determining an individual's response to xenobiotics, and recent studies have begun to suggest roles for SULT polymorphism in disease susceptibility. This review aims to provide a summary of our present understanding of the function of human cytosolic sulfotransferases.

Comparison of the binding of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines with their isosteric sulfonamides to the active site of phenylethanolamine N-methyltransferase

3-Fluoromethyl-7-(N-substituted aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines (14, 16, and 18-22) are highly potent and selective inhibitors of phenylethanolamine N-methyltransferase (PNMT). Molecular modeling studies with 3-fluoromethyl-7-(N-alkyl aminosulfonyl)-1,2,3,4-tetrahydroisoquinolines, such as 16, suggested that the sulfonamide -NH-could form a hydrogen bond with the side chain of Lys57. However, SAR studies and analysis of the crystal structure of human PNMT (hPNMT) in complex with 7 indicated that the sulfonamide oxygens, and not the sulfonamide -NH-, formed favorable interactions with the enzyme. Thus, we hypothesized that replacement of the sulfonamide -NH-with a methylene group could result in compounds that would retain potency at PNMT and that would have increased lipophilicity, thus increasing the likelihood they will cross the blood brain barrier. A series of 3-fluoromethyl-7-sulfonyl-1,2,3,4-tetrahydroisoquinolines (23-30) were synthesized and evaluated for their PNMT inhibitory potency and affinity for the R2-adrenoceptor. A comparison of these compounds with their isosteric sulfonamides (14, 16, and 18-22) showed that the sulfones were more lipophilic but less potent than their corresponding sulfonamides. Sulfone 24 (hPNMT K-i = 1.3 mu M) is the most potent compound in this series and is quite selective for PNMT versus the R2-adrenoceptor, but 24 is less potent than the corresponding sulfonamide, 16 (hPNMT K-i = 0.13 mu M). We also report the crystal structure of hPNMT in complex with sulfonamide 15, from which a potential hydrogen bond acceptor within the hPNMT active site has been identified, the main chain carbonyl oxygen of Asn39. The interaction of this residue with the sulfonamide -NH-is likely responsible for much of the enhanced inhibitory potency of the sulfonamides versus the sulfones.

Electronic energy-transfer rate constants for geometrical isomers of a bichromophoric macrocyclic complex

The rate of electronic energy transfer (EET) between a naphthalene donor and an anthracene acceptor in [ZnL3]-(ClO4)(2) and [ZnL4](ClO4)(2) was determined by time-resolved fluorescence measurements, where L 3 and L 4 are the geometrical isomers of 6-[(anthracen-9-ylmethyl)amino]-trans-6,13-dimethyl-1,4,8,11-tetraazacyclotetradecane-13-amine (L-2), substituted with either a naphthalen-1-ylmethyl or naphthalen-2-ylmethyl donor, respectively. The energy transfer rate constant, k(EET), was determined to be (0.92 +/- 0.02) x 10(9) s(-1) for the naphthalen-1-ylmethyl-substituted isomer, while that for the naphthalen-2-ylmethyl-substituted isomer is somewhat faster, with k(EET) = (1.31 +/- 0.01) x 10(9) s(-1). The solid-state structure of [(ZnLCl)-Cl-3]ClO4 has been determined, and using molecular modeling calculations, the likely distributions of solution conformations in CH3CN have been evaluated for both complexes. The calculated conformational distributions in the common trans-III N-based isomeric form gave Forster EET rate constants that account for the differences observed and are in excellent agreement with the experimental values. It is shown that the full range of conformers must be considered to accurately reproduce the observed EET kinetics.

Multi-scale approaches to Food Sciences: computational QM and MM explorations at the microscopic level

L'obiettivo principale della politica di sicurezza alimentare è quello di garantire la salute dei consumatori attraverso regole e protocolli di sicurezza specifici. Al fine di rispondere ai requisiti di sicurezza alimentare e standardizzazione della qualità, nel 2002 il Parlamento Europeo e il Consiglio dell'UE (Regolamento (CE) 178/2002 (CE, 2002)), hanno cercato di uniformare concetti, principi e procedure in modo da fornire una base comune in materia di disciplina degli alimenti e mangimi provenienti da Stati membri a livello comunitario. La formalizzazione di regole e protocolli di standardizzazione dovrebbe però passare attraverso una più dettagliata e accurata comprensione ed armonizzazione delle proprietà globali (macroscopiche), pseudo-locali (mesoscopiche), ed eventualmente, locali (microscopiche) dei prodotti alimentari. L'obiettivo principale di questa tesi di dottorato è di illustrare come le tecniche computazionali possano rappresentare un valido supporto per l'analisi e ciò tramite (i) l’applicazione di protocolli e (ii) miglioramento delle tecniche ampiamente applicate. Una dimostrazione diretta delle potenzialità già offerte dagli approcci computazionali viene offerta nel primo lavoro in cui un virtual screening basato su docking è stato applicato al fine di valutare la preliminare xeno-androgenicità di alcuni contaminanti alimentari. Il secondo e terzo lavoro riguardano lo sviluppo e la convalida di nuovi descrittori chimico-fisici in un contesto 3D-QSAR. Denominata HyPhar (Hydrophobic Pharmacophore), la nuova metodologia così messa a punto è stata usata per esplorare il tema della selettività tra bersagli molecolari strutturalmente correlati e ha così dimostrato di possedere i necessari requisiti di applicabilità e adattabilità in un contesto alimentare. Nel complesso, i risultati ci permettono di essere fiduciosi nel potenziale impatto che le tecniche in silico potranno avere nella identificazione e chiarificazione di eventi molecolari implicati negli aspetti tossicologici e nutrizionali degli alimenti.