HLA-DP2 binding prediction by molecular dynamics simulations

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

Peptide binding to HLA-DP proteins at pH 5.0 and pH 7.0:a quantitative molecular docking study

Background: HLA-DPs are class II MHC proteins mediating immune responses to many diseases. Peptides bind MHC class II proteins in the acidic environment within endosomes. Acidic pH markedly elevates association rate constants but dissociation rates are almost unchanged in the pH range 5.0 - 7.0. This pH-driven effect can be explained by the protonation/deprotonation states of Histidine, whose imidazole has a pKa of 6.0. At pH 5.0, imidazole ring is protonated, making Histidine positively charged and very hydrophilic, while at pH 7.0 imidazole is unprotonated, making Histidine less hydrophilic. We develop here a method to predict peptide binding to the four most frequent HLA-DP proteins: DP1, DP41, DP42 and DP5, using a molecular docking protocol. Dockings to virtual combinatorial peptide libraries were performed at pH 5.0 and pH 7.0. Results: The X-ray structure of the peptide - HLA-DP2 protein complex was used as a starting template to model by homology the structure of the four DP proteins. The resulting models were used to produce virtual combinatorial peptide libraries constructed using the single amino acid substitution (SAAS) principle. Peptides were docked into the DP binding site using AutoDock at pH 5.0 and pH 7.0. The resulting scores were normalized and used to generate Docking Score-based Quantitative Matrices (DS-QMs). The predictive ability of these QMs was tested using an external test set of 484 known DP binders. They were also compared to existing servers for DP binding prediction. The models derived at pH 5.0 predict better than those derived at pH 7.0 and showed significantly improved predictions for three of the four DP proteins, when compared to the existing servers. They are able to recognize 50% of the known binders in the top 5% of predicted peptides. Conclusions: The higher predictive ability of DS-QMs derived at pH 5.0 may be rationalised by the additional hydrogen bond formed between the backbone carbonyl oxygen belonging to the peptide position before p1 (p-1) and the protonated ε-nitrogen of His 79β. Additionally, protonated His residues are well accepted at most of the peptide binding core positions which is in a good agreement with the overall negatively charged peptide binding site of most MHC proteins. © 2012 Patronov et al.; licensee BioMed Central Ltd.

Tailoring heme -thiolate proteins into efficient biocatalysts with high specificity and selectivity

Cytochrome P450 monooxygenases, one of the most important classes of heme-thiolate proteins, have attracted considerable interest in the biochemical community because of its catalytic versatility, substrate diversity and great number in the superfamily. Although P450s are capable of catalyzing numerous difficult oxidation reactions, the relatively low stability, low turnover rates and the need of electron-donating cofactors have limited their practical biotechnological and pharmaceutical applications as isolated enzymes. The goal of this study is to tailor such heme-thiolate proteins into efficient biocatalysts with high specificity and selectivity by protein engineering and to better understand the structure-function relationship in cytochromes P450. In the effort to engineer P450cam, the prototype member of the P450 superfamily, into an efficient peroxygenase that utilizes hydrogen peroxide via the “peroxide-shunt” pathway, site-directed mutagenesis has been used to elucidate the critical roles of hydrophobic residues in the active site. Various biophysical, biochemical and spectroscopic techniques have been utilized to investigate the wild-type and mutant proteins. Three important P450cam variants were obtained showing distinct structural and functional features. In P450camV247H mutant, which exhibited almost identical spectral properties with the wild-type, it is demonstrated that a single amino acid switch turned the monooxygenase into an efficient preoxidase by increasing the peroxidase activity nearly one thousand folds. In order to tune the distal pocket of P450cam with polar residues, Leu 246 was replaced with a basic residue, lysine, resulting in a mutant with spectral features identical to P420, the inactive species of P450. But this inactive-species-like mutant showed catalytic activities without the facilitation of any cofactors. By substituting Gly 248 with a histidine, a novel Cys-Fe-His ligation set was obtained in P450cam which represented the very rare case of His ligation in heme-thiolate proteins. In addition to serving as a convenient model for hemoprotein structural studies, the G248H mutant also provided evidence about the nature of the axial ligand in cytochrome P420 and other engineered hemoproteins with thiolate ligations. Furthermore, attempts have been made to replace the proximal ligand in sperm whale myoglobin to construct a heme-thiolate protein model by mimicking the protein environment of cytochrome P450cam and chloroperoxidase.

Tailoring Heme-Thiolate Proteins into Efficient Biocatalysts with High Specificity and Selectivity

Cytochrome P450 monooxygenases, one of the most important classes of heme-thiolate proteins, have attracted considerable interest in the biochemical community because of its catalytic versatility, substrate diversity and great number in the superfamily. Although P450s are capable of catalyzing numerous difficult oxidation reactions, the relatively low stability, low turnover rates and the need of electron-donating cofactors have limited their practical biotechnological and pharmaceutical applications as isolated enzymes. The goal of this study is to tailor such heme-thiolate proteins into efficient biocatalysts with high specificity and selectivity by protein engineering and to better understand the structure-function relationship in cytochromes P450. In the effort to engineer P450cam, the prototype member of the P450 superfamily, into an efficient peroxygenase that utilizes hydrogen peroxide via the “peroxide-shunt” pathway, site-directed mutagenesis has been used to elucidate the critical roles of hydrophobic residues in the active site. Various biophysical, biochemical and spectroscopic techniques have been utilized to investigate the wild-type and mutant proteins. Three important P450cam variants were obtained showing distinct structural and functional features. In P450camV247H mutant, which exhibited almost identical spectral properties with the wild-type, it is demonstrated that a single amino acid switch turned the monooxygenase into an efficient preoxidase by increasing the peroxidase activity nearly one thousand folds. In order to tune the distal pocket of P450cam with polar residues, Leu 246 was replaced with a basic residue, lysine, resulting in a mutant with spectral features identical to P420, the inactive species of P450. But this inactive-species-like mutant showed catalytic activities without the facilitation of any cofactors. By substituting Gly 248 with a histidine, a novel Cys-Fe-His ligation set was obtained in P450cam which represented the very rare case of His ligation in heme-thiolate proteins. In addition to serving as a convenient model for hemoprotein structural studies, the G248H mutant also provided evidence about the nature of the axial ligand in cytochrome P420 and other engineered hemoproteins with thiolate ligations. Furthermore, attempts have been made to replace the proximal ligand in sperm whale myoglobin to construct a heme-thiolate protein model by mimicking the protein environment of cytochrome P450cam and chloroperoxidase.

7-Amino actinomycin D bound to single-stranded hairpin DNA enhanced by loop sequence-dependent luminescent Eu3+ and Tb3+ binding

Lanthanide Eu3+ and Tb3+ ions have been widely used in luminescent resonance energy transfer (LRET) for bioassays to study metal binding microenvironments. We report here that Eu3+ or Tb3+ can increase the binding affinity of antitumor antibiotic drug agent, 7-amino actinomycin D (7AACTD), binding to 5'-GT/TG-5' or 5'-GA/AG-5' mismatched stem region of the single-stranded hairpin DNA. Further studies indicate that the effect of Eu3+ or Tb3+ on 7AACTD binding is related to DNA loop sequence. Our results will provide new insights into how metal ions can enhance antitumor agents binding to their targets.

A forensic investigation of single human hair fibres using FTIR-ATR spectroscopy and chemometrics

Human hair fibres are ubiquitous in nature and are found frequently at crime scenes often as a result of exchange between the perpetrator, victim and/or the surroundings according to Locard's Principle. Therefore, hair fibre evidence can provide important information for crime investigation. For human hair evidence, the current forensic methods of analysis rely on comparisons of either hair morphology by microscopic examination or nuclear and mitochondrial DNA analyses. Unfortunately in some instances the utilisation of microscopy and DNA analyses are difficult and often not feasible. This dissertation is arguably the first comprehensive investigation aimed to compare, classify and identify the single human scalp hair fibres with the aid of FTIR-ATR spectroscopy in a forensic context. Spectra were collected from the hair of 66 subjects of Asian, Caucasian and African (i.e. African-type). The fibres ranged from untreated to variously mildly and heavily cosmetically treated hairs. The collected spectra reflected the physical and chemical nature of a hair from the near-surface particularly, the cuticle layer. In total, 550 spectra were acquired and processed to construct a relatively large database. To assist with the interpretation of the complex spectra from various types of human hair, Derivative Spectroscopy and Chemometric methods such as Principal Component Analysis (PCA), Fuzzy Clustering (FC) and Multi-Criteria Decision Making (MCDM) program; Preference Ranking Organisation Method for Enrichment Evaluation (PROMETHEE) and Geometrical Analysis for Interactive Aid (GAIA); were utilised. FTIR-ATR spectroscopy had two important advantages over to previous methods: (i) sample throughput and spectral collection were significantly improved (no physical flattening or microscope manipulations), and (ii) given the recent advances in FTIR-ATR instrument portability, there is real potential to transfer this work.s findings seamlessly to on-field applications. The "raw" spectra, spectral subtractions and second derivative spectra were compared to demonstrate the subtle differences in human hair. SEM images were used as corroborative evidence to demonstrate the surface topography of hair. It indicated that the condition of the cuticle surface could be of three types: untreated, mildly treated and treated hair. Extensive studies of potential spectral band regions responsible for matching and discrimination of various types of hair samples suggested the 1690-1500 cm^-1 IR spectral region was to be preferred in comparison with the commonly used 1750-800 cm^-1. The principal reason was the presence of the highly variable spectral profiles of cystine oxidation products (1200-1000 cm^-1), which contributed significantly to spectral scatter and hence, poor hair sample matching. In the preferred 1690-1500 cm^-1 region, conformational changes in the keratin protein attributed to the α-helical to β-sheet transitions in the Amide I and Amide II vibrations and played a significant role in matching and discrimination of the spectra and hence, the hair fibre samples. For gender comparison, the Amide II band is significant for differentiation. The results illustrated that the male hair spectra exhibit a more intense β-sheet vibration in the Amide II band at approximately 1511 cm^-1 whilst the female hair spectra displayed more intense α-helical vibration at 1520-1515cm^-1. In terms of chemical composition, female hair spectra exhibit greater intensity of the amino acid tryptophan (1554 cm^-1), aspartic and glutamic acid (1577 cm^-1). It was also observed that for the separation of samples based on racial differences, untreated Caucasian hair was discriminated from Asian hair as a result of having higher levels of the amino acid cystine and cysteic acid. However, when mildly or chemically treated, Asian and Caucasian hair fibres are similar, whereas African-type hair fibres are different. In terms of the investigation's novel contribution to the field of forensic science, it has allowed for the development of a novel, multifaceted, methodical protocol where previously none had existed. The protocol is a systematic method to rapidly investigate unknown or questioned single human hair FTIR-ATR spectra from different genders and racial origin, including fibres of different cosmetic treatments. Unknown or questioned spectra are first separated on the basis of chemical treatment i.e. untreated, mildly treated or chemically treated, genders, and racial origin i.e. Asian, Caucasian and African-type. The methodology has the potential to complement the current forensic analysis methods of fibre evidence (i.e. Microscopy and DNA), providing information on the morphological, genetic and structural levels.

Preparative deracemization of unnatural amino acids

Unnatural amino acids are a growing class of intermediates required for pharmaceuticals, agrochemicals and other industrial products. However, no single method has proven sufficiently versatile to prepare these compounds broadly at scale. To address this need, we have developed a general chemoenzymatic process to prepare enantiomerically pure L- and D-amino acids in high yield by deracemization of racemic starting materials. This method involves the concerted action of an enantioselective oxidase biocatalyst and a non-selective chemical reducing agent to effect the stereoinversion of one enantiomer and can result in an enantiomeric excess of >99% from the starting racemate, and product yields of over 90%. This approach compares very favourably with resolution processes, which have a maximum single-pass yield of 50%. We have developed efficient methods to adapt the process towards new target compounds and to optimize key factors that influence process efficiency and offer competitive economics at scale.

Novel ANKH amino terminus mutation (Pro5Ser) associated with early-onset calcium pyrophosphate disease with associated phosphaturia

This report describes a 32-year-old woman presenting since childhood with progressive calcium pyrophosphate disease (CPPD), characterized by severe arthropathy and chondrocalcinosis involving multiple peripheral joints and intervertebral disks. Because ANKH mutations have been previously described in familial CPPD, the proband's DNA was assessed at this locus by direct sequencing of promoter and coding regions and revealed 3 sequence variants in ANKH. Sequences of exon 1 revealed a novel isolated nonsynonymous mutation (c.13 C>T), altering amino acid in codon 5 from proline to serine (CCG>TCG). Sequencing of parental DNA revealed an identical mutation in the proband's father but not the mother. Subsequent clinical evaluation demonstrated extensive chondrocalcinosis and degenerative arthropathy in the proband's father. In summary, we report a novel mutation, not previously described, in ANKH exon 1, wherein serine replaces proline, in a case of early-onset severe CPPD associated with metabolic abnormalities, with similar findings in the proband's father.

Orchestration of Haemophilus influenzae RecJ Exonuclease by Interaction with Single-Stranded DNA-Binding Protein

RecJ exonuclease plays crucial roles in several DNA repair and recombination pathways, and its ubiquity in bacterial species points to its ancient origin and vital cellular function. RecJ exonuclease from Haemophilus influenzae is a 575-amino-acid protein that harbors the characteristic motifs conserved among RecJ homologs. The purified protein exhibits a process 5'-3' single-stranded-DNA-specific exonuclease activity. The exonuclease activity of H. influenzae RecJ (HiRecJ) was supported by Mg2+ or Mn2+ and inhibited by Cd2+ suggesting a different mode of metal binding in HiRecJ as compared to Escherichia coli RecJ (EcoRecJ). Site-directed mutagenesis of highly conserved residues in HiRecJ abolished enzymatic activity. Interestingly, substitution of alanine for aspartate 77 resulted in a catalytically inactive enzyme that bound to DNA with a significantly higher affinity as compared to the wild-type enzyme. Noticeably, steady-state kinetic studies showed that H. influenzae single-stranded DNA-binding protein (HiSSB) increased the affinity of HiRecJ for single-stranded DNA and stimulated its exonuclease activity. HiSSB, whose C-terminal tail had been deleted, failed to enhance RecJ exonuclease activity. More importantly, HiRecJ was found to directly associate with its cognate single-stranded DNA-binding protein (SSB), as demonstrated by various in vitro assays, Interaction studies carried out with the truncated variants of HiRecJ and HiSSB revealed that the two proteins interact via the C-terminus of SSB protein and the core-catalytic domain of RecJ. Taken together, these results emphasize direct interactio between RecJ and SSB, which confers functional cooperativity to these two proteins. In addition, these results implicate SSB as being involved in the recruitment of RecJ to DNA and provide insights into the interplay between these proteins in repair and recombination pathways.

X-ray studies on crystalline complexes involving amino acids and peptides. Part XIV: Closed conformation and head-to-tail arrangement in a new crystal form of L-histidine L-aspartate monohydrate

A new form of L-histidine L-aspartate monohydrate crystallizes in space group P22 witha = 5.131(1),b = 6.881(1),c= 18.277(2) Å,β= 97.26(1)° and Z = 2. The structure has been solved by the direct methods and refined to anR value of 0.044 for 1377 observed reflections. Both the amino acid molecules in the complex assume the energetically least favourable allowed conformation with the side chains staggered between the α-amino and α-scarboxylate groups. This results in characteristic distortions in some bond angles. The unlike molecules aggregate into alternating double layers with water molecules sandwiched between the two layers in the aspartate double layer. The molecules in each layer are arranged in a head-to-tail fashion. The aggregation pattern in the complex is fundamentally similar to that in other binary complexes involving commonly occurring L amino acids, although the molecules aggregate into single layers in them. The distribution of crystallographic (and local) symmetry elements in the old form of the complex is very different from that in the new form. So is the conformation of half the histidine molecules. Yet, the basic features of molecular aggregation, particularly the nature and the orientation of head-to-tail sequences, remain the same in both the forms. This supports the thesis that the characteristic aggregation patterns observed in crystal structures represent an intrinsic property of amino acid aggregation.

Studies on crystalline complexes involving amino acids. V. The structure of L-serine-L-ascorbic acid

L-Serine-L-ascorbic acid, C3HTNOa. C6HsO6, a 1:1 complex between the amino acid serine and the vitamin ascorbic acid, crystallizes in the orthorhombic space group P2~2~2~ with four formula units in a cell of dimensions a = 5.335(3), b = 8.769(2), c = 25.782 (5) A. The structure was solved by direct methods and refined by full-matrix least squares to an R of 0.036 for 951 observed reflections. Both molecules are neutral in the structure. The conformation of the serine molecule is different from that observed in the crystal structures of L-serine, DL-serine and L-serine monohydrate. The enediol group in the ascorbic acid molecule is planar, whereas significant departures from planarity are observed in the lactone group. The conformation of this molecule is similar to that observed in arginine ascorbate. The unlike molecules aggregate into separate columns in the crystal structure. The columns are held together by hydrogen bonds. Among these, a pair of hydrogen bonds between the enediol group of ascorbic acid and the carboxylate group of serine provides a possible model for a specific interaction between ascorbic acid and a carboxylate ion.

DNA cleavage in red light promoted by copper(II) complexes of -amino acids and photoactive phenanthroline bases

Ternary copper(II) complexes [Cu(L-trp)(B)(H2O)](NO3) ( 1–3) and [Cu(L-phe)(B)(H2O)](NO3) ( 4–6) of L-tryptophan (L-trp) and L-phenylalanine (L-phe) having phenanthroline bases (B), viz. 1,10-phenanthroline (phen, 1 and 4), dipyrido[3,2-d:2,3-f]quinoxaline (dpq, 2 and 5) and dipyrido[3,2-a:2,3-c]phenazine (dppz, 3 and 6), were prepared and characterized by physico-chemical techniques. Complexes 3 and 6 were structurally characterized by X-ray crystallography and show the presence of a square pyramidal (4 + 1) CuN3O2 coordination geometry in which the N,O-donor amino acid (L-trp or L-phe) and N,N-donor phenanthroline base bind at the equatorial plane with an aqua ligand coordinated at the elongated axial site. Complex 3 shows significant distortion from the square pyramidal geometry and a strong intramolecular – stacking interaction between the pendant indole ring of L-trp and the planar dppz aromatic moiety. All the complexes display good binding propensity to the calf thymus DNA giving an order: 3, 6 (dppz) > 2, 5 (dpq) > 1, 4 (phen). The binding constant (Kb) values are in the range of 2.1 × 104–1.1 × 106 mol-1 with the binding site size (s) values of 0.17–0.63. The phen and dpq complexes are minor groove binders while the dppz analogues bind at the DNA major groove. Theoretical DNA docking studies on 2 and 3 show the close proximity of two photosensitizers, viz. the indole moiety of L-trp and the quinoxaline/phenazine of the dpq/dppz bases, to the complementary DNA strands. Complexes 2 and 3 show oxidative DNA double strand breaks (dsb) of supercoiled (SC) DNA forming a significant quantity of linear DNA along with the nicked circular (NC) form on photoexposure to UV-A light of 365 nm and red light of 647.1 nm (Ar–Kr laser). Complexes 1, 5 and 6 show only single strand breaks (ssb) forming NC DNA. The red light induced DNA cleavage involves metal-assisted photosensitization of L-trp and dpq/dppz base resulting in the formation of a reactive singlet oxygen (1O2) species.

Catalysis of tRNA Aminoacylation: Single Turnover to Steady-State Kinetics of tRNA Synthetases

Aminoacyl-tRNA synthetases (aaRS) catalyze the bimolecular association reaction between amino acid and tRNA by specifically and unerringly choosing the cognate amino acid and tRNA. There are two classes of such synthetases that perform tRNA-aminoacylation reaction. Interestingly, these two classes of aminoacyl-tRNA synthetases differ not only in their structures but they also exhibit remarkably distinct kinetics under pre-steady-state condition. The class I synthetases show initial burst of product formation followed by a slower steady-state rate. This has been argued to represent the influence of slow product release. In contrast, there is no burst in the case of class H enzymes. The tight binding of product with enzyme for class I enzymes is correlated with the enhancement of rate in presence of elongation factor. EF-TU. In spite of extensive experimental studies, there is no detailed theoretical analysis that can provide a quantitative understanding of this important problem. In this article, we present a theoretical investigation of enzyme kinetics for both classes of aminoacyl-tRNA synthetases. We present an augmented kinetic scheme and then employ the methods of time-dependent probability statistics to obtain expressions for the first passage time distribution that gives both the time-dependent and the steady-state rates. The present study quantitatively explains all the above experimental observations. We propose an alternative path way in the case of class II enzymes showing the tRNA-dependent amino acid activation and the discrepancy between the single-turnover and steady-state rate.

Kinetic proofreading at single molecular level: aminoacylation of tRNA(Ile) and the role of water as an editor

Proofreading/editing in protein synthesis is essential for accurate translation of information from the genetic code. In this article we present a theoretical investigation of efficiency of a kinetic proofreading mechanism that employs hydrolysis of the wrong substrate as the discriminatory step in enzyme catalytic reactions. We consider aminoacylation of tRNA(Ile) which is a crucial step in protein synthesis and for which experimental results are now available. We present an augmented kinetic scheme and then employ methods of stochastic simulation algorithm to obtain time dependent concentrations of different substances involved in the reaction and their rates of formation. We obtain the rates of product formation and ATP hydrolysis for both correct and wrong substrates (isoleucine and valine in our case, respectively), in single molecular enzyme as well as ensemble enzyme kinetics. The present theoretical scheme correctly reproduces (i) the amplitude of the discrimination factor in the overall rates between isoleucine and valine which is obtained as (1.8x10(2)).(4.33x10(2)) = 7.8x10(4), (ii) the rates of ATP hydrolysis for both Ile and Val at different substrate concentrations in the aminoacylation of tRNA(Ile). The present study shows a non-michaelis type dependence of rate of reaction on tRNA(Ile) concentration in case of valine. The overall editing in steady state is found to be independent of amino acid concentration. Interestingly, the computed ATP hydrolysis rate for valine at high substrate concentration is same as the rate of formation of Ile-tRNA(Ile) whereas at intermediate substrate concentration the ATP hydrolysis rate is relatively low. We find that the presence of additional editing domain in class I editing enzyme makes the kinetic proofreading more efficient through enhanced hydrolysis of wrong product at the editing CP1 domain.

Computer simulation of human iodothyronine deiodinase single-chain antibody

Iodothyronine plays a major role in growth, basic metabolism and organ formation. It has an extremely limited source in the body. In this thesis, we designed iodothyronine(T4) as hapten. Then a single chain antibody displayed on phange was obtained from a human phage displaying a single chain antibody library. The specific genes of E3 was subcloned in P-5E vector. According to its amino acid sequences, we simulate its three dimention structure by computer. It has never been reported in PDB.