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

Advances in crosslinking chemistry and proximity-enabled strategies: deciphering protein complexes and interactions

Mass spectrometry, coupled with innovative crosslinking techniques to decode protein conformations and interactions through uninterrupted signal connections, has undergone remarkable progress in recent years. It is crucial to develop selective crosslinking reagents that minimally disrupt protein structure and dynamics, providing insights into protein network regulation and biological functions. Compared to traditional crosslinkers, new bifunctional chemical crosslinkers exhibit high selectivity and specificity in connecting proximal amino acid residues, resulting in stable molecular crosslinked products. The conjugation with specific amino acid residues like lysine, cysteine, arginine and tyrosine expands the XL-MS toolbox, enabling more precise modeling of target substrates and leading to improved data quality and reliability. Another emerging crosslinking method utilizes unnatural amino acids (UAAs) derived from proximity-enabled reactivity with specific amino acids or sulfur-fluoride exchange (SuFEx) reactions with nucleophilic residues. These UAAs are genetically encoded into proteins for the formation of specific covalent bonds. This technique combines the benefits of genetic encoding for live cell compatibility with chemical crosslinking, providing a valuable method for capturing transient and weak protein-protein interactions (PPIs) for mapping PPI coordinates and improving the pharmacological properties of proteins. With continued advancements in technology and applications, crosslinking mass spectrometry is poised to play an increasingly significant role in guiding our understanding of protein dynamics and function in the future.

Native Mass Spectrometry Reveals Binding Interactions of SARS-CoV-2 PLpro with Inhibitors and Cellular Targets

Here we used native mass spectrometry (native MS) to probe a SARS-CoV protease, PLpro, which plays critical roles in coronavirus disease by affecting viral protein production and antagonizing host antiviral responses. Ultraviolet photodissociation (UVPD) and variable temperature electrospray ionization (vT ESI) were used to localize binding sites of PLpro inhibitors and revealed the stabilizing effects of inhibitors on protein tertiary structure. We compared PLpro from SARS-CoV-1 and SARS-CoV-2 in terms of inhibitor and ISG15 interactions to discern possible differences in protease function. A PLpro mutant lacking a single cysteine was used to localize inhibitor binding, and thermodynamic measurements revealed that inhibitor PR-619 stabilized the folded PLpro structure. These results will inform further development of PLpro as a therapeutic target against SARS-CoV-2 and other emerging coronaviruses.

Partial wrapping of single-stranded DNA by replication protein A and modulation through phosphorylation

Single-stranded DNA (ssDNA) intermediates which emerge during DNA metabolic processes are shielded by replication protein A (RPA). RPA binds to ssDNA and acts as a gatekeeper to direct the ssDNA towards downstream DNA metabolic pathways with exceptional specificity. Understanding the mechanistic basis for such RPA-dependent functional specificity requires knowledge of the structural conformation of ssDNA when RPA-bound. Previous studies suggested a stretching of ssDNA by RPA. However, structural investigations uncovered a partial wrapping of ssDNA around RPA. Therefore, to reconcile the models, in this study, we measured the end-to-end distances of free ssDNA and RPA-ssDNA complexes using single-molecule FRET and double electron-electron resonance (DEER) spectroscopy and found only a small systematic increase in the end-to-end distance of ssDNA upon RPA binding. This change does not align with a linear stretching model but rather supports partial wrapping of ssDNA around the contour of DNA binding domains of RPA. Furthermore, we reveal how phosphorylation at the key Ser-384 site in the RPA70 subunit provides access to the wrapped ssDNA by remodeling the DNA-binding domains. These findings establish a precise structural model for RPA-bound ssDNA, providing valuable insights into how RPA facilitates the remodeling of ssDNA for subsequent downstream processes.

Xplor-NIH: Better parameters and protocols for NMR protein structure determination

The present work describes an update to the protein covalent geometry and atomic radii parameters in the Xplor-NIH biomolecular structure determination package. In combination with an improved treatment of selected non-bonded interactions between atoms three bonds apart, such as those involving methyl hydrogens, and a previously developed term that affects the system's gyration volume, the new parameters are tested using structure calculations on 30 proteins with restraints derived from nuclear magnetic resonance data. Using modern structure validation criteria, including several formally adopted by the Protein Data Bank, and a clear measure of structural accuracy, the results show superior performance relative to previous Xplor-NIH implementations. Additionally, the Xplor-NIH structures compare favorably against originally determined NMR models.

Wrapping of single-stranded DNA by Replication Protein A and modulation through phosphorylation

Single-stranded DNA (ssDNA) intermediates, which emerge during DNA metabolic processes are shielded by Replication Protein A (RPA). RPA binds to ssDNA and acts as a gatekeeper, directing the ssDNA towards downstream DNA metabolic pathways with exceptional specificity. Understanding the mechanistic basis for such RPA-dependent specificity requires a comprehensive understanding of the structural conformation of ssDNA when bound to RPA. Previous studies suggested a stretching of ssDNA by RPA. However, structural investigations uncovered a partial wrapping of ssDNA around RPA. Therefore, to reconcile the models, in this study, we measured the end-to-end distances of free ssDNA and RPA-ssDNA complexes using single-molecule FRET and Double Electron-Electron Resonance (DEER) spectroscopy and found only a small systematic increase in the end-to-end distance of ssDNA upon RPA binding. This change does not align with a linear stretching model but rather supports partial wrapping of ssDNA around the contour of DNA binding domains of RPA. Furthermore, we reveal how phosphorylation at the key Ser-384 site in the RPA70 subunit provides access to the wrapped ssDNA by remodeling the DNA-binding domains. These findings establish a precise structural model for RPA-bound ssDNA, providing valuable insights into how RPA facilitates the remodeling of ssDNA for subsequent downstream processes.

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.

Yeast Rad52 is a homodecamer and possesses BRCA2-like bipartite Rad51 binding modes

Homologous recombination (HR) is an essential double-stranded DNA break repair pathway. In HR, Rad52 facilitates the formation of Rad51 nucleoprotein filaments on RPA-coated ssDNA. Here, we decipher how Rad52 functions using single-particle cryo-electron microscopy and biophysical approaches. We report that Rad52 is a homodecameric ring and each subunit possesses an ordered N-terminal and disordered C-terminal half. An intrinsic structural asymmetry is observed where a few of the C-terminal halves interact with the ordered ring. We describe two conserved charged patches in the C-terminal half that harbor Rad51 and RPA interacting motifs. Interactions between these patches regulate ssDNA binding. Surprisingly, Rad51 interacts with Rad52 at two different bindings sites: one within the positive patch in the disordered C-terminus and the other in the ordered ring. We propose that these features drive Rad51 nucleation onto a single position on the DNA to promote formation of uniform pre-synaptic Rad51 filaments in HR.

A deep learning model for predicting optimal distance range in crosslinking mass spectrometry data

Macromolecular assemblies play an important role in all cellular processes. While there has recently been significant progress in protein structure prediction based on deep learning, large protein complexes cannot be predicted with these approaches. The integrative structure modeling approach characterizes multi-subunit complexes by computational integration of data from fast and accessible experimental techniques. Crosslinking mass spectrometry is one such technique that provides spatial information about the proximity of crosslinked residues. One of the challenges in interpreting crosslinking datasets is designing a scoring function that, given a structure, can quantify how well it fits the data. Most approaches set an upper bound on the distance between Cα atoms of crosslinked residues and calculate a fraction of satisfied crosslinks. However, the distance spanned by the crosslinker greatly depends on the neighborhood of the crosslinked residues. Here, we design a deep learning model for predicting the optimal distance range for a crosslinked residue pair based on the structures of their neighborhoods. We find that our model can predict the distance range with the area under the receiver-operator curve of 0.86 and 0.7 for intra- and inter-protein crosslinks, respectively. Our deep scoring function can be used in a range of structure modeling applications.

Homodecameric Rad52 promotes single-position Rad51 nucleation in homologous recombination

Homologous recombination (HR) is a pathway for the accurate repair of double-stranded DNA breaks. These breaks are resected to yield single-stranded DNA (ssDNA) that are coated by Replication Protein A (RPA). Saccharomyces cerevisiae Rad52 is a mediator protein that promotes HR by facilitating formation of Rad51 nucleoprotein filaments on RPA-coated ssDNA. Canonically, Rad52 has been described to function by displacing RPA to promote Rad51 binding. However, in vitro, Rad51 readily forms a filament by displacing RPA in the absence of Rad52. Yet, in vivo, Rad52 is essential for HR. Here, we resolve how Rad52 functions as a mediator using single-particle cryo-electron microscopy and biophysical approaches. We show that Rad52 functions as a homodecamer and catalyzes single-position nucleation of Rad51. The N-terminal half of Rad52 is a well-ordered ring, while the C-terminal half is disordered. An intrinsic asymmetry within Rad52 is observed, where one or a few of the C-terminal halves interact with the ordered N-terminal ring. Within the C-terminal half, we identify two conserved charged patches that harbor the Rad51 and RPA interacting motifs. Interactions between these two charged patches regulate a ssDNA binding. These features drive Rad51 binding to a single position on the Rad52 decameric ring. We propose a Rad52 catalyzed single-position nucleation model for the formation of pre-synaptic Rad51 filaments in HR.

Evolutionary conservation of the structure and function of meiotic Rec114-Mei4 and Mer2 complexes

Meiosis-specific Rec114-Mei4 and Mer2 complexes are thought to enable Spo11-mediated DNA double-strand break (DSB) formation through a mechanism that involves DNA-dependent condensation. However, the structure, molecular properties, and evolutionary conservation of Rec114-Mei4 and Mer2 are unclear. Here, we present AlphaFold models of Rec114-Mei4 and Mer2 complexes supported by nuclear magnetic resonance (NMR) spectroscopy, small-angle X-ray scattering (SAXS), and mutagenesis. We show that dimers composed of the Rec114 C terminus form α-helical chains that cup an N-terminal Mei4 α helix, and that Mer2 forms a parallel homotetrameric coiled coil. Both Rec114-Mei4 and Mer2 bind preferentially to branched DNA substrates, indicative of multivalent protein-DNA interactions. Indeed, the Rec114-Mei4 interaction domain contains two DNA-binding sites that point in opposite directions and drive condensation. The Mer2 coiled-coil domain bridges coaligned DNA duplexes, likely through extensive electrostatic interactions along the length of the coiled coil. Finally, we show that the structures of Rec114-Mei4 and Mer2 are conserved across eukaryotes, while DNA-binding properties vary significantly. This work provides insights into the mechanism whereby Rec114-Mei4 and Mer2 complexes promote the assembly of the meiotic DSB machinery and suggests a model in which Mer2 condensation is the essential driver of assembly, with the DNA-binding activity of Rec114-Mei4 playing a supportive role.

Small-Molecule Inhibition of Androgen Receptor Dimerization as a Strategy against Prostate Cancer

The clinically used androgen receptor (AR) antagonists for the treatment of prostate cancer (PCa) are all targeting the AR ligand binding pocket (LBP), resulting in various drug-resistant problems. Therefore, a new strategy to combat PCa is urgently needed. Enlightened by the gain-of-function mutations of androgen insensitivity syndrome, we discovered for the first time small-molecule antagonists toward a prospective pocket on the AR dimer interface named the dimer interface pocket (DIP) via molecular dynamics (MD) simulation, structure-based virtual screening, structure-activity relationship exploration, and bioassays. The first-in-class antagonist M17-B15 targeting the DIP is capable of effectively disrupting AR self-association, thereby suppressing AR signaling. Furthermore, M17-B15 exhibits extraordinary anti-PCa efficacy in vitro and also in mouse xenograft tumor models, demonstrating that AR dimerization disruption by small molecules targeting the DIP is a novel and valid strategy against PCa.

Disulfide bond and crosslinking analyses reveal inter-domain interactions that contribute to the rigidity of placental malaria VAR2CSA structure and formation of CSA binding channel

VAR2CSA, a multidomain Plasmodium falciparum protein, mediates the adherence of parasite-infected red blood cells to chondroitin 4-sulfate (C4S) in the placenta, contributing to placental malaria. Therefore, detailed understanding of VAR2CSA structure likely help developing strategies to treat placental malaria. The VAR2CSA ectodomain consists of an N-terminal segment (NTS), six Duffy binding-like (DBL) domains, and three interdomains (IDs) present in sequence NTS-DBL1x-ID1-DBL2x-ID2-DBL3x-DBL4ε-ID3-DBL5ε-DBL6ε. Recent electron microscopy studies showed that VAR2CSA is compactly organized into a globular structure containing C4S-binding channel, and that DBL5ε-DBL6ε arm is attached to the NTS-ID3 core structure. However, the structural elements involved in inter-domain interactions that stabilize the VAR2CSA structure remain largely not understood. Here, limited proteolysis and peptide mapping by mass spectrometry showed that VAR2CSA contains several inter-domain disulfide bonds that stabilize its compact structure. Chemical crosslinking-mass spectrometry showed that all IDs interact with DBL4ε; additionally, IDs interact with other DBL domains, demonstrating that IDs are the key structural scaffolds that shape the functional NTS-ID3 core. Ligand binding analysis suggested that NTS considerably restricts the C4S binding. Overall, our study revealed that inter-domain disulfide bonds and interactions between IDs and DBL domains contribute to the stability of VAR2CSA structural architecture and formation of C4S-binding channel.

Unveiling the biological interface of protein complexes by mass spectrometry-coupled methods

Most biomolecules become functional and bioactive by forming protein complexes through interaction with ligands that are diverse in size, shape, and physicochemical properties. In the complex biological milieu, the interaction is ligand-specific, driven by molecular sensing, and involves the recognition of a binding interface localized within a protein structure. Mapping interfaces of protein complexes is a highly sought area of research as it delivers fundamental insights into proteomes and pathology and hence strategies for therapeutics. While X-ray crystallography and electron microscopy remain the gold standard for structural elucidation of protein complexes, their artificial and static analytic nature often produces a non-native interface that otherwise might be negligible or non-existent in a biological environment. Recently, the mass spectrometry-coupled approaches, chemical crosslinking (CLMS) and hydrogen-deuterium exchange (HDMS) have become valuable analytic complements to the traditional techniques. These methods explicitly identify hot residues and motifs embedded in binding interfaces, especially when the interaction is predominantly dynamic, transient, and/or caused by an intrinsically disordered domain. Here, we review the principal role of CLMS and HDMS in protein structural biology with a particular emphasis on the contribution of recent examples to exploring biological interfaces. Additionally, we describe recent studies that utilized these methods to expand our understanding of protein complex formation and the related biological processes, to increase the probability of structure-based drug design.

Application of Multiple Length Cross-linkers to the Characterization of Gaseous Protein Structure

The speed, sensitivity, and tolerance of heterogeneity, as well as the kinetic trapping of solution-like states during electrospray, make native mass spectrometry an attractive method to study protein structure. Increases in the resolution of ion mobility measurements and in mass resolving power and range are leading to the increase of the information content of intact protein measurements and an expanded role of mass spectrometry in structural biology. Herein, a suite of different length noncovalent (sulfonate to positively charged side chain) cross-linkers was introduced via gas-phase ion/ion chemistry and used to determine distance restraints of kinetically trapped gas-phase structures of native-like cytochrome c ions. Electron capture dissociation allowed for the identification of cross-linked sites. Different length linkers resulted in distinct pairs of side chains being linked, supporting the ability of gas-phase cross-linking to be structurally specific. The gas-phase lengths of the cross-linkers were determined by conformational searches and density functional theory, allowing for the interpretation of the cross-links as distance restraints. These distance restraints were used to model gas-phase structures with molecular dynamics simulations, revealing a mixture of structures with similar overall shape/size but distinct features, thereby illustrating the kinetic trapping of multiple native-like solution structures in the gas phase.

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.

Leveraging crosslinking mass spectrometry in structural and cell biology

Crosslinking mass spectrometry (crosslinking-MS) is a versatile tool providing structural insights into protein conformation and protein-protein interactions. Its medium-resolution residue-residue distance restraints have been used to validate protein structures proposed by other methods and have helped derive models of protein complexes by integrative structural biology approaches. The use of crosslinking-MS in integrative approaches is underpinned by progress in estimating error rates in crosslinking-MS data and in combining these data with other information. The flexible and high-throughput nature of crosslinking-MS has allowed it to complement the ongoing resolution revolution in electron microscopy by providing system-wide residue-residue distance restraints, especially for flexible regions or systems. Here, we review how crosslinking-MS information has been leveraged in structural model validation and integrative modeling. Crosslinking-MS has also been a key technology for cell biology studies and structural systems biology where, in conjunction with cryoelectron tomography, it can provide structural and mechanistic insights directly in situ.

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.