Spondyloarthritis: Matrix Metalloproteinasesas Biomarkers of Pathogenesis and Response to Tumor Necrosis Factor (TNF) Inhibitors

The term spondyloarthritis (SpA) is used to describe a group of multifactorial chronic inflammatory diseases characterized by a predisposing genetic background and clinical manifestations typically involving the sacroiliac joint. The absence of pathognomonic clinical and/or laboratory findings generally results in a delay in diagnosis and, consequently, in treatment. In addition, 20-40% of SpA patients are non-responders to tumor necrosis factor (TNF) inhibitor therapies. Given these considerations, it is important to identify biomarkers that can facilitate the diagnosis and assessment of disease activity. As inflammation plays a key role in the pathogenesis of SpA, inflammatory mediators have been investigated as potential biomarkers for diagnosing the disease and predicting response to therapy. Some investigators have focused their attention on the role of matrix metalloproteinases (MMPs), which are known to be markers of synovial inflammation that is generated in the joint in reaction to inflammatory stimuli. Several studies have been carried out to verify if serum MMPs levels could be useful to diagnose SpA, to assess disease severity, and to predict response to TNF inhibitor therapy. The current review focuses on MMPs' role in SpA pathogenesis, diagnosis and therapeutic implications

A minimalist chemical model of matrix metalloproteinases- Can small peptides mimic the more rigid metal binding sites of proteins?

In order to develop a minimalist chemical model of matrix metalloproteinases (MMPs), we synthesized a pentadecapeptide (Ac-KAHEFGHSLGLDHSK-NH2) corresponding to the catalytic zinc(II) binding site of human MMP-13. The multi-domain structural organization of MMPs fundamentally determines their metal binding affinity, catalytic activity and selectivity. Our potentiometric, UV-VIS, CD, EPR, NMR, ESI-MS and kinetic study are aimed to explore the usefulness of flexible peptides to mimic the more rigid metal binding sites of proteins, to examine the intrinsic metal binding properties of this naked sequence, as well as to contribute the development of a minimalist, peptide-based chemical model of MMPs, including the catalytic properties. Since multiimidazole environment is also characteristic for copper(II), and recently copper(II) containing variants of MMPs have been identified, we also studied the copper(II) complexes of the above peptide. Around pH 6-7 the peptide, similarly to MMPs, offers {3Nim} coordinated binding site for both zinc(II) and copper(II). In the case of copper(II), the formation of amide coordinated species at higher pH ceased the analogy with the copper(II) containing MMP variant. On the other hand, the zinc(II)-peptide system mimics some basic features of the MMP active sites: the main species around pH 7 (ZnH2L) possesses {3Nim,H2O} coordination environment, the deprotonation of the zinc-bound water takes place near to the physiological pH, it forms relatively stable ternary complexes with hydroxamic acids, and the species ZnH2L(OH) and ZnH2L(OH)2 have notable hydrolytic activity between pH 7-9

OF MICE AND MEN: DISSECTING THE INTERACTION BETWEEN LISTERIA MONOCYTOGENES INTERNALIN A AND E-CADHERIN

We report a study of the interaction between internalin A (inlA) and human or murine E-cadherin (Ecad). inlA is used by Listeria monocytogenes to internalize itself into host cell, but the bacterium is unable to invade murine cells, which has been attributed to the difference in sequence between hEcad and mEcad. Using molecular dynamics simulations, MM/GBSA free energy calculations, hydrogen bond analysis, water characterization and umbrella sampling, we provide a complete atomistic picture of the binding between inlA and Ecad. We dissect key residues in the protein–protein interface and analyze the energetics using MM/GBSA.'From this analysis it is clear that the binding of inlA–mEcad is weaker than inlA–hEcad, on par with the experimentally observed inability of inlA to bind to mEcad.'However, extended MD simulations of 200 ns in length show no destabilization of the inlA–mEcad complex and the estimation of the potential of mean force (PMF) using umbrella sampling corroborates this conclusion. The binding strength computed from the PMFs show no significant difference between the two protein complexes. Hence, our study suggests that the inability of L. monocytogenes to invade murine cells cannot be explained by processes at the nanosecond to sub-microsecond time scale probed by the simulations performed here

Variational Methods for Biomolecular Modeling

Structure, function and dynamics of many biomolecular systems can be characterized by the energetic variational principle and the corresponding systems of partial differential equations (PDEs). This principle allows us to focus on the identification of essential energetic components, the optimal parametrization of energies, and the efficient computational implementation of energy variation or minimization. Given the fact that complex biomolecular systems are structurally non-uniform and their interactions occur through contact interfaces, their free energies are associated with various interfaces as well, such as solute-solvent interface, molecular binding interface, lipid domain interface, and membrane surfaces. This fact motivates the inclusion of interface geometry, particular its curvatures, to the parametrization of free energies. Applications of such interface geometry based energetic variational principles are illustrated through three concrete topics: the multiscale modeling of biomolecular electrostatics and solvation that includes the curvature energy of the molecular surface, the formation of microdomains on lipid membrane due to the geometric and molecular mechanics at the lipid interface, and the mean curvature driven protein localization on membrane surfaces. By further implicitly representing the interface using a phase field function over the entire domain, one can simulate the dynamics of the interface and the corresponding energy variation by evolving the phase field function, achieving significant reduction of the number of degrees of freedom and computational complexity. Strategies for improving the efficiency of computational implementations and for extending applications to coarse-graining or multiscale molecular simulations are outlined.Comment: 36 page

Identification of the ADPR binding pocket in the NUDT9 homology domain of TRPM2

Activation of the transient receptor potential melastatin 2 (TRPM2) channel occurs during the response to oxidative stress under physiological conditions as well as in pathological processes such as ischemia and diabetes. Accumulating evidence indicates that adenosine diphosphate ribose (ADPR) is the most important endogenous ligand of TRPM2. However, although it is known that ADPR binds to the NUDT9 homology (NUDT9-H) domain in the intracellular C-terminal region, the molecular mechanism underlying ADPR binding and activation of TRPM2 remains unknown. In this study, we generate a structural model of the NUDT9-H domain and identify the binding pocket for ADPR using induced docking and molecular dynamics simulation. We find a subset of 11 residues—H1346, T1347, T1349, L1379, G1389, S1391, E1409, D1431, R1433, L1484, and H1488—that are most likely to directly interact with ADPR. Results from mutagenesis and electrophysiology approaches support the predicted binding mechanism, indicating that ADPR binds tightly to the NUDT9-H domain, and suggest that the most significant interactions are the van der Waals forces with S1391 and L1484, polar solvation interaction with E1409, and electronic interactions (including π–π interactions) with H1346, T1347, Y1349, D1431, and H1488. These findings not only clarify the roles of a range of newly identified residues involved in ADPR binding in the TRPM2 channel, but also reveal the binding pocket for ADPR in the NUDT9-H domain, which should facilitate structure-based drug design for the TRPM2 channel

Transferability of conformational dependent charges from protein simulations

We have studied the transferability of atomic charges for proteins, fitted to the quantum mechanical electrostatic potential and extensively averaged over a set of structures sampled by molecular dynamics (MD) and over all residues of the same kind in the protein sequence (xAvESP). Previously, such charges were obtained for one single protein (avidin). In this study, we use five additional proteins. The aim of this study is fourfold. First, we provide xAvESP charges for all amino acids, including amino- and carboxy-terminal variants of all, as well as alternative protonation states of His, Asp, Glu, Lys, Arg, Cys, and Tyr. Second, we show that the xAvESP charges averaged over the five new proteins are similar to charges obtained in the same way for avidin, with a correlation coefficient of 0.997. This shows that the charges are transferable and system-independent. Electrostatic proteinligand interaction energies calculated with charges obtained from different proteins differ by only 13 kJ/mol on average. The xAvESP charges correlate rather well with Amber charges (except for the N atom of amino-terminal residues that are erroneous in Amber), although they are obtained in a more general way. Third, the conformational dependence of the charges is significant and gives rise to quite large differences in energies. However, these differences are to a large extent screened by solvation effects. For example, the solvent-screened electrostatic interaction energy between the protein galectin-3 and five different ligands varies with the charge sets by less than 3 kJ/mol on average. Finally, we show that the xAvESP charges give a comparable root-mean-squared deviation as the Amber charges for the MD simulations of 18 proteinligand complexes, they give comparable or slightly worse backbone N?H order parameters for two galectin-3 complexes, but they give a better correlation between calculated and experimental affinities for the binding of seven biotin analogues to avidin and for nine inhibitors of factor Xa. (c) 2011 Wiley Periodicals, Inc. Int J Quantum Chem 112:17681785, 201

The efficiency of C-4 substituents in activating the -lactam scaffold towards serine proteases and hydroxide ion

The presence of a leaving group at C-4 of monobactams is usually considered to be a requirement for mechanism-based inhibition of human leukocyte elastase by these acylating agents. We report that second-order rate constants for the alkaline hydrolysis and elastase inactivation by N-carbamoyl monobactams are independent of the pKa of the leaving group at C-4. Indeed, the effect exerted by these substituents is purely inductive: electron-withdrawing substituents at C-4 of N-carbamoyl-3,3-diethylmonobactams increase the rate of alkaline hydrolysis and elastase inactivation, with Hammett rho-I values of 3.4 and 2.5, respectively, which indicate the development of a negative charge in the transition-states. The difference in magnitude between these rho-I values is consistent with an earlier transition-state for the enzymatic reaction when compared with that for the chemical process. These results suggest that rate limiting step in elastase inactivation is the formation of the tetrahedral intermediate, and that beta-lactam ring-opening is not concerted with the departure of a leaving group from C-4. Monobactam sulfones emerged as potent elastase inhibitors even when the ethyl groups at C-3, required for interaction with the primary recognition site, are absent. For one such compound, a 1:1 enzyme-inhibitor complex involving porcine pancreatic elastase has been examined by X-ray crystallography and shown to result from serine acylation and sulfinate departure from the β-lactam C-4

Water-Membrane Partition Thermodynamics of an Amphiphilic Lipopeptide: An Enthalpy-Driven Hydrophobic Effect

To shed light on the driving force for the hydrophobic effect that partitions amphiphilic lipoproteins between water and membrane, we carried out an atomically detailed thermodynamic analysis of a triply lipid modified H-ras heptapeptide anchor (ANCH) in water and in a DMPC (1,2-dimyristoyl-sn-glycero-3-phosphocholine) bilayer. Combining molecular mechanical and continuum solvent approaches with an improved technique for solute entropy calculation, we obtained an overall transfer free energy of ∼−13 kcal mol−1. This value is in qualitative agreement with free energy changes derived from a potential of mean force calculation and indirect experimental observations. Changes in free energies of solvation and ANCH conformational reorganization are unfavorable, whereas ANCH-DMPC interactions—especially van der Waals—favor insertion. These results are consistent with an enthalpy-driven hydrophobic effect, in accord with earlier calorimetric data on the membrane partition of other amphiphiles. Furthermore, structural and entropic analysis of molecular dynamics-generated ensembles suggests that conformational selection may play a hitherto unappreciated role in membrane insertion of lipid-modified peptides and proteins

Dissection of the Inflammatory Bowel Disease Transcriptome Using Genome-Wide cDNA Microarrays

BACKGROUND: The differential pathophysiologic mechanisms that trigger and maintain the two forms of inflammatory bowel disease (IBD), Crohn disease (CD), and ulcerative colitis (UC) are only partially understood. cDNA microarrays can be used to decipher gene regulation events at a genome-wide level and to identify novel unknown genes that might be involved in perpetuating inflammatory disease progression. METHODS AND FINDINGS: High-density cDNA microarrays representing 33,792 UniGene clusters were prepared. Biopsies were taken from the sigmoid colon of normal controls (n = 11), CD patients (n = 10) and UC patients (n = 10). (33)P-radiolabeled cDNA from purified poly(A)(+) RNA extracted from biopsies (unpooled) was hybridized to the arrays. We identified 500 and 272 transcripts differentially regulated in CD and UC, respectively. Interesting hits were independently verified by real-time PCR in a second sample of 100 individuals, and immunohistochemistry was used for exemplary localization. The main findings point to novel molecules important in abnormal immune regulation and the highly disturbed cell biology of colonic epithelial cells in IBD pathogenesis, e.g., CYLD (cylindromatosis, turban tumor syndrome) and CDH11 (cadherin 11, type 2). By the nature of the array setup, many of the genes identified were to our knowledge previously uncharacterized, and prediction of the putative function of a subsection of these genes indicate that some could be involved in early events in disease pathophysiology. CONCLUSION: A comprehensive set of candidate genes not previously associated with IBD was revealed, which underlines the polygenic and complex nature of the disease. It points out substantial differences in pathophysiology between CD and UC. The multiple unknown genes identified may stimulate new research in the fields of barrier mechanisms and cell signalling in the context of IBD, and ultimately new therapeutic approaches