Direct molecular level measurements of the electrostatic properties of a protein surface

In this work, we used direct measurements with the surface force apparatus to determine the pH-dependent electrostatic charge density of a single binding face of streptavidin. Mean field calculations have been used with considerable success to model electrostatic potential fields near protein surfaces, but these models and their inherent assumptions have not been tested directly at the molecular level. Using the force apparatus and immobilized, oriented monolayers of streptavidin, we measured a pI of 5–5.5 for the biotin-binding face of the protein. This differs from the pI of 6.3 for the soluble protein and confirms that we probed the local electrostatic features of the macromolecule. With finite difference solutions of the linearized Poisson–Boltzmann equation, we then calculated the pH-dependent charge densities adjacent to the same face of the protein. These calculated values agreed quantitatively with those obtained by direct force measurements. Although our study focuses on the pH-dependence of surface electrostatics, this direct approach to probing the electrostatic features of proteins is applicable to investigations of any perturbations that alter the charge distribution of the surfaces of immobilized molecules.

Cell surface expression of the epithelial Na channel and a mutant causing Liddle syndrome: A quantitative approach

The epithelial amiloride-sensitive sodium channel (ENaC) controls transepithelial Na+ movement in Na+-transporting epithelia and is associated with Liddle syndrome, an autosomal dominant form of salt-sensitive hypertension. Detailed analysis of ENaC channel properties and the functional consequences of mutations causing Liddle syndrome has been, so far, limited by lack of a method allowing specific and quantitative detection of cell-surface-expressed ENaC. We have developed a quantitative assay based on the binding of 125I-labeled M2 anti-FLAG monoclonal antibody (M2Ab*) directed against a FLAG reporter epitope introduced in the extracellular loop of each of the α, β, and γ ENaC subunits. Insertion of the FLAG epitope into ENaC sequences did not change its functional and pharmacological properties. The binding specificity and affinity (Kd = 3 nM) allowed us to correlate in individual Xenopus oocytes the macroscopic amiloride-sensitive sodium current (INa) with the number of ENaC wild-type and mutant subunits expressed at the cell surface. These experiments demonstrate that: (i) only heteromultimeric channels made of α, β, and γ ENaC subunits are maximally and efficiently expressed at the cell surface; (ii) the overall ENaC open probability is one order of magnitude lower than previously observed in single-channel recordings; (iii) the mutation causing Liddle syndrome (β R564stop) enhances channel activity by two mechanisms, i.e., by increasing ENaC cell surface expression and by changing channel open probability. This quantitative approach provides new insights on the molecular mechanisms underlying one form of salt-sensitive hypertension.

Molecular computation: RNA solutions to chess problems

We have expanded the field of “DNA computers” to RNA and present a general approach for the solution of satisfiability problems. As an example, we consider a variant of the “Knight problem,” which asks generally what configurations of knights can one place on an n × n chess board such that no knight is attacking any other knight on the board. Using specific ribonuclease digestion to manipulate strands of a 10-bit binary RNA library, we developed a molecular algorithm and applied it to a 3 × 3 chessboard as a 9-bit instance of this problem. Here, the nine spaces on the board correspond to nine “bits” or placeholders in a combinatorial RNA library. We recovered a set of “winning” molecules that describe solutions to this problem.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Molecular cloning of a cytochrome P450 taxane 10β-hydroxylase cDNA from Taxus and functional expression in yeast

The early steps in the biosynthesis of Taxol involve the cyclization of geranylgeranyl diphosphate to taxa-4(5),11(12)-diene followed by cytochrome P450-mediated hydroxylation at C5, acetylation of this intermediate, and a second cytochrome P450-dependent hydroxylation at C10 to yield taxadien-5α-acetoxy-10β-ol. Subsequent steps of the pathway involve additional cytochrome P450 catalyzed oxygenations and CoA-dependent acylations. The limited feasibility of reverse genetic cloning of cytochrome P450 oxygenases led to the use of Taxus cell cultures induced for Taxol production and the development of an approach based on differential display of mRNA-reverse transcription-PCR, which ultimately provided full-length forms of 13 unique but closely related cytochrome P450 sequences. Functional expression of these enzymes in yeast was monitored by in situ spectrophotometry coupled to in vivo screening of oxygenase activity by feeding taxoid substrates. This strategy yielded a family of taxoid-metabolizing enzymes and revealed the taxane 10β-hydroxylase as a 1494-bp cDNA that encodes a 498-residue cytochrome P450 capable of transforming taxadienyl acetate to the 10β-hydroxy derivative; the identity of this latter pathway intermediate was confirmed by chromatographic and spectrometric means. The 10β-hydroxylase represents the initial cytochrome P450 gene of Taxol biosynthesis to be isolated by an approach that should provide access to the remaining oxygenases of the pathway.

Integrated fossil and molecular data reconstruct bat echolocation

Molecular and morphological data have important roles in illuminating evolutionary history. DNA data often yield well resolved phylogenies for living taxa, but are generally unattainable for fossils. A distinct advantage of morphology is that some types of morphological data may be collected for extinct and extant taxa. Fossils provide a unique window on evolutionary history and may preserve combinations of primitive and derived characters that are not found in extant taxa. Given their unique character complexes, fossils are critical in documenting sequences of character transformation over geologic time and may elucidate otherwise ambiguous patterns of evolution that are not revealed by molecular data alone. Here, we employ a methodological approach that allows for the integration of molecular and paleontological data in deciphering one of the most innovative features in the evolutionary history of mammals—laryngeal echolocation in bats. Molecular data alone, including an expanded data set that includes new sequences for the A2AB gene, suggest that microbats are paraphyletic but do not resolve whether laryngeal echolocation evolved independently in different microbat lineages or evolved in the common ancestor of bats and was subsequently lost in megabats. When scaffolds from molecular phylogenies are incorporated into parsimony analyses of morphological characters, including morphological characters for the Eocene taxa Icaronycteris, Archaeonycteris, Hassianycteris, and Palaeochiropteryx, the resulting trees suggest that laryngeal echolocation evolved in the common ancestor of fossil and extant bats and was subsequently lost in megabats. Molecular dating suggests that crown-group bats last shared a common ancestor 52 to 54 million years ago.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

New approach for inhibiting Rev function and HIV-1 production using the influenza virus NS1 protein.

The Rev protein of HIV-1, which facilitates the nuclear export of HIV-1 pre-mRNAs, has been a target for antiviral therapy. Here we describe a new strategy for inhibiting Rev function and HIV-1 replication. In contrast to previous approaches, we use a wild-type rather than a mutant Rev protein and covalently link this Rev sequence to the NS1 protein of influenza A virus, a protein that inhibits the nuclear export of mRNAs. The NS1 protein contains an RNA-binding domain mutation (RM), so that the only functional RNA-binding domain in the chimeric protein (NS1RM-Rev) is in the Rev protein sequence. In the presence of the NS1RM-Rev chimeric protein, HIV-1 pre-mRNAs were retained in, rather than exported from, the nucleus. In addition, this chimeric protein effectively inhibited Rev function in trans in transfection experiments and effectively inhibited the production of HIV-1 in tissue culture cells transfected with an infectious molecular clone of HIV-1 DNA. The inhibitory activities of the NS1RM-Rev chimera were at least equivalent to those of the Rev M10 mutant protein, which has been considered to be the prototype trans inhibitor of Rev function and is currently in phase I clinical trials for the treatment of AIDS patients. We discuss (i) the potential for increasing the inhibitory activity of NS1-Rev chimeras against HIV-1 and (ii) the need for additional studies to evaluate these chimeras for the treatment of AIDS.

Human immunodeficiency virus (HIV) type 2-mediated inhibition of HIV type 1: a new approach to gene therapy of HIV-infection.

Human immunodeficiency virus (HIV) type 2, the second AIDS-associated human retrovirus, differs from HIV-1 in its natural history, infectivity, and pathogenicity, as well as in details of its genomic structure and molecular behavior. We report here that HIV-2 inhibits the replication of HIV-1 at the molecular level. This inhibition was selective, dose-dependent, and nonreciprocal. The closely related simian immunodeficiency provirus also inhibited HIV-1. The selectivity of inhibition was shown by the observation that HIV-2 did not significantly downmodulate the expression of the unrelated murine leukemia virus; neither did the murine leukemia virus markedly affect HIV-1 or HIV-2 expression. Moreover, while HIV-2 potently inhibited HIV-1, the reverse did not happen, thus identifying yet another and remarkable difference between HIV-1 and HIV-2. Mutational analysis of the HIV-2 genome suggested that the inhibition follows a complex pathway, possibly involving multiple genes and redundant mechanisms. Introduction of inactivating mutations into the structural and regulatory/accessory genes did not render the HIV-2 provirus ineffective. Some of the HIV-2 gene defects, such as that of tat and rev genes, were phenotypically transcomplemented by HIV-1. The HIV-2 proviruses with deletions in the putative packaging signal and defective for virus replication were effective in inducing the suppressive phenotype. Though the exact mechanism remains to be defined, the inhibition appeared to be mainly due to an intracellular molecular event because it could not be explained solely on the basis of cell surface receptor mediated interference. The results support the notion that the inhibition likely occurred at the level of viral RNA, possibly involving competition between viral RNAs for some transcriptional factor essential for virus replication. Induction of a cytokine is another possibility. These findings might be relevant to the clinical-epidemiological data suggesting that infection with HIV-2 may offer some protection against HIV-1 infection.

Molecular biology of retinal ganglion cells.

Retinal ganglion cells are the output neurons that encode and transmit information from the eye to the brain. Their diverse physiologic and anatomic properties have been intensively studied and appear to account well for a number of psychophysical phenomena such as lateral inhibition and chromatic opponency. In this paper, we summarize our current view of retinal ganglion cell properties and pose a number of questions regarding underlying molecular mechanisms. As an example of one approach to understanding molecular mechanisms, we describe recent work on several POU domain transcription factors that are expressed in subsets of retinal ganglion cells and that appear to be involved in ganglion cell development.

Molecular organization of histidine-tagged biomolecules at self-assembled lipid interfaces using a novel class of chelator lipids.

In molecular biology, the expression of fusion proteins is a very useful and well-established technique for the identification and one-step purification of gene products. Even a short fused sequence of five or six histidines enables proteins to bind to an immobilized metal ion chelate complex. By synthesis of a class of chelator lipids, we have transferred this approach to the concept of self-assembly. The specific interaction and lateral organization of a fluorescent fusion molecule containing a C-terminal oligohistidine sequence was studied by film balance techniques in combination with epifluorescence microscopy. Due to the phase behavior of the various lipid mixtures used, the chelator lipids can be laterally structured, generating two-dimensional arrays of histidine-tagged biomolecules. Because of the large variety of fusion proteins already available, this concept represents a powerful technique for orientation and organization of proteins at lipid interfaces with applications in biosensing, biofunctionalization of nanostructured interfaces, two-dimensional crystallization, and studies of lipid-anchored proteins.

Interfacial activation-based molecular bioimprinting of lipolytic enzymes.

Interfacial activation-based molecular (bio)-imprinting (IAMI) has been developed to rationally improve the performance of lipolytic enzymes in nonaqueous environments. The strategy combinedly exploits (i) the known dramatic enhancement of the protein conformational rigidity in a water-restricted milieu and (ii) the reported conformational changes associated with the activation of these enzymes at lipid-water interfaces, which basically involves an increased substrate accessibility to the active site and/or an induction of a more competent catalytic machinery. Six model enzymes have been assayed in several model reactions in nonaqueous media. The results, rationalized in light of the present biochemical and structural knowledge, show that the IAMI approach represents a straightforward, versatile method to generate manageable, activated (kinetically trapped) forms of lipolytic enzymes, providing under optimal conditions nonaqueous rate enhancements of up to two orders of magnitude. It is also shown that imprintability of lipolytic enzymes depends not only on the nature of the enzyme but also on the "quality" of the interface used as the template.

Nutraceutical approach to the modulation of cholesterol metabolism: in vitro and in vivo studies

Negli ultimi anni, si sono diffusi nuove strategie per il trattamento delle malattie cardiovascolari, che possano supportare una terapia medica, o in alcuni casi, sostituirla. Infatti, l’abbandono delle terapie è il più importante problema di salute pubblica del mondo occidentale, soprattutto per le malattie croniche. Ciò è dovuto alla complessità delle terapie farmacologiche e ai numerosi e in alcuni casi gravi effetti collaterali dei farmaci somministrati. Di conseguenza, una riduzione di questi effetti migliorerebbe le condizioni di vita del paziente e quindi diminuirebbe il rischio di abbandono della terapia. Per ottenere ciò, è possibile affiancare al trattamento farmacologico una terapia nutraceutica, consistente nella somministrazione di complessi molecolari o microorganismi, provenienti da piante, latte o cibi funzionali. Lo scopo generale di questo studio è indagare le attività ipolipidemizzanti di un composto nutraceutico e di un ceppo batterio specifico nel modello animale che presenta elevati alti livelli plasmatici di colesterolo. Inoltre, sono stati analizzati gli effetti del trattamento nutraceutico sui meccanismi fisiologici che contrastano la creazione della placca aterosclerotica come l’efflusso di colesterolo dalle “foam cells” presenti nell’ateroma, o la riduzione dell’assorbimento intestinale di colesterolo. La presente tesi è divisa in due parti. Nella prima parte, abbiamo analizzato la capacità dei Bifidobacteria di ridurre i livelli di colesterolo nel medium di crescita. Dall’analisi, si è osservato che vari ceppi del genere Bifidobacteria presentano un’ampia capacità di assimilazione del colesterolo all’interno della cellula batterica, in particolare il Bifidobacterium bifidum PRL2010. Le analisi di trascrittomica del Bb PRL2010 incubato in presenza di colesterolo, hanno rivelato un significativo aumento dei livelli di trascrizione di geni codificanti trasportatori e riduttasi, responsabili del meccanismo di accumulo all’interno della cellula batterica e della conversione del colesterolo in coprostanolo. L’attività ipolipidemizzante del Bb PRL2010 è stata poi valutata nel modello murino, mostrando la modificazione del microbiota dei topi trattati dopo somministrazione del batterio in questione. Nella seconda parte del progetto di ricerca, abbiamo indagato sugli effetti di un composto coperto da brevetto, chiamato “Ola”, sull’efflusso di colesterolo di criceti trattati con questo composto nutraceutico. L’efflusso di colesterolo è il primo step del meccanismo fisiologico noto come Trasporto Inverso del Colesterolo, che consente l’eliminazione del colesterolo dalle placche aterosclerotiche, attraverso l’interazione fra le HDL, presenti nella circolazione sanguigna, e specifici trasportatori delle foam cells, come ABCA1/G1 e SR-BI. In seguito, le lipoproteine rilasciano il colesterolo alle cellule epatiche, dove è metabolizzato ed escreto attraverso le feci. Per valutare l’effetto dell’Ola sul profilo lipidico dei criceti, sono state condotte analisi in vitro. I risultati mostrano un aumento dell’efflusso di colesterolo in cellule che esprimono il trasportatore ABCA1, comparato con il gruppo controllo. Questi due studi mostrano come l’approccio nutraceutico può essere un importante modo per contrastare l’aterosclerosi. Come mostrato in letteratura, gli effetti dei composti nutraceutici sull’aterosclerosi e su altre malattie croniche, hanno portato a un ampio uso come supporto alle terapie farmacologiche, ed in alcuni casi hanno rimpiazzato la terapia farmacologica stessa.

Intermolecular and Intramolecular Charge Transfer in Functional Molecular Materials

This thesis is devoted to the investigation of inter and intramolecular charge transfer (CT) in molecular functional materials and specifically organic dyes and CT crystals. An integrated approach encompassing quantum-chemical calculations, semiempirical tools, theoretical models and spectroscopic measurements is applied to understand structure-property relationships governing the low-energy physics of these materials. Four main topics were addressed: 1) Spectral properties of organic dyes. Charge-transfer dyes are constituted by electron donor (D) and electron acceptor (A) units linked through bridge(s) to form molecules with different symmetry and dimensionality. Their low-energy physics is governed by the charge resonance between D and A groups and is effectively described by a family of parametric Hamiltonians known as essential-state models. These models account for few electronic states, corresponding to the main resonance structures of the relevant dye, leading to a simple picture that is completed introducing the coupling of the electronic system to molecular vibrations, treated in a non-adiabatic way, and an effective classical coordinate, describing polar solvation. In this work a specific essential-state model was proposed and parametrized for the dye Brilliant Green. The central issue in this work has been the definition of the diabatic states, a not trivial task for a multi-branched chromophore. In a second effort, we have used essential-state models for the description of the early-stage dynamics of excited states after ultrafast excitation. Crucial to this work is the fully non-adiabatic treatment of the coupled electronic and vibrational motion, allowing for a reliable description of the dynamics of systems showing a multistable, broken-symmetry excited state. 2) Mixed-stack CT salts. Mixed-stack (MS) CT crystals are an interesting class of multifunctional molecular materials, where D and A molecules arrange themselves to form stacks, leading to delocalized electrons in one dimension. The interplay between the intermolecular CT, electrostatic interactions, lattice phonons and molecular vibrations leads to intriguing physical properties that include (photoinduced) phase transitions, multistability, antiferromagnetism, ferroelectricity and potential multiferroicity. The standard microscopic model to describe this family of materials is the Modified Hubbard model accounting for electron-phonon coupling (Peierls coupling), electron-molecular vibrations coupling (Holstein coupling) and electrostatic interactions. We adopt and validate a method, based on DFT calculations on dimeric DA structures, to extract relevant model parameters. The approach offers a powerful tool to shed light on the complex physics of MS-CT salts. 3) Charge transfer in organic radical dipolar dyes. In collaboration with the group of Prof. Jaume Veciana (ICMAB- Barcellona), we have studied spectral properties of a special class of CT dyes with D-bridge-A structure where the acceptor group is a stable radical (of the perchlorotriphenylmethyl, PTM, family), leading to an open-shell CT dyes. These materials are of interest since they associate the electronic and optical properties of CT dyes with magnetic properties from the unpaired electron. The first effort was devoted to the parametrization of the relevant essential-state model. Two strategies were adopted, one based on the calculation of the low-energy spectral properties, the other based on the variation of ground state properties with an applied electric field. 4) The spectral properties of organic nanoparticles based on radical species are investigated in collaboration with Dr. I. Ratera (ICMAB- Barcellona). Intriguing spectroscopic behavior was observed pointing to the presence of excimer states. In an attempt to rationalize these findings, extensive calculations (TD-DFT and ZINDO) were performed. The results for the isolated dyes are validated against experimental spectra in solution. To address intermolecular interactions we studied dimeric structures in the gas phase, but the preliminary results obtained do not support excimer formation.