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.

Differential localization of GABAA receptor subunits in relation to rat striatopallidal and pallidopallidal synapses

As a central integrator of basal ganglia function, the external segment of the globus pallidus (GP) plays a critical role in the control of voluntary movement. The GP is composed of a network of inhibitory GABA-containing projection neurons which receive GABAergic input from axons of the striatum (Str) and local collaterals of GP neurons. Here, using electrophysiological techniques and immunofluorescent labeling we have investigated the differential cellular distribution of a1, a2 and a3 GABAA receptor subunits in relation to striatopallidal (Str-GP) and pallidopallidal (GP-GP) synapses. Electrophysiological investigations showed that zolpidem (100 nm; selective for the a1 subunit) increased the amplitude and the decay time of both Str-GP and GP-GP IPSCs, indicating the presence of the a1 subunits at both synapses. However, the application of drugs selective for the a2, a3 and a5 subunits (zolpidem at 400 nm, L-838,417 and TP003) revealed differential effects on amplitude and decay time of IPSCs, suggesting the nonuniform distribution of non-a1 subunits. Immunofluorescence revealed widespread distribution of the a1 subunit at both soma and dendrites, while double- and triple-immunofluorescent labeling for parvalbumin, enkephalin, gephyrin and the ?2 subunit indicated strong immunoreactivity for GABAAa3 subunits in perisomatic synapses, a region mainly targeted by local axon collaterals. In contrast, immunoreactivity for synaptic GABAAa2 subunits was observed in dendritic compartments where striatal synapses are preferentially located. Due to the kinetic properties which each GABAAa subunit confers, this distribution is likely to contribute differentially to both physiological and pathological patterns of activity.

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.

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.