Rational design of glycosaminoglycan binding cyclic peptides using cPEPmatch.

Cyclic peptides present a robust platform for drug design, offering high specificity and stability due to their conformationally constrained structures. In this study, we introduce an updated version of the Cyclic Peptide Matching program (cPEPmatch) tailored for the identification of cyclic peptides capable of mimicking protein-glycosaminoglycan (GAG) binding sites. We focused on engineering cyclic peptides to replicate the GAG-binding affinity of antithrombin III (ATIII), a protein that plays a crucial role in modulating anticoagulation through interaction with the GAG heparin. By integrating computational and experimental methods, we successfully identified a cyclic peptide binder with promising potential for future optimization. MD simulations and MM-GBSA calculations were used to assess binding efficacy, supplemented by umbrella sampling to approximate free energy landscapes. The binding specificity was further validated through NMR and ITC experiments. Our findings demonstrate that the computationally designed cyclic peptides effectively target GAGs, suggesting their potential as novel therapeutic agents. This study advances our understanding of peptide-GAG interactions and lays the groundwork for future development of cyclic peptide-based therapeutics.

Unraveling the mechanism of small molecule induced activation of Staphylococcus aureus signal peptidase IB.

Staphylococcus aureus signal peptidase IB (SpsB) is an essential enzyme for protein secretion. While inhibition of its activity by small molecules is a well-precedented mechanism to kill bacteria, the mode of activation is however less understood. We here investigate the activation mechanism of a recently introduced activator, the antibiotic compound PK150, and demonstrate by combined experimental and Molecular Dynamics (MD) simulation studies a unique principle of enzyme stimulation. Mass spectrometric studies with an affinity-based probe of PK150 unravel the binding site of PK150 in SpsB which is used as a starting point for MD simulations. Our model shows the localization of the molecule in an allosteric pocket next to the active site which shields the catalytic dyad from excess water that destabilizes the catalytic geometry. This mechanism is validated by the placement of mutations aligning the binding pocket of PK150. While the mutants retain turnover of the SpsB substrate, no stimulation of activity is observed upon PK150 addition. Overall, our study elucidates a previously little investigated mechanism of enzyme activation and serves as a starting point for the development of future enzyme activators.

Refinement of Docked Protein-Protein Complexes Using Repulsive Scaling Replica Exchange Simulations.

Accurate prediction and evaluation of protein-protein complex structures is of major importance to understand the cellular interactome. Predicted complex structures based on deep learning approaches or traditional docking methods require often structural refinement and rescoring for realistic evaluation. Standard molecular dynamics (MD) simulations are time-consuming and often do not structurally improve docking solutions. Better refinement can be achieved with our recently developed replica-exchange-based scheme employing different levels of repulsive biasing between proteins in each replica simulation (RS-REMD). The bias acts specifically on the intermolecular interactions based on an increase in effective pairwise van der Waals radii without changing interactions within each protein or with the solvent. It allows for an improvement of the predicted protein-protein complex structure and simultaneous realistic free energy scoring of protein-protein complexes. The setup of RS-REMD simulations is described in detail including the application on two examples (all necessary scripts and input files can be obtained from https://gitlab.com/TillCyrill/mmib ).

How arginine inhibits substrate-binding domain 2 elucidated using molecular dynamics simulations.

The substrate-binding domain 2 (SBD2) is an important part of the bacterial glutamine (GLN) transporter and mediates binding and delivery of GLN to the transporter translocation subunit. The SBD2 consists of two domains, D1 and D2, that bind GLN in the space between domains in a closed structure. In the absence of ligand, the SBD2 adopts an open conformation with larger space between domains. The GLN binding and closing are essential for the subsequent transport into the cell. Arginine (ARG) can also bind to SBD2 but does not induce closing and inhibits GLN transport. We use atomistic molecular dynamics (MD) simulations in explicit solvent to study ARG binding in the presence of the open SBD2 structure and observed reversible binding to the native GLN binding site with similar contacts but no transition to a closed SBD2 state. Absolute binding free energy simulations predict a considerable binding affinity of ARG and GLN to the binding site on the D1 domain. Free energy simulations to induce subsequent closing revealed a strong free energy penalty in case of ARG binding in contrast to GLN binding that favors the closed SBD2 state but still retains a free energy barrier for closing. The simulations allowed the identification of the molecular origin of the closing penalty in case of bound ARG and suggested a mutation of lysine at position 373 to alanine that strongly reduced the penalty and allowed closing even in the presence of bound ARG. The study offers an explanation of the molecular mechanism of how ARG competitively inhibits GLN from binding to SBD2 and from triggering the transition to a closed conformation. The proposed Lys373Ala mutation shows promise as a potential tool to validate whether a conformational mismatch between open SBD2 and the translocator is responsible for preventing ARG uptake to the cell.

Invasive Pulmonary Aspergillosis in Critically Ill Patients with Hantavirus Infection, Austria.

We investigated a cohort of 370 patients in Austria with hantavirus infections (7.8% ICU admission rate) and detected 2 cases (cumulative incidence 7%) of invasive pulmonary aspergillosis; 1 patient died. Hantavirus-associated pulmonary aspergillosis may complicate the course of critically ill patients who have hemorrhagic fever with renal syndrome.

Helical reorganization in the context of membrane protein folding: Insights from simulations with bacteriorhodopsin (BR) fragments.

Membrane protein folding is distinct from folding of soluble proteins. Conformational acquisition in major membrane protein subclasses can be delineated into insertion and folding processes. An exception to the "two stage" folding, later developed to "three stage" folding, is observed within the last two helices in bacteriorhodopsin (BR), a system that serves as a model membrane protein. We employ a reductionist approach to understand interplay of molecular factors underlying the apparent defiance. Leveraging available solution NMR structures, we construct, sample in silico, and analyze partially (PIn) and fully inserted (FIn) BR membrane states. The membrane lateral C-terminal helix (CH) in PIn is markedly prone to transient structural distortions over microsecond timescales; a disorder prone region (DPR) is thereby identified. While clear transmembrane propensities are not acquired, the distortions induce alterations in local membrane curvature and area per lipid. Importantly, energetic decompositions reveal that overall, the N-terminal helix (NH) is thermodynamically more stable in the PIn. Higher overall stability of the FIn arises from favorable interactions between the NH and the CH. Our results establish lack of spontaneous transition of the PIn to the FIn, and attributes their partitioning to barriers that exceed those accessible with thermal fluctuations. This work paves the way for further detailed studies aimed at determining the thermo-kinetic roles of the initial five helices, or complementary external factors, in complete helical folding and insertion in BR. We comment that complementing such efforts with the growing field of machine learning assisted energy landscape searches may offer unprecedented insights.

Structural Insights into Seeding Mechanisms of hIAPP Fibril Formation.

The deposition of islet amyloid polypeptide (hIAPP) fibrils is a hallmark of ß-cell death in type II diabetes. In this study, we employ state-of-the-art MAS solid-state spectroscopy to investigate the previously elusive N-terminal region of hIAPP fibrils, uncovering both rigidity and heterogeneity. Comparative analysis between wild-type hIAPP and a disulfide-deficient variant (hIAPPC2S,C7S) unveils shared fibril core structures yet strikingly distinct dynamics in the N-terminus. Specifically, the variant fibrils exhibit extended ß-strand conformations, facilitating surface nucleation. Moreover, our findings illuminate the pivotal roles of specific residues in modulating secondary nucleation rates. These results deepen our understanding of hIAPP fibril assembly and provide critical insights into the molecular mechanisms underpinning type II diabetes, holding promise for future therapeutic strategies.

Longitudinal analysis of PD-L1 expression in patients with relapsed NSCLC.

BACKGROUND: The use and approval of immune checkpoint inhibitors for the treatment of non-small cell lung cancer (NSCLC) depends on PD-L1 expression in the tumor tissue. Nevertheless, PD-L1 often fails to predict response to treatment. One possible explanation could be a change in PD-L1 expression during the course of the disease and the neglect of reassessment. The purpose of this study was a longitudinal analysis of PD-L1 expression in patients with relapsed NSCLC. METHODS: We retrospectively analyzed PD-L1 expression in patients with early-stage NSCLC and subsequent relapse in preoperative samples, matched surgical specimens and biopsy samples of disease recurrence. Ventana PD-L1 (SP263) immunohistochemistry assay was used for all samples. PD-L1 expression was scored based on clinically relevant groups (0%, 1%-49%, and ≥50%). The primary endpoint was the change in PD-L1 score group between preoperative samples, matched surgical specimens and relapsed tumor tissue. RESULTS: 395 consecutive patients with stages I-III NSCLC and 136 (34%) patients with a subsequent relapse were identified. For 87 patients at least two specimens for comparison of PD-L1 expression between early stage and relapsed disease were available. In 72 cases, a longitudinal analysis between preoperative biopsy, the surgically resected specimen and biopsy of disease recurrence was feasible. When comparing preoperative and matched surgical specimens, a treatment-relevant conversion of PD-L1 expression group was found in 25 patients (34.7%). Neoadjuvant treatment showed no significant effect on PD-L1 alteration (p=0.39). In 32 (36.8%) out of 87 cases, a change in PD-L1 group was observed when biopsies of disease relapse were compared with early-stage disease. Adjuvant treatment was not significantly associated with a change in PD-L1 expression (p=0.53). 39 patients (54.2%) showed at least 1 change into a different PD-L1 score group during the course of disease. 14 patients (19.4%) changed the PD-L1 score group twice, 5 (6.9%) of them being found in all different score groups. CONCLUSION: PD-L1 expression shows dynamic changes during the course of disease. There is an urgent need for consensus guidelines to define a PD-L1 testing strategy including time points of reassessment, the number of biopsies to be obtained and judgment of surgical specimens.

Structure-Based Protein Assembly Simulations Including Various Binding Sites and Conformations.

Many biological functions are mediated by large complexes formed by multiple proteins and other cellular macromolecules. Recent progress in experimental structure determination, as well as in integrative modeling and protein structure prediction using deep learning approaches, has resulted in a rapid increase in the number of solved multiprotein assemblies. However, the assembly process of large complexes from their components is much less well-studied. We introduce a rapid computational structure-based (SB) model, GoCa, that allows to follow the assembly process of large multiprotein complexes based on a known native structure. Beyond existing SB Go̅-type models, it distinguishes between intra- and intersubunit interactions, allowing us to include coupled folding and binding. It accounts automatically for the permutation of identical subunits in a complex and allows the definition of multiple minima (native) structures in the case of proteins that undergo global transitions during assembly. The model is successfully tested on several multiprotein complexes. The source code of the GoCa program including a tutorial is publicly available on Github: https://github.com/ZachariasLab/GoCa. We also provide a web source that allows users to quickly generate the necessary input files for a GoCa simulation: https://goca.t38webservices.nat.tum.de.

Efficient and accurate binding free energy calculation of Aß9-40 protofilament propagation.

Self-assembled aggregation of peptides and proteins into regular amyloid fibrils is associated with several neurodegenerative diseases. In case of Alzheimer's disease proteolytic cleavage products of the amyloid precursor protein form pathological amyloid-beta fibrils in a nucleation and propagation phase. The molecular details and thermodynamic driving forces of amyloid formation are not well understood, but are of high relevance for potential pharmacological interference. We used atomistic binding free energy simulations to calculate the free energy of protofilament propagation by an additional Aß9-40 peptide binding to the protofilament tip. It requires sampling of relevant conformational transitions which is challenging since the monomeric Aß9-40 peptide is intrinsically disordered. However, the convergence of umbrella simulations can be enhanced by applying additional restraining potentials on the axial, orientational and conformational degrees of freedom. The improved convergence leads to a much closer agreement with experimental binding free energy data compared to unrestrained umbrella sampling. Moreover, the restraining approach results in a separation of contributions to the total binding free energy. The calculated contributions indicate that the free energy change associated with the restriction of conformational freedom upon propagation makes a large opposing contribution of higher magnitude than the total binding free energy. Finally, optimization of the approach leads to further significant reduction of the computational demand which is crucial for systematic studies on mutations, denaturants and inhibitors in the fibril propagation step.

Liquid-Vapor Coexistence and Spontaneous Evaporation at Atmospheric Pressure of Common Rigid Three-Point Water Models in Molecular Simulations.

Molecular dynamics (MD) simulations are widely used to investigate molecular systems at atomic resolution including biomolecular structures, drug-receptor interactions, and novel materials. Frequently, MD simulations are performed in an aqueous solution with explicit models of water molecules. Commonly, such models are parameterized to reproduce the liquid phase of water under ambient conditions. However, often, simulations at significantly higher temperatures are also of interest. Hence, it is important to investigate the equilibrium of the liquid and vapor phases of molecular models of water at elevated temperatures. Here, we evaluate the behavior of 11 common rigid three-point water models over a wide range of temperatures. From liquid-vapor coexistence simulations, we estimated the critical points and studied the spontaneous evaporation of these water models. Moreover, we investigated the influence of the system size, choice of the pressure-coupling algorithm, and rate of heating on the process and compared them with the experimental data. We found that modern rigid three-point water models reproduce the critical point surprisingly well. Furthermore, we discovered that the critical temperature correlates with the quadrupole moment of the respective water model. This indicates that the spatial arrangement of the partial charges is important for reproducing the liquid-vapor phase transition. Our findings may guide the selection of water models for simulations conducted at high temperatures.

Expanding Broad Molecular Reflex Testing in Non-Small Cell Lung Cancer to Squamous Histology.

Due to the success story of biomarker-driven targeted therapy, most NSCLC guidelines agree that molecular reflex testing should be performed in all cases with non-squamous cell carcinoma (non-SCC). In contrast, testing recommendations for squamous cell carcinoma (SCC) vary considerably, specifically concerning the exclusion of patients of certain age or smoking status from molecular testing strategies. We performed a retrospective single-center study examining the value of molecular reflex testing in an unselected cohort of 316 consecutive lung SCC cases, tested by DNA- and RNA-based next-generation sequencing (NGS) at our academic institution between 2019 and 2023. Clinicopathological data from these cases were obtained from electronic medical records and correlated with sequencing results. In 21/316 (6.6%) cases, we detected an already established molecular target for an approved drug. Among these were seven cases with an EGFR mutation, seven with a KRAS G12C mutation, four with an ALK fusion, two with an EGFR fusion and one with a METex14 skipping event. All patients harboring a targetable alteration were >50 years of age and most of them had >15 pack-years, questioning restrictive molecular testing strategies. Based on our real-world data, we propose a reflex testing workflow using DNA- and RNA-based NGS that includes all newly diagnosed NSCLC cases, irrespective of histology, but also irrespective of age or smoking status.

Comprehensive Analysis of Coupled Proline Cis-Trans States in Bradykinin Using ωBP-REMD Simulations.

It is well-known that proline (Pro) cis-trans isomerization plays a decisive role in the folding and stabilization of proteins. The conformational coupling between isomerization states of different Pro residues in proteins during conformational adaptation processes is not well understood. In the present work, we investigate the coupled cis-trans isomerization of three Pro residues using bradykinin (BK), a partially unstructured nonapeptide hormone, as a model system. We use a recently developed enhanced-sampling molecular dynamics method (ω-bias potential replica exchange molecular dynamics; ωBP-REMD) that allows us to exhaustively sample all combinations of Pro isomer states and obtain converged probability densities of all eight state combinations within 885 ns ωBP-REMD simulations. In agreement with experiment, the all-trans state is seen to be the preferred isomer of zwitterionic aqueous BK. In about a third of its structures, this state presents the characteristic C-terminal ß-turn conformation; however, other isomer combinations also contribute significantly to the structural ensemble. Unbiased probabilities can be projected onto the peptide bond dihedral angles of the three Pro residues. This unveils the interdependence of the individual Pro isomerization states, i.e., a possible coupling of the different Pro isomers. The cis/trans equilibrium of a Pro residue can change by up to 2.5 kcal·mol-1, depending on the isomerization state of other Pro residues. For example, for Pro7, the simulations indicate that its cis state becomes favored compared to its trans state when Pro2 is switched from the trans state to the cis state. Our findings demonstrate the efficiency of the ωBP-REMD methodology and suggest that the coupling of Pro isomerization states may play an even more decisive role in larger folded proteins subject to more conformational restraints.

The Role of Force Fields and Water Models in Protein Folding and Unfolding Dynamics.

Protein folding is a fascinating, not fully understood phenomenon in biology. Molecular dynamics (MD) simulations are an invaluable tool to study conformational changes in atomistic detail, including folding and unfolding processes of proteins. However, the accuracy of the conformational ensembles derived from MD simulations inevitably relies on the quality of the underlying force field in combination with the respective water model. Here, we investigate protein folding, unfolding, and misfolding of fast-folding proteins by examining different force fields with their recommended water models, i.e., ff14SB with the TIP3P model and ff19SB with the OPC model. To this end, we generated long conventional MD simulations highlighting the perks and pitfalls of these setups. Using Markov state models, we defined kinetically independent conformational substates and emphasized their distinct characteristics, as well as their corresponding state probabilities. Surprisingly, we found substantial differences in thermodynamics and kinetics of protein folding, depending on the combination of the protein force field and water model, originating primarily from the different water models. These results emphasize the importance of carefully choosing the force field and the respective water model as they determine the accuracy of the observed dynamics of folding events. Thus, the findings support the hypothesis that the water model is at least equally important as the force field and hence needs to be considered in future studies investigating protein dynamics and folding in all areas of biophysics.

Characterizing ATP processing by the AAA+ protein p97 at the atomic level.

The human enzyme p97 regulates various cellular pathways by unfolding hundreds of protein substrates in an ATP-dependent manner, making it an essential component of protein homeostasis and an impactful pharmacological target. The hexameric complex undergoes substantial conformational changes throughout its catalytic cycle. Here we elucidate the molecular motions that occur at the active site in the temporal window immediately before and after ATP hydrolysis by merging cryo-EM, NMR spectroscopy and molecular dynamics simulations. p97 populates a metastable reaction intermediate, the ADP·Pi state, which is poised between hydrolysis and product release. Detailed snapshots reveal that the active site is finely tuned to trap and eventually discharge the cleaved phosphate. Signalling pathways originating at the active site coordinate the action of the hexamer subunits and couple hydrolysis with allosteric conformational changes. Our multidisciplinary approach enables a glimpse into the sophisticated spatial and temporal orchestration of ATP handling by a prototype AAA+ protein.

The mutual value of histopathology and ITS sequencing in the diagnosis of mucormycosis.

AIMS: Mucormycosis is a fast-progressing disease with a high mortality rate. The most important factor determining survival of patients is early and accurate diagnosis. Although histopathology often recognises invasive mould infections at first, histomorphology alone is insufficient in providing an accurate diagnosis. Unbiased molecular methods to detect and identify fungi are promising, yet their role in complementing routine histopathological workflows has not been studied sufficiently. METHODS AND RESULTS: We performed a retrospective single-centre study examining the clinical value of complementing histopathology with internal transcribed spacer (ITS) sequencing of fungal DNA in the routine diagnosis of mucormycosis. At our academic centre, we identified 14 consecutive mucormycosis cases diagnosed by histopathology and subsequent ITS sequencing. Using histomorphological examination, fungal hyphae could be detected in all cases; however, morphological features were unreliable regarding specifying the taxa. Subsequent ITS sequencing identified a remarkable phylogenetic diversity among Mucorales: the most common species was Rhizopus microsporus (six of 14; 42.9%), followed by Lichtheimia corymbifera (three of 14, 21.4%) and single detections of Rhizopus oryzae, Actinomucor elegans, Mucor circinelloides, Rhizomucor pusillus and Rhizomucor miehei (one of 14; 7.1%, respectively). In one case, we additionally detected Pneumocystis jirovecii in the same lung tissue specimen, suggesting a clinically relevant co-infection. Fungal culture was performed in 10 cases but yielded positive results in only two of 10 (20%), revealing its limited value in the diagnosis of mucormycosis. CONCLUSIONS: Our study demonstrates that a combination of histopathology and ITS sequencing is a practically feasible approach that outperforms fungal culture in detecting Mucorales in tissue-associated infections. Therefore, pathologists might adapt diagnostic workflows accordingly when mucormycosis is suspected.

Design of cyclic peptides as novel inhibitors of ICOS/ICOSL interaction.

Compared to small molecules and antibodies, cyclic peptides exhibit unique biochemical and therapeutic attributes in the realm of pharmaceutical applications. The interaction between the inducible costimulator (ICOS) and its ligand (ICOSL) plays a key role in T-cell differentiation and activation. ICOS/ICOSL inhibition results in a reduction in the promotion of immunosuppressive regulatory T cells (Tregs) in both hematologic malignancies and solid tumors. Herein, we implement the computational cPEPmatch approach to design the first examples of cyclic peptides that inhibit ICOS/ICOSL interaction. The top cyclic peptide from our approach possessed an IC50 value of 1.87 ± 0.15 µM as an ICOS/ICOSL inhibitor and exhibited excellent in vitro pharmacokinetic properties as a drug candidate. Our work will lay the groundwork for future endeavors in cancer drug discovery, with the goal of developing cyclic peptides that target the ICOS/ICOSL interaction.

Studying specificity in protein-glycosaminoglycan recognition with umbrella sampling.

In the past few decades, glycosaminoglycan (GAG) research has been crucial for gaining insights into various physiological, pathological, and therapeutic aspects mediated by the direct interactions between the GAG molecules and diverse proteins. The structural and functional heterogeneities of GAGs as well as their ability to bind specific proteins are determined by the sugar composition of the GAG, the size of the GAG chains, and the degree and pattern of sulfation. A deep understanding of the interactions in protein-GAG complexes is essential to explain their biological functions. In this study, the umbrella sampling (US) approach is used to pull away a GAG ligand from the binding site and then pull it back in. We analyze the binding interactions between GAGs of three types (heparin, desulfated heparan sulfate, and chondroitin sulfate) with three different proteins (basic fibroblast growth factor, acidic fibroblast growth factor, and cathepsin K). The main focus of our study was to evaluate whether the US approach is able to reproduce experimentally obtained structures, and how useful it can be for getting a deeper understanding of GAG properties, especially protein recognition specificity and multipose binding. We found that the binding free energy landscape in the proximity of the GAG native binding pose is complex and implies the co-existence of several binding poses. The sliding of a GAG chain along a protein surface could be a potential mechanism of GAG particular sequence recognition by proteins.