iBIS2Analyzer: a web server for a phylogeny-driven coevolution analysis of protein families.

Residue coevolution within and between proteins is used as a marker of physical interaction and/or residue functional cooperation. Pairs or groups of coevolving residues are extracted from multiple sequence alignments based on a variety of computational approaches. However, coevolution signals emerging in subsets of sequences might be lost if the full alignment is considered. iBIS2Analyzer is a web server dedicated to a phylogeny-driven coevolution analysis of protein families with different evolutionary pressure. It is based on the iterative version, iBIS2, of the coevolution analysis method BIS, Blocks in Sequences. iBIS2 is designed to iteratively select and analyse subtrees in phylogenetic trees, possibly large and comprising thousands of sequences. With iBIS2Analyzer, openly accessible at http://ibis2analyzer.lcqb.upmc.fr/, the user visualizes, compares and inspects clusters of coevolving residues by mapping them onto sequences, alignments or structures of choice, greatly simplifying downstream analysis steps. A rich and interactive graphic interface facilitates the biological interpretation of the results.

MyCLADE: a multi-source domain annotation server for sequence functional exploration.

The ever-increasing number of genomic and metagenomic sequences accumulating in our databases requires accurate approaches to explore their content against specific domain targets. MyCLADE is a user-friendly webserver designed for targeted functional profiling of genomic and metagenomic sequences based on a database of a few million probabilistic models of Pfam domains. It uses the MetaCLADE multi-source domain annotation strategy, modelling domains based on multiple probabilistic profiles. MyCLADE takes a list of protein sequences and possibly a target set of domains/clans as input and, for each sequence, it provides a domain architecture built from the targeted domains or from all Pfam domains. It is linked to the Pfam and QuickGO databases in multiple ways for easy retrieval of domain and clan information. E-value, bit-score, domain-dependent probability scores and logos representing the match of the model with the sequence are provided to help the user to assess the quality of each annotation. Availability and implementation: MyCLADE is freely available at http://www.lcqb.upmc.fr/myclade.

A large disordered region confers a wide spanning volume to vertebrate Suppressor of Fused as shown in a trans-species solution study.

Hedgehog (Hh) pathway inhibition by the conserved protein Suppressor of Fused (SuFu) is crucial to vertebrate development. By constrast, SuFu loss-of-function mutant has little effect in drosophila. Previous publications showed that the crystal structures of human and drosophila SuFu consist of two ordered domains that are capable of breathing motions upon ligand binding. However, the crystal structure of human SuFu does not give information about twenty N-terminal residues (IDR1) and an eighty-residue-long region predicted as disordered (IDR2) in the C-terminus, whose function is important for the pathway repression. These two intrinsically disordered regions (IDRs) are species-dependent. To obtain information about the IDR regions, we studied full-length SuFu's structure in solution, both with circular dichroism and small angle X-ray scattering, comparing drosophila, zebrafish and human species, to better understand this considerable difference. Our studies show that, in spite of similar crystal structures restricted to ordered domains, drosophila and vertebrate SuFu have very different structures in solution. The IDR2 of vertebrates spans a large area, thus enabling it to reach for partners and be accessible for post-translational modifications. Furthermore, we show that the IDR2 region is highly conserved within phyla but varies in length and sequence, with insects having a shorter disordered region while that of vertebrates is broad and mobile. This major variation may explain the different phenotypes observed upon SuFu removal.

COVTree: Coevolution in OVerlapped sequences by Tree analysis server.

Overlapping genes are commonplace in viruses and play an important role in their function and evolution. For these genes, molecular coevolution may be seen as a mechanism to decrease the evolutionary constraints of amino acid positions in the overlapping regions and to tolerate or compensate unfavorable mutations. Tracing these mutational sites, could help to gain insight on the direct or indirect effect of the mutations in the corresponding overlapping proteins. In the past, coevolution analysis has been used to identify residue pairs and coevolutionary signatures within or between proteins that served as markers of physical interactions and/or functional relationships. Coevolution in OVerlapped sequences by Tree analysis (COVTree) is a web server providing the online analysis of coevolving amino-acid pairs in overlapping genes, where residues might be located inside or outside the overlapping region. COVTree is designed to handle protein families with various characteristics, among which those that typically display a small number of highly conserved sequences. It is based on BIS2, a fast version of the coevolution analysis tool Blocks in Sequences (BIS). COVTree provides a rich and interactive graphical interface to ease biological interpretation of the results and it is openly accessible at http://www.lcqb.upmc.fr/COVTree/.

Phylogenetic Reconstruction Based on Synteny Block and Gene Adjacencies.

Gene order can be used as an informative character to reconstruct phylogenetic relationships between species independently from the local information present in gene/protein sequences. PhyChro is a reconstruction method based on chromosomal rearrangements, applicable to a wide range of eukaryotic genomes with different gene contents and levels of synteny conservation. For each synteny breakpoint issued from pairwise genome comparisons, the algorithm defines two disjoint sets of genomes, named partial splits, respectively, supporting the two block adjacencies defining the breakpoint. Considering all partial splits issued from all pairwise comparisons, a distance between two genomes is computed from the number of partial splits separating them. Tree reconstruction is achieved through a bottom-up approach by iteratively grouping sister genomes minimizing genome distances. PhyChro estimates branch lengths based on the number of synteny breakpoints and provides confidence scores for the branches. PhyChro performance is evaluated on two data sets of 13 vertebrates and 21 yeast genomes by using up to 130,000 and 179,000 breakpoints, respectively, a scale of genomic markers that has been out of reach until now. PhyChro reconstructs very accurate tree topologies even at known problematic branching positions. Its robustness has been benchmarked for different synteny block reconstruction methods. On simulated data PhyChro reconstructs phylogenies perfectly in almost all cases, and shows the highest accuracy compared with other existing tools. PhyChro is very fast, reconstructing the vertebrate and yeast phylogenies in <15 min.

Large-scale identification of coevolution signals across homo-oligomeric protein interfaces by direct coupling analysis.

Proteins have evolved to perform diverse cellular functions, from serving as reaction catalysts to coordinating cellular propagation and development. Frequently, proteins do not exert their full potential as monomers but rather undergo concerted interactions as either homo-oligomers or with other proteins as hetero-oligomers. The experimental study of such protein complexes and interactions has been arduous. Theoretical structure prediction methods are an attractive alternative. Here, we investigate homo-oligomeric interfaces by tracing residue coevolution via the global statistical direct coupling analysis (DCA). DCA can accurately infer spatial adjacencies between residues. These adjacencies can be included as constraints in structure prediction techniques to predict high-resolution models. By taking advantage of the ongoing exponential growth of sequence databases, we go significantly beyond anecdotal cases of a few protein families and apply DCA to a systematic large-scale study of nearly 2,000 Pfam protein families with sufficient sequence information and structurally resolved homo-oligomeric interfaces. We find that large interfaces are commonly identified by DCA. We further demonstrate that DCA can differentiate between subfamilies with different binding modes within one large Pfam family. Sequence-derived contact information for the subfamilies proves sufficient to assemble accurate structural models of the diverse protein-oligomers. Thus, we provide an approach to investigate oligomerization for arbitrary protein families leading to structural models complementary to often-difficult experimental methods. Combined with ever more abundant sequential data, we anticipate that this study will be instrumental to allow the structural description of many heteroprotein complexes in the future.

BIS2Analyzer: a server for co-evolution analysis of conserved protein families.

Along protein sequences, co-evolution analysis identifies residue pairs demonstrating either a specific co-adaptation, where changes in one of the residues are compensated by changes in the other during evolution or a less specific external force that affects the evolutionary rates of both residues in a similar magnitude. In both cases, independently of the underlying cause, co-evolutionary signatures within or between proteins serve as markers of physical interactions and/or functional relationships. Depending on the type of protein under study, the set of available homologous sequences may greatly differ in size and amino acid variability. BIS2Analyzer, openly accessible at http://www.lcqb.upmc.fr/BIS2Analyzer/, is a web server providing the online analysis of co-evolving amino-acid pairs in protein alignments, especially designed for vertebrate and viral protein families, which typically display a small number of highly similar sequences. It is based on BIS2, a re-implemented fast version of the co-evolution analysis tool Blocks in Sequences (BIS). BIS2Analyzer provides a rich and interactive graphical interface to ease biological interpretation of the results.

Direct coevolutionary couplings reflect biophysical residue interactions in proteins.

Coevolution of residues in contact imposes strong statistical constraints on the sequence variability between homologous proteins. Direct-Coupling Analysis (DCA), a global statistical inference method, successfully models this variability across homologous protein families to infer structural information about proteins. For each residue pair, DCA infers 21 × 21 matrices describing the coevolutionary coupling for each pair of amino acids (or gaps). To achieve the residue-residue contact prediction, these matrices are mapped onto simple scalar parameters; the full information they contain gets lost. Here, we perform a detailed spectral analysis of the coupling matrices resulting from 70 protein families, to show that they contain quantitative information about the physico-chemical properties of amino-acid interactions. Results for protein families are corroborated by the analysis of synthetic data from lattice-protein models, which emphasizes the critical effect of sampling quality and regularization on the biochemical features of the statistical coupling matrices.

Influence of the single-strand linker composition on the structural/dynamical properties of a truncated octahedral DNA nano-cage family.

The structural/dynamical properties of three truncated octahedral DNA nano-cages composed by identical double helices but single strand linkers with different composition, namely 7 thymidines, 7 adenines, and 7 alternated thymidines and adenines, have been investigated through classical molecular dynamics simulations. Trajectories have been analyzed to investigate the role of the linkers in defining nano-cages stability and flexibility, including possible influence on the internal cages motions. The data indicate that the cages behavior is almost identical and that the structural/dynamical parameters measured along the trajectories are not particularly affected by the presence of different bases. These results demonstrate that the constraints imposed by the nano-structure geometry are the main factor in modulating these properties

The weak, fluctuating, dipole moment of membrane-bound hydrogenase from Aquifex aeolicus accounts for its adaptability to charged electrodes.

[NiFe] hydrogenases from Aquifex aeolicus (AaHase) and Desulfovibrio fructosovorans (DfHase) have been mainly studied to characterize physiological electron transfer processes, or to develop biotechnological devices such as biofuel cells. In this context, it remains difficult to control the orientation of AaHases on electrodes to achieve a fast interfacial electron transfer. Here, we study the electrostatic properties of these two proteins based on microsecond-long molecular dynamics simulations that we compare to voltammetry experiments. Our calculations show weak values and large fluctuations of the dipole direction in AaHase compared to DfHase, enabling the AaHase to absorb on both negatively and positively charged electrodes, with an orientation distribution that induces a spread in electron transfer rates. Moreover, we discuss the role of the transmembrane helix of AaHase and show that it does not substantially impact the general features of the dipole moment.

Simulative and experimental investigation on the cleavage site that generates the soluble human LOX-1.

Lectin-like oxidized low-density lipoprotein receptor-1 (LOX-1) is a scavenger receptor that mediates the recognition, the binding and internalization of ox-LDL. A truncated soluble form of LOX-1 (sLOX-1) has been identified that, at elevated levels, has been associated to acute coronary syndrome. Human sLOX-1 is the extracellular part of membrane LOX-1 which is cleaved in the NECK domain with a mechanism that has not yet been identified. Purification of human sLOX-1 has been carried out to experimentally identify the cleavage site region within the NECK domain. Molecular modelling and classical molecular dynamics simulation techniques have been used to characterize the structural and dynamical properties of the LOX-1 NECK domain in the presence and absence of the CTLD recognition region, taking into account the obtained proteolysis results. The simulative data indicate that the NECK domain is stabilized by the coiled-coil heptad repeat motif along the simulations, shows a definite flexibility pattern and is characterized by specific electrostatic potentials. The detection of a mobile inter-helix space suggests an explanation for the in vivo susceptibility of the NECK domain to the proteolytic cleavage, validating the assumption that the NECK domain sequence is composed of a coiled-coil motif destabilized in specific regions of functional significance.

Characterization of a novel CYP2C9 gene mutation and structural bioinformatic protein analysis in a warfarin hypersensitive patient.

Warfarin (coumadin) is a worldwide-prescribed anticoagulant for the long-term treatment and prevention of thromboembolic events, presenting a great interindividual variability in the required dose. It is known that both environmental and genetic factors influence the dose necessary for the therapeutic effect. Herein we describe a pharmacogenetic study conducted on an Italian patient with warfarin hypersensitivity, who required a very low dosage to achieve therapeutic anticoagulation effect. We genotyped common polymorphisms in VKORC1, CYP2C9, and CYP4F2 genes, known to be involved in warfarin dosing. As the patient resulted in a mixture of low-dosing and high-dosing polymorphic variants, we searched for rare mutations by direct sequencing of the same genes. We identified in the CYP2C9 gene, a novel mutation in heterozygote status, c.374G>T, which produces the Arg125Leu substitution. We have observed, through an electrostatic analysis, that the new mutation produces an electrostatic alteration on the cytochrome surface.

Effect of the charge distribution along the "ferritin-like" pores of the proteins from the Dps family on the iron incorporation process.

DNA-binding proteins from starved cells (Dps) differ in the number and position of charged residues along the "ferritin-like" pores that are used by iron to reach the ferroxidase center and the protein cavity. These differences are shown to affect significantly the electrostatic potential at the pores, which determines the extent of cooperativity in the iron uptake kinetics and thereby the mass distribution of the ferric hydroxide micelles inside the protein cavity. These conclusions are of biotechnological value in the preparation of protein-enclosed nanomaterials and are expected to apply also to ferritins. They were reached after characterization of the Dps from Listeria innocua, Helicobacter pylori, Thermosynechococcus elongatus, Escherichia coli, and Mycobacterium smegmatis. The characterization comprised the calculation of the electrostatic potential at the pores, determination of the iron uptake kinetics in the presence of molecular oxygen or hydrogen peroxide, and analysis of the proteins by means of the sedimentation velocity after iron incorporation.

Structural-dynamical investigation of the ZnuA histidine-rich loop: involvement in zinc management and transport.

Comparative homology modelling techniques have been used to model the protein ZnuA from Salmonella enterica serovar Typhimurium using the 3D structure of the homologous protein from Escherichia coli. These two-domain proteins bind one Zn(2+) atom, with high affinity, in the inter-domain cleft and possess a histidine-rich loop in the N-terminal domain. Alternative structures of the ZnuA histidine-rich loop, never resolved by the X-ray diffraction method, have been modelled. A model of the apo form, one with the histidine-rich loop deleted and two alternative structures with a second zinc ion bound to the histidine-rich loop, have been generated. In all the modelled proteins, investigated through molecular dynamics simulation, the histidine-rich loop is highly mobile and its fluctuations are correlated to the ligand stability observed in the zinc sites. Based on the plasticity of the histidine-rich loop and its significant effects on protein mobility a possible role in the capture and/or transfer of the zinc ions has been suggested.

ADP/ATP mitochondrial carrier MD simulations to shed light on the structural-dynamical events that, after an additional mutation, restore the function in a pathological single mutant.

Molecular dynamics simulations of the wild type bovine ADP/ATP mitochondrial carrier, and of the single Ala113Pro and double Ala113Pro/Val180Met mutants, embedded in a lipid bilayer, have been carried out for 30ns to shed light on the structural-dynamical changes induced by the Val180Met mutation restoring the carrier function in the Ala113Pro pathologic mutant. Principal component analysis indicates that, for the three systems, the protein dynamics is mainly characterized by the motion of the matrix loops and of the odd-numbered helices having a conserved proline in their central region. Analysis of the motions shows a different behaviour of single pathological mutant with respect of the other two systems. The single mutation induces a regularization and rigidity of the H3 helix, lost upon the introduction of the second mutation. This is directly correlated to the salt bridge distribution involving residues Arg79, Asp134 and Arg234, hypothesized to interact with the substrate. In fact, in the wild type simulation two stable inter-helices salt bridges, crucial for substrate binding, are present almost over all the simulation time. In line with the impaired ADP transport, one salt interaction is lost in the single mutant trajectory but reappears in the double mutant simulation, where a salt bridge network matching the wild type is restored. Other important structural-dynamical properties, such as the trans-membrane helices mobility, analyzed via the principal component analysis, are similar for the wild type and double mutant while are different for the single mutant, providing a mechanistic explanation for their different functional properties.

Conjugated eicosapentaenoic acid inhibits human topoisomerase IB with a mechanism different from camptothecin.

Conjugated eicosapentaenoic acid (cEPA) has been found to have antitumor effects which has been ascribed to their ability to inhibit DNA topoisomerases and DNA polymerases. We here show that cEPA inhibits the catalytic activity of human topoisomerase I, but unlike camptothecin it does not stabilize the cleavable complex, indicating a different mechanism of action. cEPA inhibits topoisomerase by impeding the catalytic cleavage of the DNA substrate as demonstrated using specific oligonucleotide substrates, and prevents the stabilization of the cleavable complex by camptothecin. Preincubation of the inhibitor with the enzyme is required to obtain complete inhibition. Molecular docking simulations indicate that the preferred cEPA binding site is proximal to the active site with the carboxylic group strongly interacting with the positively charged K443 and K587. Taken together the results indicate that cEPA inhibitor does not prevent DNA binding but inhibits DNA cleavage, binding in a region close to the topoisomerase active site.

Residues of Alpha Helix H3 Determine Distinctive Features of Transforming Growth Factor ß3.

Transforming growth factors (TGF-ßs) are proteins that regulate cell growth by binding to their receptors. In contrast to transforming growth factor (TGF) ß1, TGF-ß3 homodimer is believed to exist also in an open conformation, in which both of its monomers are loosely packed against each other. At the origin of this difference is the H3-helix. Its sequence and degree of structuration seem to govern the outcome of TGF dimerization. We docked two monomers of TGF-ß3 with intact and altered H3 α-helix against each other using HADDOCK. TGF-ß3 monomer with an intact H3-helix exclusively forms closed conformations of homodimer, whereas the open conformation may coexist with the closed one when a part of the H3 α-helix is destabilized. We quantify the difference in its conformational preference for the open versus the closed structure by calculating the binding energy between monomers using the MMPBSA approach. We compare the wild type (wt) TGFß3/TGFß1 homodimers in the Protein Data Bank to a swapped mutant where all residues of the H3-helix were mutated to the respective TGFß1/TGFß3 sequence. Swapping stabilizes the closed conformation and destabilizes the open conformation of TGFß3. Further detailed insight is derived from molecular dynamics simulation studies suggesting that Val 61 of the H3-helix may act as an anchor residue for the closed conformation of TGFß3. Computational alanine scanning mutagenesis confirms that several residues of the H3-helix are the hot residues for the closed conformation of TGFß3. These observations may bear relevance to general conformational transitions in proteins and specifically in the TGFß superfamily.

Molecular mechanism of statin-mediated LOX-1 inhibition.

Statins are largely used in clinics in the treatment of patients with cardiovascular diseases for their effect on lowering circulating cholesterol. Lectin-like oxidized low-density lipoprotein (LOX-1), the primary receptor for ox-LDL, plays a central role in the pathogenesis of atherosclerosis and cardiovascular disorders. We have recently shown that chronic exposure of cells to lovastatin disrupts LOX-1 receptor cluster distribution in plasma membranes, leading to a marked loss of LOX-1 function. Here we investigated the molecular mechanism of statin-mediated LOX-1 inhibition and we demonstrate that all tested statins are able to displace the binding of fluorescent ox-LDL to LOX-1 by a direct interaction with LOX-1 receptors in a cell-based binding assay. Molecular docking simulations confirm the interaction and indicate that statins completely fill the hydrophobic tunnel that crosses the C-type lectin-like (CTLD) recognition domain of LOX-1. Classical molecular dynamics simulation technique applied to the LOX-1 CTLD, considered in the entire receptor structure with or without a statin ligand inside the tunnel, indicates that the presence of a ligand largely increases the dimer stability. Electrophoretic separation and western blot confirm that different statins binding stabilize the dimer assembly of LOX-1 receptors in vivo. The simulative and experimental results allow us to propose a CTLD clamp motion, which enables the receptor-substrate coupling. These findings reveal a novel and significant functional effect of statins.

Multiscale simulations give insight into the hydrogen in and out pathways of [NiFe]-hydrogenases from Aquifex aeolicus and Desulfovibrio fructosovorans.

[NiFe]-hydrogenases catalyze the cleavage of molecular hydrogen into protons and electrons and represent promising tools for H2-based technologies such as biofuel cells. However, many aspects of these enzymes remain to be understood, in particular how the catalytic center can be protected from irreversible inactivation by O2. In this work, we combined homology modeling, all-atom molecular dynamics, and coarse-grain Brownian dynamics simulations to investigate and compare the dynamic and mechanical properties of two [NiFe]-hydrogenases: the soluble O2-sensitive enzyme from Desulfovibrio fructosovorans, and the O2-tolerant membrane-bound hydrogenase from Aquifex aeolicus. We investigated the diffusion pathways of H2 from the enzyme surface to the central [NiFe] active site, and the possible proton pathways that are used to evacuate hydrogen after the oxidation reaction. Our results highlight common features of the two enzymes, such as a Val/Leu/Arg triad of key residues that controls ligand migration and substrate access in the vicinity of the active site, or the key role played by a Glu residue for proton transfer after hydrogen oxidation. We show specificities of each hydrogenase regarding the enzymes internal tunnel network or the proton transport pathways.

Mapping multiple potential ATP binding sites on the matrix side of the bovine ADP/ATP carrier by the combined use of MD simulation and docking.

The mitochondrial adenosine diphosphate/adenosine triphosphate (ADP/ATP) carrier-AAC-was crystallized in complex with its specific inhibitor carboxyatractyloside (CATR). The protein consists of a six-transmembrane helix bundle that defines the nucleotide translocation pathway, which is closed towards the matrix side due to sharp kinks in the odd-numbered helices. In this paper, we describe the interaction between the matrix side of the AAC transporter and the ATP(4-) molecule using carrier structures obtained through classical molecular dynamics simulation (MD) and a protein-ligand docking procedure. Fifteen structures were extracted from a previously published MD trajectory through clustering analysis, and 50 docking runs were carried out for each carrier conformation, for a total of 750 runs ("MD docking"). The results were compared to those from 750 docking runs performed on the X-ray structure ("X docking"). The docking procedure indicated the presence of a single interaction site in the X-ray structure that was conserved in the structures extracted from the MD trajectory. MD docking showed the presence of a second binding site that was not found in the X docking. The interaction strategy between the AAC transporter and the ATP(4-) molecule was analyzed by investigating the composition and 3D arrangement of the interaction pockets, together with the orientations of the substrate inside them. A relationship between sequence repeats and the ATP(4-) binding sites in the AAC carrier structure is proposed.