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.

Why are MD simulated protein folding times wrong?

The question of significant deviations of protein folding times simulated using molecular dynamics from experimental values is investigated. It is shown that in the framework of Markov State Model (MSM) describing the conformational dynamics of peptides and proteins, the folding time is very sensitive to the simulation model parameters, such as forcefield and temperature. Using two peptides as examples, we show that the deviations in the folding times can reach an order of magnitude for modest variations of the molecular model. We, therefore, conclude that the folding rate values obtained in molecular dynamics simulations have to be treated with care.

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.

G-protein-coupled receptors:from structural insights to functional mechanisms

The papers resulting from the recent Biochemical Society Focused Meeting 'G-Protein-Coupled Receptors: from Structural Insights to Functional Mechanisms' held in Prato in September 2012 are introduced in the present overview. A number of future goals for GPCR (G-protein-coupled receptor) research are considered, including the need to develop biophysical and computational methods to explore the full range of GPCR conformations and their dynamics, the need to develop methods to take this into account for drug discovery and the importance of relating observations on isolated receptors or receptors expressed in model systems to receptor function in vivo. © 2013 Biochemical Society.

Tear analysis and lens-tear interactions:Part I. Protein fingerprinting with microfluidic technology

The purpose of this work is to establish the application of a fully automated microfluidic chip based protein separation assay in tear analysis. It is rapid, requires small sample volumes and is vastly superior to, and more convenient than, comparable conventional gel electrophoresis assays. The protein sizing chip technology was applied to three specific fields of analysis. Firstly tear samples were collected regularly from subjects establishing the baseline effects of tear stimulation, tear state and patient health. Secondly tear samples were taken from lens wearing eyes and thirdly the use of microfluidic technology was assessed as a means to investigate a novel area of tear analysis, which we have termed the 'tear envelope'. Utilising the Agilent 2100 Bioanalyzer in combination with the Protein 200 Plus LabChip kit, these studies investigated tear proteins in the range of 14-200 kDa. Particular attention was paid to the relative concentrations of lysozyme, tear lipocalin, secretory IgA (sIgA), IgG and lactoferrin, together with the overall tear electropherogram 'fingerprint'. Furthermore, whilst lens-tear interaction studies are generally thought of as an investigation into the effects of tears components on the contact lens material, i.e. deposition studies, this report addresses the reverse phenomenon-the effect of the lens, and particularly the newly inserted lens, on the tear fluid composition and dynamics. The use of microfluidic technology provides a significant advance in tear studies and should prove invaluable in tear diagnostics and contact lens performance analysis.

Influence of downstream processing on the breakage of whey protein precipitates

Whey proteins may be fractionated by isoelectric precipitation followed by centrifugal recovery of the precipitate phase. Transport and processing of protein precipitates may alter the precipitate particle properties, which may affect how they behave in subsequent processes. For example, the transport of precipitate solution through pumps, pipes and valves and into a centrifugal separator may cause changes in particle size and density, which may affect the performance of the separator. This work investigates the effect of fluid flow intensity, flow geometry and exposure time on the breakage of whey protein precipitates: Computational fluid dynamics (CFD) was used to quantify the flow intensity in different geometries. Flow geometry can have a critical impact on particle breakage. Sharp geometrical transitions induce large increases in turbulence that can result in substantial particle breakage. As protein precipitate particles break, they tend to form denser more compact structures. The reduction in particle size and increase in compaction is due to breakage. This makes the particles become more resistant to further breakage as particle compactness increases. The effect of flow intensity on particle breakage is coupled to exposure time, with greater exposure time producing more breakage. However, it is expected that the particles will attain an equilibrium particle size and density after prolonged exposure in a constant flow field where no further breakage will occur with exposure time. © 2005 Institution of Chemical Engineers.

Large-scale molecular dynamics simulations of HLA-A*0201 complexed with a tumor-specific antigenic peptide:can the α3 and β2m domains be neglected?

Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histo-compatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVY-DGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available Large-scale massively parallel molecular dynamics (MD) simulations of the human class I major histocompatibility complex (MHC) protein HLA-A*0201 bound to a decameric tumor-specific antigenic peptide GVYDGREHTV were performed using a scalable MD code on high-performance computing platforms. Such computational capabilities put us in reach of simulations of various scales and complexities. The supercomputing resources available for this study allow us to compare directly differences in the behavior of very large molecular models; in this case, the entire extracellular portion of the peptide–MHC complex vs. the isolated peptide binding domain. Comparison of the results from the partial and the whole system simulations indicates that the peptide is less tightly bound in the partial system than in the whole system. From a detailed study of conformations, solvent-accessible surface area, the nature of the water network structure, and the binding energies, we conclude that, when considering the conformation of the α1–α2 domain, the α3 and β2m domains cannot be neglected. © 2004 Wiley Periodicals, Inc. J Comput Chem 25: 1803–1813, 2004

A proteo-liposome system for the analysis of the intracellular interactome of membrane proteins using amyloid precursor protein as a model

Transmembrane proteins play crucial roles in many important physiological processes. The intracellular domain of membrane proteins is key for their function by interacting with a wide variety of cytosolic proteins. It is therefore important to examine this interaction. A recently developed method to study these interactions, based on the use of liposomes as a model membrane, involves the covalent coupling of the cytoplasmic domains of membrane proteins to the liposome membrane. This allows for the analysis of interaction partners requiring both protein and membrane lipid binding. This thesis further establishes the liposome recruitment system and utilises it to examine the intracellular interactome of the amyloid precursor protein (APP), most well-known for its proteolytic cleavage that results in the production and accumulation of amyloid beta fragments, the main constituent of amyloid plaques in Alzheimer’s disease pathology. Despite this, the physiological function of APP remains largely unclear. Through the use of the proteo-liposome recruitment system two novel interactions of APP’s intracellular domain (AICD) are examined with a view to gaining a greater insight into APP’s physiological function. One of these novel interactions is between AICD and the mTOR complex, a serine/threonine protein kinase that integrates signals from nutrients and growth factors. The kinase domain of mTOR directly binds to AICD and the N-terminal amino acids of AICD are crucial for this interaction. The second novel interaction is between AICD and the endosomal PIKfyve complex, a lipid kinase involved in the production of phosphatidylinositol-3,5-bisphosphate (PI(3,5)P2) from phosphatidylinositol-3-phosphate, which has a role in controlling ensdosome dynamics. The scaffold protein Vac14 of the PIKfyve complex binds directly to AICD and the C-terminus of AICD is important for its interaction with the PIKfyve complex. Using a recently developed intracellular PI(3,5)P2 probe it is shown that APP controls the formation of PI(3,5)P2 positive vesicular structures and that the PIKfyve complex is involved in the trafficking and degradation of APP. Both of these novel APP interactors have important implications of both APP function and Alzheimer’s disease. The proteo-liposome recruitment method is further validated through its use to examine the recruitment and assembly of the AP-2/clathrin coat from purified components to two membrane proteins containing different sorting motifs. Taken together this thesis highlights the proteo-liposome recruitment system as a valuable tool for the study of membrane proteins intracellular interactome. It allows for the mimicking of the protein in its native configuration therefore identifying weaker interactions that are not detected by more conventional methods and also detecting interactions that are mediated by membrane phospholipids.

Histidine hydrogen bonding in MHC at pH 5 and pH 7 modeled by molecular docking and molecular dynamics simulations

Hydrogen bonds play important roles in maintaining the structure of proteins and in the formation of most biomolecular protein-ligand complexes. All amino acids can act as hydrogen bond donors and acceptors. Among amino acids, Histidine is unique, as it can exist in neutral or positively charged forms within the physiological pH range of 5.0 to 7.0. Histidine can thus interact with other aromatic residues as well as forming hydrogen bonds with polar and charged residues. The ability of His to exchange a proton lies at the heart of many important functional biomolecular interactions, including immunological ones. By using molecular docking and molecular dynamics simulation, we examine the influence of His protonation/deprotonation on peptide binding affinity to MHC class II proteins from locus HLA-DP. Peptide-MHC interaction underlies the adaptive cellular immune response, upon which the next generation of commercially-important vaccines will depend. Consistent with experiment, we find that peptides containing protonated His residues bind better to HLA-DP proteins than those with unprotonated His. Enhanced binding at pH 5.0 is due, in part, to additional hydrogen bonds formed between peptide His+ and DP proteins. In acidic endosomes, protein His79β is predominantly protonated. As a result, the peptide binding cleft narrows in the vicinity of His79β, which stabilizes the peptide - HLA-DP protein complex. © 2014 Bentham Science Publishers.

The amyloid precursor protein (APP) binds the PIKfyve complex and modulates its function

Phosphoinositides are important components of eukaryotic membranes that are required for multiple forms of membrane dynamics. Phosphoinositides are involved in defining membrane identity, mediate cell signalling and control membrane trafficking events. Due to their pivotal role in membrane dynamics, phosphoinositide de-regulation contributes to various human diseases. In this review, we will focus on the newly emerging regulation of the PIKfyve complex, a phosphoinositide kinase that converts the endosomal phosphatidylinositol-3-phosphate [PI(3)P] to phosphatidylinositol-3,5-bisphosphate [PI(3,5)P2)], a low abundance phosphoinositide of outstanding importance for neuronal integrity and function. Loss of PIKfyve function is well known to result in neurodegeneration in both mousemodels and human patients. Our recent work has surprisingly identified the amyloid precursor protein (APP), the central molecule in Alzheimer s disease aetiology, as a novel interaction partner of a subunit of the PIKfyve complex, Vac14. Furthermore, it has been shown that APP modulates PIKfyve function and PI(3,5)P2 dynamics, suggesting that the APP gene family functions as regulator of PI(3,5)P2 metabolism. The recent advances discussed in this review suggest a novel, unexpected, â-amyloid-independent mechanism for neurodegeneration in Alzheimer s disease.

Quantifying protein folding pathways

A deep understanding of the proteins folding dynamics can be get quantifying folding landscape by calculating how the number of microscopic configurations (entropy) varies with the energy of the chain, Ω=Ω(E). Because of the incredibly large number of microstates available to a protein, direct enumeration of Ω(E) is not possible on realistic computer simulations. An estimate of Ω(E) can be obtained by use of a combination of statistical mechanics and thermodynamics. By combining different definitions of entropy that are valid for a system whose probability for occupying a state is given by the canonical Boltzmann probability, computers allow the determination of Ω(E). ^ The energy landscapes of two similar, but not identical model proteins were studied. One protein contains no kinetic tracks. Results show a smooth funnel for the folding landscape. That allows the contour determination of the folding funnel. Also it was presented results for the folding landscape for a modified protein with kinetic traps. Final results show that the computational approach is able to distinguish and explore regions of the folding landscape that are due to kinetic traps from the native state folding funnel.^

Conformational dynamics associated with ligand binding to vertebrate hexa-coordinate hemoglobins

Neuroglobin (Ngb) and cytoglobin (Cygb) are two new additions to the globin family, exhibiting heme iron hexa-coordination, a disulfide bond and large internal cavities. These proteins are implicated in cytoprotection under hypoxic-ischemic conditions, but the molecular basis of their cytoprotective function is unclear. Herein, a photothermal and spectroscopic study of the interactions of diatomic ligands with Ngb, Cygb, myoglobin and hemoglobin is presented. The impact of the disulfide bond in Ngb and Cygb and role of conserved residues in Ngb His64, Val68, Cys55, Cys120 and Tyr44 on conformational dynamics associated with ligand binding/dissociation were investigated. Transient absorption and photoacoustic calorimetry studies indicate that CO photo-dissociation from Ngb leads to a volume expansion (13.4±0.9 mL mol-1), whereas a smaller volume change was determined for Ngb with reduced Cys (ΔV=4.6±0.3 mL mol-1). Furthermore, Val68 side chain regulates ligand migration between the distal pocket and internal hydrophobic cavities since Val68Phe geminate quantum yield is ∼2.7 times larger than that of WT Ngb. His64Gln and Tyr44Phe mutations alter the thermodynamic parameters associated with CO photo-release indicating that electrostatic/hydrogen binding network that includes heme propionate groups, Lys 67, His64, and Tyr 44 in Ngb modulates the energetics of CO photo-dissociation. In Cygb, CO escape from the protein matrix is fast (< 40 ns) with a ΔH of 18±2 kcal mol-1 in Cygbred, whereas disulfide bridge formation promotes a biphasic ligand escape associated with an overall enthalpy change of 9±4 kcal mol-1. Therefore, the disulfide bond modulates conformational dynamics in Ngb and Cygb. I propose that in Cygb with reduced Cys the photo-dissociated ligand escapes through the hydrophobic tunnel as occurs in Ngb, whereas the CO preferentially migrates through the His64 gate in Cygbox. To characterize Cygb surface 1,8-ANS interactions with Cygb were investigated employing fluorescence spectroscopy, ITC and docking simulations. Two 1,8-ANS binding sites were identified. One binding site is located close to the extended N-terminus of Cygb and was also identified as a binding site for oleate. Furthermore, guanidinium hydrochloride-induced unfolding studies of Cygb reveal that the disulfide bond does not impact Cygb stability, whereas binding of cyanide slightly increases the protein stability.

On Provable Algorithms for Determination of Continuous Protein Interdomain Motions from Residual Dipolar Couplings

Dynamics of biomolecules over various spatial and time scales are essential for biological functions such as molecular recognition, catalysis and signaling. However, reconstruction of biomolecular dynamics from experimental observables requires the determination of a conformational probability distribution. Unfortunately, these distributions cannot be fully constrained by the limited information from experiments, making the problem an ill-posed one in the terminology of Hadamard. The ill-posed nature of the problem comes from the fact that it has no unique solution. Multiple or even an infinite number of solutions may exist. To avoid the ill-posed nature, the problem needs to be regularized by making assumptions, which inevitably introduce biases into the result.

Here, I present two continuous probability density function approaches to solve an important inverse problem called the RDC trigonometric moment problem. By focusing on interdomain orientations we reduced the problem to determination of a distribution on the 3D rotational space from residual dipolar couplings (RDCs). We derived an analytical equation that relates alignment tensors of adjacent domains, which serves as the foundation of the two methods. In the first approach, the ill-posed nature of the problem was avoided by introducing a continuous distribution model, which enjoys a smoothness assumption. To find the optimal solution for the distribution, we also designed an efficient branch-and-bound algorithm that exploits the mathematical structure of the analytical solutions. The algorithm is guaranteed to find the distribution that best satisfies the analytical relationship. We observed good performance of the method when tested under various levels of experimental noise and when applied to two protein systems. The second approach avoids the use of any model by employing maximum entropy principles. This 'model-free' approach delivers the least biased result which presents our state of knowledge. In this approach, the solution is an exponential function of Lagrange multipliers. To determine the multipliers, a convex objective function is constructed. Consequently, the maximum entropy solution can be found easily by gradient descent methods. Both algorithms can be applied to biomolecular RDC data in general, including data from RNA and DNA molecules.

Spatiotemporal Dynamics of CaMKI During Structural Plasticity of Single Dendritic Spines

Multifunctional calcium/calmodulin dependent protein kinases (CaMKs) are key regulators of spine structural plasticity and long-term potentiation (LTP) in neurons. CaMKs have promiscuous and overlapping substrate recognition motifs, and are distinguished in their regulatory role based on differences in the spatiotemporal dynamics of activity. While the function and activity of CaMKII in synaptic plasticity has been extensively studied, that of CaMKI, another major class of CaMK required for LTP, still remain elusive.

Here, we develop a Förster’s Resonance Energy Transfer (FRET) based sensor to measure the spatiotemporal activity dynamics of CaMK1. We monitored CaMKI activity using 2-photon fluorescence lifetime imaging, while inducing LTP in single dendritic spines of rat (Rattus Norvegicus, strain Sprague Dawley) hippocampal CA1 pyramidal neurons using 2-photon glutamate uncaging. Using RNA-interference and pharmacological means, we also characterize the role of CaMKI during spine structural plasticity.

We found that CaMKI was rapidly and transiently activated with a rise time of ~0.3 s and decay time of ~1 s in response to each uncaging pulse. Activity of CaMKI spread out of the spine. Phosphorylation of CaMKI by CaMKK was required for this spreading and for the initial phase of structural LTP. Combined with previous data showing that CaMKII is restricted to the stimulated spine and required for long-term maintenance of structural LTP, these results suggest that CaMK diversity allows the same incoming signal – calcium – to independently regulate distinct phases of LTP by activating different CaMKs with distinct spatiotemporal dynamics.

Role of Vinculin in Regulating Force-Sensitive Dynamics of Adherens Junctions During Collective Cell Migration

Dynamic processes such as morphogenesis and tissue patterning require the precise control of many cellular processes, especially cell migration. Historically, these processes are thought to be mediated by genetic and biochemical signaling pathways. However, recent advances have unraveled a previously unappreciated role of mechanical forces in regulating these homeostatic processes in of multicellular systems. In multicellular systems cells adhere to both deformable extracellular matrix (ECM) and other cells, which are sources of applied forces and means of mechanical support. Cells detect and respond to these mechanical signals through a poorly understood process called mechanotransduction, which can have profound effects on processes such as cell migration. These effects are largely mediated by the sub cellular structures that link cells to the ECM, called focal adhesions (FAs), or cells to other cells, termed adherens junctions (AJs).

Overall this thesis is comprised of my work on identifying a novel force dependent function of vinculin, a protein which resides in both FAs and AJs - in dynamic process of collective migration. Using a collective migration assay as a model for collective cell behavior and a fluorescence resonance energy transfer (FRET) based molecular tension sensor for vinculin I demonstrated a spatial gradient of tension across vinculin in the direction of migration. To define this novel force-dependent role of vinculin in collective migration I took advantage of previously established shRNA based vinculin knock down Marin-Darby Canine Kidney (MDCK) epithelial cells.

The first part of my thesis comprises of my work demonstrating the mechanosensitive role of vinculin at AJ’s in collectively migrating cells. Using vinculin knockdown cells and vinculin mutants, which specifically disrupt vinculin’s ability to bind actin (VinI997A) or disrupt its ability to localize to AJs without affecting its localization at FAs (VinY822F), I establish a role of force across vinculin in E-cadherin internalization and clipping. Furthermore by measuring E-cadherin dynamics using fluorescence recovery after bleaching (FRAP) analysis I show that vinculin inhibition affects the turnover of E-cadherin at AJs. Together these data reveal a novel mechanosensitive role of vinculin in E-cadherin internalization and turnover in a migrating cell layer, which is contrary to the previously identified role of vinculin in potentiating E-cadherin junctions in a static monolayer.

For the last part of my thesis I designed a novel tension sensor to probe tension across N-cadherin (NTS). N-cadherin plays a critical role in cardiomyocytes, vascular smooth muscle cells, neurons and neural crest cells. Similar to E-cadherin, N-cadherin is also believed to bear tension and play a role in mechanotransduction pathways. To identify the role of tension across N-cadherin I designed a novel FRET-based molecular tension sensor for N-cadherin. I tested the ability of NTS to sense molecular tension in vascular smooth muscle cells, cardiomyocytes and cancer cells. Finally in collaboration with the Horwitz lab we have been able to show a role of tension across N-cadherin in synaptogenesis of neurons.