Annotating Macromolecular Complexes in the Protein Data Bank: Improving the FAIRness of Structure Data.

Macromolecular complexes are essential functional units in nearly all cellular processes, and their atomic-level understanding is critical for elucidating and modulating molecular mechanisms. The Protein Data Bank (PDB) serves as the global repository for experimentally determined structures of macromolecules. Structural data in the PDB offer valuable insights into the dynamics, conformation, and functional states of biological assemblies. However, the current annotation practices lack standardised naming conventions for assemblies in the PDB, complicating the identification of instances representing the same assembly. In this study, we introduce a method leveraging resources external to PDB, such as the Complex Portal, UniProt and Gene Ontology, to describe assemblies and contextualise them within their biological settings accurately. Employing the proposed approach, we assigned standard names to over 90% of unique assemblies in the PDB and provided persistent identifiers for each assembly. This standardisation of assembly data enhances the PDB, facilitating a deeper understanding of macromolecular complexes. Furthermore, the data standardisation improves the PDB's FAIR attributes, fostering more effective basic and translational research and scientific education.

Stages of OCP-FRP Interactions in the Regulation of Photoprotection in Cyanobacteria, Part 2: Small-Angle Neutron Scattering with Partial Deuteration.

We used small-angle neutron scattering partially coupled with size-exclusion chromatography to unravel the solution structures of two variants of the Orange Carotenoid Protein (OCP) lacking the N-terminal extension (OCP-ΔNTE) and its complex formation with the Fluorescence Recovery Protein (FRP). The dark-adapted, orange form OCP-ΔNTEO is fully photoswitchable and preferentially binds the pigment echinenone. Its complex with FRP consists of a monomeric OCP component, which closely resembles the compact structure expected for the OCP ground state, OCPO. In contrast, the pink form OCP-ΔNTEP, preferentially binding the pigment canthaxanthin, is mostly nonswitchable. The pink OCP form appears to occur as a dimer and is characterized by a separation of the N- and C-terminal domains, with the canthaxanthin embedded only into the N-terminal domain. Therefore, OCP-ΔNTEP can be viewed as a prototypical model system for the active, spectrally red-shifted state of OCP, OCPR. The dimeric structure of OCP-ΔNTEP is retained in its complex with FRP. Small-angle neutron scattering using partially deuterated OCP-FRP complexes reveals that FRP undergoes significant structural changes upon complex formation with OCP. The observed structures are assigned to individual intermediates of the OCP photocycle in the presence of FRP.

Docking-based long timescale simulation of cell-size protein systems at atomic resolution.

Computational methodologies are increasingly addressing modeling of the whole cell at the molecular level. Proteins and their interactions are the key component of cellular processes. Techniques for modeling protein interactions, thus far, have included protein docking and molecular simulation. The latter approaches account for the dynamics of the interactions but are relatively slow, if carried out at all-atom resolution, or are significantly coarse grained. Protein docking algorithms are far more efficient in sampling spatial coordinates. However, they do not account for the kinetics of the association (i.e., they do not involve the time coordinate). Our proof-of-concept study bridges the two modeling approaches, developing an approach that can reach unprecedented simulation timescales at all-atom resolution. The global intermolecular energy landscape of a large system of proteins was mapped by the pairwise fast Fourier transform docking and sampled in space and time by Monte Carlo simulations. The simulation protocol was parametrized on existing data and validated on a number of observations from experiments and molecular dynamics simulations. The simulation protocol performed consistently across very different systems of proteins at different protein concentrations. It recapitulated data on the previously observed protein diffusion rates and aggregation. The speed of calculation allows reaching second-long trajectories of protein systems that approach the size of the cells, at atomic resolution.

3D time-lapse imaging of a mouse embryo using intensity diffraction tomography embedded inside a deep learning framework.

We present a compact 3D diffractive microscope that can be inserted directly in a cell incubator for long-term observation of developing organisms. Our setup is particularly simple and robust, since it does not include any moving parts and is compatible with commercial cell culture containers. It has been designed to image large specimens (>100×100×100µm3) with subcellular resolution. The sample's optical properties [refractive index (RI) and absorption] are reconstructed in 3D from intensity-only images recorded with different illumination angles produced by an LED array. The reconstruction is performed using the beam propagation method embedded inside a deep-learning network where the layers encode the optical properties of the object. This deep neural network is trained for a given multiangle intensity acquisition. After training, the weights of the neural network deliver the 3D distribution of the optical properties of the sample. The effect of spherical aberrations due to the sample holder/air interfaces are taken into account in the forward model. Using this approach, we performed time-lapse 3D imaging of preimplantation mouse embryos over six days. Images of embryos from a single cell (low-scattering regime) to the blastocyst stage (highly scattering regime) were successfully reconstructed. Due to its subcellular resolution, our system can provide quantitative information on the embryos' development and viability. Hence, this technology opens what we believe to be novel opportunities for 3D label-free live-cell imaging of whole embryos or organoids over long observation times.

Session introduction: AI-driven Advances in Modeling of Protein Structure.

The last few years mark dramatic improvements in modeling of protein structure. Progress was initially due to breakthroughs in residue-residue contact prediction, first with global statistical models and later with deep learning. These advancements were then followed by an even broader application of the deep learning techniques to the protein structure modeling itself, first using Convolutional Neural Networks (CNNs) and then switching to Natural Language Processing (NLP), including Attention models, and to Geometric Deep Learning (GDL). The accuracy of protein structure models generated with current state-of-the-art methods rivals that of experimental structures, while models themselves are used to solve structures or to make them more accurate.Looking at the near future of machine learning applications in structural biology, we ask the following questions: Which specific problems should we expect to be solved next? Which new methods will prove to be the most effective? Which actions are likely to stimulate further progress the most? In addressing these questions, we invite the 2022 PSB attendees to actively participate in session discussions.The AI-driven Advances in Modeling of Protein Structure session includes five papers specifically dedicated to:Evaluating the significance of training data selection in machine learning.Geometric pattern transferability, from protein self-interactions to protein-ligand interactions.Supervised versus unsupervised sequence to contact learning, using attention models.Side chain packing using SE(3) transformers.Feature detection in electrostatic representations of ligand binding sites.

Protein sequence-to-structure learning: Is this the end(-to-end revolution)?

The potential of deep learning has been recognized in the protein structure prediction community for some time, and became indisputable after CASP13. In CASP14, deep learning has boosted the field to unanticipated levels reaching near-experimental accuracy. This success comes from advances transferred from other machine learning areas, as well as methods specifically designed to deal with protein sequences and structures, and their abstractions. Novel emerging approaches include (i) geometric learning, that is, learning on representations such as graphs, three-dimensional (3D) Voronoi tessellations, and point clouds; (ii) pretrained protein language models leveraging attention; (iii) equivariant architectures preserving the symmetry of 3D space; (iv) use of large meta-genome databases; (v) combinations of protein representations; and (vi) finally truly end-to-end architectures, that is, differentiable models starting from a sequence and returning a 3D structure. Here, we provide an overview and our opinion of the novel deep learning approaches developed in the last 2 years and widely used in CASP14.

HOPMA: Boosting Protein Functional Dynamics with Colored Contact Maps.

In light of the recent very rapid progress in protein structure prediction, accessing the multitude of functional protein states is becoming more central than ever before. Indeed, proteins are flexible macromolecules, and they often perform their function by switching between different conformations. However, high-resolution experimental techniques such as X-ray crystallography and cryogenic electron microscopy can catch relatively few protein functional states. Many others are only accessible under physiological conditions in solution. Therefore, there is a pressing need to fill this gap with computational approaches. We present HOPMA, a novel method to predict protein functional states and transitions by using a modified elastic network model. The method exploits patterns in a protein contact map, taking its 3D structure as input, and excludes some disconnected patches from the elastic network. Combined with nonlinear normal mode analysis, this strategy boosts the protein conformational space exploration, especially when the input structure is highly constrained, as we demonstrate on a set of more than 400 transitions. Our results let us envision the discovery of new functional conformations, which were unreachable previously, starting from the experimentally known protein structures. The method is computationally efficient and available at https://github.com/elolaine/HOPMA and https://team.inria.fr/nano-d/software/nolb-normal-modes.

VoroCNN: deep convolutional neural network built on 3D Voronoi tessellation of protein structures.

MOTIVATION: Effective use of evolutionary information has recently led to tremendous progress in computational prediction of three-dimensional (3D) structures of proteins and their complexes. Despite the progress, the accuracy of predicted structures tends to vary considerably from case to case. Since the utility of computational models depends on their accuracy, reliable estimates of deviation between predicted and native structures are of utmost importance. RESULTS: For the first time, we present a deep convolutional neural network (CNN) constructed on a Voronoi tessellation of 3D molecular structures. Despite the irregular data domain, our data representation allows us to efficiently introduce both convolution and pooling operations and train the network in an end-to-end fashion without precomputed descriptors. The resultant model, VoroCNN, predicts local qualities of 3D protein folds. The prediction results are competitive to state of the art and superior to the previous 3D CNN architectures built for the same task. We also discuss practical applications of VoroCNN, for example, in recognition of protein binding interfaces. AVAILABILITY AND IMPLEMENTATION: The model, data and evaluation tests are available at https://team.inria.fr/nano-d/software/vorocnn/. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

KORP-PL: a coarse-grained knowledge-based scoring function for protein-ligand interactions.

MOTIVATION: Despite the progress made in studying protein-ligand interactions and the widespread application of docking and affinity prediction tools, improving their precision and efficiency still remains a challenge. Computational approaches based on the scoring of docking conformations with statistical potentials constitute a popular alternative to more accurate but costly physics-based thermodynamic sampling methods. In this context, a minimalist and fast sidechain-free knowledge-based potential with a high docking and screening power can be very useful when screening a big number of putative docking conformations. RESULTS: Here, we present a novel coarse-grained potential defined by a 3D joint probability distribution function that only depends on the pairwise orientation and position between protein backbone and ligand atoms. Despite its extreme simplicity, our approach yields very competitive results with the state-of-the-art scoring functions, especially in docking and screening tasks. For example, we observed a twofold improvement in the median 5% enrichment factor on the DUD-E benchmark compared to Autodock Vina results. Moreover, our results prove that a coarse sidechain-free potential is sufficient for a very successful docking pose prediction. AVAILABILITYAND IMPLEMENTATION: The standalone version of KORP-PL with the corresponding tests and benchmarks are available at https://team.inria.fr/nano-d/korp-pl/ and https://chaconlab.org/modeling/korp-pl. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Hydroxylation of Antitubercular Drug Candidate, SQ109, by Mycobacterial Cytochrome P450.

Spreading of the multidrug-resistant (MDR) strains of the one of the most harmful pathogen Mycobacterium tuberculosis (Mtb) generates the need for new effective drugs. SQ109 showed activity against resistant Mtb and already advanced to Phase II/III clinical trials. Fast SQ109 degradation is attributed to the human liver Cytochrome P450s (CYPs). However, no information is available about interactions of the drug with Mtb CYPs. Here, we show that Mtb CYP124, previously assigned as a methyl-branched lipid monooxygenase, binds and hydroxylates SQ109 in vitro. A 1.25 Å-resolution crystal structure of the CYP124-SQ109 complex unambiguously shows two conformations of the drug, both positioned for hydroxylation of the ω-methyl group in the trans position. The hydroxylated SQ109 presumably forms stabilizing H-bonds with its target, Mycobacterial membrane protein Large 3 (MmpL3). We anticipate that Mtb CYPs could function as analogs of drug-metabolizing human CYPs affecting pharmacokinetics and pharmacodynamics of antitubercular (anti-TB) drugs.

AnAnaS: Software for Analytical Analysis of Symmetries in Protein Structures.

Symmetry is very common among proteins found in structural databases such as the Protein Data Bank (PDB). We present novel software, called AnAnaS, that finds positions and orientations of the symmetry axes in all types of symmetrical protein assemblies. It deals with five symmetry groups: cyclic, dihedral, tetrahedral, octahedral, and icosahedral. The software also assesses the quality of symmetry and can detect symmetries in incomplete cyclic assemblies. Internally, AnAnaS comprises discrete and continuous optimization steps and is applicable to assemblies with multiple chains in the asymmetric subunits or to those with pseudosymmetry. The method is very fast as most of the steps are performed analytically.

Interdomain Flexibility within NADPH Oxidase Suggested by SANS Using LMNG Stealth Carrier.

Small angle neutron scattering (SANS) provides a method to obtain important low-resolution information for integral membrane proteins (IMPs), challenging targets for structural determination. Specific deuteration furnishes a "stealth" carrier for the solubilized IMP. We used SANS to determine a structural envelope of SpNOX, the Streptococcus pneumoniae NADPH oxidase (NOX), a prokaryotic model system for exploring structure and function of eukaryotic NOXes. SpNOX was solubilized in the detergent lauryl maltose neopentyl glycol, which provides optimal SpNOX stability and activity. Using deuterated solvent and protein, the lauryl maltose neopentyl glycol was experimentally undetected in SANS. This affords a cost-effective SANS approach for obtaining novel structural information on IMPs. Combining SANS data with molecular modeling provided a first, to our knowledge, structural characterization of an entire NOX enzyme. It revealed a distinctly less compact structure than that predicted from the docking of homologous crystal structures of the separate transmembrane and dehydrogenase domains, consistent with a flexible linker connecting the two domains.

Combining molecular dynamics simulations with small-angle X-ray and neutron scattering data to study multi-domain proteins in solution.

Many proteins contain multiple folded domains separated by flexible linkers, and the ability to describe the structure and conformational heterogeneity of such flexible systems pushes the limits of structural biology. Using the three-domain protein TIA-1 as an example, we here combine coarse-grained molecular dynamics simulations with previously measured small-angle scattering data to study the conformation of TIA-1 in solution. We show that while the coarse-grained potential (Martini) in itself leads to too compact conformations, increasing the strength of protein-water interactions results in ensembles that are in very good agreement with experiments. We show how these ensembles can be refined further using a Bayesian/Maximum Entropy approach, and examine the robustness to errors in the energy function. In particular we find that as long as the initial simulation is relatively good, reweighting against experiments is very robust. We also study the relative information in X-ray and neutron scattering experiments and find that refining against the SAXS experiments leads to improvement in the SANS data. Our results suggest a general strategy for studying the conformation of multi-domain proteins in solution that combines coarse-grained simulations with small-angle X-ray scattering data that are generally most easy to obtain. These results may in turn be used to design further small-angle neutron scattering experiments that exploit contrast variation through 1H/2H isotope substitutions.

Predicting Protein Functional Motions: an Old Recipe with a New Twist.

Large macromolecules, including proteins and their complexes, very often adopt multiple conformations. Some of them can be seen experimentally, for example with x-ray crystallography or cryo-electron microscopy. This structural heterogeneity is not occasional and is frequently linked with specific biological function. Thus, the accurate description of macromolecular conformational transitions is crucial for understanding fundamental mechanisms of life's machinery. We report on a real-time method to predict such transitions by extrapolating from instantaneous eigen motions, computed using the normal mode analysis, to a series of twists. We demonstrate the applicability of our approach to the prediction of a wide range of motions, including large collective opening-closing transitions and conformational changes induced by partner binding. We also highlight particularly difficult cases of very small transitions between crystal and solution structures. Our method guarantees preservation of the protein structure during the transition and allows accessing conformations that are unreachable with classical normal mode analysis. We provide practical solutions to describe localized motions with a few low-frequency modes and to relax some geometrical constraints along the predicted transitions. This work opens the way to the systematic description of protein motions, whatever their degree of collectivity. Our method is freely available as a part of the NOn-Linear rigid Block (NOLB) package.

Docking rigid macrocycles using Convex-PL, AutoDock Vina, and RDKit in the D3R Grand Challenge 4.

The D3R Grand Challenge 4 provided a brilliant opportunity to test macrocyclic docking protocols on a diverse high-quality experimental data. We participated in both pose and affinity prediction exercises. Overall, we aimed to use an automated structure-based docking pipeline built around a set of tools developed in our team. This exercise again demonstrated a crucial importance of the correct local ligand geometry for the overall success of docking. Starting from the second part of the pose prediction stage, we developed a stable pipeline for sampling macrocycle conformers. This resulted in the subangstrom average precision of our pose predictions. In the affinity prediction exercise we obtained average results. However, we could improve these when using docking poses submitted by the best predictors. Our docking tools including the Convex-PL scoring function are available at https://team.inria.fr/nano-d/software/.

Assessment of chemical-crosslink-assisted protein structure modeling in CASP13.

With the advance of experimental procedures obtaining chemical crosslinking information is becoming a fast and routine practice. Information on crosslinks can greatly enhance the accuracy of protein structure modeling. Here, we review the current state of the art in modeling protein structures with the assistance of experimentally determined chemical crosslinks within the framework of the 13th meeting of Critical Assessment of Structure Prediction approaches. This largest-to-date blind assessment reveals benefits of using data assistance in difficult to model protein structure prediction cases. However, in a broader context, it also suggests that with the unprecedented advance in accuracy to predict contacts in recent years, experimental crosslinks will be useful only if their specificity and accuracy further improved and they are better integrated into computational workflows.

Small angle X-ray scattering-assisted protein structure prediction in CASP13 and emergence of solution structure differences.

Small angle X-ray scattering (SAXS) measures comprehensive distance information on a protein's structure, which can constrain and guide computational structure prediction algorithms. Here, we evaluate structure predictions of 11 monomeric and oligomeric proteins for which SAXS data were collected and provided to predictors in the 13th round of the Critical Assessment of protein Structure Prediction (CASP13). The category for SAXS-assisted predictions made gains in certain areas for CASP13 compared to CASP12. Improvements included higher quality data with size exclusion chromatography-SAXS (SEC-SAXS) and better selection of targets and communication of results by CASP organizers. In several cases, we can track improvements in model accuracy with use of SAXS data. For hard multimeric targets where regular folding algorithms were unsuccessful, SAXS data helped predictors to build models better resembling the global shape of the target. For most models, however, no significant improvement in model accuracy at the domain level was registered from use of SAXS data, when rigorously comparing SAXS-assisted models to the best regular server predictions. To promote future progress in this category, we identify successes, challenges, and opportunities for improved strategies in prediction, assessment, and communication of SAXS data to predictors. An important observation is that, for many targets, SAXS data were inconsistent with crystal structures, suggesting that these proteins adopt different conformation(s) in solution. This CASP13 result, if representative of PDB structures and future CASP targets, may have substantive implications for the structure training databases used for machine learning, CASP, and use of prediction models for biology.

Blind prediction of homo- and hetero-protein complexes: The CASP13-CAPRI experiment.

We present the results for CAPRI Round 46, the third joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of 20 targets including 14 homo-oligomers and 6 heterocomplexes. Eight of the homo-oligomer targets and one heterodimer comprised proteins that could be readily modeled using templates from the Protein Data Bank, often available for the full assembly. The remaining 11 targets comprised 5 homodimers, 3 heterodimers, and two higher-order assemblies. These were more difficult to model, as their prediction mainly involved "ab-initio" docking of subunit models derived from distantly related templates. A total of ~30 CAPRI groups, including 9 automatic servers, submitted on average ~2000 models per target. About 17 groups participated in the CAPRI scoring rounds, offered for most targets, submitting ~170 models per target. The prediction performance, measured by the fraction of models of acceptable quality or higher submitted across all predictors groups, was very good to excellent for the nine easy targets. Poorer performance was achieved by predictors for the 11 difficult targets, with medium and high quality models submitted for only 3 of these targets. A similar performance "gap" was displayed by scorer groups, highlighting yet again the unmet challenge of modeling the conformational changes of the protein components that occur upon binding or that must be accounted for in template-based modeling. Our analysis also indicates that residues in binding interfaces were less well predicted in this set of targets than in previous Rounds, providing useful insights for directions of future improvements.

KAP1 is an antiparallel dimer with a functional asymmetry.

KAP1 (KRAB domain-associated protein 1) plays a fundamental role in regulating gene expression in mammalian cells by recruiting different transcription factors and altering the chromatin state. In doing so, KAP1 acts both as a platform for macromolecular interactions and as an E3 small ubiquitin modifier ligase. This work sheds light on the overall organization of the full-length protein combining solution scattering data, integrative modeling, and single-molecule experiments. We show that KAP1 is an elongated antiparallel dimer with an asymmetry at the C-terminal domains. This conformation is consistent with the finding that the Really Interesting New Gene (RING) domain contributes to KAP1 auto-SUMOylation. Importantly, this intrinsic asymmetry has key functional implications for the KAP1 network of interactions, as the heterochromatin protein 1 (HP1) occupies only one of the two putative HP1 binding sites on the KAP1 dimer, resulting in an unexpected stoichiometry, even in the context of chromatin fibers.