Identification of a α‐helical molten globule intermediate and structural characterization of β‐cardiotoxin,an all β‐sheet protein isolated from the venom of Ophiophagus hannah (king cobra)

β‐Cardiotoxin is a novel member of the snake venom three‐finger toxin (3FTX) family. This is the first exogenous protein to antagonize β‐adrenergic receptors and thereby causing reduction in heart rates (bradycardia) when administered into animals, unlike the conventional cardiotoxins as reported earlier. 3FTXs are stable all β‐sheet peptides with 60–80 amino acid residues. Here, we describe the three‐dimensional crystal structure of β‐cardiotoxin together with the identification of a molten globule intermediate in the unfolding pathway of this protein. In spite of the overall structural similarity of this protein with conventional cardiotoxins, there are notable differences observed at the loop region and in the charge distribution on the surface, which are known to be critical for cytolytic activity of cardiotoxins. The molten globule intermediate state present in the thermal unfolding pathway of β‐cardiotoxin was however not observed during the chemical denaturation of the protein. Interestingly, circular dichroism (CD) and NMR studies revealed the presence of α‐helical secondary structure in the molten globule intermediate. These results point to substantial conformational plasticity of β‐cardiotoxin, which might aid the protein in responding to the sometimes conflicting demands of structure, stability, and function during its biological lifetime.     

β‐strand interactions at the domain interface critical for the stability of human lens γD‐crystallin

Human age‐onset cataracts are believed to be caused by the aggregation of partially unfolded or covalently damaged lens crystallin proteins; however, the exact molecular mechanism remains largely unknown. We have used microseconds of molecular dynamics simulations with explicit solvent to investigate the unfolding process of human lens γD‐crystallin protein and its isolated domains. A partially unfolded folding intermediate of γD‐crystallin is detected in simulations with its C‐terminal domain (C‐td) folded and N‐terminal domain (N‐td) unstructured, in excellent agreement with biochemical experiments. Our simulations strongly indicate that the stability and the folding mechanism of the N‐td are regulated by the interdomain interactions, consistent with experimental observations. A hydrophobic folding core was identified within the C‐td that is comprised of a and b strands from the Greek key motif 4, the one near the domain interface. Detailed analyses reveal a surprising non‐native surface salt‐bridge between Glu135 and Arg142 located at the end of the ab folded hairpin turn playing a critical role in stabilizing the folding core. On the other hand, an in silico single E135A substitution that disrupts this non‐native Glu135‐Arg142 salt‐bridge causes significant destabilization to the folding core of the isolated C‐td, which, in turn, induces unfolding of the N‐td interface. These findings indicate that certain highly conserved charged residues, that is, Glu135 and Arg142, of γD‐crystallin are crucial for stabilizing its hydrophobic domain interface in native conformation, and disruption of charges on the γD‐crystallin surface might lead to unfolding and subsequent aggregation.     

Cholesterol‐β1AR interaction versus cholesterol‐β2AR interaction

Two 8‐µs all‐atom molecular dynamics simulations have been performed on the two highly homologous G protein‐coupled receptor (GPCR) subtypes, β₁‐ and β₂‐adrenergic receptors, which were embedded in a lipid bilayer with randomly dispersed cholesterol molecules. During the simulations, cholesterol molecules accumulate to different surface regions of the two receptors, suggesting the subtype specificity of cholesterol–β‐adrenergic receptor interaction and providing some clues to the physiological difference of the two subtypes. Meanwhile, comparison between the two receptors in interacting with cholesterols shed some new light on general determinants of cholesterol binding to GPCRs. Our results indicate that although the concave surface, charged residues and aromatic residues are important, neither of these stabilizing factors is indispensable for a cholesterol interaction site. Different combinations of these factors lead to the diversified binding modes of cholesterol binding to the receptors. Our long‐time simulations, for the first time, revealed the pathway of a cholesterol molecule entering the consensus cholesterol motif (CCM) site, and the binding process of cholesterol to CCM is accompanied by a side chain flipping of the conserved Trp4.50. Moreover, the simulation results suggest that the I‐/V‐/L‐rich region on the extracellular parts of helix 6 might be an alternatively conserved cholesterol‐binding site for the class‐A GPCRs. Proteins 2014; 82:760–770. © 2013 Wiley Periodicals, Inc.     

Membrane dynamics of γ‐secretase with the anterior pharynx‐defective 1B subunit

The four‐subunit protease complex γ‐secretase cleaves many single‐pass transmembrane (TM) substrates, including Notch and β‐amyloid precursor protein to generate amyloid‐β (Aβ), central to Alzheimer's disease. Two of the subunits anterior pharynx‐defective 1 (APH‐1) and presenilin (PS) exist in two homologous forms APH1‐A and APH1‐B, and PS1 and PS2. The consequences of these variations are poorly understood and could affect Aβ production and γ‐secretase medicine. Here, we developed the first complete structural model of the APH‐1B subunit using the published cryo‐electron microscopy (cryo‐EM) structures of APH1‐A (Protein Data Bank: 5FN2, 5A63, and 6IYC). We then performed all‐atom molecular dynamics simulations at 303 K in a realistic bilayer system to understand both APH‐1B alone and in γ‐secretase without and with substrate C83‐bound. We show that APH‐1B adopts a 7TM topology with a water channel topology similar to APH‐1A. We demonstrate direct transport of water through this channel, mainly via Glu84, Arg87, His170, and His196. The apo and holo states closely resemble the experimental cryo‐EM structures with APH‐1A, however with subtle differences: The substrate‐bound APH‐1B γ‐secretase was quite stable, but some TM helices of PS1 and APH‐1B rearranged in the membrane consistent with the disorder seen in the cryo‐EM data. This produces different accessibility of water molecules for the catalytic aspartates of PS1, critical for Aβ production. In particular, we find that the typical distance between the catalytic aspartates of PS1 and the C83 cleavage sites are shorter in APH‐1B, that is, it represents a more closed state, due to interactions with the C‐terminal fragment of PS1. Our structural‐dynamic model of APH‐1B alone and in γ‐secretase suggests generally similar topology but some notable differences in water accessibility which may be relevant to the protein's existence in two forms and their specific function and location.     

Phosphorylation stabilizes the N‐termini of α‐helices

The role of phosphorylation in stabilizing the N‐termini of α‐helices is examined using computer simulations of model peptides. The models comprise either a phosphorylated or unphosphorylated serine at the helix N‐terminus, followed by nine alanines. Monte Carlo/stochastic Dynamics simulations were performed on the model helices. The simulations revealed a distinct stabilization of the helical conformation at the N‐terminus after phosphorylation. The stabilization was attributable to favorable electrostatic interactions between the phosphate and the helix backbone. However, direct helix capping by the phosphorylated sidechain was not observed. The results of the calculations are consistent with experimental evidence on the stabilization of helices by phosphates and other anions. © 1999 John Wiley & Sons, Inc. Biopoly 49: 225–233, 1999     

Protein prosthesis: β‐peptides as reverse‐turn surrogates

The introduction of non‐natural modules could provide unprecedented control over folding/unfolding behavior, conformational stability, and biological function of proteins. Success requires the interrogation of candidate modules in natural contexts. Here, expressed protein ligation is used to replace a reverse turn in bovine pancreatic ribonuclease (RNase A) with a synthetic β‐dipeptide: β²‐homoalanine–β³‐homoalanine. This segment is known to adopt an unnatural reverse‐turn conformation that contains a 10‐membered ring hydrogen bond, but one with a donor–acceptor pattern opposite to that in the 10‐membered rings of natural reverse turns. The RNase A variant has intact enzymatic activity, but unfolds more quickly and has diminished conformational stability relative to native RNase A. These data indicate that hydrogen‐bonding pattern merits careful consideration in the selection of beneficial reverse‐turn surrogates.     

Conserved buried water molecules enable the β‐trefoil architecture

Available high‐resolution crystal structures for the family of β‐trefoil proteins in the structural databank were queried for buried waters. Such waters were classified as either: (a) unique to a particular domain, family, or superfamily or (b) conserved among all β‐trefoil folds. Three buried waters conserved among all β‐trefoil folds were identified. These waters are related by the threefold rotational pseudosymmetry characteristic of this protein architecture (representing three instances of an identical structural environment within each repeating trefoil‐fold motif). The structural properties of this buried water are remarkable and include: residing in a cavity space no larger than a single water molecule, exhibiting a positional uncertainty (i.e., normalized B‐factor) substantially lower than the average Cα atom, providing essentially ideal H‐bonding geometry with three solvent‐inaccessible main chain groups, simultaneously serving as a bridging H‐bond for three different β‐strands at a point of secondary structure divergence, and orienting conserved hydrophobic side chains to form a nascent core‐packing group. Other published work supports an interpretation that these interactions are key to the formation of an efficient folding nucleus and folded thermostability. The fundamental threefold symmetric structural element of the β‐trefoil fold is therefore, surprisingly, a buried water molecule.     

5α‐Androst‐16‐en‐3α‐ol β‐D‐Glucuronide,Precursor of 5α‐Androst‐16‐en‐3α‐ol in Human Sweat

5α‐Androst‐16‐en‐3α‐ol (α‐androstenol) is an important contributor to human axilla sweat odor. It is assumed that α‐andostenol is excreted from the apocrine glands via a H₂O‐soluble conjugate, and this precursor was formally characterized in this study for the first time in human sweat. The possible H₂O‐soluble precursors, sulfate and glucuronide derivatives, were synthesized as analytical standards, i.e., α‐androstenol, β‐androstenol sulfates, 5α‐androsta‐5,16‐dien‐3β‐ol (β‐androstadienol) sulfate, α‐androstenol β‐glucuronide, α‐androstenol α‐glucuronide, β‐androstadienol β‐glucuronide, and α‐androstenol β‐glucuronide furanose. The occurrence of α‐androstenol β‐glucuronide was established by ultra performance liquid chromatography (UPLC)/MS (heated electrospray ionization (HESI)) in negative‐ion mode in pooled human sweat, containing eccrine and apocrine secretions and collected from 25 female and 24 male underarms. Its concentration was of 79 ng/ml in female secretions and 241 ng/ml in male secretions. The release of α‐androstenol was observed after incubation of the sterile human sweat or α‐androstenol β‐glucuronide with a commercial glucuronidase enzyme, the urine‐isolated bacteria Streptococcus agalactiae, and the skin bacteria Staphylococcus warneri DSM 20316, Staphylococcus haemolyticus DSM 20263, and Propionibacterium acnes ATCC 6919, reported to have β‐glucuronidase activities. We demonstrated that if α‐ and β‐androstenols and androstadienol sulfates were present in human sweat, their concentrations would be too low to be considered as potential precursors of malodors; therefore, the H₂O‐soluble precursor of α‐androstenol in apocrine secretion should be a β‐glucuronide.     

Structural characterization of a neurotoxic threonine‐rich peptide corresponding to the human prion protein α2‐helical 180–195 segment,and comparison with full‐length α2‐helix‐derived peptides

The 173–195 segment corresponding to the helix 2 of the globular PrP domain is a good candidate to be one of the several ‘spots’ of intrinsic structural flexibility, which might induce local destabilization and concur to protein transformation, leading to aggregation‐prone conformations. Here, we report CD and NMR studies on the α2‐helix‐derived peptide of maximal length (hPrP[180–195]) that is able to exhibit a regular structure different from the prevalently random arrangement of other α2‐helix‐derived peptides. This peptide, which has previously been shown to be affected by buffer composition via the ion charge density dependence typical of Hofmeister effects, corresponds to the C‐terminal sequence of the PrP^C full‐length α2‐helix and includes the highly conserved threonine‐rich 188–195 segment. At neutral pH, its conformation is dominated by β‐type contributions, which only very strong environmental modifications are able to modify. On TFE addition, an increase of α‐helical content can be observed, but a fully helical conformation is only obtained in neat TFE. However, linking of the 173–179 segment, as occurring in wild‐type and mutant peptides corresponding to the full‐length α2‐helix, perturbs these intrinsic structural propensities in a manner that depends on whether the environment is water or TFE. Overall, these results confirm that the 180–195 parental region in hPrP^C makes a strong contribution to the chameleon conformational behavior of the segment corresponding to the full‐length α2‐helix, and could play a role in determining structural rearrangements of the entire globular domain. Copyright © 2008 European Peptide Society and John Wiley & Sons, Ltd.     

Tumour necrosis factor‐α promotes liver ischaemia‐reperfusion injury through the PGC‐1α/Mfn2 pathway

Tumour necrosis factor (TNF)‐α has been considered to induce ischaemia‐reperfusion injury (IRI) of liver which is characterized by energy dysmetabolism. Peroxisome proliferator–activated receptor‐γ co‐activator (PGC)‐1α and mitofusion2 (Mfn2) are reported to be involved in the regulation of mitochondrial function. However, whether PGC‐1α and Mfn2 form a pathway that mediates liver IRI, and if so, what the underlying involvement is in that pathway remain unclear. In this study, L02 cells administered recombinant human TNF‐α had increased TNF‐α levels and resulted in down‐regulation of PGC‐1α and Mfn2 in a rat liver IRI model. This was associated with hepatic mitochondrial swelling, decreased adenosine triphosphate (ATP) production, and increased levels of reactive oxygen species (ROS) and alanine aminotransferase (ALT) activity as well as cell apoptosis. Inhibition of TNF‐α by neutralizing antibody reversed PGC‐1α and Mfn2 expression, and decreased hepatic injury and cell apoptosis both in cell culture and in animals. Treatment by rosiglitazone sustained PGC‐1α and Mfn2 expression both in IR livers, and L02 cells treated with TNF‐α as indicated by increased hepatic mitochondrial integrity and ATP production, reduced ROS and ALT activity as well as decreased cell apoptosis. Overexpression of Mfn2 by lentiviral‐Mfn2 transfection decreased hepatic injury in IR livers and L02 cells treated with TNF‐α. However, there was no up‐regulation of PGC‐1α. These findings suggest that PGC‐1α and Mfn2 constitute a regulatory pathway, and play a critical role in TNF‐α‐induced hepatic IRI. Inhibition of the TNF‐α or PGC‐1α/Mfn2 pathways may represent novel therapeutic interventions for hepatic IRI.     

Mitotic spindle association of TACC3 requires Aurora‐A‐dependent stabilization of a cryptic α‐helix

Aurora‐A regulates the recruitment of TACC3 to the mitotic spindle through a phospho‐dependent interaction with clathrin heavy chain (CHC). Here, we describe the structural basis of these interactions, mediated by three motifs in a disordered region of TACC3. A hydrophobic docking motif binds to a previously uncharacterized pocket on Aurora‐A that is blocked in most kinases. Abrogation of the docking motif causes a delay in late mitosis, consistent with the cellular distribution of Aurora‐A complexes. Phosphorylation of Ser558 engages a conformational switch in a second motif from a disordered state, needed to bind the kinase active site, into a helical conformation. The helix extends into a third, adjacent motif that is recognized by a helical‐repeat region of CHC, not a recognized phospho‐reader domain. This potentially widespread mechanism of phospho‐recognition provides greater flexibility to tune the molecular details of the interaction than canonical recognition motifs that are dominated by phosphate binding.     

α4 phosphoprotein interacts with EDD E3 ubiquitin ligase and poly(A)‐binding protein

Mammalian α4 phosphoprotein, the homolog of yeast Tap42, is a component of the mammalian target‐of‐rapamycin (mTOR) pathway that regulates ribogenesis, the initiation of translation, and cell‐cycle progression. α4 is known to interact with the catalytic subunit of protein phosphatase 2A (PP2Ac) and to regulate PP2A activity. Using α4 as bait in yeast two‐hybrid screening of a human K562 erythroleukemia cDNA library, EDD (E3 isolated by differential display) E3 ubiquitin ligase was identified as a new protein partner of α4. EDD is the mammalian ortholog of Drosophila hyperplastic discs gene (hyd) that controls cell proliferation during development. The EDD protein contains a PABC domain that is present in poly(A)‐binding protein (PABP), suggesting that PABP may also interact with α4. PABP recruits translation factors to the poly(A)‐tails of mRNAs. In the present study, immunoprecipitation/immunoblotting (IP/IB) analyses showed a physical interaction between α4 and EDD in rat Nb2 T‐lymphoma and human MCF‐7 breast cancer cell lines. α4 also interacted with PABP in Nb2, MCF‐7 and the human Jurkat T‐leukemic and K562 myeloma cell lines. COS‐1 cells, transfected with Flag‐tagged‐pSG5‐EDD, gave a (Flag)‐EDD–α4 immunocomplex. Furthermore, deletion mutants of α4 were constructed to determine the binding site for EDD. IP/IB analysis showed that EDD bound to the C‐terminal region of α4, independent of the α4‐PP2Ac binding site. Therefore, in addition to PP2Ac, α4 interacts with EDD and PABP, suggesting its involvement in multiple steps in the mTOR pathway that leads to translation initiation and cell‐cycle progression. J. Cell. Biochem. 110: 1123–1129, 2010. Published 2010 Wiley‐Liss, Inc.     

nNOS–CAPON interaction mediates amyloid‐β‐induced neurotoxicity,especially in the early stages

In neurons, increased protein–protein interactions between neuronal nitric oxide synthase (nNOS) and its carboxy‐terminal PDZ ligand (CAPON) contribute to excitotoxicity and abnormal dendritic spine development, both of which are involved in the development of Alzheimer's disease. In models of Alzheimer's disease, increased nNOS–CAPON interaction was detected after treatment with amyloid‐β in vitro, and a similar change was found in the hippocampus of APP/PS1 mice (a transgenic mouse model of Alzheimer's disease), compared with age‐matched background mice in vivo. After blocking the nNOS–CAPON interaction, memory was rescued in 4‐month‐old APP/PS1 mice, and dendritic impairments were ameliorated both in vivo and in vitro. Furthermore, we demonstrated that S‐nitrosylation of Dexras1 and inhibition of the ERK–CREB–BDNF pathway might be downstream of the nNOS–CAPON interaction.     

The minimal α‐crystallin domain of Mj Hsp16.5 is functional at non‐heat‐shock conditions

The small heat shock protein (sHSP) from Methanococcus jannaschii (Mj Hsp16.5) forms a monodisperse 24mer and each of its monomer contains two flexible N‐ and C‐terminals and a rigid α‐crystallin domain with an extruding β‐strand exchange loop. The minimal α‐crystallin domain with a β‐sandwich fold is conserved in sHSP family, while the presence of the β‐strand exchange loop is divergent. The function of the β‐strand exchange loop and the minimal α‐crystallin domain of Mj Hsp16.5 need further study. In the present study, we constructed two fragment‐deletion mutants of Mj Hsp16.5, one with both the N‐ and C‐terminals deleted (ΔNΔC) and the other with a further deletion of the β‐strand exchange loop (ΔNΔLΔC). ΔNΔC existed as a dimer in solution. In contrast, the minimal α‐crystallin domain ΔNΔLΔC became polydisperse in solution and exhibited more efficient chaperone‐like activities to prevent amorphous aggregation of insulin B chain and fibril formation of the amyloidogenic peptide dansyl‐SSTSAA‐W than the mutant ΔNΔC and the wild type did. The hydrophobic probe binding experiments indicated that ΔNΔLΔC exposed much more hydrophobic surface than ΔNΔC. Our study also demonstrated that Mj Hsp16.5 used different mechanisms for protecting different substrates. Though Mj Hsp16.5 formed stable complexes with substrates when preventing thermal aggregation, no complexes were detected when preventing aggregation under non‐heat‐shock conditions. Proteins 2014; 82:1156–1167. © 2013 Wiley Periodicals, Inc.     

Monte Carlo simulation of 4‐α‐glucanotransferase reaction

4‐α‐Glucanotransferase (GTase, D ‐enzyme) catalyzes disproportionation between two short polymers of maltooligosaccharides linked by α‐1,4‐glucoside bonds. Using action modes of the potato GTase for the donor and acceptor substrates, the Monte Carlo method was applied to simulate the GTase reaction. The simulation starts from a single enzyme molecule and a finite number (10⁵) of substrate molecules. All selection processes were performed using random numbers produced by computer. The initial substrates were from trimer to 10‐mer. In every case, the final stage was the steady‐state distribution of polymers. The steady‐state distribution by the potato GTase reaction was different from those by the hypothetical random disproportionation reaction. The simulated data from the reaction of potato GTase and trimer almost quantitatively agreed with experimental data. The mechanism of the GTase reaction was accumulation of probabilistic processes and was well simulated by the Monte Carlo method. GTase randomizes the overall distribution of chain length of the substrate. Therefore the GTase reaction is an entropy‐driven process. © 1999 John Wiley & Sons, Inc. Biopoly 50: 145–151, 1999     

The β‐sheet breakers and π‐stacking

Fibrillation of β‐amyloid is recognized as a key process leading to the development of Alzheimer's disease. Small peptides called β‐sheet breakers were found to inhibit the process of β‐amyloid fibrillation and to dissolve amyloid fibrils in vitro, in vivo, and in cell culture studies [1,2]. The mechanism by which peptide inhibition takes place remains elusive and a detailed model needs to be established. Here, we present new insights into the possible role of consecutive Phe residues, present in the structure of β‐sheet breakers, supported by the results obtained by means of MD simulations. We performed a 30‐ns MD of two β‐sheet breakers: iAβ5 (LPFFD) and iAβ6 (LPFFFD) which have two and three consecutive Phe residues, respectively. We have found that Phe rings in these peptides tend to form stacked conformations. For one of the peptides – iAβ6 – the calculated electrostatic contribution to free energy of one of the conformers with three rings stacked (c2) is significantly lower than that corresponding to the unstacked one (c1), two rings stacked (c0) and second conformer with three rings stacked (c3). This may favor the interaction of the c2 conformer with the target on amyloid fibril. We hypothesize that the mechanism of inhibition of amyloidogenesis by β‐sheet breaker involves competition among π‐stacked Phe residues of the inhibitor and π‐stacking within the β‐amyloid fibril. iAβ6 may be a promising candidate for a lead compound of amyloidogenesis inhibitors. Copyright © 2013 European Peptide Society and John Wiley & Sons, Ltd.