Chasing long-range evolutionary couplings in the AlphaFold era

Coevolution between protein residues is normally interpreted as direct contact. However, the evolutionary record of a protein sequence contains rich information that may include long-range functional couplings, couplings that report on homo-oligomeric states or even conformational changes. Due to the complexity of the sequence space and the lack of structural information on various members of a protein family, it has been difficult to effectively mine the additional information encoded in a multiple sequence alignment (MSA). Here, taking advantage of the recent release of the AlphaFold (AF) database we attempt to identify coevolutionary couplings that cannot be explained simply by spatial proximity. We propose a simple computational method that performs direct coupling analysis on a MSA and searches for couplings that are not satisfied in any of the AF models of members of the identified protein family. Application of this method on 2012 protein families suggests that ~12% of the total identified coevolving residue pairs are spatially distant and more likely to be disordered than their contacting counterparts. We expect that this analysis will help improve the quality of coevolutionary distance restraints used for structure determination and will be useful in identifying potentially functional/allosteric cross-talk between distant residues.     

Drosophila acidic ribosomal protein rpA2: Sequence and characterization

A cDNA encoding the Drosophila melanogaster acidic ribosomal protein rpA2 was cloned and sequenced. rpA2 is homologous to the Artemia salina acidic ribosomal protein eL12′. In situ hybridization to salivary gland polytene chromosomes localizes the rpA2 gene to band 21C. It is a single copy gene, with an mRNA of 0.8 kb. Two-dimensional gel electrophoresis of Drosophila ribosomal proteins followed by immuno-blotting showed that the rpA2 protein has an apparent relative mobility in SDS of 17 kD and an isoelectric point less than pH 5.0. Although the Drosophila gene rp21C may be the same as rpA2, the reported sequences differ. Comparisons of the aligned nucleotide sequences coding for the acidic ribosomal proteins rpA1 and rpA2 of Drosophila with those of other eukaryotes support the view of two separate, though closely related, groups of acidic proteins. Comparison with the Artemia homologues suggests that nucleotide identity may have been conserved by some constraint that acts in addition to the requirement for substantial similarity of amino acid sequences. © 1993 Wiley-Liss, Inc.     

Performance of a new family of modular,bed-supported,chromatography devices

Prepacked chromatography columns and cassette filtration units offer many advantages in bioprocessing. These include reduced labor costs and processing times, ease of storage, and enhanced process flexibility. Rectangular formats are particularly attractive as they can be easily stacked and multiplexed together for continuous processing. Cylindrical chromatography beds have dominated bioprocessing even though their bed support and pressure-flow performance vary with bed dimensions. This work presents the performance of novel, rhombohedral chromatography devices with internally supported beds. They are compatible with existing chromatography workstations and can be packed with any standard commercial resin. The devices offer pressure-flow characteristics independent of container-volume, simple multiplexing, and separation performance comparable to cylindrical columns. Their bi-planar, internal bed support allows mechanically less-rigid resins to be used at up to four times higher maximal linear velocities, and productivities approaching 200 g/L/h for affinity resins, compared to the 20 g/L/h typical of many column-based devices. Three 5 L devices should allow processing of up to 3 kg of monoclonal antibody per hour.     

Bioprocess development to produce a hyperthermostable S-methyl-5′-thioadenosine phosphorylase in Escherichia coli

Nucleoside phosphorylases are important biocatalysts for the chemo-enzymatic synthesis of nucleosides and their analogs which are, among others, used for the treatment of viral infections or cancer. S-methyl-5′-thioadenosine phosphorylases (MTAP) are a group of nucleoside phosphorylases and the thermostable MTAP of Aeropyrum pernix (ApMTAP) was described to accept a wide range of modified nucleosides as substrates. Therefore, it is an interesting biocatalyst for the synthesis of nucleoside analogs for industrial and therapeutic applications. To date, thermostable nucleoside phosphorylases were produced in shake flask cultivations using complex media. The drawback of this approach is low volumetric protein yields which hamper the wide-spread application of the thermostable nucleoside phosphorylases in large scale. High cell density (HCD) cultivations allow the production of recombinant proteins with high volumetric yields, as final optical densities >100 can be achieved. Therefore, in this study, we developed a suitable protocol for HCD cultivations of ApMTAP. Initially, optimum expression conditions were determined in 24-well plates using a fed-batch medium. Subsequently, HCD cultivations were performed using E. coli BL21-Gold cells, by employing a glucose-limited fed-batch strategy. Comparing different growth rates in stirred-tank bioreactors, cultivations revealed that growth at maximum growth rates until induction resulted in the highest yields of ApMTAP. On a 500-mL scale, final cell dry weights of 87.1–90.1 g L⁻¹ were observed together with an overproduction of ApMTAP in a 1.9%–3.8% ratio of total protein. Compared to initially applied shake flask cultivations with terrific broth (TB) medium the volumetric yield increased by a factor of 136. After the purification of ApMTAP via heat treatment and affinity chromatography, a purity of more than 90% was determined. Activity testing revealed specific activities in the range of 0.21 ± 0.11 (low growth rate) to 3.99 ± 1.02 U mg⁻¹ (growth at maximum growth rate). Hence, growth at maximum growth rate led to both an increased expression of the target protein and an increased specific enzyme activity. This study paves the way towards the application of thermostable nucleoside phosphorylases in industrial applications due to an improved heterologous expression in Escherichia coli.     

Proteolytic processing of the vitellogenin precursor in the boll weevil,Anthonomus grandis

The soluble proteins of the eggs of the coleopteran insect Anthonomus grandis Boheman, the cotton boll weevil, consist almost entirely of two vitellin types with M_rs of 160,000 and 47,000. We sequenced their N-terminal ends and one internal cyanogen bromide fragment of the large vitellin and compared these sequences with the deduced amino acid sequence from the vitellogenin gene. The results suggest that both the boll weevil vitellin proteins are products of the proteolytic cleavage of a single precursor protein. The smaller 47,000 M vitellin protein is derived from the N-terminal portion of the precursor adjacent to an 18 amino acid signal peptide. The cleavage site between the large and small vitellins at amino acid 362 is adjacent to a pentapeptide sequence containing two pairs of arginine residues. Comparison of the boll weevil sequences with limited known sequences from the single 180,000 M_r honey bee protein show that the honey bee vitellin N-terminal exhibits sequence homology to the N-terminal of the 47,000 M_r boll weevil vitellin. Treatment of the vitellins with an N-glycosidase results in a decrease in molecular weight of both proteins, from 47,000 to 39,000 and from 160,000 to 145,000, indicating that about 10–15% of the molecular weight of each vitellin consists of N-linked carbohydrate. The molecular weight of the deglycosylated large vitellin is smaller than that predicted from the gene sequence, indicating possible further proteolytic processing at the C-terminal of that protein. © 1993 Wiley-Liss, Inc.

1 This article is a US Government work and, as such, is in the public domain in the United States of America.

     

Selection of synthetic proteins to modulate the human frataxin function

Frataxin is a kinetic activator of the mitochondrial supercomplex for iron-sulfur cluster assembly. Low frataxin expression or a decrease in its functionality results in Friedreich's Ataxia (FRDA). With the aim of creating new molecular tools to study this metabolic pathway, and ultimately, to explore new therapeutic strategies, we have investigated the possibility of obtaining small proteins exhibiting a high affinity for frataxin. In this study, we applied the ribosome display approach, using human frataxin as the target. We focused on Affi_224, one of the proteins that we were able to select after five rounds of selection. We have studied the interaction between both proteins and discussed some applications of this specific molecular tutor, concerning the modulation of the supercomplex activity. Affi_224 and frataxin showed a K_D value in the nanomolar range, as judged by surface plasmon resonance analysis. Most likely, it binds to the frataxin acidic ridge, as suggested by the analysis of chemical shift perturbations (nuclear magnetic resonance) and computational simulations. Affi_224 was able to increase Cys NFS1 desulfurase activation exerted by the FRDA frataxin variant G130V. Importantly, Affi_224 interacts with frataxin in a human cellular model. Our results suggest quaternary addition may be a new tool to modulate frataxin function in vivo. Nevertheless, more functional experiments under physiological conditions should be carried out to evaluate Affi_224 effectiveness in FRDA cell models.     

Insights into the Structure of Invisible Conformations of Large Methyl Group Labeled Molecular Machines from High Pressure NMR

Most proteins are highly flexible and can adopt conformations that deviate from the energetically most favorable ground state. Structural information on these lowly populated, alternative conformations is often lacking, despite the functional importance of these states. Here, we study the pathway by which the Dcp1:Dcp2 mRNA decapping complex exchanges between an autoinhibited closed and an open conformation. We make use of methyl Carr–Purcell–Meiboom–Gill (CPMG) NMR relaxation dispersion (RD) experiments that report on the population of the sparsely populated open conformation as well as on the exchange rate between the two conformations. To obtain volumetric information on the open conformation as well as on the transition state structure we made use of RD measurements at elevated pressures. We found that the open Dcp1:Dcp2 conformation has a lower molecular volume than the closed conformation and that the transition state is close in volume to the closed state. In the presence of ATP the volume change upon opening of the complex increases and the volume of the transition state lies in-between the volumes of the closed and open state. These findings show that ATP has an effect on the volume changes that are associated with the opening-closing pathway of the complex. Our results highlight the strength of pressure dependent NMR methods to obtain insights into structural features of protein conformations that are not directly observable. As our work makes use of methyl groups as NMR probes we conclude that the applied methodology is also applicable to high molecular weight complexes.     

Coevolution-Guided Mapping of the Type VI Secretion Membrane Complex-Baseplate Interface

The type VI secretion system (T6SS) is a multiprotein weapon evolved by Gram-negative bacteria to deliver effectors into eukaryotic cells or bacterial rivals. The T6SS uses a contractile mechanism to propel an effector-loaded needle into its target. The contractile tail is built on an assembly platform, the baseplate, which is anchored to a membrane complex. Baseplate-membrane complex interactions are mainly mediated by contacts between the C-terminal domain of the TssK baseplate component and the cytoplasmic domain of the TssL inner membrane protein. Currently, the structural details of this interaction are unknown due to the marginal stability of the TssK-TssL complex. Here we conducted a mutagenesis study based on putative TssK-TssL contact pairs identified by co-evolution analyses. We then evaluated the impact of these mutations on T6SS activity, TssK-TssL interaction and sheath assembly and dynamics in enteroaggregative Escherichia coli. Finally, we probed the TssK-TssL interface by disulfide cross-linking, allowing to propose a model for the baseplate-membrane complex interface.     

The Conserved Yeast Protein Knr4 Involved in Cell Wall Integrity Is a Multi-domain Intrinsically Disordered Protein

Knr4/Smi1 proteins are specific to the fungal kingdom and their deletion in the model yeast Saccharomyces cerevisiae and the human pathogen Candida albicans results in hypersensitivity to specific antifungal agents and a wide range of parietal stresses. In S. cerevisiae, Knr4 is located at the crossroads of several signalling pathways, including the conserved cell wall integrity and calcineurin pathways. Knr4 interacts genetically and physically with several protein members of those pathways. Its sequence suggests that it contains large intrinsically disordered regions. Here, a combination of small-angle X-ray scattering (SAXS) and crystallographic analysis led to a comprehensive structural view of Knr4. This experimental work unambiguously showed that Knr4 comprises two large intrinsically disordered regions flanking a central globular domain whose structure has been established. The structured domain is itself interrupted by a disordered loop. Using the CRISPR/Cas9 genome editing technique, strains expressing KNR4 genes deleted from different domains were constructed. The N-terminal domain and the loop are essential for optimal resistance to cell wall-binding stressors. The C-terminal disordered domain, on the other hand, acts as a negative regulator of this function of Knr4. The identification of molecular recognition features, the possible presence of secondary structure in these disordered domains and the functional importance of the disordered domains revealed here designate these domains as putative interacting spots with partners in either pathway. Targeting these interacting regions is a promising route to the discovery of inhibitory molecules that could increase the susceptibility of pathogens to the antifungals currently in clinical use.     

An empirical energy function for threading protein sequence through the folding motif

In this paper we present a new residue contact potantial derived by statistical analysis of protein crystal structures. This gives mean hydrophobic and pairwise contact energies as a function of residue type and distance interval. To test the accuracy of this potential we generate model structures by “threading” different sequences through backbone folding motifs found in the structural data base. We find that conformational energies calculated by summing contact potentials show perfect specificity in matching the correct sequences with each globular folding motif in a 161-protcin data set. They also identify correct models with the core folding motifs of heme-rythrin and immunoglobulin McPC603 V₁-do- main, among millions of alternatives possible when we align subsequences with α-helices and β-strands, and allow for variation in the lengths of intervening loops. We suggest that contact potentials reflect important constraints on nonbonded interaction in native proteins, and that “threading” may be useful for structure prediction by recognition of folding motif. © 1993 Wiley-Liss, Inc.