HLA-DP2 binding prediction by molecular dynamics simulations

Major histocompatibility complex (MHC) II proteins bind peptide fragments derived from pathogen antigens and present them at the cell surface for recognition by T cells. MHC proteins are divided into Class I and Class II. Human MHC Class II alleles are grouped into three loci: HLA-DP, HLA-DQ, and HLA-DR. They are involved in many autoimmune diseases. In contrast to HLA-DR and HLA-DQ proteins, the X-ray structure of the HLA-DP2 protein has been solved quite recently. In this study, we have used structure-based molecular dynamics simulation to derive a tool for rapid and accurate virtual screening for the prediction of HLA-DP2-peptide binding. A combinatorial library of 247 peptides was built using the "single amino acid substitution" approach and docked into the HLA-DP2 binding site. The complexes were simulated for 1 ns and the short range interaction energies (Lennard-Jones and Coulumb) were used as binding scores after normalization. The normalized values were collected into quantitative matrices (QMs) and their predictive abilities were validated on a large external test set. The validation shows that the best performing QM consisted of Lennard-Jones energies normalized over all positions for anchor residues only plus cross terms between anchor-residues.

Peptide binding prediction for the human class II MHC allele HLA-DP2:a molecular docking approach

MHC class II proteins bind oligopeptide fragments derived from proteolysis of pathogen antigens, presenting them at the cell surface for recognition by CD4+ T cells. Human MHC class II alleles are grouped into three loci: HLA-DP, HLA-DQ and HLA-DR. In contrast to HLA-DR and HLA-DQ, HLA-DP proteins have not been studied extensively, as they have been viewed as less important in immune responses than DRs and DQs. However, it is now known that HLA-DP alleles are associated with many autoimmune diseases. Quite recently, the X-ray structure of the HLA-DP2 molecule (DPA*0103, DPB1*0201) in complex with a self-peptide derived from the HLA-DR a-chain has been determined. In the present study, we applied a validated molecular docking protocol to a library of 247 modelled peptide-DP2 complexes, seeking to assess the contribution made by each of the 20 naturally occurred amino acids at each of the nine binding core peptide positions and the four flanking residues (two on both sides).

Molecular basis for reduced estrone sulfate transport and altered modulator sensitivity of transmembrane helix (TM) 6 and TM17 mutants of multidrug resistance protein 1 (ABCC1)

Multidrug resistance protein 1 (MRP1) confers drug resistance and also mediates cellular efflux of many organic anions. MRP1 also transports glutathione (GSH); furthermore, this tripeptide stimulates transport of several substrates, including estrone 3-sulfate. We have previously shown that mutations of Lys(332) in transmembrane helix (TM) 6 and Trp(1246) in TM17 cause different substrate-selective losses in MRP1 transport activity. Here we have extended our characterization of mutants K332L and W1246C to further define the different roles these two residues play in determining the substrate and inhibitor specificity of MRP1. Thus, we have shown that TM17-Trp(1246) is crucial for conferring drug resistance and for binding and transport of methotrexate, estradiol glucuronide, and estrone 3-sulfate, as well as for binding of the tricyclic isoxazole inhibitor N-[3-(9-chloro-3-methyl-4-oxo-4H-isoxazolo-[4,3-c]quinolin-5-yl)-cyclohexylmethyl]-benzamide (LY465803). In contrast, TM6-Lys(332) is important for enabling GSH and GSH-containing compounds to serve as substrates (e.g., leukotriene C(4)) or modulators (e.g., S-decyl-GSH, GSH disulfide) of MRP1 and, further, for enabling GSH (or S-methyl-GSH) to enhance the transport of estrone 3-sulfate and increase the inhibitory potency of LY465803. On the other hand, both mutants are as sensitive as wild-type MRP1 to the non-GSH-containing inhibitors (E)-3-[[[3-[2-(7-chloro-2-quinolinyl)ethenyl]phenyl][[3-(dimethylamino)-3-oxopropyl]thio]methyl]thio]-propanoic acid (MK571), 1-[2-hydroxy-3-propyl-4-[4-(1H-tetrazol-5-yl)butoxy]phenyl]-ethanone (LY171883), and highly potent 6-[4'-carboxyphenylthio]-5[S]-hydroxy-7[E], 11[Z]14[Z]-eicosatetrenoic acid (BAY u9773). Finally, the differing abilities of the cysteinyl leukotriene derivatives leukotriene C(4), D(4), and F(4) to inhibit estradiol glucuronide transport by wild-type and K332L mutant MRP1 provide further evidence that TM6-Lys(332) is involved in the recognition of the gamma-Glu portion of substrates and modulators containing GSH or GSH-like moieties.

Modeling the peptide-T cell receptor interaction by the comparative molecular similarity indices analysis-soft independent modeling of class analogy technique

A set of 38 epitopes and 183 non-epitopes, which bind to alleles of the HLA-A3 supertype, was subjected to a combination of comparative molecular similarity indices analysis (CoMSIA) and soft independent modeling of class analogy (SIMCA). During the process of T cell recognition, T cell receptors (TCR) interact with the central section of the bound nonamer peptide; thus only positions 4−8 were considered in the study. The derived model distinguished 82% of the epitopes and 73% of the non-epitopes after cross-validation in five groups. The overall preference from the model is for polar amino acids with high electron density and the ability to form hydrogen bonds. These so-called “aggressive” amino acids are flanked by small-sized residues, which enable such residues to protrude from the binding cleft and take an active role in TCR-mediated T cell recognition. Combinations of “aggressive” and “passive” amino acids in the middle part of epitopes constitute a putative TCR binding motif

Structural basis for receptor activity-modifying protein-dependent selective peptide recognition by a G protein-coupled receptor

Association of receptor activity-modifying proteins (RAMP1-3) with the G protein-coupled receptor (GPCR) calcitonin receptor-like receptor (CLR) enables selective recognition of the peptides calcitonin gene-related peptide (CGRP) and adrenomedullin (AM) that have diverse functions in the cardiovascular and lymphatic systems. How peptides selectively bind GPCR:RAMP complexes is unknown. We report crystal structures of CGRP analog-bound CLR:RAMP1 and AM-bound CLR:RAMP2 extracellular domain heterodimers at 2.5 and 1.8 Å resolutions, respectively. The peptides similarly occupy a shared binding site on CLR with conformations characterized by a β-turn structure near their C termini rather than the α-helical structure common to peptides that bind related GPCRs. The RAMPs augment the binding site with distinct contacts to the variable C-terminal peptide residues and elicit subtly different CLR conformations. The structures and accompanying pharmacology data reveal how a class of accessory membrane proteins modulate ligand binding of a GPCR and may inform drug development targeting CLR:RAMP complexes.

Cellular and molecular consequences of S100A4-induced motility in rat breast tumour Rama 37 cells

Since the first discovery of S100 members in 1965, their expressions have been affiliated with numerous biological functions in all cells of the body. However, in the recent years, S100A4, a member of this superfamily has emerged as the central target in generating new avenue for cancer therapy as its overexpression has been correlated with cancer patients’ mortality as well as established roles as motility and metastasis promoter. As it has no catalytic activity, S100A4 has to interact with its target proteins to regulate such effects. Up to date, more than 10 S100A4 target proteins have been identified but the mechanical process regulated by S100A4 to induce motility remains vague. In this work, we demonstrated that S100A4 overexpression resulted in actin filaments disorganisation, reduction in focal adhesions, instability of filopodia as well as exhibiting polarised morphology. However, such effects were not observed in truncated versions of S100A4 possibly highlighting the importance of C terminus of S100A4 target recognition. In order to assess some of the intracellular mechanisms that may be involved in promoting migrations, different strategies were used, including active pharmaceutical agents, inhibitors and knockdown experiments. Treatment of S100A4 overexpressing cells with blebbistatin and Y-27632, non muscle myosin IIA (NMMIIA) inhibitors, as well as knockdown of NMMIIA, resulted in motility enhancement and focal adhesions reduction proposing that NMMIIA assisted S100A4 in regulating cell motility but its presence is not essential. Further work done using Cos 7 cell lines, naturally lacking NMMIIA, further demonstrated that S100A4 is capable of regulating cell motility independent of NMMIIA, possibly through poor maturation of focal adhesion. Given that all these experiments highlighted the independency of NMMIIA towards migration, a protein that has been put at the forefront of S100A4-induced motility, we aimed to gather further understanding regarding the other molecular mechanisms that may be at play for motility. Using high throughput imaging (HCI), 3 compounds were identified to be capable of inhibiting S100A4-mediated migration. Although we have yet to investigate the underlying mechanism for their effects, these compounds have been shown to target membrane proteins and the externalisation of S100 proteins, for at least one of the compounds, leading us to speculate that preventing externalisation of S100A4 could potentially regulate cell motility.

MHC class II binding prediction by molecular docking

Proteins of the Major Histocompatibility Complex (MHC) bind self and nonself peptide antigens or epitopes within the cell and present them at the cell surface for recognition by T cells. All T-cell epitopes are MHC binders but not all MCH binders are T-cell epitopes. The MHC class II proteins are extremely polymorphic. Polymorphic residues cluster in the peptide-binding region and largely determine the MHC's peptide selectivity. The peptide binding site on MHC class II proteins consist of five binding pockets. Using molecular docking, we have modelled the interactions between peptide and MHC class II proteins from locus DRB1. A combinatorial peptide library was generated by mutation of residues at peptide positions which correspond to binding pockets (so called anchor positions). The binding affinities were assessed using different scoring functions. The normalized scoring functions for each amino acid at each anchor position were used to construct quantitative matrices (QM) for MHC class II binding prediction. Models were validated by external test sets comprising 4540 known binders. Eighty percent of the known binders are identified in the best predicted 15% of all overlapping peptides, originating from one protein. © 2011 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim.

Functionalized carbon nanotube nanoelectrodes for biomolecular recognition

The objective of this research is to develop nanoscale ultrasensitive transducers for detection of biological species at molecular level using carbon nanotubes as nanoelectrodes. Rapid detection of ultra low concentration or even single DNA molecules are essential for medical diagnosis and treatment, pharmaceutical applications, gene sequencing as well as forensic analysis. Here the use of functionalized single walled carbon nanotubes (SWNT) as nanoscale detection platform for rapid detection of single DNA molecules is demonstrated. The detection principle is based on obtaining electrical signal from a single amine terminated DNA molecule which is covalently bridged between two ends of an SWNT separated by a nanoscale gap. The synthesis, fabrication, chemical functionalization of nanoelectrodes and DNA attachment were optimized to perform reliable electrical characterization these molecules. Using this detection system fundamental study on charge transport in DNA molecule of both genomic and non genomic sequences is performed. We measured an electrical signal of about 30 pA through a hybridized DNA molecule of 80 base pair in length which encodes a portion of sequence of H5N1 gene of avian Influenza A virus. Due the dynamic nature of the DNA molecules the local environment such as ion concentration, pH and temperature significantly influence its physical properties. We observed a decrease in DNA conductance of about 33% in high vacuum conditions. The counterion variation was analyzed by changing the buffer from sodium acetate to tris(hydroxymethyl) aminomethane, which resulted in a two orders of magnitude increase in the conductivity of the DNA. The fabrication of large array of identical SWNT nanoelectrodes was achieved by using ultralong SWNTs. Using these nanoelectrode array we have investigated the sequence dependent charge transport in DNA. A systematic study performed on PolyG - PolyC sequence with varying number of intervening PolyA - PolyT pairs showed a decrease in electrical signal from 180 pA (PolyG - PolyC) to 30 pA with increasing number of the PolyA - PolyT pairs. This work also led to the development of ultrasensitive nanoelectrodes based on enzyme functionalized vertically aligned high density multiwalled CNTs for electrochemical detection of cholesterol. The nanoelectrodes exhibited selectively detection of cholesterol in the presence of common interferents found in human blood.

Estudio en aves silvestres de canto ("Passeriformes") de las moléculas MHC y reclasificación molecular de la familia "Carduelini"

The Major Histocompatibility Complex (MHC) comprises the most polymorphic loci in animals. MHC plays an important role during the first steps of the immune response in vertebrates. In humans, MHC molecules (also named human leukocyte antigens, HLA) were initially regarded as class I or class II molecules. Each of them, presents to different T cells subsets. MHC class I molecules, are heterodimers in which the heavy chain (alpha) has three extracellular domains, two of which (alpha 1 and alpha 2) are polymorphic and conform the antigen recognition sites (ARS). The ARS is thought to be subjected to balancing selection for variability, which is the cause of the very high polymorphism of the MHC molecules. Different pathogenic epitopes would be the evolutionary force causing balancing selection. MHC class I genes have been completely sequenced (α1 and α2 protein domains) and thoroughly studied in Gallus gallus (chicken) as well as in mammals. In fact, the MHC locus was first defined in chicken, specifically in the highly consanguineous variety „Leghorn‟. It has been found that, in the case of chickens the MHC genetic region is considerably smaller than it is in mammals (remarkably shorter introns were found in chickens), and is organized quite differently. The noteworthy presence of short introns in chickens; supported the hypothesis that chicken‟s MHC represented a „minimal essential MHC‟. Until now, it has been assumed that chicken (order Galliformes) MHC was similar to all species included in the whole class Aves...

Structural Basis of RNA Recognition by Mycobacterium tuberculosis MazF Homologues

The MazEF toxin-antitoxin (TA) system consists of the antitoxin MazE and the toxin MazF. MazF is a sequence-specific endoribonuclease that upon activation causes cellular growth arrest and increass the level of persisters. Moreover, MazF-induced cells are in a quasi-dormant state that cells remain metabolically active while stop dividing. The quasi-dormancy is similar to the nonreplicating state of M. tuberculosis during latent tuberculosis, thus suggesting the role of mazEF in M. tuberculosis dormancy and persistence. M. tuberculosis has nine mazEF TA modules, each with different RNA cleavage specificities and implicated in selective gene expression during stress conditions. To date only the Bacillus subtilis MazF-RNA complex structure has been determined. As M. tuberculosis MazF homologues recognize distinct RNA sequences, their molecular mechanisms of substrate specificity remain unclear. By taking advantage of X-ray crystallography, we have determined structures of two M. tuberculosis MazF-RNA complexes, MazF-mt1 (Rv2801c) and MazF-mt3 (Rv1991c) in complex with an uncleavable RNA substrate. These structures have provided the molecular basis of sequence-specific RNA recognition and cleavage by MazF toxins.

Both MazF-mt1-RNA and MazF-mt3-RNA complexes showed similar structural organization with one molecule of RNA bound to a MazF-mt1 or MazF-mt3 dimer and occupying the same pocket within the MazF dimer interface. Similar to B. subtilis MazF-RNA complex, MazF-mt1 and MazF-mt3 displayed a conserved active site architecture, where two highly conserved residues, Arg and Thr, form hydrogen bonds with the scissile phosphate group in the cleavage site of the bound RNA. The MazF-mt1-RNA complex also showed specific interactions with its three-base RNA recognition element. Compared with the B. subtilis MazF-RNA complex, our structures showed that residues involved in sequence-specific recognition of target RNA vary between the MazF homologues, therefore explaining the molecular basis for their different RNA recognition sequences. In addition, local conformational changes of the loops in the RNA binding site of MazF-mt1 appear to play a role in MazF targeting different RNA lengths and sequences. In contrast, the MazF-mt3-RNA complex is in a non-optimal RNA binding state with a symmetry-related MazF-mt3 molecule found to make interactions with the bound RNA in the crystal. The crystal-packing interactions were further examined by isothermal titration calorimetry (ITC) studies on selected MazF-mt3 mutants. Our attempts to utilize a MazF-mt3 mutant bearing mutations involved in crystal contacts all crystallized with few nucleotides, which are still found to interact with a symmetry mate. However, these different crystal forms revealed the conformational flexibility of loops in the RNA binding interface of MazF-mt3, suggesting their role in RNA binding and recognition, which will require further studies on additional MazF-mt3-RNA complex interactions.

In conclusion, the structures of the MazF-mt1-RNA and MazF-mt3-RNA complexes provide the first structural information on any M. tuberculosis MazF homologues. Supplemented with structure-guided mutational studies on MazF toxicity in vivo, this study has addressed the structural basis of different RNA cleavage specificities among MazF homologues. Our work will guide future studies on the function of other M. tuberculosis MazF and MazE-MazF homologues, and will help delineate their physiological roles in M. tuberculosis stress responses and pathogenesis.

Programming Molecular Devices using Nucleic Acid Hairpins

Nucleic Acid hairpins have been a subject of study for the last four decades. They are composed of single strand that is

hybridized to itself, and the central section forming an unhybridized loop. In nature, they stabilize single stranded RNA, serve as nucleation

sites for RNA folding, protein recognition signals, mRNA localization and regulation of mRNA degradation. On the other hand,

DNA hairpins in biological contexts have been studied with respect to forming cruciform structures that can regulate gene expression.

The use of DNA hairpins as fuel for synthetic molecular devices, including locomotion, was proposed and experimental demonstrated in 2003. They

were interesting because they bring to the table an on-demand energy/information supply mechanism.

The energy/information is hidden (from hybridization) in the hairpin’s loop, until required.

The energy/information is harnessed by opening the stem region, and exposing the single stranded loop section.

The loop region is now free for possible hybridization and help move the system into a thermodynamically favourable state.

The hidden energy and information coupled with

programmability provides another functionality, of selectively choosing what reactions to hide and

what reactions to allow to proceed, that helps develop a topological sequence of events.

Hairpins have been utilized as a source of fuel for many different DNA devices. In this thesis, we program four different

molecular devices using DNA hairpins, and experimentally validate them in the

laboratory. 1) The first device: A

novel enzyme-free autocatalytic self-replicating system composed entirely of DNA that operates isothermally. 2) The second

device: Time-Responsive Circuits using DNA have two properties: a) asynchronous: the final output is always correct

regardless of differences in the arrival time of different inputs.

b) renewable circuits which can be used multiple times without major degradation of the gate motifs

(so if the inputs change over time, the DNA-based circuit can re-compute the output correctly based on the new inputs).

3) The third device: Activatable tiles are a theoretical extension to the Tile assembly model that enhances

its robustness by protecting the sticky sides of tiles until a tile is partially incorporated into a growing assembly.

4) The fourth device: Controlled Amplification of DNA catalytic system: a device such that the amplification

of the system does not run uncontrollably until the system runs out of fuel, but instead achieves a finite

amount of gain.

Nucleic acid circuits with the ability

to perform complex logic operations have many potential practical applications, for example the ability to achieve point of care diagnostics.

We discuss the designs of our DNA Hairpin molecular devices, the results we have obtained, and the challenges we have overcome

to make these truly functional.

Insight into the molecular and functional determinants of HP1043, the essential master regulator of Helicobacter pylori.

Helicobacter pylori is one of the most widespread and successful human pathogens, colonizing half of the population stomach mucosa and causing gastric malignancies in 1% of carriers. Due to the increasing number of antimicrobial-resistant strains, in 2017 the WHO included H. pylori among pathogens that pose a major threat for humankind. In this study, we propose as a molecular target for novel antimicrobial strategies HP1043, an orphan response regulator essential for the viability of H. pylori as it orchestrates all the most important cellular processes. Amino acids most relevant for HP1043 dimerization and target DNA recognition were identified and used to guide an in-silico protein-DNA docking and generate a high-resolution structural model of the interacting HP1043 dimer and its target DNA. The model was experimentally validated and exploited to carry out a virtual screening of small molecule libraries, identifying 8 compounds potentially able to interfere with HP1043 function and likely block H. pylori infection. A second line of research aimed at the characterization of the regulatory function of HP1043 and the tight mechanisms of regulation of hp1043 gene expression. In particular, we proved a direct interaction between HP1043 and the housekeeping sigma80 factor of the RNA polymerase. A conditional mutant H. pylori strain overexpressing a synthetic copy of the hp1043 gene altered in nucleotide sequence yet encoding the wild-type protein was generated, achieving increased intracellular levels of HP1043. However, overexpression of HP1043 did not result in an upregulation of target genes transcription nor modulation of hp1043 transcript levels, pinpointing the existence of multiple overlayed mechanisms of regulation that affect both protein levels and functionality as well as maintain steady the amount of hp1043 transcript. Finally, we proposed that a mechanism of post-transcriptional regulation could depend on an antisense transcript to the hp1043 gene which was validated in two different strains.

Photoactive (supra)molecular devices in double layer membranes

Biological systems are complex and highly organized architectures governed by non-covalent interactions responsible for the regulation of essential tasks in all living organisms. These systems are a constant source of inspiration for supramolecular chemists aiming to design multicomponent molecular assemblies able to perform elaborated tasks, thanks to the role and action of the components that constitute them. Artificial supramolecular systems exploit non-covalent interactions to mimic naturally occurring events. In this context, stimuli-responsive supramolecular systems have attracted attention due to the possibility to control macroscopic effects through modifications at the nanoscale. This thesis is divided in three experimental chapters, characterized by a progressive increase in molecular complexity. Initially, the preparation and studies of liposomes functionalized with a photoactive guest such as azobenzene in the bilayer were tackled, in order to evaluate the effect of such photochrome on the vesicle properties. Subsequently, the synthesis and studies of thread-like molecules comprising an azobenzene functionality was reported. Such molecules were conceived to be intercalated in the bilayer membrane of liposomes with the aim to be used as components for photoresponsive transmembrane molecular pumps. Finally, a [3]rotaxane was developed and studied in solution. This system is composed of two crown ether rings interlocked with an axle containing three recognition sites for the macrocycles, i.e. two pH-switchable ammonium stations and a permanent triazolium station. Such molecule was designed to achieve a change in the ratio between the recognition sites and the crown ethers as a consequence of acid-base inputs. This leads to the formation of rotaxanes containing a number of recognition sites respectively larger, equal or lower than the number of interlocked rings and connected by a network of acid-base reactions.

Molecular And Biochemical Characterization Of Caffeine Synthase And Purine Alkaloid Concentration In Guarana Fruit.

Guarana seeds have the highest caffeine concentration among plants accumulating purine alkaloids, but in contrast with coffee and tea, practically nothing is known about caffeine metabolism in this Amazonian plant. In this study, the levels of purine alkaloids in tissues of five guarana cultivars were determined. Theobromine was the main alkaloid that accumulated in leaves, stems, inflorescences and pericarps of fruit, while caffeine accumulated in the seeds and reached levels from 3.3% to 5.8%. In all tissues analysed, the alkaloid concentration, whether theobromine or caffeine, was higher in young/immature tissues, then decreasing with plant development/maturation. Caffeine synthase activity was highest in seeds of immature fruit. A nucleotide sequence (PcCS) was assembled with sequences retrieved from the EST database REALGENE using sequences of caffeine synthase from coffee and tea, whose expression was also highest in seeds from immature fruit. The PcCS has 1083bp and the protein sequence has greater similarity and identity with the caffeine synthase from cocoa (BTS1) and tea (TCS1). A recombinant PcCS allowed functional characterization of the enzyme as a bifunctional CS, able to catalyse the methylation of 7-methylxanthine to theobromine (3,7-dimethylxanthine), and theobromine to caffeine (1,3,7-trimethylxanthine), respectively. Among several substrates tested, PcCS showed higher affinity for theobromine, differing from all other caffeine synthases described so far, which have higher affinity for paraxanthine. When compared to previous knowledge on the protein structure of coffee caffeine synthase, the unique substrate affinity of PcCS is probably explained by the amino acid residues found in the active site of the predicted protein.

On The Distinct Molecular Architectures Of Dipping- And Spray-lbl Films Containing Lipid Vesicles.

The introduction of spraying procedures to fabricate layer-by-layer (LbL) films has brought new possibilities for the control of molecular architectures and for making the LbL technique compliant with industrial processes. In this study we show that significantly distinct architectures are produced for dipping and spray-LbL films of the same components, which included DODAB/DPPG vesicles. The films differed notably in their thickness and stratified nature. The electrical response of the two types of films to aqueous solutions containing erythrosin was also different. With multidimensional projections we showed that the impedance for the DODAB/DPPG spray-LbL film is more sensitive to changes in concentration, being therefore more promising as sensing units. Furthermore, with surface-enhanced Raman scattering (SERS) we could ascribe the high sensitivity of the LbL films to adsorption of erythrosin.