Purification and characterization of trypsin produced by gut bacteria from Anticarsia gemmatalis

Purification of active trypsin in the digestive process of insects is essential for the development of potent protease inhibitors (PIs) as an emerging pest control technology and research into insect adaptations to dietary PIs. An important aspect is the presence of proteolytic microorganisms, which contribute to host nutrition. Here, we purified trypsins produced by bacteria Bacillus cereus, Enterococcus mundtii, Enterococcus gallinarum, and Staphylococcus xylosus isolated from the midgut of Anticarsia gemmatalis. The trypsins had a molecular mass of approximately 25 kDa. The enzymes showed increased activity at 40°C, and they were active at pH values 7.5–10. Aprotinin, bis‐benzamidine, and soybean Kunitz inhibitor (SKTI) significantly inhibited trypsin activity. The l ‐1‐tosyl‐amido‐2‐phenylethylchloromethyl ketone (TPCK), pepstatin A, E‐64, ethylenediamine tetraacetic acid, and calcium ions did not affect the enzyme activity at the concentrations tested. We infer the purified trypsins do not require calcium ions, by which they differ from the trypsins of other microorganisms and the soluble and insoluble trypsins characterized from A. gemmatalis. These data suggest the existence of different isoforms of trypsin in the velvetbean caterpillar midguts.     

Refinement of the structure of human basic fibroblast growth factor at 1.6 A resolution and analysis of presumed heparin binding sites by selenate substitution.

The three-dimensional structure of human basic fibroblast growth factor has been refined to a crystallographic residual of 16.1% at 1.6 A resolution. The structure has a Kunitz-type fold and is composed of 12 antiparallel beta-strands, 6 of which form a beta-barrel. One bound sulfate ion has been identified in the model, hydrogen bonded to the side chains of Asn 27, Arg 120, and Lys 125. The side chain of Arg 120 has two conformations, both of which permit hydrogen bonds to the sulfate. This sulfate binding site has been suggested as the binding site for heparin (Eriksson, A.E., Cousens, L.S., Weaver, L.H., & Matthews, B.W., 1991, Proc. Natl. Acad. Sci. USA 88, 3441-3445). Two beta-mercaptoethanol (BME) molecules are also included in the model, each forming a disulfide bond to the S gamma atoms of Cys 69 and Cys 92, respectively. The side chain of Cys 92 has two conformations of which only one can bind BME. Therefore the BME molecule is half occupied at this site. The locations of possible sulfate binding sites on the protein were examined by replacing the ammonium sulfate in the crystallization medium with ammonium selenate. Diffraction data were measured to 2.2 A resolution and the structure refined to an R-factor of 13.8%. The binding of the more electron-dense selenate ion was identified at two positions. One position was identical to the sulfate binding site identified previously. The second selenate binding site, which is of lower occupancy, is situated 5.6 A from the first. This ion is hydrogen bonded by the side chain of Lys 135 and Arg 120. Thus the side chain of Arg 120 binds two selenate ions simultaneously. It is suggested that the observed second selenate binding site should also be considered as a possible binding site for heparin, or that both selenate binding sites might simultaneously contribute to the binding of heparin.     

Molecular mechanism of enzyme inhibition: prediction of the three-dimensional structure of the dimeric trypsin inhibitor from Leucaena leucocephala by homology modelling

Serine proteinase inhibitors are widely distributed in nature and inhibit the activity of enzymes like trypsin and chymotrypsin. These proteins interfere with the physiological processes such as germination, maturation and form the first line of defense against the attack of seed predator. The most thoroughly examined plant serine proteinase inhibitors are found in the species of the families Leguminosae, Graminae, and Solanaceae. Leucaena leucocephala belongs to the family Leguminosae. It is widely used both as an ornamental tree as well as cattle food. We have constructed a three-dimensional model of a serine proteinase inhibitor from L. leucocephala seeds (LTI) complexed with trypsin. The model was built based on its comparative homology with soybean trypsin inhibitor (STI) using the program, MODELLER6. The quality of the model was assessed stereochemically by PROCHECK. LTI shows structural features characteristic of the Kunitz type trypsin inhibitor and shows 39% residue identity with STI. LTI consists of 172 amino acid residues and is characterized by two disulfide bridges. The protein is a dimer with the two chains being linked by a disulfide bridge. Despite the high similarity in the overall tertiary structure, significant differences exist at the active site between STI and LTI. The present study aims at analyzing these interactions based on the available amino acid sequences and structural data. We have also studied some functional sites such as phosphorylation, myristoylation, which can influence the inhibitory activity or complexation with other molecules. Some of the differences observed at the active site and functional sites can explain the unique features of LTI.     

Loss of soybean trypsin inhibitor in callus as monitored by inhibition enzyme immunoassay

Summary A monoclonal antibody was produced against Kunitz soybean inhibitor (KSBTI) and used in an inhibition enzyme immunoassay (EIA). The inhibition EIA was as sensitive as competetive EIAs and was easily modified for other protein-antibody interactions. The KSBTI assay described detected KSBTI in complex mixtures from 100 μg/ml to 50 ng/ml and did not react with the Bowman-Birk trypsin inhibitor. The assay was used to examine levels of KSBTI inGlycine max hypocotyl-derived callus tissue. The developing hypocotyls contained 0.21 μg KSBTI per mg of fresh tissue. This level of KSBTI rapidly decreased when placed in culture and was undetectable 6 days later. The decrease in KSBTI correlated with the development of callus.     

Ligula intestinalis (Cestoda: Pseudophyllidae): Studies on the properties of proteolytic and protease inhibitor activities of plerocercoid larvae

The somatic extract of L. intestinalis plerocercoids reveals hydrolytic activity against N-Benzoyl-l-tyrosine ethyl ester (BTEE) and Azocoll, and inactivates the esterolysis by mammalian trypsin and chymotripsin. The proteolytic enzyme activity and the inhibitory effect were completely separated by Sephadex G-100 column chromatography. Gel chromatography of the somatic extract revealed two peaks of proteolytic activity : one is bound to macromolecular substances, the other appears to be in free form and has a molecular weight of approx 60,000–65,000. The proteolytic activity showed the following characteristics : Tris-HCl buffer provided the highest activity against BTEE, the pH optimum was 7·4–7·8; the enzyme was activated by 10^?5m-Ca²⁺, Mg²⁺ or Mn²⁺, it was inhibited by 10^?5m-Cu²⁺, but not by 10^?5m-Zn²⁺. 0.001% soybean trypsin inhibitor, 2 × 10^?3m-EDTA, 1 mm-tosyl-l-phenylalanyl chloromethane, 1000 KIU/ml Trasylol did not inhibit the proteolytic activity, but it was inhibited by 1 mm-phenylmethyl-sulphonyl fluoride. The enzyme activity completely ceased upon 5 % TCA treatment or incubation at 56°C for 30 min. The trypsin and chyrnotrypsin inhibitor activities were eluted from the Sephadex G-100 column in a single peak with an estimated molecular weight of 6700–7200. The inhibitory effect was not sensitive to pH changes, and treatment by 5% TCA or incubation at 80°C for 15 min was ineffective. The proteolytic activity of plerocercoid extract was not effected ‘in vitro’ by the inhibitors isolated from this parasite.     

Effects of soybean trypsin inhibitor on hypopharyngeal gland protein content, total midgut protease activity and survival of the honey bee (Apis mellifera L.)

Insecticidal properties of protease inhibitors have been established in transgenic plants. In the wake of continuous research and rapid development of protease inhibitors it is important to assess possible effects on beneficial insects like the honey bee (Apis mellifera L.). In this study, newly emerged caged bees were fed pollen diets containing three different concentrations (0.1%, 0.5% and 1% w:w) of soybean trypsin inhibitor (SBTI). Hypopharyngeal gland protein content, total midgut proteolytic enzyme activity of these bees, and survival were measured. Bees fed 1% SBTI had significantly reduced hypopharyngeal gland protein content and midgut proteolytic enzyme activity. There were no significant differences between control, 0.1% and 0.5% SBTI treatments. Bees fed a diet containing 1% SBTI had the lowest survival, followed by 0.5% and 0.1%, over a 30-day period. We concluded that nurse bees fed a pollen diet containing at least 1% SBTI would be poor producers of larval food, potentially threatening colony growth and maintenance.     

Formation of the tight-binding inhibitor, 3-ketoarabinitol-1,5-bisphosphate by ribulose-1,5-bisphosphate carboxylase/oxygenase is O2-dependent

Ribulose-1,5-bisphosphate carboxylase/oxygenase (Rubisco) (EC 4.1.1.39) not only catalyzes carboxylation and oxygenation of ribulose-1,5-bisphosphate (RuBP), but it can also act either as an epimerase or isomerase converting RuBP into xylulose-1,5-bisphosphate (XuBP) or 3-ketoarabinitol-1,5-bisphosphate (KABP), respectively, a process called misfire. XuBP is formed as a result of misprotonation at C3 of the RuBP-enediol. It is released from Rubisco active sites and accumulates in the reaction mixture. Increasing the amounts of CO₂ or O₂ decreases XuBP production. However, KABP synthesis, which has been proposed to be only a product due to C2 misprotonation of the RuBP-endiol, is dependent upon the presence of O₂. KABP remains tightly bound to Rubisco active sites after its formation, causing the loss of Rubisco activity (fallover). The results suggest that the non-stabilized form of the peroxy-intermediate in the oxygenase reaction can be converted in a backreaction to KABP and molecular oxygen. The stabilization of the peroxy-intermediate due to the presence of Mn²⁺ instead of Mg²⁺ eliminates the formation of KABP.     

Inga vera trypsin inhibitor interferes in the proteolytic activity and nutritional physiology of Ephestia kuehniella larvae

Plant peptidase inhibitors provide plants with a defense strategy to inhibit insect digestive enzymes and have been studied as an alternative strategy for pest control as they interfere in normal insect physiology. We evaluated the effects of ingestion of the trypsin inhibitor from Inga vera Willd. (Fabaceae) seeds on the nutritional and digestive physiology of Ephestia kuehniella (Zeller) (Lepidoptera: Pyralidae) larvae. Inga vera trypsin inhibitor (IVTI) reduced the efficiency of the conversion of ingested and digested food in these larvae and increased the metabolic cost, causing an anti‐nutritional effect. In both short‐ and long‐term bioassays, the ingestion of IVTI inactivated most of the insect's trypsin activity, but increased chymotrypsin activity as a compensatory response by the insect; however, protein digestion continued to be partially blocked. Consequently, chymotrypsin‐like enzymes, which were over‐produced in the gut, were excreted more into the frass of IVTI‐fed larvae. As such, the resistance of IVTI to hydrolysis by insect midgut proteases resulted in detrimental effects to larvae. These data provide support for the use of IVTI as a biotechnological tool for pest control.     

A biophysical and computational study unraveling the molecular interaction mechanism of a new Janus kinase inhibitor Tofacitinib with bovine serum albumin

The interaction of a recently certified kinase inhibitor Tofacitinib (TFB) with bovine serum albumin (BSA) has been studied, by spectroscopic and molecular docking studies. Spectrofluorimetric measurements at 3 different temperatures (288, 298, and 310 K) showed that TFB quench the intrinsic fluorescence of BSA upon forming a nonfluorescent complex. The intrinsic fluorescence data showed that TFB binds to BSA with binding constant (K _b) of approximately 10⁴M⁻¹, affirming a significant affinity of TFB with BSA. The decrease in Stern‐Volmer quenching constant with increasing temperature exhibited the static mechanism of quenching. Negative value of ΔG (−6.94 ± 0.32 kcal·mol⁻¹), ΔH (−7.87 ± 0.52 kcal·mol⁻¹), and ΔS (−3.14 ± 0.42 cal·mol⁻¹·K⁻¹) at all 3 temperatures declared the reaction between BSA and TFB to be spontaneous and exothermic. Far‐UV circular dichroism spectroscopy results demonstrated an increase in helical content of BSA in the presence of TFB. Moreover, dynamic light scattering measurements showed that TFB resulted into a decrease in the hydrodynamic radii (from 3.6 ± 0.053 to 2.9 ± 0.02 nm) of BSA. Molecular docking studies confirmed that TFB binds near site II on BSA, hydrogen bonding, and hydrophobic interaction were involved in the BSA‐TFB complex formation. The present study characterizing the BSA‐TFB interaction could be significant towards gaining an insight into the drug pharmacokinetics and pharmacodynamics and also in the direction of rational drug designing with better competence, against emerging immune‐mediated diseases, ie, alopecia and rheumatoid arthritis.     

Discovery of potential sclerostin inhibitors from plants with loop2 region of sclerostin inhibition by interacting with residues outside Pro-Asn-Ala-Ile-Gly motif

Abstract

Sclerostin, an antagonist of the Wnt/β-catenin signaling pathway, was discovered as a potential therapeutic target for stimulating bone formation in osteoporosis. In this study, molecular docking was employed to predict the binding of 29 herbal compounds, which were reported as bone formation stimulators, to the loop2 region of sclerostin. Then, the 50 ns molecular dynamics (MD) simulation of the complexes between sclerostin and the top 10 hits obtained from molecular docking were carried out. Root mean square deviations (RMSDs) analysis of MD trajectories pointed out that all ligands-complexes remain stable throughout the duration of MD simulations. In addition, the molecular mechanics/generalized born surface area (MM/GBSA) binding free energy and energy decomposition analyses were determined. The results here suggested that baicalin is the most promising inhibitor of sclerostin. Interestingly, baicalin binds to sclerostin via the hydrophobic interaction with the amino acid residues on loop2 region but outside the Pro-Asn-Ala-Ile-Gly (PNAIG) motif, particularly the Arg-Gly-Lys-Trp-Trp-Arg (RGKWWR) motif. This finding could be a novel strategy for developing new sclerostin inhibitors in the future.

Communicated by Ramaswamy H. Sarma     

Inhibition of dihydroceramide desaturase activity by the sphingosine kinase inhibitor SKI II

Sphingosine kinase inhibitor (SKI) II has been reported as a dual inhibitor of sphingosine kinases (SKs) 1 and 2 and has been extensively used to prove the involvement of SKs and sphingosine-1-phosphate (S1P) in cellular processes. Dihydroceramide desaturase (Des1), the last enzyme in the de novo synthesis of ceramide (Cer), regulates the balance between dihydroceramides (dhCers) and Cers. Both SKs and Des1 have interest as therapeutic targets. Here we show that SKI II is a noncompetitive inhibitor (K_i = 0.3 μM) of Des1 activity with effect also in intact cells without modifying Des1 protein levels. Molecular modeling studies support that the SKI II-induced decrease in Des1 activity could result from inhibition of NADH-cytochrome b5 reductase. SKI II, but not the SK1-specific inhibitor PF-543, provoked a remarkable accumulation of dhCers and their metabolites, while both SKI II and PF-543 reduced S1P to almost undetectable levels. SKI II, but not PF543, reduced cell proliferation with accumulation of cells in the G0/G1 phase. SKI II, but not PF543, induced autophagy. These overall findings should be taken into account when using SKI II as a pharmacological tool, as some of the effects attributed to decreased S1P may actually be caused by augmented dhCers and/or their metabolites.     

Designing cyclic peptide inhibitor of dengue virus NS3-NS2B protease by using molecular docking approach

Peptides are preferred for designing inhibitors because of their high activity and specificity. Seven cyclopentapeptide inhibitors were designed in this study against dengue virus type 2 (DEN-2) NS3-NS2B protease: CKRRC, CGRRC, CRGRC, CRTRC, CTRRC, CKRKC and CRRKC. Docking analysis was performed to study the enzyme-inhibitor binding interactions. The free energy binding and estimated Ki values for all the inhibitors were found to be small (within micromolar range), indicating that the inhibitors bind considerably well to the binding site. The results showed that the cyclopentapeptide CKRKC was the best peptide inhibitor candidate with estimated free binding energy of -8.39 kcal/mol and Ki of 0.707 μM when compared to the standard inhibitor Bz-Nle-Lys-Arg-Arg-H that has been experimentally tested and shown to exhibit Ki value of 5.8 μM. Several modes of weak interactions were observed between the cyclopentapeptide CKRKC and the active site of DEN-2 NS3-NS2B protease. Thus, the cyclopentapeptide is proposed as a potential inhibitor to the NS3-NS2B protease activities of DEN-2. While these preliminary results are promising, further experimental investigation is necessary to validate the results.     

Functional characterization of a slow and tight-binding inhibitor of plasmin isolated from Russell's viper venom

Background

Snake venoms are rich in Kunitz-type protease inhibitors that may have therapeutic applications. However, apart from trypsin or chymotrypsin inhibition, the functions of most of these inhibitors have not been elucidated. A detailed functional characterization of these inhibitors may lead to valuable drug candidates.

Methods

A Kunitz-type protease inhibitor, named DrKIn-II, was tested for its ability to inhibit plasmin using various approaches such as far western blotting, kinetic analyses, fibrin plate assay and euglobulin clot lysis assay. In addition, the antifibrinolytic activity of DrKIn-II was demonstrated in vivo.

Results

DrKIn-II potently decreased the amidolytic activity of plasmin in a dose-dependent manner, with a global inhibition constant of 0.2 nM. Inhibition kinetics demonstrated that the initial binding of DrKIn-II causes the enzyme to isomerize, leading to the formation of a much tighter enzyme-inhibitor complex. DrKIn-II also demonstrated antifibrinolytic activity in fibrin plate assay and significantly prolonged the lysis of the euglobulin clot. Screening of DrKIn-II against a panel of serine proteases indicated that plasmin is the preferential target of DrKIn-II. Furthermore, DrKIn-II treatment prevented the increase of FDP in coagulation-stimulated mice and significantly reduced the bleeding time in a murine tail bleeding model.

Conclusion

DrKIn-II is a potent, slow and tight-binding plasmin inhibitor that demonstrates antifibrinolytic activity both in vitro and in vivo.

General significance

This is the first in-depth functional characterization of a plasmin inhibitor from a viperid snake. The potent antifibrinolytic activity of DrKIn-II makes it a potential candidate for the development of novel antifibrinolytic agents.     

The role of factor XIa (FXIa) catalytic domain exosite residues in substrate catalysis and inhibition by the Kunitz protease inhibitor domain of protease nexin 2

To select residues in coagulation factor XIa (FXIa) potentially important for substrate and inhibitor interactions, we examined the crystal structure of the complex between the catalytic domain of FXIa and the Kunitz protease inhibitor (KPI) domain of a physiologically relevant FXIa inhibitor, protease nexin 2 (PN2). Six FXIa catalytic domain residues (Glu(98), Tyr(143), Ile(151), Arg(3704), Lys(192), and Tyr(5901)) were subjected to mutational analysis to investigate the molecular interactions between FXIa and the small synthetic substrate (S-2366), the macromolecular substrate (factor IX (FIX)) and inhibitor PN2KPI. Analysis of all six Ala mutants demonstrated normal K(m) values for S-2366 hydrolysis, indicating normal substrate binding compared with plasma FXIa; however, all except E98A and K192A had impaired values of k(cat) for S-2366 hydrolysis. All six Ala mutants displayed deficient k(cat) values for FIX hydrolysis, and all were inhibited by PN2KPI with normal values of K(i) except for K192A, and Y5901A, which displayed increased values of K(i). The integrity of the S1 binding site residue, Asp(189), utilizing p-aminobenzamidine, was intact for all FXIa mutants. Thus, whereas all six residues are essential for catalysis of the macromolecular substrate (FIX), only four (Tyr(143), Ile(151), Arg(3704), and Tyr(5901)) are important for S-2366 hydrolysis; Glu(98) and Lys(192) are essential for FIX but not S-2366 hydrolysis; and Lys(192) and Tyr(5901) are required for both inhibitor and macromolecular substrate interactions.     

Properties of a trypsin and chymotrypsin inhibitor secreted by larval Taenia pisiformis

Living metacestodes of Taenia pisiformis maintained in vitro discharge into the surrounding medium a protease inhibitor, which has been purified from the medium by affinity chromatography on bovine α-chymotrypsin immobilized to CNBr-activated Sepharose 4B. The purified inhibitor was shown to inactivate the hydrolysis of $\text{N-α-}\text{benzoyl}\text{-}\text{L}\text{-}\text{arginine}$ ethyl ester and $\text{N-}\text{benzoyl}\text{-}\text{L}\text{-}\text{tyrosine}$ ethyl ester, respectively, by trypsin and chymotrypsin of bovine, rabbit and dog origin, and also the hydrolysis of casein by both bovine trypsin and bovine α and β chymotrypsins, but it did not affect the enzymic activity of subtilisin, elastase, collagenase, pepsin, rennin and papain. The inhibitor withstood heating at 100°C for up to 30 min, was stable in the pH range of 1.5–8.0, was unaffected by 8 M-urea or 0.2 M-2-mercaptoethanol, and had a molecular weight of about 7000 as calculated from its gel chromatographic behaviour. The inhibitor specifically inhibits either trypsin or chymotrypsin with the formation of stable enzyme inhibitor complexes that are not dissociated by 4 M-KCl. Inhibition of trypsin and chymotrypsin is non-competitive and is linear with inhibitor concentration up to 70–80% inhibition. Inhibitory activities toward both enzymes are functions of the same binding site of the inhibitor molecule. Complex formation between the inhibitor and the enzymes is timedependent; it requires 3–4 min for completion.     

Discovery of octahydropyrrolo [3,2-b] pyridin derivative as a highly selective Type I inhibitor of FGFR3 over VEGFR2 by high-throughput virtual screening

Although the aberrant activity of fibroblast growth factor receptor 3 (FGFR3) is implicated in various cancers, the reported kinase inhibitors of FGFR3 tend to cause side effects resulting from the inhibitory activity on vascular endothelial growth factor receptor 2 (VEGFR2). Therefore, it is necessary to find a novel high-selective inhibitor of FGFR3 over VEGFR2 from the small-molecule compound database. In this study, integrated virtual screening protocols were established to screen for selective inhibitors of FGFR3 over VEGFR2 in Drugbank and Asinex databases by combining three-dimensional pharmacophore model, molecular docking, molecular dynamics (MD) simulation, and molecular mechanics Poisson–Boltzmann surface area (MMPBSA) calculations. Finally, it is found that Asinex-5082, as an octahydropyrrolo[3,2-b] pyridin derivative, has larger binding free energy with FGFR3 (−39.3 kcal/mol) than reference drug Erdafitinib (−29.9 kcal/mol), while cannot bind with VEGFR2, resulting in considerable inhibitory selectivity. This is because Asinex-5082, unlike Erdafitinib, has not m-dimethoxybenzene with large steric hindrance, thus can enter the larger ATP-binding pocket of FGFR3 with DFG-in conformation to form hydrophobic interaction with residues Met529, Ile539, and Tyr557 as well as hydrogen bond with Ala558. On the other hand, due to the fact that the benzodioxane and N-heterocyclic rings are connected by carbonyl (C=O), Asinex-5082 cannot rotate freely so as to enter the smaller ATP binding pocket of VEGFR2 on the DFG-out conformation. The lead molecule Asinex-5082 may facilitate the rational design and development of novel selective inhibitors of FGFR3 over VEGFR2 as anticancer drugs.     

Characterisation of the α1-protease inhibitor system in Thoroughbred horse plasma by horizontal two-dimensional (ISO-DALT) electrophoresis. 2. Protease inhibition

The protease inhibitory spectra of the eight homozygous Thoroughbred Pi types against trypsin, elastase and chymotrypsin have been determined. The α1-protease inhibitor proteins exhibit three classes of inhibitory specificity towards these enzymes. The Pi types F, I, N and U exhibit class I (trypsin, elastase and chymotrypsin) and class II (trypsin and elastase) types of inhibition and fit Juneja et al.s (1979) classification of two separate genetic systems Pi 1 and Pi 2 based on differences in the inhibitory spectra against trypsin and chymotrypsin. The remaining four Pi types are exceptions to Juneja et al.s (1979) classification. Types G, L, S₁ and S₂ possess class I but not class II proteins. A third class of proteins (class III) which exclusively inhibit chymotrypsin was detected in all eight protease inhibitor types. Type G is well represented by class III proteins because two of the three major proteins of the ISO-DALT pattern inhibit only chymotrypsin and is thus an exception to Juneja et al.s (1979) classification.     

Prediction of the substrate for nonribosomal peptide synthetase (NRPS) adenylation domains by virtual screening

Nonribosomal peptide synthetases (NRPSs) synthesize a diverse array of bioactive small peptides, many of which are used in medicine. There is considerable interest in predicting NRPS substrate specificity in order to facilitate investigation of the many “cryptic” NRPS genes that have not been linked to any known product. However, the current sequence similarity‐based methods are unable to produce reliable predictions when there is a lack of prior specificity data, which is a particular problem for fungal NRPSs. We conducted virtual screening on the specificity‐determining domain of NRPSs, the adenylation domain, and found that virtual screening using experimentally determined structures results in good enrichment of the cognate substrate. Our results indicate that the conformation of the adenylation domain and in particular the conformation of a key conserved aromatic residue is important in determining the success of the virtual screening. When homology models of NRPS adenylation domains of known specificity, rather than experimentally determined structures, were built and used for virtual screening, good enrichment of the cognate substrate was also achieved in many cases. However, the accuracy of the models was key to the reliability of the predictions and there was a large variation in the results when different models of the same domain were used. This virtual screening approach is promising and is able to produce enrichment of the cognate substrates in many cases, but improvements in building and assessing homology models are required before the approach can be reliably applied to these models. Proteins 2015; 83:2052–2066. © 2015 Wiley Periodicals, Inc.