Analysis of agonist and antagonist effects on thyroid hormone receptor conformation by hydrogen/deuterium exchange

Thyroid hormone receptors (TRs) are ligand-gated transcription factors with critical roles in development and metabolism. Although x-ray structures of TR ligand-binding domains (LBDs) with agonists are available, comparable structures without ligand (apo-TR) or with antagonists are not. It remains important to understand apo-LBD conformation and the way that it rearranges with ligands to develop better TR pharmaceuticals. In this study, we conducted hydrogen/deuterium exchange on TR LBDs with or without agonist (T 3) or antagonist (NH3). Both ligands reduce deuterium incorporation into LBD amide hydrogens, implying tighter overall folding of the domain. As predicted, mass spectroscopic analysis of individual proteolytic peptides after hydrogen/ deuterium exchange reveals that ligand increases the degree of solvent protection of regions close to the buried ligand-binding pocket. However, there is also extensive ligand protection of other regions, including the dimer surface at H10-H11, providing evidence for allosteric communication between the ligand-binding pocket and distant interaction surfaces. Surprisingly, Cterminal activation helix H12, which is known to alter position with ligand, remains relatively protected from solvent in all conditions suggesting that it is packed against the LBD irrespective of the presence or type of ligand. T 3, but not NH3, increases accessibility of the upper part of H3-H5 to solvent, and we propose that TR H12 interacts with this region in apo-TR and that this interaction is blocked by T 3 but not NH3.Wepresent data from site-directed mutagenesis experiments and molecular dynamics simulations that lend support to this structural model of apo-TR and its ligand-dependent conformational changes. (Molecular Endocrinology 25: 15-31, 2011). Copyright © 2011 by The Endocrine Society.

Pressure-assisted cold denaturation of hen egg white lysozyme: the influence of co-solvents probed by hydrogen exchange nuclear magnetic resonance

COSY proton nuclear magnetic resonance was used to measure the exchange rates of amide protons of hen egg white lysozyme (HEWL) in the pressure-assisted cold-denatured state and in the heat-denatured state. After dissolving lysozyme in deuterium oxide buffer, labile protons exchange for deuterons in such a way that exposed protons are substituted rapidly, whereas "protected" protons within structured parts of the protein are substituted slowly. The exchange rates k obs were determined for HEWL under heat treatment (80ºC) and under high pressure conditions at low temperature (3.75 kbar, -13ºC). Moreover, the influence of co-solvents (sorbitol, urea) on the exchange rate was examined under pressure-assisted cold denaturation conditions, and the corresponding protection factors, P, were determined. The exchange kinetics upon heat treatment was found to be a two-step process with initial slow exchange followed by a fast one, showing residual protection in the slow-exchange state and P-factors in the random-coil-like range for the final temperature-denatured state. Addition of sorbitol (500 mM) led to an increase of P-factors for the pressure-assisted cold denatured state, but not for the heat-denatured state. The presence of 2 M urea resulted in a drastic decrease of the P-factors of the pressure-assisted cold denatured state. For both types of co-solvents, the effect they exert appears to be cooperative, i.e., no particular regions within the protein can be identified with significantly diverse changes of P-factors.

Destabilizing Mutations Alter the Hydrogen Exchange Mechanism in Ribonuclease A

The effect of strongly destabilizing mutations, I106A and V108G of Ribonuclease A (RNase A), on its structure and stability has been determined by NMR. The solution structures of these variants are essentially equivalent to RNase A. The exchange rates of the most protected amide protons in RNase A (35ºC), the I106A variant (35ºC), and the V108G variant (10ºC) yield stability values of 9.9, 6.0, and 6.8 kcal/mol, respectively, when analyzed assuming an EX2 exchange mechanism. Thus, the destabilization induced by these mutations is propagated throughout the protein. Simulation of RNase A hydrogen exchange indicates that the most protected protons in RNase A and the V108G variant exchange via the EX2 regime, whereas those of I106A exchange through a mixed EX1 1 EX2 process. It is striking that a single point mutation can alter the overall exchange mechanism. Thus, destabilizing mutations joins high temperatures, high pH and the presence of denaturating agents as a factor that induces EX1 exchange in proteins. The calculations also indicate a shift from the EX2 to the EX1 mechanism for less protected groups within the same protein. This should be borne in mind when interpreting exchange data as a measure of local stability in less protected regions

Signal Processing for the Measurement of the Deuterium/Hydrogen Ratio in the Local Interstellar Medium

We report on a comprehensive signal processing procedure for very low signal levels for the measurement of neutral deuterium in the local interstellar medium from a spacecraft in Earth orbit. The deuterium measurements were performed with the IBEX-Lo camera on NASA’s Interstellar Boundary Explorer (IBEX) satellite. Our analysis technique for these data consists of creating a mass relation in three-dimensional time of flight space to accurately determine the position of the predicted D events, to precisely model the tail of the H events in the region where the H tail events are near the expected D events, and then to separate the H tail from the observations to extract the very faint D signal. This interstellar D signal, which is expected to be a few counts per year, is extracted from a strong terrestrial background signal, consisting of sputter products from the sensor’s conversion surface. As reference we accurately measure the terrestrial D/H ratio in these sputtered products and then discriminate this terrestrial background source. During the three years of the mission time when the deuterium signal was visible to IBEX, the observation geometry and orbit allowed for a total observation time of 115.3 days. Because of the spinning of the spacecraft and the stepping through eight energy channels the actual observing time of the interstellar wind was only 1.44 days. With the optimised data analysis we found three counts that could be attributed to interstellar deuterium. These results update our earlier work.

A general two-process model describes the hydrogen exchange behavior of RNase A in unfolding conditions.

When NMR hydrogen exchange was used previously to monitor the kinetics of RNase A unfolding, some peptide NH protons were found to show EX2 exchange (detected by base catalysis) in addition to the expected EX1 exchange, whose rate is limited by the kinetic unfolding process. In earlier work, two groups showed independently that a restricted two-process model successfully fits published hydrogen exchange rates of native RNase A in the range 0-0.7 M guanidinium chloride. We find that this model predicts properties that are very different from the observed properties of the EX2 exchange reactions of RNase A in conditions where guanidine-induced unfolding takes place. The model predicts that EX2 exchange should be too fast to measure by the technique used, whereas it is readily measurable. Possible explanations for the contradiction are considered here, and we show that removing the restriction from the earlier two-process model is sufficient to resolve the contradiction; instead of specifying that exchange caused by global unfolding occurs by the EX2 mechanism, we allow it to occur by the general mechanism, which includes both the EX1 and EX2 cases. It is logical to remove this restriction because global unfolding of RNase A is known to give rise to EX1 exchange in these unfolding conditions. Resolving the contradiction makes it possible to determine whether populated unfolding intermediates contribute to the EX2 exchange, and this question is considered elsewhere. The results and simulations indicate that moderate or high denaturant concentrations readily give rise to EX1 exchange in native proteins. Earlier studies showed that hydrogen exchange in native proteins typically occurs by the EX2 mechanism but that high temperatures or pH values above 7 may give rise to EX1 exchange. High denaturant concentrations should be added to the list of variables likely to cause EX1 exchange.

Kinetics of hydrogen bond breakage in the process of unfolding of ribonuclease A measured by pulsed hydrogen exchange.

A sensitive test for kinetic unfolding intermediates in ribonuclease A (EC 3.1.27.5) is performed under conditions where the enzyme unfolds slowly (10 degrees C, pH 8.0, 4.5 M guanidinium chloride). Exchange of peptide NH protons (2H-1H) is used to monitor structural opening of individual hydrogen bonds during unfolding, and kinetic models are developed for hydrogen exchange during the process of protein unfolding. The analysis indicates that the kinetic process of unfolding can be monitored by EX1 exchange (limited by the rate of opening) for ribonuclease A in these conditions. Of the 49 protons whose unfolding/exchange kinetics was measured, 47 have known hydrogen bond acceptor groups. To test whether exchange during unfolding follows the EX2 (base-catalyzed) or the EX1 (uncatalyzed) mechanism, unfolding/exchange was measured both at pH 8.0 and at pH 9.0. A few faster-exchanging protons were found that undergo exchange by both EX1 and EX2 processes, but the 43 slower-exchanging protons at pH 8 undergo exchange only by the EX1 mechanism, and they have closely similar rates. Thus, it is likely that all 49 protons undergo EX1 exchange at the same rate. The results indicate that a single rate-limiting step in unfolding breaks the entire network of peptide hydrogen bonds and causes the overall unfolding of ribonuclease A. The additional exchange observed for some protons that follows the EX2 mechanism probably results from equilibrium unfolding intermediates and will be discussed elsewhere.

Fractionation factors and activation energies for exchange of the low barrier hydrogen bonding proton in peptidyl trifluoromethyl ketone complexes of chymotrypsin

NMR investigations have been carried out of complexes between bovine chymotrypsin Aα and a series of four peptidyl trifluoromethyl ketones, listed here in order of increasing affinity for chymotrypsin: N-Acetyl-l-Phe-CF3, N-Acetyl-Gly-l-Phe-CF3, N-Acetyl-l-Val-l-Phe-CF3, and N-Acetyl-l-Leu-l-Phe-CF3. The D/H fractionation factors (φ) for the hydrogen in the H-bond between His 57 and Asp 102 (His 57-Hδ1) in these four complexes at 5°C were in the range φ = 0.32–0.43, expected for a low-barrier hydrogen bond. For this series of complexes, measurements also were made of the chemical shifts of His 57-Hɛ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 8.97–9.18), the exchange rate of the His 57-Hδ1 proton with bulk water protons (284–12.4 s−1), and the activation enthalpies for this hydrogen exchange (14.7–19.4 kcal⋅mol−1). It was found that the previously noted correlations between the inhibition constants (Ki 170–1.2 μM) and the chemical shifts of His 57-Hδ1 (δ2,2-dimethylsilapentane-5-sulfonic acid 18.61–18.95) for this series of peptidyl trifluoromethyl ketones with chymotrypsin [Lin, J., Cassidy, C. S. & Frey, P. A. (1998) Biochemistry 37, 11940–11948] could be extended to include the fractionation factors, hydrogen exchange rates, and hydrogen exchange activation enthalpies. The results support the proposal of low barrier hydrogen bond-facilitated general base catalysis in the addition of Ser 195 to the peptidyl carbonyl group of substrates in the mechanism of chymotrypsin-catalyzed peptide hydrolysis. Trends in the enthalpies for hydrogen exchange and the fractionation factors are consistent with a strong, double-minimum or single-well potential hydrogen bond in the strongest complexes. The lifetimes of His 57-Hδ1, which is solvent shielded in these complexes, track the strength of the hydrogen bond. Because these lifetimes are orders of magnitude shorter than those of the complexes themselves, the enzyme must have a pathway for hydrogen exchange at this site that is independent of dissociation of the complexes.

Identification of initiation sites for T4 lysozyme folding using CD and NMR spectroscopy of peptide fragments

Using CD and 2D H-1 NMR spectroscopy, we have identified potential initiation sites for the folding of T4 lysozyme by examining the conformational preferences of peptide fragments corresponding to regions of secondary structure. CD spectropolarimetry showed most peptides were unstructured in water, but adopted partial helical conformations in TFE and SDS solution. This was also consistent with the H-1 NMR data which showed that the peptides were predominantly disordered in water, although in some cases, nascent or small populations of partially folded conformations could be detected. NOE patterns, coupling constants, and deviations from random coil Her chemical shift values complemented the CD data and confirmed that many of the peptides were helical in TFE and SDS micelles. In particular, the peptide corresponding to helix E in the native enzyme formed a well-defined helix in both TFE and SDS, indicating that helix E potentially forms an initiation site for T4 lysozyme folding. The data for the other peptides indicated that helices D, F, G, and H are dependent on tertiary interactions for their folding and/or stability. Overall, the results from this study, and those of our earlier studies, are in agreement with modeling and IID-deuterium exchange experiments, and support an hierarchical model of folding for T4 lysozyme.

Molecular modeling study of the binding and signaling properties of the TCRpMHC complex

The activation of the specific immune response against tumor cells is based on the recognition by the CD8+ Cytotoxic Τ Lymphocytes (CTL), of antigenic peptides (p) presented at the surface of the cell by the class I major histocompatibility complex (MHC). The ability of the so-called T-Cell Receptors (TCR) to discriminate between self and non-self peptides constitutes the most important specific control mechanism against infected cells. The TCR/pMHC interaction has been the subject of much attention in cancer therapy since the design of the adoptive transfer approach, in which Τ lymphocytes presenting an interesting response against tumor cells are extracted from the patient, expanded in vitro, and reinfused after immunodepletion, possibly leading to cancer regression. In the last decade, major progress has been achieved by the introduction of engineered lypmhocytes. In the meantime, the understanding of the molecular aspects of the TCRpMHC interaction has become essential to guide in vitro and in vivo studies. In 1996, the determination of the first structure of a TCRpMHC complex by X-ray crystallography revealed the molecular basis of the interaction. Since then, molecular modeling techniques have taken advantage of crystal structures to study the conformational space of the complex, and understand the specificity of the recognition of the pMHC by the TCR. In the meantime, experimental techniques used to determine the sequences of TCR that bind to a pMHC complex have been used intensively, leading to the collection of large repertoires of TCR sequences that are specific for a given pMHC. There is a growing need for computational approaches capable of predicting the molecular interactions that occur upon TCR/pMHC binding without relying on the time consuming resolution of a crystal structure. This work presents new approaches to analyze the molecular principles that govern the recognition of the pMHC by the TCR and the subsequent activation of the T-cell. We first introduce TCRep 3D, a new method to model and study the structural properties of TCR repertoires, based on homology and ab initio modeling. We discuss the methodology in details, and demonstrate that it outperforms state of the art modeling methods in predicting relevant TCR conformations. Two successful applications of TCRep 3D that supported experimental studies on TCR repertoires are presented. Second, we present a rigid body study of TCRpMHC complexes that gives a fair insight on the TCR approach towards pMHC. We show that the binding mode of the TCR is correctly described by long-distance interactions. Finally, the last section is dedicated to a detailed analysis of an experimental hydrogen exchange study, which suggests that some regions of the constant domain of the TCR are subject to conformational changes upon binding to the pMHC. We propose a hypothesis of the structural signaling of TCR molecules leading to the activation of the T-cell. It is based on the analysis of correlated motions in the TCRpMHC structure. - L'activation de la réponse immunitaire spécifique dirigée contre les cellules tumorales est basée sur la reconnaissance par les Lymphocytes Τ Cytotoxiques (CTL), d'un peptide antigénique (p) présenté à la suface de la cellule par le complexe majeur d'histocompatibilité de classe I (MHC). La capacité des récepteurs des lymphocytes (TCR) à distinguer les peptides endogènes des peptides étrangers constitue le mécanisme de contrôle le plus important dirigé contre les cellules infectées. L'interaction entre le TCR et le pMHC est le sujet de beaucoup d'attention dans la thérapie du cancer, depuis la conception de la méthode de transfer adoptif: les lymphocytes capables d'une réponse importante contre les cellules tumorales sont extraits du patient, amplifiés in vitro, et réintroduits après immunosuppression. Il peut en résulter une régression du cancer. Ces dix dernières années, d'importants progrès ont été réalisés grâce à l'introduction de lymphocytes modifiés par génie génétique. En parallèle, la compréhension du TCRpMHC au niveau moléculaire est donc devenue essentielle pour soutenir les études in vitro et in vivo. En 1996, l'obtention de la première structure du complexe TCRpMHC à l'aide de la cristallographie par rayons X a révélé les bases moléculaires de l'interaction. Depuis lors, les techniques de modélisation moléculaire ont exploité les structures expérimentales pour comprendre la spécificité de la reconnaissance du pMHC par le TCR. Dans le même temps, de nouvelles techniques expérimentales permettant de déterminer la séquence de TCR spécifiques envers un pMHC donné, ont été largement exploitées. Ainsi, d'importants répertoires de TCR sont devenus disponibles, et il est plus que jamais nécessaire de développer des approches informatiques capables de prédire les interactions moléculaires qui ont lieu lors de la liaison du TCR au pMHC, et ce sans dépendre systématiquement de la résolution d'une structure cristalline. Ce mémoire présente une nouvelle approche pour analyser les principes moléculaires régissant la reconnaissance du pMHC par le TCR, et l'activation du lymphocyte qui en résulte. Dans un premier temps, nous présentons TCRep 3D, une nouvelle méthode basée sur les modélisations par homologie et ab initio, pour l'étude de propriétés structurales des répertoires de TCR. Le procédé est discuté en détails et comparé à des approches standard. Nous démontrons ainsi que TCRep 3D est le plus performant pour prédire des conformations pertinentes du TCR. Deux applications à des études expérimentales des répertoires TCR sont ensuite présentées. Dans la seconde partie de ce travail nous présentons une étude de complexes TCRpMHC qui donne un aperçu intéressant du mécanisme d'approche du pMHC par le TCR. Finalement, la dernière section se concentre sur l'analyse détaillée d'une étude expérimentale basée sur les échanges deuterium/hydrogène, dont les résultats révèlent que certaines régions clés du domaine constant du TCR sont sujettes à un changement conformationnel lors de la liaison au pMHC. Nous proposons une hypothèse pour la signalisation structurelle des TCR, menant à l'activation du lymphocyte. Celle-ci est basée sur l'analyse des mouvements corrélés observés dans la structure du TCRpMHC.

Stable water isotopes of precipitation and firn cores from the northern Antarctic Peninsula region as a proxy for climate reconstruction

In order to investigate the climate variability in the northern Antarctic Peninsula region, this paper focuses on the relationship between stable isotope content of precipitation and firn, and main meteorological variables (air temperature, relative humidity, sea surface temperature, and sea ice extent). Between 2008 and 2010, we collected precipitation samples and retrieved firn cores from several key sites in this region. We conclude that the deuterium excess oscillation represents a robust indicator of the meteorological variability on a seasonal to sub-seasonal scale. Low absolute deuterium excess values and the synchronous variation of both deuterium excess and air temperature imply that the evaporation of moisture occurs in the adjacent Southern Ocean. The delta O-18-air temperature relationship is complicated and significant only at a (multi)seasonal scale. Backward trajectory calculations show that air-parcels arriving at the region during precipitation events predominantly originate at the South Pacific Ocean and Bellingshausen Sea. These investigations will be used as a calibration for ongoing and future research in the area, suggesting that appropriate locations for future ice core research are located above 600 m a.s.l. We selected the Plateau Laclavere, Antarctic Peninsula as the most promising site for a deeper drilling campaign.

Identification of protein–protein interfaces by decreased amide proton solvent accessibility

Matrix-assisted laser desorption ionization–time-of-flight mass spectrometry was used to identify peptic fragments from protein complexes that retained deuterium under hydrogen exchange conditions due to decreased solvent accessibility at the interface of the complex. Short deuteration times allowed preferential labeling of rapidly exchanging surface amides so that primarily solvent accessibility changes and not conformational changes were detected. A single mass spectrum of the peptic digest mixture was analyzed to determine the deuterium content of all proteolytic fragments of the protein. The protein–protein interface was reliably indicated by those peptides that retained more deuterons in the complex compared with control experiments in which only one protein was present. The method was used to identify the kinase inhibitor [PKI(5–24)] and ATP-binding sites in the cyclic-AMP-dependent protein kinase. Three overlapping peptides identified the ATP-binding site, three overlapping peptides identified the glycine-rich loop, and two peptides identified the PKI(5–24)-binding site. A complex of unknown structure also was analyzed, human α-thrombin bound to an 83-aa fragment of human thrombomodulin [TMEGF(4–5)]. Five peptides from thrombin showed significantly decreased solvent accessibility in the complex. Three peptides identified the anion-binding exosite I, confirming ligand competition experiments. Two peptides identified a new region of thrombin near the active site providing a potential mechanism of how thrombomodulin alters thrombin substrate specificity.

Acyclic permutants of naturally occurring cyclic proteins - Characterization of cystine knot and beta-sheet formation in the macrocyclic polypeptide kalata B1

Kalata B1 is a prototypic member of the unique cyclotide family of macrocyclic polypeptides in which the major structural features are a circular peptide backbone, a triple stranded beta-sheet, and a cystine knot arrangement of three disulfide bonds. The cyclotides are the only naturally occurring family of circular proteins and have prompted us to explore the concept of acyclic permutation, i.e. opening the backbone of a cross-linked circular protein in topologically permuted ways. We have synthesized the complete suite of acyclic permutants of kalata B1 and examined the effect of acyclic permutation on structure and activity. Only two of six topologically distinct backbone loops are critical for folding into the native conformation, and these involve disruption of the embedded ring in the cystine knot. Surprisingly, it is possible to disrupt regions of the p-sheet and still allow folding into native-like structure, provided the cystine knot is intact. Kalata B1 has mild hemolytic activity, but despite the overall structure of the native peptide being retained in all but two cases, none of the acyclic permutants displayed hemolytic activity. This loss of activity is not localized to one particular region and suggests that cyclization is critical for hemolytic activity.

Stabilization of native protein fold by intein-mediated covalent cyclization

A mutant version of the N-terminal domain of Escherichia coli DnaB helicase was used as a model system to assess the stabilization against unfolding gained by covalent cyclization. Cyclization was achieved in vivo by formation of an amide bond between the N and C termini with the help of a split mini-intein. Linear and circular proteins were constructed to be identical in amino acid sequence. Mutagenesis of Phe102 to Glu rendered the protein monomeric even at high concentration. A difference in free energy of unfolding, DeltaDeltaG, between circular and linear protein of 2.3(+/-0.5) kcal mol(-1) was measured at 10degreesC by circular dichroism. A theoretical estimate of the difference in conformational entropy of linear and circular random chains in a three-dimensional cubic lattice model predicted DeltaDeltaG = 2.3 kcal mol(-1), suggesting that stabilization by protein cyclization is driven by the reduced conformational entropy of the unfolded state. Amide-proton exchange rates measured by NMR spectroscopy and mass spectrometry showed a uniform, approximately tenfold decrease of the exchange rates of the most slowly exchanging amide protons, demonstrating that cyclization globally decreases the unfolding rate of the protein. The amide proton exchange was found to follow EX1 kinetics at near-neutral pH, in agreement with an unusually slow refolding I measured by stopped-flow circular dichroism. rate of less than 4 min(-1) The linear and circular proteins differed more in their unfolding than in their folding rates. Global unfolding of the N-terminal domain of E. coli DnaB is thus promoted strongly by spatial separation of the N and C termini, whereas their proximity is much less important for folding. (C) 2005 Elsevier Ltd. All rights reserved.

Seasonal changes in the D/H ratio of fatty acids of pelagic microorganisms

Culture studies of microorganisms have shown that the hydrogen isotopic composition of fatty acids depends on their metabolism, but there are only few environmental studies available to confirm this observation. Here we studied the seasonal variability of the deuterium/hydrogen (D/H) ratio of fatty acids in the coastal Dutch North Sea and compared this with the diversity of the phyto- and bacterioplankton. Over the year, the stable hydrogen isotopic fractionation factor epsilon between fatty acids and water ranged between -172 per mil and -237 per mil, the algal-derived polyunsaturated fatty acid nC20:5 being the most D-depleted and nC18:0 the least D-depleted fatty acid. The D-depleted nC20:5 is in agreement with culture studies, which indicates that photoautotrophic microorganisms produce fatty acids which are significantly depleted in D relative to water. The epsilon-lipid/water of all fatty acids showed a transient shift towards increased fractionation during the spring phytoplankton bloom, indicated by increasing chlorophyll a concentrations and relative abundance of the nC20:5 PUFA, suggesting increased contributions of photoautotrophy. Time periods with decreased fractionation (less negative epsilon-lipid/water values) can be explained by an increased contribution by heterotrophy to the fatty acid pool. Our results show that the hydrogen isotopic composition of fatty acids is a useful tool to assess the community metabolism of coastal plankton.

Stability and dynamics in a hyperthermophilic protein with melting temperature close to 200°C

The rubredoxin protein from the hyperthermophilic archaebacterium Pyrococcus furiosus was examined by a hydrogen exchange method. Even though the protein does not exhibit reversible thermal unfolding, one can determine its stability parameters—free energy, enthalpy, entropy, and melting temperature—and also the distribution of stability throughout the protein, by using hydrogen exchange to measure the reversible cycling of the protein between native and unfolded states that occurs even under native conditions.