Toward an atomistic understanding of the immune synapse:large-scale molecular dynamics simulation of a membrane-embedded TCR-pMHC-CD4 complex

T-cell activation requires interaction of T-cell receptors (TCR) with peptide epitopes bound by major histocompatibility complex (MHC) proteins. This interaction occurs at a special cell-cell junction known as the immune or immunological synapse. Fluorescence microscopy has shown that the interplay among one agonist peptide-MHC (pMHC), one TCR and one CD4 provides the minimum complexity needed to trigger transient calcium signalling. We describe a computational approach to the study of the immune synapse. Using molecular dynamics simulation, we report here on a study of the smallest viable model, a TCR-pMHC-CD4 complex in a membrane environment. The computed structural and thermodynamic properties are in fair agreement with experiment. A number of biomolecules participate in the formation of the immunological synapse. Multi-scale molecular dynamics simulations may be the best opportunity we have to reach a full understanding of this remarkable supra-macromolecular event at a cell-cell junction.

Close contact fluctuations:the seeding of signalling domains in the immunological synapse

We analyse the size and density of thermally induced regions of close contact in cell : cell contact interfaces within a harmonic potential approximation, estimating these regions to be below one-tenth of a micron across. Our calculations indicate that as the distance between the close contact threshold depth and the mean membrane-membrane separation increases, the density of close contact patches decreases exponentially while there is only a minimal variation in their mean size. The technique developed can be used to calculate the probability of first crossing in reflection symmetry violating systems. © Europhysics Letters Association.

Contact time periods in immunological synapse

This paper resolves the long standing debate as to the proper time scale τ of the onset of the immunological synapse bond, the noncovalent chemical bond defining the immune pathways involving T cells and antigen presenting cells. Results from our model calculations show τ to be of the order of seconds instead of minutes. Close to the linearly stable regime, we show that in between the two critical spatial thresholds defined by the integrin:ligand pair (Δ2∼ 40-45 nm) and the T-cell receptor TCR:peptide-major-histocompatibility-complex pMHC bond (Δ1∼ 14-15 nm), τ grows monotonically with increasing coreceptor bond length separation δ (= Δ2-Δ1∼ 26-30 nm) while τ decays with Δ1 for fixed Δ2. The nonuniversal δ-dependent power-law structure of the probability density function further explains why only the TCR:pMHC bond is a likely candidate to form a stable synapse.

Temporal dynamics in an immunological synapse:role of thermal fluctuations in signaling

The article analyzes the contribution of stochastic thermal fluctuations in the attachment times of the immature T-cell receptor TCR: peptide-major-histocompatibility-complex pMHC immunological synapse bond. The key question addressed here is the following: how does a synapse bond remain stabilized in the presence of high-frequency thermal noise that potentially equates to a strong detaching force? Focusing on the average time persistence of an immature synapse, we show that the high-frequency nodes accompanying large fluctuations are counterbalanced by low-frequency nodes that evolve over longer time periods, eventually leading to signaling of the immunological synapse bond primarily decided by nodes of the latter type. Our analysis shows that such a counterintuitive behavior could be easily explained from the fact that the survival probability distribution is governed by two distinct phases, corresponding to two separate time exponents, for the two different time regimes. The relatively shorter timescales correspond to the cohesion:adhesion induced immature bond formation whereas the larger time reciprocates the association:dissociation regime leading to TCR:pMHC signaling. From an estimate of the bond survival probability, we show that, at shorter timescales, this probability PΔ(τ) scales with time τ as a universal function of a rescaled noise amplitude DΔ2, such that PΔ(τ)∼τ-(ΔD+12),Δ being the distance from the mean intermembrane (T cell:Antigen Presenting Cell) separation distance. The crossover from this shorter to a longer time regime leads to a universality in the dynamics, at which point the survival probability shows a different power-law scaling compared to the one at shorter timescales. In biological terms, such a crossover indicates that the TCR:pMHC bond has a survival probability with a slower decay rate than the longer LFA-1:ICAM-1 bond justifying its stability.

Immunological synapse: a mathematical model of the bond formation process:mathematical modelling of the immature synapse

The cell:cell bond between an immune cell and an antigen presenting cell is a necessary event in the activation of the adaptive immune response. At the juncture between the cells, cell surface molecules on the opposing cells form non-covalent bonds and a distinct patterning is observed that is termed the immunological synapse. An important binding molecule in the synapse is the T-cell receptor (TCR), that is responsible for antigen recognition through its binding with a major-histocompatibility complex with bound peptide (pMHC). This bond leads to intracellular signalling events that culminate in the activation of the T-cell, and ultimately leads to the expression of the immune eector function. The temporal analysis of the TCR bonds during the formation of the immunological synapse presents a problem to biologists, due to the spatio-temporal scales (nanometers and picoseconds) that compare with experimental uncertainty limits. In this study, a linear stochastic model, derived from a nonlinear model of the synapse, is used to analyse the temporal dynamics of the bond attachments for the TCR. Mathematical analysis and numerical methods are employed to analyse the qualitative dynamics of the nonequilibrium membrane dynamics, with the specic aim of calculating the average persistence time for the TCR:pMHC bond. A single-threshold method, that has been previously used to successfully calculate the TCR:pMHC contact path sizes in the synapse, is applied to produce results for the average contact times of the TCR:pMHC bonds. This method is extended through the development of a two-threshold method, that produces results suggesting the average time persistence for the TCR:pMHC bond is in the order of 2-4 seconds, values that agree with experimental evidence for TCR signalling. The study reveals two distinct scaling regimes in the time persistent survival probability density prole of these bonds, one dominated by thermal uctuations and the other associated with the TCR signalling. Analysis of the thermal fluctuation regime reveals a minimal contribution to the average time persistence calculation, that has an important biological implication when comparing the probabilistic models to experimental evidence. In cases where only a few statistics can be gathered from experimental conditions, the results are unlikely to match the probabilistic predictions. The results also identify a rescaling relationship between the thermal noise and the bond length, suggesting a recalibration of the experimental conditions, to adhere to this scaling relationship, will enable biologists to identify the start of the signalling regime for previously unobserved receptor:ligand bonds. Also, the regime associated with TCR signalling exhibits a universal decay rate for the persistence probability, that is independent of the bond length.

Antioxidants and protein oxidation

Proteins are susceptible to oxidation by reactive oxygen species, where the type of damage induced is characteristic of the denaturing species. The induction of protein carbonyls is a widely applied biomarker, arising from primary oxidative insult. However, when applied to complex biological and pathological conditions it can be subject to interference from lipid, carbohydrate and DNA oxidation products. More recently, interest has focused on the analysis of specific protein bound oxidised amino acids. Of the 22 amino acids, aromatic and sulphydryl containing residues have been regarded as being particularly susceptible to oxidative modification, with L-DOPA from tyrosine, ortho-tyrosine from phenylalanine; sulphoxides and disulphides from methionine and cysteine respectively; and kynurenines from tryptophan. Latterly, the identification of valine and leucine hydroxides, reduced from hydroperoxide intermediates, has been described and applied. In order to examine the nature of oxidative damage and protective efficacy of antioxidants the markers must be thoroughly evaluated for dosimetry in vitro following damage by specific radical species. Antioxidant protection against formation of the biomarker should be demonstrated in vitro. Quantification of biomarkers in proteins from normal subjects should be within the limits of detection of any analytical procedure. Further to this, the techniques for isolation and hydrolysis of specific proteins should demonstrate that in vitro oxidation is minimised. There is a need for the development of standards for quality assurance material to standardise procedures between laboratories. At present, antioxidant effects on protein oxidation in vivo are limited to animal studies, where dietary antioxidants have been reported to reduce dityrosine formation during rat exercise training. Two studies on humans have been reported last year. The further application of these methods to human studies is indicated, where the quality of the determinations will be enhanced through inter-laboratory validation.

Studies of clay minerals and their decomposition products

Mõssbauer spectroscopy and X-ray diffraction of five coals revealed the presence of pyrite, illite, kaolinite and Quartz, together with other minor phases. Analysis of the coal ashes indicated the formation of hematite and an Fe (3+) paramagnetic phase, the latter resulting from .the dehydroxylation of the clay minerals during ashing at 700 to 750 C. By using a combination of several physicochemical methods, different successive stages of dehydroxylation, structural consolidation, and recrystallisation of illite, montmorillonite and hectorite upon thermal treatment to 1300 C were investigated. Dehydroxylation of the clay minerals occurred between 450 and 750 C, the X-ray crysdallinity of illite and montmorillonite remaining until 800 C. Hectorite gradually recrystallises to enstatite at temperatures above 700°C. At 900 C the crystalline structure of all three clay minerals had totally collapsed. Solid state reactions occurred above 900 C producing such phases as spinel, hematite, enstatite, cristobalite and mullite. Illite and montmorillonite started to melt between 1200 and 1300°C, producing a silicate glass that contained Fe(3+) and Fe(2+) ions. Ortho-pnstatite, clino-enstatite and proto-enstatite were identified in the thermal products of hectorite, their relative proportions varying with temperature. Protoenstatite was stabilised with respect to metastable clinoenstatite upon cooling from 12000 C by the presence of exchanged transition metal cations. Solid state Nuclear Magnetic Resonance spectroscopy of thermally treated transition metal exchanged hectorite indicated the levels at which paramagnetic cations could be loaded on to the clay before spectral resolution is significantly diminished.

On the mechanisms of cerebellar synaptic plasticity

Changes in the strength of signalling between neurones are thought to provide a cellular substrate for learning and memory. In the cerebellar cortex, raising the frequency and the strength of parallel fibre (PF) stimulation leads to a long-term depression (LTD) of the strength of signalling at the synapse between PFs and Purkinje cells (PCs), which spreads to distant synapses to the same cell via a nitric oxide (NO) dependent mechanism. At the same synapse, but under conditions of reduced post-synaptic calcium activity, raised frequency stimulation (RFS) of PFs triggers a long-term potentiation of synaptic transmission. The aims of the work described in this thesis were to investigate the conditions necessary for LTD and LTP at this synapse following RFS and to identify the origins and second messenger cascades involved in the induction and spread of LTP and LTD. In thin, parasagittal cerebellar slices whole cell patch clamp recordings were made from PCs and the effects of RFS of one of two, independent PF inputs to the same PC were examined under a range of experimental conditions. Under conditions designed to reduce post-synaptic calcium activity, RFS to a single PF input led to LTP and a decreases in paired pulse facilitation (PPF) in both pathways. This heterosynaptic potentiation was prevented by inhibition of protein kinase A (PKA) or by inhibition of NO synthase with either 7-nitroindazole (7-NI) or NG Nitro-L-argenine methyl ester. Inhibition of guanylate cyclase (GC) or protein kinase G (PKG) had no effect. A similar potentiation was observed upon application of the adenylyl cyclase (AC) activator forskolin or the NO donor spermine NONOate. Both of these treatments also resulted in an increase in the frequency of mEPSCs, which provides further evidence for a presynaptic origin of LTP. Forskolin induced potentiation and the increase in mEPSC frequency were blocked by 7-NI. The styryl dye FM1-43, a fluorescent reporter of endo- and exocytosis, was also used to further examine the possible pre-synaptic origins of LTP. RFS or forskolin application enhanced FM1-43 de-staining and NOS inhibitors blocked this effect. Application of NONOate also enhanced FM1-43 de-staining. When post-synaptic calcium activity was less strictly buffered, RFS to a single PF input led to a transient potentiation that was succeeded by LTD in both pathways. This LTD, which resembled previously described forms, was prevented by inhibition of the NO/cGMP/PKG cascade. Modification of the AC/cAMP/PKA cascade had no effect. In summary, the direction of synaptic plasticity at the PF-PC synapse in response to RFS depends largely on the level of post-synaptic calcium activity. LTP and LTD were non-input specific and both forms of plasticity were dependent on NOS activity. Induction of LTP was mediated by a presynaptic mechanism and depended on NO and cAMP production. LTD on the other hand was a post-synaptic process and required activity of the NO/cGMP/PKG signalling cascade.

Non-classical inhibitors of dihydrofolate reductase

This thesis comprises two main objectives. The first objective involved the stereochemical studies of chiral 4,6-diamino-1-aryl-1,2-dihydro-s-triazines and an investigation on how the different conformations of these stereoisomers may affect their binding affinity to the enzyme dihydrofolate reductase (DHFR). The ortho-substituted 1-aryl-1,2-dihydro-s-triazines were synthesised by the three component method. An ortho-substitution at the C6' position was observed when meta-azidocycloguanil was decomposed in acid. The ortho-substituent restricts free rotation and this gives rise to atropisomerism. Ortho-substituted 4,6-diamino-1-aryl-2-ethyl-1,2-dihydro-2-methyl-s-triazine contains two elements of chirality and therefore exists as four stereoisomers: (S,aR), (R,aS), (R,aR) and (S,aS). The energy barriers to rotation of these compounds were calculated by a semi-empirical molecular orbital program called MOPAC and they were found to be in excess of 23 kcal/mol. The diastereoisomers were resolved and enriched by C18 reversed phase h.p.l.c. Nuclear overhauser effect experiments revealed that (S,aR) and (R,aS) were the more stable pair of stereoisomers and therefore existed as the major component. The minor diastereoisomers showed greater binding affinity for the rat liver DHFR in in vitro assay. The second objective entailed the investigation into the possibility of retaining DHFR inhibitory activity by replacing the classical diamino heterocyclic moiety with an amidinyl group. 4-Benzylamino-3-nitro-N,N-dimethyl-phenylamidine was synthesised in two steps. One of the two phenylamidines indicated weak inhibition against the rat liver DHFR. This weak activity may be due to the failure of the inhibitor molecule to form strong hydrogen bonds with residue Glu-30 at the active site of the enzyme.

T-cell epitope prediction and immune complex simulation using molecular dynamics:state of the art and persisting challenges

Atomistic Molecular Dynamics provides powerful and flexible tools for the prediction and analysis of molecular and macromolecular systems. Specifically, it provides a means by which we can measure theoretically that which cannot be measured experimentally: the dynamic time-evolution of complex systems comprising atoms and molecules. It is particularly suitable for the simulation and analysis of the otherwise inaccessible details of MHC-peptide interaction and, on a larger scale, the simulation of the immune synapse. Progress has been relatively tentative yet the emergence of truly high-performance computing and the development of coarse-grained simulation now offers us the hope of accurately predicting thermodynamic parameters and of simulating not merely a handful of proteins but larger, longer simulations comprising thousands of protein molecules and the cellular scale structures they form. We exemplify this within the context of immunoinformatics.

Glutamatergic input-output properties of thalamic astrocytes

Astrocytes in the somatosensory ventrobasal (VB) thalamus of rats respond to glutamatergic synaptic input with metabotropic glutamate receptor (mGluR) mediated intracellular calcium ([Ca²?](i)) elevations. Astrocytes in the VB thalamus also release the gliotransmitter (GT) glutamate in a Ca²?-dependent manner. The tripartite synapse hypothesis posits that astrocytic [Ca²?](i) elevations resulting from synaptic input releases gliotransmitters that then feedback to modify the synapse. Understanding the dynamics of this process and the conditions under which it occurs are therefore important steps in elucidating the potential roles and impact of GT release in particular brain activities. In this study, we investigated the relationship between VB thalamus afferent synaptic input and astrocytic glutamate release by recording N-methyl-D-aspartate (NMDA) receptor-mediated slow inward currents (SICs) elicited in neighboring neurons. We found that Lemniscal or cortical afferent stimulation, which can elicit astrocytic [Ca²?](i) elevations, do not typically result in the generation of SICs in thalamocortical (TC) neurons. Rather, we find that the spontaneous emergence of SICs is largely resistant to acute afferent input. The frequency of SICs, however, is correlated to long-lasting afferent activity. In contrast to short-term stimulus-evoked GT release effects reported in other brain areas, astrocytes in the VB thalamus do not express a straightforward input-output relationship for SIC generation but exhibit integrative characteristics.

Proteomic analysis of a noninvasive human model of acute inflammation and its resolution:the twenty-one day gingivitis model

The 21-day experimental gingivitis model, an established noninvasive model of inflammation in response to increasing bacterial accumulation in humans, is designed to enable the study of both the induction and resolution of inflammation. Here, we have analyzed gingival crevicular fluid, an oral fluid comprising a serum transudate and tissue exudates, by LC-MS/MS using Fourier transform ion cyclotron resonance mass spectrometry and iTRAQ isobaric mass tags, to establish meta-proteomic profiles of inflammation-induced changes in proteins in healthy young volunteers. Across the course of experimentally induced gingivitis, we identified 16 bacterial and 186 human proteins. Although abundances of the bacterial proteins identified did not vary temporally, Fusobacterium outer membrane proteins were detected. Fusobacterium species have previously been associated with periodontal health or disease. The human proteins identified spanned a wide range of compartments (both extracellular and intracellular) and functions, including serum proteins, proteins displaying antibacterial properties, and proteins with functions associated with cellular transcription, DNA binding, the cytoskeleton, cell adhesion, and cilia. PolySNAP3 clustering software was used in a multilayered analytical approach. Clusters of proteins that associated with changes to the clinical parameters included neuronal and synapse associated proteins.

Astrocyte plasticity:implications for synaptic and neuronal activity

Astrocytes are increasingly implicated in a range of functions in the brain, many of which were previously ascribed to neurons. Much of the prevailing interest centers on the role of astrocytes in the modulation of synaptic transmission and their involvement in the induction of forms of plasticity such as long-term potentiation and long-term depression. However, there is also an increasing realization that astrocytes themselves can undergo plasticity. This plasticity may be manifest as changes in protein expression which may modify calcium activity within the cells, changes in morphology that affect the environment of the synapse and the extracellular space, or changes in gap junction astrocyte coupling that modify the transfer of ions and metabolites through astrocyte networks. Plasticity in the way that astrocytes release gliotransmitters can also have direct effects on synaptic activity and neuronal excitability. Astrocyte plasticity can potentially have profound effects on neuronal network activity and be recruited in pathological conditions. An emerging principle of astrocyte plasticity is that it is often induced by neuronal activity, reinforcing our emerging understanding of the working brain as a constant interaction between neurons and glial cells. © The Author(s) 2013.

Direct observation of morphological development during the spin-coating of polystyrene-poly(methyl methacrylate) polymer blends

We present results of the direct observation, in real-space, of the phase separation of high molecular weight polystyrene and poly(methyl methacrylate) from ortho-xylene using our newly developed technique of high speed stroboscopic interference microscopy. Taking a fixed concentration (3 wt % in o-xylene) at a fixed composition (1:4 by weight) and by varying the rotational rate during the spin-coating process, we are able to observe the formation of a range of phase separated bicontinuous morphologies of differing length-scales. Importantly, we are able to show that the mechanism by which the final phase separated structure is formed is through domain coarsening when rich in solvent, before vitrification occurs and fixes the phase separated structure. The ability to directly observe morphological development offers a route toward controlling the length-scale of the final morphology through process control and in situ feedback, from a single stock solution. © 2013 Wiley Periodicals, Inc.