PredyFlexy: flexibility and local structure prediction from sequence.

Protein structures are necessary for understanding protein function at a molecular level. Dynamics and flexibility of protein structures are also key elements of protein function. So, we have proposed to look at protein flexibility using novel methods: (i) using a structural alphabet and (ii) combining classical X-ray B-factor data and molecular dynamics simulations. First, we established a library composed of structural prototypes (LSPs) to describe protein structure by a limited set of recurring local structures. We developed a prediction method that proposes structural candidates in terms of LSPs and predict protein flexibility along a given sequence. Second, we examine flexibility according to two different descriptors: X-ray B-factors considered as good indicators of flexibility and the root mean square fluctuations, based on molecular dynamics simulations. We then define three flexibility classes and propose a method based on the LSP prediction method for predicting flexibility along the sequence. This method does not resort to sophisticate learning of flexibility but predicts flexibility from average flexibility of predicted local structures. The method is implemented in PredyFlexy web server. Results are similar to those obtained with the most recent, cutting-edge methods based on direct learning of flexibility data conducted with sophisticated algorithms. PredyFlexy can be accessed at http://www.dsimb.inserm.fr/dsimb_tools/predyflexy/.

Knowledge graph-based recommendation framework identifies drivers of resistance in EGFR mutant non-small cell lung cancer.

Resistance to EGFR inhibitors (EGFRi) presents a major obstacle in treating non-small cell lung cancer (NSCLC). One of the most exciting new ways to find potential resistance markers involves running functional genetic screens, such as CRISPR, followed by manual triage of significantly enriched genes. This triage process to identify 'high value' hits resulting from the CRISPR screen involves manual curation that requires specialized knowledge and can take even experts several months to comprehensively complete. To find key drivers of resistance faster we build a recommendation system on top of a heterogeneous biomedical knowledge graph integrating pre-clinical, clinical, and literature evidence. The recommender system ranks genes based on trade-offs between diverse types of evidence linking them to potential mechanisms of EGFRi resistance. This unbiased approach identifies 57 resistance markers from >3,000 genes, reducing hit identification time from months to minutes. In addition to reproducing known resistance markers, our method identifies previously unexplored resistance mechanisms that we prospectively validate.

Predicting protein flexibility through the prediction of local structures.

Protein structures are valuable tools for understanding protein function. However, protein dynamics is also considered a key element in protein function. Therefore, in addition to structural analysis, fully understanding protein function at the molecular level now requires accounting for flexibility. However, experimental techniques that produce both types of information simultaneously are still limited. Prediction approaches are useful alternative tools for obtaining otherwise unavailable data. It has been shown that protein structure can be described by a limited set of recurring local structures. In this context, we previously established a library composed of 120 overlapping long structural prototypes (LSPs) representing fragments of 11 residues in length and covering all known local protein structures. On the basis of the close sequence-structure relationship observed in LSPs, we developed a novel prediction method that proposes structural candidates in terms of LSPs along a given sequence. The prediction accuracy rate was high given the number of structural classes. In this study, we use this methodology to predict protein flexibility. We first examine flexibility according to two different descriptors, the B-factor and root mean square fluctuations from molecular dynamics simulations. We then show the relevance of using both descriptors together. We define three flexibility classes and propose a method based on the LSP prediction method for predicting flexibility along the sequence. The prediction rate reaches 49.6%. This method competes rather efficiently with the most recent, cutting-edge methods based on true flexibility data learning with sophisticated algorithms. Accordingly, flexibility information should be taken into account in structural prediction assessments.

Modelling human liver fibrosis in the context of non-alcoholic steatohepatitis using a microphysiological system.

Non-alcoholic steatohepatitis (NASH) is a common form of chronic liver disease characterised by lipid accumulation, infiltration of immune cells, hepatocellular ballooning, collagen deposition and liver fibrosis. There is a high unmet need to develop treatments for NASH. We have investigated how liver fibrosis and features of advanced clinical disease can be modelled using an in vitro microphysiological system (MPS). The NASH MPS model comprises a co-culture of primary human liver cells, which were cultured in a variety of conditions including+/- excess sugar, fat, exogenous TGFß or LPS. The transcriptomic, inflammatory and fibrotic phenotype of the model was characterised and compared using a system biology approach to identify conditions that mimic more advanced clinical disease. The transcriptomic profile of the model was shown to closely correlate with the profile of patient samples and the model displayed a quantifiable fibrotic phenotype. The effects of Obeticholic acid and Elafibranor, were evaluated in the model, as wells as the effects of dietary intervention, with all able to significantly reduce inflammatory and fibrosis markers. Overall, we demonstrate how the MPS NASH model can be used to model different aspects of clinical NASH but importantly demonstrate its ability to model advanced disease with a quantifiable fibrosis phenotype.

Influence of assignment on the prediction of transmembrane helices in protein structures.

α-Helical transmembrane proteins (TMPα) are composed of a series of helices embedded in the lipid bilayer. Due to technical difficulties, few 3D structures are available. Therefore, the design of structural models of TMPα is of major interest. We study the secondary structures of TMPα by analyzing the influence of secondary structures assignment methods (SSAMs). For this purpose, a published and updated benchmark databank of TMPα is used and several SSAMs (9) are evaluated. The analysis of the results points to significant differences in SSA depending on the methods used. Pairwise comparisons between SSAMs led to more than 10% of disagreement. Helical regions corresponding to transmembrane zones are often correctly characterized. The study of the sequence-structure relationship shows very limited differences with regard to the structural disagreement. Secondary structure prediction based on Bayes' rule and using only a single sequence give correct prediction rates ranging from 78 to 81%. A structural alphabet approach gives a slightly better prediction, i.e., only 2% less than the best equivalent approach, whereas the prediction rate with a very different assignment bypasses 86%. This last result highlights the importance of the correct assignment choice to evaluate the prediction assessment.

Delivering Robust Candidates to the Drug Pipeline through Computational Analysis of Arrayed CRISPR Screens.

Genome-wide arrayed CRISPR screening is a powerful method for drug target identification as it enables exploration of the effect of individual gene perturbations using diverse highly multiplexed functional and phenotypic assays. Using high-content imaging, we can measure changes in biomarker expression, intracellular localization, and cell morphology. Here we present the computational pipeline we have developed to support the analysis and interpretation of arrayed CRISPR screens. This includes evaluating the quality of guide RNA libraries, performing image analysis, evaluating assay results quality, data processing, hit identification, ranking, visualization, and biological interpretation.

Developing an Arrayed CRISPR-Cas9 Co-Culture Screen for Immuno-Oncology Target ID.

Immunotherapies including PD-L1 blockade have shown remarkable increases in the T cell-directed antitumor response; however, efficacy is seen only in a minority of patients. Recently, pooled CRISPR-Cas9 knockout (CRISPRn) screens in tumor/immune co-culture systems have identified a number of genes that confer resistance to T cell killing in pathways including antigen presentation and cytokine signaling, providing insight into tumor mechanisms that cause resistance to immunotherapies. The development of an arrayed CRISPRn screen in a tumor/immune co-culture system would allow the identification of novel targets for immuno-oncology, characterization of hits from pooled screens, and multiple assay endpoints to be measured per gene. Here, a small-scale arrayed CRISPRn screen was successfully developed to investigate the effects on a co-culture of T cells and Cas9-expressing PC9 lung adenocarcinoma cells modified to express anti-CD3 antibody on the cell surface (PC9-OKT3 T cell system). A focused CRISPRn library was designed to target genes involved in known resistance mechanisms (including antigen presentation, cytokine signaling, and apoptosis) as well as genes involved in immune synapse interactions. The viability of PC9 cells was assessed in two-dimensional adherent co-cultures via longitudinal imaging analysis. Knockout of epidermal growth factor receptor (EGFR) and PLK1 in tumor cells cultured alone or with T cells resulted in increased tumor cell death, as expected, whereas knockout of the test gene ICAM1 showed subtle donor-specific resistance to T cell killing. Taken together, these data provide proof of concept for arrayed CRISPRn screens in tumor/immune co-culture systems and warrant further investigation of in vitro co-culture models.

Arrayed CRISPR Screening Identifies Novel Targets That Enhance the Productive Delivery of mRNA by MC3-Based Lipid Nanoparticles.

Modified messenger RNAs (mRNAs) hold great potential as therapeutics by using the body's own processes for protein production. However, a key challenge is efficient delivery of therapeutic mRNA to the cell cytosol and productive protein translation. Lipid nanoparticles (LNPs) are the most clinically advanced system for nucleic acid delivery; however, a relatively narrow therapeutic index makes them unsuitable for many therapeutic applications. A key obstacle to the development of more potent LNPs is a limited mechanistic understanding of the interaction of LNPs with cells. To address this gap, we performed an arrayed CRISPR screen to identify novel pathways important for the functional delivery of MC3 lipid-based LNP encapsulated mRNA (LNP-mRNA). Here, we have developed and validated a robust, high-throughput screening-friendly phenotypic assay to identify novel targets that modulate productive LNP-mRNA delivery. We screened the druggable genome (7795 genes) and validated 44 genes that either increased (37 genes) or inhibited (14 genes) the productive delivery of LNP-mRNA. Many of these genes clustered into families involved with host cell transcription, protein ubiquitination, and intracellular trafficking. We show that both UDP-glucose ceramide glucosyltransferase and V-type proton ATPase can significantly modulate the productive delivery of LNP-mRNA, increasing and decreasing, respectively, with both genetic perturbation and by small-molecule inhibition. Taken together, these findings shed new light into the molecular machinery regulating the delivery of LNPs into cells and improve our mechanistic understanding of the cellular processes modulating the interaction of LNPs with cells.

A new prediction strategy for long local protein structures using an original description.

A relevant and accurate description of three-dimensional (3D) protein structures can be achieved by characterizing recurrent local structures. In a previous study, we developed a library of 120 3D structural prototypes encompassing all known 11-residues long local protein structures and ensuring a good quality of structural approximation. A local structure prediction method was also proposed. Here, overlapping properties of local protein structures in global ones are taken into account to characterize frequent local networks. At the same time, we propose a new long local structure prediction strategy which involves the use of evolutionary information coupled with Support Vector Machines (SVMs). Our prediction is evaluated by a stringent geometrical assessment. Every local structure prediction with a Calpha RMSD less than 2.5 A from the true local structure is considered as correct. A global prediction rate of 63.1% is then reached, corresponding to an improvement of 7.7 points compared with the previous strategy. In the same way, the prediction of 88.33% of the 120 structural classes is improved with 8.65% mean gain. 85.33% of proteins have better prediction results with a 9.43% average gain. An analysis of prediction rate per local network also supports the global improvement and gives insights into the potential of our method for predicting super local structures. Moreover, a confidence index for the direct estimation of prediction quality is proposed. Finally, our method is proved to be very competitive with cutting-edge strategies encompassing three categories of local structure predictions.

In silico prediction of protein flexibility with local structure approach.

Flexibility is an intrinsic essential feature of protein structures, directly linked to their functions. To this day, most of the prediction methods use the crystallographic data (namely B-factors) as the only indicator of protein's inner flexibility and predicts them as rigid or flexible. PredyFlexy stands differently from other approaches as it relies on the definition of protein flexibility (i) not only taken from crystallographic data, but also (ii) from Root Mean Square Fluctuation (RMSFs) observed in Molecular Dynamics simulations. It also uses a specific representation of protein structures, named Long Structural Prototypes (LSPs). From Position-Specific Scoring Matrix, the 120 LSPs are predicted with a good accuracy and directly used to predict (i) the protein flexibility in three categories (flexible, intermediate and rigid), (ii) the normalized B-factors, (iii) the normalized RMSFs, and (iv) a confidence index. Prediction accuracy among these three classes is equivalent to the best two class prediction methods, while the normalized B-factors and normalized RMSFs have a good correlation with experimental and in silico values. Thus, PredyFlexy is a unique approach, which is of major utility for the scientific community. It support parallelization features and can be run on a local cluster using multiple cores.

Drug Screening in Human PSC-Cardiac Organoids Identifies Pro-proliferative Compounds Acting via the Mevalonate Pathway.

We have previously developed a high-throughput bioengineered human cardiac organoid (hCO) platform, which provides functional contractile tissue with biological properties similar to native heart tissue, including mature, cell-cycle-arrested cardiomyocytes. In this study, we perform functional screening of 105 small molecules with pro-regenerative potential. Our findings reveal surprising discordance between our hCO system and traditional 2D assays. In addition, functional analyses uncovered detrimental effects of many hit compounds. Two pro-proliferative small molecules without detrimental impacts on cardiac function were identified. High-throughput proteomics in hCO revealed synergistic activation of the mevalonate pathway and a cell-cycle network by the pro-proliferative compounds. Cell-cycle reentry in hCO and in vivo required the mevalonate pathway as inhibition of the mevalonate pathway with a statin attenuated pro-proliferative effects. This study highlights the utility of human cardiac organoids for pro-regenerative drug development, including identification of underlying biological mechanisms and minimization of adverse side effects.

Protein contacts, inter-residue interactions and side-chain modelling.

Three-dimensional structures of proteins are the support of their biological functions. Their folds are stabilized by contacts between residues. Inner protein contacts are generally described through direct atomic contacts, i.e. interactions between side-chain atoms, while contact prediction methods mainly used inter-Calpha distances. In this paper, we have analyzed the protein contacts on a recent high quality non-redundant databank using different criteria. First, we have studied the average number of contacts depending on the distance threshold to define a contact. Preferential contacts between types of amino acids have been highlighted. Detailed analyses have been done concerning the proximity of contacts in the sequence, the size of the proteins and fold classes. The strongest differences have been extracted, highlighting important residues. Then, we studied the influence of five different side-chain conformation prediction methods (SCWRL, IRECS, SCAP, SCATD and SCCOMP) on the distribution of contacts. The prediction rates of these different methods are quite similar. However, using a distance criterion between side chains, the results are quite different, e.g. SCAP predicts 50% more contacts than observed, unlike other methods that predict fewer contacts than observed. Contacts deduced are quite distinct from one method to another with at most 75% contacts in common. Moreover, distributions of amino acid preferential contacts present unexpected behaviours distinct from previously observed in the X-ray structures, especially at the surface of proteins. For instance, the interactions involving Tryptophan greatly decrease.

The Role of Historical Bioactivity Data in the Deconvolution of Phenotypic Screens.

A substantial challenge in phenotypic drug discovery is the identification of the molecular targets that govern a phenotypic response of interest. Several experimental strategies are available for this, the so-called target deconvolution process. Most of these approaches exploit the affinity between a small-molecule compound and its putative targets or use large-scale genetic manipulations and profiling. Each of these methods has strengths but also limitations such as bias toward high-affinity interactions or risks from genetic compensation. The use of computational methods for target and mechanism of action identification is a complementary approach that can influence each step of a phenotypic screening campaign. Here, we describe how cheminformatics and bioinformatics are embedded in the process from initial selection of a focused compound library from a large set of historical small-molecule screens through the analysis of screening results. We present a deconvolution method based on enrichment analysis and using known bioactivity data of screened compounds to infer putative targets, pathways, and biological processes that are consistent with the observed phenotypic response. As an example, the approach is applied to a cellular screen aiming at identifying inhibitors of tumor necrosis factor-α production in lipopolysaccharide-stimulated THP-1 cells. In summary, we find that the approach can contribute to solving the often very complex target deconvolution task.

Systematic exploration of dual-acting modulators from a combined medicinal chemistry and biology perspective.

The traditional drug discovery strategy of pursuing "one compound-one target" has had difficulties delivering novel therapies for complex diseases currently lacking adequate treatments. An alternative and complementary approach is the design of multitargeted modulators simultaneously addressing multiple pathological mechanisms or overcoming pathway robustness. In this study, we propose a methodology to increase the probability of success for developing dual-acting modulators by systematically and rationally evaluating all dual-acting modulator opportunities within a specific disease area. This approach employs a combination of a five-step medicinal chemistry evaluation and a two-step biological analysis to help select the optimal target combination. It provides a novel methodology suitable for widespread application across disease areas. To exemplify the power of this approach, we focus on an analysis of the gastrointestinal (GI) disease area to identify opportunities supported by current literature data.

A short survey on protein blocks.

Protein structures are classically described in terms of secondary structures. Even if the regular secondary structures have relevant physical meaning, their recognition from atomic coordinates has some important limitations such as uncertainties in the assignment of boundaries of helical and ß-strand regions. Further, on an average about 50% of all residues are assigned to an irregular state, i.e., the coil. Thus different research teams have focused on abstracting conformation of protein backbone in the localized short stretches. Using different geometric measures, local stretches in protein structures are clustered in a chosen number of states. A prototype representative of the local structures in each cluster is generally defined. These libraries of local structures prototypes are named as "structural alphabets". We have developed a structural alphabet, named Protein Blocks, not only to approximate the protein structure, but also to predict them from sequence. Since its development, we and other teams have explored numerous new research fields using this structural alphabet. We review here some of the most interesting applications.

Analysis of protein contacts into Protein Units.

Three-dimensional structures of proteins are the support of their biological functions. Their folds are maintained by inter-residue interactions which are one of the main focuses to understand the mechanisms of protein folding and stability. Furthermore, protein structures can be composed of single or multiple functional domains that can fold and function independently. Hence, dividing a protein into domains is useful for obtaining an accurate structure and function determination. In previous studies, we enlightened protein contact properties according to different definitions and developed a novel methodology named Protein Peeling. Within protein structures, Protein Peeling characterizes small successive compact units along the sequence called protein units (PUs). The cutting done by Protein Peeling maximizes the number of contacts within the PUs and minimizes the number of contacts between them. This method is so a relevant tool in the context of the protein folding research and particularly regarding the hierarchical model proposed by George Rose. Here, we accurately analyze the PUs at different levels of cutting, using a non-redundant protein databank. Distribution of PU sizes, number of PUs or their accessibility are screened to determine their common and different features. Moreover, we highlight the preferential amino acid interactions inside and between PUs. Our results show that PUs are clearly an intermediate level between secondary structures and protein structural domains.

Analysis of loop boundaries using different local structure assignment methods.

Loops connect regular secondary structures. In many instances, they are known to play important biological roles. Analysis and prediction of loop conformations depend directly on the definition of repetitive structures. Nonetheless, the secondary structure assignment methods (SSAMs) often lead to divergent assignments. In this study, we analyzed, both structure and sequence point of views, how the divergence between different SSAMs affect boundary definitions of loops connecting regular secondary structures. The analysis of SSAMs underlines that no clear consensus between the different SSAMs can be easily found. Because these latter greatly influence the loop boundary definitions, important variations are indeed observed, that is, capping positions are shifted between different SSAMs. On the other hand, our results show that the sequence information in these capping regions are more stable than expected, and, classical and equivalent sequence patterns were found for most of the SSAMs. This is, to our knowledge, the most exhaustive survey in this field as (i) various databank have been used leading to similar results without implication of protein redundancy and (ii) the first time various SSAMs have been used. This work hence gives new insights into the difficult question of assignment of repetitive structures and addresses the issue of loop boundaries definition. Although SSAMs give very different local structure assignments capping sequence patterns remain efficiently stable.

Protein short loop prediction in terms of a structural alphabet.

Loops connect regular secondary structures. In many instances, they are known to play crucial biological roles. To bypass the limitation of secondary structure description, we previously defined a structural alphabet composed of 16 structural prototypes, called Protein Blocks (PBs). It leads to an accurate description of every region of 3D protein backbones and has been used in local structure prediction. In the present study, we used our structural alphabet to predict the loops connecting two repetitive structures. Thus, we showed interest to take into account the flanking regions, leading to prediction rate improvement up to 19.8%, but we also underline the sensitivity of such an approach. This research can be used to propose different structures for the loops and to probe and sample their flexibility. It is a useful tool for ab initio loop prediction and leads to insights into flexible docking approach.

Analysis of protein chameleon sequence characteristics.

Conversion of local structural state of a protein from an alpha-helix to a beta-strand is usually associated with a major change in the tertiary structure. Similar changes were observed during the self assembly of amyloidogenic proteins to form fibrils, which are implicated in severe diseases conditions, e.g., Alzheimer disease. Studies have emphasized that certain protein sequence fragments known as chameleon sequences do not have a strong preference for either helical or the extended conformations. Surprisingly, the information on the local sequence neighborhood can be used to predict their secondary at a high accuracy level. Here we report a large scale-analysis of chameleon sequences to estimate their propensities to be associated with different local structural states such as alpha -helices, beta-strands and coils. With the help of the propensity information derived from the amino acid composition, we underline their complexity, as more than one quarter of them prefers coil state over to the regular secondary structures. About half of them show preference for both alpha-helix and beta-sheet conformations and either of these two states is favored by the rest.