Reconciling the lock-and-key and dynamic views of canonical serine protease inhibitor action

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.

Structural Basis of Selective Human Indoleamine-2,3-dioxygenase 1 (hIDO1) Inhibition

hIDO1 is a heme-dioxygenase overexpressed in the tumor microenvironment and is implicated in the survival of cancer cells. Metabolism of tryptophan to N-formyl-kynurenine by hIDO1 leads to immune suppression to result in cancer cell immune escape. In this article, we discuss the discovery of selective hIDO1 inhibitors for therapeutic intervention that have been promoted to clinical trials and for which crystallographic structural information is available for the respective inhibitor-enzyme complex. The structural insights are based on the complex crystal structures and the relative biological data profiles. The structural basis of selective hIDO1 inhibition, as discussed herein, opens new avenues to the discovery of novel inhibitors with improved activity profiles, selectivity, and distinct structure frameworks.

Bowman-Birk Inhibitors: Insights into Family of Multifunctional Proteins and Peptides with Potential Therapeutical Applications

Bowman-Birk inhibitors (BBIs) are found primarily in seeds of legumes and in cereal grains. These canonical inhibitors share a highly conserved nine-amino acids binding loop motif CTP1SXPPXC (where P1 is the inhibitory active site, while X stands for various amino acids). They are natural controllers of plants' endogenous proteases, but they are also inhibitors of exogenous proteases present in microbials and insects. They are considered as plants' protective agents, as their elevated levels are observed during injury, presence of pathogens, or abiotic stress, i.a. Similar properties are observed for peptides isolated from amphibians' skin containing 11-amino acids disulfide-bridged loop CWTP1SXPPXPC. They are classified as Bowman-Birk like trypsin inhibitors (BBLTIs). These inhibitors are resistant to proteolysis and not toxic, and they are reported to be beneficial in the treatment of various pathological states. In this review, we summarize up-to-date research results regarding BBIs' and BBLTIs' inhibitory activity, immunomodulatory and anti-inflammatory activity, antimicrobial and insecticidal strength, as well as chemopreventive properties.

Directed Evolution of Canonical Loops and Their Swapping between Unrelated Serine Proteinase Inhibitors Disprove the Interscaffolding Additivity Model

Reversible serine proteinase inhibitors comprise 18 unrelated families. Each family has a distinct representative structure but contains a surface loop that adopts the same, canonical conformation in the enzyme-inhibitor complex. The Laskowski mechanism universally applies for the action of all canonical inhibitors independent of their scaffold, but it has two nontrivial extrapolations. Intrascaffolding additivity states that all enzyme-contacting loop residues act independently of each other, while interscaffolding additivity claims that these residues act independently of the scaffold. These theories have great importance for engineering proteinase inhibitors but have not been comprehensively challenged. Therefore, we tested the interscaffolding additivity theory by hard-randomizing all enzyme-contacting canonical loop positions of a Kazal- and a Pacifastin-scaffold inhibitor, displaying the variants on M13 phage, and selecting the libraries on trypsin and chymotrypsin. Directed evolution delivered different patterns on both scaffolds against both enzymes, which contradicts interscaffolding additivity. To quantitatively assess the extent of non-additivity, we measured the affinities of the optimal binding loop variants and their binding loop-swapped versions. While optimal variants have picomolar affinities, swapping the evolved loops results in up to 200,000-fold affinity loss. To decipher the underlying causes, we characterized the stability, overall structure and dynamics of the inhibitors with differential scanning calorimetry, circular dichroism and NMR spectroscopy and molecular dynamic simulations. These studies revealed that the foreign loop destabilizes the lower-stability Pacifastin scaffold, while the higher-stability Kazal scaffold distorts the foreign loop. Our findings disprove interscaffolding additivity and show that loop and scaffold form one integrated unit that needs to be coevolved to provide high-affinity inhibition.

Fine-tuning the extent and dynamics of binding cleft opening as a potential general regulatory mechanism in parvulin-type peptidyl prolyl isomerases

Parvulins or rotamases form a distinct group within peptidyl prolyl cis-trans isomerases. Their exact mode of action as well as the role of conserved residues in the family are still not unambiguously resolved. Using backbone S² order parameters and NOEs as restraints, we have generated dynamic structural ensembles of three distinct parvulins, SaPrsA, TbPin1 and CsPinA. The resulting ensembles are in good agreement with the experimental data but reveal important differences between the three enzymes. The largest difference can be attributed to the extent of the opening of the substrate binding cleft, along which motional mode the three molecules occupy distinct regions. Comparison with a wide range of other available parvulin structures highlights structural divergence along the bottom of the binding cleft acting as a hinge during the opening-closing motion. In the prototype WW-domain containing parvulin, Pin1, this region is also important in forming contacts with the WW domain known to modulate enzymatic activity of the catalytic domain. We hypothesize that modulation of the extent and dynamics of the identified ‘breathing motion’ might be one of the factors responsible for functional differences in the distinct parvulin subfamilies.

Protein Ensembles: How Does Nature Harness Thermodynamic Fluctuations for Life? The Diverse Functional Roles of Conformational Ensembles in the Cell

All soluble proteins populate conformational ensembles that together constitute the native state. Their fluctuations in water are intrinsic thermodynamic phenomena, and the distributions of the states on the energy landscape are determined by statistical thermodynamics; however, they are optimized to perform their biological functions. In this review we briefly describe advances in free energy landscape studies of protein conformational ensembles. Experimental (nuclear magnetic resonance, small-angle X-ray scattering, single-molecule spectroscopy, and cryo-electron microscopy) and computational (replica-exchange molecular dynamics, metadynamics, and Markov state models) approaches have made great progress in recent years. These address the challenging characterization of the highly flexible and heterogeneous protein ensembles. We focus on structural aspects of protein conformational distributions, from collective motions of single- and multi-domain proteins, intrinsically disordered proteins, to multiprotein complexes. Importantly, we highlight recent studies that illustrate functional adjustment of protein conformational ensembles in the crowded cellular environment. We center on the role of the ensemble in recognition of small- and macro-molecules (protein and RNA/DNA) and emphasize emerging concepts of protein dynamics in enzyme catalysis. Overall, protein ensembles link fundamental physicochemical principles and protein behavior and the cellular network and its regulation.

Molecular dynamics simulations of the free and inhibitor-bound cruzain systems in aqueous solvent: insights on the inhibition mechanism in acidic pH

The major cysteine protease of Trypanosoma cruzi, cruzain (CRZ), has been described as a therapeutic target for Chagas' disease, which affects millions of people worldwide. Thus, a series of CRZ inhibitors has been studied, including a new competitive inhibitor, Nequimed176 (NEQ176). Nevertheless, the structural and dynamic basis for CRZ inhibition remains unclear. Hoping to contribute to this ever-growing understanding of timescale dynamics in the CRZ inhibition mechanism, we have performed the first study using 100 ns of molecular dynamics (MD) simulations of two CRZ systems in an aqueous solvent under pH 5.5: CRZ in the apo form (ligand free) and CRZ complexed to NEQ176. According to the MD simulations, the enzyme adopts an open conformation in the apo form and a closed conformation in the NEQ176-CRZ complex. We also suggest that this closed conformation is related to the hydrogen-bonding interactions between NEQ176 and CRZ, which occurs through key residues, mainly Gly66, Met68, Asn69, and Leu160. In addition, the cross-correlation analysis shows evidence of the correlated motions among Ala110-Asp140, Leu160-Gly189, and Glu190-Gly215 subdomains, as well as, the movements related to Ala1-Thr59 and Asp60-Pro90 regions seem to be crucial for CRZ activity.

Simulated unbound structures for benchmarking of protein docking in the DOCKGROUND resource

Background

Proteins play an important role in biological processes in living organisms. Many protein functions are based on interaction with other proteins. The structural information is important for adequate description of these interactions. Sets of protein structures determined in both bound and unbound states are essential for benchmarking of the docking procedures. However, the number of such proteins in PDB is relatively small. A radical expansion of such sets is possible if the unbound structures are computationally simulated.

Results

The Dockground public resource provides data to improve our understanding of protein–protein interactions and to assist in the development of better tools for structural modeling of protein complexes, such as docking algorithms and scoring functions. A large set of simulated unbound protein structures was generated from the bound structures. The modeling protocol was based on 1 ns Langevin dynamics simulation. The simulated structures were validated on the ensemble of experimentally determined unbound and bound structures. The set is intended for large scale benchmarking of docking algorithms and scoring functions.

Conclusions

A radical expansion of the unbound protein docking benchmark set was achieved by simulating the unbound structures. The simulated unbound structures were selected according to criteria from systematic comparison of experimentally determined bound and unbound structures. The set is publicly available at http://dockground.compbio.ku.edu.

Engineered protease inhibitors based on sunflower trypsin inhibitor-1 (SFTI-1) provide insights into the role of sequence and conformation in Laskowski mechanism inhibition

Laskowski inhibitors regulate serine proteases by an intriguing mode of action that involves deceiving the protease into synthesizing a peptide bond. Studies exploring naturally occurring Laskowski inhibitors have uncovered several structural features that convey the inhibitor's resistance to hydrolysis and exceptional binding affinity. However, in the context of Laskowski inhibitor engineering, the way that various modifications intended to fine-tune an inhibitor's potency and selectivity impact on its association and dissociation rates remains unclear. This information is important as Laskowski inhibitors are becoming increasingly used as design templates to develop new protease inhibitors for pharmaceutical applications. In this study, we used the cyclic peptide, sunflower trypsin inhibitor-1 (SFTI-1), as a model system to explore how the inhibitor's sequence and structure relate to its binding kinetics and function. Using enzyme assays, MD simulations and NMR spectroscopy to study SFTI variants with diverse sequence and backbone modifications, we show that the geometry of the binding loop mainly influences the inhibitor's potency by modulating the association rate, such that variants lacking a favourable conformation show dramatic losses in activity. Additionally, we show that the inhibitor's sequence (including both the binding loop and its scaffolding) influences its potency and selectivity by modulating both the association and the dissociation rates. These findings provide new insights into protease inhibitor function and design that we apply by engineering novel inhibitors for classical serine proteases, trypsin and chymotrypsin and two kallikrein-related peptidases (KLK5 and KLK14) that are implicated in various cancers and skin diseases.

"Invisible" conformers of an antifungal disulfide protein revealed by constrained cold and heat unfolding, CEST-NMR experiments, and molecular dynamics calculations

Transition between conformational states in proteins is being recognized as a possible key factor of function. In support of this, hidden dynamic NMR structures were detected in several cases up to populations of a few percent. Here, we show by two- and three-state analysis of thermal unfolding, that the population of hidden states may weight 20-40 % at 298 K in a disulfide-rich protein. In addition, sensitive (15) N-CEST NMR experiments identified a low populated (0.15 %) state that was in slow exchange with the folded PAF protein. Remarkably, other techniques failed to identify the rest of the NMR "dark matter". Comparison of the temperature dependence of chemical shifts from experiments and molecular dynamics calculations suggests that hidden conformers of PAF differ in the loop and terminal regions and are most similar in the evolutionary conserved core. Our observations point to the existence of a complex conformational landscape with multiple conformational states in dynamic equilibrium, with diverse exchange rates presumably responsible for the completely hidden nature of a considerable fraction.

Defense response in non-genomic model species: methyl jasmonate exposure reveals the passion fruit leaves' ability to assemble a cocktail of functionally diversified Kunitz-type trypsin inhibitors and recruit two of them against papain

Multiplicity of protease inhibitors induced by predators may increase the understanding of a plant's intelligent behavior toward environmental challenges. Information about defense mechanisms of non-genomic model plant passion fruit (Passiflora edulis Sims) in response to predator attack is still limited. Here, via biochemical approaches, we showed its flexibility to build-up a broad repertoire of potent Kunitz-type trypsin inhibitors (KTIs) in response to methyl jasmonate. Seven inhibitors (20-25 kDa) were purified from exposed leaves by chromatographic techniques. Interestingly, the KTIs possessed truncated Kunitz motif in their N-terminus and some of them also presented non-consensus residues. Gelatin-Native-PAGE established multiple isoforms for each inhibitor. Significant differences regarding inhibitors' activity toward trypsin and chymotrypsin were observed, indicating functional polymorphism. Despite its rarity, two of them also inhibited papain, and such bifunctionality suggests a recruiting process onto another mechanistic class of target protease (cysteine-type). All inhibitors acted strongly on midgut proteases from sugarcane borer, Diatraea saccharalis (a lepidopteran insect) while in vivo assays supported their insecticide properties. Moreover, the bifunctional inhibitors displayed activity toward midgut proteases from cowpea weevil, Callosobruchus maculatus (a coleopteran insect). Unexpectedly, all inhibitors were highly effective against midgut proteases from Aedes aegypti a dipteran insect (vector of neglected tropical diseases) opening new avenues for plant-derived PIs for vector control-oriented research. Our results reflect the KTIs' complexities in passion fruit which could be wisely exploited by influencing plant defense conditions. Therefore, the potential of passion fruit as source of bioactive compounds with diversified biotechnological application was strengthened.

Specificity of a protein-protein interface: local dynamics direct substrate recognition of effector caspases

Proteases are prototypes of multispecific protein–protein interfaces. Proteases recognize and cleave protein and peptide substrates at a well-defined position in a substrate binding groove and a plethora of experimental techniques provide insights into their substrate recognition. We investigate the caspase family of cysteine proteases playing a key role in programmed cell death and inflammation, turning caspases into interesting drug targets. Specific ligand binding to one particular caspase is difficult to achieve, as substrate specificities of caspase isoforms are highly similar. In an effort to rationalize substrate specificity of two closely related caspases, we investigate the substrate promiscuity of the effector Caspases 3 and 7 by data mining (cleavage entropy) and by molecular dynamics simulations. We find a strong correlation between binding site rigidity and substrate readout for individual caspase subpockets explaining more stringent substrate readout of Caspase 7 via its narrower conformational space. Caspase 3 subpockets S3 and S4 show elevated local flexibility explaining the more unspecific substrate readout of that isoform in comparison to Caspase 7. We show by in silico exchange mutations in the S3 pocket of the proteases that a proline residue in Caspase 7 contributes to the narrowed conformational space of the binding site. These findings explain the substrate specificities of caspases via a mechanism of conformational selection and highlight the crucial importance of binding site local dynamics in substrate recognition of proteases. Proteins 2014; 82:546–555.

Ensemble-based interpretations of NMR structural data to describe protein internal dynamics

NMR spectroscopy is the leading technique to characterize protein internal dynamics at the atomic level and on multiple time scales. However, the structural interpretation of the observables obtained by various measurements is not always straightforward and in many cases dynamics-related parameters are only used to “decorate” static structural models without offering explicit description of conformational heterogeneity. To overcome such limitations, several computational techniques have been developed to generate ensemble-based representations of protein structure and dynamics with the use of NMR-derived data. An important common aspect of the methods is that NMR observables and derived parameters are interpreted as properties of the ensemble instead of individual conformers. The resulting ensembles reflect the experimentally determined internal mobility of proteins at a given time scale and can be used to understand the role of internal motions in biological processes at atomic detail. In this review we provide an overview of the calculation methods currently available and examples of biological insights obtained by the ensemble-based models of the proteins investigated.

Applications of isothermal titration calorimetry in pure and applied research--survey of the literature from 2010

Isothermal titration calorimetry (ITC) is a biophysical technique for measuring the formation and dissociation of molecular complexes and has become an invaluable tool in many branches of science from cell biology to food chemistry. By measuring the heat absorbed or released during bond formation, ITC provides accurate, rapid, and label-free measurement of the thermodynamics of molecular interactions. In this review, we survey the recent literature reporting the use of ITC and have highlighted a number of interesting studies that provide a flavour of the diverse systems to which ITC can be applied. These include measurements of protein-protein and protein-membrane interactions required for macromolecular assembly, analysis of enzyme kinetics, experimental validation of molecular dynamics simulations, and even in manufacturing applications such as food science. Some highlights include studies of the biological complex formed by Staphylococcus aureus enterotoxin C3 and the murine T-cell receptor, the mechanism of membrane association of the Parkinson's disease-associated protein α-synuclein, and the role of non-specific tannin-protein interactions in the quality of different beverages. Recent developments in automation are overcoming limitations on throughput imposed by previous manual procedures and promise to greatly extend usefulness of ITC in the future. We also attempt to impart some practical advice for getting the most out of ITC data for those researchers less familiar with the method.

Understanding biomolecular motion, recognition, and allostery by use of conformational ensembles

We review the role conformational ensembles can play in the analysis of biomolecular dynamics, molecular recognition, and allostery. We introduce currently available methods for generating ensembles of biomolecules and illustrate their application with relevant examples from the literature. We show how, for binding, conformational ensembles provide a way of distinguishing the competing models of induced fit and conformational selection. For allostery we review the classic models and show how conformational ensembles can play a role in unravelling the intricate pathways of communication that enable allostery to occur. Finally, we discuss the limitations of conformational ensembles and highlight some potential applications for the future.

CoNSEnsX: an ensemble view of protein structures and NMR-derived experimental data

Background

In conjunction with the recognition of the functional role of internal dynamics of proteins at various timescales, there is an emerging use of dynamic structural ensembles instead of individual conformers. These ensembles are usually substantially more diverse than conventional NMR ensembles and eliminate the expectation that a single conformer should fulfill all NMR parameters originating from 10¹⁶ - 10¹⁷ molecules in the sample tube. Thus, the accuracy of dynamic conformational ensembles should be evaluated differently to that of single conformers.

Results

We constructed the web application CoNSEnsX (Consistency of NMR-derived Structural Ensembles with eXperimental data) allowing fast, simple and convenient assessment of the correspondence of the ensemble as a whole with diverse independent NMR parameters available. We have chosen different ensembles of three proteins, human ubiquitin, a small protease inhibitor and a disordered subunit of cGMP phosphodiesterase 5/6 for detailed evaluation and demonstration of the capabilities of the CoNSEnsX approach.

Conclusions

Our results present a new conceptual method for the evaluation of dynamic conformational ensembles resulting from NMR structure determination. The designed CoNSEnsX approach gives a complete evaluation of these ensembles and is freely available as a web service at http://consensx.chem.elte.hu.

Enzyme dynamics point to stepwise conformational selection in catalysis

Recent data increasingly reveal that conformational dynamics are indispensable to enzyme function throughout the catalytic cycle, in substrate recruiting, chemical transformation, and product release. Conformational transitions may involve conformational selection and induced fit, which can be viewed as a special case in the catalytic network. NMR, X-ray crystallography, single-molecule FRET, and simulations clearly demonstrate that the free enzyme dynamics already encompass all the conformations necessary for substrate binding, preorganization, transition-state stabilization, and product release. Conformational selection and substate population shift at each step of the catalytic turnover can accommodate enzyme specificity and efficiency. Within such a framework, entropy can have a larger role in conformational dynamics than in direct energy transfer in dynamically promoted catalysis.