An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structuralvariations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

An automated framework for understanding structural variations in the binding grooves of MHC class II molecules

Background: MHC/HLA class II molecules are important components of the immune system and play a critical role in processes such as phagocytosis. Understanding peptide recognition properties of the hundreds of MHC class II alleles is essential to appreciate determinants of antigenicity and ultimately to predict epitopes. While there are several methods for epitope prediction, each differing in their success rates, there are no reports so far in the literature to systematically characterize the binding sites at the structural level and infer recognition profiles from them. Results: Here we report a new approach to compare the binding sites of MHC class II molecules using their three dimensional structures. We use a specifically tuned version of our recent algorithm, PocketMatch. We show that our methodology is useful for classification of MHC class II molecules based on similarities or differences among their binding sites. A new module has been used to define binding sites in MHC molecules. Comparison of binding sites of 103 MHC molecules, both at the whole groove and individual sub-pocket levels has been carried out, and their clustering patterns analyzed. While clusters largely agree with serotypic classification, deviations from it and several new insights are obtained from our study. We also present how differences in sub-pockets of molecules associated with a pair of autoimmune diseases, narcolepsy and rheumatoid arthritis, were captured by PocketMatch(13). Conclusion: The systematic framework for understanding structural variations in MHC class II molecules enables large scale comparison of binding grooves and sub-pockets, which is likely to have direct implications towards predicting epitopes and understanding peptide binding preferences.

Prognostic significance of cyclooxygenase-2 and associated molecules in gastric cancer

Background and aims: Low stage and curative surgery are established factors for improved survival in gastric cancer. However, not all low-stage patients have a good prognosis. Cyclooxygenase-2 (COX-2) is known to associate with reduced survival in several cancers, and has been shown to play an important role in gastric carcinogenesis. Since new and better prognostic markers are needed for gastric cancer, we studied the prognostic significance of COX-2 and of markers that associate with COX-2 expression. We also studied markers reflecting proliferation and apoptosis, and evaluated their association with COX-2. Our purpose was to construct an accurate prognostic model by combining tissue markers and clinicopathogical factors. Materials and methods: Of 342 consecutive patients who underwent surgery for gastric cancer at Meilahti Hospital, Helsinki University Central Hospital, 337 were included in this study. Low stages I to II were represented by 141 (42%) patients, and high stages III to IV by 196 (58%). Curative surgery was performed on 176 (52%) patients. Survival data were obtained from the national registers. Slides from archive tissue blocks were prepared for immunohistochemistry by use of COX-2, human antigen R (HuR), cyclin A, matrix metalloproteinases 2 and 9 (MMP-2, MMP-9), and Ki-67 antibodies. Immunostainings were scored by microscopy, and scores were entered into a database. Associations of tumor markers with clinicopathological factors were calculated, as well as associations with p53, p21, and results of flow cytometry from earlier studies. Survival analysis was performed by the Kaplan-Meier method, and Cox multivariate models were reconstructed. Cell culture experiments were performed to explore the effect of small interfering (si)RNA of HuR on COX-2 expression in a TMK-1 gastric cancer cell line. Results: Overall 5-year survival was 35.1%. Study I showed that COX-2 was an independent prognostic factor, and that the prognostic impact of COX-2 was more pronounced in low-stage patients. Cytoplasmic HuR expression also associated with reduced survival in gastric cancer patients in a non-independent manner. Cell culture experiments showed that HuR can regulate COX-2 expression in TMK-1 cells in vitro, with an association also between COX-2 and HuR tissue expression in a clinical material. In Study II, cyclin A was an independent prognostic factor and was associated with HuR expression in the gastric cancer material. The results of Study III showed that epithelial MMP-2 associated with survival in univariate, but not in multivariate analysis. However, MMP-9 showed no prognostic value. MMP-2 expression was associated with COX-2 expression. In Study IV, the prognostic power of COX-2 was compared with that of all tested markers associated with survival in Studies I to III, as well as with p21, p53, and flow cytometry results. COX-2 and p53 were independent prognostic factors, and COX-2 expression was associated with that of p53 and Ki-67 and also with aneuploidy. Conclusions: COX-2 is an independent prognostic factor in gastric cancer, and its prognostic power emerges especially in low stage cancer. COX-2 is regulated by HuR, and is associated with factors reflecting invasion, proliferation, and apoptosis. In an extended multivariate model, COX-2 retained its position as an independent prognosticator. COX-2 can be considered a promising new prognostic marker in gastric cancer.

Excited-state structures of molecules with nearby electronic states

Vibronic coupling among the nearby excited electronic states via the in-plane and the out-of-plane nuclear motions is examined in benzene, pyrazine, formaldehyde and thioformaldehyde. Results reveal that in benzene the structure distorts via the most active nuclear bending (planar) motion while in the other molecules the structures distort through an out-of-plane bending motion in their respective lowest excited states.

Effects of electron correlations in conjugated organic molecules and solids

The discrepancies between the non-interacting models and experimental results for conjugated systems is highlighted in this brief review. The interacting model hamiltonians correctly give the forbidden singlet state below the optical gap in polyenes and also explain both the nonvanishing optical gap in polyacetylenes and the vanishing optical gap in symmetric cyanine dyes. The negative spin densities in polyene radicals is also understood in terms of a correlated picture. The role of electron-electron interactions in other strongly correlated systems, such as polydiacetylene and mixed and segregated stack charge transfer solids, are also briefly discussed.

A general procedure for generation of curved DNA molecules.

A general method for generation of base-pairs in a curved DNA structure, for any prescribed values of helical parameters--unit rise (h), unit twist (theta), wedge roll (theta R) and wedge tilt (theta T), propeller twist (theta p) and displacement (D) is described. Its application for generation of uniform as well curved structures is also illustrated with some representative examples. An interesting relationship is observed between helical twist (theta), base-pair parameters theta x, theta y and the wedge parameters theta R, theta T, which has important consequences for the description and estimation of DNA curvature.

Glioblastoma-Specific Protein Interaction Network Identifies PP1A and CSK21 as Connecting Molecules between Cell Cycle-Associated Genes

Glioblastoma (GBM; grade IV astrocytoma) is a very aggressive form of brain cancer with a poor survival and few qualified predictive markers. This study integrates experimentally validated genes that showed specific upregulation in GBM along with their protein-protein interaction information. A system level analysis was used to construct GBM-specific network. Computation of topological parameters of networks showed scale-free pattern and hierarchical organization. From the large network involving 1,447 proteins, we synthesized subnetworks and annotated them with highly enriched biological processes. A careful dissection of the functional modules, important nodes, and their connections identified two novel intermediary molecules CSK21 and protein phosphatase 1 alpha (PP1A) connecting the two subnetworks CDC2-PTEN-TOP2A-CAV1-P53 and CDC2-CAV1-RB-P53-PTEN, respectively. Real-time quantitative reverse transcription-PCR analysis revealed CSK21 to be moderately upregulated and PP1A to be overexpressed by 20-fold in GBM tumor samples. Immunohistochemical staining revealed nuclear expression of PP1A only in GBM samples. Thus, CSK21 and PP1A, whose functions are intimately associated with cell cycle regulation, might play key role in gliomagenesis. Cancer Res; 70(16); 6437-47. (C)2010 AACR.

Identification of molecules relevant for the invasiveness of fibrosarcomas and melanomas

Cancer is becoming the leading cause of deaths in the world. As 90% of all deaths from cancer are caused by metastasis, discovery of the mechanisms behind cancer cell invasion and metastasis is of utmost importance. Only new effective therapies targeting cancer progression can reduce cancer mortality rates. The aim of this study was to identify molecules that are relevant for tumor cell invasion and spreading in fibrosarcomas and melanomas, and to analyze their potential for cancer biomarkers or therapeutic targets. First, the gene expression changes of normal cells and transformed cells showing high invasiveness, S-adenosylmethionine decarboxylase (AdoMetDC)-transfected murine fibroblasts and human melanoma cells, were studied by microarray analyses. The function of the identified candidate molecules were then studied in detail in these cell lines. Finally, the physiological relevance of the identified changes was studied by immunohistochemical analyses of human sarcoma and melanoma specimens or by a mouse xenograft model. In fibrosarcoma cells, the most remarkable change detected was a dramatic up-regulation of the actin-sequestering molecule thymosin beta 4 (TB4), which was shown to be important for the transformed phenotype of the AdoMetDC-transfected cells (Amdc-s and -as). A sponge toxin latrunculin A, inhibiting the binding of TB4 to actin, was found to selectively inhibit the migration and invasion of these cells. Further, Amdc-s-induced mouse tumors and human high-grade sarcomas were found to show intense TB4 immunostaining. In addition to TB4, integrin subunits alfa 6 and beta 7 (ItgA6 and ItgB7) were found to be up-regulated in Amdc-s and -as cells. ItgA6 was shown to dimerize mainly with ItgB1 in Amdc-s. Inhibition of ItgA6 or ItgB1 function with neutralizing antibodies fully blocked the invasiveness of Amdc-s cells, and importantly also human HT-1080 fibrosarcoma cells, in three-dimensional (3D)-Matrigel mimicking tumor extracellular matrix (ECM). By immunohistochemical analyses, strong staining for ITGA6 was detected in human high-grade fibrosarcomas and other sarcomas, especially at the invasion fronts of the tumors. In the studied melanoma cell lines, the expression levels of the adhesion-related ECM proteins tenascin-C (TN-C), fibronectin (FN), and transforming growth factor beta-induced (TGFBI) were found to be highly up-regulated. By immunohistochemistry, intense TN-C and FN staining was detected in invasive and metastatic melanoma tumors, showing co-localization (together with procollagen-I) in tubular meshworks and channels around the invading melanoma cells. In vitro, TN-C and FN were further found to directly stimulate the migration of melanoma cells in 3D-collagen-I matrix. The third candidate protein, TGFBI, was found to be an anti-adhesive molecule for melanoma cells, and knockdown of its expression in metastatic melanoma cells (TGFBI-KD cells) led to dramatically impaired tumor growth in immunocompromized mice. Interestingly, the control tumors showed intense TGFBI immunostaining in the invasion fronts, showing partial co-localization with the fibrillar FN staining, whereas the small TGFBI-KD cell-induced tumors displayed amorphous, non-fibrillar FN staining. These data suggest an important role for TGFBI in FN fibrillogenesis and melanoma progression. In conclusion, we have identified several invasion-related molecules, which show potential for cancer diagnostic or prognostic markers, or therapeutic targets. Based on our previous and present fibrosarcoma studies, we propose the possibility of using ITGA6 antagonists (affecting tumor cell adhesion) in combination with TB4 inhibitors (affecting tumor cell migration) and cathepsin L inhibitors (affecting the degradation of basement membrane and ECM proteins) for the treatment of fibrosarcomas and other tumors overexpressing these molecules. With melanoma cells, in turn, we point to the importance of three secreted ECM proteins, TN-C, FN, and TGFBI, in melanoma progression. Of these, especially the potential of TN-C as a prognostic melanoma biomarker and TGFBI as a promising therapeutic target molecule are clearly worth additional studies.

Stabilization of the bio-membrane by small molecules: interaction of trehalose with the phospholipid bilayer

Anhydrobiotic organisms undergo periods of acute dehydration during their life cycle. It is of interest to understand how the biomembrane remains intact through such stress. A disaccharide, trehalose, which is metabolised during anhydrobiosis is found to prevent disruption of model membrane systems. Molecular modelling techniques are used to investigate the possible mode of interaction of trehalose with a model monolayer. The objective is to maximise hydrogen bonding between the two systems. A phospholipid matrix consisting of 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) is chosen to represent the monolayer. The crystal structure of DMPC reveals that there are two distinct conformers designated as A and B. An expansion of the monolayer, coplanar with its surface, results in the trehalose molecule being accommodated in a pocket formed by four B conformers. One glucose ring of the sugar rests on the hydrophobic patch provided by the choline methyls of an A conformer. Five hydrogen bonds are formed involving the phosphate oxygens of three of the surrounding B conformers. The model will be discussed with reference to relevant experimental data on the interaction.

Interaction of water molecules in non-identical protein structures

The interaction of the protein atoms with the surrounding water oxygen atoms has been computed for 392 protein chains from 369 protein structures belonging to 90% non-homologous high resolution (<= 1.5 angstrom) protein Structures with a crystallographic R-factor <= 20%. The percentage composition of the polar atoms is found to be 36.3%. An average of 82.55% of water oxygen atoms are found to be in the primary hydration shell and 15.12% in the secondary hydration shell. The average Percentage of interactions of water oxygen atoms with the polar atoms of the main chain and side chain are 54% and 46%. respectively. The interaction of the acidic residues, aspartate and glutamate, with the water oxygen atoms is more when compared to that of the other residues.

Magnetically induced currents and magnetic response properties of molecules

The magnetically induced currents in organic monoring and multiring molecules, in Möbius shaped molecules and in inorganic all-metal molecules have been investigated by means of the Gauge-including magnetically induced currents (GIMIC) method. With the GIMIC method, the ring-current strengths and the ring-current density distributions can be calculated. For open-shell molecules, also the spin current can be obtained. The ring-current pathways and ring-current strengths can be used to understand the magnetic resonance properties of the molecules, to indirectly identify the effect of non-bonded interactions on NMR chemical shifts, to design new molecules with tailored properties and to discuss molecular aromaticity. In the thesis, the magnetic criterion for aromaticity has been adopted. According to this, a molecule which has a net diatropic ring current might be aromatic. Similarly, a molecule which has a net paratropic current might be antiaromatic. If the net current is zero, the molecule is nonaromatic. The electronic structure of the investigated molecules has been resolved by quantum chemical methods. The magnetically induced currents have been calculated with the GIMIC method at the density-functional theory (DFT) level, as well as at the self-consistent field Hartree-Fock (SCF-HF), at the Møller-Plesset perturbation theory of the second order (MP2) and at the coupled-cluster singles and doubles (CCSD) levels of theory. For closed-shell molecules, accurate ring-current strengths can be obtained with a reasonable computational cost at the DFT level and with rather small basis sets. For open-shell molecules, it is shown that correlated methods such as MP2 and CCSD might be needed to obtain reliable charge and spin currents. The basis set convergence has to be checked for open-shell molecules by performing calculations with large enough basis sets. The results discussed in the thesis have been published in eight papers. In addition, some previously unpublished results on the ring currents in the endohedral fullerene Sc3C2@C80 and in coronene are presented. It is shown that dynamical effects should be taken into account when modelling magnetic resonance parameters of endohedral metallofullerenes such as Sc3C2@C80. The ring-current strengths in a series of nano-sized hydrocarbon rings are related to static polarizabilities and to H-1 nuclear magnetic resonance (NMR) shieldings. In a case study on the possible aromaticity of a Möbius-shaped [16]annulene we found that, according to the magnetic criterion, the molecule is nonaromatic. The applicability of the GIMIC method to assign the aromatic character of molecules was confirmed in a study on the ring currents in simple monocylic aromatic, homoaromatic, antiaromatic, and nonaromatic hydrocarbons. Case studies on nanorings, hexaphyrins and [n]cycloparaphenylenes show that explicit calculations are needed to unravel the ring-current delocalization pathways in complex multiring molecules. The open-shell implementation of GIMIC was applied in studies on the charge currents and the spin currents in single-ring and bi-ring molecules with open shells. The aromaticity predictions that are made based on the GIMIC results are compared to other aromaticity criteria such as H-1 NMR shieldings and shifts, electric polarizabilities, bond-length alternation, as well as to predictions provided by the traditional Hückel (4n+2) rule and its more recent extensions that account for Möbius twisted molecules and for molecules with open shells.

Hyperfine effects in the nuclear magnetic resonance of paramagnetic molecules

Paramagnetic, or open-shell, systems are often encountered in the context of metalloproteins, and they are also an essential part of molecular magnets. Nuclear magnetic resonance (NMR) spectroscopy is a powerful tool for chemical structure elucidation, but for paramagnetic molecules it is substantially more complicated than in the diamagnetic case. Before the present work, the theory of NMR of paramagnetic molecules was limited to spin-1/2 systems and it did not include relativistic corrections to the hyperfine effects. It also was not systematically expandable. --- The theory was first expanded by including hyperfine contributions up to the fourth power in the fine structure constant α. It was then reformulated and its scope widened to allow any spin state in any spatial symmetry. This involved including zero-field splitting effects. In both stages the theory was implemented into a separate analysis program. The different levels of theory were tested by demonstrative density functional calculations on molecules selected to showcase the relative strength of new NMR shielding terms. The theory was also tested in a joint experimental and computational effort to confirm assignment of 11 B signals. The new terms were found to be significant and comparable with the terms in the earlier levels of theory. The leading-order magnetic-field dependence of shielding in paramagnetic systems was formulated. The theory is now systematically expandable, allowing for higher-order field dependence and relativistic contributions. The prevailing experimental view of pseudocontact shift was found to be significantly incomplete, as it only includes specific geometric dependence, which is not present in most of the new terms introduced here. The computational uncertainty in density functional calculations of the Fermi contact hyperfine constant and zero-field splitting tensor sets a limit for quantitative prediction of paramagnetic shielding for now.