Zinc finger domains: hypotheses and current knowledge

The glycine-rich domain of GRP7 plays a crucial role in binding long RNAs and facilitating phase separation

Microscale thermophoresis (MST) is a well-established method to quantify protein-RNA interactions. In this study, we employed MST to analyze the RNA binding properties of glycine-rich RNA binding protein 7 (GRP7), which is known to have multiple biological functions related to its ability to bind different types of RNA. However, the exact mechanism of GRP7's RNA binding is not fully understood. While the RNA-recognition motif of GRP7 is known to be involved in RNA binding, the glycine-rich region (known as arginine-glycine-glycine-domain or RGG-domain) also influences this interaction. To investigate to which extend the RGG-domain of GRP7 is involved in RNA binding, mutation studies on putative RNA interacting or modulating sites were performed. In addition to MST experiments, we examined liquid-liquid phase separation of GRP7 and its mutants, both with and without RNA. Furthermore, we systemically investigated factors that might affect RNA binding selectivity of GRP7 by testing RNAs of different sizes, structures, and modifications. Consequently, our study revealed that GRP7 exhibits a high affinity for a variety of RNAs, indicating a lack of pronounced selectivity. Moreover, we established that the RGG-domain plays a crucial role in binding longer RNAs and promoting phase separation.

Decoding the Parkinson's Symphony: PARIS, Maestro of Transcriptional Regulation and Metal Coordination for Dopamine Release

Parkin interacting substrate (PARIS) is a pivotal transcriptional regulator in the brain that orchestrates the activity of various enzymes through its intricate interactions with biomolecules, including nucleic acids. Notably, the binding of PARIS to insulin response sequences (IRSs) triggers a cascade of events that results in the functional loss in the substantia nigra, which impairs dopamine release and, subsequently, exacerbates the relentless neurodegeneration. Here, we report the details of the interactions of PARIS with IRSs via classical zinc finger (ZF) domains in PARIS, namely, PARIS(ZF2-4). Our biophysical studies with purified PARIS(ZF2-4) elucidated the binding partner of PARIS, which generates specific interactions with the IRS1 (5'-TATTTTT, Kd = 38.9 ± 2.4 nM) that is positioned in the promoter region of peroxisome proliferator-activated receptor gamma coactivator-1α (PGC-1α). Mutational and metal-substitution studies demonstrated that Zn(II)-PARIS(ZF2-4) could recognize its binding partner selectively. Overall, our work provides submolecular details regarding PARIS and shows that it is a transcriptional factor that regulates dopamine release. Thus, PARIS could be a crucial target for therapeutic applications.

The past, present, and future of artificial zinc finger proteins: design strategies and chemical and biological applications

Zinc finger proteins are abundant in the human proteome and are responsible for a variety of functions. The domains that constitute zinc finger proteins are compact spherical structures, each comprising approximately 30 amino acid residues, but they also have precise molecular factor functions: zinc binding and DNA recognition. Due to the biological importance of zinc finger proteins and their unique structural and functional properties, many artificial zinc finger proteins have been created and are expected to improve their functions and biological applications. In this study, we review previous studies on the redesign and application of artificial zinc finger proteins, focusing on the experimental results obtained by our research group. In addition, we systematically review various design strategies used to construct artificial zinc finger proteins and discuss in detail their potential biological applications, including gene editing. This review will provide relevant information to researchers involved or interested in the field of artificial zinc finger proteins as a potential new treatment for various diseases.

Quantum chemistry rules retinoid biology

This Perspective discusses how retinol catalyzes resonance energy transfer (RET) reactions pivotally important for mitochondrial energy homeostasis by protein kinase C δ (PKCδ). PKCδ signals to the pyruvate dehydrogenase complex, controlling oxidative phosphorylation. The PKCδ-retinol complex reversibly responds to the redox potential of cytochrome c, that changes with the electron transfer chain workload. In contrast, the natural retinoid anhydroretinol irreversibly activates PKCδ. Its elongated conjugated-double-bond system limits the energy quantum absorbed by RET. Consequently, while capable of triggering the exergonic activating pathway, anhydroretinol fails to activate the endergonic silencing path, trapping PKCδ in the ON position and causing harmful levels of reactive oxygen species. However, physiological retinol levels displace anhydroretinol, buffer cyotoxicity and potentially render anhydroretinol useful for rapid energy generation. Intriguingly, apocarotenoids, the primary products of the mitochondrial β-carotene,9'-10'-oxygenase, have all the anhydroretinol-like features, including modulation of energy homeostasis. We predict significant conceptual advances to stem from further understanding of the retinoid-catalyzed RET.

Alignment-free metal ion-binding site prediction from protein sequence through pretrained language model and multi-task learning

More than one-third of the proteins contain metal ions in the Protein Data Bank. Correct identification of metal ion-binding residues is important for understanding protein functions and designing novel drugs. Due to the small size and high versatility of metal ions, it remains challenging to computationally predict their binding sites from protein sequence. Existing sequence-based methods are of low accuracy due to the lack of structural information, and time-consuming owing to the usage of multi-sequence alignment. Here, we propose LMetalSite, an alignment-free sequence-based predictor for binding sites of the four most frequently seen metal ions in BioLiP (Zn2+, Ca2+, Mg2+ and Mn2+). LMetalSite leverages the pretrained language model to rapidly generate informative sequence representations and employs transformer to capture long-range dependencies. Multi-task learning is adopted to compensate for the scarcity of training data and capture the intrinsic similarities between different metal ions. LMetalSite was shown to surpass state-of-the-art structure-based methods by more than 19.7, 14.4, 36.8 and 12.6% in area under the precision recall on the four independent tests, respectively. Further analyses indicated that the self-attention modules are effective to learn the structural contexts of residues from protein sequence. We provide the data sets, source codes and trained models of LMetalSite at https://github.com/biomed-AI/LMetalSite.

Molecular Dynamics Simulations of Reduced and Oxidized TFIIIA Zinc Fingers Free and Interacting with 5S RNA

Interactions of zinc finger (ZF) proteins with nucleic acids and proteins play an important role in DNA transcription and repair, biochemical recognition, and protein regulation. The release of Zn²⁺ through oxidation of cysteine thiolates is associated with disruption of gene expression and DNA repair, preventing tumor growth. Multi-microsecond molecular dynamics (MD) simulations were carried out to examine the effect of Cys oxidation on the ZF456 fragment of transcription factor III A (TFIIIA) and its complex with 5S RNA. In the absence of 5S RNA, the reduced ZF456 peptide undergoes conformational changes in the secondary structure due to the reorientation of the intact ZF domains. Upon oxidation, the individual ZF domains unfold to various degrees, yielding a globular ZF456 peptide with ZF4 and ZF6, responsible for base-specific hydrogen bonds with 5S RNA, losing their ββα-folds. ZF5, on the other hand, participates in nonspecific interactions through its α-helix that conditionally unravels early in the simulation. In the presence of RNA, oxidation of the ZF456 peptide disrupts the key hydrogen bonding interactions between ZF5/ZF6 and 5S RNA. However, interactions with ZF4 are dependent on the protonation state of His119.

Chemoproteomic Profiling by Cysteine Fluoroalkylation Reveals Myrocin G as an Inhibitor of the Nonhomologous End Joining DNA Repair Pathway

Chemoproteomic profiling of cysteines has emerged as a powerful method for screening the proteome-wide targets of cysteine-reactive fragments, drugs, and natural products. Herein, we report the development and an in-depth evaluation of a tetrafluoroalkyl benziodoxole (TFBX) as a cysteine-selective chemoproteomic probe. We show that this probe features numerous key improvements compared to the traditionally used cysteine-reactive probes, including a superior target occupancy, faster labeling kinetics, and broader proteomic coverage, thus enabling profiling of cysteines directly in live cells. In addition, the fluorine "signature" of probe 7 constitutes an additional advantage resulting in a more confident adduct-amino acid site assignment in mass-spectrometry-based identification workflows. We demonstrate the utility of our new probe for proteome-wide target profiling by identifying the cellular targets of (-)-myrocin G, an antiproliferative fungal natural product with a to-date unknown mechanism of action. We show that this natural product and a simplified analogue target the X-ray repair cross-complementing protein 5 (XRCC5), an ATP-dependent DNA helicase that primes DNA repair machinery for nonhomologous end joining (NHEJ) upon DNA double-strand breaks, making them the first reported inhibitors of this biomedically highly important protein. We further demonstrate that myrocins disrupt the interaction of XRCC5 with DNA leading to sensitization of cancer cells to the chemotherapeutic agent etoposide as well as UV-light-induced DNA damage. Altogether, our next-generation cysteine-reactive probe enables broader and deeper profiling of the cysteinome, rendering it a highly attractive tool for elucidation of targets of electrophilic small molecules.

DeepCys: Structure-based multiple cysteine function prediction method trained on deep neural network: Case study on domains of unknown functions belonging to COX2 domains

Cysteine (Cys) is the most reactive amino acid participating in a wide range of biological functions. In-silico predictions complement the experiments to meet the need of functional characterization. Multiple Cys function prediction algorithm is scarce, in contrast to specific function prediction algorithms. Here we present a deep neural network-based multiple Cys function prediction, available on web-server (DeepCys) (https://deepcys.herokuapp.com/). DeepCys model was trained and tested on two independent datasets curated from protein crystal structures. This prediction method requires three inputs, namely, PDB identifier (ID), chain ID and residue ID for a given Cys and outputs the probabilities of four cysteine functions, namely, disulphide, metal-binding, thioether and sulphenylation and predicts the most probable Cys function. The algorithm exploits the local and global protein properties, like, sequence and secondary structure motifs, buried fractions, microenvironments and protein/enzyme class. DeepCys outperformed most of the multiple and specific Cys function algorithms. This method can predict maximum number of cysteine functions. Moreover, for the first time, explicitly predicts thioether function. This tool was used to elucidate the cysteine functions on domains of unknown functions belonging to cytochrome C oxidase subunit-II like transmembrane domains. Apart from the web-server, a standalone program is also available on GitHub (https://github.com/vam-sin/deepcys).

Two-Dimensional Coordination Networks from Cyclic Dipeptides

Peptide-based biomimetic nanostructures and metal-organic coordination networks on surfaces are two promising classes of hybrid materials which have been explored recently. However, despite the great versatility and structural variability of natural and synthetic peptides, the two directions have so far not been merged in fabrication of metal-organic coordination networks using peptides as building blocks. Here we demonstrate that cyclic peptides can be used as ligands to form highly ordered, two-dimensional, peptide-based metal-organic coordination networks. The networks are formed on a Au(111) surface through coadsorption of cyclic dialanine with Cu-adatoms under Ultra-High Vacuum (UHV) conditions. Scanning Tunneling Microscopy (STM) in combination with X-ray Photoelectron spectroscopy (XPS) has been utilized to characterize the network structures at submolecular resolution and expound the chemical changes involved in network coordination. The networks involve a motif of three cyclic dialanine molecules coordinating to a central Cu-adatom. Interestingly the networks expose pores functionalized by the side chain of the cyclic peptide, suggesting a general method to form functionalized porous metal-organic networks on surfaces.

Global pentapeptide statistics are far away from expected distributions

The relationships between polypeptide composition, sequence, structure and function have been puzzling biologists ever since first protein sequences were determined. Here, we study the statistics of occurrence of all possible pentapeptide sequences in known proteins. To compensate for the non-uniform distribution of individual amino acid residues in protein sequences, we investigate separately all possible permutations of every given amino acid composition. For the majority of permutation groups we find that pentapeptide occurrences deviate strongly from the expected binomial distributions, and that the observed distributions are also characterized by high numbers of outlier sequences. An analysis of identified outliers shows they often contain known motifs and rare amino acids, suggesting that they represent important functional elements. We further compare the pentapeptide composition of regions known to correspond to protein domains with that of non-domain regions. We find that a substantial number of pentapeptides is clearly strongly favored in protein domains. Finally, we show that over-represented pentapeptides are significantly related to known functional motifs and to predicted ancient structural peptides.

In Vivo Genome Editing as a Therapeutic Approach

Genome editing has been well established as a genome engineering tool that enables researchers to establish causal linkages between genetic mutation and biological phenotypes, providing further understanding of the genetic manifestation of many debilitating diseases. More recently, the paradigm of genome editing technologies has evolved to include the correction of mutations that cause diseases via the use of nucleases such as zinc-finger nucleases (ZFN), transcription activator-like effector nucleases (TALENs), and more recently, Cas9 nuclease. With the aim of reversing disease phenotypes, which arise from somatic gene mutations, current research focuses on the clinical translatability of correcting human genetic diseases in vivo, to provide long-term therapeutic benefits and potentially circumvent the limitations of in vivo cell replacement therapy. In this review, in addition to providing an overview of the various genome editing techniques available, we have also summarized several in vivo genome engineering strategies that have successfully demonstrated disease correction via in vivo genome editing. The various benefits and challenges faced in applying in vivo genome editing in humans will also be discussed.

Human PRDM2: Structure, function and pathophysiology

PRDM2/RIZ is a member of a superfamily of histone/protein methyltransferases (PRDMs), which are characterized by the conserved N-terminal PR domain, with methyltransferase activity and zinc finger arrays at the C-terminus. Similar to other family members, two main protein types, known as RIZ1 and RIZ2, are produced from the PRDM2 locus differing by the presence or absence of the PR domain. The imbalance in their respective amounts may be an important cause of malignancy, with the PR-positive isoform commonly lost or downregulated and the PR-negative isoform always being present at higher levels in cancer cells. Interestingly, the RIZ1 isoform also represents an important target of estradiol action downstream of the interaction with hormone receptor. Furthermore, the imbalance between the two products could also be a molecular basis for other human diseases. Thus, understanding the molecular mechanisms underlying PRDM2 function could be useful in the pathophysiological context, with a potential to exploit this information in clinical practice.

Cysteine-reactive probes and their use in chemical proteomics

Proteomic profiling using bioorthogonal chemical probes that selectively react with certain amino acids is now a widely used method in life sciences to investigate enzymatic activities, study posttranslational modifications and discover novel covalent inhibitors. Over the past two decades, researchers have developed selective probes for several different amino acids, including lysine, serine, cysteine, threonine, tyrosine, aspartate and glutamate. Among these amino acids, cysteines are particularly interesting due to their highly diverse and complex biochemical role in our cells. In this feature article, we focus on the chemical probes and methods used to study cysteines in complex proteomes.

Applications of Reactive Cysteine Profiling

Cysteine thiols are involved in a diverse set of biological transformations, including nucleophilic and redox catalysis, metal coordination and formation of both dynamic and structural disulfides. Often posttranslationally modified, cysteines are also frequently alkylated by electrophilic compounds, including electrophilic metabolites, drugs, and natural products, and are attractive sites for covalent probe and drug development. Quantitative proteomics combined with activity-based protein profiling has been applied to annotate cysteine reactivity, susceptibility to posttranslational modifications, and accessibility to chemical probes, uncovering thousands of functional and small-molecule targetable cysteines across a diverse set of proteins, proteome-wide in an unbiased manner. Reactive cysteines have been targeted by high-throughput screening and fragment-based ligand discovery efforts. New cysteine-reactive electrophiles and compound libraries have been synthesized to enable inhibitor discovery broadly and to minimize nonspecific toxicity and off-target activity of compounds. With the recent blockbuster success of several covalent inhibitors, and the development of new chemical proteomic strategies to broadly identify reactive, ligandable and posttranslationally modified cysteines, cysteine profiling is poised to enable the development of new potent and selective chemical probes and even, in some cases, new drugs.

The 91-205 amino acid region of AcMNPV ORF34 (Ac34), which comprises a potential C3H zinc finger, is required for its nuclear localization and optimal virus multiplication

During baculovirus infection, most viral proteins must be imported to the nucleus to support virus multiplication. Autographa californica multiple nucleopolyhedrovirus (AcMNPV) orf34 (ac34) is an alphabaculovirus unique gene that is required for optimal virus production. Ac34 distributes in both the cytoplasm and the nuclei of virus-infected Sf9 cells, but contains no conventional nuclear localization signal (NLS). In this study, we investigated the nuclear targeting domains in Ac34. Transient expression assays showed that Ac34 localized in both the cytoplasm and the nuclei of Sf9 cells, indicating that no viral protein is required for Ac34 nuclear localization. Subcellular localization analysis of Ac34 truncations and internal deletions fused with green fluorescent protein in plasmid-transfected Sf9 cells identified that the 91-205 amino acid (aa) region is required for Ac34 nuclear localization. Mutations in a potential C3H zinc finger (aa 116-131) in Ac34 resulted in exclusive cytoplasmic distribution of GFP:Ac34, suggesting that the zinc finger is required for Ac34 nuclear localization. To assess the functional importance of Ac34 in the nucleus during virus replication, recombinant AcMNPV bacmids containing a series of Ac34 truncations, internal deletions, or site mutations fused with HA tags were constructed. Subcellular localization analysis showed that Ac34 with internal deletions in aa 91-205 or site mutations in the potential zinc finger was predominantly distributed in the cytoplasm. Viral plaque assays and virus growth curves indicated that disruption of Ac34 nuclear localization significantly impaired virus replication. Taken together, our findings demonstrated that the nuclear localization of Ac34 requires the 91-205 aa region and its nuclear localization is essential for optimal virus replication.

Recognizing metal and acid radical ion-binding sites by integrating ab initio modeling with template-based transferals

MOTIVATION

More than half of proteins require binding of metal and acid radical ions for their structure and function. Identification of the ion-binding locations is important for understanding the biological functions of proteins. Due to the small size and high versatility of the metal and acid radical ions, however, computational prediction of their binding sites remains difficult.

RESULTS

We proposed a new ligand-specific approach devoted to the binding site prediction of 13 metal ions (Zn²⁺, Cu²⁺, Fe²⁺, Fe³⁺, Ca²⁺, Mg²⁺, Mn²⁺, Na⁺, K⁺) and acid radical ion ligands (CO3^2-, NO2^-, SO4^2-, PO4^3-) that are most frequently seen in protein databases. A sequence-based ab initio model is first trained on sequence profiles, where a modified AdaBoost algorithm is extended to balance binding and non-binding residue samples. A composite method IonCom is then developed to combine the ab initio model with multiple threading alignments for further improving the robustness of the binding site predictions. The pipeline was tested using 5-fold cross validations on a comprehensive set of 2,100 non-redundant proteins bound with 3,075 small ion ligands. Significant advantage was demonstrated compared with the state of the art ligand-binding methods including COACH and TargetS for high-accuracy ion-binding site identification. Detailed data analyses show that the major advantage of IonCom lies at the integration of complementary ab initio and template-based components. Ion-specific feature design and binding library selection also contribute to the improvement of small ion ligand binding predictions.

AVAILABILITY AND IMPLEMENTATION

http://zhanglab.ccmb.med.umich.edu/IonCom CONTACT: hxz@imut.edu.cn or zhng@umich.eduSupplementary information: Supplementary data are available at Bioinformatics online.

XLOS-observed mutations of MID1 Bbox1 domain cause domain unfolding

MID1 catalyzes the ubiquitination of the protein alpha4 and the catalytic subunit of protein phosphatase 2A. Mutations within the MID1 Bbox1 domain are associated with X-linked Opitz G syndrome (XLOS). Our functional assays have shown that mutations of Ala130 to Val or Thr, Cys142 to Ser and Cys145 to Thr completely disrupt the polyubiquitination of alpha4. Using NMR spectroscopy, we characterize the effect of these mutations on the tertiary structure of the Bbox1 domain by itself and in tandem with the Bbox2 domain. The mutation of either Cys142 or Cys145, each of which is involved in coordinating one of the two zinc ions, results in the collapse of signal dispersion in the HSQC spectrum of the Bbox1 domain indicating that the mutant protein structure is unfolded. Each mutation caused the coordination of both zinc ions, which are ∼ 13 Å apart, to be lost. Although Ala130 is not involved in the coordination of a zinc ion, the Ala130Thr mutant Bbox1 domain yields a poorly dispersed HSQC spectrum similar to those of the Cys142Ser and Cys145Thr mutants. Interestingly, neither cysteine mutation affects the structure of the adjacent Bbox2 domain when the two Bbox domains are engineered in their native tandem Bbox1-Bbox2 protein construct. Dynamic light scattering measurements suggest that the mutant Bbox1 domain has an increased propensity to form aggregates compared to the wild type Bbox1 domain. These studies provide insight into the mechanism by which mutations observed in XLOS affect the structure and function of the MID1 Bbox1 domain.

Fluxes in "free" and total zinc are essential for progression of intraerythrocytic stages of Plasmodium falciparum

Dynamic fluxes in the concentration of ions and small molecules are fundamental features of cell signaling, differentiation, and development. Similar roles for fluxes in transition metal concentrations are less well established. Here, we show that massive zinc fluxes are essential in the infection cycle of an intracellular eukaryotic parasite. Using single-cell quantitative imaging, we show that growth of the blood-stage Plasmodium falciparum parasite requires acquisition of 30 million zinc atoms per erythrocyte before host cell rupture, corresponding to a 400% increase in total zinc concentration. Zinc accumulates in a freely available form in parasitophorous compartments outside the food vacuole, including mitochondria. Restriction of zinc availability via small molecule treatment causes a drop in mitochondrial membrane potential and severely inhibits parasite growth. Thus, extraordinary zinc acquisition and trafficking are essential for parasite development.

Taking into Account the Ion-induced Dipole Interaction in the Nonbonded Model of Ions

Metal ions exist in almost half of the proteins in the protein databank and they serve as structural, electron-transfer and catalytic elements in the metabolic processes of organisms. Molecular Dynamics (MD) simulation is a powerful tool that provides information about biomolecular systems at the atomic level. Coupled with the growth in computing power, algorithms like the Particle Mesh Ewald (PME) method have become the accepted standard when dealing with long-range interactions in MD simulations. The nonbonded model of metal ions consists of an electrostatic plus 12-6 Lennard Jones (LJ) potential and is used largely because of its speed relative to more accurate models. In previous work we found that ideal parameters do not exist that reproduce several experimental properties for M(II) ions simultaneously using the nonbonded model coupled with the PME method due to the underestimation of metal ion-ligand interactions. Via a consideration of the nature of the nonbonded model, we proposed that the observed error largely arises from overlooking charge-induced dipole interactions. The electrostatic plus 12-6 LJ potential model works reasonably well for neutral systems but does struggle with more highly charged systems. In the present work we designed and parameterized a new nonbonded model for metal ions by adding a 1/r4 term to the 12-6 model. We call it the 12-6-4 LJ-type nonbonded model due to its mathematical construction. Parameters were determined for 16 +2 metal ions for the TIP3P, SPC/E and TIP4PEW water models. The final parameters reproduce the experimental hydration free energies (HFE), ion-oxygen distances (IOD) in the first solvation shell and coordination numbers (CN) accurately for the metal ions investigated. Preliminary tests on MgCl2 at different concentrations in aqueous solution and Mg2+--nucleic acid systems show reasonable results suggesting that the present parameters can work in mixed systems. The 12-6-4 LJ-type nonbonded model is readily adopted into standard force fields like AMBER, CHARMM and OPLS-AA with only a modest computational overhead. The new nonbonded model doesn't consider charge-transfer effects explicitly and, hence, may not suitable for the simulation of systems where charge-transfer effects play a decisive role.

Rational Design of Particle Mesh Ewald Compatible Lennard-Jones Parameters for +2 Metal Cations in Explicit Solvent

Metal ions play significant roles in biological systems. Accurate molecular dynamics (MD) simulations on these systems require a validated set of parameters. Although there are more detailed ways to model metal ions, the nonbonded model, which employs a 12-6 Lennard-Jones (LJ) term plus an electrostatic potential is still widely used in MD simulations today due to its simple form. However, LJ parameters have limited transferability due to different combining rules, various water models and diverse simulation methods. Recently, simulations employing a Particle Mesh Ewald (PME) treatment for long-range electrostatics have become more and more popular owing to their speed and accuracy. In the present work we have systematically designed LJ parameters for 24 +2 metal (M(II)) cations to reproduce different experimental properties appropriate for the Lorentz-Berthelot combining rules and PME simulations. We began by testing the transferability of currently available M(II) ion LJ parameters. The results showed that there are differences between simulations employing Ewald summation with other simulation methods and that it was necessary to design new parameters specific for PME based simulations. Employing the thermodynamic integration (TI) method and performing periodic boundary MD simulations employing PME, allowed for the systematic investigation of the LJ parameter space. Hydration free energies (HFEs), the ion-oxygen distance in the first solvation shell (IOD) and coordination numbers (CNs) were obtained for various combinations of the parameters of the LJ potential for four widely used water models (TIP3P, SPC/E, TIP4P and TIP4PEW). Results showed that the three simulated properties were highly correlated. Meanwhile, M(II) ions with the same parameters in different water models produce remarkably different HFEs but similar structural properties. It is difficult to reproduce various experimental values simultaneously because the nonbonded model underestimates the interaction between the metal ions and water molecules at short range. Moreover, the extent of underestimation increases successively for the TIP3P, SPC/E, TIP4PEW and TIP4P water models. Nonetheless, we fitted a curve to describe the relationship between ε (the well depth) and radius (Rmin/2) from experimental data on noble gases to facilitate the generation of the best possible compromise models. Hence, by targeting different experimental values, we developed three sets of parameters for M(II) cations for three different water models (TIP3P, SPC/E and TIP4PEW). These parameters we feel represent the best possible compromise that can be achieved using the nonbonded model for the ions in combination with simple water models. From a computational uncertainty analysis we estimate that the uncertainty in our computed HFEs is on the order of ±1kcal/mol. Further improvements will require more advanced non-bonded models likely with inclusion of polarization.

Computationally characterizing and comprehensive analysis of zinc-binding sites in proteins

Zinc is one of the most essential metals utilized by organisms, and zinc-binding proteins play an important role in a variety of biological processes such as transcription regulation, cell metabolism and apoptosis. Thus, characterizing the precise zinc-binding sites is fundamental to an elucidation of the biological functions and molecular mechanisms of zinc-binding proteins. Using systematic analyses of structural characteristics, we observed that 4-residue and 3-residue zinc-binding sites have distinctly specific geometric features. Based on the results, we developed the novel computational program Geometric REstriction for Zinc-binding (GRE4Zn) to characterize the zinc-binding sites in protein structures, by restricting the distances between zinc and its coordinating atoms. The comparison between GRE4Zn and analogous tools revealed that it achieved a superior performance. A large-scale prediction for structurally characterized proteins was performed with this powerful predictor, and statistical analyses for the results indicated zinc-binding proteins have come to be significantly involved in more complicated biological processes in higher species than simpler species during the course of evolution. Further analyses suggested that zinc-binding proteins are preferentially implicated in a variety of diseases and highly enriched in known drug targets, and the prediction of zinc-binding sites can be helpful for the investigation of molecular mechanisms. In this regard, these prediction and analysis results should prove to be highly useful be helpful for further biomedical study and drug design. The online service of GRE4Zn is freely available at: http://biocomp.ustc.edu.cn/gre4zn/. This article is part of a Special Issue entitled: Computational Proteomics, Systems Biology & Clinical Implications. Guest Editor: Yudong Cai.

EENdb: a database and knowledge base of ZFNs and TALENs for endonuclease engineering

We report here the construction of engineered endonuclease database (EENdb) (http://eendb.zfgenetics.org/), a searchable database and knowledge base for customizable engineered endonucleases (EENs), including zinc finger nucleases (ZFNs) and transcription activator-like effector nucleases (TALENs). EENs are artificial nucleases designed to target and cleave specific DNA sequences. EENs have been shown to be a very useful genetic tool for targeted genome modification and have shown great potentials in the applications in basic research, clinical therapies and agricultural utilities, and they are specifically essential for reverse genetics research in species where no other gene targeting techniques are available. EENdb contains over 700 records of all the reported ZFNs and TALENs and related information, such as their target sequences, the peptide components [zinc finger protein-/transcription activator-like effector (TALE)-binding domains, FokI variants and linker peptide/framework], the efficiency and specificity of their activities. The database also lists EEN engineering tools and resources as well as information about forms and types of EENs, EEN screening and construction methods, detection methods for targeting efficiency and many other utilities. The aim of EENdb is to represent a central hub for EEN information and an integrated solution for EEN engineering. These studies may help to extract in-depth properties and common rules regarding ZFN or TALEN efficiency through comparison of the known ZFNs or TALENs.

Zinc and zinc related enzymes in precancerous and cancerous tissue in the colon of dimethyl hydrazine treated rats

Trace element zinc deficiency or excess is implicated in the development or progression of some cancers. The exact role of zinc in the etiology of colon cancer is unclear. To cast light on this question, an experimental model of colon carcinogenesis was applied here. Six week old rats were given sub cutaneous injections of DMH (30 mg/kg body weight) twice a week for three months and sacrificed after 4 months (precancer model) and 6 months (cancer model). Plasma zinc levels showed a significant decrease (p<0.05) at 4 months and a greater significant decrease at 6 months (p<0.01) as compared with controls. In the large intestine there was a significant decrease in tissue zinc levels (p<0.005) and in CuZnSOD, and alkaline phosphatase activity (p<0.05) in the pre-cancerous model and a greater significant decrease in tissue zinc (p<0.0001), and in CuZnSOD and alkaline phosphatase activity (p<0.001), in the carcinoma model. The tissue zinc levels showed a significant decrease in the small intestine and stomach (p<0.005) and in liver (p<0.05) in the cancer model. 87% of the rats in the precancer group and 92% rats in the cancer group showed histological evidence of precancerous lesions and carcinomas respectively in the colon mucosa. This study suggests that the decrease in plasma zinc, tissue zinc and activity of zinc related enzymes are associated with the development of preneoplastic lesions and these biochemical parameters further decrease with progression to carcinoma in the colon.

Reactive cysteine in the structural Zn(2+) site of the C1B domain from PKCα

Structural cysteine-rich Zn(2+) sites that stabilize protein folds are considered to be unreactive. In this article, we identified a reactive cysteine residue, Cys151, in a treble-clef zinc finger with a Cys(3)His coordination sphere. The protein in question is the C1B domain of Protein Kinase Cα (PKCα). Mass-tagging cysteine assays of several C1B variants were employed to ascertain the site specificity of the covalent modification. The reactivity of Cys151 in C1B also manifests itself in the structural dynamics of the Zn(2+) coordination sphere where the Sγ of Cys151 alternates between the Zn(2+)-bound thiolate and free thiol states. We used NMR-detected pH titrations, ZZ-exchange spectroscopy, and residual dipolar coupling (RDC)-driven structure refinement to characterize the two exchanging conformations of C1B that differ in zinc coordination. Our data suggest that Cys151 serves as an entry point for the reactive oxygen species that activate PKCα in a process involving Zn(2+) release.

A BEN-domain-containing protein associates with heterochromatin and represses transcription

In eukaryotes, higher order chromatin structure governs crucial cellular processes including DNA replication, transcription and post-transcriptional gene regulation. Specific chromatin-interacting proteins play vital roles in the maintenance of chromatin structure. We have identified BEND3, a quadruple BEN domain-containing protein that is highly conserved amongst vertebrates. BEND3 colocalizes with HP1 and H3 trimethylated at K9 at heterochromatic regions in mammalian cells. Using an in vivo gene locus, we have been able to demonstrate that BEND3 associates with the locus only when it is heterochromatic and dissociates upon activation of transcription. Furthermore, tethering BEND3 inhibits transcription from the locus, indicating that BEND3 is involved in transcriptional repression through its interaction with histone deacetylases and Sall4, a transcription repressor. We further demonstrate that BEND3 is SUMOylated and that such modifications are essential for its role in transcriptional repression. Finally, overexpression of BEND3 causes premature chromatin condensation and extensive heterochromatinization, resulting in cell cycle arrest. Taken together, our data demonstrate the role of a novel heterochromatin-associated protein in transcriptional repression.

Using surface-assisted laser desorption/ionization mass spectrometry to detect proteins and protein-protein complexes

In this study, we combined surface-assisted laser desorption/ionization mass spectrometry (SALDI-MS) with HgTe nanostructures as matrix for the detection of several proteins (α1-antitrypsin, trypsin, IgG, protein G) and their complexes. We investigated the effects of several parameters (the concentration and nature of surfactants and metal ions, the pH, and concentration of the analytes in the sample matrixes) on the sensitivity of the detection of these proteins and their complexes. The presence of stabilizing Brij 76 surfactant and Zn(II) ions allowed the detection of weak protein complexes, such as α1-antitrypsin-trypsin and IgG-protein G complexes, at the picomole level. We observed multiply charged states at m/z 72,160 ([α1-antitrypsin + trypsin + H](+)) and 86,585 ([IgG + protein G + 2H](2+)) for the α1-antitrypsin-trypsin and IgG-protein G complexes, respectively. To the best of our knowledge, detection of weak protein complexes and determination of their stoichiometry have not been demonstrated previously when a combination of SALDI-MS and nanostructures were used. This simple and reproducible SALDI-MS approach using HgTe nanostructures holds great potential for the detection of other proteins and their complexes.

Multiple genes in the 15q13-q14 chromosomal region are associated with schizophrenia

OBJECTIVE

The chromosomal region, 15q13-q14, including the α7 nicotinic acetylcholine receptor gene, CHRNA7, is a replicated region for schizophrenia. This study fine-mapped genes at 15q13-q14 to determine whether the association is unique to CHRNA7.

METHODS

Family-based and case-control association studies were performed on Caucasian-non-Hispanic and African-American individuals from 120 families as well as 468 individual patients with schizophrenia and 144 well-characterized controls. Single-nucleotide polymorphism (SNP) markers were genotyped, and association analyses carried out for the outcomes of schizophrenia, smoking, and smoking in schizophrenia.

RESULTS

Three genes were associated with schizophrenia in both ethnic populations: TRPM1, KLF13, and RYR3. Two SNPs in CHRNA7 were associated with schizophrenia in African-Americans, and a second SNP in CHRNA7 was significant for an association with smoking and smoking in schizophrenia in Caucasians.

CONCLUSION

Results of these studies support association of the 15q13-q14 region with schizophrenia. The broad positive association suggests that more than one 15q gene may be contributing to the disorder, either in combination or through a regulatory mechanism.

Misexpression of testicular microRNA in sterile Xenopus hybrids points to tetrapod-specific microRNAs associated with male fertility

Spermatogenesis is one of the most complex biological processes undergone by any organism, making it susceptible to perturbations that result in male sterility. Research has demonstrated that mutant phenotypes can be obtained from the disruption of epigenetic modifications, which are commonly microRNA guided. Employing the Xenopus system, whereby homogametic interspecies males are always sterile, thus violating Haldane's Rule, we deep-sequenced testes-specific small-RNAs to identify microRNAs most frequently misexpressed between sterile hybrids and their fertile parental taxa. Using these data, we cross-referenced our expression information with previously published mouse (Mus musculus) data and identified a subset of seven microRNAs common to both (miR-338, miR-222, miR-18, miR-30, miR-10, miR-196, and miR-365). We propose that these microRNAs are likely critical for spermatogenesis in all tetrapods, having retained testicular expression across ~350 million years of evolution (Amphibian-Mammal split). Gene targets of six of these microRNAs are known, and all the six associate with zinc and zinc finger proteins (both previously found critical in male fertility), and three with Hox genes (some of which have also previously been deemed critical for testicular development and male fertility). Expression information for these targets revealed that all those associated with zinc have previously been found to express in mammalian testes. One Hox target has known mammalian testicular expression, two have close relatives with known mammalian testicular expression, and two more are associated with proteins known to have mammalian testicular expression. In addition, miR-222 has prior association with spermatogenesis, and miR-30 has been found to be abundantly expressed in both mouse and human testes.

Cysteine and histidine shuffling: mixing and matching cysteine and histidine residues in zinc finger proteins to afford different folds and function

Zinc finger proteins utilize zinc for structural purposes: zinc binds to a combination of cysteine and histidine ligands in a tetrahedral coordination geometry facilitating protein folding and function. While much is known about the classical zinc finger proteins, which utilize a Cys(2)His(2) ligand set to coordinate zinc and fold into an anti-parallel beta sheet/alpha helical fold, there are thirteen other families of 'non-classical' zinc finger proteins for which relationships between metal coordination and protein structure/function are less defined. This 'Perspective' article focuses on two classes of these non-classical zinc finger proteins: Cys(3)His type zinc finger proteins and Cys(2)His(2)Cys type zinc finger proteins. These proteins bind zinc in a tetrahedral geometry, like the classical zinc finger proteins, yet they adopt completely different folds and target different oligonucleotides. Our current understanding of the relationships between ligand set, metal ion, fold and function for these non-classical zinc fingers is discussed.

Are zinc-finger domains of protein kinase C dynamic structures that unfold by lipid or redox activation?

Protein kinase C (PKC) is activated by lipid second messengers or redox action, raising the question whether these activation modes involve the same or alternate mechanisms. Here we show that both lipid activators and oxidation target the zinc-finger domains of PKC, suggesting a unifying activation mechanism. We found that lipid agonist-binding or redox action leads to zinc release and disassembly of zinc fingers, thus triggering large-scale unfolding that underlies conversion to the active enzyme. These results suggest that PKC zinc fingers, originally considered purely structural devices, are in fact redox-sensitive flexible hinges, whose conformation is controlled both by redox conditions and lipid agonists.

Molecular dynamics investigation of dipeptide-transition metal salts in aqueous solutions

Molecular dynamics (MD) simulations of glycylglycine dipeptide with transition metal ions (Mn(2+), Fe(2+), Co(2+), Ni(2+), Cu(2+), and Zn(2+)) in aqueous solutions have been carried out to get an insight into the solvation structure, intermolecular interactions, and salt effects in these systems. The solvation structure and hydrogen bonding were described in terms of radial distribution function (RDF) and spatial distribution function (SDF). The dynamical properties of the solvation structure were also analyzed in terms of diffusion and residence times. The simulation results show the presence of a well-defined first hydration shell around the dipeptide, with water molecules forming hydrogen bonds to the polar groups of the dipeptide. This shell is, however, affected by the strong electric field of divalent metal ions, which at higher ion concentrations lead to the shift in the dipeptide-water RDFs. Higher salt concentrations lead also to increased residence times and slower diffusion rates. In general, smaller ions (Cu(2+), Zn(2+)) demonstrate stronger binding to dipeptide than the larger ones (Fe(2+), Mn(2+)). Simulations do not show any stronger association of peptide molecules indicating their dissolution in water. The above results may be of potential interest to future researchers on these molecular interactions.

Energy variational analysis of ions in water and channels: Field theory for primitive models of complex ionic fluids

Ionic solutions are mixtures of interacting anions and cations. They hardly resemble dilute gases of uncharged noninteracting point particles described in elementary textbooks. Biological and electrochemical solutions have many components that interact strongly as they flow in concentrated environments near electrodes, ion channels, or active sites of enzymes. Interactions in concentrated environments help determine the characteristic properties of electrodes, enzymes, and ion channels. Flows are driven by a combination of electrical and chemical potentials that depend on the charges, concentrations, and sizes of all ions, not just the same type of ion. We use a variational method EnVarA (energy variational analysis) that combines Hamilton's least action and Rayleigh's dissipation principles to create a variational field theory that includes flow, friction, and complex structure with physical boundary conditions. EnVarA optimizes both the action integral functional of classical mechanics and the dissipation functional. These functionals can include entropy and dissipation as well as potential energy. The stationary point of the action is determined with respect to the trajectory of particles. The stationary point of the dissipation is determined with respect to rate functions (such as velocity). Both variations are written in one Eulerian (laboratory) framework. In variational analysis, an "extra layer" of mathematics is used to derive partial differential equations. Energies and dissipations of different components are combined in EnVarA and Euler-Lagrange equations are then derived. These partial differential equations are the unique consequence of the contributions of individual components. The form and parameters of the partial differential equations are determined by algebra without additional physical content or assumptions. The partial differential equations of mixtures automatically combine physical properties of individual (unmixed) components. If a new component is added to the energy or dissipation, the Euler-Lagrange equations change form and interaction terms appear without additional adjustable parameters. EnVarA has previously been used to compute properties of liquid crystals, polymer fluids, and electrorheological fluids containing solid balls and charged oil droplets that fission and fuse. Here we apply EnVarA to the primitive model of electrolytes in which ions are spheres in a frictional dielectric. The resulting Euler-Lagrange equations include electrostatics and diffusion and friction. They are a time dependent generalization of the Poisson-Nernst-Planck equations of semiconductors, electrochemistry, and molecular biophysics. They include the finite diameter of ions. The EnVarA treatment is applied to ions next to a charged wall, where layering is observed. Applied to an ion channel, EnVarA calculates a quick transient pile-up of electric charge, transient and steady flow through the channel, stationary "binding" in the channel, and the eventual accumulation of salts in "unstirred layers" near channels. EnVarA treats electrolytes in a unified way as complex rather than simple fluids. Ad hoc descriptions of interactions and flow have been used in many areas of science to deal with the nonideal properties of electrolytes. It seems likely that the variational treatment can simplify, unify, and perhaps derive and improve those descriptions.

Mipu1, a novel rat zinc-finger protein, inhibits transcriptional activities of AP-1 and SRE in mitogen-activated protein kinase signaling pathway

Mipu1 is a novel rat gene recently identified in our lab. Mipu1 cDNA contains a 1,824 bp open reading frame (ORF) and encoded a 608 amino acid protein with an N-terminal Krüppel-associated box (KRAB) domain and classical zinc finger C(2)H(2) motifs in the C-terminus. Mipu1 protein is located in the nuclei. Fused to Gal-4 DNA-binding domain and cotransfected with pG5-luc, Mipu1 played a transcriptional suppressive effect. Deletion analysis with a series of truncated fusion proteins indicated that the KRAB motif was a basal repression domain. Overexpression of Mipu1 in H9c2 myogenic cells inhibited the transcriptional activities of SRE and AP-1. RNAi of Mipu1 in H9c2 myogenic cells activated the transcriptional activities of SRE and AP-1. These results suggested that Mipu1 protein might act as a transcriptional repressor in mitogen-activated protein kinase (MAPK) signaling pathway to mediate cellular functions.

Mapping of nucleocapsid residues important for HIV-1 genomic RNA dimerization and packaging

Retroviral genomic RNA (gRNA) dimerization appears essential for viral infectivity, and the nucleocapsid protein (NC) of human immunodeficiency virus type 1 (HIV-1) facilitates HIV-1 gRNA dimerization. To identify the relevant and dispensable positions of NC, 34 of its 55 residues were mutated, individually or in small groups, in a panel of 40 HIV-1 mutants prepared by site-directed mutagenesis. It was found that the amino-terminus, the proximal zinc finger, the linker, and the distal zinc finger of NC each contributed roughly equally to efficient HIV-1 gRNA dimerization. The N-terminal and linker segments appeared to play predominantly electrostatic and steric roles, respectively. Mutating the hydrophobic patch of either zinc finger, or substituting alanines for their glycine doublet, was as disabling as deleting the corresponding finger. Replacing the CysX(2)CysX(4)HisX(4)Cys motif of either finger by CysX(2)CysX(4)CysX(4)Cys or CysX(2)CysX(4)HisX(4)His, interchanging the zinc fingers or, replacing one zinc finger by a copy of the other one, had generally intermediate effects; among these mutations, the His23-->Cys substitution in the N-terminal zinc finger had the mildest effect. The charge of NC could be increased or decreased by up to 18%, that of the linker could be reduced by 75% or increased by 50%, and one or two electric charges could be added or subtracted from either zinc finger, without affecting gRNA dimerization. Shortening, lengthening, or making hydrophobic the linker was as disabling as deleting the N-terminal or the C-terminal zinc finger, but a neutral and polar linker was innocuous. The present work multiplies by 4 and by 33 the number of retroviral and lentiviral NC mutations known to inhibit gRNA dimerization, respectively. It shows the first evidence that gRNA dimerization can be inhibited by: 1) mutations in the N-terminus or the linker of retroviral NC; 2) mutations in the proximal zinc finger of lentiviral NC; 3) mutations in the hydrophobic patch or the conserved glycines of the proximal or the distal retroviral zinc finger. Some NC mutations impaired gRNA dimerization more than mutations inactivating the viral protease, indicating that gRNA dimerization may be stimulated by the NC component of the Gag polyprotein. Most, but not all, mutations inhibited gRNA packaging; some had a strong effect on virus assembly or stability.