Structure-based protein design with deep learning

Since the first revelation of proteins functioning as macromolecular machines through their three dimensional structures, researchers have been intrigued by the marvelous ways the biochemical processes are carried out by proteins. The aspiration to understand protein structures has fueled extensive efforts across different scientific disciplines. In recent years, it has been demonstrated that proteins with new functionality or shapes can be designed via structure-based modeling methods, and the design strategies have combined all available information — but largely piece-by-piece — from sequence derived statistics to the detailed atomic-level modeling of chemical interactions. Despite the significant progress, incorporating data-derived approaches through the use of deep learning methods can be a game changer. In this review, we summarize current progress, compare the arc of developing the deep learning approaches with the conventional methods, and describe the motivation and concepts behind current strategies that may lead to potential future opportunities.     

Dynamic and Reversible Aggregation of the Human CAP Superfamily Member GAPR-1 in Protein Inclusions in Saccharomyces cerevisiae

Many proteins that can assemble into higher order structures termed amyloids can also concentrate into cytoplasmic inclusions via liquid–liquid phase separation. Here, we study the assembly of human Golgi-Associated plant Pathogenesis Related protein 1 (GAPR-1), an amyloidogenic protein of the Cysteine-rich secretory proteins, Antigen 5, and Pathogenesis-related 1 proteins (CAP) protein superfamily, into cytosolic inclusions in Saccharomyces cerevisiae. Overexpression of GAPR-1-GFP results in the formation GAPR-1 oligomers and fluorescent inclusions in yeast cytosol. These cytosolic inclusions are dynamic and reversible organelles that gradually increase during time of overexpression and decrease after promoter shut-off. Inclusion formation is, however, a regulated process that is influenced by factors other than protein expression levels. We identified N-myristoylation of GAPR-1 as an important determinant at early stages of inclusion formation. In addition, mutations in the conserved metal-binding site (His54 and His103) enhanced inclusion formation, suggesting that these residues prevent uncontrolled protein sequestration. In agreement with this, we find that addition of Zn²⁺ metal ions enhances inclusion formation. Furthermore, Zn²⁺ reduces GAPR-1 protein degradation, which indicates stabilization of GAPR-1 in inclusions. We propose that the properties underlying both the amyloidogenic properties and the reversible sequestration of GAPR-1 into inclusions play a role in the biological function of GAPR-1 and other CAP family members.     

A mini-Tn5-derived transposon with reportable and selectable markers enables rapid generation and screening of insertional mutants in Gram-negative bacteria

We re-engineered a classic tool for mutagenesis and gene expression studies in Gram-negative bacteria. Our modified Tn5-based transposon contains multiple features that allow rapid selection for mutants, direct quantification of gene expression and straightforward cloning of the inactivated gene. The promoter-less gfp-km cassette provides selection and reporter assay depending on the activity of the promoter upstream of the transposon insertion site. The cat gene facilitates positive antibiotic selection for mutants, while the narrow R6Kγ replication origin forces transposition in recipient strains lacking the pir gene and enables cloning of the transposon flanked with the disrupted gene from the chromosome. The suicide vector pCKD100, a plasmid that could be delivered into recipient cells through biparental mating or electroporation, harbours the modified transposon. We used the transposon to mutagenize Pectobacterium versatile KD100, Pseudumonas coronafaciens PC27R and Escherichia coli 35150N. The fluorescence intensities of mutants expressing high GFP could be quantified and detected qualitatively. Transformation efficiency from conjugation ranged from 1600 to 1900 CFU per ml. We sequenced the upstream flanking regions, identified the putative truncated genes and demonstrated the restoration of the GFP phenotype through marker exchange. The mini-Tn5 transposon was also utilized to construct mutant a library of P. versatile for forward genetic screens.     

Protistan microbial observatory in the Cariaco Basin,Caribbean. I. Pyrosequencing vs Sanger insights into species richness

Microbial diversity and distribution are topics of intensive research. In two companion papers in this issue, we describe the results of the Cariaco Microbial Observatory (Caribbean Sea, Venezuela). The Basin contains the largest body of marine anoxic water, and presents an opportunity to study protistan communities across biogeochemical gradients. In the first paper, we survey 18S ribosomal RNA (rRNA) gene sequence diversity using both Sanger- and pyrosequencing-based approaches, employing multiple PCR primers, and state-of-the-art statistical analyses to estimate microbial richness missed by the survey. Sampling the Basin at three stations, in two seasons, and at four depths with distinct biogeochemical regimes, we obtained the largest, and arguably the least biased collection of over 6000 nearly full-length protistan rRNA gene sequences from a given oceanographic regime to date, and over 80 000 pyrosequencing tags. These represent all major and many minor protistan taxa, at frequencies globally similar between the two sequence collections. This large data set provided, via the recently developed parametric modeling, the first statistically sound prediction of the total size of protistan richness in a large and varied environment, such as the Cariaco Basin: over 36 000 species, defined as almost full-length 18S rRNA gene sequence clusters sharing over 99% sequence homology. This richness is a small fraction of the grand total of known protists (over 100 000–500 000 species), suggesting a degree of protistan endemism.     

Limited proteolysis patterns of the B chain of insulin.

The oxidized B chain of insulin was used as a simple model for further consideration of limited proteolysis with low substrate:enzyme ratios. With low B chain:trypsin ratios, the ordinarily slower cleavage rate of the -Lys₂₉-Ala₃₀ bond essentially equaled the cleavage saturation rate of the -Arg₂₂-Gly₂₃ bond. This led to the disappearance of octapeptide which ordinarily forms most rapidly. Heptapeptide and alanine, formed mainly by cleavage of the octapeptide, decreased somewhat at high enzyme relative levels. Trypsin added to B chain formed a single chromatographic peak.     

DNA Sequences from Formalin-Fixed Nematodes: Integrating Molecular and Morphological Approaches to Taxonomy

To effectively integrate DNA sequence analysis and classical nematode taxonomy, we must be able to obtain DNA sequences from formalin-fixed specimens. Microdissected sections of nematodes were removed from specimens fixed in formalin, using standard protocols and without destroying morphological features. The fixed sections provided sufficient template for multiple polymerase chain reaction-based DNA sequence analyses.     

Trichoderma reesei sequences that bind to the nuclear matrix enhance transformation frequency

Three DNA fragments, trs1, 2 and 3, were isolated from the Trichoderma reesei genome on the basis of their ability to promote autonomous replication of plasmids in Saccharomyces cerevisiae. Each trs element bound specifically to the isolated T. reesei nuclear matrix in vitro, and two of them bound in vivo, indicating that they are matrix attachment regions (MARs). A similar sequence previously isolated from Aspergillus nidulans (ans1) was also shown to bind specifically to the T. reesei nuclear matrix in vitro. The T. reesei MARs are AT-rich sequences containing 70%, 86% and 73% A+T over 2.9, 0.8 and 3.7 kb, respectively for trs1, 2 and 3. They exhibited no significant sequence homology, but were shown to contain a number of sequence motifs that occur frequently in many MARs identified in other eukaryotes. However, these motifs occurred as frequently in the trs elements as in randomly generated sequences with the same A+T content. trs1 and 3 were shown to be present as single copies in the T. reesei genome. The presence of the trs elements in transforming plasmids enhanced the frequency of integrative transformation of T. reesei up to five fold over plasmids without a trs. No evidence was obtained to suggest that the trs elements promoted efficient replication of plasmids in T. reseei. A mechanism for the enhancement of transformation frequency by the trs elements is proposed. Received: 1 March 1997 / Accepted: 13 May 1997     

Dissecting the Conformational Dynamics of the Bile Acid Transporter Homologue ASBTNM

Apical sodium-dependent bile acid transporter (ASBT) catalyses uphill transport of bile acids using the electrochemical gradient of Na⁺ as the driving force. The crystal structures of two bacterial homologues ASBT_NM and ASBT_Yf have previously been determined, with the former showing an inward-facing conformation, and the latter adopting an outward-facing conformation accomplished by the substitution of the critical Na⁺-binding residue glutamate-254 with an alanine residue. While the two crystal structures suggested an elevator-like movement to afford alternating access to the substrate binding site, the mechanistic role of Na⁺ and substrate in the conformational isomerization remains unclear. In this study, we utilized site-directed alkylation monitored by in-gel fluorescence (SDAF) to probe the solvent accessibility of the residues lining the substrate permeation pathway of ASBT_NM under different Na⁺ and substrate conditions, and interpreted the conformational states inferred from the crystal structures. Unexpectedly, the crosslinking experiments demonstrated that ASBT_NM is a monomer protein, unlike the other elevator-type transporters, usually forming a homodimer or a homotrimer. The conformational dynamics observed by the biochemical experiments were further validated using DEER measuring the distance between the spin-labelled pairs. Our results revealed that Na⁺ ions shift the conformational equilibrium of ASBT_NM toward the inward-facing state thereby facilitating cytoplasmic uptake of substrate. The current findings provide a novel perspective on the conformational equilibrium of secondary active transporters.