Molecular Modelling of the Ni(II)-Responsive Synechocystis PCC 6803 Transcriptional Regulator InrS in the Metal Bound Form

InrS (internal nickel-responsive sensor) is a transcriptional regulator found in cyanobacteria that represses the transcription of the nickel exporter NrsD in the apo form and de-represses expression of the exporter upon Ni(II) binding. Although a crystal structure of apo-InrS from Synechocystis PCC 6803 has been reported, no structure of the protein with metal ions bound is available. Here we report the results of a computational study aimed to reconstruct the metal binding site by taking advantage of recent X-ray absorption spectroscopy (XAS) data and to envisage the structural rearrangements occurring upon Ni(II) binding. The modelled Ni(II) binding site shows a square planar geometry consistent with experimental data. The structural details of the conformational changes occurring upon metal binding are also discussed in the framework of trying to rationalize the different affnity of the apo- and holo-forms of the protein for DNA

Structure of the Intermolecular Complex between Plastocyanin and Cytochrome f from Spinach

In oxygenic photosynthesis, plastocyanin shuttles electrons between the membrane-bound complexes cytochrome b6f and photosystem I. The homologous complex between cytochrome f and plastocyanin, both from spinach, is the object of this study. The solution structure of the reduced spinach plastocyanin was determined using high field NMR spectroscopy, whereas the model structure of oxidized cytochrome f was obtained by homology modeling calculations and molecular dynamics. The model structure of the intermolecular complex was calculated using the program AUTODOCK, taking into account biological information obtained from mutagenesis experiments. The best electron transfer pathway from the heme group of cytochrome f to the copper ion of plastocyanin was calculated using the program HARLEM, obtaining a coupling decay value of 1.8 x 10(-4). Possible mechanisms of interaction and electron transfer between plastocyanin and cytochrome f were discussed considering the possible formation of a supercomplex that associates one cytochrome b6f, one photosystem I, and one plastocyanin

PacDOCK: A Web Server for Positional Distance-Based and Interaction-Based Analysis of Docking Results

Molecular docking is a key method for structure-based drug design used to predict the conformations assumed by small drug-like ligands when bound to their target. However, the evaluation of molecular docking studies can be hampered by the lack of a free and easy to use platform for the complete analysis of results obtained by the principal docking programs. To this aim, we developed PacDOCK, a freely available and user-friendly web server that comprises a collection of tools for positional distance-based and interaction-based analysis of docking results, which can be provided in several file formats. PacDOCK allows a complete analysis of molecular docking results through root mean square deviation (RMSD) calculation, molecular visualization, and cluster analysis of docked poses. The RMSD calculation compares docked structures with a reference structure, also when atoms are randomly labelled, and their conformational and positional differences can be visualised. In addition, it is possible to visualise a ligand into the target binding pocket and investigate the key receptor-ligand interactions. Moreover, PacDOCK enables the clustering of docking results by identifying a restrained number of clusters from many docked poses. We believe that PacDOCK will contribute to facilitating the analysis of docking results to improve the efficiency of computer-aided drug design

Computational Study of Helicase from SARS-CoV-2 in RNA-Free and Engaged Form

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the pandemic that broke out in 2020 and continues to be the cause of massive global upheaval. Coronaviruses are positive-strand RNA viruses with a genome of ~30 kb. The genome is replicated and transcribed by RNA-dependent RNA polymerase together with accessory factors. One of the latter is the protein helicase (NSP13), which is essential for viral replication. The recently solved helicase structure revealed a tertiary structure composed of five domains. Here, we investigated NSP13 from a structural point of view, comparing its RNA-free form with the RNA-engaged form by using atomistic molecular dynamics (MD) simulations at the microsecond timescale. Structural analyses revealed conformational changes that provide insights into the contribution of the different domains, identifying the residues responsible for domain–domain interactions in both observed forms. The RNA-free system appears to be more flexible than the RNA-engaged form. This result underlies the stabilizing role of the nucleic acid and the functional core role of these domains

Modelling eNvironment for Isoforms (MoNvIso): A general platform to predict structural determinants of protein isoforms in genetic diseases

: The seamless integration of human disease-related mutation data into protein structures is an essential component of any attempt to correctly assess the impact of the mutation. The key step preliminary to any structural modelling is the identification of the isoforms onto which mutations should be mapped due to there being several functionally different protein isoforms from the same gene. To handle large sets of data coming from omics techniques, this challenging task needs to be automatized. Here we present the MoNvIso (Modelling eNvironment for Isoforms) code, which identifies the most useful isoform for computational modelling, balancing the coverage of mutations of interest and the availability of templates to build a structural model of both the wild-type isoform and the related variants

Conformational equilibria in monomeric alpha-synuclein at the single molecule level

Natively unstructured proteins defy the classical "one sequence-one structure" paradigm of protein science. Monomers of these proteins in pathological conditions can aggregate in the cell, a process that underlies socially relevant neurodegenerative diseases such as Alzheimer and Parkinson. A full comprehension of the formation and structure of the so-called misfolded intermediates from which the aggregated states ensue is still lacking. We characterized the folding and the conformational diversity of alpha-synuclein (aSyn), a natively unstructured protein involved in Parkinson disease, by mechanically stretching single molecules of this protein and recording their mechanical properties. These experiments permitted us to directly observe directly and quantify three main classes of conformations that, under in vitro physiological conditions, exist simultaneously in the aSyn sample, including disordered and "beta-like" structures. We found that this class of "beta-like" structures is directly related to aSyn aggregation. In fact, their relative abundance increases drastically in three different conditions known to promote the formation of aSyn fibrils: the presence of Cu2+, the occurrence of the pathogenic A30P mutation, and high ionic strength. We expect that a critical concentration of aSyn with a "beta-like" structure must be reached to trigger fibril formation. This critical concentration is therefore controlled by a chemical equilibrium. Novel pharmacological strategies can now be tailored to act upstream, before the aggregation process ensues, by targeting this equilibrium. To this end, Single Molecule Force Spectroscopy can be an effective tool to tailor and test new pharmacological agents.Comment: 37 pages, 9 figures (including supplementary material

Cloning the barley nec3 disease lesion mimic mutant using complementation by sequencing

Disease lesion mimic (DLM) or necrotic mutants display necrotic lesions in the absence of pathogen infections. They can show improved resistance to some pathogens and their molecular dissection can contribute to revealing components of plant defense pathways. Although forward-genetics strategies to find genes causal to mutant phenotypes are available in crops, these strategies require the production of experimental cross populations, mutagenesis, or gene editing and are time- and resource-consuming or may have to deal with regulated plant materials. In this study, we described a collection of 34 DLM mutants in barley (Hordeum vulgare L.) and applied a novel method called complementation by sequencing (CBS), which enables the identification of the gene responsible for a mutant phenotype given the availability of two or more chemically mutagenized individuals showing the same phenotype. Complementation by sequencing relies on the feasibility to obtain all induced mutations present in chemical mutants and on the low probability that different individuals share the same mutated genes. By CBS, we identified a cytochrome P450 CYP71P1 gene as responsible for orange blotch DLM mutants, including the historical barley nec3 locus. By comparative phylogenetic analysis we showed that CYP71P1 gene family emerged early in angiosperm evolution but has been recurrently lost in some lineages including Arabidopsis thaliana (L.) Heynh. Complementation by sequencing is a straightforward cost-effective approach to clone genes controlling phenotypes in a chemically mutagenized collection. The TILLMore (TM) collection will be instrumental for understanding the molecular basis of DLM phenotypes and to contribute knowledge about mechanisms of host-pathogen interaction

Deciphering OPA1 mutations pathogenicity by combined analysis of human, mouse and yeast cell models

OPA1 is the major gene responsible for Dominant Optic Atrophy (DOA) and the syndromic form DOA “plus”. Over 370 OPA1 mutations have been identified so far, although their pathogenicity is not always clear. We have analyzed one novel and a set of known OPA1 mutations to investigate their impact on protein functions in primary skin fibroblasts and in two “ad hoc” generated cell systems: the MGM1/OPA1 chimera yeast model and the Opa1−/− MEFs model expressing the mutated human OPA1 isoform 1. The yeast model allowed us to confirm the deleterious effects of these mutations and to gain information on their dominance/recessivity. The MEFs model enhanced the phenotypic alteration caused by mutations, nicely correlating with the clinical severity observed in patients, and suggested that the DOA “plus” phenotype could be induced by the combinatorial effect of mitochondrial network fragmentation with variable degrees of mtDNA depletion. Overall, the two models proved to be valuable tools to functionally assess and define the deleterious mechanism and the pathogenicity of novel OPA1 mutations, and useful to testing new therapeutic interventions

Deciphering OPA1 mutations pathogenicity by combined analysis of human, mouse and yeast cell models

OPA1 is the major gene responsible for Dominant Optic Atrophy (DOA) and the syndromic form DOA “plus”. Over 370 OPA1 mutations have been identified so far, although their pathogenicity is not always clear. We have analyzed one novel and a set of known OPA1 mutations to investigate their impact on protein functions in primary skin fibroblasts and in two “ad hoc” generated cell systems: the MGM1/OPA1 chimera yeast model and the Opa1−/− MEFs model expressing the mutated human OPA1 isoform 1. The yeast model allowed us to confirm the deleterious effects of these mutations and to gain information on their dominance/recessivity. The MEFs model enhanced the phenotypic alteration caused by mutations, nicely correlating with the clinical severity observed in patients, and suggested that the DOA “plus” phenotype could be induced by the combinatorial effect of mitochondrial network fragmentation with variable degrees of mtDNA depletion. Overall, the two models proved to be valuable tools to functionally assess and define the deleterious mechanism and the pathogenicity of novel OPA1 mutations, and useful to testing new therapeutic interventions