Contribution of Noncovalent Recognition and Reactivity to the Optimization of Covalent Inhibitors: A Case Study on KRasG12C.

Covalent drugs might bear electrophiles to chemically modify their targets and have the potential to target previously undruggable proteins with high potency. Covalent binding of drug-size molecules includes a noncovalent recognition provided by secondary interactions and a chemical reaction leading to covalent complex formation. Optimization of their covalent mechanism of action should involve both types of interactions. Noncovalent and covalent binding steps can be characterized by an equilibrium dissociation constant (KI) and a reaction rate constant (kinact), respectively, and they are affected by both the warhead and the scaffold of the ligand. The relative contribution of these two steps was investigated on a prototypic drug target KRASG12C, an oncogenic mutant of KRAS. We used a synthetically more accessible nonchiral core derived from ARS-1620 that was equipped with four different warheads and a previously described KRAS-specific basic side chain. Combining these structural changes, we have synthesized novel covalent KRASG12C inhibitors and tested their binding and biological effect on KRASG12C by various biophysical and biochemical assays. These data allowed us to dissect the effect of scaffold and warhead on the noncovalent and covalent binding event. Our results revealed that the atropisomeric core of ARS-1620 is not indispensable for KRASG12C inhibition, the basic side chain has little effect on either binding step, and warheads affect the covalent reactivity but not the noncovalent binding. This type of analysis helps identify structural determinants of efficient covalent inhibition and may find use in the design of covalent agents.

Opening Amyloid-Windows to the secondary structure of proteins: The amyloidogenecity increases tenfold inside beta-sheets.

Methods from artificial intelligence (AI), in general, and machine learning, in particular, have kept conquering new territories in numerous areas of science. Most of the applications of these techniques are restricted to the classification of large data sets, but new scientific knowledge can seldom be inferred from these tools. Here we show that an AI-based amyloidogenecity predictor can strongly differentiate the border- and the internal hexamers of ß-pleated sheets when screening all the Protein Data Bank-deposited homology-filtered protein structures. Our main result shows that more than 30% of internal hexamers of ß sheets are predicted to be amyloidogenic, while just outside the border regions, only 3% are predicted as such. This result may elucidate a general protection mechanism of proteins against turning into amyloids: if the borders of ß-sheets were amyloidogenic, then the whole ß sheet could turn more easily into an insoluble amyloid-structure, characterized by periodically repeated parallel ß-sheets. We also present that no analogous phenomenon exists on the borders of α-helices or randomly chosen subsequences of the studied protein structures.

LoCoHD: a metric for comparing local environments of proteins.

Protein folds and the local environments they create can be compared using a variety of differently designed measures, such as the root mean squared deviation, the global distance test, the template modeling score or the local distance difference test. Although these measures have proven to be useful for a variety of tasks, each fails to fully incorporate the valuable chemical information inherent to atoms and residues, and considers these only partially and indirectly. Here, we develop the highly flexible local composition Hellinger distance (LoCoHD) metric, which is based on the chemical composition of local residue environments. Using LoCoHD, we analyze the chemical heterogeneity of amino acid environments and identify valines having the most conserved-, and arginines having the most variable chemical environments. We use LoCoHD to investigate structural ensembles, to evaluate critical assessment of structure prediction (CASP) competitors, to compare the results with the local distance difference test (lDDT) scoring system, and to evaluate a molecular dynamics simulation. We show that LoCoHD measurements provide unique information about protein structures that is distinct from, for example, those derived using the alignment-based RMSD metric, or the similarly distance matrix-based but alignment-free lDDT metric.

Unveiling the Oxazolidine Character of Pseudoproline Derivatives by Automated Flow Peptide Chemistry.

Pseudoproline derivatives such as Thr(ΨPro)-OH are commonly used in peptide synthesis to reduce the likelihood of peptide aggregation and to prevent aspartimide (Asi) formation during the synthesis process. In this study, we investigate notable by-products such as aspartimide formation and an imine derivative of the Thr(ΨPro) moiety observed in flow peptide chemistry synthesis. To gain insight into the formation of these unexpected by-products, we design a series of experiments. Furthermore, we demonstrate the oxazolidine character of the pseudoproline moiety and provide plausible mechanisms for the two-way ring opening of oxazolidine leading to these by-products. In addition, we present evidence that Asi formation appears to be catalyzed by the presence of the pseudoproline moiety. These observed side reactions are attributed to elevated temperature and pressure; therefore, caution is advised when using ΨPro derivatives under such harsh conditions. In addition, we propose a solution whereby thermodynamically controlled Asi formation can be kinetically prevented.

Amide isomerization pathways: Electronic and structural background of protonation- and deprotonation-mediated cis-trans interconversions.

The cis-trans isomerization of amide bonds leads to wide range of structural and functional changes in proteins and can easily be the rate-limiting step in folding. The trans isomer is thermodynamically more stable than the cis, nevertheless the cis form can play a role in biopolymers' function. The molecular system of N-methylacetamide · 2H2O is complex enough to reveal energetics of the cis-trans isomerization at coupled cluster single-double and coupled cluster single-double and perturbative triple [CCSD(T)] levels of theory. The cis-trans isomerization cannot be oversimplified by a rotation along ω, since this rotation is coupled with the N-atom pyramidal inversion, requesting the introduction of a second dihedral angle "α." Full f(ω,α) potential energy surfaces of the different amide protonation states, critical points and isomerization reaction paths were determined, and the barriers of the neutral, O-protonated and N-deprotonated amides were found too high to allow cis-trans interconversion at room temperature: ∼85, ∼140, and ∼110 kJ mol-1, respectively. For the N-protonated amide bond, the cis form (ω = 0°) is a maximum rather than a minimum, and each ω state is accessible for less than ∼10 kJ mol-1. Here we outline a cis-trans isomerization pathway with a previously undescribed low energy transition state, which suggests that the proton is transferred from the more favorable O- to the N-protonation site with the aid of nearby water molecules, allowing the trans → cis transition to occur at an energy cost of ≤11.6 kJ mol-1. Our results help to explain why isomerase enzymes operate via protonated amide bonds and how N-protonation of the peptide bond occurs via O-protonation.

N6-Methyladenine Progressively Accumulates in Mitochondrial DNA during Aging.

N6-methyladenine (6mA) in the DNA is a conserved epigenetic mark with various cellular, physiological and developmental functions. Although the presence of 6mA was discovered a few years ago in the nuclear genome of distantly related animal taxa and just recently in mammalian mitochondrial DNA (mtDNA), accumulating evidence at present seriously questions the presence of N6-adenine methylation in these genetic systems, attributing it to methodological errors. In this paper, we present a reliable, PCR-based method to determine accurately the relative 6mA levels in the mtDNA of Caenorhabditis elegans, Drosophila melanogaster and dogs, and show that these levels gradually increase with age. Furthermore, daf-2(-)-mutant worms, which are defective for insulin/IGF-1 (insulin-like growth factor) signaling and live twice as long as the wild type, display a half rate at which 6mA progressively accumulates in the mtDNA as compared to normal values. Together, these results suggest a fundamental role for mtDNA N6-adenine methylation in aging and reveal an efficient diagnostic technique to determine age using DNA.

Polymorphic amyloid nanostructures of hormone peptides involved in glucose homeostasis display reversible amyloid formation.

A large group of hormones are stored as amyloid fibrils in acidic secretion vesicles before they are released into the bloodstream and readopt their functional state. Here, we identify an evolutionarily conserved hexapeptide sequence as the major aggregation-prone region (APR) of gastrointestinal peptides of the glucagon family: xFxxWL. We determine nine polymorphic crystal structures of the APR segments of glucagon-like peptides 1 and 2, and exendin and its derivatives. We follow amyloid formation by CD, FTIR, ThT assays, and AFM. We propose that the pH-dependent changes of the protonation states of glutamate/aspartate residues of APRs initiate switching between the amyloid and the folded, monomeric forms of the hormones. We find that pH sensitivity diminishes in the absence of acidic gatekeepers and amyloid formation progresses over a broad pH range. Our results highlight the dual role of short aggregation core motifs in reversible amyloid formation and receptor binding.

NMR-Chemical-Shift-Driven Protocol Reveals the Cofactor-Bound, Complete Structure of Dynamic Intermediates of the Catalytic Cycle of Oncogenic KRAS G12C Protein and the Significance of the Mg2+ Ion.

In this work, catalytically significant states of the oncogenic G12C variant of KRAS, those of Mg2+-free and Mg2+-bound GDP-loaded forms, have been determined using CS-Rosetta software and NMR-data-driven molecular dynamics simulations. There are several Mg2+-bound G12C KRAS/GDP structures deposited in the Protein Data Bank (PDB), so this system was used as a reference, while the structure of the Mg2+-free but GDP-bound state of the RAS cycle has not been determined previously. Due to the high flexibility of the Switch-I and Switch-II regions, which also happen to be the catalytically most significant segments, only chemical shift information could be collected for the most important regions of both systems. CS-Rosetta was used to derive an "NMR ensemble" based on the measured chemical shifts, which, however, did not contain the nonprotein components of the complex. We developed a torsional restraint set for backbone torsions based on the CS-Rosetta ensembles for MD simulations, overriding the force-field-based parametrization in the presence of the reinserted cofactors. This protocol (csdMD) resulted in complete models for both systems that also retained the structural features and heterogeneity defined by the measured chemical shifts and allowed a detailed comparison of the Mg2+-bound and Mg2+-free states of G12C KRAS/GDP.

Inhibitor Design Strategy for Myostatin: Dynamics and Interaction Networks Define the Affinity and Release Mechanisms of the Inhibited Complexes.

Myostatin, an important negative regulator of muscle mass, is a therapeutic target for muscle atrophic disorders such as muscular dystrophy. Thus, the inhibition of myostatin presents a strategy to treat these disorders. It has long been established that the myostatin prodomain is a strong inhibitor of the mature myostatin, and the minimum peptide of the prodomain-corresponding to the α1-helix of its lasso-region-responsible for the inhibitory efficiency was defined and characterized as well. Here we show that the minimum peptide segment based on the growth differentiation factor 11 (GDF11), which we found to be more helical in its stand-alone solvated stfate than the similar segment of myostatin, is a promising new base scaffold for inhibitor design. The proposed inhibitory peptides in their solvated state and in complex with the mature myostatin were analyzed by in silico molecule modeling supplemented with the electronic circular dichroism spectroscopy measurements. We defined the Gaussian-Mahalanobis mean score to measure the fraction of dihedral angle-pairs close to the desired helical region of the Ramachandran-plot, carried out RING analysis of the peptide-protein interaction networks and characterized the internal motions of the complexes using our rigid-body segmentation protocol. We identified a variant-11m2-that is sufficiently ordered both in solvent and within the inhibitory complex, forms a high number of contacts with the binding-pocket and induces such changes in its internal dynamics that lead to a rigidified, permanently locked conformation that traps this peptide in the binding site. We also showed that the naturally evolved α1-helix has been optimized to simultaneously fulfill two very different roles: to function as a strong binder as well as a good leaving group. It forms an outstanding number of non-covalent interactions with the mature core of myostatin and maintains the most ordered conformation within the complex, while it induces independent movement of the gate-keeper ß-hairpin segment assisting the dissociation and also results in the least-ordered solvated form which provides extra stability for the dissociated state and discourages rebinding.

Downregulation of transposable elements extends lifespan in Caenorhabditis elegans.

Mobility of transposable elements (TEs) frequently leads to insertional mutations in functional DNA regions. In the potentially immortal germline, TEs are effectively suppressed by the Piwi-piRNA pathway. However, in the genomes of ageing somatic cells lacking the effects of the pathway, TEs become increasingly mobile during the adult lifespan, and their activity is associated with genomic instability. Whether the progressively increasing mobilization of TEs is a cause or a consequence of ageing remains a fundamental problem in biology. Here we show that in the nematode Caenorhabditis elegans, the downregulation of active TE families extends lifespan. Ectopic activation of Piwi proteins in the soma also promotes longevity. Furthermore, DNA N6-adenine methylation at TE stretches gradually rises with age, and this epigenetic modification elevates their transcription as the animal ages. These results indicate that TEs represent a novel genetic determinant of ageing, and that N6-adenine methylation plays a pivotal role in ageing control.

Acetyl group for proper protection of ß-sugar-amino acids used in SPPS.

The synthesis of D-glucosamine-1-carboxylic acid based ß-sugar amino acids (ß-SAAs) is typically performed in nine consecutive steps via an inefficient OAc → Br → CN conversion protocol with low overall yield. Here, we present the improved and more efficient synthesis of both Fmoc-GlcAPC-OH and Fmoc-GlcAPC(Ac)-OH, ß-SAAs consisting of only 4-5 synthetic steps. Their active ester and amide bond formation with glycine methyl ester (H-Gly-OMe) was completed and monitored by 1H NMR. The stability of the pyranoid OHs protecting the acetyl groups was investigated under three different Fmoc cleavage conditions and was found to be satisfactory even at high piperidine concentration (e.g. 40%). We designed a SPPS protocol using Fmoc-GlcAPC(Ac)-OH to produce model peptides Gly-ß-SAA-Gly as well as Gly-ß-SAA-ß-SAA-Gly with high coupling efficiency. The products were deacetylated using the Zemplén method, which allows the hydrophilicity of a building block and/or chimera to be fine-tuned, even after the polypeptide chain has already been synthesized.

A Novel Fusion Protein System for the Production of Nanobodies and the SARS-CoV-2 Spike RBD in a Bacterial System.

Antibodies are key proteins of the immune system, and they are widely used for both research and theragnostic applications. Among them, camelid immunoglobulins (IgG) differ from the canonical human IgG molecules, as their light chains are completely missing; thus, they have only variable domains on their heavy chains (VHHs). A single VHH domain, often called a nanobody, has favorable structural, biophysical, and functional features compared to canonical antibodies. Therefore, robust and efficient production protocols relying on recombinant technologies are in high demand. Here, by utilizing ecotin, an Escherichia coli protein, as a fusion partner, we present a bacterial expression system that allows an easy, fast, and cost-effective way to prepare nanobodies. Ecotin was used here as a periplasmic translocator and a passive refolding chaperone, which allowed us to reach high-yield production of nanobodies. We also present a new, easily applicable prokaryotic expression and purification method of the receptor-binding domain (RBD) of the SARS-CoV-2 S protein for interaction assays. We demonstrate using ECD spectroscopy that the bacterially produced RBD is well-folded. The bacterially produced nanobody was shown to bind strongly to the recombinant RBD, with a Kd of 10 nM. The simple methods presented here could facilitate rapid interaction measurements in the event of the appearance of additional SARS-CoV-2 variants.

Covalent fragment mapping of KRasG12C revealed novel chemotypes with in vivo potency.

G12C mutant KRas is considered druggable by allele-specific covalent inhibitors due to the nucleophilic character of the oncogenic mutant cysteine at position 12. Discovery of these inhibitors requires the optimization of both covalent and noncovalent interactions. Here, we report covalent fragment screening of our electrophilic fragment library of diverse non-covalent scaffolds equipped with 40 different electrophilic functionalities to identify fragments as suitable starting points targeting Cys12. Screening the library against KRasG12C using Ellman's free thiol assay, followed by protein NMR and cell viability assays, resulted in two potential inhibitor chemotypes. Characterization of these scaffolds in in vitro cellular- and in vivo xenograft models revealed them as promising starting points for covalent drug discovery programs.

[Farnesyltransferase inhibitors have antitumoral effects in mutant KRAS containing cancer cells in preclinical models]. / Farnezil-transzferáz-inhibitorok (FTI) szinergista módon potencírozzák a mutáns KRAS inhibitorai daganatellenes hatásait preklinikai modellekben.

In silico studies raised the possibility that farnesyltransferase inhibitors (FTIs) may have antitumoral effects on KRAS mutant cancer cells. Accordingly, we have tested FTIs (tipifarnib and lonafarnib) in G12C mutant human cancer cell lines in vitro and in vivo. We have discovered that the combination of the two drugs has a synergistic antitumoral effect. Next, we have tested FTIs on G12D mutant human cancer cell lines and found that the combination has antitumoral effect in various preclinical cancer models. At last, we have also tested FTIs on G12V mutant human cancer cells and again we have detected antitumoral effects. We suggest that FTIs may have clinical relevance outside the HRAS mutant cancers.

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.

Succinct Amyloid and Nonamyloid Patterns in Hexapeptides.

Hexapeptides are widely applied as a model system for studying the amyloid-forming properties of polypeptides, including proteins. Recently, large experimental databases have become publicly available with amyloidogenic labels. Using these data sets for training and testing purposes, one may build artificial intelligence (AI)-based classifiers for predicting the amyloid state of peptides. In our previous work (Biomolecules 2021, 11, 500), we described the Support Vector Machine (SVM)-based Budapest Amyloid Predictor (https://pitgroup.org/bap). Here, we apply the Budapest Amyloid Predictor for discovering numerous amyloidogenic and nonamyloidogenic hexapeptide patterns with accuracy between 80% and 84%, as surprising and succinct novel rules for further understanding the amyloid state of peptides. For example, we have shown that for any independently mutated residue (position marked by "x"), the patterns CxFLWx, FxFLFx, or xxIVIV are predicted to be amyloidogenic, while those of PxDxxx, xxKxEx, and xxPQxx are nonamyloidogenic. We note that each amyloidogenic pattern with two x's (e.g.,CxFLWx) describes succinctly 202 = 400 hexapeptides, while the nonamyloidogenic patterns comprising four point mutations (e.g.,PxDxxx) give 204 = 160 000 hexapeptides in total. We also examine the restricted substitutions for positions "x" from subclasses of proteinogenic amino acid residues; for example, if "x" is substituted with hydrophobic amino acids, then there exist patterns containing three x's, like MxVVxx, predicted to be amyloidogenic. If we can choose for the x positions any hydrophobic amino acids, except the "structure breaker" proline, then we get amyloid patterns with five x positions, for example, xxxFxx, each corresponding to 32 768 hexapeptides. To our knowledge, no similar applications of artificial intelligence tools or succinct amyloid patterns were described before the present work.

D27-LIKE1 isomerase has a preference towards trans/cis and cis/cis conversions of carotenoids in Arabidopsis.

Carotenoids contribute to a variety of physiological processes in plants, functioning also as biosynthesis precursors of ABA and strigolactones (SLs). SL biosynthesis starts with the enzymatic conversion of all-trans-ß-carotene to 9-cis-ß-carotene by the DWARF27 (D27) isomerase. In Arabidopsis, D27 has two closely related paralogs, D27-LIKE1 and D27-LIKE2, which were predicted to be ß-carotene-isomerases. In the present study, we characterised D27-LIKE1 and identified some key aspects of its physiological and enzymatic functions in Arabidopsis. d27-like1-1 mutant does not display any strigolactone-deficient traits and exhibits a substantially higher 9-cis-violaxanthin content, which is accompanied by a slightly higher ABA level. In vitro feeding assays with recombinant D27-LIKE1 revealed that the protein exhibits affinity to all ß-carotene isoforms but with an exclusive preference towards trans/cis conversions and the interconversion between 9-cis, 13-cis and 15-cis-ß-carotene forms, and accepts zeaxanthin and violaxanthin as substrates. Finally, we present evidence showing that D27-LIKE1 mRNA is phloem mobile and D27-LIKE1 is an ancient isomerase with a long evolutionary history. In summary, we demonstrate that D27-LIKE1 is a carotenoid isomerase with multi-substrate specificity and has a characteristic preference towards the catalysation of cis/cis interconversion of carotenoids. Therefore, D27-LIKE1 is a potential regulator of carotenoid cis pools and, eventually, SL and ABA biosynthesis pathways.

Synthesis of the extracellular domain of GLP-1R by chemical and biotechnological approaches.

The extracellular domain of the glucagon-like peptide-1 receptor, GLP-1R, is responsible for the binding of GLP-1, and a handful of additional agonists (such as exenatide, lixisenatide, and liraglutide) used daily for treating type II diabetes mellitus. Lead discovery and optimization, however, require binding studies, which, in turn, necessitate the total synthesis of GLP-1R, comprising 108 residues. A protein domain of 10-15 kDa size could be obtained either by expression in E. coli or by ligating solid-phase peptide synthesis (SPPS)-made fragments. However, direct overexpression fails to give a properly folded protein, as GLP-1R forms an inclusion body, which fails to refold due to improper disulfide pairing. Several bacterial strains, constructs, and fusion partners were probed and it was found that only co-expression with MBP gave a 3D-fold allowing the native disulfide bond pattern formation. Some fusion partners can act as covalently linked or in situ chaperones for guiding the refolding of GLP-1R toward success. Therefore, the bottleneck to preparing GPCR extracellular domains is the correct pairing of the Cys residues. As a proof-of-concept model, nGLP1-R was made by SPPS to form the purified full-length polypeptide chain, subjected to self-guided or spontaneous Cys pairing. However, the formation of correct SS-pairs was lagging behind any protocol in use support, and the bottleneck of large-scale protein production relies on the risky step of proper refolding, which is sometimes possible only if a suitable fusion partner effectively helps and catalysis of the correct disulfide formation.

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.