-In silico functional characterization of a double histone fold domain from the Heliothis zea virus 1

BACKGROUND: Histones are short proteins involved in chromatin packaging; in eukaryotes, two H2a-H2b and H3-H4 histone dimers form the nucleosomal core, which acts as the fundamental DNA-packaging element. The double histone fold is a rare globular protein fold in which two consecutive regions characterized by the typical structure of histones assemble together, thus originating a histone pseudodimer. This fold is included in a few prokaryotic histones and in the regulatory region of guanine nucleotide exchange factors of the Sos family. For the prokaryotic histones, there is no direct structural counterpart in the nucleosomal core particle, while the pseudodimer from Sos proteins is very similar to the dimer formed by histones H2a and H2b RESULTS: The absence of a H3-H4-like histone pseudodimer in the available structural databases prompted us to search for proteins that could assume such fold. The application of several secondary structure prediction and fold recognition methods allowed to show that the viral protein gi|22788712 is compatible with the structure of a H3-H4-like histone pseudodimer. Further in silico analyses revealed that this protein module could retain the ability of mediating protein-DNA interactions, and could consequently act as a DNA-binding domain. CONCLUSION: Our results suggest a possible functional role in viral pathogenicity for this novel double histone fold domain; thus, the computational analyses here reported will be helpful in directing future biochemical studies on gi|22788712 protein

GPR17: Molecular modeling and dynamics studies of the 3-D structure and purinergic ligand binding features in comparison with P2Y receptors

<p>Abstract</p> <p>Background</p> <p>GPR17 is a G-protein-coupled receptor located at intermediate phylogenetic position between two distinct receptor families: the P2Y and CysLT receptors for extracellular nucleotides and cysteinyl-LTs, respectively. We previously showed that GPR17 can indeed respond to both classes of endogenous ligands and to synthetic compounds active at the above receptor families, thus representing the first fully characterized non-peptide "hybrid" GPCR. In a rat brain focal ischemia model, the selective <it>in vivo </it>knock down of GPR17 by anti-sense technology or P2Y/CysLT antagonists reduced progression of ischemic damage, thus highlighting GPR17 as a novel therapeutic target for stroke. Elucidation of the structure of GPR17 and of ligand binding mechanisms are the necessary steps to obtain selective and potent drugs for this new potential target. On this basis, a 3-D molecular model of GPR17 embedded in a solvated phospholipid bilayer and refined by molecular dynamics simulations has been the first aim of this study. To explore the binding mode of the "purinergic" component of the receptor, the endogenous agonist UDP and two P2Y receptor antagonists demonstrated to be active on GPR17 (MRS2179 and cangrelor) were then modeled on the receptor.</p> <p>Results</p> <p>Molecular dynamics simulations suggest that GPR17 nucleotide binding pocket is similar to that described for the other P2Y receptors, although only one of the three basic residues that have been typically involved in ligand recognition is conserved (Arg255). The binding pocket is enclosed between the helical bundle and covered at the top by EL2. Driving interactions are H-bonds and salt bridges between the 6.55 and 6.52 residues and the phosphate moieties of the ligands. An "accessory" binding site in a region formed by the EL2, EL3 and the Nt was also found.</p> <p>Conclusion</p> <p>Nucleotide binding to GPR17 occurs on the same receptor regions identified for already known P2Y receptors. Agonist/antagonist binding mode are similar, but not identical. An accessory external binding site could guide small ligands to the deeper principal binding site in a multi-step mechanism of activation. The nucleotide binding pocket appears to be unable to allocate the leukotrienic type ligands in the same effective way.</p

Fast Generation of Broken-Symmetry States in a Large System Including Multiple Iron-Sulfur Assemblies: Investigation of QM/MM Energies, Clusters Charges, and Spin Populations

A density functional theory study is presented regarding the energetics and the Mulliken population analyses of a quantum mechanical/molecular mechanical (QM/MM) system including multiple iron-sulfur clusters in the QM region. The [FeFe]-hydrogenase from Desulfovibrio desulfuricans was studied, and both the active site (an Fe6S6 assembly generally referred to as the H-cluster) and an ancillary Fe4S4 site were treated at the BP86-RI/TZVP level. The antiferromagnetic coupling that characterizes both sites was modeled using the broken-symmetry (BS) approach. For such a QM system, 36 different BS couplings can be defined, depending on the localization of spin excess on the various spin centers. All the BS states were obtained by means of an effective and simple method for spin localization, that is here described and compared with more sophisticated approaches already available in literature. The variation of the QM/MM energy among the various geometry-optimized protein models was found to be less than 25 kJ mol(-1). This energy variation almost doubles if no geometry optimization is performed. A detailed analysis of the additive nature of these variations in QM/MM energy is reported. The Mulliken charges show very small variations among the 36 BS states, whereas the Mulliken spin populations were found to be somewhat more variable. The relevance of such variations is discussed in light of the available Mossbauer and Electron Paramagnetic Resonance (EPR) spectroscopic data for the enzyme. Finally, the influence of the basis set on the spin populations, charges, and structural parameters of the models was investigated, by means of QM/MM computations on the same system at the BP86-RI/SVP level. (C) 2010 Wiley Periodicals, Inc. Int J Quantum Chem 111: 3949-3960, 201

Mechanistic and Physiological Implications of the Interplay among Iron-Sulfur Clusters in [FeFe]-Hydrogenases. A QM/MM Perspective

Key stereoelectronic properties of Desulfovibrio desulfuricans [FeFe]-hydrogenase (DdH) were investigated by quantum mechanical description of its complete inorganic core, which includes a Fe6S6 active site (the H-cluster), as well as two ancillary Fe4S4 assemblies (the F and F' clusters). The partially oxidized, active-ready form of DdH is able to efficiently bind dihydrogen, thus starting H-2 oxidation catalysis. The calculations allow us to unambiguously assign a mixed Fe(H)Fe(I) state to the catalytic core of the active-ready enzyme and show that H-2 uptake exerts subtle, yet crucial influences on the redox properties of DdH. In fact, H-2 binding can promote electron transfer from the H-cluster to the solvent-exposed F'-cluster, thanks to a 50% decrease of the energy gap between the HOMO (that is localized on the H-cluster) and the LUMO (which is centered on the F'-cluster). Our results also indicate that the binding of the redox partners of DdH in proximity of its F'-cluster can trigger one-electron oxidation of the H-2-bound enzyme, a process that is expected to have an important role in H-2 activation. Our findings are analyzed not only from a mechanistic perspective, but also in consideration of the physiological role of DdH. In fact, this enzyme is known to be able to catalyze both the oxidation and the evolution of H-2, depending on the cellular metabolic requirements. Hints for the design of targeted mutations that could lead to the enhancement of the oxidizing properties of DdH are proposed and discussed

Validation of protein models by a neural network approach

Background: The development and improvement of reliable computational methods designed to evaluate the quality of protein models is relevant in the context of protein structure refinement, which has been recently identified as one of the bottlenecks limiting the quality and usefulness of protein structure prediction. Results: In this contribution, we present a computational method (Artificial Intelligence Decoys Evaluator: AIDE) which is able to consistently discriminate between correct and incorrect protein models. In particular, the method is based on neural networks that use as input 15 structural parameters, which include energy, solvent accessible surface, hydrophobic contacts and secondary structure content. The results obtained with AIDE on a set of decoy structures were evaluated using statistical indicators such as Pearson correlation coefficients, Znat, fraction enrichment, as well as ROC plots. It turned out that AIDE performances are comparable and often complementary to available state-of-the-art learning-based methods. Conclusion: In light of the results obtained with AIDE, as well as its comparison with available learning-based methods, it can be concluded that AIDE can be successfully used to evaluate the quality of protein structures. The use of AIDE in combination with other evaluation tools is expected to further enhance protein refinement effort

Catalytic Mechanism of Fungal Lytic Polysaccharide Monooxygenases Investigated by First-Principles Calculations

Lytic polysaccharide monooxygenases (LPMOs) are Cu-containing enzymes that facilitate the degradation of recalcitrant polysaccharides by the oxidative cleavage of glycosidic bonds. They are gaining rapidly increasing attention as key players in biomass conversion, especially for the production of second-generation biofuels. Elucidation of the detailed mechanism of the LPMO reaction is a major step toward the assessment and optimization of LPMO efficacy in industrial biotechnology, paving the way to utilization of sustainable fuel sources. Here, we used density functional theory calculations to study the reaction pathways suggested to date, exploiting a very large active-site model for a fungal AA9 LPMO and using a celloheptaose unit as a substrate mimic. We identify a copper oxyl intermediate as being responsible for H-atom abstraction from the substrate, followed by a rapid, water-assisted hydroxyl rebound, leading to substrate hydroxylation

An Acidic Loop and Cognate Phosphorylation Sites Define a Molecular Switch That Modulates Ubiquitin Charging Activity in Cdc34-Like Enzymes

E2 ubiquitin-conjugating enzymes are crucial mediators of protein ubiquitination, which strongly influence the ultimate fate of the target substrates. Recently, it has been shown that the activity of several enzymes of the ubiquitination pathway is finely tuned by phosphorylation, an ubiquitous mechanism for cellular regulation, which modulates protein conformation. In this contribution, we provide the first rationale, at the molecular level, of the regulatory mechanism mediated by casein kinase 2 (CK2) phosphorylation of E2 Cdc34-like enzymes. In particular, we identify two co-evolving signature elements in one of the larger families of E2 enzymes: an acidic insertion in β4α2 loop in the proximity of the catalytic cysteine and two conserved key serine residues within the catalytic domain, which are phosphorylated by CK2. Our investigations, using yeast Cdc34 as a model, through 2.5 µs molecular dynamics simulations and biochemical assays, define these two elements as an important phosphorylation-controlled switch that modulates opening and closing of the catalytic cleft. The mechanism relies on electrostatic repulsions between a conserved serine phosphorylated by CK2 and the acidic residues of the β4α2 loop, promoting E2 ubiquitin charging activity. Our investigation identifies a new and unexpected pivotal role for the acidic loop, providing the first evidence that this loop is crucial not only for downstream events related to ubiquitin chain assembly, but is also mandatory for the modulation of an upstream crucial step of the ubiquitin pathway: the ubiquitin charging in the E2 catalytic cleft