Cation-pi interactions at the active site of factor Xa: dramatic enhancement upon stepwise N-alkylation of ammonium ions

Trimethyllysine Reader Proteins Exhibit Widespread Charge-Agnostic Binding via Different Mechanisms to Cationic and Neutral Ligands

In the last 40 years, cation-π interactions have become part of the lexicon of noncovalent forces that drive protein binding. Indeed, tetraalkylammoniums are universally bound by aromatic cages in proteins, suggesting that cation-π interactions are a privileged mechanism for binding these ligands. A prominent example is the recognition of histone trimethyllysine (Kme3) by the conserved aromatic cage of reader proteins, dictating gene expression. However, two proteins have recently been suggested as possible exceptions to the conventional understanding of tetraalkylammonium recognition. To broadly interrogate the role of cation-π interactions in protein binding interactions, we report the first large-scale comparative evaluation of reader proteins for a neutral Kme3 isostere, experimental and computational mechanistic studies, and structural analysis. We find unexpected widespread binding of readers to a neutral isostere with the first examples of readers that bind the neutral isostere more tightly than Kme3. We find that no single factor dictates the charge selectivity, demonstrating the challenge of predicting such interactions. Further, readers that bind both cationic and neutral ligands differ in mechanism: binding Kme3 via cation-π interactions and the neutral isostere through the hydrophobic effect in the same aromatic cage. This discovery explains apparently contradictory results in previous studies, challenges traditional understanding of molecular recognition of tetraalkylammoniums by aromatic cages in myriad protein-ligand interactions, and establishes a new framework for selective inhibitor design by exploiting differences in charge dependence.

Systematic Study of Heteroarene Stacking Using a Congeneric Set of Molecular Glues for Procaspase-6

Heteroaromatic stacking interactions are important in drug binding, supramolecular chemistry, and materials science, making protein-ligand model systems of these interactions of considerable interest. Here we studied 30 congeneric ligands that each present a distinct heteroarene for stacking between tyrosine residues at the dimer interface of procaspase-6. Complex X-ray crystal structures of 10 analogs showed that stacking geometries were well conserved, while high-accuracy computations showed that heteroarene stacking energy was well correlated with predicted overall ligand binding energies. Empirically determined KD values in this system thus provide a useful measure of heteroarene stacking with tyrosine. Stacking energies are discussed in the context of torsional strain, the number and positioning of heteroatoms, tautomeric state, and coaxial orientation of heteroarene in the stack. Overall, this study provides an extensive data set of empirical and high-level computed binding energies in a versatile new protein-ligand system amenable to studies of other intermolecular interactions.

Investigation of the tight binding mechanism of a new anticoagulant DD217 to factor Xa by means of molecular docking and molecular dynamics

A new promising drug candidate DD217 has been proposed recently as a potent anticoagulant acting on factor Xa (fXa) target. It exhibits the lowest concentration of doubling the prothrombin time among the known anticoagulants. In order to explain the efficacy of DD217 in terms of molecular interactions with its target we studied the hypothesis of the tight binding mechanism by means of molecular dynamics simulations and statistical analysis of the trajectory. The conducted analysis confirms the significant contributions to the MM/GBSA estimated binding free energy of the S4 pocket residues as well the crucial role of establishing the hydrogen bonds between the ligand and the backbone amides of Gly216 and Gly218 of the target. The simulation results support the hypothesis of the tight binding mechanism of DD217 to fXa.Communicated by Ramaswamy H. Sarma.

A systematic approach for evaluating the role of surface-exposed loops in trypsin-like serine proteases applied to the 170 loop in coagulation factor VIIa

Proteases play a major role in many vital physiological processes. Trypsin-like serine proteases (TLPs), in particular, are paramount in proteolytic cascade systems such as blood coagulation and complement activation. The structural topology of TLPs is highly conserved, with the trypsin fold comprising two β-barrels connected by a number of variable surface-exposed loops that provide a surprising capacity for functional diversity and substrate specificity. To expand our understanding of the roles these loops play in substrate and co-factor interactions, we employ a systematic methodology akin to the natural truncations and insertions observed through evolution of TLPs. The approach explores a larger deletion space than classical random or directed mutagenesis. Using FVIIa as a model system, deletions of 1–7 amino acids through the surface exposed 170 loop, a vital allosteric regulator, was introduced. All variants were extensively evaluated by established functional assays and computational loop modelling with Rosetta. The approach revealed detailed structural and functional insights recapitulation and expanding on the main findings in relation to 170 loop functions elucidated over several decades using more cumbersome crystallization and single deletion/mutation methodologies. The larger deletion space was key in capturing the most active variant, which unexpectedly had a six-amino acid truncation. This variant would have remained undiscovered if only 2–3 deletions were considered, supporting the usefulness of the methodology in general protease engineering approaches. Our findings shed further light on the complex role that surface-exposed loops play in TLP function and supports the important role of loop length in the regulation and fine-tunning of enzymatic function throughout evolution.

Identification of Spiro-Fused Pyrrolo[3,4-a]pyrrolizines and Tryptanthrines as Potential Antitumor Agents: Synthesis and In Vitro Evaluation

A series of heterocyclic compounds containing a spiro-fused pyrrolo[3,4-a]pyrrolizine and tryptanthrin framework have been synthesized and studied as potential antitumor agents. Cytotoxicity of products was screened against human erythroleukemia (K562) and human cervical carcinoma (HeLa) cell lines. Among the screened compounds. 4a, 4b and 5a were active against human erythroleukemia (K562) cell line, while 4a and 5a were active against cervical carcinoma (HeLa) cell line. In agreement with the DNA cytometry studies, the tested compounds have achieved significant cell-cycle perturbation with higher accumulation of cells in G2/M phase and induced apoptosis. Using confocal microscopy, we found that with 4a and 5a treatment of HeLa cells, actin filaments disappeared, and granular actin was distributed diffusely in the cytoplasm in 76–91% of cells. We discovered that HeLa cells after treatment with compounds 4a and 5a significantly reduced the number of cells with filopodium-like membrane protrusions (from 63 % in control cells to 29% after treatment) and a decrease in cell motility.

Aromatic interactions in Glycochemistry: from molecular recognition to catalysis

Aromatic platforms are ubiquitous recognition motifs occurring in protein carbohydrate binding domains (CBDs), RNA receptors and enzymes. They stabilize the glycoside/receptor complexes by participating in stacking CH/π interactions with either the α- or β- face of the corresponding pyranose units. In addition, the role played by aromatic units in the stabilization of glycoside cationic transition states has started being recognized in recent years. Extensive studies carried out during the last decade have allowed to dissect the main contributing forces that stabilize the carbohydrate/aromatic complexes, while helping delineate not only the standing relationship between the glycoside/aromatic chemical structures and the strength of this interaction, but also their potential influence on glycoside reactivity.

Machine Learning-Based Scoring Functions, Development and Applications with SAnDReS

BACKGROUND

Analysis of atomic coordinates of protein-ligand complexes can provide three-dimensional data to generate computational models to evaluate binding affinity and thermodynamic state functions. Application of machine learning techniques can create models to assess protein-ligand potential energy and binding affinity. These methods show superior predictive performance when compared with classical scoring functions available in docking programs.

OBJECTIVE

Our purpose here is to review the development and application of the program SAnDReS. We describe the creation of machine learning models to assess the binding affinity of protein-ligand complexes.

METHODS

SAnDReS implements machine learning methods available in the scikit-learn library. This program is available for download at https://github.com/azevedolab/sandres. SAnDReS uses crystallographic structures, binding and thermodynamic data to create targeted scoring functions.

RESULTS

Recent applications of the program SAnDReS to drug targets such as Coagulation factor Xa, cyclin-dependent kinases and HIV-1 protease were able to create targeted scoring functions to predict inhibition of these proteins. These targeted models outperform classical scoring functions.

CONCLUSION

Here, we reviewed the development of machine learning scoring functions to predict binding affinity through the application of the program SAnDReS. Our studies show the superior predictive performance of the SAnDReS-developed models when compared with classical scoring functions available in the programs such as AutoDock4, Molegro Virtual Docker and AutoDock Vina.

Homology Modeling of the Human P-glycoprotein (ABCB1) and Insights into Ligand Binding through Molecular Docking Studies

The ABCB1 transporter also known as P-glycoprotein (P-gp) is a transmembrane protein belonging to the ATP binding cassette super-family of transporters; it is a xenobiotic efflux pump that limits intracellular drug accumulation by pumping the compounds out of cells. P-gp contributes to a decrease of toxicity and possesses broad substrate specificity. It is involved in the failure of numerous anticancer and antiviral chemotherapies due to the multidrug resistance (MDR) phenomenon, where it removes the chemotherapeutics out of the targeted cells. Understanding the details of the ligand–P-gp interaction is therefore crucial for the development of drugs that might overcome the MRD phenomenon and for obtaining a more effective prediction of the toxicity of certain compounds. In this work, an in silico modeling was performed using homology modeling and molecular docking methods with the aim of better understanding the ligand–P-gp interactions. Based on different mouse P-gp structural templates from the PDB repository, a 3D model of the human P-gp (hP-gp) was constructed by means of protein homology modeling. The homology model was then used to perform molecular docking calculations on a set of thirteen compounds, including some well-known compounds that interact with P-gp as substrates, inhibitors, or both. The sum of ranking differences (SRD) was employed for the comparison of the different scoring functions used in the docking calculations. A consensus-ranking scheme was employed for the selection of the top-ranked pose for each docked ligand. The docking results showed that a high number of π interactions, mainly π–sigma, π–alkyl, and π–π type of interactions, together with the simultaneous presence of hydrogen bond interactions contribute to the stability of the ligand–protein complex in the binding site. It was also observed that some interacting residues in hP-gp are the same when compared to those observed in a co-crystallized ligand (PBDE-100) with mouse P-gp (PDB ID: 4XWK). Our in silico approach is consistent with available experimental results regarding P-gp efflux transport assay; therefore it could be useful in the prediction of the role of new compounds in systemic toxicity.

Can glycine betaine denature proteins?

Glycine betaine (GB) is a naturally occurring osmolyte that has been widely recognized as a protein protectant. Since GB consists of a methylated ammonium moiety, it can engage in strong cation-π interactions with aromatic amino acid sidechains. We hypothesize that such specific binding interactions would allow GB to decrease the stability of proteins that are predominantly stabilized by a cluster of aromatic amino acids. To test this hypothesis, we investigate the effect of GB on the stability of two β-hairpins (or mini-proteins) that contain such a cluster. We find that for both systems the stability of the folded state first decreases and then increases with increasing GB concentration. Such non-monotonic dependence not only confirms that GB can act as a protein denaturant, but also underscores the complex interplay between GB's stabilizing and destabilizing forces toward a given protein. While stabilizing osmolytes all have the tendency to be excluded from the protein surface which is the action underlying their stabilizing effect, our results suggest that in order to quantitatively assess the effect of GB on the stability of any given protein, specific cation-π binding interactions need to be explicitly considered. Moreover, our results show, consistent with other studies, that cation methylation can strengthen the respective cation-π interactions. Taken together, these findings provide new insight into the mechanism by which amino acid-based osmolytes interact with proteins.

Urea-aromatic interactions in biology

Noncovalent interactions are key determinants in both chemical and biological processes. Among such processes, the hydrophobic interactions play an eminent role in folding of proteins, nucleic acids, formation of membranes, protein-ligand recognition, etc.. Though this interaction is mediated through the aqueous solvent, the stability of the above biomolecules can be highly sensitive to any small external perturbations, such as temperature, pressure, pH, or even cosolvent additives, like, urea-a highly soluble small organic molecule utilized by various living organisms to regulate osmotic pressure. A plethora of detailed studies exist covering both experimental and theoretical regimes, to understand how urea modulates the stability of biological macromolecules. While experimentalists have been primarily focusing on the thermodynamic and kinetic aspects, theoretical modeling predominantly involves mechanistic information at the molecular level, calculating atomistic details applying the force field approach to the high level electronic details using the quantum mechanical methods. The review focuses mainly on examples with biological relevance, such as (1) urea-assisted protein unfolding, (2) urea-assisted RNA unfolding, (3) urea lesion interaction within damaged DNA, (4) urea conduction through membrane proteins, and (5) protein-ligand interactions those explicitly address the vitality of hydrophobic interactions involving exclusively the urea-aromatic moiety.

Characterization of the enzymatic activity of the serine protease domain of Factor VII activating protease (FSAP)

Factor VII (FVII) activating protease (FSAP) is a circulating serine protease. Human genetic studies, based on the Marburg I (MI) (Gly221Glu, chymotrypsin numbering system) polymorphism, implicate FSAP in the pathogenesis of many diseases. Here, we describe the molecular and functional changes caused by the Gly221Glu substitution in the 220 loop using recombinant proteins expressed in E. coli. The serine protease domain (SPD) of wild type (WT) FSAP displayed auto-catalytic activation whereas the MI isoform displayed very low autocatalytic activation and low proteolytic activity against the chromogenic substrate S-2288, Factor VII, tissue factor pathway inhibitor as well as pro-urokinase. Introduction of a thermolysin cleavage site in the activation position (Arg15Gln) led to cleavage of both WT- and MI-SPD and the resulting WT-SPD, but not the MI-SPD, was active. Mutating the Gly221 position to Asp, Gln and Leu led to a loss of activity whereas the Ala substitution was partially active. These results suggest a disturbance of the active site, or non-accessibility of the substrate to the active site in MI-SPD. With respect to regulation with metal ions, calcium, more than sodium, increased the enzymatic activity of WT-SPD. Thus, we describe a novel method for the production of recombinant FSAP-SPD to understand the role of the MI-single nucleotide polymorphism (SNP) in the regulation of its activity.

An Empirical Study of Amide-Heteroarene π-Stacking Interactions Using Reversible Inhibitors of a Bacterial Serine Hydrolase

Compared to aryl-aryl π-stacking interactions, the analogous stacking of heteroarenes on amide π systems is less well understood and vastly underutilized in structure-based drug design. Recent theoretical studies have delineated the important geometric coordinates of the interaction, some of which have been confirmed with synthetic model systems based on Rebek imides. Unfortunately, a broadly useful and tractable protein-ligand model system of this interaction has remained elusive. Here we employed a known inhibitor scaffold to study π-stacking of diverse heteroarene substituents on the amide face of Gly238 in the cephalosporinases CTX-M-14 and CTX-M-27. Biochemical inhibition constants (K _i) and biophysical binding constants (K _d) were determined for nineteen new analogues against both enzymes, while multiple high-resolution co-crystal structures revealed remarkably consistent placement of the probe heteroarene on Gly238. The data presented support the predicted importance of opposing dipoles in amide-heteroarene interactions and should be useful for evaluating other theoretical predictions concerning these interactions.

Sodium-site in serine protease domain of human coagulation factor IXa: evidence from the crystal structure and molecular dynamics simulations study

Essentials Consensus sequence and biochemical data suggest a Na+ -site in the factor (F) IXa protease domain. X-ray structure of the FIXa EGF2/protease domain at 1.37 Å reveals a Na+ -site not observed earlier. Molecular dynamics simulations data support that Na+  ± Ca2+ promote FIXa protease domain stability. Sulfate ions found in the protease domain mimic heparin sulfate binding mode in FIXa. SUMMARY: Background Activated coagulation factor IX (FIXa) consists of a γ-carboxyglutamic acid domain, two epidermal growth factor-like (EGF) domains, and a C-terminal protease domain. Consensus sequence and biochemical data support the existence of a Na+ -site in the FIXa protease domain. However, soaking experiments or crystals grown in high concentration of ammonium sulfate did not reveal a Na+ -site in wild-type or mutant FIXa EGF2/protease domain structure. Objective Determine the structure of the FIXa EGF2/protease domain in the presence of Na+ ; perform molecular dynamics (MD) simulations to explore the role of Na+ in stabilizing FIXa structure. Methods Crystallography, MD simulations, and modeling heparin binding to FIXa. Results Crystal structure at 1.37-Å resolution revealed that Na+ is coordinated to carbonyl groups of residues 184A, 185, 221A, and 224 in the FIXa protease domain. The Na+ -site in FIXa is similar to that of FXa and is linked to the Asp189 S1-site. In MD simulations, Na+ reduced fluctuations in residues 217-225 (Na+ -loop) and 70-80 (Ca2+ -loop), whereas Ca2+ reduced fluctuations only in residues of the Ca2+ -loop. Ca2+ and Na+ together reduced fluctuations in residues of the Ca2+ -loop and Na+ -loop (residues 70-80, 183-194, and 217-225). Moreover, we observed four sulfate ions that make salt bridges with FIXa protease domain Arg/Lys residues, which have been implicated in heparin binding. Based upon locations of the sulfate ions, we modeled heparin binding to FIXa, which is similar to the heparin binding in thrombin. Conclusions The FIXa Na+ -site in association with Ca2+ contributes to stabilization of the FIXa protease domain. The heparin binding mode in FIXa is similar to that in thrombin.

N+-C-H···O Hydrogen bonds in protein-ligand complexes

In the context of drug design, C-H···O hydrogen bonds have received little attention so far, mostly because they are considered weak relative to other noncovalent interactions such as O-H···O hydrogen bonds, π/π interactions, and van der Waals interactions. Herein, we demonstrate the significance of hydrogen bonds between C-H groups adjacent to an ammonium cation and an oxygen atom (N⁺-C-H···O hydrogen bonds) in protein-ligand complexes. Quantum chemical calculations revealed details on the strength and geometrical requirements of these N⁺-C-H···O hydrogen bonds, and a subsequent survey of the Protein Data Bank (PDB) based on these criteria suggested that numerous protein-ligand complexes contain such N⁺-C-H···O hydrogen bonds. An ensuing experimental investigation into the G9a-like protein (GLP)-inhibitor complex demonstrated that N⁺-C-H···O hydrogen bonds affect the activity of the inhibitors against the target enzyme. These results should provide the basis for the use of N⁺-C-H···O hydrogen bonds in drug discovery.

Discovery of a Novel Highly Selective Histamine H4 Receptor Antagonist for the Treatment of Atopic Dermatitis

The histamine H4 receptor (H4R), a member of the G-protein coupled receptor family, has been considered as a potential therapeutic target for treating atopic dermatitis (AD). A large number of H4R antagonists have been disclosed, but no efficient agents controlling both pruritus and inflammation in AD have been developed yet. Here, we have discovered a novel class of orally available H4R antagonists showing strong anti-itching and anti-inflammation activity as well as excellent selectivity against off-targets. A pharmacophore-based virtual screening system constructed in-house successfully identified initial hit compound 9, and the subsequent homology model-guided optimization efficiently led us to discover pyrido[2,3- e]tetrazolo[1,5- a]pyrazine analogue 48 as a novel chemotype of a potent and highly selective H4R antagonist. Importantly, orally administered compound 48 exhibits remarkable efficacy on antipruritus and anti-inflammation with a favorable pharmacokinetic (PK) profile in several mouse models of AD. Thus, these data strongly suggest that our compound 48 is a promising clinical candidate for treatment of AD.

Recognition of shorter and longer trimethyllysine analogues by epigenetic reader proteins

Combined thermodynamic data, molecular dynamics simulations, and quantum chemical studies reveal that epigenetic reader proteins efficiently bind trimethylornithine and trimethylhomolysine.

A systematic analysis of atomic protein-ligand interactions in the PDB

As the protein databank (PDB) recently passed the cap of 123 456 structures, it stands more than ever as an important resource not only to analyze structural features of specific biological systems, but also to study the prevalence of structural patterns observed in a large body of unrelated structures, that may reflect rules governing protein folding or molecular recognition. Here, we compiled a list of 11 016 unique structures of small-molecule ligands bound to proteins - 6444 of which have experimental binding affinity - representing 750 873 protein-ligand atomic interactions, and analyzed the frequency, geometry and impact of each interaction type. We find that hydrophobic interactions are generally enriched in high-efficiency ligands, but polar interactions are over-represented in fragment inhibitors. While most observations extracted from the PDB will be familiar to seasoned medicinal chemists, less expected findings, such as the high number of C-H···O hydrogen bonds or the relatively frequent amide-π stacking between the backbone amide of proteins and aromatic rings of ligands, uncover underused ligand design strategies.

Amine promiscuity and toxicology analysis

Drug discovery programs often face challenges to obtain sufficient duration of action of the drug (i.e. seek longer half-lives). If the pharmacodynamic response is driven by free plasma concentration of the drug then extending the plasma drug concentration is a valid approach. Half-life is dependent on the volume of distribution, which in turn can be dependent upon the ionization state of the molecule. Basic compounds tend to have a higher volume of distribution leading to longer half-lives. However, it has been shown that bases may also have higher promiscuity. In this work, we describe an analysis of in vitro pharmacological profiling and toxicology data investigating the role of primary, secondary, and tertiary amines in imparting promiscuity and thus off-target toxicity. Primary amines are found to be less promiscuous in in vitro assays and have improved profiles in in vivo toxicology studies compared to secondary and tertiary amines.

Experimental Binding Energies in Supramolecular Complexes

On the basis of many literature measurements, a critical overview is given on essential noncovalent interactions in synthetic supramolecular complexes, accompanied by analyses with selected proteins. The methods, which can be applied to derive binding increments for single noncovalent interactions, start with the evaluation of consistency and additivity with a sufficiently large number of different host-guest complexes by applying linear free energy relations. Other strategies involve the use of double mutant cycles, of molecular balances, of dynamic combinatorial libraries, and of crystal structures. Promises and limitations of these strategies are discussed. Most of the analyses stem from solution studies, but a few also from gas phase. The empirically derived interactions are then presented on the basis of selected complexes with respect to ion pairing, hydrogen bonding, electrostatic contributions, halogen bonding, π-π-stacking, dispersive forces, cation-π and anion-π interactions, and contributions from the hydrophobic effect. Cooperativity in host-guest complexes as well as in self-assembly, and entropy factors are briefly highlighted. Tables with typical values for single noncovalent free energies and polarity parameters are in the Supporting Information.

Regulation of protein-ligand binding affinity by hydrogen bond pairing

Hydrogen (H)-bonds potentiate diverse cellular functions by facilitating molecular interactions. The mechanism and the extent to which H-bonds regulate molecular interactions are a largely unresolved problem in biology because the H-bonding process continuously competes with bulk water. This interference may significantly alter our understanding of molecular function, for example, in the elucidation of the origin of enzymatic catalytic power. We advance this concept by showing that H-bonds regulate molecular interactions via a hitherto unappreciated donor-acceptor pairing mechanism that minimizes competition with water. On the basis of theoretical and experimental correlations between H-bond pairings and their effects on ligand binding affinity, we demonstrate that H-bonds enhance receptor-ligand interactions when both the donor and acceptor have either significantly stronger or significantly weaker H-bonding capabilities than the hydrogen and oxygen atoms in water. By contrast, mixed strong-weak H-bond pairings decrease ligand binding affinity due to interference with bulk water, offering mechanistic insight into why indiscriminate strengthening of receptor-ligand H-bonds correlates poorly with experimental binding affinity. Further support for the H-bond pairing principle is provided by the discovery and optimization of lead compounds targeting dietary melamine and Clostridium difficile toxins, which are not realized by traditional drug design methods. Synergistic H-bond pairings have therefore evolved in the natural design of high-affinity binding and provide a new conceptual framework to evaluate the H-bonding process in biological systems. Our findings may also guide wider applications of competing H-bond pairings in lead compound design and in determining the origin of enzymatic catalytic power.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Molecular Basis of Enhanced Activity in Factor VIIa-Trypsin Variants Conveys Insights into Tissue Factor-mediated Allosteric Regulation of Factor VIIa Activity

The complex of coagulation factor VIIa (FVIIa), a trypsin-like serine protease, and membrane-bound tissue factor (TF) initiates blood coagulation upon vascular injury. Binding of TF to FVIIa promotes allosteric conformational changes in the FVIIa protease domain and improves its catalytic properties. Extensive studies have revealed two putative pathways for this allosteric communication. Here we provide further details of this allosteric communication by investigating FVIIa loop swap variants containing the 170 loop of trypsin that display TF-independent enhanced activity. Using x-ray crystallography, we show that the introduced 170 loop from trypsin directly interacts with the FVIIa active site, stabilizing segment 215-217 and activation loop 3, leading to enhanced activity. Molecular dynamics simulations and novel fluorescence quenching studies support that segment 215-217 conformation is pivotal to the enhanced activity of the FVIIa variants. We speculate that the allosteric regulation of FVIIa activity by TF binding follows a similar path in conjunction with protease domain N terminus insertion, suggesting a more complete molecular basis of TF-mediated allosteric enhancement of FVIIa activity.

New Infestin-4 Mutants with Increased Selectivity against Factor XIIa

Factor XIIa (fXIIa) is a serine protease that triggers the coagulation contact pathway and plays a role in thrombosis. Because it interferes with coagulation testing, the need to inhibit fXIIa exists in many cases. Infestin-4 (Inf4) is a Kazal-type inhibitor of fXIIa. Its specificity for fXIIa can be enhanced by point mutations in the protease-binding loop. We attempted to adapt Inf4 for the selective repression of the contact pathway under various in vitro conditions, e.g., during blood collection and in ‘global’ assays of tissue factor (TF)-dependent coagulation. First, we designed a set of new Inf4 mutants that, in contrast to wt-Inf4, had stabilized canonical conformations during molecular dynamics simulation. Off-target activities against factor Xa (fXa), plasmin, and other coagulation proteases were either reduced or eliminated in these recombinant mutants, as demonstrated by chromogenic assays. Interactions with fXIIa and fXa were also analyzed using protein-protein docking. Next, Mutant B, one of the most potent mutants (its K_i for fXIIa is 0.7 nM) was tested in plasma. At concentrations 5–20 μM, this mutant delayed the contact-activated generation of thrombin, as well as clotting in thromboelastography and thrombodynamics assays. In these assays, Mutant B did not affect coagulation initiated by TF, thus demonstrating sufficient selectivity and its potential practical significance as a reagent for coagulation diagnostics.

A theoretical study of ternary indole-cation-anion complexes

The simultaneous interactions of an anion and a cation with a π system were investigated by MP2 and M06-2X theoretical calculations. Indole was chosen as a model π system for its relevance in biological environments. Two different orientations of the anion, interacting with the N-H and with the C-H groups of indole, were considered. The four cations (Na(+), NH4(+), C(NH2)3(+) and N(CH3)4(+)) and the four anions (Cl(-), NO3(-), HCOO(-) and BF4(-)) included in the study are of biological interest. The total interaction energy of the ternary complexes was calculated and separated into its two- and three-body components and all of them are further divided into their electrostatic, exchange, repulsion, polarization and dispersion contributions using the local molecular orbital-energy decomposition analysis (LMO-EDA) methodology. The binding energy of the indole-cation-anion complexes depends on both ions, with the cation having the strongest effect. The intense cation-anion attraction determines the geometric and energetic features in all ternary complexes. These structures, with both ions on the same side of the π system, show an anti-cooperative interaction. However, the interaction is not only determined by electrostatics, but also the polarization contribution is important. Specific interactions like the one established between the anion and the N-H group of indole or the proton transfer between an acidic cation and a basic anion play a significant role in the energetics and the structure of particular complexes. The presence of the polar solvent as modelled with the polarizable continuum model (PCM) does not seem to have a significant effect on the geometry of the ternary complexes, but drastically weakens the interaction energy. Also, the strength of the interaction is reduced at a faster rate when the anion is pushed away, compared to the results obtained in the gas phase. The combination of PCM with the addition of one water molecule indicates that the PCM method properly reproduces the main energetic and geometrical changes, even at the quantitative level, but the explicit hydration allows refining the solvent effect and detecting cases that do not follow the general trend.

Correlating structure and ligand affinity in drug discovery: a cautionary tale involving second shell residues

A high-resolution crystallographic structure determination of a protein-ligand complex is generally accepted as the 'gold standard' for structure-based drug design, yet the relationship between structure and affinity is neither obvious nor straightforward. Here we analyze the interactions of a series of serine proteinase inhibitors with trypsin variants onto which the ligand-binding site of factor Xa has been grafted. Despite conservative mutations of only two residues not immediately in contact with ligands (second shell residues), significant differences in the affinity profiles of the variants are observed. Structural analyses demonstrate that these are due to multiple effects, including differences in the structure of the binding site, differences in target flexibility and differences in inhibitor binding modes. The data presented here highlight the myriad competing microscopic processes that contribute to protein-ligand interactions and emphasize the difficulties in predicting affinity from structure.

Molecular recognition in chemical and biological systems

Structure-based ligand design in medicinal chemistry and crop protection relies on the identification and quantification of weak noncovalent interactions and understanding the role of water. Small-molecule and protein structural database searches are important tools to retrieve existing knowledge. Thermodynamic profiling, combined with X-ray structural and computational studies, is the key to elucidate the energetics of the replacement of water by ligands. Biological receptor sites vary greatly in shape, conformational dynamics, and polarity, and require different ligand-design strategies, as shown for various case studies. Interactions between dipoles have become a central theme of molecular recognition. Orthogonal interactions, halogen bonding, and amide⋅⋅⋅π stacking provide new tools for innovative lead optimization. The combination of synthetic models and biological complexation studies is required to gather reliable information on weak noncovalent interactions and the role of water.

A series of potent CREBBP bromodomain ligands reveals an induced-fit pocket stabilized by a cation-π interaction

The benzoxazinone and dihydroquinoxalinone fragments were employed as novel acetyl lysine mimics in the development of CREBBP bromodomain ligands. While the benzoxazinone series showed low affinity for the CREBBP bromodomain, expansion of the dihydroquinoxalinone series resulted in the first potent inhibitors of a bromodomain outside the BET family. Structural and computational studies reveal that an internal hydrogen bond stabilizes the protein-bound conformation of the dihydroquinoxalinone series. The side chain of this series binds in an induced-fit pocket forming a cation–π interaction with R1173 of CREBBP. The most potent compound inhibits binding of CREBBP to chromatin in U2OS cells.

Cation-π interactions of methylated ammonium ions: a quantum mechanical study

Cation-π interactions of methylated ammonium ions play a key role in a broad range of biochemical systems. These include methyl-lysine binding proteins which bind to methylated sites on histone proteins, lysine demethylase enzymes which demethylate these sites, and neurotransmitter receptor complexes which bind choline-derived ligands. Recognition in these systems is achieved through an 'aromatic cage' motif in the binding site. Here we use high-level quantum mechanical calculations to address how cation-π interactions of methylated ammonium ions are modulated by a change in methylation state and interaction geometry. We survey methyl-lysine and choline-derived complexes in the Protein Databank to validate our results against available structural data. A quantitative description of cation-π interactions of methylated ammonium systems is critical to structure-based efforts to target methyl-lysine binding proteins and demethylase enzymes in the treatment of unregulated transcriptional control, and neurotransmitter receptors in the treatment of neurological disease. It is our hope that our work will serve as a benchmark for the development of physical chemistry based force fields that can accurately model the contribution of cation-π interactions to binding and specificity in these systems.

Quantifying additive interactions of the osmolyte proline with individual functional groups of proteins: comparisons with urea and glycine betaine, interpretation of m-values

To quantify interactions of the osmolyte l-proline with protein functional groups and predict their effects on protein processes, we use vapor pressure osmometry to determine chemical potential derivatives dμ2/dm3 = μ23, quantifying the preferential interactions of proline (component 3) with 21 solutes (component 2) selected to display different combinations of aliphatic or aromatic C, amide, carboxylate, phosphate or hydroxyl O, and amide or cationic N surface. Solubility data yield μ23 values for four less-soluble solutes. Values of μ23 are dissected using an ASA-based analysis to test the hypothesis of additivity and obtain α-values (proline interaction potentials) for these eight surface types and three inorganic ions. Values of μ23 predicted from these α-values agree with the experiment, demonstrating additivity. Molecular interpretation of α-values using the solute partitioning model yields partition coefficients (Kp) quantifying the local accumulation or exclusion of proline in the hydration water of each functional group. Interactions of proline with native protein surfaces and effects of proline on protein unfolding are predicted from α-values and ASA information and compared with experimental data, with results for glycine betaine and urea, and with predictions from transfer free energy analysis. We conclude that proline stabilizes proteins because of its unfavorable interactions with (exclusion from) amide oxygens and aliphatic hydrocarbon surfaces exposed in unfolding and that proline is an effective in vivo osmolyte because of the osmolality increase resulting from its unfavorable interactions with anionic (carboxylate and phosphate) and amide oxygens and aliphatic hydrocarbon groups on the surface of cytoplasmic proteins and nucleic acids.

Crystal structures of a plant trypsin inhibitor from Enterolobium contortisiliquum (EcTI) and of its complex with bovine trypsin

A serine protease inhibitor from Enterolobium contortisiliquum (EcTI) belongs to the Kunitz family of plant inhibitors, common in plant seeds. It was shown that EcTI inhibits the invasion of gastric cancer cells through alterations in integrin-dependent cell signaling pathway. We determined high-resolution crystal structures of free EcTI (at 1.75 Å) and complexed with bovine trypsin (at 2 Å). High quality of the resulting electron density maps and the redundancy of structural information indicated that the sequence of the crystallized isoform contained 176 residues and differed from the one published previously. The structure of the complex confirmed the standard inhibitory mechanism in which the reactive loop of the inhibitor is docked into trypsin active site with the side chains of Arg64 and Ile65 occupying the S1 and S1′ pockets, respectively. The overall conformation of the reactive loop undergoes only minor adjustments upon binding to trypsin. Larger deviations are seen in the vicinity of Arg64, driven by the needs to satisfy specificity requirements. A comparison of the EcTI-trypsin complex with the complexes of related Kunitz inhibitors has shown that rigid body rotation of the inhibitors by as much as 15° is required for accurate juxtaposition of the reactive loop with the active site while preserving its conformation. Modeling of the putative complexes of EcTI with several serine proteases and a comparison with equivalent models for other Kunitz inhibitors elucidated the structural basis for the fine differences in their specificity, providing tools that might allow modification of their potency towards the individual enzymes.

Maleimide-functionalized photochromic spirodihydroindolizines

Two photochromic spirodihydroindolizine/betaine systems for tethering to peptides and proteins via a maleimide function have been prepared. The absorption spectra of the betaines are in the red region of the visible spectrum and in the near-IR spectral domain, which are suitable energies of light for future in vivo applications. The half-times of cyclization have been determined for both DHI/betaine systems. The findings are consistent with a thermal barrier of varying size between the transoid and cisoid conformers of the betaines.