Amino acid positions subject to multiple co-evolutionary constraints can be robustly identified by their eigenvector network centrality scores

As proteins evolve, amino acid positions key to protein structure or function are subject to mutational constraints. These positions can be detected by analyzing sequence families for amino acid conservation or for co-evolution between pairs of positions. Co-evolutionary scores are usually rank-ordered and thresholded to reveal the top pairwise scores, but they also can be treated as weighted networks. Here, we used network analyses to bypass a major complication of co-evolution studies: For a given sequence alignment, alternative algorithms usually identify different, top pairwise scores. We reconciled results from five commonly-used, mathematically divergent algorithms (ELSC, McBASC, OMES, SCA, and ZNMI), using the LacI/GalR and 1,6-bisphosphate aldolase protein families as models. Calculations used unthresholded co-evolution scores from which column-specific properties such as sequence entropy and random noise were subtracted; “central” positions were identified by calculating various network centrality scores. When compared among algorithms, network centrality methods, particularly eigenvector centrality, showed markedly better agreement than comparisons of the top pairwise scores. Positions with large centrality scores occurred at key structural locations and/or were functionally sensitive to mutations. Further, the top central positions often differed from those with top pairwise co-evolution scores: Instead of a few strong scores, central positions often had multiple, moderate scores. We conclude that eigenvector centrality calculations reveal a robust evolutionary pattern of constraints – detectable by divergent algorithms – that occur at key protein locations. Finally, we discuss the fact that multiple patterns co-exist in evolutionary data that, together, give rise to emergent protein functions

A Coevolutionary Residue Network at the Site of a Functionally Important Conformational Change in a Phosphohexomutase Enzyme Family

Coevolution analyses identify residues that co-vary with each other during evolution, revealing sequence relationships unobservable from traditional multiple sequence alignments. Here we describe a coevolutionary analysis of phosphomannomutase/phosphoglucomutase (PMM/PGM), a widespread and diverse enzyme family involved in carbohydrate biosynthesis. Mutual information and graph theory were utilized to identify a network of highly connected residues with high significance. An examination of the most tightly connected regions of the coevolutionary network reveals that most of the involved residues are localized near an interdomain interface of this enzyme, known to be the site of a functionally important conformational change. The roles of four interface residues found in this network were examined via site-directed mutagenesis and kinetic characterization. For three of these residues, mutation to alanine reduces enzyme specificity to ∼10% or less of wild-type, while the other has ∼45% activity of wild-type enzyme. An additional mutant of an interface residue that is not densely connected in the coevolutionary network was also characterized, and shows no change in activity relative to wild-type enzyme. The results of these studies are interpreted in the context of structural and functional data on PMM/PGM. Together, they demonstrate that a network of coevolving residues links the highly conserved active site with the interdomain conformational change necessary for the multi-step catalytic reaction. This work adds to our understanding of the functional roles of coevolving residue networks, and has implications for the definition of catalytically important residues

Peptidyl-Prolyl cis/trans Isomerase-Independent Functional NifH Mutant of Azotobacter vinelandii

Peptidyl-prolyl cis/trans isomerases (PPIases) play a pivotal role in catalyzing the correct folding of many prokaryotic and eukaryotic proteins that are implicated in a variety of biological functions, ranging from cell cycle regulation to bacterial infection. The nif accessory protein NifM, which is essential for the biogenesis of a functional NifH component of nitrogenase, is a PPIase. To understand the nature of the molecular signature that defines the NifM dependence of NifH, we screened a library of nifH mutants in the nitrogen-fixing bacterium Azotobacter vinelandii for mutants that acquired NifM independence. Here, we report that NifH can acquire NifM independence when the conserved Pro258 located in the C-terminal region of NifH, which wraps around the other subunit in the NifH dimer, is replaced by serine

In vitro transcription accurately predicts lac repressor phenotype in vivo in Escherichia coli

A multitude of studies have looked at the in vivo and in vitro behavior of the lac repressor binding to DNA and effector molecules in order to study transcriptional repression, however these studies are not always reconcilable. Here we use in vitro transcription to directly mimic the in vivo system in order to build a self consistent set of experiments to directly compare in vivo and in vitro genetic repression. A thermodynamic model of the lac repressor binding to operator DNA and effector is used to link DNA occupancy to either normalized in vitro mRNA product or normalized in vivo fluorescence of a regulated gene, YFP. An accurate measurement of repressor, DNA and effector concentrations were made both in vivo and in vitro allowing for direct modeling of the entire thermodynamic equilibrium. In vivo repression profiles are accurately predicted from the given in vitro parameters when molecular crowding is considered. Interestingly, our measured repressor–operator DNA affinity differs significantly from previous in vitro measurements. The literature values are unable to replicate in vivo binding data. We therefore conclude that the repressor-DNA affinity is much weaker than previously thought. This finding would suggest that in vitro techniques that are specifically designed to mimic the in vivo process may be necessary to replicate the native system

Impaired Mitophagy Plays a Role in Denervation of Neuromuscular Junctions in ALS Mice

Motor neurons in amyotrophic lateral sclerosis (ALS) patients and animal models show degeneration from the nerve terminal, known as dying-back neuropathy. To investigate the mechanism underlying this neuropathy, we analyzed the neuromuscular junctions (NMJs) and motor neuron cell bodies in SOD1G93A mice using electron microscopy. NMJs of SOD1G93A mice exhibited significantly higher numbers of autophagosomes and degenerated mitochondria compared to wild-type controls. Mitophagosomes were identified in the NMJ presynaptic terminals of wild-type mice and SOD1G93A mice. However, the number of mitophagosomes did not increase significantly in SOD1G93A NMJs indicating a defect in mitophagy, the autophagic process to degrade mitochondria. Consistent with this, proteins essential for mitophagy, p62/SQSTM1, Bnip3, Pink1, and Parkin were down-regulated in motor neurons in SOD1G93A mice. Importantly, SQSTM1 is one of the genes mutated in familial ALS patients. We evaluated the effect of impaired mitophagy on motor neurons by analyzing the double knockout mice of Pink1 and Parkin, two genes responsible for sensing depolarized mitochondria and delivering degenerated mitochondria to mitophagosomes. The double knockout mice exhibited NMJ degeneration, including axon swelling and NMJ fragmentation at 4 months of age. These phenotypes were rarely observed in wild-type control mice of the same age. The protein level of ATP synthase β subunit increased in the NMJ presynaptic terminals, suggesting the accumulation of mitochondria at NMJs of the double knockout mice. Importantly, NMJ denervation was observed in the double knockout mice. These data suggest that the reduced mitophagy function in motor neurons of SOD1G93A mice is one of the mechanisms causing degeneration of ALS NMJs