Evaluation of an Affinity-Enhanced Anti-SARS-CoV2 Nanobody Design Workflow Using Machine Learning and Molecular Dynamics.

In silico optimization of protein binding has received a great deal of attention in the recent years. Since in silico prefiltering of strong binders is fast and cheap compared to in vitro library screening methods, the advent of powerful hardware and advanced machine learning algorithms has made this strategy more accessible and preferred. These advances have already impacted the global response to pandemic threats. In this study, we proposed and tested a workflow for designing nanobodies targeting the SARS-CoV-2 spike protein receptor binding domain (S-RBD) using machine learning techniques complemented by molecular dynamics simulations. We evaluated the feasibility of this workflow using a test set of 3 different nanobodies and 2 different S-RBD variants, from in silico design and bacterial expression to binding assays of the designed nanobody mutants. We successfully designed nanobodies that were subsequently tested against both the wild-type (Wuhan type) and the delta variant S-RBD and found 2 of them to be stronger binders compared to the wild-type nanobody. We use this case study to describe both the strengths and weaknesses of this in silico assisted nanobody design strategy.

Unravelling the Complexity of Amyloid Peptide Core Interfaces.

Amyloids, large intermolecular sandwiched ß-sheet structures, underlie several protein misfolding diseases but have also been shown to have functional roles and can be a basis for designing smart and responsive nanomaterials. Short segments of proteins, called aggregation-prone regions (APRs), have been identified that nucleate amyloid formation. Here we present the database of 173 APR crystal structures currently available in the PDB, and a tool, ACW, for analyzing their topologies and the 267 inter-ß-sheet interfaces of zipper regions assigned in these structures. We defined a new descriptor of zipper interfaces, the surface detail index (SDi), which quantifies the intertwining between the side chains of both ß-sheets of the zipper, an important factor for the molecular recognition and self-assembly of these mesostructures. This allowed a comparative analysis of the zipper interfaces and identification of 6 clusters with different intertwining, steric fit, and size characteristics using three complementary descriptors, SDi, shape complementarity, and buried surface area. 60% of the APR structures are formed by parallel ß-sheets, of which 52% belong to the topological class 1. This could be explained by the better fit and a deeper entanglement of the zipper regions of the parallel structures than of the antiparallel structures, as the analysis showed that both their shape complementarity (0.79 vs 0.70) and SDi (1.53 vs 1.32) were higher. The higher abundance of certain residues (Asn and Gln in parallel and Leu and Ala in antiparallel ß-sheets) can be explained by their ability to form different ladder-like secondary interaction patterns within ß-sheets. Analogous to the hierarchy of protein structure, we interpreted the primary, secondary, tertiary, and quaternary structure levels of APRs revealing different characteristics of the zipper regions for both parallel and antiparallel ß-sheet structures, which may provide clues to the structural conditions of amyloid core formation and the rational design of amyloid polymorphs.

Influence of Trp-Cage on the Function and Stability of GLP-1R Agonist Exenatide Derivatives.

Exenatide (Ex4), a GLP-1 incretin mimetic polypeptide, is an effective therapeutic agent against diabetes and obesity. We highlight the indirect role of Ex4's structure-stabilizing Trp-cage (Tc) motif in governing GLP-1 receptor (GLP-1R) signal transduction. We use various Ex4 derivatives to explore how Tc compactness influences thermal stability, aggregation, enhancement of insulin secretion, and GLP-1R binding. We found that Ex4 variants decorated with fortified Tc motifs exhibit increased resistance to unfolding and aggregation but show an inverse relationship between the bioactivity and stability. Molecular dynamics simulations coupled with a rigid-body segmentation protocol to analyze dynamic interconnectedness revealed that the constrained Tc motifs remain intact within the receptor-ligand complexes but interfere with one of the major stabilizing contacts and recognition loci on the extracellular side of GLP-1R, dislodging the N-terminal activating region of the hormone mimetics, and restrict the free movement of TM6, the main signal transduction device of GLP-1R.

Target-Templated Construction of Functional Proteomimetics Using Photo-Foldamer Libraries.

Current methods for proteomimetic engineering rely on structure-based design. Here we describe a design strategy that allows the construction of proteomimetics against challenging targets without a priori characterization of the target surface. Our approach employs (i) a 100-membered photoreactive foldamer library, the members of which act as local surface mimetics, and (ii) the subsequent affinity maturation of the primary hits using systems chemistry. Two surface-oriented proteinogenic side chains drove the interactions between the short helical foldamer fragments and the proteins. Diazirine-based photo-crosslinking was applied to sensitively detect and localize binding even to shallow and dynamic patches on representatively difficult targets. Photo-foldamers identified functionally relevant protein interfaces, allosteric and previously unexplored targetable regions on the surface of STAT3 and an oncogenic K-Ras variant. Target-templated dynamic linking of foldamer hits resulted in two orders of magnitude affinity improvement in a single step. The dimeric K-Ras ligand mimicked protein-like catalytic functions. The photo-foldamer approach thus enables the highly efficient mapping of protein-protein interaction sites and provides a viable starting point for proteomimetic ligand development without a priori structural hypotheses.

Solvent induced amyloid polymorphism and the uncovering of the elusive class 3 amyloid topology.

Aggregation-prone-motifs (APRs) of proteins are short segments, which - as isolated peptides - form diverse amyloid-like crystals. We introduce two APRs - designed variants of the incretin mimetic Exendin-4 - that both display crystal-phase polymorphism. Crystallographic and spectroscopic analysis revealed that a single amino-acid substitution can greatly reduce topological variability: while LYIQWL can form both parallel and anti-parallel ß-sheets, LYIQNL selects only the former. We also found that the parallel/anti-parallel switch of LYIQWL can be induced by simply changing the crystallization temperature. One crystal form of LYIQNL was found to belong to the class 3 topology, an arrangement previously not encountered among proteinogenic systems. We also show that subtle environmental changes lead to crystalline assemblies with different topologies, but similar interfaces. Spectroscopic measurements showed that polymorphism is already apparent in the solution state. Our results suggest that the temperature-, sequence- and environmental sensitivity of physiological amyloids is reflected in assemblies of the APR segments, which, complete with the new class 3 crystal form, effectively sample all the originally proposed basic topologies of amyloid-like aggregates.

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.

Mapping protein binding sites by photoreactive fragment pharmacophores.

Fragment screening is a popular strategy of generating viable chemical starting points especially for challenging targets. Although fragments provide a better coverage of chemical space and they have typically higher chance of binding, their weak affinity necessitates highly sensitive biophysical assays. Here, we introduce a screening concept that combines evolutionary optimized fragment pharmacophores with the use of a photoaffinity handle that enables high hit rates by LC-MS-based detection. The sensitivity of our screening protocol was further improved by a target-conjugated photocatalyst. We have designed, synthesized, and screened 100 diazirine-tagged fragments against three benchmark and three therapeutically relevant protein targets of different tractability. Our therapeutic targets included a conventional enzyme, the first bromodomain of BRD4, a protein-protein interaction represented by the oncogenic KRasG12D protein, and the yet unliganded N-terminal domain of the STAT5B transcription factor. We have discovered several fragment hits against all three targets and identified their binding sites via enzymatic digestion, structural studies and modeling. Our results revealed that this protocol outperforms screening traditional fully functionalized and photoaffinity fragments in better exploration of the available binding sites and higher hit rates observed for even difficult targets.

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.

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.

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.

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.

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.

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.

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.