Institute for Bioscience and Biotechnology Research

This article is an introduction to the special issue of the journal PROTEINS, dedicated to the tenth Critical Assessment of Structure Prediction (CASP) experiment to assess the state of the art in protein structure modeling. The article describes the conduct of the experiment, the categories of prediction included, and outlines the evaluation and assessment procedures. The 10 CASP experiments span almost 20 years of progress in the field of protein structure modeling, and there have been enormous advances in methods and model accuracy in that period. Notable in this round is the first sustained improvement of models with refinement methods, using molecular dynamics. For the first time, we tested the ability of modeling methods to make use of sparse experimental three-dimensional contact information, such as may be obtained from new experimental techniques, with encouraging results. On the other hand, new contact prediction methods, though holding considerable promise, have yet to make an impact in CASP testing. The nature of CASP targets has been changing in recent CASPs, reflecting shifts in experimental structural biology, with more irregular structures, more multi-domain and multi-subunit structures, and less standard versions of known folds. When allowance is made for these factors, we continue to see steady progress in the overall accuracy of models, particularly resulting from improvement of non-template regions.

-

CASP (critical assessment of structure prediction) assesses the state of the art in modeling protein structure from amino acid sequence. The most recent experiment (CASP13 held in 2018) saw dramatic progress in structure modeling without use of structural templates (historically "ab initio" modeling). Progress was driven by the successful application of deep learning techniques to predict inter-residue distances. In turn, these results drove dramatic improvements in three-dimensional structure accuracy: With the proviso that there are an adequate number of sequences known for the protein family, the new methods essentially solve the long-standing problem of predicting the fold topology of monomeric proteins. Further, the number of sequences required in the alignment has fallen substantially. There is also substantial improvement in the accuracy of template-based models. Other areas-model refinement, accuracy estimation, and the structure of protein assemblies-have again yielded interesting results. CASP13 placed increased emphasis on the use of sparse data together with modeling and chemical crosslinking, SAXS, and NMR all yielded more mature results. This paper summarizes the key outcomes of CASP13. The special issue of PROTEINS contains papers describing the CASP13 assessments in each modeling category and contributions from the participants.

Critical assessment of structure prediction (CASP) is a community experiment to advance methods of computing three-dimensional protein structure from amino acid sequence. Core components are rigorous blind testing of methods and evaluation of the results by independent assessors. In the most recent experiment (CASP14), deep-learning methods from one research group consistently delivered computed structures rivaling the corresponding experimental ones in accuracy. In this sense, the results represent a solution to the classical protein-folding problem, at least for single proteins. The models have already been shown to be capable of providing solutions for problematic crystal structures, and there are broad implications for the rest of structural biology. Other research groups also substantially improved performance. Here, we describe these results and outline some of the many implications. Other related areas of CASP, including modeling of protein complexes, structure refinement, estimation of model accuracy, and prediction of inter-residue contacts and distances, are also described.

High-resolution experimental structural determination of protein-protein interactions has led to valuable mechanistic insights, yet due to the massive number of interactions and experimental limitations there is a need for computational methods that can accurately model their structures. Here we explore the use of the recently developed deep learning method, AlphaFold, to predict structures of protein complexes from sequence. With a benchmark of 152 diverse heterodimeric protein complexes, multiple implementations and parameters of AlphaFold were tested for accuracy. Remarkably, many cases (43%) had near-native models (medium or high critical assessment of predicted interactions accuracy) generated as top-ranked predictions by AlphaFold, greatly surpassing the performance of unbound protein-protein docking (9% success rate for near-native top-ranked models), however AlphaFold modeling of antibody-antigen complexes within our set was unsuccessful. We identified sequence and structural features associated with lack of AlphaFold success, and we also investigated the impact of multiple sequence alignment input. Benchmarking of a multimer-optimized version of AlphaFold (AlphaFold-Multimer) with a set of recently released antibody-antigen structures confirmed a low rate of success for antibody-antigen complexes (11% success), and we found that T cell receptor-antigen complexes are likewise not accurately modeled by that algorithm, showing that adaptive immune recognition poses a challenge for the current AlphaFold algorithm and model. Overall, our study demonstrates that end-to-end deep learning can accurately model many transient protein complexes, and highlights areas of improvement for future developments to reliably model any protein-protein interaction of interest.

This article reports the outcome of the 12th round of Critical Assessment of Structure Prediction (CASP12), held in 2016. CASP is a community experiment to determine the state of the art in modeling protein structure from amino acid sequence. Participants are provided sequence information and in turn provide protein structure models and related information. Analysis of the submitted structures by independent assessors provides a comprehensive picture of the capabilities of current methods, and allows progress to be identified. This was again an exciting round of CASP, with significant advances in 4 areas: (i) The use of new methods for predicting three-dimensional contacts led to a two-fold improvement in contact accuracy. (ii) As a consequence, model accuracy for proteins where no template was available improved dramatically. (iii) Models based on a structural template showed overall improvement in accuracy. (iv) Methods for estimating the accuracy of a model continued to improve. CASP continued to develop new areas: (i) Assessing methods for building quaternary structure models, including an expansion of the collaboration between CASP and CAPRI. (ii) Modeling with the aid of experimental data was extended to include SAXS data, as well as again using chemical cross-linking information. (iii) A team of assessors evaluated the suitability of models for a range of applications, including mutation interpretation, analysis of ligand binding properties, and identification of interfaces. This article describes the experiment and summarizes the results. The rest of this special issue of PROTEINS contains papers describing CASP12 results and assessments in more detail.

Genome editing promises giant leaps forward in advancing biotechnology, agriculture, and basic research. The process relies on the use of sequence specific nucleases (SSNs) to make DNA double stranded breaks at user defined genomic loci, which are subsequently repaired by two main DNA repair pathways: non-homologous end joining (NHEJ) and homology directed repair (HDR). NHEJ can result in frameshift mutations that often create genetic knockouts. These knockout lines are useful for functional and reverse genetic studies but also have applications in agriculture. HDR has a variety of applications as it can be used for gene replacement, gene stacking, and for creating various fusion proteins. In recent years, transcription activator-like effector nucleases and clustered regularly interspaced palindromic repeats (CRISPR) and CRISPR associated protein 9 or CRISPR from Prevotella and Francisella 1 have emerged as the preferred SSNs for research purposes. Here, we review their applications in plant research, discuss current limitations, and predict future research directions in plant genome editing.

ABSTRACT The cell wall of Mycobacterium tuberculosis is composed of unique lipids that are important for pathogenesis. Indeed, the first-ever genetic screen in M. tuberculosis identified genes involved in the biosynthesis and transport of the cell wall lipid PDIM ( p hthiocerol dim ycocerosates) as crucial for the survival of M. tuberculosis in mice. Here we show evidence for a novel molecular mechanism of the PDIM-mediated virulence in M. tuberculosis . We characterized the DNA interaction and the regulon of Rv3167c, a transcriptional repressor that is involved in virulence regulation of M. tuberculosis , and discovered that it controls the PDIM operon. A loss-of-function genetic approach showed that PDIM levels directly correlate with the capacity of M. tuberculosis to escape the phagosome and induce host cell necrosis and macroautophagy. In conclusion, our study attributes a novel role of the cell wall lipid PDIM in intracellular host cell modulation, which is important for host cell exit and dissemination of M. tuberculosis . IMPORTANCE Mycobacterium tuberculosis is a major human pathogen that has coevolved with its host for thousands of years. The complex and unique cell wall of M. tuberculosis contains the lipid PDIM ( p hthiocerol dim ycocerosates), which is crucial for virulence of the bacterium, but its function is not well understood. Here we show that PDIM expression by M. tuberculosis is negatively regulated by a novel transcriptional repressor, Rv3167c. In addition, we discovered that the escape of M. tuberculosis from its intracellular vacuole was greatly augmented by the presence of PDIM. The increased release of M. tuberculosis into the cytosol led to increased host cell necrosis. The discovery of a link between the cell wall lipid PDIM and a major pathogenesis pathway of M. tuberculosis provides important insights into the molecular mechanisms of host cell manipulation by M. tuberculosis .

Modeling of protein structure from amino acid sequence now plays a major role in structural biology. Here we report new developments and progress from the CASP11 community experiment, assessing the state of the art in structure modeling. Notable points include the following: (1) New methods for predicting three dimensional contacts resulted in a few spectacular template free models in this CASP, whereas models based on sequence homology to proteins with experimental structure continue to be the most accurate. (2) Refinement of initial protein models, primarily using molecular dynamics related approaches, has now advanced to the point where the best methods can consistently (though slightly) improve nearly all models. (3) The use of relatively sparse NMR constraints dramatically improves the accuracy of models, and another type of sparse data, chemical crosslinking, introduced in this CASP, also shows promise for producing better models. (4) A new emphasis on modeling protein complexes, in collaboration with CAPRI, has produced interesting results, but also shows the need for more focus on this area. (5) Methods for estimating the accuracy of models have advanced to the point where they are of considerable practical use. (6) A first assessment demonstrates that models can sometimes successfully address biological questions that motivate experimental structure determination. (7) There is continuing progress in accuracy of modeling regions of structure not directly available by comparative modeling, while there is marginal or no progress in some other areas. Proteins 2016; 84(Suppl 1):4-14. © 2016 Wiley Periodicals, Inc.

CRISPR/Cas9 and Cas12a (Cpf1) nucleases are two of the most powerful genome editing tools in plants. In this work, we compared their activities by targeting maize glossy2 gene coding region that has overlapping sequences recognized by both nucleases. We introduced constructs carrying SpCas9-guide RNA (gRNA) and LbCas12a-CRISPR RNA (crRNA) into maize inbred B104 embryos using Agrobacterium-mediated transformation. On-target mutation analysis showed that 90%-100% of the Cas9-edited T0 plants carried indel mutations and 63%-77% of them were homozygous or biallelic mutants. In contrast, 0%-60% of Cas12a-edited T0 plants had on-target mutations. We then conducted CIRCLE-seq analysis to identify genome-wide potential off-target sites for Cas9. A total of 18 and 67 potential off-targets were identified for the two gRNAs, respectively, with an average of five mismatches compared to the target sites. Sequencing analysis of a selected subset of the off-target sites revealed no detectable level of mutations in the T1 plants, which constitutively express Cas9 nuclease and gRNAs. In conclusion, our results suggest that the CRISPR/Cas9 system used in this study is highly efficient and specific for genome editing in maize, while CRISPR/Cas12a needs further optimization for improved editing efficiency.

Many models of evolution calculate the rate of evolution by multiplying the rate at which new mutations originate within a population by a probability of fixation. Here we review the historical origins, contemporary applications, and evolutionary implications of these "origin-fixation" models, which are widely used in evolutionary genetics, molecular evolution, and phylogenetics. Origin-fixation models were first introduced in 1969, in association with an emerging view of "molecular" evolution. Early origin-fixation models were used to calculate an instantaneous rate of evolution across a large number of independently evolving loci; in the 1980s and 1990s, a second wave of origin-fixation models emerged to address a sequence of fixation events at a single locus. Although origin fixation models have been applied to a broad array of problems in contemporary evolutionary research, their rise in popularity has not been accompanied by an increased appreciation of their restrictive assumptions or their distinctive implications. We argue that origin-fixation models constitute a coherent theory of mutation-limited evolution that contrasts sharply with theories of evolution that rely on the presence of standing genetic variation. A major unsolved question in evolutionary biology is the degree to which these models provide an accurate approximation of evolution in natural populations.

In bacterial biofilms, high molecular weight, secreted exopolysaccharides can serve as a scaffold to which additional carbohydrates, proteins, lipids, and nucleic acids adhere, forming the matrix of the developing biofilm. Here we report methods to extract and purify high molecular weight (>15 kDa) exopolysaccharides from biofilms of eight human pathogens, including species of Staphylcococcus, Klebsiella, Acinetobacter, Pseudomonas, and a toxigenic strain of Escherichia coli O157:H7. Glycosyl composition analysis indicated a high total mannose content across all strains with P. aeruginosa and A. baumannii exopolysaccharides comprised of 80-90% mannose, K. pneumoniae and S. epidermidis strains containing 40-50% mannose, and E. coli with ∼10% mannose. Galactose and glucose were also present in all eight strains, usually as the second and third most abundant carbohydrates. N-acetyl-glucosamine and galacturonic acid were found in 6 of 8 strains, while arabinose, fucose, rhamnose, and xylose were found in 5 of 8 strains. For linkage analysis, 33 distinct residue-linkage combinations were detected with the most abundant being mannose-linked moieties, in line with the composition analysis. The exopolysaccharides of two P. aeruginosa strains analyzed were consistent with the Psl carbohydrate, but not Pel or alginate. The S. epidermidis strain had a composition rich in mannose and glucose, which is consistent with the previously described slime associated antigen (SAA) and the extracellular slime substance (ESS), respectively, but no polysaccharide intracellular adhesion (PIA) was detected. The high molecular weight exopolysaccharides from E. coli, K. pneumoniae, and A. baumannii appear to be novel, based on composition and/or ratio analysis of carbohydrates.

Significance The exact function of gasdermin-B, a protein involved in epithelial cell development, is unknown. We provide insights into gasdermin-B function and how it may contribute to cancer progression and genetic susceptibility to asthma and inflammatory bowel disease (IBD). In contrast to other gasdermins, which bind phosphoinositides and cardiolipin only upon cleavage between their N- and C-terminal domains, intact gasdermin-B binds phosphoinositides and, uniquely, sulfatide, a component of the apical membrane of epithelial cells. Polymorphism residues in the C-terminal domain, associated with asthma and IBD, induce structural changes that may affect protein activity. Components of the apical plasma membrane maintain the cell barrier integrity; thus, aberrant sulfatide levels due to changes in the cellular gasdermin-B concentration or activity could affect disease risk.

This article is an introduction to the special issue of the journal PROTEINS, dedicated to the ninth Critical Assessment of Structure Prediction (CASP) experiment to assess the state of the art in protein structure modeling. The article describes the conduct of the experiment, the categories of prediction included, and outlines the evaluation and assessment procedures. Methods for modeling protein structure continue to advance, although at a more modest pace than in the early CASP experiments. CASP developments of note are indications of improvement in model accuracy for some classes of target, an improved ability to choose the most accurate of a set of generated models, and evidence of improvement in accuracy for short "new fold" models. In addition, a new analysis of regions of models not derivable from the most obvious template structure has revealed better performance than expected.

Sphingolipids comprise a major class of structural materials and lipid signaling molecules in all eukaryotic cells. Over the past two decades, there has been a phenomenal growth in the study of sphingolipids (i.e., sphingobiology) at an average rate of ∼1000 research articles per year. Sphingolipid studies in plants, though accounting for only a small fraction (∼6%) of the total number of publications, have also enjoyed proportionally rapid growth in the past decade. Concomitant with the growth of sphingobiology, there has also been tremendous progress in our understanding of the molecular mechanisms of plant innate immunity. In this review, we (i) cross examine and analyze the major findings that establish and strengthen the intimate connections between sphingolipid metabolism and plant programmed cell death (PCD) associated with plant defense or disease; (ii) highlight and compare key bioactive sphingolipids involved in the regulation of plant PCD and possibly defense; (iii) discuss the potential role of sphingolipids in polarized membrane/protein trafficking and formation of lipid rafts as subdomains of cell membranes in relation to plant defense; and (iv) where possible, attempt to identify potential parallels for immunity-related mechanisms involving sphingolipids across kingdoms.

MicroRNAs (miRNAs) are small non-coding RNAs that play important roles in plant development and stress responses. Loss-of-function analysis of miRNA genes has been traditionally challenging due to lack of appropriate knockout tools. In this study, single miRNA genes (OsMIR408 and OsMIR528) and miRNA gene families (miR815a/b/c and miR820a/b/c) in rice were targeted by CRISPR-Cas9. We showed single strand conformation polymorphism (SSCP) is a more reliable method than restriction fragment length polymorphism (RFLP) for identifying CRISPR-Cas9 generated mutants. Frequencies of targeted mutagenesis among regenerated T0 lines ranged from 48% to 89% at all tested miRNA target sites. In the case of miRNA528, three independent guide RNAs (gRNAs) all generated biallelic mutations among confirmed mutant lines. When targeted by two gRNAs, miRNA genes were readily to be deleted at a frequency up to 60% in T0 rice lines. Thus, we demonstrate CRISPR-Cas9 is an effective tool for knocking out plant miRNAs. Single-base pair (bp) insertion/deletion mutations (indels) in mature miRNA regions can lead to the generation of functionally redundant miRNAs. Large deletions at either the mature miRNA or the complementary miRNA* were found to readily abolish miRNA function. Utilizing mutants of OsMIR408 and OsMIR528, we find that knocking out a single miRNA can result in expression profile changes of many other seemingly unrelated miRNAs. In a case study on OsMIR528, we reveal it is a positive regulator in salt stress. Our work not only provides empirical guidelines on targeting miRNAs with CRISPR-Cas9, but also brings new insights into miRNA function and complex cross-regulation in rice.

De novo protein design has been successful in expanding the natural protein repertoire. However, most de novo proteins lack biological function, presenting a major methodological challenge. In vaccinology, the induction of precise antibody responses remains a cornerstone for next-generation vaccines. Here, we present a protein design algorithm called TopoBuilder, with which we engineered epitope-focused immunogens displaying complex structural motifs. In both mice and nonhuman primates, cocktails of three de novo-designed immunogens induced robust neutralizing responses against the respiratory syncytial virus. Furthermore, the immunogens refocused preexisting antibody responses toward defined neutralization epitopes. Overall, our design approach opens the possibility of targeting specific epitopes for the development of vaccines and therapeutic antibodies and, more generally, will be applicable to the design of de novo proteins displaying complex functional motifs.

Colloidal shear thickening presents a significant challenge because the macroscopic rheology becomes increasingly controlled by the microscopic details of short ranged particle interactions in the shear thickening regime. Our measurements here of the first normal stress difference over a wide range of particle volume fractions elucidate the relative contributions from hydrodynamic lubrication and frictional contact forces, which have been debated. At moderate volume fractions we find N_{1}<0, consistent with hydrodynamic models; however, at higher volume fractions and shear stresses these models break down and we instead observe dilation (N_{1}>0), indicating frictional contact networks. Remarkably, there is no signature of this transition in the viscosity; instead, this change in the sign of N_{1} occurs while the shear thickening remains continuous. These results suggest a scenario where shear thickening is driven primarily by the formation of frictional contacts, with hydrodynamic forces playing a supporting role at lower concentrations. Motivated by this picture, we introduce a simple model that combines these frictional and hydrodynamic contributions and accurately fits the measured viscosity over a wide range of particle volume fractions and shear stress.

CRISPR-Cas12a is a promising genome editing system for targeting AT-rich genomic regions. Comprehensive genome engineering requires simultaneous targeting of multiple genes at defined locations. Here, to expand the targeting scope of Cas12a, we screen nine Cas12a orthologs that have not been demonstrated in plants, and identify six, ErCas12a, Lb5Cas12a, BsCas12a, Mb2Cas12a, TsCas12a and MbCas12a, that possess high editing activity in rice. Among them, Mb2Cas12a stands out with high editing efficiency and tolerance to low temperature. An engineered Mb2Cas12a-RVRR variant enables editing with more relaxed PAM requirements in rice, yielding two times higher genome coverage than the wild type SpCas9. To enable large-scale genome engineering, we compare 12 multiplexed Cas12a systems and identify a potent system that exhibits nearly 100% biallelic editing efficiency with the ability to target as many as 16 sites in rice. This is the highest level of multiplex edits in plants to date using Cas12a. Two compact single transcript unit CRISPR-Cas12a interference systems are also developed for multi-gene repression in rice and Arabidopsis. This study greatly expands the targeting scope of Cas12a for crop genome engineering.

There has been a long-standing controversy regarding the effect of chemical denaturants on the dimensions of unfolded and intrinsically disordered proteins: A wide range of experimental techniques suggest that polypeptide chains expand with increasing denaturant concentration, but several studies using small-angle X-ray scattering (SAXS) have reported no such increase of the radius of gyration (Rg). This inconsistency challenges our current understanding of the mechanism of chemical denaturants, which are widely employed to investigate protein folding and stability. Here, we use a combination of single-molecule Förster resonance energy transfer (FRET), SAXS, dynamic light scattering (DLS), and two-focus fluorescence correlation spectroscopy (2f-FCS) to characterize the denaturant dependence of the unfolded state of the spectrin domain R17 and the intrinsically disordered protein ACTR in two different denaturants. Standard analysis of the primary data clearly indicates an expansion of the unfolded state with increasing denaturant concentration irrespective of the protein, denaturant, or experimental method used. This is the first case in which SAXS and FRET have yielded even qualitatively consistent results regarding expansion in denaturant when applied to the same proteins. To more directly illustrate this self-consistency, we used both SAXS and FRET data in a Bayesian procedure to refine structural ensembles representative of the observed unfolded state. This analysis demonstrates that both of these experimental probes are compatible with a common ensemble of protein configurations for each denaturant concentration. Furthermore, the resulting ensembles reproduce the trend of increasing hydrodynamic radius with denaturant concentration obtained by 2f-FCS and DLS. We were thus able to reconcile the results from all four experimental techniques quantitatively, to obtain a comprehensive structural picture of denaturant-induced unfolded state expansion, and to identify the most likely sources of earlier discrepancies.