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.

Target Molecular Size and Sodium Dodecyl Sulfate-Polyacrylamide Gel Electrophoresis Analysis of the ATP-and Pyrophosphate-Dependent Proton Pumps from Maize Root Tonoplast.

Tonoplast-enriched membranes were prepared from maize (Zea mays L. cv LG 11) primary roots, using sucrose nonlinear gradients. The functional molecular size of the tonoplast ATP-and PPi-dependent proton pumps were analyzed by radiation inactivation. Glucose-6-phosphate dehydrogenase (G6PDH) was added as an internal standard. Frozen samples (-196 degrees C) of the membranes were irradiated with (60)Co for different periods of time. After thawing the samples, the activities of G6PDH, ATPase, and PPase were tested. By applying target theory, the functional sizes of the ATPase and PPase in situ were found to be around 540 and 160 kilodaltons, respectively. The two activities were solubilized and separated by gel filtration chromatography. The different polypeptides copurifying with the two pumps were analyzed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Two bands (around 59 and 65 kilodaltons) were associated with the ATPase activity, whereas a double band (around 40 kilodaltons) was recovered with the PPase activity.

Molecular mechanism of myosin Va recruitment to dense core secretory granules.

The brain-spliced isoform of Myosin Va (BR-MyoVa) plays an important role in the transport of dense core secretory granules (SGs) to the plasma membrane in hormone and neuropeptide-producing cells. The molecular composition of the protein complex that recruits BR-MyoVa to SGs and regulates its function has not been identified to date. We have identified interaction between SG-associated proteins granuphilin-a/b (Gran-a/b), BR-MyoVa and Rab27a, a member of the Rab family of GTPases. Gran-a/b-BR-MyoVa interaction is direct, involves regions downstream of the Rab27-binding domain, and the C-terminal part of Gran-a determines exon specificity. MyoVa and Gran-a/b are partially colocalised on SGs and disruption of Gran-a/b-BR-MyoVa binding results in a perinuclear accumulation of SGs which augments nutrient-stimulated hormone secretion in pancreatic beta-cells. These results indicate the existence of at least another binding partner of BR-MyoVa that was identified as rabphilin-3A (Rph-3A). BR-MyoVa-Rph-3A interaction is also direct and enhanced when secretion is activated. The BR-MyoVa-Rph-3A and BR-MyoVa-Gran-a/b complexes are linked to a different subset of SGs, and simultaneous inhibition of these complexes nearly completely blocks stimulated hormone release. This study demonstrates that multiple binding partners of BR-MyoVa regulate SG transport, and this molecular mechanism is universally used by neuronal, endocrine and neuroendocrine cells.

Molecular insights for optimizing T cell receptor specificity against cancer.

Cytotoxic CD8 T cells mediate immunity to pathogens and they are able to eliminate malignant cells. Immunity to viruses and bacteria primarily involves CD8 T cells bearing high affinity T cell receptors (TCRs), which are specific to pathogen-derived (non-self) antigens. Given the thorough elimination of high affinity self/tumor-antigen reactive T cells by central and peripheral tolerance mechanisms, anti-cancer immunity mostly depends on TCRs with intermediate-to-low affinity for self-antigens. Because of this, a promising novel therapeutic approach to increase the efficacy of tumor-reactive T cells is to engineer their TCRs, with the aim to enhance their binding kinetics to pMHC complexes, or to directly manipulate the TCR-signaling cascades. Such manipulations require a detailed knowledge on how pMHC-TCR and co-receptors binding kinetics impact the T cell response. In this review, we present the current knowledge in this field. We discuss future challenges in identifying and targeting the molecular mechanisms to enhance the function of natural or TCR-affinity optimized T cells, and we provide perspectives for the development of protective anti-tumor T cell responses.

Identification, using a siRNA screen approach, and molecular characterization of a new Fas pathway modulator: STK11/

Abstract : Apoptosis is an evolutionarily conserved cellular suicide mechanism that can be triggered by activation of various pathways, such as the Fas-Pathway. Upon stimulation by its specific ligand (FasL), present at the surface of Cytotoxic Τ lymphocytes, the death receptor Fas initiates a signaling cascade culminating in the activation of cellular caspases, leading thus to cell death of the target cell (e.g. transformed cell). Dysregulation of apoptosis in general, and of Fas pathway in particular, was shown to contribute to pathogenesis of cancers and many human diseases. Even though, during the last decades the molecular mechanisms of apoptosis have been widely studied, it is important to better understand the mechanisms leading to apoptosis, to improve our understanding of pathological processes, and generate more subtle apoptosis-modulating therapies to fight cancer and other diseases. In order to identify new components of the Fas signaling pathway, a screen based on the mechanism of RNA interference was undertaken. After a first and a second manual whole-kinome screen, we identified several strong positive hits that showed a protection against Fas ligand-induced apoptosis with distinct siRNAs, notably STK11, an interesting tumor suppressor mutated in several sporadic and inherited cancers. The STK11 functional characterization reveals that this kinase represents an apically acting general pro-apoptotic modulator of the extrinsic pathway (FasL, TRAIL, TNF-induced apoptosis), but not of the intrinsic apoptotic pathway. The STK11 action on the Fas pathway was shown to be dependent on its kinase activity, but independent of AMPK, a well-characterized STK11 downstream substrate. Furthermore, STK11 was shown to interact with caspase-8, a major mediator of the extrinsic pathway, and modulate its activity through an unclear mechanism that may involve an STK11-dependant caspase-8 phosphorylation. This modification may allow a proper caspase-8 polyubiquitination and activation in p62 sequestosmes aggregates, but may also increase the activation of caspase-8 at the DISC level. In addition, we observed that STK11 modulate not only the apoptotic pathway induced by Fas engagement, but also FasL-induced JNK and NF- KB, sustaining an upstream role of this kinase in the pathway. In conclusion, our report reveals that STK11 is an important pro-apoptotic modulator of the Fas pathway in particular, and extrinsic pathway in general. Our finding could explain, at least partially, why inactivating mutations of the kinase leads to cancer, by allowing resistance to apoptosis and accordingly evasion of immune surveillance. Résumé : L'apoptose est un mécanisme de suicide cellulaire, conservé dans diverses espèces, et qui au niveau moléculaire est déclenché par différentes voies de signalisation, comme par exemple lors de l'activation du récepteur Fas. La liaison du ligand FasL au récepteur de la mort Fas, induit une cascade de signalisation qui conduit à l'activation des caspases. Les lymphocytes Τ cytotoxiques peuvent utiliser la voie Fas pour induire la mort et se débarrasser de cellules dangereuses pour le reste de l'organisme, tel que les cellules transformées. La dysrégulation de l'apoptose en général, et de la voie Fas en particulier, peut contribuer à diverses maladies telles que le cancer. Même si ces dernières décennies, les mécanismes moléculaires conduisant à l'apoptose ont été extensivement étudiés, il reste néanmoins important de mieux comprendre le phénomène d'apoptose, pour améliorer notre compréhension des processus pathologiques, mais surtout dans le but de développer de nouvelles thérapies ciblant l'apoptose contre le cancer et d'autres pathologies. Pour identifier de nouveau constituants de la voie Fas, un criblage génétique basé sur l'interférence à l'ARN a été entrepris. Après un premier et un deuxième criblage d'une librairie du kinome, nous avons identifié différentes protéines qui pourraient jouer un rôle positif dans la voie Fas, et en particulier la protéine suppresseur de tumeur STK11, qui est fréquemment mutée dans divers cancers sporadiques et héréditaires. La caractérisation fonctionnelle de STK11 a révélé que cette kinase était un modulateur apical de la voie extrinsèque de l'apoptose en général (Fas, TNF, TRAIL), mais pas de la voie intrinsèque. L'action de STK11 sur la voie Fas est dépendante de sa fonction kinase, mais indépendante de l'AMPK, un substrat bien caractérisé de STK11. De plus, STK11 interagît avec la caspase-8, un constituant majeur de la voie Fas, et module son activité, par un mécanisme encore peu clair qui pourrait impliquer une phosphorylation de la caspase-8 par STK11. Cette modification pourrait permettre une activation optimale de la caspase-8 en jouant un rôle dans le processus de polyubiquitination de la caspase-8, phénomène qui semble être important pour l'activation de la caspase-8 dans des agrégats protéiques avec p62, mais qui pourrait aussi augmenter son activation au niveau du DISC. Finalement, nous avons observé que STK11 modulait non seulement la voie apoptotique déclenchée par l'activation de Fas, mais aussi les voies non-apoptotiques de Fas, comme JNK et NF-KB. En conclusion notre étude, révèle que STK11 est un important modulateur pro- apoptotique de la voie Fas, et de la voie extrinsèque en général. Cette découverte pourrait expliquer, du moins partiellement, pourquoi les mutations inactivatrices de STK11 conduisent au cancer, par une augmentation de la résistance à l'apoptose et donc par l'évasion de la surveillance immunitaire.

Cytogenetic and molecular relationships between Zarudny's rock shrew (Crocidura zarudnyi; Mammalia : Soricomorpha) and Eurasian taxa

We karyotyped and sequenced 1,140 base pairs of the mitochondrial DNA cytochrome b of a specimen of Zarudny's rock shrew (Crocidura zarudnyi) from Baluchestan, southeastern Iran, to clarify its cytogenetic and molecular relationships with other Eurasian species of Crocidura. According to the karyotype (2N = 40, FN = 50), Zarudny's rock shrew belongs to the group of the lesser white-toothed shrew (C. suaveolens), which is different from other known crocidurine karyotypes, considering the combination of the diploid and fundamental number of chromosomes. Molecular results revealed that C. zarudnyi is included in a monophyletic clade with the C. suaveolens group, where it is a sister taxon to the others (mean Kimura 2-parameter distance = 9.7%).

Molecular changes underlying mammalian phenotypic innovation

Mammals are characterized by specific phenotypic traits that include lactation, hair, and relatively large brains with unique structures. Individual mammalian lineages have, in turn, evolved characteristic traits that distinguish them from others. These include obvious anatom¬ical differences but also differences related to reproduction, life span, cognitive abilities, be¬havior. and disease susceptibility. However, the molecular basis of the diverse mammalian phenotypes and the selective pressures that shaped their evolution remain largely unknown. In the first part of my thesis, I analyzed the genetic factors associated with the origin of a unique mammalian phenotype lactation and I studied the selective pressures that forged the transition from oviparity to viviparity. Using a comparative genomics approach and evolutionary simulations, I showed that the emergence of lactation, as well as the appear¬ance of the casein gene family, significantly reduced selective pressure on the major egg-yolk proteins (the vitellogenin family). This led to a progressive loss of vitellogenins, which - in oviparous species - act as storage proteins for lipids, amino acids, phosphorous and calcium in the isolated egg. The passage to internal fertilization and placentation in therian mam¬mals rendered vitellogenins completely dispensable, which ended in the loss of the whole gene family in this lineage. As illustrated by the vitellogenin study, changes in gene content are one possible underlying factor for the evolution of mammalian-specific phenotypes. However, more subtle genomic changes, such as mutations in protein-coding sequences, can also greatly affect the phenotypes. In particular, it was proposed that changes at the level of gene reg¬ulation could underlie many (or even most) phenotypic differences between species. In the second part of my thesis, I participated in a major comparative study of mammalian tissue transcriptomes, with the goal of understanding how evolutionary forces affected expression patterns in the past 200 million years of mammalian evolution. I showed that, while com¬parisons of gene expressions are in agreement with the known species phylogeny, the rate of expression evolution varies greatly among lineages. Species with low effective population size, such as monotremes and hominoids, showed significantly accelerated rates of gene expression evolution. The most likely explanation for the high rate of gene expression evolution in these lineages is the accumulation of mildly deleterious mutations in regulatory regions, due to the low efficiency of purifying selection. Thus, our observations are in agreement with the nearly neutral theory of molecular evolution. I also describe substantial differences in evolutionary rates between tissues, with brain being the most constrained (especially in primates) and testis significantly accelerated. The rate of gene expression evolution also varies significantly between chromosomes. In particular, I observed an acceleration of gene expression changes on the X chromosome, probably as a result of adaptive processes associated with the origin of therian sex chromosomes. Lastly, I identified several individual genes as well as co-regulated expression modules that have undergone lineage specific expression changes and likely under¬lie various phenotypic innovations in mammals. The methods developed during my thesis, as well as the comprehensive gene content analyses and transcriptomics datasets made available by our group, will likely prove to be useful for further exploratory analyses of the diverse mammalian phenotypes.