Voltage-Gated Proton Channels in the Tree of Life.

With a single gene encoding HV1 channel, proton channel diversity is particularly low in mammals compared to other members of the superfamily of voltage-gated ion channels. Nonetheless, mammalian HV1 channels are expressed in many different tissues and cell types where they exert various functions. In the first part of this review, we regard novel aspects of the functional expression of HV1 channels in mammals by differentially comparing their involvement in (1) close conjunction with the NADPH oxidase complex responsible for the respiratory burst of phagocytes, and (2) in respiratory burst independent functions such as pH homeostasis or acid extrusion. In the second part, we dissect expression of HV channels within the eukaryotic tree of life, revealing the immense diversity of the channel in other phylae, such as mollusks or dinoflagellates, where several genes encoding HV channels can be found within a single species. In the last part, a comprehensive overview of the biophysical properties of a set of twenty different HV channels characterized electrophysiologically, from Mammalia to unicellular protists, is given.

Proton channels in molluscs: A new bivalvian-specific minimal HV 4 channel.

Recently, three proton channels (HV ) have been identified and characterized in Aplysia californica (AcHV 1-3). Focusing on AcHV 1 and AcHV 2, analysis of Transcriptome Shotgun Assembly and genomic databases of 91 molluscs identified HV homologous channels in other molluscs: channels homologous to AcHV 1 and to AcHV 2 were found in 90 species (56 full-length sequences) and in 33 species (18 full-length sequences), respectively. Here, we report the discovery of a fourth distinct proton channel family, HV 4. This new family has high homology to AcHV 1 and AcHV 2 and was identified only in bivalvian molluscs (13 species, 12 full-length sequences). Typically, these channels possess an extracellular S1-S2 loop of intermediate size (~ 20 amino acids) compared to the shorter loops of molluscan HV 1 channels (~ 13 amino acids) and the much larger loops of molluscan HV 2 channels (> 65 amino acids). The characteristic voltage-sensor motif in S4 possesses only two arginine residues with the common third arginine being replaced by a lysine. Moreover, HV 4 channels are much smaller with only around 200 amino acids in total length. The smallest functional channel found so far in nature (189 amino acids) is expressed in the pacific oyster Crassostrea gigas (CgHV 4) and might be considered an archetypical minimal proton channel. Functional expression and electrophysiological characterization demonstrated that CgHV 4 shares distinctive hallmarks of other investigated proton channels as high proton selectivity, slow activation, and pH- and voltage-regulated gating. This work is the first description of a HV 4 type channel, adding a new member to the recently expanded family of proton channels.

Unexpected expansion of the voltage-gated proton channel family.

Voltage-gated ion channels, whose first identified function was to generate action potentials, are divided into subfamilies with numerous members. The family of voltage-gated proton channels (HV ) is tiny. To date, all species found to express HV have exclusively one gene that codes for this unique ion channel. Here we report the discovery and characterization of three proton channel genes in the classical model system of neural plasticity, Aplysia californica. The three channels (AcHV 1, AcHV 2, and AcHV 3) are distributed throughout the whole animal. Patch-clamp analysis confirmed proton selectivity of these channels but they all differed markedly in gating. AcHV 1 gating resembled HV in mammalian cells where it is responsible for proton extrusion and charge compensation. AcHV 2 activates more negatively and conducts extensive inward proton current, properties likely to acidify the cytosol. AcHV 3, which differs from AcHV 1 and AcHV 2 in lacking the first arginine in the S4 helix, exhibits proton selective leak currents and weak voltage dependence. We report the expansion of the proton channel family, demonstrating for the first time the expression of three functionally distinct proton channels in a single species.

The pH-dependent gating of the human voltage-gated proton channel from computational simulations.

Gating of the voltage-gated proton channel HV1 is strongly controlled by pH. There is evidence that this involves the sidechains of titratable amino acids that change their protonation state with changes of the pH. Despite experimental investigations to identify the amino acids involved in pH sensing only few progress has been made, including one histidine at the cytoplasmic side of the channel that is involved in sensing cellular pH. We have used constant pH molecular dynamics simulations in symmetrical and asymmetrical pH conditions across the membrane to investigate the pH- and ΔpH-dependent gating of the human HV1 channel. Therefore, the pKa of every titratable amino acids has been assessed in single simulations. Our simulations captured initial conformational changes between a deactivated and an activated state of the channel induced solely by changes of the pH. The pH-dependent gating is accompanied by an outward displacement of the three S4 voltage sensing arginines that moves the second arginine past the hydrophobic gasket (HG) which separates the inner and outer pores of the channel. HV1 activation, when outer pH increases, involves amino acids at the extracellular entrance of the channel that extend the network of interactions from the external solution down to the HG. Whereas, amino acids at the cytoplasmic entrance of the channel are involved in activation, when inner pH decreases, and in a network of interactions that extend from the cytoplasm up to the HG.

Voltage-gated proton channels in polyneopteran insects.

Voltage-gated proton channels (HV 1) are expressed in eukaryotes, including basal hexapods and polyneopteran insects. However, currently, there is little known about HV 1 channels in insects. A characteristic aspartate (Asp) that functions as the proton selectivity filter (SF) and the RxWRxxR voltage-sensor motif are conserved structural elements in HV 1 channels. By analysing Transcriptome Shotgun Assembly (TSA) databases, we found 33 polyneopteran species meeting these structural requirements. Unexpectedly, an unusual natural variation Asp to glutamate (Glu) at SF was found in Phasmatodea and Mantophasmatodea. Additionally, we analysed the expression and function of HV 1 in the phasmatodean stick insect Extatosoma tiaratum (Et). EtHV 1 is strongly expressed in nervous tissue and shows pronounced inward proton conduction. This is the first study of a natural occurring Glu within the SF of a functional HV 1 and might be instrumental in uncovering the physiological function of HV 1 in insects.

Assessing Structural Determinants of Zn2+ Binding to Human HV1 via Multiple MD Simulations.

Voltage-gated proton channels (HV1) are essential for various physiological tasks but are strongly inhibited by Zn2+ cations. Some determinants of Zn2+ binding have been elucidated experimentally and in computational studies. However, the results have always been interpreted under the assumption that Zn2+ binds to monomeric HV1 despite evidence that HV1 expresses as a dimer and that the dimer has a higher affinity for zinc than the monomer and experimental data that suggest coordination in the dimer interface. The results of former studies are also controversial, e.g., supporting either one single or two binding sites. Some structural determinants of the binding are still elusive. We performed a series of molecular dynamics simulations to address different structures of the human proton channel, the monomer and two plausible dimer conformations, to compare their respective potential to interact with and bind Zn2+ via the essential histidines. The series consisted of several copies of the system to generate independent trajectories and increase the significance compared to a single simulation. The amount of time simulated totals 29.9 µs for 126 simulations of systems comprising ∼59,000 to ∼187,000 atoms. Our approach confirms the existence of two binding sites in monomeric and dimeric human HV1. The dimer interface is more efficient for attracting and binding Zn2+ via the essential histidines than the monomer or a dimer with the histidines in the periphery. The higher affinity is due to the residues in the dimer interface that create an attractive electrostatic potential funneling the zinc cations toward the binding sites.

Binding properties of SUMO-interacting motifs (SIMs) in yeast.

Small ubiquitin-like modifier (SUMO) conjugation and interaction play an essential role in many cellular processes. A large number of yeast proteins is known to interact non-covalently with SUMO via short SUMO-interacting motifs (SIMs), but the structural details of this interaction are yet poorly characterized. In the present work, sequence analysis of a large dataset of 148 yeast SIMs revealed the existence of a hydrophobic core binding motif and a preference for acidic residues either within or adjacent to the core motif. Thus the sequence properties of yeast SIMs are highly similar to those described for human. Molecular dynamics simulations were performed to investigate the binding preferences for four representative SIM peptides differing in the number and distribution of acidic residues. Furthermore, the relative stability of two previously observed alternative binding orientations (parallel, antiparallel) was assessed. For all SIMs investigated, the antiparallel binding mode remained stable in the simulations and the SIMs were tightly bound via their hydrophobic core residues supplemented by polar interactions of the acidic residues. In contrary, the stability of the parallel binding mode is more dependent on the sequence features of the SIM motif like the number and position of acidic residues or the presence of additional adjacent interaction motifs. This information should be helpful to enhance the prediction of SIMs and their binding properties in different organisms to facilitate the reconstruction of the SUMO interactome.

An information-theoretic classification of amino acids for the assessment of interfaces in protein-protein docking.

Docking represents a versatile and powerful method to predict the geometry of protein-protein complexes. However, despite significant methodical advances, the identification of good docking solutions among a large number of false solutions still remains a difficult task. We have previously demonstrated that the formalism of mutual information (MI) from information theory can be adapted to protein docking, and we have now extended this approach to enhance its robustness and applicability. A large dataset consisting of 22,934 docking decoys derived from 203 different protein-protein complexes was used for an MI-based optimization of reduced amino acid alphabets representing the protein-protein interfaces. This optimization relied on a clustering analysis that allows one to estimate the mutual information of whole amino acid alphabets by considering all structural features simultaneously, rather than by treating them individually. This clustering approach is fast and can be applied in a similar fashion to the generation of reduced alphabets for other biological problems like fold recognition, sequence data mining, or secondary structure prediction. The reduced alphabets derived from the present work were converted into a scoring function for the evaluation of docking solutions, which is available for public use via the web service score-MI: http://score-MI.biochem.uni-erlangen.de.

Nuclear import of isoforms of the cytomegalovirus kinase pUL97 is mediated by differential activity of NLS1 and NLS2 both acting through classical importin-α binding.

The multifunctional protein kinase pUL97 of human cytomegalovirus (HCMV) strongly determines the efficiency of virus replication. Previously, the existence of two pUL97 isoforms that arise from alternative translational initiation and show a predominant nuclear localization was described. Two bipartite nuclear localization sequences, NLS1 and NLS2, were identified in the N terminus of the large isoform, whilst the small isoform exclusively contained NLS2. The current study found the following: (i) pUL97 nuclear localization in HCMV-infected primary fibroblasts showed accumulations in virus replication centres and other nuclear sections; (ii) in a quantitative evaluation system for NLS activity, the large isoform showed higher efficiency of nuclear translocation than the small isoform; (iii) NLS1 was mapped to aa 6-35 and NLS2 to aa 190-213; (iv) using surface plasmon resonance spectroscopy, the binding of both NLS1 and NLS2 to human importin-α was demonstrated, stressing the importance of individual arginine residues in the bipartite consensus motifs; (v) nuclear magnetic resonance spectroscopy of pUL97 peptides confirmed an earlier statement about the functional requirement of NLS1 embedding into an intact α-helical structure; and (vi) a bioinformatics investigation of the solvent-accessible surface suggested a high accessibility of NLS1 and an isoform-specific, variable accessibility of NLS2 for interaction with importin-α. Thus, the nucleocytoplasmic transport mechanism of the isoforms appeared to be differentially regulated, and this may have consequences for isoform-dependent functions of pUL97 during virus replication.

A molecular model for the differential activation of STAT3 and STAT6 by the herpesviral oncoprotein tip.

Constitutive STAT signaling provides growth promoting signals in many forms of malignancy. We performed molecular modeling and molecular dynamics studies of the interaction between the regulatory Src homology 2 (SH2) domains of STAT3 and 6 with phosphorylated peptides of the herpesviral oncoprotein Tip, which facilitates Src kinase mediated STAT-activation and T cell proliferation. The studies give insight into the ligand binding specificity of the STAT SH2 domains and provide the first model for the differential activation of STAT3 or STAT6 by two distinct regions of the viral Tip protein. The biological relevance of the modeled interactions was then confirmed by activation studies using corresponding recombinant oncoproteins, and finally by respective recombinant viruses. The functional data give experimental validation of the molecular dynamics study, and provide evidence for the involvement of STAT6 in the herpesvirus induced T cell proliferation.

Identification of the structural features that mediate binding specificity in the recognition of STAT proteins by dual-specificity phosphatases.

Inactivation of signal transducers and activators of transcription (STAT) proteins is regulated by dual-specificity phosphatases (DSPs) with high substrate specificity. Although experiments have provided useful information about the phosphatase activity and the specificity for STATs, there is up-to-date no data at a molecular level to explain the specific recognition of STAT substrates by this subfamily of phosphatases. Here, a combined approach of molecular modeling, docking and molecular dynamics simulations was used to address the binding between DSPs and their STAT substrates. We identified a binding interface at the protein tyrosine phosphatase (PTP) domain of the DSP VHR that interacts with the SH2-domain of STAT5. This finding is consistent with previous mutational data and supports a "two-step" mechanism for the dephosphorylation event. Application of the same approach suggests the presence of a similar interface between the viral DSP VH1 and STAT1. Furthermore, the interaction network at this interface provides an explanation for the specificity of the DSP-STAT recognition.

Two isoforms of the protein kinase pUL97 of human cytomegalovirus are differentially regulated in their nuclear translocation.

The pUL97 protein kinase encoded by human cytomegalovirus is a multifunctional determinant of the efficiency of viral replication and phosphorylates viral as well as cellular substrate proteins. Here, we report that pUL97 is expressed in two isoforms with molecular masses of approximately 90 and 100 kDa. ORF UL97 comprises an unusual coding strategy in that five in-frame ATG start codons are contained within the N-terminal 157 aa. Site-directed mutagenesis, transient expression of point and deletion mutants and proteomic analyses accumulated evidence that the formation of the large and small isoforms result from alternative initiation of translation, with the start points being at amino acids 1 and 74, respectively. In vitro kinase assays demonstrated that catalytic activity, in terms of autophosphorylation and histone substrate phosphorylation, was indistinguishable for the two isoforms. An analysis of the intracellular distribution of pUL97 by confocal laser-scanning microscopy demonstrated that both isoforms have a pronounced nuclear localization. Surprisingly, mapping experiments performed to identify the nuclear localization signal (NLS) of pUL97 strongly suggest that the mechanism of nuclear transport is distinct for the two isoforms. While the extreme N terminus (large isoform) comprises a highly efficient, bipartite NLS (amino acids 6-35), a second sequence apparently conferring a less efficient mode of nuclear translocation was identified downstream of amino acid 74 (small and large isoforms). Taken together, the findings argue for a complex mechanism of nuclear translocation for pUL97 which might be linked with fine-regulatory differences between the two isoforms.

DNA binding by Corynebacterium glutamicum TetR-type transcription regulator AmtR.

BACKGROUND: The TetR family member AmtR is the central regulator of nitrogen starvation response in Corynebacterium glutamicum. While the AmtR regulon was physiologically characterized in great detail up to now, mechanistic questions of AmtR binding were not addressed. This study presents a characterization of functionally important amino acids in the DNA binding domain of AmtR and of crucial nucleotides in the AmtR recognition motif. RESULTS: Site-directed mutagenesis, the characterization of corresponding mutant proteins by gel retardation assays and surface plasmon resonance and molecular modelling revealed several amino acids, which are directly involved in DNA binding, while others have more structural function. Furthermore, we could show that the spacing of the binding motif half sites is crucial for repression of transcription by AmtR. CONCLUSION: Although the DNA binding domain of TetR-type repressors is highly conserved and a core binding motif was identified for AmtR and TetR(D), the AmtR binding domain shows individual properties compared to other TetR proteins. Besides by distinct amino acids of AmtR, DNA binding is influenced by nucleotides not only of the conserved binding motif but also by spacing nucleotides in C. glutamicum.

Insight into the phosphoryl transfer of the Escherichia coli glucose phosphotransferase system from QM/MM simulations.

Phosphoryl transfer is a key reaction in many aspects of metabolism, gene regulation, and signal transduction. One prominent example is the phosphoenolpyruvate:sugar phosphotransferase system (PTS), which represents an integral part of the bacterial sugar metabolism. The transfer between the enzymes IIA (Glc) and IIB (Glc) in the glucose-specific branch of the PTS is of particular interest due to the unusual combination of donor and acceptor residues involved in phosphoryl transfer: The phosphoryl group is initially attached to the Nepsilon2 atom of His 90 in IIA (Glc) and then transferred to the S gamma atom of Cys 35 in IIB (Glc). To gain insight into the details of the transfer mechanism, we have performed a QM/MM simulation which treats the entire active site quantum-mechanically. The transfer has a high dissociative character, and the Nepsilon2-P bond gets immediately destabilized after complex formation by numerous interactions formed between residues of IIB (Glc) and the phosphoryl group. The final formation of a tight S gamma-P bond is accompanied by a reorientation of the side chain of the phosphoryl donor. This reorientation results in the loss of interaction between the imidazole ring and the phosphate group thus hindering the reverse transfer. A comparison to the transfer in protein tyrosine phosphatases, which also use a cysteine as acceptor of the phosphoryl group, reveals significant similarities in the conformation of the active site, the energy profile of the reaction, and in the pattern of interactions that stabilize the phosphoryl group during the transfer.

Synthesis, crystal structures, linear and nonlinear optical properties, and theoretical studies of (p-R-phenyl)-, (p-R-phenylethynyl)-, and (E)-[2-(p-R-phenyl)ethenyl]dimesitylboranes and related compounds.

The (p-R-phenyl)dimesitylboranes (R=Me(2)N, MeO, MeS, Br, I), (p-R-phenylethynyl)dimesitylboranes (R=Me(2)N, MeO, MeS, H), (E)-[2-(p-R-phenyl)ethenyl]dimesitylboranes (R=Me(2)N, H(2)N, MeO, MeS, H, CN, NO(2)), (E)-[2-(2-thienyl)ethenyl]dimesitylborane, and (E)-[2-(o-carboranyl)ethenyl]dimesitylborane have been prepared through the reaction of the appropriate p-R-phenyl- and p-R-phenylethynyllithium reagents with dimesitylboron fluoride and by hydroboration of the appropriate p-R-phenylacetylene, 2-ethynylthiophene, and o-ethynylcarborane with dimesitylborane. Their UV/Vis absorption and emission spectra have been recorded in a range of solvents with the fluorescence maxima of the donor-substituted compounds in particular exhibiting large bathochromic shifts in highly polar solvents, indicative of charge transfer leading to large dipole moments in the excited state. The molecular structures of the (p-R-phenyl)dimesitylboranes (R=Me(2)N, MeO, MeS, Br, I), the (E)-[2-(p-R-phenyl)ethenyl]dimesitylboranes (R=Me(2)N, H(2)N MeO, MeS, H), (p-R-phenylethynyl)dimesitylborane (R=Me(2)N), and (E)-[2-(2-thienyl)ethenyl]dimesitylborane, which have been determined from single-crystal X-ray diffraction measurements, offer evidence of increased conjugation in the ground state with increased donor strength of the R substituent. Their first- and second-order molecular hyperpolarizabilities have been obtained from EFISH and THG measurements, the first-order hyperpolarizabilities being largest for the strongest R-substituent donors. AM1 calculations have been performed on these compounds, showing reasonable agreement with the experimentally obtained bond lengths and hyperpolarizabilities, as well as on several related hypothetical compounds containing multiple C==C bonds, most of which are proposed to have even larger hyperpolarizabilities.