Evolutionary diversification of protein-protein interactions by interface add-ons

Effect of Direct and Water-Mediated Interactions on the Identification of Hotspots in Biomolecular Complexes with Multiple Subsystems

Identification of important residues in biochemical complexes is often a crucial step for many problems in molecular biology and biochemistry. A method is proposed to identify hotspots in biomolecular complexes based on a new metric, derived from networks representing molecular subunits (residues, bridging solvent molecules, ligands etc.) connected by interactions. A singular value decomposition of the weighted adjacency matrix is used to construct a scalar rank for each subunit that reflects its importance in the residue interaction network. This metric is called the singular value centrality. In addition, a new formalism is proposed to account for water-mediated interactions in the ranking of residues. Interactions for a residue network can be provided by various computational methods. In this work interactions are obtained from full quantum-mechanical calculations of protein-protein complexes using the fragment molecular orbital method. The ranking results are shown to be in good agreement with earlier computational and experimental studies. The developed method can be used to gain a deeper insight into the role of subunits in complex biomolecular systems.

Validation of aminodeoxychorismate synthase and anthranilate synthase as novel targets for bispecific antibiotics inhibiting conserved protein-protein interactions

Multi-resistant bacteria are a rapidly emerging threat to modern medicine. It is thus essential to identify and validate novel antibacterial targets that promise high robustness against resistance-mediating mutations. This can be achieved by simultaneously targeting several conserved function-determining protein-protein interactions in enzyme complexes from prokaryotic primary metabolism. Here, we selected two evolutionary related glutamine amidotransferase complexes, aminodeoxychorismate synthase and anthranilate synthase, that are required for the biosynthesis of folate and tryptophan in most prokaryotic organisms. Both enzymes rely on the interplay of a glutaminase and a synthase subunit that is conferred by a highly conserved subunit interface. Consequently, inhibiting subunit association in both enzymes by one competing bispecific inhibitor has the potential to suppress bacterial proliferation. We comprehensively verified two conserved interface hot-spot residues as potential inhibitor-binding sites in vitro by demonstrating their crucial role in subunit association and enzymatic activity. For in vivo target validation, we generated genomically modified Escherichia coli strains in which subunit association was disrupted by modifying these central interface residues. The growth of such strains was drastically retarded on liquid and solid minimal medium due to a lack of folate and tryptophan. Remarkably, the bacteriostatic effect was observed even in the presence of heat-inactivated human plasma, demonstrating that accessible host metabolite concentrations do not compensate for the lack of folate and tryptophan within the tested bacterial cells. We conclude that a potential inhibitor targeting both enzyme complexes will be effective against a broad spectrum of pathogens and offer increased resilience against antibiotic resistance.

IMPORTANCE

Antibiotics are indispensable for the treatment of bacterial infections in human and veterinary medicine and are thus a major pillar of modern medicine. However, the exposure of bacteria to antibiotics generates an unintentional selective pressure on bacterial assemblies that over time promotes the development or acquisition of resistance mechanisms, allowing pathogens to escape the treatment. In that manner, humanity is in an ever-lasting race with pathogens to come up with new treatment options before resistances emerge. In general, antibiotics with novel modes of action require more complex pathogen adaptations as compared to chemical derivates of existing entities, thus delaying the emergence of resistance. In this contribution, we use modified Escherichia coli strains to validate two novel targets required for folate and tryptophan biosynthesis that can potentially be targeted by one and the same bispecific protein-protein interaction inhibitor and promise increased robustness against bacterial resistances.

Role of tungsten disulfide quantum dots in specific protein-protein interactions at air-water interface

The intriguing network of antibody-antigen (Ab-Ag) interactions is highly governed by environmental perturbations and the nature of biomolecular interaction. Protein-protein interactions (PPIs) have potential applications in developing protein-adsorption-based sensors and nano-scale materials. Therefore, characterizing PPIs in the presence of a nanomaterial at the molecular level becomes imperative. The present work involves the investigation of antiferritin-ferritin (Ab-Ag) protein interactions under the influence of tungsten disulfide quantum dots (WS2 QDs). Isothermal calorimetry and contact angle measurements validated the strong influence of WS2 QDs on Ab-Ag interactions. The interfacial signatures of nano-bio-interactions were evaluated using sum frequency generation vibration spectroscopy (SFG-VS) at the air-water interface. Our SFG results reveal a variation in the tilt angle of methyl groups by ∼12° ± 2° for the Ab-Ag system in the presence of WS2 QDs. The results illustrated an enhanced ordering of water molecules in the presence of QDs, which underpins the active role of interfacial water molecules during nano-bio-interactions. We have also witnessed a differential impact of QDs on Ab-Ag by raising the concentration of the Ab-Ag combination, which showcased an increased inter-molecular interaction among the Ab and Ag molecules and a minimal influence on the methyl tilt angle. These findings suggest the formation of stronger and ordered Ab-Ag complexes upon introducing WS2 QDs in the aqueous medium and signify the potentiality of WS2 QDs relevant to protein-based sensing assays.

Structural basis for the allosteric pathway of 4-amino-4-deoxychorismate synthase

4-Amino-4-deoxychorismate synthase (ADCS), a chorismate-utilizing enzyme, is composed of two subunits: PabA and PabB. PabA is a glutamine amidotransferase that hydrolyzes glutamine into glutamate and ammonia. PabB is an aminodeoxychorismate synthase that converts chorismate to 4-amino-4-deoxychorismate (ADC) using the ammonia produced by PabA. ADCS functions under allosteric regulation between PabA and PabB. However, the allosteric mechanism remains unresolved because the structure of the PabA-PabB complex has not been determined. Here, the crystal structure and characterization of PapA from Streptomyces venezuelae (SvPapA), a bifunctional enzyme comprising the PabA and PabB domains, is reported. SvPapA forms a unique dimer in which PabA and PabB domains from different monomers complement each other and form an active structure. The chorismate-bound structure revealed that recognition of the C1 carboxyl group by Thr501 and Gly502 of the 498-PIKTG-502 motif in the PabB domain is essential for the catalytic Lys500 to reach the C2 atom, a reaction-initiation site. SvPapA demonstrated ADCS activity in the presence of Mg2+ when glutamate or NH+4 was used as the amino donor. The crystal structure indicated that the Mg2+-binding position changed depending on the binding of chorismate. In addition, significant structural changes were observed in the PabA domain depending on the presence or absence of chorismate. This study provides insights into the structural factors that are involved in the allosteric regulation of ADCS.

Designing Cyclic-Constrained Peptides to Inhibit Human Phosphoglycerate Dehydrogenase

Although loop epitopes at protein-protein binding interfaces often play key roles in mediating oligomer formation and interaction specificity, their binding sites are underexplored as drug targets owing to their high flexibility, relatively few hot spots, and solvent accessibility. Prior attempts to develop molecules that mimic loop epitopes to disrupt protein oligomers have had limited success. In this study, we used structure-based approaches to design and optimize cyclic-constrained peptides based on loop epitopes at the human phosphoglycerate dehydrogenase (PHGDH) dimer interface, which is an obligate homo-dimer with activity strongly dependent on the oligomeric state. The experimental validations showed that these cyclic peptides inhibit PHGDH activity by directly binding to the dimer interface and disrupting the obligate homo-oligomer formation. Our results demonstrate that loop epitope derived cyclic peptides with rationally designed affinity-enhancing substitutions can modulate obligate protein homo-oligomers, which can be used to design peptide inhibitors for other seemingly intractable oligomeric proteins.

Structural and Functional Characterization of the Ureidoacrylate Amidohydrolase RutB from Escherichia coli

We present a detailed structure-function analysis of the ureidoacrylate amidohydrolase RutB from Eschericha coli, which is an essential enzyme of the Rut pathway for pyrimidine utilization. Crystals of selenomethionine-labeled RutB were produced, which allowed us to determine the first structure of the enzyme at a resolution of 1.9 Å and to identify it as a new member of the isochorismatase-like hydrolase family. RutB was co-crystallized with the substrate analogue ureidopropionate, revealing the mode of substrate binding. Mutation of residues constituting the catalytic triad (D24A, D24N, K133A, C166A, C166S, C166T, C166Y) resulted in complete inactivation of RutB, whereas mutation of other residues close to the active site (Y29F, Y35F, N72A, W74A, W74F, E80A, E80D, S92A, S92T, S92Y, Q105A, Y136A, Y136F) leads to distinct changes of the turnover number (k_cat) and/or the Michaelis constant (K_M). The results of our structural and mutational studies allowed us to assign specific functions to individual residues and to formulate a plausible reaction mechanism for RutB.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

Photoswitching of Feedback Inhibition by Tryptophan in Anthranilate Synthase

The artificial regulation of enzymatic activity by light is an important goal of synthetic biology that can be achieved by the incorporation of light-responsive noncanonical amino acids via genetic code expansion. Here, we apply this concept to anthranilate synthase from Salmonella typhimurium (stTrpE). This enzyme catalyzes the first step of tryptophan biosynthesis, and its activity is feedback-inhibited by the binding of the end-product of the pathway to an allosteric site. To put this feedback inhibition of stTrpE by tryptophan under the control of light, we individually replaced 15 different amino acid residues with the photosensitive noncanonical amino acid o-nitrobenzyl-O-tyrosine (ONBY). ONBY contains a sterically demanding caging group that was meant to cover the allosteric site. Steady-state enzyme kinetics showed that the negative effect of tryptophan on the catalytic activity of the two variants stTrpE-K50ONBY and stTrpE-Y455ONBY was diminished compared to the wild-type enzyme by 1 to 2 orders of magnitude. Upon light-induced decaging of ONBY to the less space-consuming tyrosine residue, tryptophan binding to the allosteric site was restored and catalytic activity was inhibited almost as efficiently as observed for wild-type stTrpE. Based on these results, direct photocontrol of feedback inhibition of stTrpE-K50ONBY and stTrpE-Y455ONBY could be achieved by irradiation during the reaction. Molecular modeling studies allowed us to rationalize the observed functional conversion from the noninhibited caged to the tryptophan-inhibited decaged states. Our study shows that feedback inhibition, which is an important mechanism to regulate key metabolic enzymes, can be efficiently controlled by the purposeful use of light-responsive noncanonical amino acids.

Auxiliary interfaces support the evolution of specific toxin-antitoxin pairing

Toxin-antitoxin (TA) systems are a large family of genes implicated in the regulation of bacterial growth and its arrest in response to attacks. These systems encode nonsecreted toxins and antitoxins that specifically pair, even when present in several paralogous copies per genome. Salmonella enterica serovar Typhimurium contains three paralogous TacAT systems that block bacterial translation. We determined the crystal structures of the three TacAT complexes to understand the structural basis of specific TA neutralization and the evolution of such specific pairing. In the present study, we show that alteration of a discrete structural add-on element on the toxin drives specific recognition by their cognate antitoxin underpinning insulation of the three pairs. Similar to other TA families, the region supporting TA-specific pairing is key to neutralization. Our work reveals that additional TA interfaces beside the main neutralization interface increase the safe space for evolution of pairing specificity.

An improved statistical method to identify chemical-genetic interactions by exploiting concentration-dependence

Chemical-genetics (C-G) experiments can be used to identify interactions between inhibitory compounds and bacterial genes, potentially revealing the targets of drugs, or other functionally interacting genes and pathways. C-G experiments involve constructing a library of hypomorphic strains with essential genes that can be knocked-down, treating it with an inhibitory compound, and using high-throughput sequencing to quantify changes in relative abundance of individual mutants. The hypothesis is that, if the target of a drug or other genes in the same pathway are present in the library, such genes will display an excessive fitness defect due to the synergy between the dual stresses of protein depletion and antibiotic exposure. While assays at a single drug concentration are susceptible to noise and can yield false-positive interactions, improved detection can be achieved by requiring that the synergy between gene and drug be concentration-dependent. We present a novel statistical method based on Linear Mixed Models, called CGA-LMM, for analyzing C-G data. The approach is designed to capture the dependence of the abundance of each gene in the hypomorph library on increasing concentrations of drug through slope coefficients. To determine which genes represent candidate interactions, CGA-LMM uses a conservative population-based approach in which genes with negative slopes are considered significant only if they are outliers with respect to the rest of the population (assuming that most genes in the library do not interact with a given inhibitor). We applied the method to analyze 3 independent hypomorph libraries of M. tuberculosis for interactions with antibiotics with anti-tubercular activity, and we identify known target genes or expected interactions for 7 out of 9 drugs where relevant interacting genes are known.

Reprogramming the Specificity of a Protein Interface by Computational and Data-Driven Design

The formation of specific protein complexes in a cell is a non-trivial problem given the co-existence of thousands of different polypeptide chains. A particularly difficult case are two glutamine amidotransferase complexes (anthranilate synthase [AS] and aminodeoxychorismate synthase [ADCS]), which are composed of homologous pairs of synthase and glutaminase subunits. We have attempted to identify discriminating interface residues of the glutaminase subunit TrpG from AS, which are responsible for its specific interaction with the synthase subunit TrpEx and prevent binding to the closely related synthase subunit PabB from ADCS. For this purpose, TrpG-specific interface residues were grafted into the glutaminase subunit PabA from ADCS by two different approaches, namely a computational and a data-driven one. Both approaches resulted in PabA variants that bound TrpEx with higher affinity than PabB. Hence, we have accomplished a reprogramming of protein-protein interaction specificity that provides insights into the evolutionary adaptation of protein interfaces.

Conserved genomic neighborhood is a strong but no perfect indicator for a direct interaction of microbial gene products

BACKGROUND

The order of genes in bacterial genomes is not random; for example, the products of genes belonging to an operon work together in the same pathway. The cotranslational assembly of protein complexes is deemed to conserve genomic neighborhoods even stronger than a common function. This is why a conserved genomic neighborhood can be utilized to predict, whether gene products form protein complexes.

RESULTS

We were interested to assess the performance of a neighborhood-based classifier that analyzes a large number of genomes. Thus, we determined for the genes encoding the subunits of 494 experimentally verified hetero-dimers their local genomic context. In order to generate phylogenetically comprehensive genomic neighborhoods, we utilized the tools offered by the Enzyme Function Initiative. For each subunit, a sequence similarity network was generated and the corresponding genome neighborhood network was analyzed to deduce the most frequent gene product. This was predicted as interaction partner, if its abundance exceeded a threshold, which was the frequency giving rise to the maximal Matthews correlation coefficient. For the threshold of 16%, the true positive rate was 45%, the false positive rate 0.06%, and the precision 55%. For approximately 20% of the subunits, the interaction partner was not found in a neighborhood of ± 10 genes.

CONCLUSIONS

Our phylogenetically comprehensive analysis confirmed that complex formation is a strong evolutionary factor that conserves genome neighborhoods. On the other hand, for 55% of the cases analyzed here, classification failed. Either, the interaction partner was not present in a ± 10 gene window or was not the most frequent gene product.

Asymmetrical diversification of the receptor-ligand interaction controlling self-incompatibility in Arabidopsis

How two-component genetic systems accumulate evolutionary novelty and diversify in the course of evolution is a fundamental problem in evolutionary systems biology. In the Brassicaceae, self-incompatibility (SI) is a spectacular example of a diversified allelic series in which numerous highly diverged receptor-ligand combinations are segregating in natural populations. However, the evolutionary mechanisms by which new SI specificities arise have remained elusive. Using in planta ancestral protein reconstruction, we demonstrate that two allelic variants segregating as distinct receptor-ligand combinations diverged through an asymmetrical process whereby one variant has retained the same recognition specificity as their (now extinct) putative ancestor, while the other has functionally diverged and now represents a novel specificity no longer recognized by the ancestor. Examination of the structural determinants of the shift in binding specificity suggests that qualitative rather than quantitative changes of the interaction are an important source of evolutionary novelty in this highly diversified receptor-ligand system.

Protein-Protein Interactions of Phosphodiesterases

BACKGROUND

Phosphodiesterases (PDEs) are enzymes that play a key role in terminating cyclic nucleotides signalling by catalysing the hydrolysis of 3', 5'- cyclic adenosine monophosphate (cAMP) and/or 3', 5' cyclic guanosine monophosphate (cGMP), the second messengers within the cell that transport the signals produced by extracellular signalling molecules which are unable to get into the cells. However, PDEs are proteins which do not operate alone but in complexes that made up of a many proteins.

OBJECTIVE

This review highlights some of the general characteristics of PDEs and focuses mainly on the Protein-Protein Interactions (PPIs) of selected PDE enzymes. The objective is to review the role of PPIs in the specific mechanism for activation and thereby regulation of certain biological functions of PDEs.

METHODS

The article discusses some of the PPIs of selected PDEs as reported in recent scientific literature. These interactions are critical for understanding the biological role of the target PDE.

RESULTS

The PPIs have shown that each PDE has a specific mechanism for activation and thereby regulation a certain biological function.

CONCLUSION

Targeting of PDEs to specific regions of the cell is based on the interaction with other proteins where each PDE enzyme binds with specific protein(s) via PPIs.

'Democratized' genomic enzymology web tools for functional assignment

The protein databases contain an exponentially growing number of sequences as a result of the recent increase in ease and decrease in cost of genome sequencing. The rate of data accumulation far exceeds the rate of functional studies, producing an increase in genomic 'dark matter', sequences for which no precise and validated function is defined. Publicly accessible, that is 'democratized,' genomic enzymology web tools are essential to leverage the protein and genome databases for discovery of the in vitro activities and in vivo functions of novel enzymes and proteins belonging to the dark matter. In this review, we discuss the use of web tools that have proven successful for functional assignment. We also describe a mechanism for ensuring the capture of published functional data so that the quality of both curated and automated annotations transfer can be improved.

"Pinching" the ammonia tunnel of CTP synthase unveils coordinated catalytic and allosteric-dependent control of ammonia passage

Molecular gates within enzymes often play important roles in synchronizing catalytic events. We explored the role of a gate in cytidine-5'-triphosphate synthase (CTPS) from Escherichia coli. This glutamine amidotransferase catalyzes the biosynthesis of CTP from UTP using either l-glutamine or exogenous NH₃ as a substrate. Glutamine is hydrolyzed in the glutaminase domain, with GTP acting as a positive allosteric effector, and the nascent NH₃ passes through a gate located at the end of a ~25-Å tunnel before entering the synthase domain where CTP is generated. Substitution of the gate residue Val 60 by Ala, Cys, Asp, Trp, or Phe using site-directed mutagenesis and subsequent kinetic analyses revealed that V60-substitution impacts glutaminase activity, nucleotide binding, salt-dependent inhibition, and inter-domain NH₃ transport. Surprisingly, the increase in steric bulk present in V60F perturbed the local structure consistent with "pinching" the tunnel, thereby revealing processes that synchronize the transfer of NH₃ from the glutaminase domain to the synthase domain. V60F had a slightly reduced coupling efficiency at maximal glutaminase activity that was ameliorated by slowing down the glutamine hydrolysis reaction, consistent with a "bottleneck" effect. The inability of V60F to use exogenous NH₃ was overcome in the presence of GTP, and more so if CTPS was covalently modified by 6-diazo-5-oxo-l-norleucine. Use of NH₂OH by V60F as an alternative bulkier substrate occurred most efficiently when it was concomitant with the glutaminase reaction. Thus, the glutaminase activity and GTP-dependent activation act in concert to open the NH₃ gate of CTPS to mediate inter-domain NH₃ transport.

Characteristics of a PHD Finger Subtype

Although the plant homeodomain (PHD) finger superfamily is known for its site-specific readouts of histone tails, the origins of the mechanistic differences in histone H3 readout by different PHD subtypes remain less clear. We show that sequences containing the xCDxCDx motif in the PHD treble clef (xCDxCDx-PHD) constitute a distinct subtype, based on the following observations: (i) the amino acid composition of the binding site is strikingly different from other subtypes due to position-specific enrichment of negatively charged and bulky nonpolar residues, (ii) the binding site positions are mutually and positively correlated, and this correlation is absent in other subtypes, and (iii) there are only small structural deviations, despite low sequence similarity. The xCDxCDx-PHD constitutes ∼20% of the PHD family, and the double PHD fingers (DPFs) are 10% of the total number of xCDxCDx-PHDs. This subtype originated early in the evolution of eukaryotes but has diversified within the metazoan lineage. Despite sequence diversification, the positions of the enriched nonpolar residues, in particular, show very small structural deviations, suggesting critical contributions of nonpolar residues in the binding mechanism of this subtype. Using mutagenesis, we probed the contributions of the binding-site positions enriched in nonpolar residues in four xCDxCDx-PHD proteins and found that they contribute to the tight packing of the H3 residues. This effect may potentially be exploited, as we observed affinity enhancement upon substituting a bulky nonpolar residue at the same binding site in another histone reader. Overall, we present a detailed characterization of PHD subtypes.