MHC–Peptide Binding is Assisted by Bound Water Molecules

Impact of water models on the structure and dynamics of enzyme tunnels

Protein hydration plays a vital role in many biological functions, and molecular dynamics simulations are frequently used to study it. However, the accuracy of these simulations is often sensitive to the water model used, a phenomenon particularly evident in intrinsically disordered proteins. Here, we investigated the extent to which the choice of water model alters the behavior of complex networks of tunnels within proteins. Tunnels are essential because they allow the exchange of substrates and products between buried enzyme active sites and the bulk solvent, directly affecting enzyme efficiency and selectivity, making the study of tunnels crucial for a holistic understanding of enzyme function at the molecular level. By performing simulations of haloalkane dehalogenase LinB and its two variants with engineered tunnels using TIP3P and OPC models, we investigated their effects on the overall tunnel topology. We also analyzed the properties of the primary tunnels, including their conformation, bottleneck dimensions, sampling efficiency, and the duration of tunnel openings. Our data demonstrate that all three proteins exhibited similar conformational behavior in both models but differed in the geometrical characteristics of their auxiliary tunnels, consistent with experimental observations. Interestingly, the results indicate that the stability of the open tunnels might be sensitive to the water model used. Because TIP3P can provide comparable data on the overall tunnel network, it is a valid choice when computational resources are limited or compatibility issues impede the use of OPC. However, OPC seems preferable for calculations requiring an accurate description of transport kinetics.

Allopurinol non-covalently facilitates binding of unconventional peptides to HLA-B*58:01

Allopurinol, widely used in gout treatment, is the most common cause of severe cutaneous adverse drug reactions. The risk of developing such life-threatening reactions is increased particularly for HLA-B*58:01 positive individuals. However the mechanism of action between allopurinol and HLA remains unknown. We demonstrate here that a Lamin A/C peptide KAGQVVTI which is unable to bind HLA-B*58:01 on its own, is enabled to form a stable peptide-HLA complex only in the presence of allopurinol. Crystal structure analysis reveal that allopurinol non-covalently facilitated KAGQVVTI to adopt an unusual binding conformation, whereby the C-terminal isoleucine does not engage as a PΩ that typically fit deeply in the binding F-pocket. A similar observation, though to a lesser degree was seen with oxypurinol. Presentation of unconventional peptides by HLA-B*58:01 aided by allopurinol contributes to our fundamental understanding of drug-HLA interactions. The binding of peptides from endogenously available proteins such as self-protein lamin A/C and viral protein EBNA3B suggest that aberrant loading of unconventional peptides in the presence of allopurinol or oxypurinol may be able to trigger anti-self reactions that can lead to Stevens-Johnson syndrome/toxic epidermal necrolysis (SJS/TEN) and Drug Reaction with Eosinophilia and Systemic Symptoms (DRESS).

Evolutionary Role of Water-Accessible Cavities in Src Homology 2 (SH2) Domains

Protein excited states are fundamental in the understanding of biological function, despite the fact they are hardly observed using traditional biophysical methodologies. Pressure perturbation coupled with nuclear magnetic resonance (NMR) spectroscopy is a powerful physicochemical tool to glance at these low-populated high-energy states on a residue-by-residue basis and underpin mechanistic insights into protein functionalities. Here we performed pressure titrations using NMR spectroscopy and relaxation dispersion experiments to identify the low-lying energetic states of the c-Abl SH2 domain. By showing that the SH2 excited state contains a hydrated hydrophobic cavity, fast-exchange motions, and highly conserved residues facing the water-accessible hole, we discuss the implications of water-protein interactions in SH2 modules achieving high-affinity binding and promiscuous phospho-Tyr peptide recognition.

The staphylococcal inhibitory protein SPIN binds to human myeloperoxidase with picomolar affinity but only dampens halide oxidation

The heme enzyme myeloperoxidase (MPO) is one of the key players in the neutrophil-mediated killing of invading pathogens as part of the innate immune system. MPO generates antimicrobial oxidants which indiscriminately and effectively kill phagocytosed pathogens. Staphylococcus aureus however is able to escape this fate, in part by secreting a small protein called SPIN (Staphylococcal Peroxidase Inhibitor), which specifically targets and inhibits MPO in a structurally complex manner. Here we present the first crystal structures of the complex of SPIN-aureus and a truncated version (SPIN-truncated) with mature dimeric leukocyte MPO. We unravel the contributions of the two domains to the kinetics and thermodynamics of SPIN-aureus binding to MPO by using a broad array of complementary biochemical and biophysical methods. The C-terminal "recognition" domain is shown to mediate specific binding to MPO, while interaction of the N-terminal "inhibitory" domain is guided mainly by hydrophobic effects and thus is less sequence-dependent. We found that inhibition of MPO is achieved by reducing substrate migration, but SPIN-aureus cannot completely block MPO activity. Its' effectiveness is inversely related to substrate size, with no discernible dependence on other factors. Thus, SPIN-aureus is an extremely high-affinity inhibitor and highly efficient for substrates larger than halogens. As aberrant MPO activity is implicated in a number of chronic inflammatory diseases, SPIN-aureus is the first promising protein inhibitor for specific inhibition of human MPO.

Homotypic and heterotypic in cis associations of MHC class I molecules at the cell surface

Through the presentation of peptide antigens to cytotoxic T lymphocytes, major histocompatibility complex (MHC) class I molecules mediate the adaptive immune response against tumors and viruses. Additional non-immunological functions include the heterotypic association of class I molecules with cell surface receptors, regulating their activities by unknown mechanisms. Also, homotypic associations resulting in class I dimers and oligomers - of unknown function - have been related to pathological outcomes. In this review, we provide an overview of the current knowledge about the occurrence, biochemical nature, and dynamics of homotypic and heterotypic associations of class I molecules at the cell surface with special focus on the molecular species that take part in the complexes and on the evidence that supports novel biological roles for class I molecules. We show that both heterotypic and homotypic class I associations reported in the literature describe not one but several kinds of oligomers with distinctive stoichiometry and biochemical properties.

•

Major histocompatibility complex class I molecules form homotypic and heterotypic associations at the cell surface.

•

Associations show distinctive stoichiometry and biochemical properties.

•

Associations might regulate immunological and non-immunological processes.

•

Heterotypic association with cell surface receptors might regulate receptor's activity.

•

Homotypic associations have been related to pathological outcomes.

Dynamics-Based Peptide-MHC Binding Optimization by a Convolutional Variational Autoencoder: A Use-Case Model for CASTELO

An unsolved challenge in the development of antigen-specific immunotherapies is determining the optimal antigens to target. Comprehension of antigen-major histocompatibility complex (MHC) binding is paramount toward achieving this goal. Here, we apply CASTELO, a combined machine learning-molecular dynamics (ML-MD) approach, to identify per-residue antigen binding contributions and then design novel antigens of increased MHC-II binding affinity for a type 1 diabetes-implicated system. We build upon a small-molecule lead optimization algorithm by training a convolutional variational autoencoder (CVAE) on MD trajectories of 48 different systems across four antigens and four HLA serotypes. We develop several new machine learning metrics including a structure-based anchor residue classification model as well as cluster comparison scores. ML-MD predictions agree well with experimental binding results and free energy perturbation-predicted binding affinities. Moreover, ML-MD metrics are independent of traditional MD stability metrics such as contact area and root-mean-square fluctuations (RMSF), which do not reflect binding affinity data. Our work supports the role of structure-based deep learning techniques in antigen-specific immunotherapy design.

Toward Guided Mutagenesis: Gaussian Process Regression Predicts MHC Class II Antigen Mutant Binding

Antigen-specific immunotherapies (ASI) require successful loading and presentation of antigen peptides into the major histocompatibility complex (MHC) binding cleft. One route of ASI design is to mutate native antigens for either stronger or weaker binding interaction to MHC. Exploring all possible mutations is costly both experimentally and computationally. To reduce experimental and computational expense, here we investigate the minimal amount of prior data required to accurately predict the relative binding affinity of point mutations for peptide-MHC class II (pMHCII) binding. Using data from different residue subsets, we interpolate pMHCII mutant binding affinities by Gaussian process (GP) regression of residue volume and hydrophobicity. We apply GP regression to an experimental data set from the Immune Epitope Database, and theoretical data sets from NetMHCIIpan and Free Energy Perturbation calculations. We find that GP regression can predict binding affinities of nine neutral residues from a six-residue subset with an average R² coefficient of determination value of 0.62 ± 0.04 (±95% CI), average error of 0.09 ± 0.01 kcal/mol (±95% CI), and with an receiver operating characteristic (ROC) AUC value of 0.92 for binary classification of enhanced or diminished binding affinity. Similarly, metrics increase to an R² value of 0.69 ± 0.04, average error of 0.07 ± 0.01 kcal/mol, and an ROC AUC value of 0.94 for predicting seven neutral residues from an eight-residue subset. Our work finds that prediction is most accurate for neutral residues at anchor residue sites without register shift. This work holds relevance to predicting pMHCII binding and accelerating ASI design.

Complement activation by drug carriers and particulate pharmaceuticals: Principles, challenges and opportunities

Considering the multifaceted protective and homeostatic roles of the complement system, many consequences arise when drug carriers, and particulate pharmaceutical formulations clash with complement proteins, and trigger complement cascade. Complement activation may induce formulation destabilization, promote opsonization, and affect biological and therapeutic performance of pharmaceutical nano- and micro-particles. In some cases, complement activation is beneficial, where complement may play a role in prophylactic protection, whereas uncontrolled complement activation is deleterious, and contributes to disease progression. Accordingly, design initiatives with particulate medicines should consider complement activation properties of the end formulation within the context of administration route, dosing, systems biology, and therapeutic perspective. Here we examine current progress in mechanistic processes underlying complement activation by pre-clinical and clinical particles, identify opportunities and challenges ahead, and suggest future directions in nanomedicine-complement interface research.

Exploring interfacial water trapping in protein-ligand complexes with multithermal titration calorimetry

In this work, we examine the hypothesis about how trapped water molecules at the interface between triosephosphate isomerase (TIM) and either of two phosphorylated inhibitors, 2-phosphoglycolate (2PG) or phosphoglycolohydroxamate (PGH), can explain the anomalous highly negative binding heat capacities (ΔC_p,b) of both complexes, TIM-2PG and TIM-PGH. We performed fluorimetric titrations of the enzyme with PGH inhibitor under osmotic stress conditions, using various concentrations of either osmolyte: sucrose, ethylene glycol or glycine betaine. We also analyze the binding processes under various stressor concentrations using a novel calorimetric methodology that allows ΔC_p,b determinations in single experiments: Multithermal Titration Calorimetry. The binding constant of the TIM-PGH complex decreased gradually with the concentration of all osmolytes, but at diverse extents depending on the osmolyte nature. According to the osmotic stress theory, this decrease indicates that the number of water molecules associated with the enzyme increases with inhibitor binding, i.e. some solvent molecules became trapped. Additionally, the binding heat capacities became less negative at higher osmolyte concentrations, their final values depending on the osmolyte. These effects were also observed in the TIM-2PG complex using sucrose as stressor. Our results strongly suggest that some water molecules became immobilized when the TIM-inhibitor complexes were formed. A computational analysis of the hydration state of the binding site of TIM in both its free state and its complexed form with 2PG or PGH, based on molecular dynamics (MD) simulations in explicit solvent, showed that the binding site effectively immobilized additional water molecules after binding these inhibitors.

Unraveling the structural basis for the unusually rich association of human leukocyte antigen DQ2.5 with class-II-associated invariant chain peptides

Human leukocyte antigen (HLA)-DQ2.5 (DQA1*05/DQB1*02) is a class-II major histocompatibility complex protein associated with both type 1 diabetes and celiac disease. One unusual feature of DQ2.5 is its high class-II-associated invariant chain peptide (CLIP) content. Moreover, HLA-DQ2.5 preferentially binds the non-canonical CLIP2 over the canonical CLIP1. To better understand the structural basis of HLA-DQ2.5's unusual CLIP association characteristics, better insight into the HLA-DQ2.5·CLIP complex structures is required. To this end, we determined the X-ray crystal structure of the HLA-DQ2.5· CLIP1 and HLA-DQ2.5·CLIP2 complexes at 2.73 and 2.20 Å, respectively. We found that HLA-DQ2.5 has an unusually large P4 pocket and a positively charged peptide-binding groove that together promote preferential binding of CLIP2 over CLIP1. An α9-α22-α24-α31-β86-β90 hydrogen bond network located at the bottom of the peptide-binding groove, spanning from the P1 to P4 pockets, renders the residues in this region relatively immobile. This hydrogen bond network, along with a deletion mutation at α53, may lead to HLA-DM insensitivity in HLA-DQ2.5. A molecular dynamics simulation experiment reported here and recent biochemical studies by others support this hypothesis. The diminished HLA-DM sensitivity is the likely reason for the CLIP-rich phenotype of HLA-DQ2.5.

In silico-based vaccine design against Ebola virus glycoprotein

Ebola virus (EBOV) is one of the lethal viruses, causing more than 24 epidemic outbreaks to date. Despite having available molecular knowledge of this virus, no definite vaccine or other remedial agents have been developed yet for the management and avoidance of EBOV infections in humans. Disclosing this, the present study described an epitope-based peptide vaccine against EBOV, using a combination of B-cell and T-cell epitope predictions, followed by molecular docking and molecular dynamics simulation approach. Here, protein sequences of all glycoproteins of EBOV were collected and examined via in silico methods to determine the most immunogenic protein. From the identified antigenic protein, the peptide region ranging from 186 to 220 and the sequence HKEGAFFLY from the positions of 154-162 were considered the most potential B-cell and T-cell epitopes, correspondingly. Moreover, this peptide (HKEGAFFLY) interacted with HLA-A*32:15 with the highest binding energy and stability, and also a good conservancy of 83.85% with maximum population coverage. The results imply that the designed epitopes could manifest vigorous enduring defensive immunity against EBOV.

Water Determines the Structure and Dynamics of Proteins

Water is an essential participant in the stability, structure, dynamics, and function of proteins and other biomolecules. Thermodynamically, changes in the aqueous environment affect the stability of biomolecules. Structurally, water participates chemically in the catalytic function of proteins and nucleic acids and physically in the collapse of the protein chain during folding through hydrophobic collapse and mediates binding through the hydrogen bond in complex formation. Water is a partner that slaves the dynamics of proteins, and water interaction with proteins affect their dynamics. Here we provide a review of the experimental and computational advances over the past decade in understanding the role of water in the dynamics, structure, and function of proteins. We focus on the combination of X-ray and neutron crystallography, NMR, terahertz spectroscopy, mass spectroscopy, thermodynamics, and computer simulations to reveal how water assist proteins in their function. The recent advances in computer simulations and the enhanced sensitivity of experimental tools promise major advances in the understanding of protein dynamics, and water surely will be a protagonist.

Stay Wet, Stay Stable? How Internal Water Helps the Stability of Thermophilic Proteins

We present a systematic computational investigation of the internal hydration of a set of homologous proteins of different stability content and molecular complexities. The goal of the study is to verify whether structural water can be part of the molecular mechanisms ensuring enhanced stability in thermophilic enzymes. Our free-energy calculations show that internal hydration in the thermophilic variants is generally more favorable, and that the cumulated effect of wetting multiple sites results in a meaningful contribution to stability. Moreover, thanks to a more effective capability to retain internal water, some thermophilic proteins benefit by a systematic gain from internal wetting up to their optimal working temperature. Our work supports the idea that internal wetting can be viewed as an alternative molecular variable to be tuned for increasing protein stability.

Dynamical insights into the differential characteristics of Mycobacterium avium subsp. paratuberculosis peptide binding to HLA-DRB1 proteins associated with multiple sclerosis

Differential properties of MAP binding to HLA proteins in Sardinian MS patients.

Cancer nanomedicine and the complement system activation paradigm: anaphylaxis and tumour growth

A wide variety of nanocarriers and particularly cancer nanomedicines activate the complement system, which is the first line of the innate immune defence mechanism. Complement activation may induce inflammatory responses, but such responses arising from uncontrolled complement activation could be life threatening. Accordingly, the role of complement in initiation of adverse reactions to particulate and polymer therapeutics is receiving increasing attention. Furthermore, the involvement of complement-activation products in promoting tumour growth has also been indicated. This could be of serious concern for development of cancer nanomedicines and cancer nanotechnology initiatives. These concepts are reviewed with preliminary evidence that intra-tumoural accumulation of model long circulating nanoparticles could promote tumour growth.

Altered backbone and side-chain interactions result in route heterogeneity during the folding of interleukin-1β (IL-1β)

Deletion of the β-bulge trigger-loop results in both a switch in the preferred folding route, from the functional loop packing folding route to barrel closure, as well as conversion of the agonist activity of IL-1β into antagonist activity. Conversely, circular permutations of IL-1β conserve the functional folding route as well as the agonist activity. These two extremes in the folding-functional interplay beg the question of whether mutations in IL-1β would result in changes in the populations of heterogeneous folding routes and the signaling activity. A series of topologically equivalent water-mediated β-strand bridging interactions within the pseudosymmetric β-trefoil fold of IL-1β highlight the backbone water interactions that stabilize the secondary and tertiary structure of the protein. Additionally, conserved aromatic residues lining the central cavity appear to be essential for both stability and folding. Here, we probe these protein backbone-water molecule and side chain-side chain interactions and the role they play in the folding mechanism of this geometrically stressed molecule. We used folding simulations with structure-based models, as well as a series of folding kinetic experiments to examine the effects of the F42W core mutation on the folding landscape of IL-1β. This mutation alters water-mediated backbone interactions essential for maintaining the trefoil fold. Our results clearly indicate that this perturbation in the primary structure alters a structural water interaction and consequently modulates the population of folding routes accessed during folding and signaling activity.

Quantitative prediction of peptide binding to HLA-DP1 protein

The exogenous proteins are processed by the host antigen-processing cells. Peptidic fragments of them are presented on the cell surface bound to the major hystocompatibility complex (MHC) molecules class II and recognized by the CD4+ T lymphocytes. The MHC binding is considered as the crucial prerequisite for T-cell recognition. Only peptides able to form stable complexes with the MHC proteins are recognized by the T-cells. These peptides are known as T-cell epitopes. All T-cell epitopes are MHC binders, but not all MHC binders are T-cell epitopes. The T-cell epitope prediction is one of the main priorities of immunoinformatics. In the present study, three chemometric techniques are combined to derive a model for in silico prediction of peptide binding to the human MHC class II protein HLA-DP1. The structures of a set of known peptide binders are described by amino acid z-descriptors. Data are processed by an iterative self-consisted algorithm using the method of partial least squares, and a quantitative matrix (QM) for peptide binding prediction to HLA-DP1 is derived. The QM is validated by two sets of proteins and showed an average accuracy of 86 percent.

Using databases and data mining in vaccinology

Throughout time functional immunology has accumulated vast amounts of quantitative and qualitative data relevant to the design and discovery of vaccines. Such data includes, but is not limited to, components of the host and pathogen genome (including antigens and virulence factors), T- and B-cell epitopes and other components of the antigen presentation pathway and allergens. In this review the authors discuss a range of databases that archive such data. Built on such information, increasingly sophisticated data mining techniques have developed that create predictive models of utilitarian value. With special reference to epitope data, the authors discuss the strengths and weaknesses of the available techniques and how they can aid computer-aided vaccine design deliver added value for vaccinology.

T-cell epitope vaccine design by immunoinformatics

Vaccination is generally considered to be the most effective method of preventing infectious diseases. All vaccinations work by presenting a foreign antigen to the immune system in order to evoke an immune response. The active agent of a vaccine may be intact but inactivated ('attenuated') forms of the causative pathogens (bacteria or viruses), or purified components of the pathogen that have been found to be highly immunogenic. The increased understanding of antigen recognition at molecular level has resulted in the development of rationally designed peptide vaccines. The concept of peptide vaccines is based on identification and chemical synthesis of B-cell and T-cell epitopes which are immunodominant and can induce specific immune responses. The accelerating growth of bioinformatics techniques and applications along with the substantial amount of experimental data has given rise to a new field, called immunoinformatics. Immunoinformatics is a branch of bioinformatics dealing with in silico analysis and modelling of immunological data and problems. Different sequence- and structure-based immunoinformatics methods are reviewed in the paper.

Biomolecular electrostatics and solvation: a computational perspective

An understanding of molecular interactions is essential for insight into biological systems at the molecular scale. Among the various components of molecular interactions, electrostatics are of special importance because of their long-range nature and their influence on polar or charged molecules, including water, aqueous ions, proteins, nucleic acids, carbohydrates, and membrane lipids. In particular, robust models of electrostatic interactions are essential for understanding the solvation properties of biomolecules and the effects of solvation upon biomolecular folding, binding, enzyme catalysis, and dynamics. Electrostatics, therefore, are of central importance to understanding biomolecular structure and modeling interactions within and among biological molecules. This review discusses the solvation of biomolecules with a computational biophysics view toward describing the phenomenon. While our main focus lies on the computational aspect of the models, we provide an overview of the basic elements of biomolecular solvation (e.g. solvent structure, polarization, ion binding, and non-polar behavior) in order to provide a background to understand the different types of solvation models.

Explicit treatment of water molecules in data-driven protein-protein docking: the solvated HADDOCKing approach

Water molecules are active components in, literally, every biochemical event, forming hydrogen bonds, filling cavities, and mediating interactions with other (bio)molecules. Therefore, solvent drastically affects the kinetics and thermodynamics of numerous cellular events, including protein-protein interactions. While docking techniques are becoming successful in predicting the three-dimensional structure of protein-protein complexes, they are still limited in accounting explicitly for water in the binding process. HADDOCK is one of the few programs so far that can explicitly deal with water molecules during docking. Its solvated docking protocol starts from hydrated molecules, and a fraction of the interfacial water is subsequently removed from the docked models in a biased Monte Carlo procedure. The Monte Carlo-based removal is based on interfacial amino acid-water contact propensities derived from a dataset of high-resolution crystal structures of protein-protein complexes. In this chapter, this solvated docking protocol is described and associated methodological aspects are illustrated through an application example. It is shown that, although docking results do not always improve when the solvated docking protocol is applied, scoring is improved and the positions of buried water molecules in an interface are correctly predicted. Therefore, by identifying important water molecules, solvated docking can aid the development of novel inhibitors of protein-protein complexes that might be better accommodated at an interface.

MOLECULAR DYNAMICS SIMULATION OF POINT MUTATIONS OF ARG 221 IN THE ACTIVE SITE OF PROTEIN TYROSINE PHOSPHATASE 1B

Protein tyrosine phosphatase 1B (PTP1B) is one of the important regulators of signal transduction pathways. The present study aims to investigate the effect of Arg 221 on the active site of PTP1B. Six mutants were carried out using Schrödinger Suite 2007 and molecular dynamics simulation was performed by using the Tinker package. Results show that point mutations at position 221 have great influence on shape of active site, backbone movement of active site, and interaction between substrate and PTP1B. R221H and R221K lead to increased total interaction energies. R221G, R221F and R221T cause increase in total interaction energies, but decrease in interaction energies between pTyr 4 and P loop (catalytic residues). R221E results in both decreased total interaction energies and interaction energies between pTyr 4 and P loop. This indicates that Arg 221 mutated to basic residues can lead to enhanced binding affinity between substrate and protein; when mutated to acidic residues it will decrease binding affinity and catalytic activity; other kinds of mutations result in increased binding affinity but decreased catalytic activity.

Loss of T cell antigen recognition arising from changes in peptide and major histocompatibility complex protein flexibility: implications for vaccine design

Modification of the primary anchor positions of antigenic peptides to improve binding to major histocompatibility complex (MHC) proteins is a commonly used strategy for engineering peptide-based vaccine candidates. However, such peptide modifications do not always improve antigenicity, complicating efforts to design effective vaccines for cancer and infectious disease. Here we investigated the MART-1(27-35) tumor antigen, for which anchor modification (replacement of the position two alanine with leucine) dramatically reduces or ablates antigenicity with a wide range of T cell clones despite significantly improving peptide binding to MHC. We found that anchor modification in the MART-1(27-35) antigen enhances the flexibility of both the peptide and the HLA-A*0201 molecule. Although the resulting entropic effects contribute to the improved binding of the peptide to MHC, they also negatively impact T cell receptor binding to the peptide·MHC complex. These results help explain how the "anchor-fixing" strategy fails to improve antigenicity in this case, and more generally, may be relevant for understanding the high specificity characteristic of the T cell repertoire. In addition to impacting vaccine design, modulation of peptide and MHC flexibility through changes to antigenic peptides may present an evolutionary strategy for the escape of pathogens from immune destruction.

Structural insights into the binding of hepatitis B virus core peptide to HLA-A2 alleles: towards designing better vaccines

Binding of specific antigenic peptides with human leukocyte antigen (HLA) molecules is a prerequisite for the initiation of T-cell responses and structural information about the peptide-HLA complex is essential for the detailed understanding of such interactions. HLA-A2 is the most prevalent HLA allele globally but aside from A*02:01 there is a significant lack of crystal structures, particularly for alleles that occur in high frequencies among Asian populations. Here, we report three HLA-A2 structures with the immunodominant hepatitis B core antigen 18-27 (HBcAg18-27) epitope, namely A*02:03, A*02:06, and A*02:07 at resolutions of 2.16, 1.70, and 1.75 Å respectively. This comparative analysis reveals that minor polymorphic residue changes between different HLA alleles can induce significant alterations in the major histocompatibility complex-peptide interface, and introduce conformational changes in the p3-p8 peptide region. Circular dichroism analysis demonstrated the HLA-A2-peptide complexes to have a hierarchy of thermostability and binding affinity in the order of A*02:06>A*02:07>A*02:01>A*02:03. Our findings provide structural insights into the varied HLA-A2 allele binding of the hepatitis B core antigen 18-27 epitope and the data suggest that chemical modifications of the peptide side chains could be a promising strategy to modulate and improve HLA-A2-peptide binding affinity for vaccine design.

Simulations of the confinement of ubiquitin in self-assembled reverse micelles

We describe the effects of confinement on the structure, hydration, and the internal dynamics of ubiquitin encapsulated in reverse micelles (RM). We performed molecular dynamics simulations of the encapsulation of ubiquitin into self-assembled protein/surfactant reverse micelles to study the positioning and interactions of the protein with the RM and found that ubiquitin binds to the RM interface at low salt concentrations. The same hydrophobic patch that is recognized by ubiquitin binding domains in vivo is found to make direct contact with the surfactant head groups, hydrophobic tails, and the iso-octane solvent. The fast backbone N-H relaxation dynamics show that the fluctuations of the protein encapsulated in the RM are reduced when compared to the protein in bulk. This reduction in fluctuations can be explained by the direct interactions of ubiquitin with the surfactant and by the reduced hydration environment within the RM. At high concentrations of excess salt, the protein does not bind strongly to the RM interface and the fast backbone dynamics are similar to that of the protein in bulk. Our simulations demonstrate that the confinement of protein can result in altered protein dynamics due to the interactions between the protein and the surfactant.

Water in the polar and nonpolar cavities of the protein interleukin-1β

Water in the protein interior serves important structural and functional roles and is also increasingly recognized as a relevant factor in drug binding. The nonpolar cavity in the protein interleukin-1β has been reported to be filled by water on the basis of some experiments and simulations and to be empty on the basis of others. Here we study the thermodynamics of filling the central nonpolar cavity and the four polar cavities of interleukin-1β by molecular dynamics simulation. We use different water models (TIP3P and SPC/E) and protein force fields (amber94 and amber03) to calculate the semigrand partition functions term by term that quantify the hydration equilibria. We consistently find that water in the central nonpolar cavity is thermodynamically unstable, independent of force field and water model. The apparent reason is the relatively small size of the cavity, with a volume less than ∼80 Å(3). Our results are consistent with the most recent X-ray crystallographic and simulation studies but disagree with an earlier interpretation of nuclear magnetic resonance (NMR) experiments probing protein-water interactions. We show that, at least semiquantitatively, the measured nuclear Overhauser effects indicating the proximity of water to the methyl groups lining the nonpolar cavity can, in all likelihood, be attributed to interactions with buried and surface water molecules near the cavity. The same methods applied to determine the occupancy of the polar cavities show that they are filled by the same number of water molecules observed in crystallography, thereby validating the theoretical and simulation methods used to study the water occupancy in the nonpolar protein cavity.

vmdICE: a plug-in for rapid evaluation of molecular dynamics simulations using VMD

Molecular dynamics (MD) is a powerful in silico method to investigate the interactions between biomolecules. It solves Newton's equations of motion for atoms over a specified period of time and yields a trajectory file, containing the different spatial arrangements of atoms during the simulation. The movements and energies of each single atom are recorded. For evaluating of these simulation trajectories with regard to biomedical implications, several methods are available. Three well-known ones are the root mean square deviation (RMSD), the root mean square fluctuation (RMSF) and solvent accessible surface area (SASA). Herein, we present a novel plug-in for the software "visual molecular dynamics" (VMD) that allows an interactive 3D representation of RMSD, RMSF, and SASA, directly on the molecule. On the one hand, our plug-in is easy to handle for inexperienced users, and on the other hand, it provides a fast and flexible graphical impression of the spatial dynamics of a system for experts in the field.

Recognition and binding of a helix-loop-helix peptide to carbonic anhydrase occurs via partly folded intermediate structures

We have studied the association of a helix-loop-helix peptide scaffold carrying a benzenesulfonamide ligand to carbonic anhydrase using steady-state and time-resolved fluorescence spectroscopy. The helix-loop-helix peptide, developed for biosensing applications, is labeled with the fluorescent probe dansyl, which serves as a polarity-sensitive reporter of the binding event. Using maximum entropy analysis of the fluorescence lifetime of dansyl at 1:1 stoichiometry reveals three characteristic fluorescence lifetime groups, interpreted as differently interacting peptide/protein structures. We characterize these peptide/protein complexes as mostly bound but unfolded, bound and partly folded, and strongly bound and folded. Furthermore, analysis of the fluorescence anisotropy decay resulted in three different dansyl rotational correlation times, namely 0.18, 1.2, and 23 ns. Using the amplitudes of these times, we can correlate the lifetime groups with the corresponding fluorescence anisotropy component. The 23-ns rotational correlation time, which appears with the same amplitude as a 17-ns fluorescence lifetime, shows that the dansyl fluorophore follows the rotational diffusion of carbonic anhydrase when it is a part of the folded peptide/protein complex. A partly folded and partly hydrated interfacial structure is manifested in an 8-ns dansyl fluorescence lifetime and a 1.2-ns rotational correlation time. This structure, we believe, is similar to a molten-globule-like interfacial structure, which allows segmental movement and has a higher degree of solvent exposure of dansyl. Indirect excitation of dansyl on the helix-loop-helix peptide through Förster energy transfer from one or several tryptophans in the carbonic anhydrase shows that the helix-loop-helix scaffold binds to a tryptophan-rich domain of the carbonic anhydrase. We conclude that binding of the peptide to carbonic anhydrase involves a transition from a disordered to an ordered structure of the helix-loop-helix scaffold.

A molecular dynamics approach to study the importance of solvent in protein interactions

Water constitutes the cellular environment for biomolecules to interact. Solvent is important for protein folding and stability, and it is also known to actively participate in many catalytic processes in the cell. However, solvent is often ignored in molecular recognition and not taken into account in protein-protein interaction studies and rational design. Previously we developed SCOWLP, a database and its web application (http://www.scowlp.org), to perform studies on the contribution of solvent to protein interface definition in all protein complexes of the PDB. We introduced the concept of wet spots, interfacial residues interacting only through one water molecule, which were shown to considerably enrich protein interface descriptions. Analysis of interfacial solvent in a nonredundant dataset of protein complexes suggested the importance of including interfacial water molecules in protein interaction studies. In this work we use a molecular dynamics approach to gain deeper insights into solvent contribution to protein interfaces. We characterize the dynamic and energetic properties of water-mediated protein interactions by comparing different interfacial interaction types (direct, dual and wet spot) at residue and solvent level. For this purpose, we perform an analysis of 17 representative complexes from two protein families of different interface nature. Energetically wet spots are quantitatively comparable to other residues in interfaces, and their mobility is shown to be lower than protein surface residues. The residence time of water molecules in wet spots sites is higher than of those on the surface of the protein. In terms of free energy, though wet-spots-forming water molecules are very heterogeneous, their contribution to the free energy of complex formation is considerable. We find that water molecules can play an important role in interaction conservation in protein interfaces by allowing sequence variability in the corresponding binding partner, and we discuss the important implications of our observations related to the use of the correlated mutations concept in protein interactions studies. The results obtained in this work help to deepen our understanding of the physico-chemical nature underlying protein-protein interactions and strengthen the idea of using the wet spots concept to qualitatively improve the accuracy of folding, docking and rational design algorithms.

Non-canonical anchor motif peptides bound to MHC class I induce cellular responses

The major histocompatibility complex (MHC) on the surface of antigen presenting cells functions to display peptides to the T cell receptor (TCR). Recognition of peptide-MHC by T cells initiates a cascade of signals, which results in the initiation of a T cell dependent immune response. An understanding of how peptides bind to MHC molecules is important for determining the structural basis for T cell dependent immune responses and facilitates the structure-based design of peptides as candidate vaccines to elicit a specific immune response. To date, crystal structures, immunogenicity and in vivo biological relevance have mainly been characterized for high affinity peptide-MHC interactions. From the crystal structures of numerous peptide-MHC complexes it became apparent what canonical sequence features were required for high affinity binding, which led to the ability to predict in most instances peptides with high affinity for MHC. We previously identified the crystal structures of non-canonical peptides in complex with MHC class I (one bound with low affinity and the other with high affinity, but utilizing novel peptide anchors and MHC pockets). It is becoming increasingly evident that other non-canonical peptides can also bind, such as long-, short- and glyco-peptides. However, the in vivo role of non-canonical peptides is not clear and we present here the immunogenicity of two non-canonical peptides and their affinity when bound to MHC class I, H2K(b). Comparison of the three-dimensional structures in complex with MHC suggests major differences in hydrogen bonding patterns with H2K(b), despite sharing similar binding modes, which may account for the differences in affinity and immunogenicity. These studies provide further evidence for the diverse range of peptide ligands that can bind to MHC and be recognized by the TCR, which will facilitate approaches to peptide-based vaccine design.

Protein-protein interaction investigated by steered molecular dynamics: the TCR-pMHC complex

We present a novel steered molecular dynamics scheme to induce the dissociation of large protein-protein complexes. We apply this scheme to study the interaction of a T cell receptor (TCR) with a major histocompatibility complex (MHC) presenting a peptide (p). Two TCR-pMHC complexes are considered, which only differ by the mutation of a single amino acid on the peptide; one is a strong agonist that produces T cell activation in vivo, while the other is an antagonist. We investigate the interaction mechanism from a large number of unbinding trajectories by analyzing van der Waals and electrostatic interactions and by computing energy changes in proteins and solvent. In addition, dissociation potentials of mean force are calculated with the Jarzynski identity, using an averaging method developed for our steering scheme. We analyze the convergence of the Jarzynski exponential average, which is hampered by the large amount of dissipative work involved and the complexity of the system. The resulting dissociation free energies largely underestimate experimental values, but the simulations are able to clearly differentiate between wild-type and mutated TCR-pMHC and give insights into the dissociation mechanism.

Crystallographic study of hydration of an internal cavity in engineered proteins with buried polar or ionizable groups

Although internal water molecules are essential for the structure and function of many proteins, the structural and physical factors that govern internal hydration are poorly understood. We have examined the molecular determinants of internal hydration systematically, by solving the crystal structures of variants of staphylococcal nuclease with Gln-66, Asn-66, and Tyr-66 at cryo (100 K) and room (298 K) temperatures, and comparing them with existing cryo and room temperature structures of variants with Glu-66, Asp-66, Lys-66, Glu-92 or Lys-92 obtained under conditions of pH where the internal ionizable groups are in the neutral state. At cryogenic temperatures the polar moieties of all these internal side chains are hydrated except in the cases of Lys-66 and Lys-92. At room temperature the internal water molecules were observed only in variants with Glu-66 and Tyr-66; water molecules in the other variants are probably present but they are disordered and therefore undetectable crystallographically. Each internal water molecule establishes between 3 and 5 hydrogen bonds with the protein or with other internal water molecules. The strength of interactions between internal polar side chains and water molecules seems to decrease from carboxylic acids to amides to amines. Low temperature, low cavity volume, and the presence of oxygen atoms in the cavity increase the positional stability of internal water molecules. This set of structures and the physical insight they contribute into internal hydration will be useful for the development and benchmarking of computational methods for artificial hydration of pockets, cavities, and active sites in proteins.

A balancing act between net uptake of water during dihydrofolate binding and net release of water upon NADPH binding in R67 dihydrofolate reductase

R67 dihydrofolate reductase (DHFR) catalyzes the reduction of dihydrofolate (DHF) to tetrahydrofolate using NADPH as a cofactor. This enzyme is a homotetramer possessing 222 symmetry, and a single active site pore traverses the length of the protein. A promiscuous binding surface can accommodate either DHF or NADPH, thus two nonproductive complexes can form (2NADPH or 2DHF) as well as a productive complex (NADPH.DHF). The role of water in binding was monitored using a number of different osmolytes. From isothermal titration calorimetry (ITC) studies, binding of NADPH is accompanied by the net release of 38 water molecules. In contrast, from both steady state kinetics and ITC studies, binding of DHF is accompanied by the net uptake of water. Although different osmolytes have similar effects on NADPH binding, variable results are observed when DHF binding is probed. Sensitivity to water activity can also be probed by an in vivo selection using the antibacterial drug, trimethoprim, where the water content of the media is decreased by increasing concentrations of sorbitol. The ability of wild type and mutant clones of R67 DHFR to allow host Escherichia coli to grow in the presence of trimethoprim plus added sorbitol parallels the catalytic efficiency of the DHFR clones, indicating water content strongly correlates with the in vivo function of R67 DHFR.

Water penetration in the low and high pressure native states of ubiquitin

Theoretical studies on the solvation of methane molecules in water have shown that the effect of increased pressure is to stabilize solvent separated contacts relative to direct contacts. This suggests that high pressure stabilizes waters that have penetrated into a protein's core, indicating a mechanism for the high pressure denaturation of proteins. We test this theory on a folded protein by studying the penetration of water into the native state of ubiquitin at low and high pressures, using molecular dynamics. An ensemble of conformations sampled in the folded state of ubiquitin has been determined by NMR at two pressures below the protein's denaturation pressure, 30 atm and 3000 atm. We find that 1-5 more waters penetrate the high pressure conformations than the low pressure conformations. Low volume configurations of the system are favored at high pressures, but different components of the system may experience increases or decreases in their specific volumes. We find that penetrating waters have a higher volume per water than bulk waters, but that the volume per protein residue may be lowered by solvation. Furthermore, we find that penetration of the protein by water at high pressures is driven by the difference in the pressure dependence of the probability of cavity opening in the protein and pressure dependence of the probability of cavity opening in the bulk solvent. The volume changes associated with cavity opening and closing indicate that each penetrating water reduces the volume of the system by about 12 mL/mol. The experimental volume change going from the low pressure to the high pressure native state of ubiquitin is 24 mL/mol. Our results indicate that this volume change can be explained by penetration of the protein by two water molecules.

Role of flexibility and polarity as determinants of the hydration of internal cavities and pockets in proteins

Molecular dynamics simulations of Staphylococcal nuclease and of 10 variants with internal polar or ionizable groups were performed to investigate systematically the molecular determinants of hydration of internal cavities and pockets in proteins. In contrast to apolar cavities in rigid carbon structures, such as nanotubes or buckeyballs, internal cavities in proteins that are large enough to house a few water molecules will most likely be dehydrated unless they contain a source of polarity. The water content in the protein interior can be modulated by the flexibility of protein elements that interact with water, which can impart positional disorder to water molecules, or bias the pattern of internal hydration that is stabilized. This might explain differences in the patterns of hydration observed in crystal structures obtained at cryogenic and room temperature conditions. The ability of molecular dynamics simulations to determine the most likely sites of water binding in internal pockets and cavities depends on its efficiency in sampling the hydration of internal sites and alternative protein and water conformations. This can be enhanced significantly by performing multiple molecular dynamics simulations as well as simulations started from different initial hydration states.

Comparative molecular dynamics analysis of tapasin-dependent and -independent MHC class I alleles

MHC class I molecules load antigenic peptides in the endoplasmic reticulum and present them at the cell surface. Efficiency of peptide loading depends on the class I allele and can involve interaction with tapasin and other proteins of the loading complex. Allele HLA-B*4402 (Asp at position 116) depends on tapasin for efficient peptide loading, whereas HLA-B*4405 (identical to B*4402 except for Tyr116) can efficiently load peptides in the absence of tapasin. Both alleles adopt very similar structures in the presence of the same peptide. Comparative unrestrained molecular dynamics simulations on the alpha(1)/alpha(2) peptide binding domains performed in the presence of bound peptides resulted in structures in close agreement with experiments for both alleles. In the absence of peptides, allele-specific conformational changes occurred in the first segment of the alpha(2)-helix that flanks the peptide C-terminal binding region (F-pocket) and contacts residue 116. This segment is also close to the proposed tapasin contact region. For B*4402, a shift toward an altered F-pocket structure deviating significantly from the bound form was observed. Subsequent free energy simulations on induced F-pocket opening in B*4402 confirmed a conformation that deviated significantly from the bound structure. For B*4405, a free energy minimum close to the bound structure was found. The simulations suggest that B*4405 has a greater tendency to adopt a peptide receptive conformation in the absence of peptide, allowing tapasin-independent peptide loading. A possible role of tapasin could be the stabilization of a peptide-receptive class I conformation for HLA-B*4402 and other tapasin-dependent alleles.