CO rebinding kinetics and molecular dynamics simulations highlight dynamic regulation of internal cavities in human cytoglobin

Cytoglobin (Cygb) was recently discovered in the human genome and localized in different tissues. It was suggested to play tissue-specific protective roles, spanning from scavenging of reactive oxygen species in neurons to supplying oxygen to enzymes in fibroblasts. To shed light on the functioning of such versatile machinery, we have studied the processes supporting transport of gaseous heme ligands in Cygb. Carbon monoxide rebinding shows a complex kinetic pattern with several distinct reaction intermediates, reflecting rebinding from temporary docking sites, second order recombination, and formation (and dissociation) of a bis-histidyl heme hexacoordinated reaction intermediate. Ligand exit to the solvent occurs through distinct pathways, some of which exploit temporary docking sites. The remarkable change in energetic barriers, linked to heme bis-histidyl hexacoordination by HisE7, may be responsible for active regulation of the flux of reactants and products to and from the reaction site on the distal side of the heme. A substantial change in both protein dynamics and inner cavities is observed upon transition from the CO-liganded to the pentacoordinated and bis-histidyl hexacoordinated species, which could be exploited as a signalling state. These findings are consistent with the expected versatility of the molecular activity of this protein.

Identification of the molecular site of ivabradine binding to HCN4 channels

Ivabradine is a specific heart rate-reducing agent approved as a treatment of chronic stable angina. Its mode of action involves a selective and specific block of HCN channels, the molecular components of sinoatrial "funny" (f)-channels. Different studies suggest that the binding site of ivabradine is located in the inner vestibule of HCN channels, but the molecular details of ivabradine binding are unknown. We thus sought to investigate by mutagenesis and in silico analysis which residues of the HCN4 channel, the HCN isoform expressed in the sinoatrial node, are involved in the binding of ivabradine. Using homology modeling, we verified the presence of an inner cavity below the channel pore and identified residues lining the cavity; these residues were replaced with alanine (or valine) either alone or in combination, and WT and mutant channels were expressed in HEK293 cells. Comparison of the block efficiency of mutant vs WT channels, measured by patch-clamp, revealed that residues Y506, F509 and I510 are involved in ivabradine binding. For each mutant channel, docking simulations correctly explain the reduced block efficiency in terms of proportionally reduced affinity for ivabradine binding. In summary our study shows that ivabradine occupies a cavity below the channel pore, and identifies specific residues facing this cavity that interact and stabilize the ivabradine molecule. This study provides an interpretation of known properties of f/HCN4 channel block by ivabradine such as the “open channel block”, the current-dependence of block and the property of "trapping" of drug molecules in the closed configuration.

Structural insight into inhibitor of apoptosis proteins recognition by a potent divalent smac-mimetic

Genetic alterations enhancing cell survival and suppressing apoptosis are hallmarks of cancer that significantly reduce the efficacy of chemotherapy or radiotherapy. The Inhibitor of Apoptosis Protein (IAP) family hosts conserved proteins in the apoptotic pathway whose over-expression, frequently found in tumours, potentiates survival and resistance to anticancer agents. In humans, IAPs comprise eight members hosting one or more structural Baculoviral IAP Repeat (BIR) domains. Cellular IAPs (cIAP1 and 2) indirectly inhibit caspase-8 activation, and regulate both the canonical and the non-canonical NF-κB signaling pathways. In contrast to cIAPs, XIAP (X chromosome-linked Inhibitor of Apoptosis Protein) inhibits directly the effector caspases-3 and -7 through its BIR2 domain, and initiator caspase-9 through its BIR3 domain; molecular docking studies suggested that Smac/DIABLO antagonizes XIAP by simultaneously targeting both BIR2 and BIR3 domains. Here we report analytical gel filtration, crystallographic and SAXS experiments on cIAP1-BIR3, XIAP-BIR3 and XIAP-BIR2BIR3 domains, alone and in the presence of compound 9a, a divalent homodimeric Smac mimetic. 9a is shown to bind two BIR domains inter- (in the case of two BIR3) and intra-molecularly (in the case of XIAP-BIR2BIR3), with higher affinity for cIAP1-BIR3, relative to XIAP-BIR3. Despite the different crystal lattice packing, 9a maintains a right handed helical conformation in both cIAP1-BIR3 and XIAP-BIR3 crystals, that is likely conserved in solution as shown by SAXS data. Our structural results demonstrate that the 9a linker length, its conformational degrees of freedom and its hydrophobicity, warrant an overall compact structure with optimal solvent exposure of its two active moieties for IAPs binding. Our results show that 9a is a good candidate for pre-clinical and clinical studies, worth of further investigations in the field of cancer therapy.

A structure-based strategy for epitope discovery in Burkholderia pseudomallei OppA antigen

We present an approach integrating structural and computational biology with immunological tests to identify epitopes in the OppA antigen from the Gram-negative pathogen Burkholderia pseudomallei, the etiological agent of melioidosis. The crystal structure of OppA(Bp), reported here at 2.1 Å resolution, was the basis for a computational analysis that identified three potential epitopes. In parallel, antigen proteolysis and immunocapturing allowed us to identify three additional peptides. All six potential epitopes were synthesized as free peptides and tested for their immunoreactivity against sera from healthy seronegative, healthy seropositive, and recovered melioidosis patients. Three synthetic peptides allowed the different patient groups to be distinguished, underlining the potential of this approach. Extension of the computational analysis, including energy-based decomposition methods, allowed rationalizing results of the predictive analyses and the immunocapture epitope mapping. Our results illustrate a structure-based epitope discovery process, whose application may expand our perspectives in the diagnostic and vaccine design fields.

Quantitative expression of the mutated lamin A/C gene in patients with cardiolaminopathy

OBJECTIVES

The authors sought to investigate the gene and protein expression in Lamin A/C (LMNA)-mutated dilated cardiolaminopathy (DCM) patients (DCM(LMNAMut)) versus LMNA-wild-type DCM (DCM(LMNAWT)), and normal controls (CTRL(LMNAWT)).

BACKGROUND

Dilated cardiolaminopathies are clinically characterized by high arrhythmogenic risk and caused by LMNA mutations. Little is known regarding quantitative gene expression (QGE) of the LMNA gene in blood and myocardium, as well as regarding myocardial expression of the lamin A/C protein.

METHODS

Using the comparative ΔΔCT method, we evaluated the QGE of LMNA (QGE(LMNA)) in peripheral blood and myocardial RNA from carriers of LMNA mutations, versus blood and myocardial samples from DCM(LMNAWT) patients and CTRL(LMNAWT) individuals. After generating reference values in normal controls, QGE(LMNA) was performed in 311 consecutive patients and relatives, blind to genotype, to assess the predictive value of QGE(LMNA) for the identification of mutation carriers. In parallel, Lamin A/C was investigated in myocardial samples from DCM(LMNAMut) versus DCM(LMNAWT) versus normal hearts (CTRL(LMNAWT)).

RESULTS

LMNA was significantly underexpressed in mRNA from peripheral blood and myocardium of DCM(LMNAMut) patients versus DCM(LMNAWT) and CTRL(LMNAWT). In 311 individuals, blind to genotype, the QGE(LMNA) showed 100% sensitivity and 87% specificity as a predictor of LMNA mutations. The receiver-operating characteristic curve analysis yielded an area under the curve of 0.957 (p < 0.001). Loss of protein in cardiomyocytes' nuclei was documented in DCM(LMNAMut) patients.

CONCLUSIONS

The reduced expression of LMNA gene in blood is a novel potential predictive biomarker for dilated cardiolaminopathies with parallel loss of protein expression in cardiomyocyte nuclei.

Determination of ligand pathways in globins: apolar tunnels versus polar gates

Although molecular dynamics simulations suggest multiple interior pathways for O(2) entry into and exit from globins, most experiments indicate well defined single pathways. In 2001, we highlighted the effects of large-to-small amino acid replacements on rates for ligand entry and exit onto the three-dimensional structure of sperm whale myoglobin. The resultant map argued strongly for ligand movement through a short channel from the heme iron to solvent that is gated by the distal histidine (His-64(E7)) near the solvent edge of the porphyrin ring. In this work, we have applied the same mutagenesis mapping strategy to the neuronal mini-hemoglobin from Cerebratulus lacteus (CerHb), which has a large internal tunnel from the heme iron to the C-terminal ends of the E and H helices, a direction that is 180° opposite to the E7 channel. Detailed comparisons of the new CerHb map with expanded results for Mb show unambiguously that the dominant (>90%) ligand pathway in CerHb is through the internal tunnel, and the major (>75%) ligand pathway in Mb is through the E7 gate. These results demonstrate that: 1) mutagenesis mapping can identify internal pathways when they exist; 2) molecular dynamics simulations need to be refined to address discrepancies with experimental observations; and 3) alternative pathways have evolved in globins to meet specific physiological demands.

Hereditary systemic amyloidosis due to Asp76Asn variant β2-microglobulin

We describe a kindred with slowly progressive gastrointestinal symptoms and autonomic neuropathy caused by autosomal dominant, hereditary systemic amyloidosis. The amyloid consists of Asp76Asn variant β(2)-microglobulin. Unlike patients with dialysis-related amyloidosis caused by sustained high plasma concentrations of wild-type β(2)-microglobulin, the affected members of this kindred had normal renal function and normal circulating β(2)-microglobulin values. The Asp76Asn β(2)-microglobulin variant was thermodynamically unstable and remarkably fibrillogenic in vitro under physiological conditions. Previous studies of β(2)-microglobulin aggregation have not shown such amyloidogenicity for single-residue substitutions. Comprehensive biophysical characterization of the β(2)-microglobulin variant, including its 1.40-Å, three-dimensional structure, should allow further elucidation of fibrillogenesis and protein misfolding.

DIO2 modifies inflammatory responses in chondrocytes

OBJECTIVE

Selenium neutralizes interleukin-1β (IL-1β) induced inflammatory responses in chondrocytes. We investigated potential mechanisms for this through in vitro knock down of three major selenoproteins, Iodothyronine Deiodinase-2 (DIO2), Glutathione Peroxidase-1 (GPX1), and Thioredoxin Reductase-1 (TR1) in primary human chondrocytes.

METHODS

Primary human chondrocytes were transfected with scrambled small interfering ribonucleic acid (siRNA) or siRNA specific for DIO2, GPX1 and TR1. After 48 h, transfected cells were cultured in serum free media for 48 h, with or without 10 pg/ml IL-1β for the final 24h. The efficiency of siRNAs was confirmed by quantitative Real Time-Polymerase Chain Reaction (qRT-PCR) and Western blot analysis. The gene expression, by qRT-PCR, of cyclooxygenase-2 (COX2), IL-1β, and Liver X receptor (LXR) alpha and beta was evaluated to determine the impact of selenoprotein knockdown on inflammatory responses in chondrocytes.

RESULTS

The messenger RNA (mRNA) expression of DIO2, GPX1, and TR1 was significantly decreased by the specific siRNAs (reduced 56%, P=0.0004; 96%, P<0.0001; and 66%, P<0.0001, respectively). Suppression of DIO2, but not GPX1 or TR1, significantly increased (~2-fold) both basal (P=0.0005) and IL-1β induced (P<0.0001) COX2 gene expression. Similarly, suppression of DIO2 significantly increased (∼9-fold) IL-1β induced IL-1β gene expression (P=0.0056) and resulted in a 32% (P=0.0044) decrease in LXRα gene expression but no effect on LXRβ.

CONCLUSIONS

Suppression of the selenoprotein DIO2 resulted in strong pro-inflammatory effects with increased expression of inflammatory mediators, IL-1β and COX2, and decreased expression of LXRα suggesting that this may be the upstream target through which the anti-inflammatory effects of DIO2 are mediated.

Ivermectin is a potent inhibitor of flavivirus replication specifically targeting NS3 helicase activity: new prospects for an old drug

OBJECTIVES

Infection with yellow fever virus (YFV), the prototypic mosquito-borne flavivirus, causes severe febrile disease with haemorrhage, multi-organ failure and a high mortality. Moreover, in recent years the Flavivirus genus has gained further attention due to re-emergence and increasing incidence of West Nile, dengue and Japanese encephalitis viruses. Potent and safe antivirals are urgently needed.

METHODS

Starting from the crystal structure of the NS3 helicase from Kunjin virus (an Australian variant of West Nile virus), we identified a novel, unexploited protein site that might be involved in the helicase catalytic cycle and could thus in principle be targeted for enzyme inhibition. In silico docking of a library of small molecules allowed us to identify a few selected compounds with high predicted affinity for the new site. Their activity against helicases from several flaviviruses was confirmed in in vitro helicase/enzymatic assays. The effect on the in vitro replication of flaviviruses was then evaluated.

RESULTS

Ivermectin, a broadly used anti-helminthic drug, proved to be a highly potent inhibitor of YFV replication (EC₅₀ values in the sub-nanomolar range). Moreover, ivermectin inhibited, although less efficiently, the replication of several other flaviviruses, i.e. dengue fever, Japanese encephalitis and tick-borne encephalitis viruses. Ivermectin exerts its effect at a timepoint that coincides with the onset of intracellular viral RNA synthesis, as expected for a molecule that specifically targets the viral helicase.

CONCLUSIONS

The well-tolerated drug ivermectin may hold great potential for treatment of YFV infections. Furthermore, structure-based optimization may result in analogues exerting potent activity against flaviviruses other than YFV.

Structure, stability, and aggregation of β-2 microglobulin mutants: insights from a Fourier transform infrared study in solution and in the crystalline state

β-2 microglobulin (β2m) is an amyloidogenic protein involved in dialysis-related amyloidosis. We report here the study of the structural properties of the protein in solution and in the form of single crystals by Fourier transform infrared (FTIR) spectroscopy and microspectroscopy. The investigation has been extended to four β2m mutants previously characterized by x-ray crystallography: Asp(53)Pro, Asp(59)Pro, Trp(60)Gly, and Trp(60)Val. These variants displayed very similar three-dimensional structures but different thermal stability and aggregation propensity, investigated here by FTIR spectroscopy. For each variant, appreciable spectral differences were found between the protein in solution and in single crystals, consisting in a downshift of the main β-sheet band and in better resolved turn and loop bands, indicative of reduced protein secondary structure dynamics in the crystalline state. Notably, the well-resolved spectra of the β2m crystalline variants enabled us to identify structural differences induced by the single amino acid mutations. Such differences encompass turn and loop structures that might affect the stability and aggregation propensity of the investigated β2m variants. This study highlights the potential of FTIR microspectroscopy to acquire useful structural information on protein crystals, complementary to the crystallographic analyses.

Ligation tunes protein reactivity in an ancient haemoglobin: kinetic evidence for an allosteric mechanism in Methanosarcina acetivorans protoglobin

Protoglobin from Methanosarcina acetivorans (MaPgb) is a dimeric globin with peculiar structural properties such as a completely buried haem and two orthogonal tunnels connecting the distal cavity to the solvent. CO binding to and dissociation from MaPgb occur through a biphasic kinetics. We show that the heterogenous kinetics arises from binding to (and dissociation from) two tertiary conformations in ligation-dependent equilibrium. Ligation favours the species with high binding rate (and low dissociation rate). The equilibrium is shifted towards the species with low binding (and high dissociation) rates for the unliganded molecules. A quantitative model is proposed to describe the observed carbonylation kinetics.

The tempered polymerization of human neuroserpin

Neuroserpin, a member of the serpin protein superfamily, is an inhibitor of proteolytic activity that is involved in pathologies such as ischemia, Alzheimer's disease, and Familial Encephalopathy with Neuroserpin Inclusion Bodies (FENIB). The latter belongs to a class of conformational diseases, known as serpinopathies, which are related to the aberrant polymerization of serpin mutants. Neuroserpin is known to polymerize, even in its wild type form, under thermal stress. Here, we study the mechanism of neuroserpin polymerization over a wide range of temperatures by different techniques. Our experiments show how the onset of polymerization is dependent on the formation of an intermediate monomeric conformer, which then associates with a native monomer to yield a dimeric species. After the formation of small polymers, the aggregation proceeds via monomer addition as well as polymer-polymer association. No further secondary mechanism takes place up to very high temperatures, thus resulting in the formation of neuroserpin linear polymeric chains. Most interesting, the overall aggregation is tuned by the co-occurrence of monomer inactivation (i.e. the formation of latent neuroserpin) and by a mechanism of fragmentation. The polymerization kinetics exhibit a unique modulation of the average mass and size of polymers, which might suggest synchronization among the different processes involved. Thus, fragmentation would control and temper the aggregation process, instead of enhancing it, as typically observed (e.g.) for amyloid fibrillation.

On the molecular structure of human neuroserpin polymers

The polymerization of serpins is at the root of a large class of diseases; the molecular structure of serpin polymers has been recently debated. In this work, we study the polymerization kinetics of human neuroserpin by Fourier Transform Infra Red spectroscopy and by time-lapse Size Exclusion Chromatography. First, we show that two distinct neuroserpin polymers, formed at 45 and 85°C, display the same isosbestic points in the Amide I' band, and therefore share common secondary structure features. We also find a concentration independent polymerization rate at 45°C suggesting that the polymerization rate-limiting step is the formation of an activated monomeric species. The polymer structures are consistent with a model that predicts the bare insertion of portions of the reactive center loop into the A β-sheet of neighboring serpin molecule, although with different extents at 45 and 85°C.

Tetramerization dynamics of C-terminal domain underlies isoform-specific cAMP gating in hyperpolarization-activated cyclic nucleotide-gated channels

Hyperpolarization-activated cyclic nucleotide-gated (HCN) channels are dually activated by hyperpolarization and binding of cAMP to their cyclic nucleotide binding domain (CNBD). HCN isoforms respond differently to cAMP; binding of cAMP shifts activation of HCN2 and HCN4 by 17 mV but shifts that of HCN1 by only 2-4 mV. To explain the peculiarity of HCN1, we solved the crystal structures and performed a biochemical-biophysical characterization of the C-terminal domain (C-linker plus CNBD) of the three isoforms. Our main finding is that tetramerization of the C-terminal domain of HCN1 occurs at basal cAMP concentrations, whereas those of HCN2 and HCN4 require cAMP saturating levels. Therefore, HCN1 responds less markedly than HCN2 and HCN4 to cAMP increase because its CNBD is already partly tetrameric. This is confirmed by voltage clamp experiments showing that the right-shifted position of V(&frac12;) in HCN1 is correlated with its propensity to tetramerize in vitro. These data underscore that ligand-induced CNBD tetramerization removes tonic inhibition from the pore of HCN channels.

Characterization of β2-microglobulin conformational intermediates associated to different fibrillation conditions

β2-Microglobulin (β2m) is the light chain of the class-I major histocompatibility complex, being also the causing agent of dialysis-related amyloidosis, which results from its accumulation as amyloid material in the skeletal joints. This study describes conformational properties of β2m under two distinct, in vitro amyloidogenic conditions: neutral pH in the presence of 20% 2,2,2-trifluoroethanol (TFE) and acidic pH in the absence of TFE. Species distribution analysis by electrospray ionization-mass spectrometry (ESI-MS) is combined with information obtained by ion mobility-mass spectrometry (IM-MS), fluorescence and circular dichroism (CD) spectroscopy. It is shown that β2m populates quite different conformational ensembles under the two conditions, but both ensembles display a minor fraction of the population in a partially folded state. In spite of similar compactness, these two partially folded forms display different conformations: helical secondary structure is predominant in the species at pH 7.4, 20% TFE, while the low-pH form is mainly random coil. As temperature is increased, the TFE intermediate looses helical structure becoming more similar to the low-pH intermediate. The existence of different conformational ensembles may rationalize the different aggregation propensity displayed by β2m under the two fibrillation conditions analyzed here.

The two tryptophans of β2-microglobulin have distinct roles in function and folding and might represent two independent responses to evolutionary pressure

Background

We have recently discovered that the two tryptophans of human β2-microglobulin have distinctive roles within the structure and function of the protein. Deeply buried in the core, Trp95 is essential for folding stability, whereas Trp60, which is solvent-exposed, plays a crucial role in promoting the binding of β2-microglobulin to the heavy chain of the class I major histocompatibility complex (MHCI). We have previously shown that the thermodynamic disadvantage of having Trp60 exposed on the surface is counter-balanced by the perfect fit between it and a cavity within the MHCI heavy chain that contributes significantly to the functional stabilization of the MHCI. Therefore, based on the peculiar differences of the two tryptophans, we have analysed the evolution of β2-microglobulin with respect to these residues.

Results

Having defined the β2-microglobulin protein family, we performed multiple sequence alignments and analysed the residue conservation in homologous proteins to generate a phylogenetic tree. Our results indicate that Trp60 is highly conserved, whereas some species have a Leu in position 95; the replacement of Trp95 with Leu destabilizes β2-microglobulin by 1 kcal/mol and accelerates the kinetics of unfolding. Both thermodynamic and kinetic data fit with the crystallographic structure of the Trp95Leu variant, which shows how the hydrophobic cavity of the wild-type protein is completely occupied by Trp95, but is only half filled by Leu95.

Conclusions

We have established that the functional Trp60 has been present within the sequence of β2-microglobulin since the evolutionary appearance of proteins responsible for acquired immunity, whereas the structural Trp95 was selected and stabilized, most likely, for its capacity to fully occupy an internal cavity of the protein thereby creating a better stabilization of its folded state.

D-strand perturbation and amyloid propensity in beta-2 microglobulin

Proteins hosting main β-sheets adopt specific strategies to avoid intermolecular interactions leading to aggregation and amyloid deposition. Human beta-2 microglobulin (β2m) displays a typical immunoglobulin fold and is known to be amyloidogenic in vivo. Upon severe kidney deficiency, β2m accumulates in the bloodstream, triggering, over the years, pathological deposition of large amyloid aggregates in joints and bones. A β-bulge observed on the edge D β-strand of some β2m crystal structures has been suggested to be crucial in protecting the protein from amyloid aggregation. Conversely, a straight D-strand, observed in different crystal structures of monomeric β2m, could promote amyloid aggregation. More recently, the different conformations observed for the β2m D-strand have been interpreted as the result of intrinsic flexibility, rather than being assigned to a functional protective role against aggregation. To shed light on such contrasting picture, the mutation Asp53→Pro was engineered in β2m, aiming to impair the formation of a regular/straight D-strand. Such a mutant was characterized structurally and biophysically by CD, X-ray crystallography and MS, in addition to an assessment of its amyloid aggregation trends in vitro. The results reported in the present study highlight the conformational plasticity of the edge D-strand, and show that even perturbing the D-strand structure through a Pro residue has only marginal effects on protecting β2m from amyloid aggregation in vitro.