Risk assessment of recurrence in sporadic retinoblastoma using a molecular-based algorithm.

PURPOSE: Most RB1 mutations are unique and distributed throughout the RB1 gene. Their detection can be time-consuming and the yield especially low in cases of conservatively-treated sporadic unilateral retinoblastoma (Rb) patients. In order to identify patients with true risk of developing Rb, and to reduce the number of unnecessary examinations under anesthesia in all other cases, we developed a universal sensitive, efficient and cost-effective strategy based on intragenic haplotype analysis. METHODS: This algorithm allows the calculation of the a posteriori risk of developing Rb and takes into account (a) RB1 loss of heterozygosity in tumors, (b) preferential paternal origin of new germline mutations, (c) a priori risk derived from empirical data by Vogel, and (d) disease penetrance of 90% in most cases. We report the occurrence of Rb in first degree relatives of patients with sporadic Rb who visited the Jules Gonin Eye Hospital, Lausanne, Switzerland, from January 1994 to December 2006 compared to expected new cases of Rb using our algorithm. RESULTS: A total of 134 families with sporadic Rb were enrolled; testing was performed in 570 individuals and 99 patients younger than 4 years old were identified. We observed one new case of Rb. Using our algorithm, the cumulated total a posteriori risk of recurrence was 1.77. CONCLUSIONS: This is the first time that linkage analysis has been validated to monitor the risk of recurrence in sporadic Rb. This should be a useful tool in genetic counseling, especially when direct RB1 screening for mutations leaves a negative result or is unavailable.

Molecular phylogenetics of shrews (Mammalia: Soricidae) reveal timing of transcontinental colonizations.

We sequenced 2167 base pairs (bp) of mitochondrial DNA cytochrome b and 16S, and 1390 bp of nuclear genes BRCA1 and ApoB in shrews taxa (Eulipotyphla, family Soricidae). The aim was to study the relationships at higher taxonomic levels within this family, and in particular the position of difficult clades such as Anourosorex and Myosorex. The data confirmed two monophyletic subfamilies, Soricinae and Crocidurinae. In the former, the tribes Anourosoricini, Blarinini, Nectogalini, Notiosoricini, and Soricini were supported. The latter was formed by the tribes Myosoricini and Crocidurini. The genus Suncus appeared to be paraphyletic and included Sylvisorex. We further suggest a biogeographical hypothesis, which shows that North America was colonized by three independent lineages of Soricinae during middle Miocene. Our hypothesis is congruent with the first fossil records for these taxa. Using molecular dating, the first exchanges between Africa and Eurasia occurred during the middle Miocene. The last one took place in the Late Miocene, with the dispersion of the genus Crocidura through the old world.

Exome sequencing of index patients with retinal dystrophies as a tool for molecular diagnosis.

BACKGROUND: Retinal dystrophies (RD) are a group of hereditary diseases that lead to debilitating visual impairment and are usually transmitted as a Mendelian trait. Pathogenic mutations can occur in any of the 100 or more disease genes identified so far, making molecular diagnosis a rather laborious process. In this work we explored the use of whole exome sequencing (WES) as a tool for identification of RD mutations, with the aim of assessing its applicability in a diagnostic context. METHODOLOGY/PRINCIPAL FINDINGS: We ascertained 12 Spanish families with seemingly recessive RD. All of the index patients underwent mutational pre-screening by chip-based sequence hybridization and resulted to be negative for known RD mutations. With the exception of one pedigree, to simulate a standard diagnostic scenario we processed by WES only the DNA from the index patient of each family, followed by in silico data analysis. We successfully identified causative mutations in patients from 10 different families, which were later verified by Sanger sequencing and co-segregation analyses. Specifically, we detected pathogenic DNA variants (∼50% novel mutations) in the genes RP1, USH2A, CNGB3, NMNAT1, CHM, and ABCA4, responsible for retinitis pigmentosa, Usher syndrome, achromatopsia, Leber congenital amaurosis, choroideremia, or recessive Stargardt/cone-rod dystrophy cases. CONCLUSIONS/SIGNIFICANCE: Despite the absence of genetic information from other family members that could help excluding nonpathogenic DNA variants, we could detect causative mutations in a variety of genes known to represent a wide spectrum of clinical phenotypes in 83% of the patients analyzed. Considering the constant drop in costs for human exome sequencing and the relative simplicity of the analyses made, this technique could represent a valuable tool for molecular diagnostics or genetic research, even in cases for which no genotypes from family members are available.

Molecular characterization of Unc5CL/ZUD

Résumé : Au cours de l'évolution, les organismes multicellulaires ont développé le système immunitaire afin de pouvoir se défendre contre les pathogènes tel que les bactéries, les virus, et les parasites. La réponse immunitaire doit être finement régulée par différentes voies de signalisation moléculaire, afin d'assurer une efficacité optimale, et d'éviter des dommages tissulaires indésirables. Les résultats expérimentaux décrits dans ce manuscrit, mettent en évidence que la protéine Unc5CL, qui contient un death domain (DD), est impliquée dans la régulation de la réponse immunitaire des muqueuses. Il a été démontré que cette protéine contient aussi un domaine transmembranaire de type III dans sa partie N-terminale, permettant ainsi de l'ancrer et d'exposer sa partie C-terminale dans le cytosol, un prérequis pour la signalisation dans ce compartiment cellulaire. De plus, cette protéine a la capacité d'activer le facteur de transcription NFxB, qui joue un rôle important dans le système immunitaire, ainsi que dans d'autres processus cellulaires essentiels. Le profil transcriptionnel révèle que l'activation de NF-κB induite par Unc5CL conduit principalement à une réponse inflammatoire, qui se caractérise par la production de diverses chimiokines (e.g. CXCL-1, IL-8 et CCL20). Il a également été démontré que Unc5CL requiert les mêmes molécules qui sont utilisées dans la voie de signalisation des récepteurs de la famille toll et de l'interleukine-1. De manière similaire à leur protéine adaptatrice MyD88, Unc5CL a la capacité de recruter, via une interaction homotypique DD-DD, les kinases IRAK1 et IRAK4 qui contiennent elles aussi un DD, permettant ainsi au signal d'être transmis. La production d'un anticorps polyclonal contre le DD de Unc5CL a permis d'identifier des lignées cellulaires et des tissus exprimant cette protéine, ainsi que de déterminer sa localisation sub-cellulaire. Unc5CL a été détecté dans les cellules de la muqueuse utérine et intestinale, ainsi que dans une lignée cellulaire issue d'un adénocarcinome colorectal humain, les CaCo-2. Dans chacun de ces cas, Unc5CL a été principalement détectée au niveau apical des cellules épithéliales polarisées. De manière similaire à PIDD, une protéine impliquée dans la réponse aux dommages à l'ADN, et au constituant des pores nucléaires Nup98, Unc5CL est constitutivement clivé de manière autoprotéolytique, au niveau d'un site HFS. Il est intéressant d'observer que les deux fragments ainsi générés restent fortement associés l'un à l'autre après clivage. Finalement, un criblage protéomique pour identifier un partenaire d'interaction, a mis en évidence l'ubiquitin ligase E3 ITCH, qui régule de manière négative Unc5CL en augmentant sa dégradation. Summary : Multicellular organisms have evolved the immune system in order to defend themselves against pathogens such as bacteria, viruses and eukaryotic parasites. Immune responses have to be tightly orchestrated by signaling mechanisms to achieve optimal effectiveness and minimal tissue damage. The experimental results in this thesis manuscript provide evidence that the death domain (DD)-containing protein Unc5CL might be involved in the regulation of mucosal immune responses. It could be shown that the protein contains an N-terminal type-III transmembrane domain that anchors the protein with its C-terminus exposed to the cytosol, a prerequisite for signaling events in this compartment. Furthermore, the protein has the capacity to activate the transcription factor NF-κB, which plays an important role in the immune system as well as in other essential cellular processes. Transcriptional profiling revealed that Unc5CL-mediated activation of NF-κB mainly leads to an inflammatory response, characterized by the production of chemokines (e.g. CXCL-l, IL-8 and CCL20). Furthermore, it could be shown that Unc5CL requires the same downstream signaling molecules as the evolutionarily ancient tolUinterleukin-1 receptor family. Similar to their adapter protein MyD88, Unc5CL has the capacity to recruit the DD-containing kinases IRAKI and IRAK4 for signaling and can interact with these proteins via homotypic DD-DD interactions. Generation of polyclonal antibodies raised against the DD of Unc5CL allowed the identification of cell lines and tissues that express the endogenous protein as well as to confine its subcellular localization. Unc5CL was detected in primary mucosal uterine and intestinal epithelial cells as well as in the human colorectal adenocarcinoma cell line CaCo-2. In all cases, the protein was mainly localized to the apical face of these polarized epithelial cells. Similar to PIDD, a protein critically involved in responses to DNA damage, and the nuclear pore component Nup98, Unc5CL is constitutively autoproteolytically processed at an HFS site. Interestingly, the two generated cleavage fragments remain tightly associated after processing. Finally, a proteomics screen for interaction partners identified the E3 ubiquitin ligase ITCH as a negative regulator of Unc5CL by targeting the protein for degradation.

A comparative analysis of two methods for the calculation of electric-field-induced perturbations to molecular vibration

Two common methods of accounting for electric-field-induced perturbations to molecular vibration are analyzed and compared. The first method is based on a perturbation-theoretic treatment and the second on a finite-field treatment. The relationship between the two, which is not immediately apparent, is made by developing an algebraic formalism for the latter. Some of the higher-order terms in this development are documented here for the first time. As well as considering vibrational dipole polarizabilities and hyperpolarizabilities, we also make mention of the vibrational Stark effec

Monoclonal antibody against carcinoembryonic antigen (CEA) identifies two new forms of crossreacting antigens of molecular weight 90,000 and 160,000 in normal granulocytes.

Two new forms of non-specific crossreacting antigens (NCAs) were identified in the Nonidet P40 (NP-40) extracts of normal granulocytes by precipitation with the monoclonal antibody (MAb) 192 directed against carcinoembryonic antigen (CEA) and already known to crossreact with the perchloric acid soluble NCA-55. The NP-40 soluble NCAs recognized by MAb 192 have apparent mol. wts of 90,000 and 160,000 in sodium dodecyl sulfate-polyacrylamide gel electrophoresis (SDS-PAGE). Both NCAs appear to consist of a single monomeric polypeptide chain, since they have the same electrophoretic mobility in SDS-PAGE under reduced and non-reduced conditions. When granulocytes were extracted with perchloric acid instead of NP-40, only the 55,000 mol. wt antigen, corresponding to the previously described NCA-55, was precipitated by MAb 192. Furthermore, it was shown that NCA-55 is not a degradation product of NCA-90 or NCA-160 due to the perchloric acid treatment because exposure to perchloric acid of NCA preparations purified from NP-40 extracts did not change their apparent mol. wts in SDS-PAGE. It was also shown that NCA-160 is not a granulocytic form of CEA because it was not precipitated by the MAb 35 reacting exclusively with CEA. Immunocytochemical studies of granulocytes and macrophages showed that MAb 192 stained both types of cells whereas MAb 47 stained only the granulocytes and MAb 35 none of these cells. In granulocytes both MAbs reacted with antigens associated with granules and also present at the periphery of the nucleus as well as in the Golgi apparatus. The NCA-90 identified by MAb 192 was found by sequential immunodepletion to be antigenically distinct from the NCA-95 precipitated by MAb 47. The epitope recognized by MAb 192 on CEA and NCA molecules appears to be on the peptidic moiety because the antigens deglycosylated by the enzyme Endo F were still precipitated by this MAb. Taken together, the results indicate that MAb 192 identifies two novel forms of NCA (NCA-90 and NCA-160) in NP-40 extracts of granulocytes, which are distinct from CEA and the previously described NCA-55 and NCA-95 identified by MAbs 192 and 47, respectively, in perchloric acid extracts of granulocytes.

Effect of basis set superposition error on the electron density of molecular complexes

The effect of basis set superposition error (BSSE) on molecular complexes is analyzed. The BSSE causes artificial delocalizations which modify the first order electron density. The mechanism of this effect is assessed for the hydrogen fluoride dimer with several basis sets. The BSSE-corrected first-order electron density is obtained using the chemical Hamiltonian approach versions of the Roothaan and Kohn-Sham equations. The corrected densities are compared to uncorrected densities based on the charge density critical points. Contour difference maps between BSSE-corrected and uncorrected densities on the molecular plane are also plotted to gain insight into the effects of BSSE correction on the electron density

A quantum molecular similarity analysis of changes in molecular electron density caused by basis set flotation and electric field application

Quantum molecular similarity (QMS) techniques are used to assess the response of the electron density of various small molecules to application of a static, uniform electric field. Likewise, QMS is used to analyze the changes in electron density generated by the process of floating a basis set. The results obtained show an interrelation between the floating process, the optimum geometry, and the presence of an external field. Cases involving the Le Chatelier principle are discussed, and an insight on the changes of bond critical point properties, self-similarity values and density differences is performed

A comparative analysis by means of quantum molecular similarity measures of density distributions derived from conventional ab initio and density functional methods

A procedure based on quantum molecular similarity measures (QMSM) has been used to compare electron densities obtained from conventional ab initio and density functional methodologies at their respective optimized geometries. This method has been applied to a series of small molecules which have experimentally known properties and molecular bonds of diverse degrees of ionicity and covalency. Results show that in most cases the electron densities obtained from density functional methodologies are of a similar quality than post-Hartree-Fock generalized densities. For molecules where Hartree-Fock methodology yields erroneous results, the density functional methodology is shown to yield usually more accurate densities than those provided by the second order Møller-Plesset perturbation theory

Electronic couplings and on-site energies for hole transfer in DNA: systematic quantum mechanical/molecular dynamic study

The electron hole transfer (HT) properties of DNA are substantially affected by thermal fluctuations of the π stack structure. Depending on the mutual position of neighboring nucleobases, electronic coupling V may change by several orders of magnitude. In the present paper, we report the results of systematic QM/molecular dynamic (MD) calculations of the electronic couplings and on-site energies for the hole transfer. Based on 15 ns MD trajectories for several DNA oligomers, we calculate the average coupling squares 〈 V2 〉 and the energies of basepair triplets X G+ Y and X A+ Y, where X, Y=G, A, T, and C. For each of the 32 systems, 15 000 conformations separated by 1 ps are considered. The three-state generalized Mulliken-Hush method is used to derive electronic couplings for HT between neighboring basepairs. The adiabatic energies and dipole moment matrix elements are computed within the INDO/S method. We compare the rms values of V with the couplings estimated for the idealized B -DNA structure and show that in several important cases the couplings calculated for the idealized B -DNA structure are considerably underestimated. The rms values for intrastrand couplings G-G, A-A, G-A, and A-G are found to be similar, ∼0.07 eV, while the interstrand couplings are quite different. The energies of hole states G+ and A+ in the stack depend on the nature of the neighboring pairs. The X G+ Y are by 0.5 eV more stable than X A+ Y. The thermal fluctuations of the DNA structure facilitate the HT process from guanine to adenine. The tabulated couplings and on-site energies can be used as reference parameters in theoretical and computational studies of HT processes in DNA

Two complementary molecular energy decomposition schemes: the Mayer and Ziegler-Rauk methods in comparison

In the present paper we discuss and compare two different energy decomposition schemes: Mayer's Hartree-Fock energy decomposition into diatomic and monoatomic contributions [Chem. Phys. Lett. 382, 265 (2003)], and the Ziegler-Rauk dissociation energy decomposition [Inorg. Chem. 18, 1558 (1979)]. The Ziegler-Rauk scheme is based on a separation of a molecule into fragments, while Mayer's scheme can be used in the cases where a fragmentation of the system in clearly separable parts is not possible. In the Mayer scheme, the density of a free atom is deformed to give the one-atom Mulliken density that subsequently interacts to give rise to the diatomic interaction energy. We give a detailed analysis of the diatomic energy contributions in the Mayer scheme and a close look onto the one-atom Mulliken densities. The Mulliken density ρA has a single large maximum around the nuclear position of the atom A, but exhibits slightly negative values in the vicinity of neighboring atoms. The main connecting point between both analysis schemes is the electrostatic energy. Both decomposition schemes utilize the same electrostatic energy expression, but differ in how fragment densities are defined. In the Mayer scheme, the electrostatic component originates from the interaction of the Mulliken densities, while in the Ziegler-Rauk scheme, the undisturbed fragment densities interact. The values of the electrostatic energy resulting from the two schemes differ significantly but typically have the same order of magnitude. Both methods are useful and complementary since Mayer's decomposition focuses on the energy of the finally formed molecule, whereas the Ziegler-Rauk scheme describes the bond formation starting from undeformed fragment densities

Molecular chaperones and associated cellular clearance mechanisms against toxic protein conformers in Parkinson's disease.

Parkinson's disease (PD) is a slowly progressive neurodegenerative disorder marked by the loss of dopaminergic neurons (in particular in the substantia nigra) causing severe impairment of movement coordination and locomotion, associated with the accumulation of aggregated α-synuclein (α-Syn) into proteinaceous inclusions named Lewy bodies. Various early forms of misfolded α-Syn oligomers are cytotoxic. Their formation is favored by mutations and external factors, such as heavy metals, pesticides, trauma-related oxidative stress and heat shock. Here, we discuss the role of several complementing cellular defense mechanisms that may counteract PD pathogenesis, especially in youth, and whose effectiveness decreases with age. Particular emphasis is given to the 'holdase' and 'unfoldase' molecular chaperones that provide cells with potent means to neutralize and scavenge toxic protein conformers. Because chaperones can specifically recognize misfolded proteins, they are key specificity factors for other cellular defenses, such as proteolysis by the proteasome and autophagy. The efficiency of the cellular defenses decreases in stressed or aging neurons, leading to neuroinflammation, apoptosis and tissue loss. Thus, drugs that can upregulate the molecular chaperones, the ubiquitin-proteasome system and autophagy in brain tissues are promising avenues for therapies against PD and other mutation-, stress- or age-dependent protein-misfolding diseases.

Multifunctional ligands for application in Alzheimer’s Disease: molecular modelling and biological evaluation

Projecte de recerca elaborat a partir d’una estada a la University of British Columbia, Canadà, entre 2010 i 2012 La malaltia d'Alzheimer (MA) representa avui la forma més comuna de demència en la població envellida. Malgrat fa 100 anys que va ser descoberta, encara avui no existeix cap tractament preventiu i/o curatiu ni cap agent de diagnòstic que permeti valorar quantitativament l'evolució d'aquesta malaltia. L'objectiu en el que s'emmarca aquest treball és contribuir a aportar solucions al problema de la manca d'agents terapèutics i de diagnosi, unívocs i rigorosos, per a la MA. Des del camp de la química bioinorgànica és fàcil fixar-se en l'excessiva concentració d'ions Zn(II) i Cu(II) en els cervells de malalts de MA, plantejar-se la seva utilització com a dianes terapèutica i, en conseqüència, cercar agents quelants que evitin la formació de plaques senils o contribueixin a la seva dissolució. Si bé aquest va ser el punt de partida d’aquest projecte, els múltiples factors implicats en la patogènesi de la MA fan que el clàssic paradigma d’ ¨una molècula, una diana¨ limiti la capacitat de la molècula de combatre aquesta malaltia tan complexa. Per tant, un esforç considerable s’ha dedicat al disseny d’agentsmultifuncionals que combatin els múltiples factors que caracteritzen el desenvolupament de la MA. En el present treball s’han dissenyat agents multifuncionals inspirats en dos esquelets moleculars ben establers i coneguts en el camp de la química medicinal: la tioflavina-T (ThT) i la deferiprona (DFP). La utilització de tècniques in silico que inclouen càlculs farmacocinètics i modelatge molecular ha estat un procés cabdal per a l’avaluació dels millors candidats en base als següents requeriments: (a) compliment de determinades propietats farmacocinètiques que estableixin el seu possible ús com a fàrmac (b) hidrofobicitat adequada per travessar la BBB i (c) interacció amb el pèptid Aen solució.

Exploring unfrequent events with accelerated molecular dynamics simulations: from photochemistry to drug design

La meva incorporació al grup de recerca del Prof. McCammon (University of California San Diego) en qualitat d’investigador post doctoral amb una beca Beatriu de Pinós, va tenir lloc el passat 1 de desembre de 2010; on vaig dur a terme les meves tasques de recerca fins al darrer 1 d’abril de 2012. El Prof. McCammon és un referent mundial en l’aplicació de simulacions de dinàmica molecular (MD) en sistemes biològics d’interès humà. La contribució més important del Prof. McCammon en la simulació de sistemes biològics és el desenvolupament del mètode de dinàmiques moleculars accelerades (AMD). Les simulacions MD convencionals, les quals estan limitades a l’escala de temps del nanosegon (~10-9s), no son adients per l’estudi de sistemes biològics rellevants a escales de temps mes llargues (μs, ms...). AMD permet explorar fenòmens moleculars poc freqüents però que son clau per l’enteniment de molts sistemes biològics; fenòmens que no podrien ser observats d’un altre manera. Durant la meva estada a la “University of California San Diego”, vaig treballar en diferent aplicacions de les simulacions AMD, incloent fotoquímica i disseny de fàrmacs per ordinador. Concretament, primer vaig desenvolupar amb èxit una combinació dels mètodes AMD i simulacions Car-Parrinello per millorar l’exploració de camins de desactivació (interseccions còniques) en reaccions químiques fotoactivades. En segon lloc, vaig aplicar tècniques estadístiques (Replica Exchange) amb AMD en la descripció d’interaccions proteïna-lligand. Finalment, vaig dur a terme un estudi de disseny de fàrmacs per ordinador en la proteïna-G Rho (involucrada en el desenvolupament de càncer humà) combinant anàlisis estructurals i simulacions AMD. Els projectes en els quals he participat han estat publicats (o estan encara en procés de revisió) en diferents revistes científiques, i han estat presentats en diferents congressos internacionals. La memòria inclosa a continuació conté més detalls de cada projecte esmentat.