Carving out a Glycoside Hydrolase Active Site for Incorporation into a New Protein Scaffold Using Deep Network Hallucination.

Enzymes are indispensable biocatalysts for numerous industrial applications, yet stability, selectivity, and restricted substrate recognition present limitations for their use. Despite the importance of enzyme engineering in overcoming these limitations, success is often challenged by the intricate architecture of enzymes derived from natural sources. Recent advances in computational methods have enabled the de novo design of simplified scaffolds with specific functional sites. Such scaffolds may be advantageous as platforms for enzyme engineering. Here, we present a strategy for the de novo design of a simplified scaffold of an endo-α-N-acetylgalactosaminidase active site, a glycoside hydrolase from the GH101 enzyme family. Using a combination of trRosetta hallucination, iterative cycles of deep-learning-based structure prediction, and ProteinMPNN sequence design, we designed proteins with 290 amino acids incorporating the active site while reducing the molecular weight by over 100 kDa compared to the initial endo-α-N-acetylgalactosaminidase. Of 11 tested designs, six were expressed as soluble monomers, displaying similar or increased thermostabilities compared to the natural enzyme. Despite lacking detectable enzymatic activity, the experimentally determined crystal structures of a representative design closely matched the design with a root-mean-square deviation of 1.0 Å, with most catalytically important side chains within 2.0 Å. The results highlight the potential of scaffold hallucination in designing proteins that may serve as a foundation for subsequent enzyme engineering.

Checkpoint activation by Spd1: a competition-based system relying on tandem disordered PCNA binding motifs.

DNA regulation, replication and repair are processes fundamental to all known organisms and the sliding clamp proliferating cell nuclear antigen (PCNA) is central to all these processes. S-phase delaying protein 1 (Spd1) from S. pombe, an intrinsically disordered protein that causes checkpoint activation by inhibiting the enzyme ribonucleotide reductase, has one of the most divergent PCNA binding motifs known. Using NMR spectroscopy, in vivo assays, X-ray crystallography, calorimetry, and Monte Carlo simulations, an additional PCNA binding motif in Spd1, a PIP-box, is revealed. The two tandemly positioned, low affinity sites exchange rapidly on PCNA exploiting the same binding sites. Increasing or decreasing the binding affinity between Spd1 and PCNA through mutations of either motif compromised the ability of Spd1 to cause checkpoint activation in yeast. These results pinpoint a role for PCNA in Spd1-mediated checkpoint activation and suggest that its tandemly positioned short linear motifs create a neatly balanced competition-based system, involving PCNA, Spd1 and the small ribonucleotide reductase subunit, Suc22R2. Similar mechanisms may be relevant in other PCNA binding ligands where divergent binding motifs so far have gone under the PIP-box radar.

The structure of pathogenic huntingtin exon 1 defines the bases of its aggregation propensity.

Huntington's disease is a neurodegenerative disorder caused by a CAG expansion in the first exon of the HTT gene, resulting in an extended polyglutamine (poly-Q) tract in huntingtin (httex1). The structural changes occurring to the poly-Q when increasing its length remain poorly understood due to its intrinsic flexibility and the strong compositional bias. The systematic application of site-specific isotopic labeling has enabled residue-specific NMR investigations of the poly-Q tract of pathogenic httex1 variants with 46 and 66 consecutive glutamines. Integrative data analysis reveals that the poly-Q tract adopts long α-helical conformations propagated and stabilized by glutamine side chain to backbone hydrogen bonds. We show that α-helical stability is a stronger signature in defining aggregation kinetics and the structure of the resulting fibrils than the number of glutamines. Our observations provide a structural perspective of the pathogenicity of expanded httex1 and pave the way to a deeper understanding of poly-Q-related diseases.

Mechanism of human Lig1 regulation by PCNA in Okazaki fragment sealing.

During lagging strand synthesis, DNA Ligase 1 (Lig1) cooperates with the sliding clamp PCNA to seal the nicks between Okazaki fragments generated by Pol δ and Flap endonuclease 1 (FEN1). We present several cryo-EM structures combined with functional assays, showing that human Lig1 recruits PCNA to nicked DNA using two PCNA-interacting motifs (PIPs) located at its disordered N-terminus (PIPN-term) and DNA binding domain (PIPDBD). Once Lig1 and PCNA assemble as two-stack rings encircling DNA, PIPN-term is released from PCNA and only PIPDBD is required for ligation to facilitate the substrate handoff from FEN1. Consistently, we observed that PCNA forms a defined complex with FEN1 and nicked DNA, and it recruits Lig1 to an unoccupied monomer creating a toolbelt that drives the transfer of DNA to Lig1. Collectively, our results provide a structural model on how PCNA regulates FEN1 and Lig1 during Okazaki fragments maturation.

A glutamine-based single α-helix scaffold to target globular proteins.

The binding of intrinsically disordered proteins to globular ones can require the folding of motifs into α-helices. These interactions offer opportunities for therapeutic intervention but their modulation with small molecules is challenging because they bury large surfaces. Linear peptides that display the residues that are key for binding can be targeted to globular proteins when they form stable helices, which in most cases requires their chemical modification. Here we present rules to design peptides that fold into single α-helices by instead concatenating glutamine side chain to main chain hydrogen bonds recently discovered in polyglutamine helices. The resulting peptides are uncharged, contain only natural amino acids, and their sequences can be optimized to interact with specific targets. Our results provide design rules to obtain single α-helices for a wide range of applications in protein engineering and drug design.

Mg2+-dependent conformational equilibria in CorA and an integrated view on transport regulation.

The CorA family of proteins regulates the homeostasis of divalent metal ions in many bacteria, archaea, and eukaryotic mitochondria, making it an important target in the investigation of the mechanisms of transport and its functional regulation. Although numerous structures of open and closed channels are now available for the CorA family, the mechanism of the transport regulation remains elusive. Here, we investigated the conformational distribution and associated dynamic behaviour of the pentameric Mg2+ channel CorA at room temperature using small-angle neutron scattering (SANS) in combination with molecular dynamics (MD) simulations and solid-state nuclear magnetic resonance spectroscopy (NMR). We find that neither the Mg2+-bound closed structure nor the Mg2+-free open forms are sufficient to explain the average conformation of CorA. Our data support the presence of conformational equilibria between multiple states, and we further find a variation in the behaviour of the backbone dynamics with and without Mg2+. We propose that CorA must be in a dynamic equilibrium between different non-conducting states, both symmetric and asymmetric, regardless of bound Mg2+ but that conducting states become more populated in Mg2+-free conditions. These properties are regulated by backbone dynamics and are key to understanding the functional regulation of CorA.

Identification of a Catalytic Nucleophile-Activating Network in the endo-α-N-Acetylgalactosaminidase of Family GH101.

Bifidobacterium longum endo-α-N-acetylgalactosaminidase (GH101), EngBF, is highly specific toward the mucin Core 1 glycan, Galß1-3GalNAc. Apart from the side chains involved in the retaining mechanism of EngBF, Asp-682 is important for the activity. In the crystal structures of both EngBF and EngSP (from Streptococcus pneumoniae), we identified a conserved water molecule in proximity to Asp-682 and the homologue residue in EngSP. The water molecule also coordinates the catalytic nucleophile and three other residues conserved in GH101 enzymes; in EngBF, these residues are His-685, His-718, and Asn-720. With casein-glycomacropeptide as the substrate, the importance of Asp-682 was confirmed by the lack of a detectable activity for the D682N enzyme. The enzyme variants, H685A, H718A, H685Q, and H718Q, all displayed only a modestly reduction in kcat of up to 15 fold for the H718A variant. However, the double-substituted variants, H685A/H718A and H685Q/H718Q, had a greatly reduced kcat value by about 200 fold compared to that of wild-type EngBF. With the synthetic substrate, Galß(1-3)GalNAcα1-para-nitrophenol, kcat of the double-substituted variants was only up to 30-fold reduced and was found to increase with pH. Compared to the pre-steady-state kinetics of wild-type EngBF, a burst of about the size of the enzyme concentration was absent with the double-substituted EngBF variants, indicating that the nucleophilic attack had become at least as slow as the hydrolysis of the enzyme intermediate. Together, the results indicate that not only Asp-682 but also the entire conserved network of His-685, His-718, and what we suggest is a catalytic water molecule is important in the activation of the catalytic nucleophile.

Cryo-EM structure of human Pol κ bound to DNA and mono-ubiquitylated PCNA.

Y-family DNA polymerase κ (Pol κ) can replicate damaged DNA templates to rescue stalled replication forks. Access of Pol κ to DNA damage sites is facilitated by its interaction with the processivity clamp PCNA and is regulated by PCNA mono-ubiquitylation. Here, we present cryo-EM reconstructions of human Pol κ bound to DNA, an incoming nucleotide, and wild type or mono-ubiquitylated PCNA (Ub-PCNA). In both reconstructions, the internal PIP-box adjacent to the Pol κ Polymerase-Associated Domain (PAD) docks the catalytic core to one PCNA protomer in an angled orientation, bending the DNA exiting the Pol κ active site through PCNA, while Pol κ C-terminal domain containing two Ubiquitin Binding Zinc Fingers (UBZs) is invisible, in agreement with disorder predictions. The ubiquitin moieties are partly flexible and extend radially away from PCNA, with the ubiquitin at the Pol κ-bound protomer appearing more rigid. Activity assays suggest that, when the internal PIP-box interaction is lost, Pol κ is retained on DNA by a secondary interaction between the UBZs and the ubiquitins flexibly conjugated to PCNA. Our data provide a structural basis for the recruitment of a Y-family TLS polymerase to sites of DNA damage.

Accurate model of liquid-liquid phase behavior of intrinsically disordered proteins from optimization of single-chain properties.

Many intrinsically disordered proteins (IDPs) may undergo liquid-liquid phase separation (LLPS) and participate in the formation of membraneless organelles in the cell, thereby contributing to the regulation and compartmentalization of intracellular biochemical reactions. The phase behavior of IDPs is sequence dependent, and its investigation through molecular simulations requires protein models that combine computational efficiency with an accurate description of intramolecular and intermolecular interactions. We developed a general coarse-grained model of IDPs, with residue-level detail, based on an extensive set of experimental data on single-chain properties. Ensemble-averaged experimental observables are predicted from molecular simulations, and a data-driven parameter-learning procedure is used to identify the residue-specific model parameters that minimize the discrepancy between predictions and experiments. The model accurately reproduces the experimentally observed conformational propensities of a set of IDPs. Through two-body as well as large-scale molecular simulations, we show that the optimization of the intramolecular interactions results in improved predictions of protein self-association and LLPS.

Effect of the Water Model in Simulations of Protein-Protein Recognition and Association.

We study self-association of ubiquitin and the disordered protein ACTR using the most commonly used water models. We find that dissociation events are found only with TIP4P-EW and TIP4P/2005, while the widely used TIP3P water model produces straightforward aggregation of the molecules due to the absence of dissociation events. We also find that TIP4P/2005 is the only water model that reproduces the fast association/dissociation dynamics of ubiquitin and best identifies its binding surface. Our results show the critical role of the water model in the description of protein-protein interactions and binding.

DEER-PREdict: Software for efficient calculation of spin-labeling EPR and NMR data from conformational ensembles.

Owing to their plasticity, intrinsically disordered and multidomain proteins require descriptions based on multiple conformations, thus calling for techniques and analysis tools that are capable of dealing with conformational ensembles rather than a single protein structure. Here, we introduce DEER-PREdict, a software program to predict Double Electron-Electron Resonance distance distributions as well as Paramagnetic Relaxation Enhancement rates from ensembles of protein conformations. DEER-PREdict uses an established rotamer library approach to describe the paramagnetic probes which are bound covalently to the protein.DEER-PREdict has been designed to operate efficiently on large conformational ensembles, such as those generated by molecular dynamics simulation, to facilitate the validation or refinement of molecular models as well as the interpretation of experimental data. The performance and accuracy of the software is demonstrated with experimentally characterized protein systems: HIV-1 protease, T4 Lysozyme and Acyl-CoA-binding protein. DEER-PREdict is open source (GPLv3) and available at github.com/KULL-Centre/DEERpredict and as a Python PyPI package pypi.org/project/DEERPREdict.

Two-step reaction mechanism reveals new antioxidant capability of cysteine disulfides against hydroxyl radical attack.

Cysteine disulfides, which constitute an important component in biological redox buffer systems, are highly reactive toward the hydroxyl radical (•OH). The mechanistic details of this reaction, however, remain unclear, largely due to the difficulty in characterizing unstable reaction products. Herein, we have developed a combined approach involving mass spectrometry (MS) and theoretical calculations to investigate reactions of •OH with cysteine disulfides (Cys-S-S-R) in the gas phase. Four types of first-generation products were identified: protonated ions of the cysteine thiyl radical (+Cys-S•), cysteine (+Cys-SH), cysteine sulfinyl radical (+Cys-SO•), and cysteine sulfenic acid (+Cys-SOH). The relative reaction rates and product branching ratios responded sensitively to the electronic property of the R group, providing key evidence to deriving a two-step reaction mechanism. The first step involved •OH conducting a back-side attack on one of the sulfur atoms, forming sulfenic acid (-SOH) and thiyl radical (-S•) product pairs. A subsequent H transfer step within the product complex was favored for protonated systems, generating sulfinyl radical (-SO•) and thiol (-SH) products. Because sulfenic acid is a potent scavenger of peroxyl radicals, our results implied that cysteine disulfide can form two lines of defense against reactive oxygen species, one using the cysteine disulfide itself and the other using the sulfenic acid product of the conversion of cysteine disulfide. This aspect suggested that, in a nonpolar environment, cysteine disulfides might play a more active role in the antioxidant network than previously appreciated.

Computing, Analyzing, and Comparing the Radius of Gyration and Hydrodynamic Radius in Conformational Ensembles of Intrinsically Disordered Proteins.

The level of compaction of an intrinsically disordered protein may affect both its physical and biological properties, and can be probed via different types of biophysical experiments. Small-angle X-ray scattering (SAXS) probe the radius of gyration (Rg) whereas pulsed-field-gradient nuclear magnetic resonance (NMR) diffusion, fluorescence correlation spectroscopy, and dynamic light scattering experiments can be used to determine the hydrodynamic radius (Rh). Here we show how to calculate Rg and Rh from a computationally generated conformational ensemble of an intrinsically disordered protein. We further describe how to use a Bayesian/Maximum Entropy procedure to integrate data from SAXS and NMR diffusion experiments, so as to derive conformational ensembles in agreement with those experiments.

Side chain to main chain hydrogen bonds stabilize a polyglutamine helix in a transcription factor.

Polyglutamine (polyQ) tracts are regions of low sequence complexity frequently found in transcription factors. Tract length often correlates with transcriptional activity and expansion beyond specific thresholds in certain human proteins is the cause of polyQ disorders. To study the structural basis of the association between tract length, transcriptional activity and disease, we addressed how the conformation of the polyQ tract of the androgen receptor, associated with spinobulbar muscular atrophy (SBMA), depends on its length. Here we report that this sequence folds into a helical structure stabilized by unconventional hydrogen bonds between glutamine side chains and main chain carbonyl groups, and that its helicity directly correlates with tract length. These unusual hydrogen bonds are bifurcate with the conventional hydrogen bonds stabilizing α-helices. Our findings suggest a plausible rationale for the association between polyQ tract length and androgen receptor transcriptional activity and have implications for establishing the mechanistic basis of SBMA.

p15PAF binding to PCNA modulates the DNA sliding surface.

p15PAF is an oncogenic intrinsically disordered protein that regulates DNA replication and lesion bypass by interacting with the human sliding clamp PCNA. In the absence of DNA, p15PAF traverses the PCNA ring via an extended PIP-box that contacts the sliding surface. Here, we probed the atomic-scale structure of p15PAF-PCNA-DNA ternary complexes. Crystallography and MD simulations show that, when p15PAF occupies two subunits of the PCNA homotrimer, DNA within the ring channel binds the unoccupied subunit. The structure of PCNA-bound p15PAF in the absence and presence of DNA is invariant, and solution NMR confirms that DNA does not displace p15PAF from the ring wall. Thus, p15PAF reduces the available sliding surfaces of PCNA, and may function as a belt that fastens the DNA to the clamp during synthesis by the replicative polymerase (pol δ). This constraint, however, may need to be released for efficient DNA lesion bypass by the translesion synthesis polymerase (pol η). Accordingly, our biochemical data show that p15PAF impairs primer synthesis by pol η-PCNA holoenzyme against both damaged and normal DNA templates. In light of our findings, we discuss the possible mechanistic roles of p15PAF in DNA replication and suppression of DNA lesion bypass.

Exploring PAZ/3'-overhang interaction to improve siRNA specificity. A combined experimental and modeling study.

The understanding of the dynamical and mechanistic aspects that lie behind siRNA-based gene regulation is a requisite to boost the performance of siRNA therapeutics. A systematic experimental and computational study on the 3'-overhang structural requirements for the design of more specific and potent siRNA molecules was carried out using nucleotide analogues differing in structural parameters, such as sugar constraint, lack of nucleobase, distance between the phosphodiester backbone and nucleobase, enantioselectivity, and steric hindrance. The results established a set of rules governing the siRNA-mediated silencing, indicating that the thermodynamic stability of the 5'-end is a crucial determinant for antisense-mediated silencing but is not sufficient to avoid sense-mediated silencing. Both theoretical and experimental approaches consistently evidence the existence of a direct connection between the PAZ/3'-overhang binding affinity and siRNA's potency and specificity. An overall description of the systems is thus achieved by atomistic simulations and free energy calculations that allow us to propose a robust and self-contained procedure for studying the factors implied in PAZ/3'-overhang siRNA interactions. A higher RNAi activity is associated with a moderate-to-strong PAZ/3'-overhang binding. Contrarily, lower binding energies compromise siRNA potency, increase specificity, and favor siRNA downregulation by Ago2-independent mechanisms. This work provides in-depth details for the design of powerful and safe synthetic nucleotide analogues for substitution at the 3'-overhang, enabling some of the intrinsic siRNA disadvantages to be overcome.

Reactivity of hydropersulfides toward the hydroxyl radical unraveled: disulfide bond cleavage, hydrogen atom transfer, and proton-coupled electron transfer.

Hydropersulfides (RSSH) are highly reactive as nucleophiles and hydrogen atom transfer reagents. These chemical properties are believed to be key for them to act as antioxidants in cells. The reaction involving the radical species and the disulfide bond (S-S) in RSSH, a known redox-active group, however, has been scarcely studied, resulting in an incomplete understanding of the chemical nature of RSSH. We have performed a high-level theoretical investigation on the reactions of the hydroxyl radical (˙OH) toward a set of RSSH (R = -H, -CH3, -NH2, -C(O)OH, -CN, and -NO2). The results show that S-S cleavage and H-atom abstraction are the two competing channels. The electron inductive effect of R induces selective ˙OH substitution at one sulfur atom upon S-S cleavage, forming RSOH and ˙SH for the electron donating groups (EDGs), whereas producing HSOH and ˙SR for the electron withdrawing groups (EWGs). The H-Atom abstraction by ˙OH follows a classical hydrogen atom transfer (hat) mechanism, producing RSS˙ and H2O. Surprisingly, a proton-coupled electron transfer (pcet) process also occurs for R being an EDG. Although for RSSH having EWGs hat is the leading channel, S-S cleavage can be competitive or even dominant for the EDGs. The overall reactivity of RSSH toward ˙OH attack is greatly enhanced with the presence of an EDG, with CH3SSH being the most reactive species found in this study (overall rate constant: 4.55 × 1012 M-1 s-1). Our results highlight the complexity in RSSH reaction chemistry, the extent of which is closely modulated by the inductive effect of the substituents in the case of the oxidation by hydroxyl radicals.

The atmospheric oxidation of CH3OOH by the OH radical: the effect of water vapor.

The atmospheric oxidation of methyl hydroperoxide by the hydroxyl radical has been investigated employing high level theoretical methods. This reaction is important in the chemistry of the troposphere because these species contribute to the oxidizing capacity of the atmosphere and therefore we have studied the bare reaction and the effect of the relative humidity as well. In both cases the reaction can proceed either by abstraction of the terminal hydrogen atom of the OH group, producing CH3O2 + H2O, or by abstraction of one hydrogen atom of the CH3 group, forming H2CO + OH + H2O. We have employed BH&HLYP, QCISD and CCSD(T) theoretical methods along with 6-311+G(2df,2p), aug-cc-pVTZ, aug-cc-pVQZ and CBS basis sets to investigate the reaction mechanism, and conventional and variational transition state theory to study the kinetics of the reaction. For the bare reaction we have computed at room temperature, a rate constant of 3.59 × 10-12 cm3 molecule-1 s-1 for the formation of CH3O2 + H2O and of 1.68 × 10-12 cm3 molecule-1 s-1 for the production of H2CO + OH + H2O, with branching ratios of 68% and 32% respectively. Water vapor enhances the rate constant for the formation of CH3O2 + H2O between 2 and 19%, depending on the temperature and relative humidity, whereas the rate constant for the production of H2CO + OH + H2O is enhanced between 0.3 and 5% by the effect of water vapor under the same conditions, which means that the branching ratio for the formation of CH3O2 + H2O is increased up to 2.5%.

Disentangling polydispersity in the PCNA-p15PAF complex, a disordered, transient and multivalent macromolecular assembly.

The intrinsically disordered p15PAF regulates DNA replication and repair when interacting with the Proliferating Cell Nuclear Antigen (PCNA) sliding clamp. As many interactions between disordered proteins and globular partners involved in signaling and regulation, the complex between p15PAF and trimeric PCNA is of low affinity, forming a transient complex that is difficult to characterize at a structural level due to its inherent polydispersity. We have determined the structure, conformational fluctuations, and relative population of the five species that coexist in solution by combining small-angle X-ray scattering (SAXS) with molecular modelling. By using explicit ensemble descriptions for the individual species, built using integrative approaches and molecular dynamics (MD) simulations, we collectively interpreted multiple SAXS profiles as population-weighted thermodynamic mixtures. The analysis demonstrates that the N-terminus of p15PAF penetrates the PCNA ring and emerges on the back face. This observation substantiates the role of p15PAF as a drag regulating PCNA processivity during DNA repair. Our study reveals the power of ensemble-based approaches to decode structural, dynamic, and thermodynamic information from SAXS data. This strategy paves the way for deciphering the structural bases of flexible, transient and multivalent macromolecular assemblies involved in pivotal biological processes.