Bovine Pancreatic Trypsin Inhibitor Research Articles

Multilevel Framework for Analysis of Protein Folding Involving Disulfide Bond Formation.

In this study, a three-layered multicenter ONIOM approach is implemented to characterize the naive folding pathway of bovine pancreatic trypsin inhibitor (BPTI). Each layer represents a distinct level of theory, where the initial layer, encompassing the entire protein, is modeled by a general all-atom force-field GFN-FF. An intermediate electronic structure layer consisting of three multicenter fragments is introduced with the state-of-the-art semiempirical tight-binding method GFN2-xTB. Higher accuracy, specifically addressing the breaking and formation of the three disulfide bonds, is achieved at the innermost layer using the composite DFT method r2SCAN-3c. Our analysis sheds light on the structural stability of BPTI, particularly the significance of interlinking disulfide bonds. The accuracy and efficiency of the multicenter QM/SQM/MM approach are benchmarked using the oxidative formation of cystine. For the folding pathway of BPTI, relative stabilities are investigated through the calculation of free energy contributions for selected intermediates, focusing on the impact of the disulfide bond. Our results highlight the intricate trade-off between accuracy and computational cost, demonstrating that the multicenter ONIOM approach provides a well-balanced and comprehensive solution to describe electronic structure effects in biomolecular systems. We conclude that multiscale energy landscape exploration provides a robust methodology for the study of intriguing biological targets.

Role of water in the pressure-induced unfolding of Bovine Pancreatic Trypsin Inhibitor (BPTI) protein: A molecular dynamics study

The mechanism of protein denaturation under pressure is still an active area of research. Using extensive molecular dynamics simulations, for a partially unfolded structures of BPTI protein, we have studied the role of water in the unfolding of proteins at pressures of 0.001, 5, and 10 kbar at 298K. Our results show that high enough pressure leads to insertion of water into the protein’s interior, causing structural perturbation to proteins, facilitated by a decrease in the intact time of intra-protein hydrogen bond acceptor–donor distances or an increase in the residence times of shell waters or by a combination of both.

In silico study of solvation effects in solutions of biomolecules: possibilities of an approach based on the 3d-distribution of solvent atomic density

Biomolecular solvation plays one of the key roles in nature. The biological activity of molecules and the performance of their target functions depend on the features of this process. However, the study of the biomolecule hydration is a non-trivial task for both experimental methods and computer simulations. The paper demonstrates the possibilities of the non-empirical 3D-SDFT/3D-RISM approach based on the 3D-distribution of the solvent atomic density to study the features of biomolecule hydration using the example of a number of amino acids such as Gly-ZW, L-Ala-ZW, L-Val-ZW, L -Pro-ZW, two model proteins such as BP-TI (bovine pancreatic trypsin inhibitor) and PTP1B (protein tyrosine phosphatase 1B), as well as complexes of the PTP1B protein with inhibitors. The presented results show that the approach allows one to describe in detail and at the same time a holistic description of the hydration shell structure of biomolecules.

Free-Energy Landscape and Rate Estimation of the Aromatic Ring Flips in Basic Pancreatic Trypsin Inhibitors Using Metadynamics.

Aromatic side chains (phenylalanine and tyrosine) of a protein flip by 180° around the Cβ-Cγ axis (χ2 dihedral of the side chain), producing two symmetry-equivalent states. The study of ring flip dynamics with nuclear magnetic resonance (NMR) experiments helps to understand local conformational fluctuations. Ring flips are categorized as slow (milliseconds and onward) or fast (nanoseconds to near milliseconds) based on timescales accessible to NMR experiments. In this study, we investigated the ability of the infrequent metadynamics approach to estimate the flip rate and discriminate between slow and fast ring flips for eight individual aromatic side chains (F4, Y10, Y21, F22, Y23, F33, Y35, and F45) of the basic pancreatic trypsin inhibitor. Well-tempered metadynamics simulations were performed to estimate the ring-flipping free-energy surfaces for all eight aromatic residues. The results indicate that χ2 as a standalone collective variable (CV) is not sufficient to obtain computationally consistent results. Inclusion of a complementary CV, such as χ1(Cα-Cβ), solved the problem for most residues and enabled us to classify fast and slow ring flips. This indicates the importance of librational motions in ring flips. Multiple pathways and mechanisms were observed for residues F4, Y10, and F22. Recrossing events were observed for residues F22 and F33, indicating a possible role of friction effects in ring flipping. The results demonstrate the successful application of infrequent metadynamics to estimate ring flip rates and identify certain limitations of the approach.

Abstract 138: Treatment Initiation During The First Month Following Hypertension Diagnosis Improves Blood Pressure Control Up To 30 Months

Background: Adults with hypertension (HTN) have fewer cardiovascular events if controlled within the first 6 months after diagnosis, during which time they are excluded from many HTN control metrics. We compared BP control rates (BP <140/90) from 6 to 42 months in adults with HTN who did or did not have treatment initiated during the first month following diagnosis, irrespective of subsequent treatment initiation or intensification (TI). Methods: A patient cohort from 5 diverse healthcare organizations, engaged in a quality improvement program, was identified with previously undiagnosed and untreated HTN. TI and BP control were assessed by measuring the period prevalence of visit BP control and TI on uncontrolled visits, then compared with the time since diagnosis (TSD) up to 42 months. Groups were stratified by treatment initiation status during the first month. Results: A prospective cohort of 14,328 newly diagnosed and untreated persons with HTN was identified, mean age 56.2±14.7 years. Patients with TI during the first month had higher rates of BP control at 6 (56.2% vs 46.2%, p<0.001) through 30 months (65.9% vs 61.4%, p<0.001) than those initially untreated. Patients initiated during the first month were less likely to receive additional intensifications (OR 0.68, 95% CI [0.63-0.74]). Conclusion: TI in the first month after diagnosis results in better BP control at 6 months, which has prognostic importance and coincides with the time when these individuals impact standard HTN control metrics. Of concern, TI during the first month after diagnosis was insufficient to achieve control rates ≥70% and lowered the probability of subsequent TI.

Probing protein stability: towards a computational atomistic, reliable, affordable, and improvable model.

We present an improved application of a recently proposed computational method designed to evaluate the change of free energy as a function of the average value of a suitably chosen collective variable in proteins. The method is based on a full atomistic description of the protein and its environment. The goal is to understand how the protein melting temperature changes upon single-point mutations, because the sign of the temperature variation will allow us to discriminate stabilizing vs. destabilizing mutations in protein sequences. In this refined application the method is based on altruistic well-tempered metadynamics, a variant of multiple-walkers metadynamics. The resulting metastatistics is then modulated by the maximal constrained entropy principle. The latter turns out to be especially helpful in free-energy calculations as it is able to alleviate the severe limitations of metadynamics in properly sampling folded and unfolded configurations. In this work we apply the computational strategy outlined above in the case of the bovine pancreatic trypsin inhibitor, a well-studied small protein, which is a reference for computer simulations since decades. We compute the variation of the melting temperature characterizing the folding-unfolding process between the wild-type protein and two of its single-point mutations that are seen to have opposite effect on the free energy changes. The same approach is used for free energy difference calculations between a truncated form of frataxin and a set of five of its variants. Simulation data are compared to in vitro experiments. In all cases the sign of the change of melting temperature is reproduced, under the further approximation of using an empirical effective mean-field to average out protein-solvent interactions.

Accounting for Solvation Correlation Effects on the Thermodynamics of Water Networks in Protein Cavities

Macromolecular recognition and ligand binding are at the core of biological function and drug discovery efforts. Water molecules play a significant role in mediating the protein-ligand interaction, acting as more than just the surrounding medium by affecting the thermodynamics and thus the outcome of the binding process. As individual water contributions are impossible to measure experimentally, a range of computational methods have emerged to identify hydration sites in protein pockets and characterize their energetic contributions for drug discovery applications. Even though several methods model solvation effects explicitly, they focus on determining the stability of specific water sites independently and neglect solvation correlation effects upon replacement of clusters of water molecules, which typically happens in hit-to-lead optimization. In this work, we rigorously determine the conjoint effects of replacing all combinations of water molecules in protein binding pockets through the use of the RE-EDS multistate free-energy method, which combines Hamiltonian replica exchange (RE) and enveloping distribution sampling (EDS). Applications on the small bovine pancreatic trypsin inhibitor and four proteins of the bromodomain family illustrate the extent of solvation correlation effects on water thermodynamics, with the favorability of replacement of the water sites by pharmacophore probes highly dependent on the composition of the water network and the pocket environment. Given the ubiquity of water networks in biologically relevant protein targets, we believe our approach can be helpful for computer-aided drug discovery by providing a pocket-specific and a priori systematic consideration of solvation effects on ligand binding and selectivity.

Correlation between protein conformations and water structure and thermodynamics at high pressure: A molecular dynamics study of the Bovine Pancreatic Trypsin Inhibitor (BPTI) protein.

Pressure-induced perturbation of a protein structure leading to its folding-unfolding mechanism is an important yet not fully understood phenomenon. The key point here is the role of water and its coupling with protein conformations as a function of pressure. In the current work, using extensive molecular dynamics simulation at 298K, we systematically examine the coupling between protein conformations and water structures of pressures of 0.001, 5, 10, 15, 20 kbar, starting from (partially) unfolded structures of the protein Bovine Pancreatic Trypsin Inhibitor (BPTI). We also calculate localized thermodynamics at those pressures as a function of protein-water distance. Our findings show that both protein-specific and generic effects of pressure are operating. In particular, we found that (1) the amount of increase in water density near the protein depends on the protein structural heterogeneity; (2) the intra-protein hydrogen bond decreases with pressure, while the water-water hydrogen bond per water in the first solvation shell (FSS) increases; protein-water hydrogen bonds also found to increase with pressure, (3) with pressure hydrogen bonds of waters in the FSS getting twisted; and (4) water's tetrahedrality in the FSS decreases with pressure, but it is dependent on the local environment. Thermodynamically, at higher pressure, the structural perturbation of BPTI is due to the pressure-volume work, while the entropy decreases with the increase of pressure due to the higher translational and rotational rigidity of waters in the FSS. The local and subtle effects of pressure, found in this work, are likely to be typical of pressure-induced protein structure perturbation.

Oxidative Protein Folding Promotion by Imidazoyl-conjugated Thiol

The promotion of protein folding is an important subject that can contribute to the efficient production of biological medicines. Here we report the promotion of disulfide bond-coupled protein folding by imidazoyl-conjugated thiol (ImdSH). ImdSH accelerated the disulfide bond formation and folding of reduced ribonuclease A and bovine pancreatic trypsin inhibitor relative to glutathione as a conventionally used additive.

The Non-native Disulfide-Bond-Containing Landscape Orthogonal to the Oxidative Protein-Folding Trajectory: A Necessary Evil?

Oxidative protein folding describes the process by which disulfide-bond-containing proteins mature from their ribosomal, fully reduced and unfolded, origins. Over the past 40 years, a number of exemplar proteins including bovine pancreatic ribonuclease A (RNaseA), bovine pancreatic trypsin inhibitor (BPTI), and hen egg-white lysozyme (HEWL), among others, have provided rich insight into the nature of the intermolecular interactions that drive the formation of the native, biologically active fold. In this Review Article, we revisit the oxidative folding process of RNaseA with a focus on reconciling the role of non-native disulfide-bond-containing species that populate the oxidative folding landscape. Toward gaining such an understanding, we project the regeneration pathway onto a Cartesian coordinate system. This helps not only to recognize the magnitude of the seemingly "fruitless", non-native disulfide-bond-containing species that lie orthogonal to the "native-protein-forming" reaction progress but also to reconcile a role for their existence in the regenerative trajectory. Finally, we superimpose the folding funnel onto the regeneration trajectory to draw parallels between oxidative folders and conformational folders (proteins that lack disulfide bonds). The overall objective is to provide the reader with a semi-quantitative description of oxidative protein folding and the barriers to successful regeneration while underscoring a role of seemingly fruitless intermediates.

Protein 3D Hydration: A Case of Bovine Pancreatic Trypsin Inhibitor

Characterization of the hydrated state of a protein is crucial for understanding its structural stability and function. In the present study, we have investigated the 3D hydration structure of the protein BPTI (bovine pancreatic trypsin inhibitor) by molecular dynamics (MD) and the integral equation method in the three-dimensional reference interaction site model (3D-RISM) approach. Both methods have found a well-defined hydration layer around the protein and revealed the localization of BPTI buried water molecules corresponding to the X-ray crystallography data. Moreover, under 3D-RISM calculations, the obtained positions of waters bound firmly to the BPTI sites are in reasonable agreement with the experimental results mentioned above for the BPTI crystal form. The analysis of the 3D hydration structure (thickness of hydration shell and hydration numbers) was performed for the entire protein and its polar and non-polar parts using various cut-off distances taken from the literature as well as by a straightforward procedure proposed here for determining the thickness of the hydration layer. Using the thickness of the hydration shell from this procedure allows for calculating the total hydration number of biomolecules properly under both methods. Following this approach, we have obtained the thickness of the BPTI hydration layer of 3.6 Å with 369 water molecules in the case of MD simulation and 3.9 Å with 333 water molecules in the case of the 3D-RISM approach. The above procedure was also applied for a more detailed description of the BPTI hydration structure near the polar charged and uncharged radicals as well as non-polar radicals. The results presented for the BPTI as an example bring new knowledge to the understanding of protein hydration.

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 emergence of protein dynamics simulations: how computational statistical mechanics met biochemistry

In this essay, we aim to illustrate how Martin Karplus and his research group effectively set in motion the engine of molecular dynamics (MD) simulations of biomolecules. This process saw its prodromes between 1969 and the early 1970s with Karplus’ landing in biology, a transition that came to fruition with the treatment of 11-cis-retinal photoisomerization and the development of an allosteric model to account for the mechanism of cooperativity in hemoglobin. In 1977, J. Andrew McCammon, Bruce Gelin, and Martin Karplus published an article in Nature reporting the MD simulation of bovine pancreatic trypsin inhibitor (BPTI). This publication helped initiate the merger of computational statistical mechanics and biochemistry, a process that Karplus undertook at a later stage and whose beginnings we propose to reconstruct in this article through unpublished accounts of the key people who participated in this endeavor.

Characterizing Metastable States with the Help of Machine Learning

Present-day atomistic simulations generate long trajectories of ever more complex systems. Analyzing these data, discovering metastable states, and uncovering their nature are becoming increasingly challenging. In this paper, we first use the variational approach to conformation dynamics to discover the slowest dynamical modes of the simulations. This allows the different metastable states of the system to be located and organized hierarchically. The physical descriptors that characterize metastable states are discovered by means of a machine learning method. We show in the cases of two proteins, chignolin and bovine pancreatic trypsin inhibitor, how such analysis can be effortlessly performed in a matter of seconds. Another strength of our approach is that it can be applied to the analysis of both unbiased and biased simulations.

Bovine pancreatic trypsin inhibitor and soybean Kunitz trypsin inhibitor: Differential effects on proteases and larval development of the soybean pest Anticarsia gemmatalis (Lepidoptera: Noctuidae)

Pest management is challenged with resistant herbivores and problems regarding human health and environmental issues. Indeed, the greatest challenge to modern agriculture is to protect crops from pests and still maintain environmental quality. This study aimed to analyze by in silico, in vitro, and in vivo approaches to the feasibility of using the inhibitory protein extracted from mammals – Bovine Pancreatic Trypsin Inhibitor (BPTI) as a potential inhibitor of digestive trypsins from the pest Anticarsia gemmatalis and comparing the results with the host-plant inhibitor – Soybean Kunitz Trypsin Inhibitor (SKTI). BPTI and SKTI interacts with A. gemmatalis trypsin-like enzyme competitively, through hydrogen and hydrophobic bonds. A. gemmatalis larvae exposed to BPTI did not show two common adaptative mechanisms i.e., proteolytic degradation and overproduction of proteases, presenting highly reduced trypsin-like activity. On the other hand, SKTI-fed larvae did not show reduced trypsin-like activity, presenting overproduction of proteases and SKTI digestion. In addition, the larval survival was reduced by BPTI similarly to SKTI, and additionally caused a decrease in pupal weight. The non-plant protease inhibitor BPTI presents intriguing element to compose biopesticide formulations to help decrease the use of conventional refractory pesticides into integrated pest management programs.

Py-MCMD: Python Software for Performing Hybrid Monte Carlo/Molecular Dynamics Simulations with GOMC and NAMD.

py-MCMD, an open-source Python software, provides a robust workflow layer that manages communication of relevant system information between the simulation engines NAMD and GOMC and generates coherent thermodynamic properties and trajectories for analysis. To validate the workflow and highlight its capabilities, hybrid Monte Carlo/molecular dynamics (MC/MD) simulations are performed for SPC/E water in the isobaric-isothermal (NPT) and grand canonical (GC) ensembles as well as with Gibbs ensemble Monte Carlo (GEMC). The hybrid MC/MD approach shows close agreement with reference MC simulations and has a computational efficiency that is 2 to 136 times greater than traditional Monte Carlo simulations. MC/MD simulations performed for water in a graphene slit pore illustrate significant gains in sampling efficiency when the coupled-decoupled configurational-bias MC (CD-CBMC) algorithm is used compared with simulations using a single unbiased random trial position. Simulations using CD-CBMC reach equilibrium with 25 times fewer cycles than simulations using a single unbiased random trial position, with a small increase in computational cost. In a more challenging application, hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations are used to hydrate a buried binding pocket in bovine pancreatic trypsin inhibitor. Water occupancies produced by GCMC/MD simulations are in close agreement with crystallographically identified positions, and GCMC/MD simulations have a computational efficiency that is 5 times better than MD simulations. py-MCMD is available on GitHub at https://github.com/GOMC-WSU/py-MCMD.

One-Pot Chemical Protein Synthesis Utilizing Fmoc-Masked Selenazolidine to Address the Redox Functionality of Human Selenoprotein F.

Human SELENOF is an endoplasmic reticulum (ER) selenoprotein that contains the redox active motif CXU (C is cysteine and U is selenocysteine), resembling the redox motif of thiol‐disulfide oxidoreductases (CXXC). Like other selenoproteins, the challenge in accessing SELENOF has somewhat limited its full biological characterization thus far. Here we present the one‐pot chemical synthesis of the thioredoxin‐like domain of SELENOF, highlighted by the use of Fmoc‐protected selenazolidine, native chemical ligations and deselenization reactions. The redox potential of the CXU motif, together with insulin turbidimetric assay suggested that SELENOF may catalyze the reduction of disulfides in misfolded proteins. Furthermore, we demonstrate that SELENOF is not a protein disulfide isomerase (PDI)‐like enzyme, as it did not enhance the folding of the two protein models; bovine pancreatic trypsin inhibitor and hirudin. These studies suggest that SELENOF may be responsible for reducing the non‐native disulfide bonds of misfolded glycoproteins as part of the quality control system in the ER.

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.

Microscopic insights into dynamic disorder in the isomerization dynamics of the protein BPTI.

Understanding the dynamic disorder behind a process, i.e., the dynamic effect of fluctuations that occur on a timescale slower or comparable with the timescale of the process, is essential for elucidating the dynamics and kinetics of complicated molecular processes in biomolecules and liquids. Despite numerous theoretical studies of single-molecule kinetics, our microscopic understanding of dynamic disorder remains limited. In the present study, we investigate the microscopic aspects of dynamic disorder in the isomerization dynamics of the Cys14-Cys38 disulfide bond in the protein bovine pancreatic trypsin inhibitor, which has been observed by nuclear magnetic resonance. We use a theoretical model with a stochastic transition rate coefficient, which is calculated from the 1-ms-long time molecular dynamics trajectory obtained by Shaw et al. [Science 330, 341-346 (2010)]. The isomerization dynamics are expressed by the transitions between coarse-grained states consisting of internal states, i.e., conformational sub-states. In this description, the rate for the transition from the coarse-grained states is stochastically modulated due to fluctuations between internal states. We examine the survival probability for the conformational transitions from a coarse-grained state using a theoretical model, which is a good approximation to the directly calculated survival probability. The dynamic disorder changes from a slow modulation limit to a fast modulation limit depending on the aspects of the coarse-grained states. Our analysis of the rate modulations behind the survival probability, in relation to the fluctuations between internal states, reveals the microscopic origin of dynamic disorder.

Amblyomma sculptum Salivary Protease Inhibitors as Potential Anti-Tick Vaccines.

Amblyomma sculptum is the main tick associated with human bites in Brazil and the main vector of Rickettsia rickettsii, the causative agent of the most severe form of Brazilian spotted fever. Molecules produced in the salivary glands are directly related to feeding success and vector competence. In the present study, we identified sequences of A. sculptum salivary proteins that may be involved in hematophagy and selected three proteins that underwent functional characterization and evaluation as vaccine antigens. Among the three proteins selected, one contained a Kunitz_bovine pancreatic trypsin inhibitor domain (named AsKunitz) and the other two belonged to the 8.9 kDa and basic tail families of tick salivary proteins (named As8.9kDa and AsBasicTail). Expression of the messenger RNA (mRNA) encoding all three proteins was detected in the larvae, nymphs, and females at basal levels in unfed ticks and the expression levels increased after the start of feeding. Recombinant proteins rAs8.9kDa and rAsBasicTail inhibited the enzymatic activity of factor Xa, thrombin, and trypsin, whereas rAsKunitz inhibited only thrombin activity. All three recombinant proteins inhibited the hemolysis of both the classical and alternative pathways; this is the first description of tick members of the Kunitz and 8.9kDa families being inhibitors of the classical complement pathway. Mice immunization with recombinant proteins caused efficacies against A. sculptum females from 59.4% with rAsBasicTail immunization to more than 85% by immunization with rAsKunitz and rAs8.9kDa. The mortality of nymphs fed on immunized mice reached 70–100%. Therefore, all three proteins are potential antigens with the possibility of becoming a new tool in the control of A. sculptum.