Ligand-Induced Variations in Structural and Dynamical Properties Within an Enzyme Superfamily

Ancestral sequence reconstruction dissects structural and functional differences among eosinophil ribonucleases

Evolutionarily conserved structural folds can give rise to diverse biological functions, yet predicting atomic-scale interactions that contribute to the emergence of novel activities within such folds remains challenging. Pancreatic-type ribonucleases illustrate this complexity, sharing a core structure that has evolved to accommodate varied functions. In this study, we used ancestral sequence reconstruction to probe evolutionary and molecular determinants that distinguish biological activities within eosinophil members of the RNase 2/3 subfamily. Our investigation unveils functional, structural, and dynamical behaviors that differentiate the evolved ancestral ribonuclease (AncRNase) from its contemporary eosinophil RNase orthologs. Leveraging the potential of ancestral reconstruction for protein engineering, we used AncRNase predictions to design a minimal 4-residue variant that transforms human RNase 2 into a chimeric enzyme endowed with the antimicrobial and cytotoxic activities of RNase 3 members. This work provides unique insights into mutational and evolutionary pathways governing structure, function, and conformational states within the eosinophil RNase subfamily, offering potential for targeted modulation of RNase-associated functions.

Structural determinants for tRNA selective cleavage by RNase 2/EDN

tRNA-derived fragments (tRFs) have emerged as key players of immunoregulation. Some RNase A superfamily members participate in the shaping of the tRFs population. By comparing wild-type and knockout macrophage cell lines, our previous work revealed that RNase 2 can selectively cleave tRNAs. Here, we confirm the in vitro protein cleavage pattern by screening of synthetic tRNAs, single-mutant variants, and anticodon-loop DNA/RNA hairpins. By sequencing of tRF products, we identified the cleavage selectivity of recombinant RNase 2 with base specificity at B1 (U/C) and B2 (A) sites, consistent with a previous cellular study. Lastly, protein-hairpin complexes were predicted by MD simulations. Results reveal the contribution of the α1, loop 3 and loop 4, and β6 RNase 2 regions, where residues Arg36/Asn39/Gln40/Asn65/Arg68/Arg132 provide interactions, spanning from P-1 to P2 sites that are essential for anticodon loop recognition. Knowledge of RNase 2-specific tRFs generation might guide new therapeutic approaches for infectious and immune-related diseases.

Biochemical Discovery, Intracellular Evaluation, and Crystallographic Characterization of Synthetic and Natural Product Adenosine 3',5'-Cyclic Monophosphate-Dependent Protein Kinase A (PKA) Inhibitors

The recent demonstration that adenosine 3',5'-cyclic monophosphate (cAMP)-dependent protein kinase A (PKA) plays an oncogenic role in a number of important cancers has led to a renaissance in drug development interest targeting this kinase. We therefore have established a suite of biochemical, cell-based, and structural biology assays for identifying and evaluating new pharmacophores for PKA inhibition. This discovery process started with a 384-well high-throughput screen of more than 200,000 substances, including fractionated natural product extracts. Identified active compounds were further prioritized in biochemical, biophysical, and cell-based assays. Priority lead compounds were assessed in detail to fully characterize several previously unrecognized PKA pharmacophores including the generation of new X-ray crystallography structures demonstrating unique interactions between PKA and bound inhibitor molecules.

Conformational exchange divergence along the evolutionary pathway of eosinophil-associated ribonucleases

The evolutionary role of conformational exchange in the emergence and preservation of function within structural homologs remains elusive. While protein engineering has revealed the importance of flexibility in function, productive modulation of atomic-scale dynamics has only been achieved on a finite number of distinct folds. Allosteric control of unique members within dynamically diverse structural families requires a better appreciation of exchange phenomena. Here, we examined the functional and structural role of conformational exchange within eosinophil-associated ribonucleases. Biological and catalytic activity of various EARs was performed in parallel to mapping their conformational behavior on multiple timescales using NMR and computational analyses. Despite functional conservation and conformational seclusion to a specific domain, we show that EARs can display similar or distinct motional profiles, implying divergence rather than conservation of flexibility. Comparing progressively more distant enzymes should unravel how this subfamily has evolved new functions and/or altered their behavior at the molecular level.

pH-Triggered Recovery of Organic Polymer Photocatalytic Particles for the Production of High Value Compounds and Enhanced Recyclability

Pseudo-homogeneous polymeric photocatalysts are an emerging class of highly efficient and tunable photocatalytic materials, where the photocatalytic centers are easily accessible. The creation of highly efficient photocatalytic materials that can be rapidly separated and recovered is one of the critical challenges in photocatalytic chemistry. Here, we describe pH-responsive photocatalytic nanoparticles that are active and well-dispersed under acidic conditions but aggregate instantly upon elevation of pH, enabling easy recovery. These responsive photocatalytic polymers can be used in various photocatalytic transformations, including Cr^VI reduction and photoredox alkylation of indole derivative. Notably, the cationic nature of the photocatalyst accelerates reaction rate of an anionic substrate compared to uncharged species. These photocatalytic particles could be readily recycled allowing multiple successive photocatalytic reactions with no clear loss in activity.

Analysis of coordinated NMR chemical shifts to map allosteric regulatory networks in proteins

The exquisite sensitivity of the NMR chemical shift to local environment makes it an ideal probe to assess atomic level perturbations in proteins of all sizes and structural compositions. Recent advances in solution and solid-state NMR spectroscopy of biomolecules have leveraged the chemical shift to report on short- and long-range couplings between individual amino acids to establish "networks" of residues that form the basis of allosteric pathways that transmit chemical signals through the protein matrix to induce functional responses. The simple premise that thermodynamically and functionally coupled regions of a protein (i.e. active and allosteric sites) should be reciprocally sensitive to structural or dynamic perturbations has enabled NMR spectroscopy, the premier method for molecular resolution of protein structural fluctuations, to occupy a place at the forefront of investigations into protein allostery. Here, we detail several key methods of NMR chemical shift analysis to extract mechanistic information about long-range chemical signaling in a protein, focusing on practical methodological aspects and the circumstances under which a given approach would be relevant. We also detail some of the experimental considerations that should be made when applying these methods to specific protein systems.

Divergent allostery reveals critical differences between structurally homologous regulatory domains of Plasmodium falciparum and human protein kinase G

Malaria is a life-threatening infectious disease primarily caused by the Plasmodium falciparum parasite. The increasing resistance to current antimalarial drugs and their side effects has led to an urgent need for novel malaria drug targets, such as the P. falciparum cGMP-dependent protein kinase (pfPKG). However, PKG plays an essential regulatory role also in the human host. Human PKG (hPKG) and pfPKG are controlled by structurally homologs cGMP-binding domains (CBDs). Here, we show that despite the structural similarities between the essential CBDs in pfPKG and hPKG, their respective allosteric networks differ significantly. Through comparative analyses of CHESCA, molecular dynamics simulations, and backbone internal dynamics measurements, we found that conserved allosteric elements within the essential CBDs are wired differently in pfPKG and hPKG to implement cGMP-dependent kinase activation. Such pfPKG vs. hPKG rewiring of allosteric networks was unexpected due to the structural similarity between the two essential CBDs. Yet, such finding provides crucial information on which elements to target for selective inhibition of pfPKG vs. hPKG, which may potentially reduce undesired side-effects in malaria treatments.

In Silico Engineering of Enzyme Access Tunnels

Enzyme engineering is a tailoring process that allows the modification of naturally occurring enzymes to provide them with improved catalytic efficiency, stability, or specificity. By introducing partial modifications to their sequence and to their structural features, enzyme engineering can transform natural enzymes into more efficient, specific and resistant biocatalysts and render them suitable for virtually countless industrial processes. Current enzyme engineering methods mostly target the active site of the enzyme, where the catalytic reaction takes place. Nonetheless, the tunnel that often connects the surface of an enzyme with its buried active site plays a key role in the activity of the enzyme as it acts as a gatekeeper and regulates the access of the substrate to the catalytic pocket. Hence, there is an increasing interest in targeting the sequence and the structure of substrate entrance tunnels in order to fine-tune enzymatic activity, regulate substrate specificity, or control reaction promiscuity.In this chapter, we describe the use of a rational in silico design and screening method to engineer the access tunnel of a fructosyl peptide oxidase with the aim to facilitate access to its catalytic site and to expand its substrate range. Our goal is to engineer this class of enzymes in order to utilize them for the direct detection of glycated proteins in diabetes monitoring devices. The design strategy involves remodeling of the backbone structure of the enzyme , a feature that is not possible with conventional enzyme engineering techniques such as single-point mutagenesis and that is highly unlikely to occur using a directed evolution approach.The proposed strategy, which results in a significant reduction in cost and time for the experimental production and characterization of candidate enzyme variants, represents a promising approach to the expedited identification of novel and improved enzymes. Rational enzyme design aims to provide in silico strategies for the fast, accurate, and inexpensive development of biocatalysts that can meet the needs of multiple industrial sectors, thus ultimately promoting the use of green chemistry and improving the efficiency of chemical processes.

Convenient and accurate insight into solution-phase equilibria from FlowNMR titrations

Solution phase titrations are made easy by multi-nuclear FlowNMR spectroscopy with automated, continuous titre addition to give accurate insights into Brønsted acid/base, hydrogen bonding, Lewis acid/base and metal/ligand binding equilibria under native conditions.

State-selective frustration as a key driver of allosteric pluripotency

Allosteric pluripotency arises when an allosteric effector switches from agonist to antagonist depending on the experimental conditions. For example, the Rp-cAMPS ligand of Protein Kinase A (PKA) switches from agonist to antagonist as the MgATP concentration increases and/or the kinase substrate affinity or concentration decreases. Understanding allosteric pluripotency is essential to design effective allosteric therapeutics with minimal side effects. Allosteric pluripotency of PKA arises from divergent allosteric responses of two homologous tandem cAMP-binding domains, resulting in a free energy landscape for the Rp-cAMPS-bound PKA regulatory subunit R1a in which the ground state is kinase inhibition-incompetent and the kinase inhibition-competent state is excited. The magnitude of the free energy difference between the ground non-inhibitory and excited inhibitory states (ΔG_R,Gap) relative to the effective free energy of R1a binding to the catalytic subunit of PKA (ΔG_R:C) dictates whether the antagonism-to-agonism switch occurs. However, the key drivers of ΔG_R,Gap are not fully understood. Here, by analyzing an R1a mutant that selectively silences allosteric pluripotency, we show that a major determinant of ΔG_R,Gap unexpectedly arises from state-selective frustration in the ground inhibition-incompetent state of Rp-cAMPS-bound R1a. Such frustration is caused by steric clashes between the phosphate-binding cassette and the helices preceding the lid, which interact with the phosphate and base of Rp-cAMPS, respectively. These clashes are absent in the excited inhibitory state, thus reducing the ΔG_R,Gap to values comparable to ΔG_R:C, as needed for allosteric pluripotency to occur. The resulting model of allosteric pluripotency is anticipated to assist the design of effective allosteric modulators.

The Rp-cAMPS ligand of protein kinase A switches from agonist to antagonist depending on metabolite and proteomic contexts. We show that the state-selective frustration is a key driver of this allosteric pluripotency phenomenon.

Insights into Structural and Dynamical Changes Experienced by Human RNase 6 upon Ligand Binding

Ribonuclease 6 (RNase 6) is one of eight catalytically active human pancreatic-type RNases that belong to a superfamily of rapidly evolving enzymes. Like some of its human homologues, RNase 6 exhibits host defense properties such as antiviral and antibacterial activities. Recently solved crystal structures of this enzyme in its nucleotide-free form show the conservation of the prototypical kidney-shaped fold preserved among vertebrate RNases, in addition to revealing the presence of a unique secondary active site. In this study, we determine the structural and conformational properties experienced by RNase 6 upon binding to substrate and product analogues. We present the first crystal structures of RNase 6 bound to a nucleotide ligand (adenosine 5'-monophosphate), in addition to RNase 6 bound to phosphate ions. While the enzyme preserves B2 subsite ligand preferences, our results show a lack of typical B2 subsite interactions normally observed in homologous ligand-bound RNases. A comparison of the dynamical properties of RNase 6 in its apo-, substrate-, and product-bound states highlight the unique dynamical properties experienced on time scales ranging from nano- to milliseconds. Overall, our results confirm the specific evolutionary adaptation of RNase 6 relative to its unique catalytic and biological activities.

Evolutionary Trends in RNA Base Selectivity Within the RNase A Superfamily

There is a growing interest in the pharmaceutical industry to design novel tailored drugs for RNA targeting. The vertebrate-specific RNase A superfamily is nowadays one of the best characterized family of enzymes and comprises proteins involved in host defense with specific cytotoxic and immune-modulatory properties. We observe within the family a structural variability at the substrate-binding site associated to a diversification of biological properties. In this work, we have analyzed the enzyme specificity at the secondary base binding site. Towards this end, we have performed a kinetic characterization of the canonical RNase types together with a molecular dynamic simulation of selected representative family members. The RNases' catalytic activity and binding interactions have been compared using UpA, UpG and UpI dinucleotides. Our results highlight an evolutionary trend from lower to higher order vertebrates towards an enhanced discrimination power of selectivity for adenine respect to guanine at the secondary base binding site (B2). Interestingly, the shift from guanine to adenine preference is achieved in all the studied family members by equivalent residues through distinct interaction modes. We can identify specific polar and charged side chains that selectively interact with donor or acceptor purine groups. Overall, we observe selective bidentate polar and electrostatic interactions: Asn to N1/N6 and N6/N7 adenine groups in mammals versus Glu/Asp and Arg to N1/N2, N1/O6 and O6/N7 guanine groups in non-mammals. In addition, kinetic and molecular dynamics comparative results on UpG versus UpI emphasize the main contribution of Glu/Asp interactions to N1/N2 group for guanine selectivity in lower order vertebrates. A close inspection at the B2 binding pocket also highlights the principal contribution of the protein ß6 and L4 loop regions. Significant differences in the orientation and extension of the L4 loop could explain how the same residues can participate in alternative binding modes. The analysis suggests that within the RNase A superfamily an evolution pressure has taken place at the B2 secondary binding site to provide novel substrate-recognition patterns. We are confident that a better knowledge of the enzymes' nucleotide recognition pattern would contribute to identify their physiological substrate and eventually design applied therapies to modulate their biological functions.

Nucleotide substrate binding characterization in human pancreatic-type ribonucleases

Human genome contains a group of more than a dozen similar genes with diverse biological functions including antiviral, antibacterial and angiogenesis activities. The characterized gene products of this group show significant sequence similarity and a common structural fold associated with binding and cleavage of ribonucleic acid (RNA) substrates. Therefore, these proteins have been categorized as members of human pancreatic-type ribonucleases (hRNases). hRNases differ in cell/tissue localization and display distinct substrate binding preferences and a wide range of ribonucleolytic catalytic efficiencies. Limited information is available about structural and dynamical properties that influence this diversity among these homologous RNases. Here, we use computer simulations to characterize substrate interactions, electrostatics and dynamical properties of hRNases 1-7 associated with binding to two nucleotide substrates (ACAC and AUAU). Results indicate that even with complete conservation of active-site catalytic triad associated with ribonucleolytic activity, these enzymes show significant differences in substrate interactions. Detailed characterization suggests that in addition to binding site electrostatic and van der Waals interactions, dynamics of distal regions may also play a role in binding. Another key insight is that a small difference in temperature of 300 K (used in experimental studies) and 310 K (physiological temperature) shows significant changes in enzyme-substrate interactions.

Engineering enzyme access tunnels

Enzymes are efficient and specific catalysts for many essential reactions in biotechnological and pharmaceutical industries. Many times, the natural enzymes do not display the catalytic efficiency, stability or specificity required for these industrial processes. The current enzyme engineering methods offer solutions to this problem, but they mainly target the buried active site where the chemical reaction takes place. Despite being many times ignored, the tunnels and channels connecting the environment with the active site are equally important for the catalytic properties of enzymes. Changes in the enzymatic tunnels and channels affect enzyme activity, specificity, promiscuity, enantioselectivity and stability. This review provides an overview of the emerging field of enzyme access tunnel engineering with case studies describing design of all the aforementioned properties. The software tools for the analysis of geometry and function of the enzymatic tunnels and channels and for the rational design of tunnel modifications will also be discussed. The combination of new software tools and enzyme engineering strategies will provide enzymes with access tunnels and channels specifically tailored for individual industrial processes.