Dynamics of Protein-RNA Interfaces Using All-Atom Molecular Dynamics Simulations.

Facing the current challenges posed by human health diseases requires the understanding of cell machinery at a molecular level. The interplay between proteins and RNA is key for any physiological phenomenon, as well protein-RNA interactions. To understand these interactions, many experimental techniques have been developed, spanning a very wide range of spatial and temporal resolutions. In particular, the knowledge of tridimensional structures of protein-RNA complexes provides structural, mechanical, and dynamical pieces of information essential to understand their functions. To get insights into the dynamics of protein-RNA complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on nine different protein-RNA complexes with different functions and interface size by taking into account the bound and unbound forms. First, we characterized structural changes upon binding and, for the RNA part, the change in the puckering. Second, we extensively analyzed the interfaces, their dynamics and structural properties, and the structural waters involved in the binding, as well as the contacts mediated by them. Based on our analysis, the interfaces rearranged during the simulation time showing alternative and stable residue-residue contacts with respect to the experimental structure.

AlphaFold2 Predicts Whether Proteins Interact Amidst Confounding Structural Compatibility.

Predicting whether two proteins physically interact is one of the holy grails of computational biology, galvanized by rapid advancements in deep learning. AlphaFold2, although not developed with this goal, is promising in this respect. Here, I test the prediction capability of AlphaFold2 on a very challenging data set, where proteins are structurally compatible, even when they do not interact. AlphaFold2 achieves high discrimination between interacting and non-interacting proteins, and the cases of misclassifications can either be rescued by revisiting the input sequences or can suggest false positives and negatives in the data set. AlphaFold2 is thus not impaired by the compatibility between protein structures and has the potential to be applied on a large scale.

Values for the Digestibility of Pea Protein Isolate or Casein Amino Acids Determined using the Dual Isotope Method Are Not Similar to Those Derived with the Standard Ileal Balance Method in Healthy Volunteers.

BACKGROUND: The measurement of ileal amino acid (AA) digestibility is invasive and inappropriate when applied to vulnerable populations. The dual isotope method has been developed over the past 5 y as an alternative method. OBJECTIVE: The aim of this work was to compare the indispensable amino acid (IAA) digestibility values of 2 different proteins obtained using the dual isotope and the standard ileal balance methods in the same subjects. METHODS: Fifteen healthy adults completed the study. Over 4 h, they ingested 9 successive portions of mashed potatoes containing the test protein (pea protein or casein) labeled intrinsically with 15N and 2H, and a 13C-free AA mixture as a reference for the dual isotope method. Plasma was sampled regularly over the 8-h postprandial period, whereas the ileal digesta was collected continuously via a naso-ileal tube. Isotopic enrichments (15N and 13C) were measured in the digesta for the direct determination of ileal IAA digestibility, whereas plasma enrichments (2H and 13C) were measured to determine IAA digestibility using the dual isotope method. RESULTS: The 4-h repeated meal procedure enabled the almost complete digestion of test proteins at 8 h and the attainment of a plasma isotopic plateau between 2.5 and 4 h. These conditions were necessary to perform the ileal balance and dual isotope methods simultaneously. For pea protein, the mean IAA digestibility was similar between the 2 methods, but significant differences (from 10% to 20%) were observed for individual IAA values. For casein, IAA digestibility was significantly lower with the dual isotope method for all the IAA analyzed. CONCLUSIONS: Under our experimental conditions, the degree of agreement between the dual isotope and ileal balance methods varied among AAs and depended on the protein source. Further research is needed to validate the dual isotope method. This study was registered at clinicaltrials.gov as NCT04072770.

Rational Prediction of PROTAC-Compatible Protein-Protein Interfaces by Molecular Docking.

Proteolysis targeting chimeras (PROTACs) are heterobifunctional ligands that mediate the interaction between a protein target and an E3 ligase, resulting in a ternary complex, whose interaction with the ubiquitination machinery leads to target degradation. This technology is emerging as an exciting new avenue for therapeutic development, with several PROTACs currently undergoing clinical trials targeting cancer. Here, we describe a general and computationally efficient methodology combining restraint-based docking, energy-based rescoring, and a filter based on the minimal solvent-accessible surface distance to produce PROTAC-compatible PPIs suitable for when there is no a priori known PROTAC ligand. In a benchmark employing a manually curated data set of 13 ternary complex crystals, we achieved an accuracy of 92% when starting from bound structures and 77% when starting from unbound structures, respectively. Our method only requires that the ligand-bound structures of the monomeric forms of the E3 ligase and target proteins be given to run, making it general, accurate, and highly efficient, with the ability to impact early-stage PROTAC-based drug design campaigns where no structural information about the ternary complex structure is available.

Discriminating physiological from non-physiological interfaces in structures of protein complexes: A community-wide study.

Reliably scoring and ranking candidate models of protein complexes and assigning their oligomeric state from the structure of the crystal lattice represent outstanding challenges. A community-wide effort was launched to tackle these challenges. The latest resources on protein complexes and interfaces were exploited to derive a benchmark dataset consisting of 1677 homodimer protein crystal structures, including a balanced mix of physiological and non-physiological complexes. The non-physiological complexes in the benchmark were selected to bury a similar or larger interface area than their physiological counterparts, making it more difficult for scoring functions to differentiate between them. Next, 252 functions for scoring protein-protein interfaces previously developed by 13 groups were collected and evaluated for their ability to discriminate between physiological and non-physiological complexes. A simple consensus score generated using the best performing score of each of the 13 groups, and a cross-validated Random Forest (RF) classifier were created. Both approaches showed excellent performance, with an area under the Receiver Operating Characteristic (ROC) curve of 0.93 and 0.94, respectively, outperforming individual scores developed by different groups. Additionally, AlphaFold2 engines recalled the physiological dimers with significantly higher accuracy than the non-physiological set, lending support to the reliability of our benchmark dataset annotations. Optimizing the combined power of interface scoring functions and evaluating it on challenging benchmark datasets appears to be a promising strategy.

Transformative adaptation through nature-based solutions: a comparative case study analysis in China, Italy, and Germany.

This paper explores how claims for transformative adaptation toward more equitable and sustainable societies can be assessed. We build on a theoretical framework describing transformative adaptation as it manifests across four core elements of the public-sector adaptation lifecycle: vision, planning, institutional frameworks, and interventions. For each element, we identify characteristics that can help track adaptation as transformative. Our purpose is to identify how governance systems can constrain or support transformative choices and thus enable targeted interventions. We demonstrate and test the usefulness of the framework with reference to three government-led adaptation projects of nature-based solutions (NBS): river restoration (Germany), forest conservation (China), and landslide risk reduction (Italy). Building on a desktop study and open-ended interviews, our analysis adds evidence to the view that transformation is not an abrupt system change, but a dynamic complex process that evolves over time. While each of the NBS cases fails to fulfill all the transformation characteristics, there are important transformative elements in their visions, planning, and interventions. There is a deficit, however, in the transformation of institutional frameworks. The cases show institutional commonalities in multi-scale and cross-sectoral (polycentric) collaboration as well as innovative processes for inclusive stakeholder engagement; yet, these arrangements are ad hoc, short-term, dependent on local champions, and lacking the permanency needed for upscaling. For the public sector, this result highlights the potential for establishing cross-competing priorities among agencies, cross-sectoral formal mechanisms, new dedicated institutions, and programmatic and regulatory mainstreaming. Supplementary Information: The online version contains supplementary material available at 10.1007/s10113-023-02066-7.

Internal Normal Mode Analysis Applied to RNA Flexibility and Conformational Changes.

We investigated the capability of internal normal modes to reproduce RNA flexibility and predict observed RNA conformational changes and, notably, those induced by the formation of RNA-protein and RNA-ligand complexes. Here, we extended our iNMA approach developed for proteins to study RNA molecules using a simplified representation of the RNA structure and its potential energy. Three data sets were also created to investigate different aspects. Despite all the approximations, our study shows that iNMA is a suitable method to take into account RNA flexibility and describe its conformational changes opening the route to its applicability in any integrative approach where these properties are crucial.

Conformational space exploration of cryo-EM structures by variability refinement.

Cryo-EM observation of biological samples enables visualization of sample heterogeneity, in the form of discrete states that are separable, or continuous heterogeneity as a result of local protein motion before flash freezing. Variability analysis of this continuous heterogeneity describes the variance between a particle stack and a volume, and results in a map series describing the various steps undertaken by the sample in the particle stack. While this observation is absolutely stunning, it is very hard to pinpoint structural details to elements of the maps. In order to bridge the gap between observation and explanation, we designed a tool that refines an ensemble of structures into all the maps from variability analysis. Using this bundle of structures, it is easy to spot variable parts of the structure, as well as the parts that are not moving. Comparison with molecular dynamics simulations highlights the fact that the movements follow the same directions, albeit with different amplitudes. Ligand can also be investigated using this method. Variability refinement is available in the Phenix software suite, accessible under the program name phenix.varref.

Specific Xray diffraction patterns of membrane proteins caused by secondary structure collinearity.

Diffraction anisotropy is a phenomenon that impacts more specifically membrane proteins, compared to soluble ones, but the reasons for this discrepancy remained unclear. Often, it is referred to a difference in resolution limits between highest and lowest diffraction limits as a signature for anisotropy. We show in this article that there is no single correlation between anisotropy and difference in resolution limits, with notably a substantial number of structures displaying various anisotropy with no difference in resolution limits. We further investigated diffraction intensity profiles, and observed a peak centred on 4.9 Å resolution more predominant in membrane proteins. Since this peak is in the region corresponding to secondary structures, we investigated the influence of secondary structure ratio. We showed that secondary structure content has little influence on this profile, while secondary structure collinearity in membrane proteins correlate with a stronger peak. Finally, we could further show that the presence of this peak is linked to higher diffraction anisotropy. These results bring to light a specific diffraction of membrane protein crystals, which calls for a specific handling by crystallographic software. It also brings an explanation for investigators struggling with their anisotropic data.

A dynamical view of protein-protein complexes: Studies by molecular dynamics simulations.

Protein-protein interactions are at the basis of many protein functions, and the knowledge of 3D structures of protein-protein complexes provides structural, mechanical and dynamical pieces of information essential to understand these functions. Protein-protein interfaces can be seen as stable, organized regions where residues from different partners form non-covalent interactions that are responsible for interaction specificity and strength. They are commonly described as a peripheral region, whose role is to protect the core region that concentrates the most contributing interactions, from the solvent. To get insights into the dynamics of protein-protein complexes, we carried out all-atom molecular dynamics simulations in explicit solvent on eight different protein-protein complexes of different functional class and interface size by taking into account the bound and unbound forms. On the one hand, we characterized structural changes upon binding of the proteins, and on the other hand we extensively analyzed the interfaces and the structural waters involved in the binding. Based on our analysis, in 6 cases out of 8, the interfaces rearranged during the simulation time, in stable and long-lived substates with alternative residue-residue contacts. These rearrangements are not restricted to side-chain fluctuations in the periphery but also affect the core interface. Finally, the analysis of the waters at the interface and involved in the binding pointed out the importance to take into account their role in the estimation of the interaction strength.

Pharmacomodulation of a ligand targeting the HBV capsid hydrophobic pocket.

Hepatitis B virus (HBV) is a small enveloped retrotranscribing DNA virus and an important human pathogen. Its capsid-forming core protein (Cp) features a hydrophobic pocket proposed to be central notably in capsid envelopment. Indeed, mutations in and around this pocket can profoundly modulate, and even abolish, secretion of enveloped virions. We have recently shown that Triton X-100, a detergent used during Cp purification, binds to the hydrophobic pocket with micromolar affinity. We here performed pharmacomodulation of pocket binders through systematic modifications of the three distinct chemical moieties composing the Triton X-100 molecule. Using NMR and ITC, we found that the flat aromatic moiety is essential for binding, while the number of atoms of the aliphatic chain modulates binding affinity. The hydrophilic tail, in contrast, is highly tolerant to changes in both length and type. Our data provide essential information for designing a new class of HBV antivirals targeting capsid-envelope interactions.

Tertiary and Quaternary Structure Organization in GMP Synthetases: Implications for Catalysis.

Glutamine amidotransferases, enzymes that transfer nitrogen from Gln to various cellular metabolites, are modular, with the amidotransferase (GATase) domain hydrolyzing Gln, generating ammonia and the acceptor domain catalyzing the addition of nitrogen onto its cognate substrate. GMP synthetase (GMPS), an enzyme in the de novo purine nucleotide biosynthetic pathway, is a glutamine amidotransferase that catalyzes the synthesis of GMP from XMP. The reaction involves activation of XMP though adenylation by ATP in the ATP pyrophosphatase (ATPPase) active site, followed by channeling and attack of NH3 generated in the GATase pocket. This complex chemistry entails co-ordination of activity across the active sites, allosteric activation of the GATase domain to modulate Gln hydrolysis and channeling of ammonia from the GATase to the acceptor active site. Functional GMPS dimers associate through the dimerization domain. The crystal structure of the Gln-bound complex of Plasmodium falciparum GMPS (PfGMPS) for the first time revealed large-scale domain rotation to be associated with catalysis and leading to the juxtaposition of two otherwise spatially distal cysteinyl (C113/C337) residues. In this manuscript, we report on an unusual structural variation in the crystal structure of the C89A/C113A PfGMPS double mutant, wherein a larger degree of domain rotation has led to the dissociation of the dimeric structure. Furthermore, we report a hitherto overlooked signature motif tightly related to catalysis.

When Alphafold2 predictions go wrong for protein-protein complexes, is there something to be learnt?

In this short communication, I analyze cases of failed predictions for protein-protein complexes with Alphafold2, and show that they either point to erroneous annotation in the PDB or correct binding site regions.

Identification of shared tumor epitopes from endogenous retroviruses inducing high-avidity cytotoxic T cells for cancer immunotherapy.

Human endogenous retroviruses (HERVs) represent 8% of the human genome. HERV products may represent tumor antigens relevant for cancer immunotherapy. We developed a bioinformatic approach to identify shared CD8+ T cell epitopes derived from cancer-associated HERVs in solid tumors. Six candidates among the most commonly shared HLA-A2 epitopes with evidence of translation were selected for immunological evaluation. In vitro priming assays confirmed the immunogenicity of these epitopes, which induced high-avidity CD8+ T cell clones. These T cells specifically recognize and kill HLA-A2+ tumor cells presenting HERV epitopes on HLA molecules, as demonstrated by mass spectrometry. Furthermore, epitope-specific CD8+ T cells were identified by dextramer staining among tumor-infiltrating lymphocytes from HLA-A2+ patients with breast cancer. Last, we showed that HERV-specific T cells lyse patient-derived organoids. These shared virus-like epitopes are of major interest for the development of cancer vaccines or T cell-based immunotherapies, especially in tumors with low/intermediate mutational burden.

Substrate-bound and substrate-free outward-facing structures of a multidrug ABC exporter.

Multidrug ABC transporters translocate drugs across membranes by a mechanism for which the molecular features of drug release are so far unknown. Here, we resolved three ATP-Mg2+-bound outward-facing conformations of the Bacillus subtilis (homodimeric) BmrA by x-ray crystallography and single-particle cryo-electron microscopy (EM) in detergent solution, one of them with rhodamine 6G (R6G), a substrate exported by BmrA when overexpressed in B. subtilis. Two R6G molecules bind to the drug-binding cavity at the level of the outer leaflet, between transmembrane (TM) helices 1-2 of one monomer and TM5'-6' of the other. They induce a rearrangement of TM1-2, highlighting a local flexibility that we confirmed by hydrogen/deuterium exchange and molecular dynamics simulations. In the absence of R6G, simulations show a fast postrelease occlusion of the cavity driven by hydrophobicity, while when present, R6G can move within the cavity, maintaining it open.

Towards design of drugs and delivery systems with the Martini coarse-grained model.

Coarse-grained (CG) modelling with the Martini force field has come of age. By combining a variety of bead types and sizes with a new mapping approach, the newest version of the model is able to accurately simulate large biomolecular complexes at millisecond timescales. In this perspective, we discuss possible applications of the Martini 3 model in drug discovery and development pipelines and highlight areas for future development. Owing to its high simulation efficiency and extended chemical space, Martini 3 has great potential in the area of drug design and delivery. However, several aspects of the model should be improved before Martini 3 CG simulations can be routinely employed in academic and industrial settings. These include the development of automatic parameterisation protocols for a variety of molecule types, the improvement of backmapping procedures, the description of protein flexibility and the development of methodologies enabling efficient sampling. We illustrate our view with examples on key areas where Martini could give important contributions such as drugs targeting membrane proteins, cryptic pockets and protein-protein interactions and the development of soft drug delivery systems.

Real ileal amino acid digestibility of pea protein compared to casein in healthy humans: a randomized trial.

BACKGROUND: It is necessary to propose plant alternatives to animal proteins that are of good nutritional quality. Pea is a good candidate owing to its high protein content and its well-balanced amino acid (AA) profile. OBJECTIVES: This study aimed to assess the real ileal AA and nitrogen digestibility (RIDAA and RIDN) of pea protein isolate as compared to milk casein in humans. It also aimed to evaluate their nutritional quality through calculation of the digestible indispensable amino acid score (DIAAS) and to determine the net postprandial protein utilization (NPPU). METHODS: Fifteen healthy volunteers were included in a randomized, single-blinded, 2-arm, parallel-design trial. They were equipped with a naso-ileal tube. They ingested the test meals, which consisted of 9 successive portions of mashed potatoes containing either pea protein or casein, intrinsically labeled with nitrogen 15. Ileal content, plasma, and urine samples were collected regularly over an 8-h postprandial period. RESULTS: The mean RIDAA values were 93.6% ± 2.9% for pea protein and 96.8% ± 1.0% for casein, with no difference between the sources (P = 0.22). Leucine, valine, lysine, and phenylalanine were significantly less digestible in pea than in casein. The RIDN values were 92.0% ± 2.7% and 94.0% ± 1.7% for pea protein and casein, respectively, and were not different (P = 0.11). The DIAAS was 1.00 for pea protein and 1.45 for casein. The NPPU was 71.6% ± 6.2% and 71.2% ± 4.9% for pea protein and casein, respectively (P = 0.88). CONCLUSIONS: Although some AAs are less digestible in pea protein than in casein, the real ileal digestibility and the NPPU were not different. The DIAAS of 1.00 obtained for pea protein demonstrated its ability to meet all AA requirements. This study shows the potential of pea isolate as a high-quality protein. This study was registered at clinicaltrials.gov as NCT04072770.

The Candida glabrata glycogen branching enzyme structure reveals unique features of branching enzymes of the Saccharomycetaceae phylum.

Branching enzymes (BE) are responsible for the formation of branching points at the 1,6 position in glycogen and starch, by catalyzing the cleavage of α-1,4-linkages and the subsequent transfer by introducing α-1,6-linked glucose branched points. BEs are found in the large GH13 family, eukaryotic BEs being mainly classified in the GH13_8 subfamily, GH13_9 grouping almost exclusively prokaryotic enzymes. With the aim of contributing to the understanding of the mode of recognition and action of the enzymes belonging to GH13_8, and to the understanding of features distinguishing these enzymes from those belonging to subfamily 13_9, we solved the crystal structure of the glycogen branching enzyme (GBE) from the yeast Candida glabrata, CgGBE, in ligand-free forms and in complex with a maltotriose. The structures revealed the presence of a domain already observed in Homo sapiens and Oryza sativa BEs that we named α-helical N-terminal domain, in addition to the three conserved domains found in BE. We confirmed by phylogenetic analysis that this α-helical N-terminal domain is always present in the GH13_8 enzymes suggesting that it could actually present a signature for this subfamily. We identified two binding sites in the α-helical N-terminal domain and in the carbohydrate binding module 48 (CBM48), respectively, which show a unique structural organization only present in the Saccharomycotina phylum. Our structural and phylogenetic investigation provides new insight into the structural characterization of GH13_8 GBE revealing that unique structural features only present in the Saccharomycotina phylum thereby conferring original properties to this group of enzymes.

The Det.Belt Server: A Tool to Visualize and Estimate Amphipathic Solvent Belts around Membrane Proteins.

Detergents wrap around membrane proteins to form a belt covering the hydrophobic part of the protein serving for membrane insertion and interaction with lipids. The number of detergent monomers forming this belt is usually unknown to investigators, unless dedicated detergent quantification is undertaken, which for many projects is difficult to setup. Yet, having an approximate knowledge of the amount of detergent forming the belt is extremely useful, to better grasp the protein of interest in interaction with its direct environment rather than picturing the membrane protein "naked". We created the Det.Belt server to dress up membrane proteins and represent in 3D the bulk made by detergent molecules wrapping in a belt. Many detergents are included in a database, allowing investigators to screen in silico the effect of different detergents around their membrane protein. The input number of detergents is changeable with fast recomputation of the belt for interactive usage. Metrics representing the belt are readily available together with scripts to render quality 3D images for publication. The Det.Belt server is a tool for biochemists to better grasp their sample.