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

Assessment of Two Restraint Potentials for Coarse-Grained Chemical-Cross-Link-Assisted Modeling of Protein Structures

The influence of distance restraints from chemical cross-link mass spectroscopy (XL-MS) on the quality of protein structures modeled with the coarse-grained UNRES force field was assessed by using a protocol based on multiplexed replica exchange molecular dynamics, in which both simulated and experimental cross-link restraints were employed, for 23 small proteins. Six cross-links with upper distance boundaries from 4 Å to 12 Å (azido benzoic acid succinimide (ABAS), triazidotriazine (TATA), succinimidyldiazirine (SDA), disuccinimidyl adipate (DSA), disuccinimidyl glutarate (DSG), and disuccinimidyl suberate (BS3)) and two types of restraining potentials ((i) simple flat-bottom Lorentz-like potentials dependent on side chain distance (all cross-links) and (ii) distance- and orientation-dependent potentials determined based on molecular dynamics simulations of model systems (DSA, DSG, BS3, and SDA)) were considered. The Lorentz-like potentials with properly set parameters were found to produce a greater number of higher-quality models compared to unrestrained simulations than the MD-based potentials, because the latter can force too long distances between side chains. Therefore, the flat-bottom Lorentz-like potentials are recommended to represent cross-link restraints. It was also found that significant improvement of model quality upon the introduction of cross-link restraints is obtained when the sum of differences of indices of cross-linked residues exceeds 150.

Protein structure prediction with in-cell photo-crosslinking mass spectrometry and deep learning

While AlphaFold2 can predict accurate protein structures from the primary sequence, challenges remain for proteins that undergo conformational changes or for which few homologous sequences are known. Here we introduce AlphaLink, a modified version of the AlphaFold2 algorithm that incorporates experimental distance restraint information into its network architecture. By employing sparse experimental contacts as anchor points, AlphaLink improves on the performance of AlphaFold2 in predicting challenging targets. We confirm this experimentally by using the noncanonical amino acid photo-leucine to obtain information on residue-residue contacts inside cells by crosslinking mass spectrometry. The program can predict distinct conformations of proteins on the basis of the distance restraints provided, demonstrating the value of experimental data in driving protein structure prediction. The noise-tolerant framework for integrating data in protein structure prediction presented here opens a path to accurate characterization of protein structures from in-cell data.

Pragmatic Coarse-Graining of Proteins: Models and Applications

The molecular details involved in the folding, dynamics, organization, and interaction of proteins with other molecules are often difficult to assess by experimental techniques. Consequently, computational models play an ever-increasing role in the field. However, biological processes involving large-scale protein assemblies or long time scale dynamics are still computationally expensive to study in atomistic detail. For these applications, employing coarse-grained (CG) modeling approaches has become a key strategy. In this Review, we provide an overview of what we call pragmatic CG protein models, which are strategies combining, at least in part, a physics-based implementation and a top-down experimental approach to their parametrization. In particular, we focus on CG models in which most protein residues are represented by at least two beads, allowing these models to retain some degree of chemical specificity. A description of the main modern pragmatic protein CG models is provided, including a review of the most recent applications and an outlook on future perspectives in the field.

Hybrid computational methods combining experimental information with molecular dynamics

A goal of structural biology is to understand how macromolecules carry out their biological roles by identifying their metastable states, mechanisms of action, pathways leading to conformational changes, and the thermodynamic and kinetic relationships between those states. Integrative modeling brings structural insights into systems where traditional structure determination approaches cannot help. We focus on the synergies and challenges of integrative modeling combining experimental data with molecular dynamics simulations.

The Cdc48 N-terminal domain has a molecular switch that mediates the Npl4-Ufd1-Cdc48 complex formation

Cdc48 (VCP/p97) is a major AAA-ATPase involved in protein quality control, along with its canonical cofactors Ufd1 and Npl4 (UN). Here, we present novel structural insights into the interactions within the Cdc48-Npl4-Ufd1 ternary complex. Using integrative modeling, we combine subunit structures with crosslinking mass spectrometry (XL-MS) to map the interaction between Npl4 and Ufd1, alone and in complex with Cdc48. We describe the stabilization of the UN assembly upon binding with the N-terminal-domain (NTD) of Cdc48 and identify a highly conserved cysteine, C115, at the Cdc48-Npl4-binding interface which is central to the stability of the Cdc48-Npl4-Ufd1 complex. Mutation of Cys115 to serine disrupts the interaction between Cdc48-NTD and Npl4-Ufd1 and leads to a moderate decrease in cellular growth and protein quality control in yeast. Our results provide structural insight into the architecture of the Cdc48-Npl4-Ufd1 complex as well as its in vivo implications.

Improvements and new functionalities of UNRES server for coarse-grained modeling of protein structure, dynamics, and interactions

In this paper we report the improvements and extensions of the UNRES server (https://unres-server.chem.ug.edu.pl) for physics-based simulations with the coarse-grained UNRES model of polypeptide chains. The improvements include the replacement of the old code with the recently optimized one and adding the recent scale-consistent variant of the UNRES force field, which performs better in the modeling of proteins with the β and the α+β structures. The scope of applications of the package was extended to data-assisted simulations with restraints from nuclear magnetic resonance (NMR) and chemical crosslink mass-spectroscopy (XL-MS) measurements. NMR restraints can be input in the NMR Exchange Format (NEF), which has become a standard. Ambiguous NMR restraints are handled without expert intervention owing to a specially designed penalty function. The server can be used to run smaller jobs directly or to prepare input data to run larger production jobs by using standalone installations of UNRES.

Modeling the Structure, Dynamics, and Transformations of Proteins with the UNRES Force Field

The physics-based united-residue (UNRES) model of proteins ( www.unres.pl ) has been designed to carry out large-scale simulations of protein folding. The force field has been derived and parameterized based on the principles of statistical-mechanics, which makes it independent of structural databases and applicable to treat nonstandard situations such as, proteins that contain D-amino-acid residues. Powered by Langevin dynamics and its replica-exchange extensions, UNRES has found a variety of applications, including ab initio and database-assisted protein-structure prediction, simulating protein-folding pathways, exploring protein free-energy landscapes, and solving biological problems. This chapter provides a summary of UNRES and a guide for potential users regarding the application of the UNRES package in a variety of research tasks.

Recent Developments in Data-Assisted Modeling of Flexible Proteins

Many proteins can fold into well-defined conformations. However, intrinsically-disordered proteins (IDPs) do not possess a defined structure. Moreover, folded multi-domain proteins often digress into alternative conformations. Collectively, the conformational dynamics enables these proteins to fulfill specific functions. Thus, most experimental observables are averaged over the conformations that constitute an ensemble. In this article, we review the recent developments in the concept and methods for the determination of the dynamic structures of flexible peptides and proteins. In particular, we describe ways to extract information from nuclear magnetic resonance small-angle X-ray scattering (SAXS), and chemical cross-linking coupled with mass spectroscopy (XL-MS) measurements. All these techniques can be used to obtain ensemble-averaged restraints or to re-weight the simulated conformational ensembles.

Critical assessment of methods of protein structure prediction (CASP)-Round XIV

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.

Pseudopotentials for coarse-grained cross-link-assisted modeling of protein structures

Pseudopotentials for the chemical cross-links comprising the glutamic- and aspartic-acid side chains bridged with adipic- (ADH) or pimelic-acid hydrazide (PDH), and the lysine side chains bridged with glutaric (BS² G) or suberic acid (BS³ ) for coarse-grained cross-link-assisted simulations were determined by canonical molecular dynamics with the Amber14sb force field. The potentials depend on the distance between side-chain ends and on side-chain orientation, this preventing from making cross-link contacts across the globule in simulations. The potentials were implemented in the UNRES coarse-grained force field and their effect on the quality of models was assessed with 11 monomeric and 1 dimeric proteins, using synthetic or experimental cross-link data. Simulations with the new potentials resulted in improvement of the generated models compared to unrestrained simulations in more instances compared to those with the statistical potentials.

ESCASA: Analytical estimation of atomic coordinates from coarse-grained geometry for nuclear-magnetic-resonance-assisted protein structure modeling. I. Backbone and Hβ protons

A method for the estimation of coordinates of atoms in proteins from coarse-grained geometry by simple analytical formulas (ESCASA), for use in nuclear-magnetic-resonance (NMR) data-assisted coarse-grained simulations of proteins is proposed. In this paper, the formulas for the backbone H^α and amide (H^N ) protons, and the side-chain H^β protons, given the C^α -trace, have been derived and parameterized, by using the interproton distances calculated from a set of 140 high-resolution non-homologous protein structures. The mean standard deviation over all types of proton pairs in the set was 0.44 Å after fitting. Validation against a set of 41 proteins with NMR-determined structures, which were not considered in parameterization, resulted in average standard deviation from average proton-proton distances of the NMR-determined structures of 0.25 Å, compared to 0.21 Å obtained with the PULCHRA all-atom-chain reconstruction algorithm and to the 0.12 Å standard deviation of the average-structure proton-proton distance of NMR-determined ensembles. The formulas provide analytical forces and can, therefore, be used in coarse-grained molecular dynamics.

Structural interactions define assembly adapter function of a type II secretion system pseudopilin

The type IV filament superfamily comprises widespread membrane-associated polymers in prokaryotes. The type II secretion system (T2SS), a virulence pathway in many pathogens, belongs to this superfamily. A knowledge gap in understanding of the T2SS is the molecular role of a small "pseudopilin" protein. Using multiple biophysical techniques, we have deciphered how this missing component of the Xcp T2SS architecture is structurally integrated, and thereby unlocked its function. We demonstrate that low-abundance XcpH is the adapter that bridges a trimeric initiating tip complex, XcpIJK, with a periplasmic filament of XcpG subunits. Each pseudopilin protein caps an XcpG protofilament in an overall pseudopilus compatible with dimensions of the periplasm and the outer membrane-spanning secretin through which substrates pass. Unexpectedly, to fulfill its adapter function, the XcpH N-terminal helix must be unwound, a property shared with XcpG subunits. We provide an experimentally validated three-dimensional structural model of a complete type IV filament.

Decellularized nerve extracellular matrix/chitosan crosslinked by genipin to prepare a moldable nerve repair material

Decellularized nerve extracellular matrix (NECM) composited with chitosan are moldable materials suitable for spinal cord repair. But the rapid biodegradation of the materials may interrupt neural tissue reconstruction in vivo. To improve the stability of the materials, the materials produced by NECM and chitosan hydrogels were crosslinked by genipine, glutaraldehyde or ultraviolet ray. Physicochemical property, degradation and biocompatibility of materials crosslinked by genipin, glutaraldehyde or ultraviolet ray were evaluated. The scaffold crosslinked by genipin possessed a porous structure, and the porosity ratio was 89.07 + 4.90%, the average diameter of pore was 85.32 + 5.34 μm. The crosslinked degree of the scaffold crosslinked by genipin and glutaraldehyde was 75.13 ± 4.87%, 71.25 ± 5.06% respectively; Uncrosslinked scaffold disintegrated when immerged in distilled water while the scaffold crosslinked by genipin and glutaraldehyde group retained their integrity. The scaffold crosslinked by genipin has better water absorption, water retention and anti-enzymatic hydrolysis ability than the other three groups. Cell cytotoxicity showed that the cytotoxicity of scaffold crosslinked by genipin was lower than that crosslinked by glutaraldehyde. The histocompatibility of scaffold crosslinked by genipin was also better than glutaraldehyde group. More cells grew well in the scaffold crosslinked by genipin when co-cultured with L929 cells. The decellularized nerve extracellular matrix/chitosan scaffold crosslinked by the genipin has good mechanical properties, micro structure and biocompatibility, which is an ideal scaffold for the spinal cord tissue engineering.

John PC, Bomble YJ. Advances in integrative structural biology: Towards understanding protein complexes in their cellular context

Microorganisms rely on protein interactions to transmit signals, react to stimuli, and grow. One of the best ways to understand these protein interactions is through structural characterization. However, in the past, structural knowledge was limited to stable, high-affinity complexes that could be crystallized. Recent developments in structural biology have revolutionized how protein interactions are characterized. The combination of multiple techniques, known as integrative structural biology, has provided insight into how large protein complexes interact in their native environment. In this mini-review, we describe the past, present, and potential future of integrative structural biology as a tool for characterizing protein interactions in their cellular context.

Tightening the Crosslinking Distance Restraints for Better Resolution of Protein Structure and Dynamics

Chemical crosslinking coupled with mass spectrometry (CXMS) has been increasingly used in structural biology. CXMS distance restraints are usually applied to Cα or Cβ atoms of the crosslinked residues, with upper bounds typically over 20 Å. The incorporation of loose CXMS restraints only marginally improves the resolution of the calculated structures. Here, we present a revised format of CXMS distance restraints, which works by first modifying the crosslinked residue with a rigid extension derived from the crosslinker. With the flexible side chain explicitly represented, the reformatted restraint can be applied to the modification group instead, with an upper bound of 6 Å or less. The short distance restraint can be represented and back-calculated simply with a straight line. The use of tighter restraints not only afford better-resolved structures but also uncover protein dynamics. Together, our approach enables more information extracted from the CXMS data.

The Conformational Preference of Chemical Cross-linkers Determines the Cross-linking Probability of Reactive Protein Residues

Chemical cross-linking mass spectrometry (XLMS) is an emerging technique in structural biology. Providing the cross-linked peptides are identified by mass spectrometry with high confidence, a distance restraint can be applied between the two reactive protein residues, with the upper bound corresponding to the maximal span of the cross-linker. However, as the upper bound is typically over 20 Å, cross-link distance restraints are unrestrictive and provide a marginal improvement in protein structural refinement. Here we analyze the experimental cross-links for lysine or acidic residues and show that the distribution of Cβ-Cβ' distances can be described with two overlapping Gaussian species. In addition to the pairwise occurrence probability of the reactive protein residues, we show that the distribution profile of the cross-link distances is determined by the intrinsic conformational propensity of the cross-linker. The cross-linker prefers either a compact or extended conformation and, once attached to a reactive protein residue, predominantly an extended conformation. Consequently, the long-distance Gaussian species occurs at a much higher probability than the short-distance species in the observed cross-links. Together, the probabilistic distribution of the cross-link distance allows the construction of a more restrictive restraint for structural modeling and better use of the XLMS data.

Protein Footprinting: Auxiliary Engine to Power the Structural Biology Revolution

Structural biology is entering an exciting time where many new high-resolution structures of large complexes and membrane proteins are determined regularly. These advances have been driven by over fifteen years of technology advancements, first in macromolecular crystallography, and recently in Cryo-electron microscopy. These structures are allowing detailed questions about functional mechanisms of the structures, and the biology enabled by these structures, to be addressed for the first time. At the same time, mass spectrometry technologies for protein structure analysis, "footprinting" studies, have improved their sensitivity and resolution dramatically and can provide detailed sub-peptide and residue level information for validating structures and interactions or understanding the dynamics of structures in the context of ligand binding or assembly. In this perspective, we review the use of protein footprinting to extend our understanding of macromolecular systems, particularly for systems challenging for analysis by other techniques, such as intrinsically disordered proteins, amyloidogenic proteins, and other proteins/complexes so far recalcitrant to existing methods. We also illustrate how the availability of high-resolution structural information can be a foundation for a suite of hybrid approaches to divine structure-function relationships beyond what individual techniques can deliver.

Improved Consensus-Fragment Selection in Template-Assisted Prediction of Protein Structures with the UNRES Force Field in CASP13

The method for protein-structure prediction, which combines the physics-based coarse-grained UNRES force field with knowledge-based modeling, has been developed further and tested in the 13th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP13). The method implements restraints from the consensus fragments common to server models. In this work, the server models to derive fragments have been chosen on the basis of quality assessment; a fully automatic fragment-selection procedure has been introduced, and Dynamic Fragment Assembly pseudopotentials have been fully implemented. The Global Distance Test Score (GDT_TS), averaged over our “Model 1” predictions, increased by over 10 units with respect to CASP12 for the free-modeling category to reach 40.82. Our “Model 1” predictions ranked 20 and 14 for all and free-modeling targets, respectively (upper 20.2% and 14.3% of all models submitted to CASP13 in these categories, respectively), compared to 27 (upper 21.1%) and 24 (upper 18.9%) in CASP12, respectively. For oligomeric targets, the Interface Patch Similarity (IPS) and Interface Contact Similarity (ICS) averaged over our best oligomer models increased from 0.28 to 0.36 and from 12.4 to 17.8, respectively, from CASP12 to CASP13, and top-ranking models of 2 targets (H0968 and T0997o) were obtained (none in CASP12). The improvement of our method in CASP13 over CASP12 was ascribed to the combined effect of the overall enhancement of server-model quality, our success in selecting server models and fragments to derive restraints, and improvements of the restraint and potential-energy functions.

Critical assessment of methods of protein structure prediction (CASP)-Round XIII

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.

Assessment of protein assembly prediction in CASP13

We present the assembly category assessment in the 13th edition of the CASP community-wide experiment. For the second time, protein assemblies constitute an independent assessment category. Compared to the last edition we see a clear uptake in participation, more oligomeric targets released, and consistent, albeit modest, improvement of the predictions quality. Looking at the tertiary structure predictions, we observe that ignoring the oligomeric state of the targets hinders modeling success. We also note that some contact prediction groups successfully predicted homomeric interfacial contacts, though it appears that these predictions were not used for assembly modeling. Homology modeling with sizeable human intervention appears to form the basis of the assembly prediction techniques in this round of CASP. Future developments should see more integrated approaches where subunits are modeled in the context of the assemblies they form.