Structural determinants of trypsin affinity and specificity for cationic inhibitors

Primary structural features, physicochemical and biological properties of two water-soluble polysaccharides extracted from the brown Tunisian seaweed Halopteris scoparia

Marine algae are the most abundant resource in the marine environment and are still a promising source of bioactive compounds including hydrocolloids. This study contributes to the evaluation of the biological and biotechnological potentials of two water soluble polysaccharides, namely alginates (AHS) and fucoidan (FHS), extracted and purified from Halopteris scoparia, an abundant Tunisian brown macroalgae collected in Tunisia (Tabarka region). The total sugars, neutral monosaccharides, uronic acids, proteins, polyphenols, and sulfate groups contents were quantified for both fractions, as well as their functional groups and primary structural features by Fourier transform infrared spectroscopy, ionic and/or gas chromatography and nuclear magnetic resonance analyses. AHS and FHS showed significant anti-inflammatory (IC50 ≈ 1 mg/mL), anticoagulant (e.g., 27-61.7 for the activated partial thromboplastin time), antihyperglycemic (0.1-40 μg/mL) and anti-trypsin (IC50 ≈ 0.3-0.4 mg/mL) effects. FHS and a hydrolyzed fraction showed a very promising potential against herpes viruses (HSV-1) (IC50 < 28 μg/mL). Besides, AHS and two hydrolyzed fractions were able to stimulate the natural defenses of tomato seedlings, assessing their elicitor capacity, by increasing the activity of phenylalanine ammonia-lyase (66-422 %) but also significantly changing the polyphenol content in the leaves (121-243 %) and roots (30-104 %) of tomato plants.

Interactive molecular dynamics in virtual reality for accurate flexible protein-ligand docking

Simulating drug binding and unbinding is a challenge, as the rugged energy landscapes that separate bound and unbound states require extensive sampling that consumes significant computational resources. Here, we describe the use of interactive molecular dynamics in virtual reality (iMD-VR) as an accurate low-cost strategy for flexible protein-ligand docking. We outline an experimental protocol which enables expert iMD-VR users to guide ligands into and out of the binding pockets of trypsin, neuraminidase, and HIV-1 protease, and recreate their respective crystallographic protein-ligand binding poses within 5-10 minutes. Following a brief training phase, our studies shown that iMD-VR novices were able to generate unbinding and rebinding pathways on similar timescales as iMD-VR experts, with the majority able to recover binding poses within 2.15 Å RMSD of the crystallographic binding pose. These results indicate that iMD-VR affords sufficient control for users to carry out the detailed atomic manipulations required to dock flexible ligands into dynamic enzyme active sites and recover crystallographic poses, offering an interesting new approach for simulating drug docking and generating binding hypotheses.

Influential Mutations in the SMAD4 Trimer Complex Can Be Detected from Disruptions of Electrostatic Complementarity

This article examines three techniques for rapidly assessing the electrostatic contribution of individual amino acids to the stability of protein-protein complexes. Whereas the energetic minimization of modeled oligomers may yield more accurate complexes, we examined the possibility that simple modeling may be sufficient to identify amino acids that add to or detract from electrostatic complementarity. The three methods evaluated were (a) the elimination of entire side chains (e.g., glycine scanning), (b) the elimination of the electrostatic contribution from the atoms of a side chain, called nullification, and (c) side chain structure prediction using SCWRL4. These techniques generate models in seconds, enabling large-scale mutational scanning. We evaluated these techniques on the SMAD2/SMAD4 heterotrimer, whose formation plays a crucial role in antitumor pathways. Many studies have documented the clinical and structural effect of specific mutations on trimer formation. Our results describe how glycine scanning yields more specific predictions, although nullification may be more sensitive, and how side chain structure prediction enables the identification of uncharged-to-charge mutations.

Different disulfide bridge connectivity drives alternative folds in highly homologous Brassicaceae trypsin inhibitors

Low-molecular-mass trypsin inhibitors from Arabidopsis thaliana, Brassica napus var. oleifera, and Sinapis alba L. (ATTI, RTI, and MTI, respectively) display more than 69% amino acid sequence identity. Among others, the amino acid sequence Cys-Ala-Pro-Arg-Ile building up the inhibitor reactive site, and the eight Cys residues forming four disulfide bridges are conserved. However, the disulfide bridge connectivity of RTI and MTI (C1-C3, C2-C4, C5-C6, and C7-C8) is different from that of ATTI Cys (C1-C8, C2-C5, C3-C6, and C4-C7). Despite the different disulfide bridge connectivity, the reactive site loop of ATTI, RTI, and MTI is solvent exposed permitting trypsin recognition. Structural considerations here reported suggest that proteins showing high amino acid sequence identity and common functional properties could display different three-dimensional structures. This may reflect high inhibitor plasticity in relation to plant-pathogen interactions, plant tissue development as well as the different redox potential of cell compartments.

VASP-E: specificity annotation with a volumetric analysis of electrostatic isopotentials

Algorithms for comparing protein structure are frequently used for function annotation. By searching for subtle similarities among very different proteins, these algorithms can identify remote homologs with similar biological functions. In contrast, few comparison algorithms focus on specificity annotation, where the identification of subtle differences among very similar proteins can assist in finding small structural variations that create differences in binding specificity. Few specificity annotation methods consider electrostatic fields, which play a critical role in molecular recognition. To fill this gap, this paper describes VASP-E (Volumetric Analysis of Surface Properties with Electrostatics), a novel volumetric comparison tool based on the electrostatic comparison of protein-ligand and protein-protein binding sites. VASP-E exploits the central observation that three dimensional solids can be used to fully represent and compare both electrostatic isopotentials and molecular surfaces. With this integrated representation, VASP-E is able to dissect the electrostatic environments of protein-ligand and protein-protein binding interfaces, identifying individual amino acids that have an electrostatic influence on binding specificity. VASP-E was used to examine a nonredundant subset of the serine and cysteine proteases as well as the barnase-barstar and Rap1a-raf complexes. Based on amino acids established by various experimental studies to have an electrostatic influence on binding specificity, VASP-E identified electrostatically influential amino acids with 100% precision and 83.3% recall. We also show that VASP-E can accurately classify closely related ligand binding cavities into groups with different binding preferences. These results suggest that VASP-E should prove a useful tool for the characterization of specific binding and the engineering of binding preferences in proteins.

Proteins, the ubiquitous worker molecules of the cell, are a diverse class of molecules that perform very specific tasks. Understanding how proteins achieve specificity is a critical step towards understanding biological systems and a key prerequisite for rationally engineering new proteins. To examine electrostatic influences on specificity in proteins, this paper presents VASP-E, a software tool that generates solid representations of the electrostatic potential fields that surround proteins. VASP-E compares solids with constructive solid geometry, a class of techniques developed first for modeling complex machine parts. We observed that solid representations could quantify the degree of charge complementarity in protein-protein interactions and identify key residues that strengthen or weaken them. VASP-E correctly identified amino acids with established experimental influences on protein-protein binding specificity. We also observed that solid representations of electrostatic fields could identify electrostatic conservations and variations that relate to similarities and differences in binding specificity between proteins and small molecules.

Isolating influential regions of electrostatic focusing in protein and DNA structure

Electrostatic focusing is a general phenomenon that occurs in cavities and grooves on the molecular surface of biomolecules. Narrow surface features can partially shield charged atoms from the high-dielectric solvent, enhancing electrostatic potentials inside the cavity and projecting electric field lines outward into the solvent. This effect has been observed in many instances and is widely considered in the human examination of molecular structure, but it is rarely integrated into the digital representations used in protein structure comparison software. To create a computational representation of electrostatic focusing, that is compatible with structure comparison algorithms, this paper presents an approach that generates three-dimensional solids that approximate regions where focusing occurs. We verify the accuracy of this representation against instances of focusing in proteins and DNA. Noting that this representation also identifies thin focusing regions on the molecular surface that are unlikely to affect binding, we describe a second algorithm that conservatively isolates larger focusing regions. The resulting 3D solids can be compared with Boolean set operations, permitting a new range of analyses on the regions where electrostatic focusing occurs. They also represent a novel integration of molecular shape and electrostatic focusing into the same structure comparison framework.

On the importance of polar interactions for complexes containing intrinsically disordered proteins

There is a growing recognition for the importance of proteins with large intrinsically disordered (ID) segments in cell signaling and regulation. ID segments in these proteins often harbor regions that mediate molecular recognition. Coupled folding and binding of the recognition regions has been proposed to confer high specificity to interactions involving ID segments. However, researchers recently questioned the origin of the interaction specificity of ID proteins because of the overrepresentation of hydrophobic residues in their interaction interfaces. Here, we focused on the role of polar and charged residues in interactions mediated by ID segments. Making use of the extended nature of most ID segments when in complex with globular proteins, we first identified large numbers of complexes between globular proteins and ID segments by using radius-of-gyration-based selection criteria. Consistent with previous studies, we found the interfaces of these complexes to be enriched in hydrophobic residues, and that these residues contribute significantly to the stability of the interaction interface. However, our analyses also show that polar interactions play a larger role in these complexes than in structured protein complexes. Computational alanine scanning and salt-bridge analysis indicate that interfaces in ID complexes are highly complementary with respect to electrostatics, more so than interfaces of globular proteins. Follow-up calculations of the electrostatic contributions to the free energy of binding uncovered significantly stronger Coulombic interactions in complexes harbouring ID segments than in structured protein complexes. However, they are counter-balanced by even higher polar-desolvation penalties. We propose that polar interactions are a key contributing factor to the observed high specificity of ID segment-mediated interactions.

Protein-protein interactions are essential to communication and signal integration in cells. For these processes to be precise, interactions between proteins have to be specific and well coordinated. In order to understand the specificity in protein interactions, researches have focused on interfaces between two or more folded proteins. It has been shown that specificity in interactions between folded proteins relies on shape complementarity, hydrogen bonding, and salt-bridge formation. However, many proteins lack a unique folded structure; the so-called intrinsically disordered proteins. These proteins fluctuate between different conformations in isolation but often adopt a single structure when interacting with partner proteins. As many intrinsically disordered proteins are involved in signaling and regulation, their interactions have to be highly specific. The finding that the interaction interfaces of intrinsically disordered proteins are enriched in hydrophobic residues has led to questions regarding the specificity of interactions mediated by this group of proteins. Here, we show that polar and charged residues play a larger role in interfaces that involve intrinsically disordered proteins compared to interfaces that involve only folded proteins. Our results suggest that polar interactions are key contributors to the specificity of interactions that involve intrinsically disordered proteins.

Possible conformational change within the desolvated and cationized sBBI/trypsin non-covalent complex during the collision-induced dissociation process

Electrospray ionization mass spectrometry (ESI-MS) has become an analytical technique widely used for the investigation of non-covalent protein-protein and protein-ligand complexes due to the soft desolvation conditions that preserve the stoichiometry of the interacting partners. Dissociation studies of solvated or desolvated complexes (in the source and in the collision cell, respectively) allow access to information on protein conformation and localization of the metal ions involved in protein structure stabilization and biological activity. The complex of bovine trypsin and small soybean Bowman-Birk inhibitor (sBBI) was studied by ESI-MS to determine changes occurring within the complex during its transfer from droplets to the gas phase independently of the ion polarity. Under collision-induced dissociation (CID) conditions, unexpected binding of the Ca(2+) ion (cofactor of native trypsin) to the inhibitor molecule was observed within the desolvated sBBI/trypsin/Ca(2+) complex (with a 1:1:1 stoichiometry). This formal gas-phase migration of the calcium ion from trypsin to the inhibitor may be related to conformational rearrangements in the solvent-free and likely collapsed complex. However, under conditions leading to the increase in complex charge state, the appearance of the cationized trypsin molecule was detected during complex dissociation, thus reflecting different pathways of the evolution of complex conformation.

Spermine oxidase (SMO) activity in breast tumor tissues and biochemical analysis of the anticancer spermine analogues BENSpm and CPENSpm

BACKGROUND

Polyamine metabolism has a critical role in cell death and proliferation representing a potential target for intervention in breast cancer (BC). This study investigates the expression of spermine oxidase (SMO) and its prognostic significance in BC. Biochemical analysis of Spm analogues BENSpm and CPENSpm, utilized in anticancer therapy, was also carried out to test their property in silico and in vitro on the recombinant SMO enzyme.

METHODS

BC tissue samples were analyzed for SMO transcript level and SMO activity. Student's t test was applied to evaluate the significance of the differences in value observed in T and NT samples. The structure modeling analysis of BENSpm and CPENSpm complexes formed with the SMO enzyme and their inhibitory activity, assayed by in vitro experiments, were examined.

RESULTS

Both the expression level of SMO mRNA and SMO enzyme activity were significantly lower in BC samples compared to NT samples. The modeling of BENSpm and CPENSpm complexes formed with SMO and their inhibition properties showed that both were good inhibitors.

CONCLUSIONS

This study shows that underexpression of SMO is a negative marker in BC. The SMO induction is a remarkable chemotherapeutical target. The BENSpm and CPENSpm are efficient SMO inhibitors. The inhibition properties shown by these analogues could explain their poor positive outcomes in Phases I and II of clinical trials.

Thrombin inhibitors with novel P1 binding pocket functionality: free energy of binding analysis

The high incidence of thrombembolic diseases justifies the development of new antithrombotics. The search for a direct inhibitor has resulted in the synthesis of a considerable number of low molecular weight molecules that inhibit human alpha-thrombin potently. However, efforts to develop an orally active drug remain in progress as the most active inhibitors with a highly basic P1 moiety exhibit an unsatisfactory bioavailability profile. In our previous work we solved several X-ray structures of human alpha-thrombin in complexes with (1) novel bicyclic arginine mimetics attached to the glycylproline amide and pyridinone acetamide scaffold and (2) inhibitors with a novel aza scaffold and with charged or neutral P1 moieties. In the present contribution, we correlate the structures of the complex between these inhibitors and the protein with the calculated free energy of binding. The energy of solvation was calculated using the Poisson-Boltzmann approach. In particular, the requirements for successful recognition of an inhibitor at the protein's active site pocket S1 are discussed.

Functional divergence in tandemly duplicated Arabidopsis thaliana trypsin inhibitor genes

In multigene families, variation among loci and alleles can contribute to trait evolution. We explored patterns of functional and genetic variation in six duplicated Arabidopsis thaliana trypsin inhibitor (ATTI) loci. We demonstrate significant variation in constitutive and herbivore-induced transcription among ATTI loci that show, on average, 65% sequence divergence. Significant variation in ATTI expression was also found between two molecularly defined haplotype classes. Population genetic analyses for 17 accessions of A. thaliana showed that six ATTI loci arranged in tandem within 10 kb varied 10-fold in nucleotide diversity, from 0.0009 to 0.0110, and identified a minimum of six recombination events throughout the tandem array. We observed a significant peak in nucleotide and indel polymorphism spanning ATTI loci in the interior of the array, due primarily to divergence between the two haplotype classes. Significant deviation from the neutral equilibrium model for individual genes was interpreted within the context of intergene linkage disequilibrium and correlated patterns of functional differentiation. In contrast to the outcrosser Arabidopsis lyrata for which recombination is observed even within ATTI loci, our data suggest that response to selection was slowed in the inbreeding, annual A. thaliana because of interference among functionally divergent ATTI loci.

Analysing six types of protein-protein interfaces

Non-covalent residue side-chain interactions occur in many different types of proteins and facilitate many biological functions. Are these differences manifested in the sequence compositions and/or the residue-residue contact preferences of the interfaces? Previous studies analysed small data sets and gave contradictory answers. Here, we introduced a new data-mining method that yielded the largest high-resolution data set of interactions analysed. We introduced an information theory-based analysis method. On the basis of sequence features, we were able to differentiate six types of protein interfaces, each corresponding to a different functional or structural association between residues. Particularly, we found significant differences in amino acid composition and residue-residue preferences between interactions of residues within the same structural domain and between different domains, between permanent and transient interfaces, and between interactions associating homo-oligomers and hetero-oligomers. The differences between the six types were so substantial that, using amino acid composition alone, we could predict statistically to which of the six types of interfaces a pool of 1000 residues belongs at 63-100% accuracy. All interfaces differed significantly from the background of all residues in SWISS-PROT, from the group of surface residues, and from internal residues that were not involved in non-trivial interactions. Overall, our results suggest that the interface type could be predicted from sequence and that interface-type specific mean-field potentials may be adequate for certain applications.

Can the calculation of ligand binding free energies be improved with continuum solvent electrostatics and an ideal-gas entropy correction?

The prediction of a ligand binding constant requires generating three-dimensional structures of the complex concerned and reliably scoring these structures. Here, the scoring problem is investigated by examining benzamidine-like inhibitors of trypsin, a system for which errors in the structures are small. Precise and consistent binding free energies for the inhibitors are determined experimentally for this test system. To examine possible improvement of scoring methods, we test the suitability of continuum electrostatics to account for solvation effects and use an ideal-gas entropy correction to account for the changes in the degrees of freedom of the ligand. The small observed root-mean-square deviation of 0.55 kcal/mol of the calculated relative to the experimental values indicates that the essentials of the binding process have been captured. Even though all six ligands make the same salt bridge and H-bonds to the protein, the electrostatic contribution varies among the ligands by as much as 2 kcal/mol. Moreover, although the ligands are rigid and similar in size, the entropic terms also significantly affect the relative binding affinities (by up to 2.7 kcal/mol). The present approach to solvation and entropy may allow the ranking of the ligands to be considerably improved at a cost that makes the method applicable to the optimization of lead compounds or to the screening of small collections of ligands.

Flexibility in monomeric Cu,Zn superoxide dismutase detected by limited proteolysis and molecular dynamics simulation

Limited proteolysis by trypsin of monomeric Cu,Zn superoxide dismutase from Escherichia coli induces a specific cleavage of the polypeptide chain at the level of Lys60 located in the S-S subloop of loop 6,5 where, when compared to the eukaryotic enzyme, a seven-residues insertion, completely exposed to the solvent, is observed. This result suggests that this subloop is disordered and flexible, thus enabling binding and adaptation to the active site of the proteolytic enzyme. Indeed, molecular dynamics simulation indicates that the S-S subloop undergoes high fluctuations and that its high flexibility coupled to an high solvent accessibility can explain the specific bond selection of the protease. As a matter of fact, of the possible 14 solvent accessible proteolytic sites only the Lys60 flexible site is cleaved. High flexibility and solvent exposure are confirmed by the short water residence time for the residues corresponding to the cleavage site evaluated by molecular dynamics simulation. These experiments demonstrate that molecular dynamics simulation and limited proteolysis are complementary and unambiguous tools to identify flexible sites in proteins.

Contributions to the binding free energy of ligands to avidin and streptavidin

The free energy of binding of a ligand to a macromolecule is here formally decomposed into the (effective) energy of interaction, reorganization energy of the ligand and the macromolecule, conformational entropy change of the ligand and the macromolecule, and translational and rotational entropy loss of the ligand. Molecular dynamics simulations with implicit solvation are used to evaluate these contributions in the binding of biotin, biotin analogs, and two peptides to avidin and streptavidin. We find that the largest contribution opposing binding is the protein reorganization energy, which is calculated to be from 10 to 30 kcal/mol for the ligands considered here. The ligand reorganization energy is also significant for flexible ligands. The translational/rotational entropy is 4.5-6 kcal/mol at 1 M standard state and room temperature. The calculated binding free energies are in the correct range, but the large statistical uncertainty in the protein reorganization energy precludes precise predictions. For some complexes, the simulations show multiple binding modes, different from the one observed in the crystal structure. This finding is probably due to deficiencies in the force field but may also reflect considerable ligand flexibility.

Regulation of serine-type exoproteinases by endogenous inhibitors present in exoantigens of the mycelial form of Paracoccidioides brasiliensis

We have partially characterized some biochemical properties of exoproteinases secreted into culture medium by the mycelial form of Paracoccidioides brasiliensis, a dimorphic fungus that causes human disease in Latin America. Proteinase activity was analyzed in solid- and liquid-phase systems using zymography and Azocoll, respectively. Minimal or no gelatinase activity was observed by zymography in the crude filtrates among proteins with a relative mobility greater than 200 kDa. When the crude filtrate was fractionated by isoelectric focusing or ion exchange chromatography, we observed striking activation of gelatinases, both those of high apparent molecular mass and alkaline isoelectric points (pI), as well as those of lower molecular mass and acidic pI. The apparent high molecular mass gelatinases, pI 10, showed optimal activity at pH 7.0. They were totally inhibited by phenylmethylsulfonylfluoride and partially inhibited by incubation with previously neutralized fractions of pI 5.4 and 6.1. The latter inhibition could be reversed by exposure to 10% isopropanol. These results provide evidence of regulatory mechanisms controlling proteinase activity in secreted proteins. The principal mechanism appears to be the formation of reversible complexes with endogenous inhibitors.

Free energy decomposition of protein-protein interactions

A free energy decomposition scheme has been developed and tested on antibody-antigen and protease-inhibitor binding for which accurate experimental structures were available for both free and bound proteins. Using the x-ray coordinates of the free and bound proteins, the absolute binding free energy was computed assuming additivity of three well-defined, physical processes: desolvation of the x-ray structures, isomerization of the x-ray conformation to a nearby local minimum in the gas-phase, and subsequent noncovalent complex formation in the gas phase. This free energy scheme, together with the Generalized Born model for computing the electrostatic solvation free energy, yielded binding free energies in remarkable agreement with experimental data. Two assumptions commonly used in theoretical treatments; viz., the rigid-binding approximation (which assumes no conformational change upon complexation) and the neglect of vdW interactions, were found to yield large errors in the binding free energy. Protein-protein vdW and electrostatic interactions between complementary surfaces over a relatively large area (1400--1700 A(2)) were found to drive antibody-antigen and protease-inhibitor binding.

Statistical potentials and scoring functions applied to protein-ligand binding

In virtual screening, small-molecule ligands are docked into protein binding sites and their binding affinity is predicted. Knowledge-based, regression-based and first-principle-based methods have been developed to rank computer-generated binding modes. As a result of still existing deficiencies, a best compromise might be the combination of several scoring schemes into a consensus scoring approach.

Serine proteinase inhibition by the active site titrant N alpha-(N, N-dimethylcarbamoyl)-alpha-azaornithine p-nitrophenyl ester. A comparative study

Kinetics for the hydrolysis of the chromogenic active-site titrant N alpha-(N,N-dimethylcarbamoyl)-alpha-azaornithine p-nitrophenyl ester (Dmc-azaOrn-ONp) catalysed by bovine beta-trypsin, bovine alpha-thrombin, bovine Factor Xa, human alpha-thrombin, human Factor Xa, human Lys77-plasmin, human urinary kallikrein, Mr 33 000 and Mr 54 000 species of human urokinase, porcine pancreatic beta-kallikrein-A and -B and Ancrod (the coagulating serine proteinase from the Malayan pit viper Agkistrodon rhodostoma venom) have been obtained between pH 6.0 and 8.0, at 21.0 degrees C, and analysed in parallel with those for the enzymatic cleavage of N alpha-(N,N-dimethylcarbamoyl)-alpha-azalysine p-nitrophenyl ester (Dmc-azaLys-ONp). The enzyme kinetics are consistent with the minimum three-step catalytic mechanism of serine proteinases, the rate-limiting step being represented by the deacylation process. Bovine beta-trypsin kinetics are modulated by the acid-base equilibrium of the His57 catalytic residue (pKa approximately 6.9). Dmc-azaOrn-ONp and Dmc-azaLys-ONp bind stoichiometrically to the serine proteinase active site, and allow the reliable determination of the active enzyme concentration between 1.0 x 10-6 M and 3.0 x 10-4 M. The affinity and the reactivity for Dmc-azaOrn-ONp (expressed by Ks and k+2/Ks, respectively) of the serine proteinases considered are much lower than those for Dmc-azaLys-ONp. The very different affinity and reactivity properties for Dmc-azaOrn-ONp and Dmc-azaLys-ONp have been related to the different size of the ornithine/lysine side chains, and to the ensuing different positioning of the active-site titrants upon binding to the enzyme catalytic centre (i.e. to P1-S1 recognition). These data represent the first detailed comparative investigation on the catalytic properties of serine proteinases towards an ornithine derivative (i. e. Dmc-azaOrn-ONp).