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.

Diverse functional motifs within the three intracellular loops of the CGRP1 receptor

The CGRP1 receptor exists as a heterodimeric complex between a single-pass transmembrane accessory protein (RAMP1) and a family B G-protein-coupled receptor (GPCR) called the calcitonin receptor-like receptor (CLR). This study investigated the structural motifs found in the intracellular loops (ICLs) of this receptor. Molecular modeling was used to predict active and inactive conformations of each ICL. Conserved residues were altered to alanine by site-directed mutagenesis. cAMP accumulation, cell-surface expression, agonist affinity, and CGRP-stimulated receptor internalization were characterized. Within ICL1, L147 and particularly R151 were important for coupling to Gs. R151 may interact directly with the G-protein, accessing it following conformational changes involving ICL2 and ICL3. At the proximal end of ICL3, I290 and L294, probably lying on the same face of an α helix, formed a G-protein coupling motif. The largest effects on coupling were observed with I290A; additionally, it reduced CGRP affinity and impaired internalization. 1290 may interact with TM6 to stabilize the conformation of ICL3, but it could also interact directly with Gs. R314, at the distal end of ICL3, impaired G-protein coupling and to a lesser extent reduced CGRP affinity; it may stabilize the TM6-ICL3 junction by interacting with the polar headgroups of membrane phospholipids. Y215 and L214 in ICL2 are required for cell-surface expression; they form a microdomain with H216 which has the same function. This study reveals similarities between the activation of CLR and other GPCRs in the role of TM6 and ICL3 but shows that other conserved motifs differ in their function. © 2006 American Chemical Society.

The design of novel antifolates against pneumocystis carinii

The pneumonia caused by Pneumocystis carinii is ultimately responsible for the death of many acquired immunodeficiency syndrome (AIDS) patients. Large doses of trimethoprim and pyrimethamine in combination with a sulphonamide and/or pentamidine suppress the infection but produce serious side-effects and seldom prevent recurrence after treatment withdrawal. However, the partial success of the aforementioned antifolates, and also trimetrexate used alone, does suggest dihydrofolate reductase (DHFR) as a target for the development of antipneumocystis agents. From the DHFR inhibitory activities of 3'-substituted pyrimethamine analogues it was suggested that the 3'-(3'',3''-dimethyltriazen-1''-yl) substituent may be responsible for the greater activity for the P.carinii over the mammalian enzyme. Crystallographic and molecular modeling studies revealed considerable geometrical and electronic differences between the triazene and the chemically related formamidine functions that may account for the differences in DHFR inhibitory profiles. Structural and electronic parameters calculated for a series of 3'-(3'',3''-disubstitutedtriazen-1''-yl) pyrimethamine analogues did not correlate with the DHFR inhibitory activities. However, the in vitro screening against P.carinii DHFR revealed that the 3''-hydroxyethyl-3''-benzyl analogue was the most active and selective. Models of the active sites of human and P.carinii DHFRs were constructed using DHFR sequence and structural homology data which had identified key residues involved in substrate and cofactor binding. Low energy conformations of the 3'',3''-dimethyl and 3''-hydroxyethyl-3''-benzyle analogues, determined from nuclear magnetic resonance studies and theoretical calculations, were docked by superimposing the diaminopyrimidine fragment onto a previously docked pyrimethamine analogue. Enzyme kinetic data supported the 3''-hydroxyethyl-3''-benzyl moiety being located in the NADPH binding groove. The 3''-benzyl substituent was able to locate to within 3 AA of a valine residue in the active site of P.carinii DHFR thereby producing a hydrophobic contact. The equivalent residue in human DHFR is threonine, more hydrophilic and less likely to be involved in such a contact. This difference may account for the greater inhibitory activity this analogue has for P.carinii DHFR and provide a basis for future drug design. From an in vivo model of PCP in immunosuppressed rats it was established that the 3"-hydroxyethyl-3"-benzyl analogue was able to reduce the.P.carinii burden more effectively with increasing doses, without causmg any visible signs of toxicity. However, equivalent doses were not as effective as pentamidine, a current treatment of choice for Pneumocystis carinii pneumonia.

Static energy analysis of MHC class I and class II peptide-binding affinity

Antigenic peptide is presented to a T-cell receptor (TCR) through the formation of a stable complex with a major histocompatibility complex (MHC) molecule. Various predictive algorithms have been developed to estimate a peptide's capacity to form a stable complex with a given MHC class II allele, a technique integral to the strategy of vaccine design. These have previously incorporated such computational techniques as quantitative matrices and neural networks. A novel predictive technique is described, which uses molecular modeling of predetermined crystal structures to estimate the stability of an MHC class II-peptide complex. The structures are remodeled, energy minimized, and annealed before the energetic interaction is calculated.

Toward the prediction of class I and II mouse major histocompatibility complex-peptide-binding affinity:in silico bioinformatic step-by-step guide using quantitative structure-activity relationships

Quantitative structure-activity relationship (QSAR) analysis is a cornerstone of modern informatics. Predictive computational models of peptide-major histocompatibility complex (MHC)-binding affinity based on QSAR technology have now become important components of modern computational immunovaccinology. Historically, such approaches have been built around semiqualitative, classification methods, but these are now giving way to quantitative regression methods. We review three methods--a 2D-QSAR additive-partial least squares (PLS) and a 3D-QSAR comparative molecular similarity index analysis (CoMSIA) method--which can identify the sequence dependence of peptide-binding specificity for various class I MHC alleles from the reported binding affinities (IC50) of peptide sets. The third method is an iterative self-consistent (ISC) PLS-based additive method, which is a recently developed extension to the additive method for the affinity prediction of class II peptides. The QSAR methods presented here have established themselves as immunoinformatic techniques complementary to existing methodology, useful in the quantitative prediction of binding affinity: current methods for the in silico identification of T-cell epitopes (which form the basis of many vaccines, diagnostics, and reagents) rely on the accurate computational prediction of peptide-MHC affinity. We have reviewed various human and mouse class I and class II allele models. Studied alleles comprise HLA-A*0101, HLA-A*0201, HLA-A*0202, HLA-A*0203, HLA-A*0206, HLA-A*0301, HLA-A*1101, HLA-A*3101, HLA-A*6801, HLA-A*6802, HLA-B*3501, H2-K(k), H2-K(b), H2-D(b) HLA-DRB1*0101, HLA-DRB1*0401, HLA-DRB1*0701, I-A(b), I-A(d), I-A(k), I-A(S), I-E(d), and I-E(k). In this chapter we show a step-by-step guide into predicting the reliability and the resulting models to represent an advance on existing methods. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made are freely available online at the URL http://www.jenner.ac.uk/MHCPred.

Toward prediction of class II mouse major histocompatibility complex peptide binding affinity:in silico bioinformatic evaluation using partial least squares, a robust multivariate statistical technique

The accurate identification of T-cell epitopes remains a principal goal of bioinformatics within immunology. As the immunogenicity of peptide epitopes is dependent on their binding to major histocompatibility complex (MHC) molecules, the prediction of binding affinity is a prerequisite to the reliable prediction of epitopes. The iterative self-consistent (ISC) partial-least-squares (PLS)-based additive method is a recently developed bioinformatic approach for predicting class II peptide−MHC binding affinity. The ISC−PLS method overcomes many of the conceptual difficulties inherent in the prediction of class II peptide−MHC affinity, such as the binding of a mixed population of peptide lengths due to the open-ended class II binding site. The method has applications in both the accurate prediction of class II epitopes and the manipulation of affinity for heteroclitic and competitor peptides. The method is applied here to six class II mouse alleles (I-Ab, I-Ad, I-Ak, I-As, I-Ed, and I-Ek) and included peptides up to 25 amino acids in length. A series of regression equations highlighting the quantitative contributions of individual amino acids at each peptide position was established. The initial model for each allele exhibited only moderate predictivity. Once the set of selected peptide subsequences had converged, the final models exhibited a satisfactory predictive power. Convergence was reached between the 4th and 17th iterations, and the leave-one-out cross-validation statistical terms - q2, SEP, and NC - ranged between 0.732 and 0.925, 0.418 and 0.816, and 1 and 6, respectively. The non-cross-validated statistical terms r2 and SEE ranged between 0.98 and 0.995 and 0.089 and 0.180, respectively. The peptides used in this study are available from the AntiJen database (http://www.jenner.ac.uk/AntiJen). The PLS method is available commercially in the SYBYL molecular modeling software package. The resulting models, which can be used for accurate T-cell epitope prediction, will be made freely available online (http://www.jenner.ac.uk/MHCPred).