Fluorine teams up with water to restore inhibitor activity to mutant BPTI

Computational screening of the effects of mutations on protein-protein off-rates and dissociation mechanisms by τRAMD

The dissociation rate, or its reciprocal, the residence time (τ), is a crucial parameter for understanding the duration and biological impact of biomolecular interactions. Accurate prediction of τ is essential for understanding protein-protein interactions (PPIs) and identifying potential drug targets or modulators for tackling diseases. Conventional molecular dynamics simulation techniques are inherently constrained by their limited timescales, making it challenging to estimate residence times, which typically range from minutes to hours. Building upon its successful application in protein-small molecule systems, τ-Random Acceleration Molecular Dynamics (τRAMD) is here investigated for estimating dissociation rates of protein-protein complexes. τRAMD enables the observation of unbinding events on the nanosecond timescale, facilitating rapid and efficient computation of relative residence times. We tested this methodology for three protein-protein complexes and their extensive mutant datasets, achieving good agreement between computed and experimental data. By combining τRAMD with MD-IFP (Interaction Fingerprint) analysis, dissociation mechanisms were characterized and their sensitivity to mutations investigated, enabling the identification of molecular hotspots for selective modulation of dissociation kinetics. In conclusion, our findings underscore the versatility of τRAMD as a simple and computationally efficient approach for computing relative protein-protein dissociation rates and investigating dissociation mechanisms, thereby aiding the design of PPI modulators.

Prebound State Discovered in the Unbinding Pathway of Fluorinated Variants of the Trypsin-BPTI Complex Using Random Acceleration Molecular Dynamics Simulations

The serine protease trypsin forms a tightly bound inhibitor complex with the bovine pancreatic trypsin inhibitor (BPTI). The complex is stabilized by the P1 residue Lys15, which interacts with negatively charged amino acids at the bottom of the S1 pocket. Truncating the P1 residue of wildtype BPTI to α-aminobutyric acid (Abu) leaves a complex with moderate inhibitor strength, which is held in place by additional hydrogen bonds at the protein-protein interface. Fluorination of the Abu residue partially restores the inhibitor strength. The mechanism with which fluorination can restore the inhibitor strength is unknown, and accurate computational investigation requires knowledge of the binding and unbinding pathways. The preferred unbinding pathway is likely to be complex, as encounter states have been described before, and unrestrained umbrella sampling simulations of these complexes suggest additional energetic minima. Here, we use random acceleration molecular dynamics to find a new metastable state in the unbinding pathway of Abu-BPTI variants and wildtype BPTI from trypsin, which we call the prebound state. The prebound state and the fully bound state differ by a substantial shift in the position, a slight shift in the orientation of the BPTI variants, and changes in the interaction pattern. Particularly important is the breaking of three hydrogen bonds around Arg17. Fluorination of the P1 residue lowers the energy barrier of the transition between the fully bound state and prebound state and also lowers the energy minimum of the prebound state. While the effect of fluorination is in general difficult to quantify, here, it is in part caused by favorable stabilization of a hydrogen bond between Gln194 and Cys14. The interaction pattern of the prebound state offers insights into the inhibitory mechanism of BPTI and might add valuable information for the design of serine protease inhibitors.

Fluorinated Protein-Ligand Complexes: A Computational Perspective

Fluorine is an element renowned for its unique properties. Its powerful capability to modulate molecular properties makes it an attractive substituent for protein binding ligands; however, the rational design of fluorination can be challenging with effects on interactions and binding energies being difficult to predict. In this Perspective, we highlight how computational methods help us to understand the role of fluorine in protein-ligand binding with a focus on molecular simulation. We underline the importance of an accurate force field, present fluoride channels as a showcase for biomolecular interactions with fluorine, and discuss fluorine specific interactions like the ability to form hydrogen bonds and interactions with aryl groups. We put special emphasis on the disruption of water networks and entropic effects.

Non-Canonical Amino Acids in Analyses of Protease Structure and Function

All known organisms encode 20 canonical amino acids by base triplets in the genetic code. The cellular translational machinery produces proteins consisting mainly of these amino acids. Several hundred natural amino acids serve important functions in metabolism, as scaffold molecules, and in signal transduction. New side chains are generated mainly by post-translational modifications, while others have altered backbones, such as the β- or γ-amino acids, or they undergo stereochemical inversion, e.g., in the case of D-amino acids. In addition, the number of non-canonical amino acids has further increased by chemical syntheses. Since many of these non-canonical amino acids confer resistance to proteolytic degradation, they are potential protease inhibitors and tools for specificity profiling studies in substrate optimization and enzyme inhibition. Other applications include in vitro and in vivo studies of enzyme kinetics, molecular interactions and bioimaging, to name a few. Amino acids with bio-orthogonal labels are particularly attractive, enabling various cross-link and click reactions for structure-functional studies. Here, we cover the latest developments in protease research with non-canonical amino acids, which opens up a great potential, e.g., for novel prodrugs activated by proteases or for other pharmaceutical compounds, some of which have already reached the clinical trial stage.

Structural and functional properties of rBmTI-A: A Kunitz-BPTI serine protease inhibitor with therapeutical potential

Serine proteases are an important group of enzymes present in several organisms such as viruses, bacteria and eukaryotes involved in several physiological and pathological processes such as cancer, neurodegeneration, tissue inflammation and infections. Kunitz-type serine protease inhibitors have been studied as therapeutical targets with positive results in many of these diseases. rBmTI-A (recombinant B. microplus Trypsin Inhibitor A) is a Kunitz-BPTI type inhibitor based on the native protein BmTI-A. BmTI-A was extracted from tick larvae and presented inhibitory activity against trypsin, human plasma kallikrein (HuPK), human neutrophil elastase (HNE) and human plasmin. rBmTI-A presented the same inhibitory activities of the BmTI-A and its thermostability has already been demonstrated. In emphysema induced by porcine pancreatic elastase (PPE) and by cigarette smoke animal models, the treatment using rBmTI-A showed a protective effect against the development of pulmonary emphysema and prevented the increase of inflammatory cells. In chronic allergic animal model, rBmTI-A treatment resulted in attenuated bronchial hyperresponsiveness, inflammation, remodeling. These are important physiological results in emphysema and lung inflammatory animal models with rBmTI-A treatment confirming its therapeutical potential.

Water Network in the Binding Pocket of Fluorinated BPTI-Trypsin Complexes─Insights from Simulation and Experiment

Structural waters in the S1 binding pocket of β-trypsin are critical for the stabilization of the complex of β-trypsin with its inhibitor bovine pancreatic trypsin inhibitor (BPTI). The inhibitor strength of BPTI can be modulated by replacing the critical lysine residue at the P1 position by non-natural amino acids. We study BPTI variants in which the critical Lys15 in BPTI has been replaced by α-aminobutyric acid (Abu) and its fluorinated derivatives monofluoroethylglycine (MfeGly), difluoroethylglycine (DfeGly), and trifluoroethylglycine (TfeGly). We investigate the hypothesis that additional water molecules in the binding pocket can form specific noncovalent interactions with the fluorinated side chains and thereby act as an extension of the inhibitors. We report potentials of mean force (PMF) of the unbinding process for all four complexes and enzyme activity inhibition assays. Additionally, we report the protein crystal structure of the Lys15MfeGly-BPTI-β-trypsin complex (pdb: 7PH1). Both experimental and computational data show a stepwise increase in inhibitor strength with increasing fluorination of the Abu side chain. The PMF additionally shows a minimum for the encounter complex and an intermediate state just before the bound state. In the bound state, the computational analysis of the structure and dynamics of the water molecules in the S1 pocket shows a highly dynamic network of water molecules that does not indicate a rigidification or stabilizing trend in regard to energetic properties that could explain the increase in inhibitor strength. The analysis of the energy and the entropy of the water molecules in the S1 binding pocket using grid inhomogeneous solvation theory confirms this result. Overall, fluorination systematically changes the binding affinity, but the effect cannot be explained by a persistent water network in the binding pocket. Other effects, such as the hydrophobicity of fluorinated amino acids and the stability of the encounter complex as well as the additional minimum in the potential of mean force in the bound state, likely influence the affinity more directly.

The Nϵ-Rule for Serine, but Not Cysteine Catalytic Triads

Catalytic triads, composed of a serine or cysteine nucleophile, a histidine, and a third triad residue (typically Asp/Glu/Asn), are common in enzyme active sites and catalyze a wide variety of chemical reactions. Two types of triads can be distinguished: We refer to them as Nδ- or Nϵ-configured, depending on whether the histidine imidazole Nδ or Nϵ atom is close to the nucleophile Oγ/Sγ. In this study, we have analyzed triad configuration. In structural triads, the more stable Nδ-configuration predominates. For catalytic triads, the configuration depends on the nucleophile. When it is a cysteine residue, both configuration types occur, depending on the family. However, when the nucleophile is a serine residue, the less stable Nϵ-configuration is almost exclusively found. We posit that the energetically less favored conformation is selected for in serine triads to facilitate the otherwise difficult proton transfer from the nucleophile to the histidine residue.

Rational design of amphiphilic fluorinated peptides: evaluation of self-assembly properties and hydrogel formation

Advanced peptide-based nanomaterials composed of self-assembling peptides (SAPs) are of emerging interest in pharmaceutical and biomedical applications. The introduction of fluorine into peptides, in fact, offers unique opportunities to tune their biophysical properties and intermolecular interactions. In particular, the degree of fluorination plays a crucial role in peptide engineering as it can be used to control the characteristics of fluorine-specific interactions and, thus, peptide conformation and self-assembly. Here, we designed and explored a series of amphipathic peptides by incorporating the fluorinated amino acids (2S)-4-monofluoroethylglycine (MfeGly), (2S)-4,4-difluoroethylglycine (DfeGly) and (2S)-4,4,4-trifluoroethylglycine (TfeGly) as hydrophobic components. This approach enabled studying the impact of fluorination on secondary structure formation and peptide self-assembly on a systematic basis. We show that the interplay between polarity and hydrophobicity, both induced differentially by varying degrees of side chain fluorination, does affect peptide folding significantly. A greater degree of fluorination promotes peptide fibrillation and subsequent formation of physical hydrogels in physiological conditions. Molecular simulations revealed the key role played by electrostatically driven intra-chain and inter-chain contact pairs that are modulated by side chain fluorination and give insights into the different self-organization behaviour of selected peptides. Our study provides a systematic report about the distinct features of fluorinated oligomeric peptides with potential applications as peptide-based biomaterials.

Molecular Characterization of Kunitz-Type Protease Inhibitors from Blister Beetles (Coleoptera, Meloidae)

Protease inhibitors are widely studied since the unrestricted activity of proteases can cause extensive organ lesions. In particular, elastase activity is involved in the pathophysiology of acute lung injury, for example during SARS-CoV-2 infection, while serine proteases and thrombin-like proteases are involved in the development and/or pathology of the nervous system. Natural protease inhibitors have the advantage to be reversible and with few side effects and thus are increasingly considered as new drugs. Kunitz-type protease inhibitors (KTPIs), reported in the venom of various organisms, such as wasps, spiders, scorpions, and snakes, have been studied for their potent anticoagulant activity and widespread protease inhibitor activity. Putative KTPI anticoagulants have been identified in transcriptomic resources obtained for two blister beetle species, Lydus trimaculatus and Mylabris variabilis. The KTPIs of L. trimaculatus and M. variabilis were characterized by combined transcriptomic and bioinformatics methodologies. The full-length mRNA sequences were divided on the base of the sequence of the active sites of the putative proteins. In silico protein structure analyses of each group of translational products show the biochemical features of the active sites and the potential protease targets. Validation of these genes is the first step for considering these molecules as new drugs for use in medicine.

Fluorine-induced polarity increases inhibitory activity of BPTI towards chymotrypsin

Substituting the P1 position in bovine pancreatic trypsin inhibitor (BPTI) is known to heavily influence its inhibitory activity towards serine proteases. Side-chain fluorinated aliphatic amino acids have been shown to alter numerous properties of peptides and proteins and thus are of interest in the context of BPTI. In our study, we systematically investigated the site-specific incorporation of non-canonical amino acids into BPTI by microwave-assisted solid-phase peptide synthesis (SPPS). Inhibitor activity of the variants was tested towards the serine protease α-chymotrypsin. We observed enhanced inhibition of two fluorinated BPTIs compared to wild type and hydrocarbon variants. To further investigate the complexes, we performed X-ray structure analysis. Our findings underline the power fluorine offers as a tool in protein engineering to beneficially alter the effects on phenomena as protein-protein interactions.

Identification of 13 Guanidinobenzoyl- or Aminidinobenzoyl-Containing Drugs to Potentially Inhibit TMPRSS2 for COVID-19 Treatment

Positively charged groups that mimic arginine or lysine in a natural substrate of trypsin are necessary for drugs to inhibit the trypsin-like serine protease TMPRSS2 that is involved in the viral entry and spread of coronaviruses, including SARS-CoV-2. Based on this assumption, we identified a set of 13 approved or clinically investigational drugs with positively charged guanidinobenzoyl and/or aminidinobenzoyl groups, including the experimentally verified TMPRSS2 inhibitors Camostat and Nafamostat. Molecular docking using the C-I-TASSER-predicted TMPRSS2 catalytic domain model suggested that the guanidinobenzoyl or aminidinobenzoyl group in all the drugs could form putative salt bridge interactions with the side-chain carboxyl group of Asp435 located in the S1 pocket of TMPRSS2. Molecular dynamics simulations further revealed the high stability of the putative salt bridge interactions over long-time (100 ns) simulations. The molecular mechanics/generalized Born surface area-binding free energy assessment and per-residue energy decomposition analysis also supported the strong binding interactions between TMPRSS2 and the proposed drugs. These results suggest that the proposed compounds, in addition to Camostat and Nafamostat, could be effective TMPRSS2 inhibitors for COVID-19 treatment by occupying the S1 pocket with the hallmark positively charged groups.

Molecular mechanism of inhibiting the SARS-CoV-2 cell entry facilitator TMPRSS2 with camostat and nafamostat

The entry of the coronavirus SARS-CoV-2 into human lung cells can be inhibited by the approved drugs camostat and nafamostat. Here we elucidate the molecular mechanism of these drugs by combining experiments and simulations. In vitro assays confirm that both drugs inhibit the human protein TMPRSS2, a SARS-Cov-2 spike protein activator. As no experimental structure is available, we provide a model of the TMPRSS2 equilibrium structure and its fluctuations by relaxing an initial homology structure with extensive 330 microseconds of all-atom molecular dynamics (MD) and Markov modeling. Through Markov modeling, we describe the binding process of both drugs and a metabolic product of camostat (GBPA) to TMPRSS2, reaching a Michaelis complex (MC) state, which precedes the formation of a long-lived covalent inhibitory state. We find that nafamostat has a higher MC population than camostat and GBPA, suggesting that nafamostat is more readily available to form the stable covalent enzyme–substrate intermediate, effectively explaining its high potency. This model is backed by our in vitro experiments and consistent with previous virus cell entry assays. Our TMPRSS2–drug structures are made public to guide the design of more potent and specific inhibitors.

The authors unravel the molecular action principle of nafamostat and camostat, two potential COVID-19 drugs targeting the human protein TMPRSS2.

Quantifying how step-wise fluorination tunes local solute hydrophobicity, hydration shell thermodynamics and the quantum mechanical contributions of solute–water interactions

Locally tuning solute–water interactions with fluorination.

Polyarylene polyimides with hydrocarbon and semi-fluorinated backbones: synthesis, characterization, and properties

A highly phenylated diamine was synthesized for the microwave preparation of six new polyimides.

Overview of Recent Strategic Advances in Medicinal Chemistry

Introducing novel strategies, concepts, and technologies that speed up drug discovery and the drug development cycle is of great importance both in the highly competitive pharmaceutical industry as well as in academia. This Perspective aims to present a "big-picture" overview of recent strategic innovations in medicinal chemistry and drug discovery.

Molecular insight into chymotrypsin inhibitor 2 resisting proteolytic degradation

Chymotrypsin inhibitor 2 (CI2) is a special serine protease inhibitor which can resist hydrolysis for several days with a rapid equilibrium between the Michaelis complex and acyl-enzyme intermediate. The energies and conformational changes for subtilisin-catalyzed proteolysis of CI2 were examined in this paper for the first time by employing pseudo bond ab initio QM/MM MD simulations. In the acylation reaction, a low-barrier hydrogen bond between His64 and Asp32 in the transition state together with the lack of covalent backbone constraints makes the peptide bonds of CI2 break more easily than in other serine protease inhibitors. After acyl-enzyme formation, molecular dynamics simulations showed that the access of hydrolytic water to the active site requires partial dissociation of the leaving group. However, retention of the leaving group mainly by the intra- and inter-molecular H-bonding networks hinders the access of water and retards the deacylation reaction. Instead of the dissociation constant of inhibitors, we suggest employing the free energy at the acyl-enzyme state to predict the relative hydrolysis rates of CI2 mutants, which are testified by the experimental relative hydrolysis rates.

Unexpected trends in the hydrophobicity of fluorinated amino acids reflect competing changes in polarity and conformation

The hydration free energy of fluorinated amino acids is calculated with molecular simulations and explained with an analytical model.

Synthesis and physical chemical properties of 2-amino-4-(trifluoromethoxy)butanoic acid - a CF3O-containing analogue of natural lipophilic amino acids

2-Amino-2-(trifluoromethoxy)butanoic acid (O-trifluoromethyl homoserine) was synthesized as a racemate and in both enantiomeric forms. The measured pK_a and log D values establish the compound as a promising analogue of natural aliphatic amino acids.

Position-dependent impact of hexafluoroleucine and trifluoroisoleucine on protease digestion

Rapid digestion by proteases limits the application of peptides as therapeutics. One strategy to increase the proteolytic stability of peptides is the modification with fluorinated amino acids. This study presents a systematic investigation of the effects of fluorinated leucine and isoleucine derivatives on the proteolytic stability of a peptide that was designed to comprise substrate specificities of different proteases. Therefore, leucine, isoleucine, and their side-chain fluorinated variants were site-specifically incorporated at different positions of this peptide resulting in a library of 13 distinct peptides. The stability of these peptides towards proteolysis by α-chymotrypsin, pepsin, proteinase K, and elastase was studied, and this process was followed by an FL-RP-HPLC assay in combination with mass spectrometry. In a few cases, we observed an exceptional increase in proteolytic stability upon introduction of the fluorine substituents. The opposite phenomenon was observed in other cases, and this may be explained by specific interactions of fluorinated residues with the respective enzyme binding sites. Noteworthy is that 5,5,5-trifluoroisoleucine is able to significantly protect peptides from proteolysis by all enzymes included in this study when positioned N-terminal to the cleavage site. These results provide valuable information for the application of fluorinated amino acids in the design of proteolytically stable peptide-based pharmaceuticals.

Deciphering the Fluorine Code-The Many Hats Fluorine Wears in a Protein Environment

Deciphering the fluorine code is how we describe not only the focus of this Account, but also the systematic approach to studying the impact of fluorine's incorporation on the properties of peptides and proteins used by our groups and others. The introduction of fluorine has been shown to impart favorable, but seldom predictable, properties to peptides and proteins, but up until about two decades ago the outcomes of fluorine modification of peptides and proteins were largely left to chance. Driven by the motivation to extend the application of the unique properties of the element fluorine from medicinal and agro chemistry to peptide and protein engineering we have established extensive research programs that enable the systematic investigation of effects that accompany the introduction of fluorine into this class of biopolymers. The introduction of fluorine into amino acids offers a universe of options for modifications with regard to number and position of fluorine substituents in the amino acid side chain. Moreover, it is important to emphasize that the consequences of incorporating the C-F bond into a biopolymer can be attributed to two distinct yet related phenomena: (i) the fluorine substituent can directly engage in intermolecular interactions with its environment and/or (ii) the other functional groups present in the molecule can be influenced by the electron withdrawing nature of this element (intramolecular) and in turn interact differently with their immediate environment (intermolecular). Based on our studies, we have shown that a change in number and/or position of as subtle as one single fluorine substituent has the power to considerably modify key properties of amino acids such as hydrophobicity, polarity, and secondary structure propensity. These properties are crucial factors in peptide and protein engineering, and thus, fluorinated amino acids can be applied to fine-tune properties such as protein folding, proteolytic stability, and protein-protein interactions provided we understand and become able to predict the outcome of a fluorine substitution in this context. With this Account, we attempt to analyze information we gained from our recent projects on how the nature of the fluorine atom and C-F bond influence four key properties of peptides and proteins: peptide folding, protein-protein interactions, ribosomal translation, and protease stability. These results impressively show why the introduction of fluorine creates a new class of amino acids with a repertoire of functionalities that is unique to the world of proteins and in some cases orthogonal to the set of canonical and natural amino acids. Our concluding statements aim to offer a few conserved design principles that have emerged from systematic studies over the last two decades; in this way, we hope to advance the field of peptide and protein engineering based on the judicious introduction of fluorinated building blocks.

Discovery and Investigation of Natural Editing Function against Artificial Amino Acids in Protein Translation

Fluorine being not substantially present in the chemistry of living beings is an attractive element in tailoring novel chemical, biophysical, and pharmacokinetic properties of peptides and proteins. The hallmark of ribosome-mediated artificial amino acid incorporation into peptides and proteins is a broad substrate tolerance, which is assumed to rely on the absence of evolutionary pressure for efficient editing of artificial amino acids. We used the well-characterized editing proficient isoleucyl-tRNA synthetase (IleRS) from Escherichia coli to investigate the crosstalk of aminoacylation and editing activities against fluorinated amino acids. We show that translation of trifluoroethylglycine (TfeGly) into proteins is prevented by hydrolysis of TfeGly-tRNAIle in the IleRS post-transfer editing domain. The remarkable observation is that dissociation of TfeGly-tRNAIle from IleRS is significantly slowed down. This finding is in sharp contrast to natural editing reactions by tRNA synthetases wherein fast editing rates for the noncognate substrates are essential to outcompete fast aa-tRNA dissociation rates. Using a post-transfer editing deficient mutant of IleRS (IleRSAla10), we were able to achieve ribosomal incorporation of TfeGly in vivo. Our work expands the knowledge of ribosome-mediated artificial amino acid translation with detailed analysis of natural editing function against an artificial amino acid providing an impulse for further systematic investigations and engineering of the translation and editing of unusual amino acids.

Targeted Fluoro Positioning for the Discovery of a Potent and Highly Selective Matrix Metalloproteinase Inhibitor

The incorporation of fluorine atoms into functional molecules is of wide interest in synthetic organic chemistry as well as cognate disciplines. In particular, in medicinal chemistry, there is a strong desire to positively influence the physicochemical molecular properties of drug compounds by introducing fluorine into biologically active molecules. Here, we present targeted fluoro positioning as the key design principle of converting a weak matrix metalloproteinase‐13 (MMP‐13) inhibitor into a very potent (IC₅₀=6 nm) and highly selective (selectivity factors of >1000 over MMP‐1, 2, 3, 7, 8, 9, 10, 12, 14) inhibitor with excellent plasma and microsomal stability, and no binding to the hERG channel (hERG: human ether‐a‐go‐go related gene).

An Acrobatic Substrate Metamorphosis Reveals a Requirement for Substrate Conformational Dynamics in Trypsin Proteolysis

The molecular basis of enzyme catalytic power and specificity derives from dynamic interactions between enzyme and substrate during catalysis. Although considerable effort has been devoted to understanding how conformational dynamics within enzymes affect catalysis, the role of conformational dynamics within protein substrates has not been addressed. Here, we examine the importance of substrate dynamics in the cleavage of Kunitz-bovine pancreatic trypsin inhibitor protease inhibitors by mesotrypsin, finding that the varied conformational dynamics of structurally similar substrates can profoundly impact the rate of catalysis. A 1.4-Å crystal structure of a mesotrypsin-product complex formed with a rapidly cleaved substrate reveals a dramatic conformational change in the substrate upon proteolysis. By using long all-atom molecular dynamics simulations of acyl-enzyme intermediates with proteolysis rates spanning 3 orders of magnitude, we identify global and local dynamic features of substrates on the nanosecond-microsecond time scale that correlate with enzymatic rates and explain differential susceptibility to proteolysis. By integrating multiple enhanced sampling methods for molecular dynamics, we model a viable conformational pathway between substrate-like and product-like states, linking substrate dynamics on the nanosecond-microsecond time scale with large collective substrate motions on the much slower time scale of catalysis. Our findings implicate substrate flexibility as a critical determinant of catalysis.