Bovine Pancreatic Trypsin Inhibitor Research Articles

Chemical Synthesis and Structure-Activity Relationship Studies of the Coagulation Factor Xa Inhibitor Tick Anticoagulant Peptide from the Hematophagous Parasite Ornithodoros moubata.

Tick Anticoagulant Peptide (TAP), a 60-amino acid protein from the soft tick Ornithodoros moubata, inhibits activated coagulation factor X (fXa) with almost absolute specificity. Despite TAP and Bovine Pancreatic Trypsin Inhibitor (BPTI) (i.e., the prototype of the Kunitz-type protease inhibitors) sharing a similar 3D fold and disulphide bond topology, they have remarkably different amino acid sequence (only ~24% sequence identity), thermal stability, folding pathways, protease specificity, and even mechanism of protease inhibition. Here, fully active and correctly folded TAP was produced in reasonably high yields (~20%) by solid-phase peptide chemical synthesis and thoroughly characterised with respect to its chemical identity, disulphide pairing, folding kinetics, conformational dynamics, and fXa inhibition. The versatility of the chemical synthesis was exploited to perform structure-activity relationship studies on TAP by incorporating non-coded amino acids at positions 1 and 3 of the inhibitor. Using Hydrogen-Deuterium Exchange Mass Spectrometry, we found that TAP has a remarkably higher conformational flexibility compared to BPTI, and propose that these different dynamics could impact the different folding pathway and inhibition mechanisms of TAP and BPTI. Hence, the TAP/BPTI pair represents a nice example of divergent evolution, while the relative facility of TAP synthesis could represent a good starting point to design novel synthetic analogues with improved pharmacological profiles.

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.

A production platform for disulfide-bonded peptides in the periplasm of Escherichia coli

BackgroundRecombinant peptide production in Escherichia coli provides a sustainable alternative to environmentally harmful and size-limited chemical synthesis. However, in-vivo production of disulfide-bonded peptides at high yields remains challenging, due to degradation by host proteases/peptidases and the necessity of translocation into the periplasmic space for disulfide bond formation.ResultsIn this study, we established an expression system for efficient and soluble production of disulfide-bonded peptides in the periplasm of E. coli. We chose model peptides with varying complexity (size, structure, number of disulfide bonds), namely parathyroid hormone 1–84, somatostatin 1–28, plectasin, and bovine pancreatic trypsin inhibitor (aprotinin). All peptides were expressed without and with the N-terminal, low molecular weight CASPON™ tag (4.1 kDa), with the expression cassette being integrated into the host genome. During BioLector™ cultivations at microliter scale, we found that most of our model peptides can only be sufficiently expressed in combination with the CASPON™ tag, otherwise expression was only weak or undetectable on SDS-PAGE. Undesired degradation by host proteases/peptidases was evident even with the CASPON™ tag. Therefore, we investigated whether degradation happened before or after translocation by expressing the peptides in combination with either a co- or post-translational signal sequence. Our results suggest that degradation predominantly happened after the translocation, as degradation fragments appeared to be identical independent of the signal sequence, and expression was not enhanced with the co-translational signal sequence. Lastly, we expressed all CASPON™-tagged peptides in two industry-relevant host strains during C-limited fed-batch cultivations in bioreactors. We found that the process performance was highly dependent on the peptide-host-combination. The titers that were reached varied between 0.6–2.6 g L−1, and exceeded previously published data in E. coli. Moreover, all peptides were shown by mass spectrometry to be expressed to completion, including full formation of disulfide bonds.ConclusionIn this work, we demonstrated the potential of the CASPON™ technology as a highly efficient platform for the production of soluble peptides in the periplasm of E. coli. The titers we show here are unprecedented whenever parathyroid hormone, somatostatin, plectasin or bovine pancreatic trypsin inhibitor were produced in E. coli, thus making our proposed upstream platform favorable over previously published approaches and chemical synthesis.

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.

A growth type pathway for improving the folding of BPTI.

The in vitro oxidative folding of the protein bovine pancreatic trypsin inhibitor (BPTI) with oxidized dithiothreitol or glutathione has served as a paradigm for protein folding but could take weeks at physiological pH because of the need to escape from kinetic traps via a rearrangement type pathway. The two major kinetic traps are called N' and N* and contain two of the three native disulfide bonds, which occur between residues 5 and 55, 30 and 51, and 14 and 38. N' is missing the disulfide bond between residues 5 and 55 while N* is missing the disulfide bond between residues 30 and 51. By determining rate constant for the reactions of the kinetic traps N* and N' and their mixed disulfides with glutathione and glutathione disulfide, many for the first time, we demonstrate that growth type pathways are feasible and could even be more efficient than rearrangement type pathways. Thus, formally unproductive pathways became productive. Interestingly, under physiological redox conditions both rearrangement and growth type pathways are important highlighting the redundancy of oxidative protein folding. With the new set of rate constants, modeling indicated that in vitro oxidative protein folding of BPTI via a growth type pathway using an oxidation, reduction and oxidation cycle would significantly improve protein folding efficiency, albeit under non-physiological redox conditions. With these changing conditions 91 ± 2% of native BPTI was achieved in 12 h compared to 83% native protein in 24 h using our previous best conditions of 5 mM GSSG and 5 mM GSH. Therefore, changing redox conditions via an oxidation, reduction and oxidation cycle may become an additional methodology for enhancing in vitro protein folding in aqueous solution.

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.