Pseudouridylation defect due to DKC1 and NOP10 mutations causes nephrotic syndrome with cataracts, hearing impairment, and enterocolitis.

RNA modifications play a fundamental role in cellular function. Pseudouridylation, the most abundant RNA modification, is catalyzed by the H/ACA small ribonucleoprotein (snoRNP) complex that shares four core proteins, dyskerin (DKC1), NOP10, NHP2, and GAR1. Mutations in DKC1, NOP10, or NHP2 cause dyskeratosis congenita (DC), a disorder characterized by telomere attrition. Here, we report a phenotype comprising nephrotic syndrome, cataracts, sensorineural deafness, enterocolitis, and early lethality in two pedigrees: males with DKC1 p.Glu206Lys and two children with homozygous NOP10 p.Thr16Met. Females with heterozygous DKC1 p.Glu206Lys developed cataracts and sensorineural deafness, but nephrotic syndrome in only one case of skewed X-inactivation. We found telomere attrition in both pedigrees, but no mucocutaneous abnormalities suggestive of DC. Both mutations fall at the dyskerin-NOP10 binding interface in a region distinct from those implicated in DC, impair the dyskerin-NOP10 interaction, and disrupt the catalytic pseudouridylation site. Accordingly, we found reduced pseudouridine levels in the ribosomal RNA (rRNA) of the patients. Zebrafish dkc1 mutants recapitulate the human phenotype and show reduced 18S pseudouridylation, ribosomal dysregulation, and a cell-cycle defect in the absence of telomere attrition. We therefore propose that this human disorder is the consequence of defective snoRNP pseudouridylation and ribosomal dysfunction.

Off-pathway 3D-structure provides protection against spontaneous Asn/Asp isomerization: shielding proteins Achilles heel.

Spontaneous deamidation prompted backbone isomerization of Asn/Asp residues resulting in - most cases - the insertion of an extra methylene group into the backbone poses a threat to the structural integrity of proteins. Here we present a systematical analysis of how temperature, pH, presence of charged residues, but most importantly backbone conformation and dynamics affect isomerization rates as determined by nuclear magnetic resonance in the case of designed peptide-models. We demonstrate that restricted mobility (such as being part of a secondary structural element) may safeguard against isomerization, but this protective factor is most effective in the case of off-pathway folds which can slow the reaction by several magnitudes compared to their on-pathway counterparts. We show that the geometric descriptors of the initial nucleophilic attack of the isomerization can be used to classify local conformation and contribute to the design of stable protein drugs, antibodies or the assessment of the severity of mutations.At any ­Asn/AspGly­ sites in proteins a spontaneous backbone isomerization occurs within days under physiological conditions leading to various forms of proteopathy. This unwanted transformation especially harmful to long-lived proteins (e.g. hemoglobin and crystallins), can be slowed down, though never stopped, by a rigid three-dimensional protein fold, if it can delay in the conformational maze, on-pathway intermediates from occurring.

The Importance of Mg2+ -Free State in Nucleotide Exchange of Oncogenic K-Ras Mutants.

For efficient targeting of oncogenic K-Ras interaction sites, a mechanistic picture of the Ras-cycle is necessary. Herein, we used NMR relaxation techniques and molecular dynamics simulations to decipher the role of slow dynamics in wild-type and three oncogenic P-loop mutants of K-Ras. Our measurements reveal a dominant two-state conformational exchange on the ms timescale in both GDP- and GTP-bound K-Ras. The identified low-populated higher energy state in GDP-loaded K-Ras has a conformation reminiscent of a nucleotide-bound/Mg2+ -free state characterized by shortened ß2/ß3-strands and a partially released switch-I region preparing K-Ras for the interaction with the incoming nucleotide exchange factor and subsequent reactivation. By providing insight into mutation-specific differences in K-Ras structural dynamics, our systematic analysis improves our understanding of prolonged K-Ras signaling and may aid the development of allosteric inhibitors targeting nucleotide exchange in K-Ras.

Omicron Binding Mode: Contact Analysis and Dynamics of the Omicron Receptor-Binding Domain in Complex with ACE2.

On 26 November 2021, the WHO classified the Omicron variant of the SARS-CoV-2 virus (B.1.1.529 lineage) as a variant of concern (VOC) (COVID-19 Variant Data, Department of Health, 2022). The Omicron variant contains as many as 26 unique mutations of effects not yet determined (Venkatakrishnan, A., Open Science Framework, 2021). Out of its total of 34 Spike protein mutations, 15 are located on the receptor-binding domain (S-RBD) (Stanford Coronavirus Antiviral & Resistance Database, 2022) that directly contacts the angiotensin-converting enzyme 2 (ACE2) host receptor and is also a primary target for antibodies. Here, we studied the binding mode of the S-RBD domain of the Spike protein carrying the Omicron mutations and the globular domain of human ACE2 using molecular dynamics (MD) simulations. We identified new and key Omicron-specific interactions such as R493 (of mutation Q493R), which forms salt bridges both with E35 and D38 of ACE2, Y501 (N501Y), which forms an edge-to-face aromatic interaction with Y41, and Y505 (Y505H), which makes an H-bond with E37 and K353. The glycan chains of ACE2 also bind differently in the WT and Omicron variants in response to different charge distributions on the surface of Spike proteins. However, while the Omicron mutations considerably improve the overall electrostatic fit of the two interfaces, the total number of specific and favorable interactions between the two does not increase. The dynamics of the complexes are highly affected too, making the Omicron S-RBD:ACE2 complex more rigid; the two main interaction sites, Patches I and II, isolated in the WT complex, become connected in the Omicron complex through the alternating interaction of R493 and R498 with E35 and D38.

Dynamically encoded reactivity of Ras enzymes: opening new frontiers for drug discovery.

Decoding molecular flexibility in order to understand and predict biological processes-applying the principles of dynamic-structure-activity relationships (DSAR)-becomes a necessity when attempting to design selective and specific inhibitors of a protein that has overlapping interaction surfaces with its upstream and downstream partners along its signaling cascade. Ras proteins are molecular switches that meet this definition perfectly. The close-lying P-loop and the highly flexible switch I and switch II regions are the site of nucleotide-, assisting-, and effector-protein binding. Oncogenic mutations that also appear in this region do not cause easily characterized overall structural changes, due partly to the inherent conformational heterogeneity and pliability of these segments. In this review, we present an overview of the results obtained using approaches targeting Ras dynamics, such as nuclear magnetic resonance (NMR) measurements and experiment-based modeling calculations (mostly molecular dynamics (MD) simulations). These methodologies were successfully used to decipher the mutant- and isoform-specific nature of certain transient states, far-lying allosteric sites, and the internal interaction networks, as well as the interconnectivity of the catalytic and membrane-binding regions. This opens new therapeutic potential: the discovered interaction hotspots present hitherto not targeted, selective sites for drug design efforts in diverse locations of the protein matrix.

Tumor-Suppressor p53TAD1-60 Forms a Fuzzy Complex with Metastasis-Associated S100A4: Structural Insights and Dynamics by an NMR/MD Approach.

Conformationally flexible protein complexes represent a major challenge for structural and dynamical studies. We present herein a method based on a hybrid NMR/MD approach to characterize the complex formed between the disordered p53TAD1-60 and the metastasis-associated S100A4. Disorder-to-order transitions of both TAD1 and TAD2 subdomains upon interaction is detected. Still, p53TAD1-60 remains highly flexible in the bound form, with residues L26, M40, and W53 being anchored to identical hydrophobic pockets of the S100A4 monomer chains. In the resulting "fuzzy" complex, the clamp-like binding of p53TAD1-60 relies on specific hydrophobic anchors and on the existence of extended flexible segments. Our results demonstrate that structural and dynamical NMR parameters (cumulative Δδ, SSP, temperature coefficients, relaxation time, hetNOE) combined with MD simulations can be used to build a structural model even if, due to high flexibility, the classical solution structure calculation is not possible.

Configuration-Controlled Crystal and/or Gel Formation of Protected d-Glucosamines Supported by Promiscuous Interaction Surfaces and a Conformationally Heterogeneous Solution State.

The configuration-dependent self-association mode of the two anomers of O-Ac,N-Fmoc-d-glucosamine, a foldamer building block, leading to gel and/or single crystal formation is described. The ß-anomer of the sugar amino acid (2) forms a gel from various solvents (confirmed by SEM, rheology measurements, NMR, and ECD spectroscopy), whereas the α-anomer (1) does not form a gel with any solvent tested. Transition from the solution state to a gel is coupled to a concurrent shift of the Fmoc-groups: from a freely rotating (almost symmetrical) to a specific, asymmetric orientation. Whereas the crystal structure of the α-anomer is built as an evenly packed 3D system, the ß-anomer forms a looser superstructure of well-packed 2D layers. Modeling indicates that in the lowest energy, but scarcely sampled conformer of the ß-anomer, the Fmoc-group bends above the sugar moiety, stabilized by intramolecular CH↔π interactions between the aromatic rings. It is concluded that possessing an extended and promiscuous interaction surface and a conformationally heterogeneous solution state are among the basic requirements of gel formation for a candidate molecule.

The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein.

Invited for the cover of this issue is the group of András Perczel at Eötvös Loránd University, Budapest, Hungary and colleagues from Osaka University, Japan. The image depicts the amyloid buildup of an Exenatide derivate miniprotein (E5) monitored on a simplified hyperspace. Read the full text of the article at 10.1002/chem.201903826.

The Route from the Folded to the Amyloid State: Exploring the Potential Energy Surface of a Drug-Like Miniprotein.

The amyloid formation of the folded segment of a variant of Exenatide (a marketed drug for type-2 diabetes mellitus) was studied by electronic circular dichroism (ECD) and NMR spectroscopy. We found that the optimum temperature for E5 protein amyloidosis coincides with body temperature and requires well below physiological salt concentration. Decomposition of the ECD spectra and its barycentric representation on the folded-unfolded-amyloid potential energy surface allowed us to monitor the full range of molecular transformation of amyloidogenesis. We identified points of no return (e.g.; T=37 °C, pH 4.1, cE5 =250 µm, cNaCl =50 mm, t>4-6 h) that will inevitably gravitate into the amyloid state. The strong B-type far ultraviolet (FUV)-ECD spectra and an unexpectedly strong near ultraviolet (NUV)-ECD signal (Θ≈275-285   nm ) indicate that the amyloid phase of E5 is built from monomers of quasi-elongated backbone structure (φ≈-145°, ψ≈+145°) with strong interstrand Tyr↔Trp interaction. Misfolded intermediates and the buildup of "toxic" early-stage oligomers leading to self-association were identified and monitored as a function of time. Results indicate that the amyloid transition is triggered by subtle misfolding of the α-helix, exposing aromatic and hydrophobic side chains that may provide the first centers for an intermolecular reorganization. These initial clusters provide the spatial closeness and sufficient time for a transition to the ß-structured amyloid nucleus, thus the process follows a nucleated growth mechanism.

Gas Sensing by Bacterial H-NOX Proteins: An MD Study.

Gas sensing is crucial for both prokaryotes and eukaryotes and is primarily performed by heme-based sensors, including H-NOX domains. These systems may provide a new, alternative mode for transporting gaseous molecules in higher organisms, but for the development of such systems, a detailed understanding of the ligand-binding properties is required. Here, we focused on ligand migration within the protein matrix: we performed molecular dynamics simulations on three bacterial (Ka, Ns and Cs) H-NOX proteins and studied the kinetics of CO, NO and O2 diffusion. We compared the response of the protein structure to the presence of ligands, diffusion rate constants, tunnel systems and storage pockets. We found that the rate constant for diffusion decreases in the O2 > NO > CO order in all proteins, and in the Ns > Ks > Cs order if single-gas is considered. Competition between gases seems to seriously influence the residential time of ligands spent in the distal pocket. The channel system is profoundly determined by the overall fold, but the sidechain pattern has a significant role in blocking certain channels by hydrophobic interactions between bulky groups, cation-π interactions or hydrogen bonding triads. The majority of storage pockets are determined by local sidechain composition, although certain functional cavities, such as the distal and proximal pockets are found in all systems. A major guideline for the design of gas transport systems is the need to chemically bind the gas molecule to the protein, possibly joining several proteins with several heme groups together.

Assignment of Vibrational Circular Dichroism Cross-Referenced Electronic Circular Dichroism Spectra of Flexible Foldamer Building Blocks: Towards Assigning Pure Chiroptical Properties of Foldamers.

Assignment of the most established electronic circular dichroism (ECD) spectra of polypeptides and foldamers is either "evidence based" or relies on the 3D structures of longer oligomers of limited internal dynamics, which are derived from NMR spectroscopy (or X-ray) data. Critics warn that the use of NMR spectroscopy and ECD side by side has severe limitations for flexible molecules because explicit knowledge of conformational ensembles is a challenge. Herein, an old-new method of comparing ab initio computed and measured vibrational circular dichroism (VCD) data is presented to validate both the structures (conf(i)) and their relative weights (c(i)) that make up the conformational ensemble. Based on the array of {conf(i), c(i)}, the pure ECD spectra, g(i)conf(i) , can be ab initio calculated. The reconstructed spectrum Σc(i)g(i)conf(i) can thus help to assign any experimental ECD counterparts. Herein, such a protocol is successfully applied to flexible foldamer building blocks of sugar ß-amino acid diamides. The epimeric pair of the model system was selected because these molecules were conformationally tunable by simple chemical modification, and thus, the robustness of the current approach could be probed. The initial hydrogen bond (NH⋅⋅⋅O) eliminated by N-methylation reorients the amide plain, which influences the chiroptical properties of the foldamer building block; this structural change is successfully monitored by changes to the VCD and ECD transitions, which are now assigned to pure conformers. The current method seems to be general and effective without requiring extensive CPU and spectroscopic resources.

First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO.

Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O2 , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.

Predictable Conformational Diversity in Foldamers of Sugar Amino Acids.

A systematic conformational search was carried out for monomers and homohexamers of furanoid ß-amino acids: cis-(S,R) and trans-(S,S) stereoisomers of aminocyclopentane carboxylic acid (ACPC), two different aminofuranuronic acids (AFUα and AFUß), their isopropylidene derivatives (AFU(ip)), and the key intermediate ß-aminotetrahydrofurancarboxylic acid (ATFC). The stereochemistry of the building blocks was chosen to match that of the natural sugar amino acid (xylose and ribose) precursors (XylAFU and RibAFU). The results show that hexamers of cis-furanoid ß-amino acids show great variability: while hydrophobic cyclopentane (cis-ACPC)6 and hydrophilic (XylAFUα/ß)6 foldamers favor two different zigzagged conformation as hexamers, the backbone fold turns into a helix in the case of (cis-ATFC)6 (10-helix) and (XylAFU(ip))6 (14-helix). Trans stereochemistry resulted in hexamers exclusively with the right-handed helix conformation, (H12P)6, regardless of their polarity. We found that the preferred oligomeric structure of XylAFUα/ß is conformationally compatible with ß-pleated sheets, while that of the trans/(S,S) units matches with α-helices of proteins.

A self-compartmentalizing hexamer serine protease from Pyrococcus horikoshii: substrate selection achieved through multimerization.

Oligopeptidases impose a size limitation on their substrates, the mechanism of which has long been under debate. Here we present the structure of a hexameric serine protease, an oligopeptidase from Pyrococcus horikoshii (PhAAP), revealing a complex, self-compartmentalized inner space, where substrates may access the monomer active sites passing through a double-gated "check-in" system, first passing through a pore on the hexamer surface and then turning to enter through an even smaller opening at the monomers' domain interface. This substrate screening strategy is unique within the family. We found that among oligopeptidases, a residue of the catalytic apparatus is positioned near an amylogenic ß-edge, which needs to be protected to prevent aggregation, and we found that different oligopeptidases use different strategies to achieve such an end. We propose that self-assembly within the family results in characteristically different substrate selection mechanisms coupled to different multimerization states.

Structure and catalysis of acylaminoacyl peptidase: closed and open subunits of a dimer oligopeptidase.

Acylaminoacyl peptidase from Aeropyrum pernix is a homodimer that belongs to the prolyl oligopeptidase family. The monomer subunit is composed of one hydrolase and one propeller domain. Previous crystal structure determinations revealed that the propeller domain obstructed the access of substrate to the active site of both subunits. Here we investigated the structure and the kinetics of two mutant enzymes in which the aspartic acid of the catalytic triad was changed to alanine or asparagine. Using different substrates, we have determined the pH dependence of specificity rate constants, the rate-limiting step of catalysis, and the binding of substrates and inhibitors. The catalysis considerably depended both on the kind of mutation and on the nature of the substrate. The results were interpreted in terms of alterations in the position of the catalytic histidine side chain as demonstrated with crystal structure determination of the native and two mutant structures (D524N and D524A). Unexpectedly, in the homodimeric structures, only one subunit displayed the closed form of the enzyme. The other subunit exhibited an open gate to the catalytic site, thus revealing the structural basis that controls the oligopeptidase activity. The open form of the native enzyme displayed the catalytic triad in a distorted, inactive state. The mutations affected the closed, active form of the enzyme, disrupting its catalytic triad. We concluded that the two forms are at equilibrium and the substrates bind by the conformational selection mechanism.

Chemotaxis induced by SXWS tetrapeptides in Tetrahymena--overlapping chemotactic effects of SXWS sequences and their identical amino acids.

The chemotactic potential of SXWS peptides and the components of the extracellular domain of cytokine receptors were investigated in Tetrahymena as a functional index of substitution with different amino acids in the position 'X' of the tetrapeptide. Data obtained demonstrate that position X plays a special determining role in the ligand, SEWS and STWS possess extremely strong chemoattractant ability, and aromatic amino acids result in chemorepellent ligands. Diverse effects of structurally related molecules, for example, SNWS-SDWS, demonstrate a highly sensitive discrimination potential in the applied model system. Physicochemical characteristics (hydropathy, residue size, and solvent-exposed area) of the amino acids were correlated with the chemotactic activity. Data obtained by computer-assisted conformation analysis of SXWS peptides and the highly overlapping chemotactic effects of the investigated SXWS peptides as well as the presence of the amino acids in the 'X' position indicate that member 'X' of the SXWS sequence performs a special role in interactions with the chemotaxis receptors in the membrane.

Cryo-EM structure of acylpeptide hydrolase reveals substrate selection by multimerization and a multi-state serine-protease triad.

The first structure of tetrameric mammalian acylaminoacyl peptidase, an enzyme that functions as an upstream regulator of the proteasome through the removal of terminal N-acetylated residues from its protein substrates, was determined by cryo-EM and further elucidated by MD simulations. Self-association results in a toroid-shaped quaternary structure, guided by an amyloidogenic ß-edge and unique inserts. With a Pro introduced into its central ß-sheet, sufficient conformational freedom is awarded to the segment containing the catalytic Ser587 that the serine protease catalytic triad alternates between active and latent states. Active site flexibility suggests that the dual function of catalysis and substrate selection are fulfilled by a novel mechanism: substrate entrance is regulated by flexible loops creating a double-gated channel system, while binding of the substrate to the active site is required for stabilization of the catalytic apparatus - as a second filter before hydrolysis. The structure not only underlines that within the family of S9 proteases homo-multimerization acts as a crucial tool for substrate selection, but it will also allow drug design targeting of the ubiquitin-proteasome system.

A carbapenem antibiotic inhibiting a mammalian serine protease: structure of the acylaminoacyl peptidase-meropenem complex.

The structure of porcine AAP (pAAP) in a covalently bound complex with meropenem was determined by cryo-EM to 2.1 Å resolution, showing the mammalian serine-protease inhibited by a carbapenem antibiotic. AAP is a modulator of the ubiquitin-proteasome degradation system and the site of a drug-drug interaction between the widely used antipsychotic, valproate and carbapenems. The active form of pAAP - a toroidal tetramer - binds four meropenem molecules covalently linked to the catalytic Ser587 of the serine-protease triad, in an acyl-enzyme state. AAP is hindered from fully processing the antibiotic by the displacement and protonation of His707 of the catalytic triad. We show that AAP is made susceptible to the association by its unusually sheltered active pockets and flexible catalytic triads, while the carbapenems possess sufficiently small substituents on their ß-lactam rings to fit into the shallow substrate-specificity pocket of the enzyme.

Directed Evolution-Driven Increase of Structural Plasticity Is a Prerequisite for Binding the Complement Lectin Pathway Blocking MASP-Inhibitor Peptides.

MASP-1 and MASP-2 are key activator proteases of the complement lectin pathway. The first specific mannose-binding lectin-associated serine protease (MASP) inhibitors had been developed from the 14-amino-acid sunflower trypsin inhibitor (SFTI) peptide by phage display, yielding SFTI-based MASP inhibitors, SFMIs. Here, we present the crystal structure of the MASP-1/SFMI1 complex that we analyzed in comparison to other existing MASP-1/2 structures. Rigidified backbone structure has long been accepted as a structural prerequisite for peptide inhibitors of proteases. We found that a hydrophobic cluster organized around the P2 Thr residue is essential for the structural stability of wild-type SFTI. We also found that the same P2 Thr prevents binding of the rigid SFTI-like peptides to the substrate-binding cleft of both MASPs as the cleft is partially blocked by large gatekeeper enzyme loops. Directed evolution removed this obstacle by replacing the P2 Thr with a Ser, providing the SFMIs with high-degree structural plasticity, which proved to be essential for MASP inhibition. To gain more insight into the structural criteria for SFMI-based MASP-2 inhibition, we systematically modified MASP-2-specific SFMI2 by capping its two termini and by replacing its disulfide bridge with varying length thioether linkers. By doing so, we also aimed to generate a versatile scaffold that is resistant to reducing environment and has increased stability in exopeptidase-containing biological environments. We found that the reduction-resistant disulfide-substituted l-2,3-diaminopropionic acid (Dap) variant possessed near-native potency. As MASP-2 is involved in the life-threatening thrombosis in COVID-19 patients, our synthetic, selective MASP-2 inhibitors could be relevant coronavirus drug candidates.

Conformational selection mechanism governs oxygen ligation to H-NOX proteins.

H-NOX proteins are at present the only structural models available for the study of the ligand binding affinity and selectivity of soluble guanylate cyclase, the physiological receptor of nitric oxide. The oxy complex and resting state structures of two bacterial H-NOX proteins of markedly different oxygen affinity, but of quite similar sequence, were studied by molecular dynamics simulations at 300 K and 400 K. Unexpectedly, the different O(2) affinity was found to be reflected in differences of the resting states structures. A conformation containing a pre-formed oxygen-binding cage is the most populated in the resting state equilibrium ensemble of the only successful O(2) binder, Tt H-NOX at 300 K, suggesting that conformational selection governs the interaction.