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.

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.

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.

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.

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.

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.

Tissue-Specific Requirement for the GINS Complex During Zebrafish Development.

Efficient and accurate DNA replication is particularly critical in stem and progenitor cells for successful proliferation and survival. The replisome, an amalgam of protein complexes, is responsible for binding potential origins of replication, unwinding the double helix, and then synthesizing complimentary strands of DNA. According to current models, the initial steps of DNA unwinding and opening are facilitated by the CMG complex, which is composed of a GINS heterotetramer that connects Cdc45 with the mini-chromosome maintenance (Mcm) helicase. In this work, we provide evidence that in the absence of GINS function DNA replication is cell autonomously impaired, and we also show that gins1 and gins2 mutants exhibit elevated levels of apoptosis restricted to actively proliferating regions of the central nervous system (CNS). Intriguingly, our results also suggest that the rapid cell cycles during early embryonic development in zebrafish may not require the function of the canonical GINS complex as neither zygotic Gins1 nor Gins2 isoforms seem to be present during these stages.

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.

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.

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.

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 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.

Structural impact of GTP binding on downstream KRAS signaling.

Oncogenic RAS proteins, involved in ∼30% of human tumors, are molecular switches of various signal transduction pathways. Here we apply a new protocol for the NMR study of KRAS in its (inactive) GDP- and (activated) GTP-bound form, allowing a comprehensive analysis of the backbone dynamics of its WT-, G12C- and G12D variants. We found that Tyr32 shows opposite mobility with respect to the backbone of its surroundings: it is more flexible in the GDP-bound form while more rigid in GTP-complexes (especially in WT- and G12D-GTP). Using the G12C/Y32F double mutant, we showed that the presence of the hydroxyl group of Tyr32 has a marked effect on the G12C-KRAS-GTP system as well. Molecular dynamics simulations indicate that Tyr32 is linked to the γ-phosphate of GTP in the activated states - an arrangement shown, using QM/MM calculations, to support catalysis. Anchoring Tyr32 to the γ-phosphate contributes to the capture of the catalytic waters participating in the intrinsic hydrolysis of GTP and supports a simultaneous triple proton transfer step (catalytic water → assisting water → Tyr32 → O1G of the γ-phosphate) leading to straightforward product formation. The coupled flip of negatively charged residues of switch I toward the inside of the effector binding pocket potentiates ligand recognition, while positioning of Thr35 to enter the coordination sphere of the Mg2+ widens the pocket. Position 12 mutations do not disturb the capture of Tyr32 by the γ-phosphate, but (partially) displace Gln61, which opens up the catalytic pocket and destabilizes catalytic water molecules thus impairing intrinsic hydrolysis.

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.

Protein Aggregation in a Nutshell: The Splendid Molecular Architecture of the Dreaded Amyloid Fibrils.

The recent high-resolution structures of amyloid fibrils show that the organization of peptide segments into amyloid aggregate architecture is a general process, though the morphology is more complex and intricate than suspected previously. The amyloid fibrils are often cytotoxic, accumulating as intracellular inclusions or extracellular plaques and have the ability to interfere with cellular physiology causing various cellular malfunctions. At the same time, the highly ordered amyloid structures also present an opportunity for nature to store and protect peptide chains under extreme conditions - something that might be used for designing storage, formulation, and delivery of protein medications or for contriving bio-similar materials of great resistance or structure-ordering capacity. Here we summarize amyloid characteristics; discussing the basic morphologies, sequential requirements and 3D-structure that are required for the understanding of this newly (re)discovered protein structure - a prerequisite for developing either inhibitors or promoters of amyloid-forming processes.

Achieving Functionality Through Modular Build-up: Structure and Size Selection of Serine Oligopeptidases.

Enzymes of the prolyl oligopeptidase family (S9 family) recognize their substrates not only by the specificity motif to be cleaved but also by size - they hydrolyze oligopeptides smaller than 30 amino acids. They belong to the serine-protease family, but differ from classical serine-proteases in size (80 kDa), structure (two domains) and regulation system (size selection of substrates). This group of enzymes is an important target for drug design as they are linked to amnesia, schizophrenia, type 2 diabetes, trypanosomiasis, periodontitis and cell growth. By comparing the structure of various members of the family we show that the most important features contributing to selectivity and efficiency are: (i) whether the interactions weaving the two domains together play a role in stabilizing the catalytic triad and thus their absence may provide for its deactivation: these oligopeptidases can screen their substrates by opening up, and (ii) whether the interaction-prone ß-edge of the hydrolase domain is accessible and thus can guide a multimerization process that creates shielded entrance or intricate inner channels for the size-based selection of substrates. These cornerstones can be used to estimate the multimeric state and selection strategy of yet undetermined structures.

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.

Hydration shell differentiates folded and disordered states of a Trp-cage miniprotein, allowing characterization of structural heterogeneity by wide-line NMR measurements.

Hydration properties of folded and unfolded/disordered miniproteins were monitored in frozen solutions by wide-line 1H-NMR. The amount of mobile water as function of T (-80 °C < T < 0 °C) was found characteristically different for folded (TC5b), semi-folded (pH < 3, TCb5(H+)) and disordered (TC5b_N1R) variants. Comparing results of wide-line 1H-NMR and molecular dynamics simulations we found that both the amount of mobile water surrounding proteins in ice, as well as their thaw profiles differs significantly as function of the compactness and conformational heterogeneity of their structure. We found that (i) at around -50 °C ~50 H2Os/protein melt (ii) if the protein is well-folded then this amount of mobile water remains quasi-constant up to -20 °C, (iii) if disordered then the quantity of the lubricating mobile water increases with T in a constant manner up to ~200 H2Os/protein by reaching -20 °C. Especially in the -55 °C ↔ -15 °C temperature range, wide-line 1H-NMR detects the heterogeneity of protein fold, providing the size of the hydration shell surrounding the accessible conformers at a given temperature. Results indicate that freezing of protein solutions proceeds by the gradual selection of the enthalpically most favored states that also minimize the number of bridging waters.

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.