Distinct conformational and energetic features define the specific recognition of (di)aromatic peptide motifs by PEX14.

The cycling import receptor PEX5 and its membrane-located binding partner PEX14 are key constituents of the peroxisomal import machinery. Upon recognition of newly synthesized cargo proteins carrying a peroxisomal targeting signal type 1 (PTS1) in the cytosol, the PEX5/cargo complex docks at the peroxisomal membrane by binding to PEX14. The PEX14 N-terminal domain (NTD) recognizes (di)aromatic peptides, mostly corresponding to Wxxx(F/Y)-motifs, with nano-to micromolar affinity. Human PEX5 possesses eight of these conserved motifs distributed within its 320-residue disordered N-terminal region. Here, we combine biophysical (ITC, NMR, CD), biochemical and computational methods to characterize the recognition of these (di)aromatic peptides motifs and identify key features that are recognized by PEX14. Notably, the eight motifs present in human PEX5 exhibit distinct affinities and energetic contributions for the interaction with the PEX14 NTD. Computational docking and analysis of the interactions of the (di)aromatic motifs identify the specific amino acids features that stabilize a helical conformation of the peptide ligands and mediate interactions with PEX14 NTD. We propose a refined consensus motif ExWΦxE(F/Y)Φ for high affinity binding to the PEX14 NTD and discuss conservation of the (di)aromatic peptide recognition by PEX14 in other species.

Benchmarking biomolecular force field-based Zn2+ for mono- and bimetallic ligand binding sites.

Zn2+ is one of the most versatile biologically available metal ions, but accurate modeling of Zn2+ -containing metalloproteins at the biomolecular force field level can be challenging. Since most Zn2+ models are parameterized in bulk solvent, in-depth knowledge about their performance in a protein environment is limited. Thus, we systematically investigate here the behavior of non-polarizable Zn2+ models for their ability to reproduce experimentally determined metal coordination and ligand binding in metalloproteins. The benchmarking is performed in challenging environments, including mono- (carbonic anhydrase II) and bimetallic (metallo-ß-lactamase VIM-2) ligand binding sites. We identify key differences in the performance between the Zn2+ models with regard to the preferred ligating atoms (charged/non-charged), attraction of water molecules, and the preferred coordination geometry. Based on these results, we suggest suitable simulation conditions for varying Zn2+ site geometries that could guide the further development of biomolecular Zn2+ models.

Comparison of allosteric signaling in DnaK and BiP using mutual information between simulated residue conformations.

The heat shock protein 70 kDa (Hsp70) chaperone system serves as a critical component of protein quality control across a wide range of prokaryotic and eukaryotic organisms. Divergent evolution and specialization to particular organelles have produced numerous Hsp70 variants which share similarities in structure and general function, but differ substantially in regulatory aspects, including conformational dynamics and activity modulation by cochaperones. The human Hsp70 variant BiP (also known as GRP78 or HSPA5) is of therapeutic interest in the context of cancer, neurodegenerative diseases, and viral infection, including for treatment of the pandemic virus SARS-CoV-2. Due to the complex conformational rearrangements and high sequential variance within the Hsp70 protein family, it is in many cases poorly understood which amino acid mutations are responsible for biochemical differences between protein variants. In this study, we predicted residues associated with conformational regulation of human BiP and Escherichia coli DnaK. Based on protein structure networks obtained from molecular dynamics simulations, we analyzed the shared information between interaction timelines to highlight residue positions with strong conformational coupling to their environment. Our predictions, which focus on the binding processes of the chaperone's substrate and cochaperones, indicate residues filling potential signaling roles specific to either DnaK or BiP. By combining predictions of individual residues into conformationally coupled chains connecting ligand binding sites, we predict a BiP specific secondary signaling pathway associated with substrate binding. Our study sheds light on mechanistic differences in signaling and regulation between Hsp70 variants, which provide insights relevant to therapeutic applications of these proteins.

Dynamic Docking of Macrocycles in Bound and Unbound Protein Structures with DynaDock.

Macrocycles are interesting molecules with unique features due to their conformationally constrained yet flexible ring structure. This characteristic poses a difficult challenge for computational modeling studies since they rely on accurate structural descriptions. In particular, molecular docking calculations suffer from the lack of ring flexibility during pose generation, which is often compensated by using pregenerated ligand conformer ensembles. Moreover, receptor structures are mainly treated rigidly, which limits the use of many docking tools. In this study, we optimized our previous molecular dynamics-based sampling and docking pipeline specifically designed for the accurate prediction of macrocyclic compounds. We developed a dihedral classification procedure for in-depth conformational analysis of the macrocyclic rings and extracted structural ensembles that were subsequently docked in both bound and unbound protein structures employing a fully flexible approach. Our results suggest that including a ring conformer close to the bound state in the starting ensemble increases the chance of successful docking. The bioactive conformations of a diverse set of ligands could be predicted with high and decent accuracy in bound and unbound protein structures, respectively, due to the incorporation of full molecular flexibility in our approach. The remaining unsuccessful docking calculations were mainly caused by large flexible substituents that bind to surface-exposed binding sites, rather than the macrocyclic ring per se and could be further improved by explicit molecular dynamics simulations of the docked complex.

Investigation of the pH-dependent aggregation mechanisms of GCSF using low resolution protein characterization techniques and advanced molecular dynamics simulations.

Granulocyte-colony stimulating factor (GCSF) is a widely used therapeutic protein to treat neutropenia. GCSF has an increased propensity to aggregate if the pH is increased above 5.0. Although GCSF is very well experimentally characterized, the exact pH-dependent aggregation mechanism of GCSF is still under debate. This study aimed to model the complex pH-dependent aggregation behavior of GCSF using state-of-the-art simulation techniques. The conformational stability of GCSF was investigated by performing metadynamics simulations, while the protein-protein interactions were investigated using coarse-grained (CG) simulations of multiple GCSF monomers. The CG simulations were directly compared with small-angle X-ray (SAXS) data. The metadynamics simulations demonstrated that the orientations of Trp residues in GCSF are dependent on pH. The conformational change of Trp residues is due to the loss of Trp-His interactions at the physiological pH, which in turn may increase protein flexibility. The helical structure of GCSF was not affected by the pH conditions of the simulations. Our CG simulations indicate that at pH 4.0, the colloidal stability may be more important than the conformational stability of GCSF. The electrostatic potential surface and CG simulations suggested that the basic residues are mainly responsible for colloidal stability as deprotonation of these residues causes a reduction of the highly positively charged electrostatic barrier close to the aggregation-prone long loop regions.

SenseNet, a tool for analysis of protein structure networks obtained from molecular dynamics simulations.

Computational methods play a key role for investigating allosteric mechanisms in proteins, with the potential of generating valuable insights for innovative drug design. Here we present the SenseNet ("Structure ENSEmble NETworks") framework for analysis of protein structure networks, which differs from established network models by focusing on interaction timelines obtained by molecular dynamics simulations. This approach is evaluated by predicting allosteric residues reported by NMR experiments in the PDZ2 domain of hPTP1e, a reference system for which previous computational predictions have shown considerable variance. We applied two models based on the mutual information between interaction timelines to estimate the conformational influence of each residue on its local environment. In terms of accuracy our prediction model is comparable to the top performing model published for this system, but by contrast benefits from its independence from NMR structures. Our results are complementary to experimental data and the consensus of previous predictions, demonstrating the potential of our new analysis tool SenseNet. Biochemical interpretation of our model suggests that allosteric residues in the PDZ2 domain form two distinct clusters of contiguous sidechain surfaces. SenseNet is provided as a plugin for the network analysis software Cytoscape, allowing for ease of future application and contributing to a system of compatible tools bridging the fields of system and structural biology.

Structure-Guided Modulation of the Catalytic Properties of [2Fe-2S]-Dependent Dehydratases.

The FeS cluster-dependent dihydroxyacid dehydratases (DHADs) and sugar acid-specific dehydratases (DHTs) from the ilvD/EDD superfamily are key enzymes in the bioproduction of a wide variety of chemicals. We analyzed [2Fe-2S]-dependent dehydratases in silico and in vitro, deduced functionally relevant sequence, structure, and activity relationships within the ilvD/EDD superfamily, and we propose a new classification based on their evolutionary relationships and substrate profiles. In silico simulations and analyses identified several key positions for specificity, which were experimentally investigated with site-directed and saturation mutagenesis. We thus increased the promiscuity of DHAD from Fontimonas thermophila (FtDHAD), showing >10-fold improved activity toward D-gluconate, and shifted the substrate preference of DHT from Paralcaligenes ureilyticus (PuDHT) toward shorter sugar acids (recording a six-fold improved activity toward the non-natural substrate D-glycerate). The successful elucidation of the role of important active site residues of the ilvD/EDD superfamily will further guide developments of this important biocatalyst for industrial applications.

The Q41R mutation in the HCV-protease enhances the reactivity towards MAVS by suppressing non-reactive pathways.

Recent experimental findings pointed out a new mutation in the HCV protease, Q41R, responsible for a significant enhancement of the enzyme's reactivity towards the mitochondrial antiviral-signaling protein (MAVS). The Q41R mutation is located rather far from the active site, and its involvement in the overall reaction mechanism is thus unclear. We used classical molecular dynamics and QM/MM to study the acylation reaction of HCV NS3/4A protease variants bound to MAVS and the NS4A/4B substrate and uncovered the indirect mechanism by which the Q41R mutation plays a critical role in the efficient cleavage of the substrate. Our simulations reveal that there are two major conformations of the MAVS H1'(p) residue for the wild type protease and only one conformation for the Q41R mutant. The conformational space of H1'(p) is restricted by the Q41R mutation due to a π-π stacking between H1'(p) and R41 as well as a strong hydrogen bond between the backbone of H57 and the side chain of R41. Further QM/MM calculations indicate that the complex with the conformation ruled out by the Q41R substitution is a non-reactive species due to its higher free energy barrier for the acylation reaction. Based on our calculations, we propose a kinetic mechanism that explains experimental data showing an increase of apparent rate constants for MAVS cleavage in Q41R mutants. Our model predicts that the non-reactive conformation of the enzyme-substrate complex modulates reaction kinetics like an uncompetitive inhibitor.

DynaBiS: A hierarchical sampling algorithm to identify flexible binding sites for large ligands and peptides.

Knowing the ligand or peptide binding site in proteins is highly important to guide drug discovery, but experimental elucidation of the binding site is difficult. Therefore, various computational approaches have been developed to identify potential binding sites in protein structures. However, protein and ligand flexibility are often neglected in these methods due to efficiency considerations despite the recognition that protein-ligand interactions can be strongly affected by mutual structural adaptations. This is particularly true if the binding site is unknown, as the screening will typically be performed based on an unbound protein structure. Herein we present DynaBiS, a hierarchical sampling algorithm to identify flexible binding sites for a target ligand with explicit consideration of protein and ligand flexibility, inspired by our previously presented flexible docking algorithm DynaDock. DynaBiS applies soft-core potentials between the ligand and the protein, thereby allowing a certain protein-ligand overlap resulting in efficient sampling of conformational adaptation effects. We evaluated DynaBiS and other commonly used binding site identification algorithms against a diverse evaluation set consisting of 26 proteins featuring peptide as well as small ligand binding sites. We show that DynaBiS outperforms the other evaluated methods for the identification of protein binding sites for large and highly flexible ligands such as peptides, both with a holo or apo structure used as input.

Turning the Actin Nucleating Compound Miuraenamide into Nucleation Inhibitors.

Natural compounds that either increase or decrease polymerization of actin into filaments have become indispensable tools for cell biology. However, to date, it was not possible to use them as therapeutics due to their overall cytotoxicity and their unfavorable pharmacokinetics. Furthermore, their synthesis is in general quite complicated. In an attempt to find simplified analogues of miuraenamide, an actin nucleating compound, we identified derivatives with a paradoxical inversion of the mode of action: instead of increased nucleation, they caused an inhibition. Using an extensive computational approach, we propose a binding mode and a mode of action for one of these derivatives. Based on our findings, it becomes feasible to tune actin-binding compounds to one or the other direction and to generate new synthetic actin binders with increased functional selectivity.

Functional analysis of peripheral and intratumoral neoantigen-specific TCRs identified in a patient with melanoma.

BACKGROUND: Neoantigens derived from somatic mutations correlate with therapeutic responses mediated by treatment with immune checkpoint inhibitors. Neoantigens are therefore highly attractive targets for the development of therapeutic approaches in personalized medicine, although many aspects of their quality and associated immune responses are not yet well understood. In a case study of metastatic malignant melanoma, we aimed to perform an in-depth characterization of neoantigens and respective T-cell responses in the context of immune checkpoint modulation. METHODS: Three neoantigens, which we identified either by immunopeptidomics or in silico prediction, were investigated using binding affinity analyses and structural simulations. We isolated seven T-cell receptors (TCRs) from the patient's immune repertoire recognizing these antigens. TCRs were compared in vitro by multiparametric analyses including functional avidity, multicytokine secretion, and cross-reactivity screenings. A xenograft mouse model served to study in vivo functionality of selected TCRs. We investigated the patient's TCR repertoire in blood and different tumor-related tissues over 3 years using TCR beta deep sequencing. RESULTS: Selected mutated peptide ligands with proven immunogenicity showed similar binding affinities to the human leukocyte antigen complex and comparable disparity to their wild-type counterparts in molecular dynamic simulations. Nevertheless, isolated TCRs recognizing these antigens demonstrated distinct patterns in functionality and frequency. TCRs with lower functional avidity showed at least equal antitumor immune responses in vivo. Moreover, they occurred at high frequencies and particularly demonstrated long-term persistence within tumor tissues, lymph nodes and various blood samples associated with a reduced activation pattern on primary in vitro stimulation. CONCLUSIONS: We performed a so far unique fine characterization of neoantigen-specific T-cell responses revealing defined reactivity patterns of neoantigen-specific TCRs. Our data highlight qualitative differences of these TCRs associated with function and longevity of respective T cells. Such features need to be considered for further optimization of neoantigen targeting including adoptive T-cell therapies using TCR-transgenic T cells.

Extended interaction networks with HCV protease NS3-4A substrates explain the lack of adaptive capability against protease inhibitors.

Inhibitors against the NS3-4A protease of hepatitis C virus (HCV) have proven to be useful drugs in the treatment of HCV infection. Although variants have been identified with mutations that confer resistance to these inhibitors, the mutations do not restore replicative fitness and no secondary mutations that rescue fitness have been found. To gain insight into the molecular mechanisms underlying the lack of fitness compensation, we screened known resistance mutations in infectious HCV cell culture with different genomic backgrounds. We observed that the Q41R mutation of NS3-4A efficiently rescues the replicative fitness in cell culture for virus variants containing mutations at NS3-Asp168 To understand how the Q41R mutation rescues activity, we performed protease activity assays complemented by molecular dynamics simulations, which showed that protease-peptide interactions far outside the targeted peptide cleavage sites mediate substrate recognition by NS3-4A and support protease cleavage kinetics. These interactions shed new light on the mechanisms by which NS3-4A cleaves its substrates, viral polyproteins and a prime cellular antiviral adaptor protein, the mitochondrial antiviral signaling protein MAVS. Peptide binding is mediated by an extended hydrogen-bond network in NS3-4A that was effectively optimized for protease-MAVS binding in Asp168 variants with rescued replicative fitness from NS3-Q41R. In the protease harboring NS3-Q41R, the N-terminal cleavage products of MAVS retained high affinity to the active site, rendering the protease susceptible for potential product inhibition. Our findings reveal delicately balanced protease-peptide interactions in viral replication and immune escape that likely restrict the protease adaptive capability and narrow the virus evolutionary space.

Systematic analysis of the binding behaviour of UHRF1 towards different methyl- and carboxylcytosine modification patterns at CpG dyads.

The multi-domain protein UHRF1 is essential for DNA methylation maintenance and binds DNA via a base-flipping mechanism with a preference for hemi-methylated CpG sites. We investigated its binding to hemi- and symmetrically modified DNA containing either 5-methylcytosine (mC), 5-hydroxymethylcytosine (hmC), 5-formylcytosine (fC), or 5-carboxylcytosine (caC). Our experimental results indicate that UHRF1 binds symmetrically carboxylated and hybrid methylated/carboxylated CpG dyads in addition to its previously reported substrates. Complementary molecular dynamics simulations provide a possible mechanistic explanation of how the protein could differentiate between modification patterns. First, we observe different local binding modes in the nucleotide binding pocket as well as the protein's NKR finger. Second, both DNA modification sites are coupled through key residues within the NKR finger, suggesting a communication pathway affecting protein-DNA binding for carboxylcytosine modifications. Our results suggest a possible additional function of the hemi-methylation reader UHRF1 through binding of carboxylated CpG sites. This opens the possibility of new biological roles of UHRF1 beyond DNA methylation maintenance and of oxidised methylcytosine derivates in epigenetic regulation.

A modular approach to the bisbenzylisoquinoline alkaloids tetrandrine and isotetrandrine.

An efficient racemic total synthesis of the bisbenzylisoquinoline alkaloids tetrandrine and isotetrandrine in four different routes is reported herein. Key steps of the synthesis include N-acyl Pictet-Spengler condensations to access the tetrahydroisoquinoline moieties, as well as copper-catalyzed Ullmann couplings for diaryl ether formation. Starting from commercially available building blocks tetrandrine and isotetrandrine are accessed in 12 steps. Depending on the sequence of the four central condensation steps, equimolar mixtures of both diastereomers or predominantly tetrandrine or its diastereomer isotetrandrine are obtained. Through computational analysis we were able to rationalize the differences in the observed diastereomeric specificities.

Acyldepsipeptide Probes Facilitate Specific Detection of Caseinolytic Protease†P Independent of Its Oligomeric and Activity State.

Caseinolytic protease P (ClpP) is a tetradecameric peptidase that assembles with chaperones such as ClpX to gain proteolytic activity. Acyldepsipeptides (ADEPs) are small-molecule mimics of ClpX that bind into hydrophobic pockets on the apical site of the complex, thereby activating ClpP. Detection of ClpP has so far been facilitated with active-site-directed probes which depend on the activity and oligomeric state of the complex. To expand the scope of ClpP labeling, we took a stepwise synthetic approach toward customized ADEP photoprobes. Structure-activity relationship studies with small fragments and ADEP derivatives paired with modeling studies revealed the design principles for suitable probe molecules. The derivatives were tested for activation of ClpP and subsequently applied in labeling studies of the wild-type peptidase as well as enzymes bearing mutations at the active site and an oligomerization sensor. Satisfyingly, the ADEP photoprobes provided a labeling readout of ClpP independent of its activity and oligomeric state.

Repurposing human kinase inhibitors to create an antibiotic active against drug-resistant Staphylococcus aureus, persisters and biofilms.

New drugs are desperately needed to combat methicillin-resistant Staphylococcus aureus (MRSA) infections. Here, we report screening commercial kinase inhibitors for antibacterial activity and found the anticancer drug sorafenib as major hit that effectively kills MRSA strains. Varying the key structural features led to the identification of a potent analogue, PK150, that showed antibacterial activity against several pathogenic strains at submicromolar concentrations. Furthermore, this antibiotic eliminated challenging persisters as well as established biofilms. PK150 holds promising therapeutic potential as it did not induce in vitro resistance, and shows oral bioavailability and in vivo efficacy. Analysis of the mode of action using chemical proteomics revealed several targets, which included interference with menaquinone biosynthesis by inhibiting demethylmenaquinone methyltransferase and the stimulation of protein secretion by altering the activity of signal peptidase IB. Reduced endogenous menaquinone levels along with enhanced levels of extracellular proteins of PK150-treated bacteria support this target hypothesis. The associated antibiotic effects, especially the lack of resistance development, probably stem from the compound's polypharmacology.

Key Features Relevant to Select Antigens and TCR From the MHC-Mismatched Repertoire to Treat Cancer.

Adoptive transfer of T cells transgenic for tumor-reactive T-cell receptors (TCR) is an attractive immunotherapeutic approach. However, clinical translation is so far limited due to challenges in the identification of suitable target antigens as well as TCRs that are concurrent safe and efficient. Definition of key characteristics relevant for effective and specific tumor rejection is essential to improve current TCR-based adoptive T-cell immunotherapies. We here characterized in-depth two TCRs derived from the human leukocyte antigen (HLA)-mismatched allogeneic repertoire targeting two different myeloperoxidase (MPO)-derived peptides presented by the same HLA-restriction element side by side comprising state of the art biochemical and cellular in vitro, in vivo, and in silico experiments. In vitro experiments reveal comparable functional avidities, off-rates, and cytotoxic activities for both TCRs. However, we observed differences especially with respect to cytokine secretion and cross-reactivity as well as in vivo activity. Biochemical and in silico analyses demonstrate different binding qualities of MPO-peptides to the HLA-complex determining TCR qualities. We conclude from our biochemical and in silico analyses of peptide-HLA-binding that rigid and high-affinity binding of peptides is one of the most important factors for isolation of TCRs with high specificity and tumor rejection capacity from the MHC-mismatched repertoire. Based on our results, we developed a workflow for selection of such TCRs with high potency and safety profile suitable for clinical translation.

Actin stabilizing compounds show specific biological effects due to their binding mode.

Actin binding compounds are widely used tools in cell biology. We compare the biological and biochemical effects of miuraenamide A and jasplakinolide, a structurally related prototypic actin stabilizer. Though both compounds have similar effects on cytoskeletal morphology and proliferation, they affect migration and transcription in a distinctive manner, as shown by a transcriptome approach in endothelial cells. In vitro, miuraenamide A acts as an actin nucleating, F-actin polymerizing and stabilizing compound, just like described for jasplakinolide. However, in contrast to jasplakinolide, miuraenamide A competes with cofilin, but not gelsolin or Arp2/3 for binding to F-actin. We propose a binding mode of miuraenamide A, explaining both its similarities and its differences to jasplakinolide. Molecular dynamics simulations suggest that the bromophenol group of miurenamide A interacts with residues Tyr133, Tyr143, and Phe352 of actin. This shifts the D-loop of the neighboring actin, creating tighter packing of the monomers, and occluding the binding site of cofilin. Since relatively small changes in the molecular structure give rise to this selectivity, actin binding compounds surprisingly are promising scaffolds for creating actin binders with specific functionality instead of just "stabilizers".

Thiazoline-Specific Amidohydrolase PurAH Is the Gatekeeper of Bottromycin Biosynthesis.

The ribosomally synthesized and post-translationally modified peptide (RiPP) bottromycin A2 possesses potent antimicrobial activity. Its biosynthesis involves the enzymatic formation of a macroamidine, a process previously suggested to require the concerted efforts of a YcaO enzyme (PurCD) and an amidohydrolase (PurAH) in vivo. In vitro, PurCD alone is sufficient to catalyze formation of the macroamidine, but the process is reversible. We set out to probe the role of PurAH in macroamidine formation in vitro. We demonstrate that PurAH is highly selective for macroamidine-containing precursor peptides and cleaves C-terminal of a thiazoline, thus removing the follower peptide. After follower cleavage, macroamidine formation is irreversible, indicating PurAH as the gatekeeper of bottromycin biosynthesis. The structure of PurAH suggests residues involved in catalysis, which were probed through mutagenesis.