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.

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.

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.

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.

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.

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.

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.

Phosphorylation as conformational switch from the native to amyloid state: Trp-cage as a protein aggregation model.

The 20 residue long Trp-cage miniprotein is an excellent model for both computational and experimental studies of protein folding and stability. Recently, great attention emerged to study disease-related protein misfolding, aggregation, and amyloid formation, with the aim of revealing their structural and thermodynamic background. Trp-cage is sensitive to both environmental and structure-modifying effects. It aggregates with ease upon structure destabilization, and thus it is suitable for modeling aggregation and amyloid formation. Here, we characterize the amyloid formation of several sequence modified and side-chain phosphorylated Trp-cage variants. We applied NMR, circular dichroism (CD) and Fourier transform infrared (FTIR) spectroscopies, molecular dynamics (MD) simulations, and transmission electron microscopy (TEM) in conjunction with thioflavin-T (ThT) fluorescence measurements to reveal the structural consequences of side-chain phosphorylation. We demonstrate that the native fold is destabilized upon serine phosphorylation, and the resultant highly dynamic structures form amyloid-like ordered aggregates with high intermolecular ß-structure content. The only exception is the D9S(P) variant, which follows an alternative aggregation process by forming thin fibrils, presenting a CD spectrum of PPII helix, and showing low ThT binding capability. We propose a complex aggregation model for these Trp-cage miniproteins. This model assumes an additional aggregated state, a collagen triple helical form that can precede amyloid formation. The phosphorylation of a single serine residue serves as a conformational switch, triggering aggregation, otherwise mediated by kinases in cell. We show that Trp-cage miniprotein is indeed a realistic model of larger globular systems of composite folding and aggregation landscapes and helps us to understand the fundamentals of deleterious protein aggregation and amyloid formation.

Probing of primed and unprimed sites of calpains: Design, synthesis and evaluation of epoxysuccinyl-peptide derivatives as selective inhibitors.

Calpains are intracellular cysteine proteases with important physiological functions. Up- or downregulation of their expression can be responsible for several diseases, therefore specific calpain inhibitors may be considered as promising candidates for drug discovery. In this paper we describe the synthesis and characterization of a new class of inhibitors derived from the analysis of amino acid preferences in primed and unprimed sites of calpains by incorporation of l- or d-epoxysuccinyl group (Eps). Amino acids for replacement were chosen by considering the substrate preference of calpain 1 and 2 enzymes. The compounds were characterized by RP-HPLC, amino acid analysis and ESI-MS. Selectivity of the compounds was studied by using calpain 1 and 2; and cathepsin B. We have identified five calpain specific inhibitors with different extent of selectivity. Two of these also exhibited isoform selectivity. Compound NH2-Thr-Pro-Leu-(d-Eps)-Thr-Pro-Pro-Pro-Ser-NH2 proved to be a calpain 2 enzyme inhibitor with at least 11.8-fold selectivity, while compound NH2-Thr-Pro-Leu-(l-Eps)-Ser-Pro-Pro-Pro-Ser-NH2 possesses calpain 1 enzyme inhibition with at least 4-fold selectivity. The results of molecular modeling calculations suggest that the orientation of the bound inhibitor in the substrate binding cleft is markedly dependent on the stereochemistry of the epoxysuccinyl group.

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.

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.