van der Waals coefficients for positronium interactions with atoms

The random-phase approximation with exchange (RPAE) is used with a B-spline basis to compute dynamic dipole polarizabilities of noble-gas atoms and several other closed-shell atoms (Be, Mg, Ca, Zn, Sr, Cd, and Ba). From these, values of the van der Waals C₆ constants for positronium interactions with these atoms are determined and compared with existing data. After correcting the RPAE polarizabilities to fit the most accurate static polarizability data, our best predictions of C₆ for Ps–noble-gas pairs are expected to be accurate to within 1%, and to within a few percent for the alkaline-earth metals. We also used accurate dynamic dipole polarizabilities from the literature to compute the C₆ coefficients for the alkali-metal atoms. Implications of increased C₆ values for Ps scattering from more polarizable atoms are discussed.

Exploring molecular variation in Schistosoma japonicum in China

Schistosomiasis is a neglected tropical disease that affects more than 200 million people worldwide. The main disease-causing agents, Schistosoma japonicum, S. mansoni and S. haematobium, are blood flukes that have complex life cycles involving a snail intermediate host. In Asia, S. japonicum causes hepatointestinal disease (schistosomiasis japonica) and is challenging to control due to a broad distribution of its snail hosts and range of animal reservoir hosts. In China, extensive efforts have been underway to control this parasite, but genetic variability in S. japonicum populations could represent an obstacle to eliminating schistosomiasis japonica. Although a draft genome sequence is available for S. japonicum, there has been no previous study of molecular variation in this parasite on a genome-wide scale. In this study, we conducted the first deep genomic exploration of seven S. japonicum populations from mainland China, constructed phylogenies using mitochondrial and nuclear genomic data sets, and established considerable variation between some of the populations in genes inferred to be linked to key cellular processes and/or pathogen-host interactions. Based on the findings from this study, we propose that verifying intraspecific conservation in vaccine or drug target candidates is an important first step toward developing effective vaccines and chemotherapies against schistosomiasis.

An Immunogenetic Signature of Ongoing Antigen Interactions in Splenic Marginal Zone Lymphoma Expressing IGHV1-2*04 Receptors.

Purpose: Prompted by the extensive biases in the immunoglobulin (IG) gene repertoire of splenic marginal-zone lymphoma (SMZL), supporting antigen selection in SMZL ontogeny, we sought to investigate whether antigen involvement is also relevant post-transformation.

Experimental Design: We conducted a large-scale subcloning study of the IG rearrangements of 40 SMZL cases aimed at assessing intraclonal diversification (ID) due to ongoing somatic hypermutation (SHM).

Results: ID was identified in 17 of 21 (81%) rearrangements using the immunoglobulin heavy variable (IGHV)1-2*04 gene versus 8 of 19 (40%) rearrangements utilizing other IGHV genes (P= 0.001). ID was also evident in most analyzed IG light chain gene rearrangements, albeit was more limited compared with IG heavy chains. Identical sequence changes were shared by subclones from different patients utilizing the IGHV1-2*04 gene, confirming restricted ongoing SHM profiles. Non-IGHV1-2*04 cases displayed both a lower number of ongoing SHMs and a lack of shared mutations (per group of cases utilizing the same IGHV gene).

Conclusions: These findings support ongoing antigen involvement in a sizable portion of SMZL and further argue that IGHV1-2*04 SMZL may represent a distinct molecular subtype of the disease.

CTCF modulates Estrogen Receptor function through specific chromatin and nuclear matrix interactions

Enhancer regions and transcription start sites of estrogen-target regulated genes are connected by means of Estrogen Receptor long-range chromatin interactions. Yet, the complete molecular mechanisms controlling the transcriptional output of engaged enhancers and subsequent activation of coding genes remain elusive. Here, we report that CTCF binding to enhancer RNAs is enriched when breast cancer cells are stimulated with estrogen. CTCF binding to enhancer regions results in modulation of estrogen-induced gene transcription by preventing Estrogen Receptor chromatin binding and by hindering the formation of additional enhancer-promoter ER looping. Furthermore, the depletion of CTCF facilitates the expression of target genes associated with cell division and increases the rate of breast cancer cell proliferation. We have also uncovered a genomic network connecting loci enriched in cell cycle regulator genes to nuclear lamina that mediates the CTCF function. The nuclear lamina and chromatin interactions are regulated by estrogen-ER. We have observed that the chromatin loops formed when cells are treated with estrogen establish contacts with the nuclear lamina. Once there, the portion of CTCF associated with the nuclear lamina interacts with enhancer regions, limiting the formation of ER loops and the induction of genes present in the loop. Collectively, our results reveal an important, unanticipated interplay between CTCF and nuclear lamina to control the transcription of ER target genes, which has great implications in the rate of growth of breast cancer cells.

Interactions of Cytochrome P450 3A4 with Phospholipid Bilayer Nanodiscs

Thesis (Ph.D.)--University of Washington, 2016-08

Renal Drug Transporters – Clinical Significance and Role in Drug-Drug Interactions

Thesis (Ph.D.)--University of Washington, 2016-07

Studies of Protein-Protein and Protein-Water Interactions by Small Angle X-Ray Scattering, Terahertz Spectroscopy, ASMOS, And Computer Simulation

The protein folding problem has been one of the most challenging subjects in biological physics due to its complexity. Energy landscape theory based on statistical mechanics provides a thermodynamic interpretation of the protein folding process. We have been working to answer fundamental questions about protein-protein and protein-water interactions, which are very important for describing the energy landscape surface of proteins correctly. At first, we present a new method for computing protein-protein interaction potentials of solvated proteins directly from SAXS data. An ensemble of proteins was modeled by Metropolis Monte Carlo and Molecular Dynamics simulations, and the global X-ray scattering of the whole model ensemble was computed at each snapshot of the simulation. The interaction potential model was optimized and iterated by a Levenberg-Marquardt algorithm. Secondly, we report that terahertz spectroscopy directly probes hydration dynamics around proteins and determines the size of the dynamical hydration shell. We also present the sequence and pH-dependence of the hydration shell and the effect of the hydrophobicity. On the other hand, kinetic terahertz absorption (KITA) spectroscopy is introduced to study the refolding kinetics of ubiquitin and its mutants. KITA results are compared to small angle X-ray scattering, tryptophan fluorescence, and circular dichroism results. We propose that KITA monitors the rearrangement of hydrogen bonding during secondary structure formation. Finally, we present development of the automated single molecule operating system (ASMOS) for a high throughput single molecule detector, which levitates a single protein molecule in a 10 µm diameter droplet by the laser guidance. I also have performed supporting calculations and simulations with my own program codes.

Human Ezrin: Prediction of the 3D structure and interactions with mTOR

The protein Ezrin, is a member of the ERM family (Ezrin, Radixin and Moesin) that links the F-actin to the plasma membrane. The protein is made of three domains namely the FERM domain, a central α-helical domain and the CERMAD domain. The residues in Ezrin such as Ser66, Tyr145, Tyr353 and Tyr477 regulate the function of the protein through phosphorylation. The protein is found in two distinct conformations of active and dormant (inactive) state. The initial step during the conformation change is the breakage of intramolecular interaction in dormant Ezrin by phosphorylation of residue Thr567. The dormant structure of human Ezrin was predicted computationally since only partial active form structure was available. The validation analysis showed that 99.7% residues were positioned in favored, allowed and generously allowed regions of the Ramachandran plot. The Z-score of Ezrin was −7.36, G-factor was 0.1, and the QMEAN score of the model was 0.61 indicating a good model for human Ezrin. The comparison of the conformations of the activated and dormant Ezrin showed a major shift in the F2 lobe (residues 142-149 and 161-177) while changes in the conformation induced mobility shifts in lobe F3 (residues 261 to 267). The 3D positions of the phosphorylation sites Tyr145, Tyr353, Tyr477, Tyr482 and Thr567 were also located. Using targeted molecular dynamic simulation, the molecular movements during conformational change from active to dormant were visualized. The dormant Ezrin auto-inhibits itself by a head-to-tail interaction of the N-terminal and C-terminal residues. The trajectory shows the breakage of the interactions and mobility of the CERMAD domain away from the FERM domain. Protein docking and clustering analysis were used to predict the residues involved in the interaction between dormant Ezrin and mTOR. Residues Tyr477 and Tyr482 were found to be involved in interaction with mTOR.

Exploring protein-RNA interactions with site-directed mutagenesis and phage display

The importance of RNA as a mediator of genetic information is widely appreciated. RNA molecules also participate in the regulation of various post-transcriptional activities, such as mRNA splicing, editing, RNA stability and transport. Their regulatory roles for these activities are highly dependent on finely tuned associations with cognate proteins. The RNA recognition motif (RRM) is an ancient RNA binding module that participates in hundreds of essential activities where specific RNA recognition is required. We have applied phage display and site-directed mutagenesis to dissect principles of RRM-controlled RNA recognition. The model systems we are investigating are U1A and CUG-BP1. In this dissertation, the molecular basis of the binding affinity of U1A-RNA beyond individual contacts was investigated. We have identified and evaluated the contributions of the local cooperativity formed by three neighboring residues (Asn15, Asn16 and Glu19) to the stability of the U1A-RNA complex. The localized cooperative network was mapped by double-mutant cycles and explored using phage display. We also showed that a cluster of these residues forms a “hot spot” on the surface of U1A; a single substitution at position 19 with Gln or His can alter the binding properties of U1A to recognize a non-cognate G4U RNA. Finally, we applied a deletion analysis of CUG-BP1 to define the contributions of individual RRMs and RRM combinations to the stability of the complex formed between CUG-BP1 and the GRE sequence. The preliminary results showed RRM3 of CUG-BP1 is a key domain for RNA binding. It possibly binds to the GRE sequence cooperatively with RRM2 of CUG-BP1. RRM1 of CUG-BP1 is not required for GRE recognition, but may be important for maintaining the stability of the full-length CUG-BP1.

Novas estratégias para entender o sistema de influxo de antibióticos e a resistência antimicrobiana a nível molecular

As fluoroquinolonas são antibióticos que têm um largo espectro de ação contra bactérias, especialmente Gram-negativas. O seu mecanismo de ação assenta na inibição de enzimas responsáveis pela replicação do DNA. Porém, devido ao seu uso indevido, o surgimento de resistência bacteriana a estes antibióticos tem-se tornado um grave problema de saúde pública. Uma vez que os seus alvos de ação se situam no meio intracelular, a redução da permeabilidade da membrana externa de bactérias Gram-negativas constitui um dos mecanismos de resistência mais conhecidos. Esta redução é associada à baixa expressão ou mutações em porinas necessárias para permitir o seu transporte, mais concretamente, da OmpF. Estudos prévios demonstraram que a coordenação de fluoroquinolonas com iões metálicos divalentes e 1,10-fenantrolina (genericamente designados metaloantibióticos) são potenciais candidatos como alternativa às fluoroquinolonas convencionais. Estes metaloantibióticos exibem um efeito antimicrobiano comparável ou superior à fluoroquinolona na forma livre, mas parecem ter uma via de translocação diferente, independente de porinas. Estas diferenças no mecanismo de captura podem ser fundamentais para contornar a resistência bacteriana. De forma a compreender o papel dos lípidos no mecanismo de entrada dos metaloantibióticos, estudou-se a interação e localização dos metaloantibióticos da Ciprofloxacina (2ª geração), da Levofloxacina (3ª geração) e Moxifloxacina (4ª geração) com um modelo de membranas de Escherichia coli desprovido de porinas. Estes estudos foram realizados através de técnicas de espectroscopia de fluorescência, por medições em modo estacionário e resolvida no tempo. Os coeficientes de partição determinados demonstraram uma interação mais elevada dos metaloantibióticos relativamente às respetivas fluoroquinolonas na forma livre, um facto que está diretamente relacionado com as espécies existentes em solução a pH fisiológico. Os estudos de localização mostraram que estes metaloantibióticos devem estar inseridos na membrana bacteriana, confirmando a sua entrada independente de porinas. Este mecanismo de entrada, pela via hidrofóbica, é potenciado por interações eletrostáticas entre as espécies catiónicas de metaloantibiótico que existem a pH 7,4 e os grupos carregados negativamente dos fosfolípidos da membrana. Desta forma, os resultados obtidos neste estudo sugerem que a via de entrada dos metaloantibióticos e das respetivas fluoroquinolonas deve ser diferente. Os metaloantibióticos são candidatos adequados para a realização de mais testes laboratoriais e uma alternativa promissora para substituir as fluoroquinolonas convencionais, uma vez que parecem ultrapassar um dos principais mecanismos de resistência bacteriana a esta classe de antibióticos.

Whole Genome Analysis of Gene Expression Reveals Coordinated Activation of Signaling and Metabolic Pathways during Pollen-Pistil Interactions in Arabidopsis

Plant reproduction depends on the concerted activation of many genes to ensure correct communication between pollen and pistil. Here, we queried the whole transcriptome of Arabidopsis (Arabidopsis thaliana) in order to identify genes with specific reproductive functions. We used the Affymetrix ATH1 whole genome array to profile wild-type unpollinated pistils and unfertilized ovules. By comparing the expression profile of pistils at 0.5, 3.5, and 8.0 h after pollination and applying a number of statistical and bioinformatics criteria, we found 1,373 genes differentially regulated during pollen-pistil interactions. Robust clustering analysis grouped these genes in 16 time-course clusters representing distinct patterns of regulation. Coregulation within each cluster suggests the presence of distinct genetic pathways, which might be under the control of specific transcriptional regulators. A total of 78% of the regulated genes were expressed initially in unpollinated pistil and/or ovules, 15% were initially detected in the pollen data sets as enriched or preferentially expressed, and 7% were induced upon pollination. Among those, we found a particular enrichment for unknown transcripts predicted to encode secreted proteins or representing signaling and cell wall-related proteins, which may function by remodeling the extracellular matrix or as extracellular signaling molecules. A strict regulatory control in various metabolic pathways suggests that fine-tuning of the biochemical and physiological cellular environment is crucial for reproductive success. Our study provides a unique and detailed temporal and spatial gene expression profile of in vivo pollen-pistil interactions, providing a framework to better understand the basis of the molecular mechanisms operating during the reproductive process in higher plants.

Host cellular regulatory networks in dengue virus-human interactions

Dengue fever is one of the most important mosquito-borne diseases worldwide and is caused by infection with dengue virus (DENV). The disease is endemic in tropical and sub-tropical regions and has increased remarkably in the last few decades. At present, there is no antiviral or approved vaccine against the virus. Treatment of dengue patients is usually supportive, through oral or intravenous rehydration, or by blood transfusion for more severe dengue cases. Infection of DENV in humans and mosquitoes involves a complex interplay between the virus and host factors. This results in regulation of numerous intracellular processes, such as signal transduction and gene transcription which leads to progression of disease. To understand the mechanisms underlying the disease, the study of virus and host factors is therefore essential and could lead to the identification of human proteins modulating an essential step in the virus life cycle. Knowledge of these human proteins could lead to the discovery of potential new drug targets and disease control strategies in the future. Recent advances of high throughput screening technologies have provided researchers with molecular tools to carry out investigations on a large scale. Several studies have focused on determination of the host factors during DENV infection in human and mosquito cells. For instance, a genome-wide RNA interference (RNAi) screen has identified host factors that potentially play an important role in both DENV and West Nile virus replication (Krishnan et al. 2008). In the present study, a high-throughput yeast two-hybrid screen has been utilised in order to identify human factors interacting with DENV non-structural proteins. From the screen, 94 potential human interactors were identified. These include proteins involved in immune signalling regulation, potassium voltage-gated channels, transcriptional regulators, protein transporters and endoplasmic reticulum-associated proteins. Validation of fifteen of these human interactions revealed twelve of them strongly interacted with DENV proteins. Two proteins of particular interest were selected for further investigations of functional biological systems at the molecular level. These proteins, including a nuclear-associated protein BANP and a voltage-gated potassium channel Kv1.3, both have been identified through interaction with the DENV NS2A. BANP is known to be involved in NF-kB immune signalling pathway, whereas, Kv1.3 is known to play an important role in regulating passive flow of potassium ions upon changes in the cell transmembrane potential. This study also initiated a construction of an Aedes aegypti cDNA library for use with DENV proteins in Y2H screen. However, several issues were encountered during the study which made the library unsuitable for protein interaction analysis. In parallel, innate immune signalling was also optimised for downstream analysis. Overall, the work presented in this thesis, in particular the Y2H screen provides a number of human factors potentially targeted by DENV during infection. Nonetheless, more work is required to be done in order to validate these proteins and determine their functional properties, as well as testing them with infectious DENV to establish a biological significance. In the long term, data from this study will be useful for investigating potential human factors for development of antiviral strategies against dengue.

Coarse Graining to Invesitigate Peptide-Lipid Interactions

Experimental characterization of molecular details is challenging, and although single molecule experiments have gained prominence, oligomer characterization remains largely unexplored. The ability to monitor the time evolution of individual molecules while they self assemble is essential in providing mechanistic insights about biological events. Molecular dynamics (MD) simulations can fill the gap in knowledge between single molecule experiments and ensemble studies like NMR, and are increasingly used to gain a better understanding of microscopic properties. Coarse-grained (CG) models aid in both exploring longer length and time scale molecular phenomena, and narrowing down the key interactions responsible for significant system characteristics. Over the past decade, CG techniques have made a significant impact in understanding physicochemical processes. However, the realm of peptide-lipid interfacial interactions, primarily binding, partitioning and folding of amphipathic peptides, remains largely unexplored compared to peptide folding in solution. The main drawback of existing CG models is the inability to capture environmentally sensitive changes in dipolar interactions, which are indigenous to protein folding, and lipid dynamics. We have used the Drude oscillator approach to incorporate structural polarization and dipolar interactions in CG beads to develop a minimalistic peptide model, WEPPROM (Water Explicit Polarizable PROtein Model), and a lipid model WEPMEM (Water Explicit Polarizable MEmbrane Model). The addition of backbone dipolar interactions in a CG model for peptides enabled us to achieve alpha-beta secondary structure content de novo, without any added bias. As a prelude to studying amphipathic peptide-lipid membrane interactions, the balance between hydrophobicity and backbone dipolar interactions in driving ordered peptide aggregation in water and at a hydrophobic-hydrophilic interface, was explored. We found that backbone dipole interactions play a crucial role in driving ordered peptide aggregation, both in water and at hydrophobic-hydrophilic interfaces; while hydrophobicity is more relevant for aggregation in water. A zwitterionic (POPC: 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine) and an anionic lipid (POPS: 1-palmitoyl-2-oleoyl-sn-glycero-3-phospho-L-serine) are used as model lipids for WEPMEM. The addition of head group dipolar interactions in lipids significantly improved structural, dynamic and dielectric properties of the model bilayer. Using WEPMEM and WEPPROM, we studied membrane-induced peptide folding of a cationic antimicrobial peptide with anticancer activity, SVS-1. We found that membrane-induced peptide folding is driven by both (a) cooperativity in peptide self interaction and (b) cooperativity in membrane-peptide interactions. The dipolar interactions between the peptide and the lipid head-groups contribute to stabilizing folded conformations. The role of monovalent ion size and peptide concentration in driving lipid domain formation in anionic/zwitterionic lipid mixtures was also investigated. Our study suggest monovalent ion size to be a crucial determinant of interaction with lipid head groups, and hence domain formation in lipid mixtures. This study reinforces the role of dipole interactions in protein folding, lipid membrane properties, membrane induced peptide folding and lipid domain formation. Therefore, the models developed in this thesis can be used to explore a multitude of biomolecular processes, both at longer time-scales and larger system sizes.

Caractérisation des interactions établies par la région riche en prolines de la ligase de l’ubiquitine Itch

Itch est un membre de la famille des ligases de l’ubiquitine de type CWH (C2-WW- HECT) impliqué dans le contrôle de la signalisation inflammatoire, des facteurs de transcription et le tri des récepteurs membranaires. La fonction d’Itch implique généralement sa capacité à induire la dégradation de ses substrats. Pour accomplir cette fonction, Itch doit d’abord interagir avec ses cibles. Itch possède quatre domaines WW lui permettant d’accomplir la majorité de ses fonctions. En plus de ces domaines, Itch possède une PRR (région riche en prolines) unique parmi les ligases CWH. Cette région est bien conservée chez les vertébrés, ce qui suggère son importance. Cette région permet à Itch d’interagir avec des protéines contenant un domaine SH3 (Src homology 3). Plusieurs partenaires SH3 furent identifiés, cependant l’on connait peu de choses concernant la fonction et l’établissement de ces complexes. Dans ce projet, nous avons analysé les propriétés de liaison d’un sous-groupe de protéines à domaine SH3 impliquées dans l’endocytose et la signalisation cellulaire. Nos travaux ont permis d’identifier de nouveaux partenaires et aussi de déterminer que différents domaines SH3 ciblent la même région riche en prolines, mais impliquent des résidus distincts. Ces résultats démontrent la variété des propriétés de liaison démontrées par la PRR d’Itch et sa préférence marquée pour l’Endophiline. Parmi les partenaires identifiés, Grb2 (Growth factor receptor-bound protein 2) est particulièrement intéressant en raison de son rôle crucial dans la signalisation cellulaire. Nous avons démontré ici qu’Itch ubiquityle Grb2, mais ne cause pas sa dégradation, contrairement à l’Endophiline. Nos travaux démontrent que la PRR d’Itch est versatile quant à ses interactions et leurs conséquences.

Insights into the Specificity of Thioredoxin Reductase-Thioredoxin Interactions. A Structural and Functional Investigation of the Yeast Thioredoxin System

The enzymatic activity of thioredoxin reductase enzymes is endowed by at least two redox centers: a flavin and a dithiol/disulfide CXXC motif. The interaction between thioredoxin reductase and thioredoxin is generally species-specific, but the molecular aspects related to this phenomenon remain elusive. Here, we investigated the yeast cytosolic thioredoxin system, which is composed of NADPH, thioredoxin reductase (ScTrxR1), and thioredoxin 1 (ScTrx1) or thioredoxin 2 (ScTrx2). We showed that ScTrxR1 was able to efficiently reduce yeast thioredoxins (mitochondrial and cytosolic) but failed to reduce the human and Escherichia coli thioredoxin counterparts. To gain insights into this specificity, the crystallographic structure of oxidized ScTrxR1 was solved at 2.4 angstrom resolution. The protein topology of the redox centers indicated the necessity of a large structural rearrangement for FAD and thioredoxin reduction using NADPH. Therefore, we modeled a large structural rotation between the two ScTrxR1 domains (based on the previously described crystal structure, PDB code 1F6M). Employing diverse approaches including enzymatic assays, site-directed mutagenesis, amino acid sequence alignment, and structure comparisons, insights were obtained about the features involved in the species-specificity phenomenon, such as complementary electronic parameters between the surfaces of ScTrxR1 and yeast thioredoxin enzymes and loops and residues (such as Ser(72) in ScTrx2). Finally, structural comparisons and amino acid alignments led us to propose a new classification that includes a larger number of enzymes with thioredoxin reductase activity, neglected in the low/high molecular weight classification.