Understanding binding selectivity toward trypsin and factor Xa: the role of aromatic interactions

Cloning and characterization of Aedes aegypti blood downregulated chymotrypsin II

Aedes aegypti adult and larval blood downregulated chymotrypsin II was cloned, sequenced and its 3D conformation modeled. Cloning of the enzymes from adult and larval guts indicated that both genes sit at the same location on Chromosome 2. Genomic analyses showed that larval and adult genes are the same and both have four exons and three introns that are located on an 8.32 Kb DNA in direction with the Ae. aegypti genome. The adult and larval transcript synthesis is controlled by alternative splicing explaining small difference in the amino acids sequences. Chymotrypsin II that was extracted from guts of sugar-fed and at 48 after blood feeding showed a pH optimum of 4-5 with a broad shoulder of activity from pH 6 to 10. Dot blot analyses show that the enzyme's transcript is downregulated after females take a blood meal and upregulated at 48 h after the blood meal. A Chymotrypsin II transcript was also detected in the larval gut during different times of larval developmental stages, indication that Ae. aegypti chymotrypsin II is synthesized by adults and larval guts. The possibility that JH III and 20HE play an active role in the regulation is discussed.

Fast and Accurate Quantum Mechanical Modeling of Large Molecular Systems Using Small Basis Set Hartree-Fock Methods Corrected with Atom-Centered Potentials

There has been significant interest in developing fast and accurate quantum mechanical methods for modeling large molecular systems. In this work, by utilizing a machine learning regression technique, we have developed new low-cost quantum mechanical approaches to model large molecular systems. The developed approaches rely on using one-electron Gaussian-type functions called atom-centered potentials (ACPs) to correct for the basis set incompleteness and the lack of correlation effects in the underlying minimal or small basis set Hartree-Fock (HF) methods. In particular, ACPs are proposed for ten elements common in organic and bioorganic chemistry (H, B, C, N, O, F, Si, P, S, and Cl) and four different base methods: two minimal basis sets (MINIs and MINIX) plus a double-ζ basis set (6-31G*) in combination with dispersion-corrected HF (HF-D3/MINIs, HF-D3/MINIX, HF-D3/6-31G*) and the HF-3c method. The new ACPs are trained on a very large set (73 832 data points) of noncovalent properties (interaction and conformational energies) and validated additionally on a set of 32 048 data points. All reference data are of complete basis set coupled-cluster quality, mostly CCSD(T)/CBS. The proposed ACP-corrected methods are shown to give errors in the tenths of a kcal/mol range for noncovalent interaction energies and up to 2 kcal/mol for molecular conformational energies. More importantly, the average errors are similar in the training and validation sets, confirming the robustness and applicability of these methods outside the boundaries of the training set. In addition, the performance of the new ACP-corrected methods is similar to complete basis set density functional theory (DFT) but at a cost that is orders of magnitude lower, and the proposed ACPs can be used in any computational chemistry program that supports effective-core potentials without modification. It is also shown that ACPs improve the description of covalent and noncovalent bond geometries of the underlying methods and that the improvement brought about by the application of the ACPs is directly related to the number of atoms to which they are applied, allowing the treatment of systems containing some atoms for which ACPs are not available. Overall, the ACP-corrected methods proposed in this work constitute an alternative accurate, economical, and reliable quantum mechanical approach to describe the geometries, interaction energies, and conformational energies of systems with hundreds to thousands of atoms.

Potent inhibition of enterovirus D68 and human rhinoviruses by dipeptidyl aldehydes and α-ketoamides

Enterovirus D68 (EV-D68) is an emerging pathogen responsible for mild to severe respiratory infections that occur mostly in infants, children and teenagers. EV-D68, one of more than 100 non-polio enteroviruses, is acid-labile and biologically similar to human rhinoviruses (HRV) (originally classified as HRV87). However, there is no approved preventive or therapeutic measure against EV-D68, HRV, or other enteroviruses. In this study, we evaluated the antiviral activity of series of dipeptidyl compounds against EV-D68 and HRV strains, and demonstrated that several peptidyl aldehyde and α-ketoamide peptidyl compounds are potent inhibitors of EV-D68 and HRV strains with high in-vitro therapeutic indices (>1000). One of the α-ketoamide compounds is shown to have favorable pharmacokinetics profiles, including a favorable oral bioavailability in rats. Recent successful development of α-ketoamide protease inhibitors against hepatitis C virus suggests these compounds may have a high potential for further optimization and development against emerging EV-D68, as well as HRV.

•

Series of dipeptidyl aldehyde or ketomide compounds were highly effective against enterovirus-D68 and human rhinoviruses.

•

The highly effective ketoamide compound is shown to have favorable pharmacokinetics profiles.

•

These may have a high potential for further antiviral development against emerging enterovirus-D68 and human rhinoviruses.

Development of Orally Active Thrombin Inhibitors for the Treatment of Thrombotic Disorder Diseases

Thrombotic disorders represent the major share of the various cardiovascular diseases, and significant progress has been made in the development of synthetic thrombin inhibitors as new anticoagulants. In addition to the development of highly potent and selective inhibitors with improved safety and suitable half-life, several allosteric inhibitors have been designed and synthesized, that did not fully nullify the procoagulant signal and thus could result in reduced bleeding complications. Furthermore, natural products with thrombin inhibitory activity have been isolated, and some natural products have been modified in order to improve their inhibitory activity and metabolic stability. This review summarizes the development of orally active thrombin inhibitors for the treatment of thrombotic disorder diseases, which could serve as a reference for the interested researchers.

Structure-guided design and optimization of dipeptidyl inhibitors of norovirus 3CL protease. Structure-activity relationships and biochemical, X-ray crystallographic, cell-based, and in vivo studies

Norovirus infection constitutes the primary cause of acute viral gastroenteritis. There are currently no vaccines or norovirus-specific antiviral therapeutics available for the management of norovirus infection. Norovirus 3C-like protease is essential for viral replication, consequently, inhibition of this enzyme is a fruitful avenue of investigation that may lead to the emergence of antinorovirus therapeutics. We describe herein the optimization of dipeptidyl inhibitors of norovirus 3C-like protease using iterative SAR, X-ray crystallographic, and enzyme and cell-based studies. We also demonstrate herein in vivo efficacy of an inhibitor using the murine model of norovirus infection.

Hitting the target: fragment screening with acoustic in situ co-crystallization of proteins plus fragment libraries on pin-mounted data-collection micromeshes

Acoustic droplet ejection (ADE) is a powerful technology that supports crystallographic applications such as growing, improving and manipulating protein crystals. A fragment-screening strategy is described that uses ADE to co-crystallize proteins with fragment libraries directly on MiTeGen MicroMeshes. Co-crystallization trials can be prepared rapidly and economically. The high speed of specimen preparation and the low consumption of fragment and protein allow the use of individual rather than pooled fragments. The Echo 550 liquid-handling instrument (Labcyte Inc., Sunnyvale, California, USA) generates droplets with accurate trajectories, which allows multiple co-crystallization experiments to be discretely positioned on a single data-collection micromesh. This accuracy also allows all components to be transferred through small apertures. Consequently, the crystallization tray is in equilibrium with the reservoir before, during and after the transfer of protein, precipitant and fragment to the micromesh on which crystallization will occur. This strict control of the specimen environment means that the crystallography experiments remain identical as the working volumes are decreased from the few microlitres level to the few nanolitres level. Using this system, lysozyme, thermolysin, trypsin and stachydrine demethylase crystals were co-crystallized with a small 33-compound mini-library to search for fragment hits. This technology pushes towards a much faster, more automated and more flexible strategy for structure-based drug discovery using as little as 2.5 nl of each major component.

Potent inhibition of norovirus by dipeptidyl α-hydroxyphosphonate transition state mimics

The design, synthesis, and evaluation of a series of dipeptidyl α-hydroxyphosphonates is reported. The synthesized compounds displayed high anti-norovirus activity in a cell-based replicon system, as well as high enzyme selectivity.

Quantifying protein-ligand binding constants using electrospray ionization mass spectrometry: a systematic binding affinity study of a series of hydrophobically modified trypsin inhibitors

NanoESI-MS is used for determining binding strengths of trypsin in complex with two different series of five congeneric inhibitors, whose binding affinity in solution depends on the size of the P3 substituent. The ligands of the first series contain a 4-amidinobenzylamide as P1 residue, and form a tight complex with trypsin. The inhibitors of the second series have a 2-aminomethyl-5-chloro-benzylamide as P1 group, and represent a model system for weak binders. The five different inhibitors of each group are based on the same scaffold and differ only in the length of the hydrophobic side chain of their P3 residue, which modulates the interactions in the S3/4 binding pocket of trypsin. The dissociation constants (K(D)) for high affinity ligands investigated by nanoESI-MS ranges from 15 nM to 450 nM and decreases with larger hydrophobic P3 side chains. Collision-induced dissociation (CID) experiments of five trypsin and benzamidine-based complexes show a correlation between trends in K(D) and gas-phase stability. For the second inhibitor series we could show that the effect of imidazole, a small stabilizing additive, can avoid the dissociation of the complex ions and as a result increases the relative abundance of weakly bound complexes. Here the K(D) values ranging from 2.9 to 17.6 μM, some 1-2 orders of magnitude lower than the first series. For both ligand series, the dissociation constants (K(D)) measured via nanoESI-MS were compared with kinetic inhibition constants (K(i)) in solution.

Prediction of trypsin/molecular fragment binding affinities by free energy decomposition and empirical scores

Two families of binding affinity estimation methodologies are described which were utilized in the SAMPL3 trypsin/fragment binding affinity challenge. The first is a free energy decomposition scheme based on a thermodynamic cycle, which included separate contributions from enthalpy and entropy of binding as well as a solvent contribution. Enthalpic contributions were estimated with PM6-DH2 semiempirical quantum mechanical interaction energies, which were modified with a statistical error correction procedure. Entropic contributions were estimated with the rigid-rotor harmonic approximation, and solvent contributions to the free energy were estimated with several different methods. The second general methodology is the empirical score LISA, which contains several physics-based terms trained with the large PDBBind database of protein/ligand complexes. Here we also introduce LISA+, an updated version of LISA which, prior to scoring, classifies systems into one of four classes based on a ligand's hydrophobicity and molecular weight. Each version of the two methodologies (a total of 11 methods) was trained against a compiled set of known trypsin binders available in the Protein Data Bank to yield scaling parameters for linear regression models. Both raw and scaled scores were submitted to SAMPL3. Variants of LISA showed relatively low absolute errors but also low correlation with experiment, while the free energy decomposition methods had modest success when scaling factors were included. Nonetheless, re-scaled LISA yielded the best predictions in the challenge in terms of RMS error, and six of these models placed in the top ten best predictions by RMS error. This work highlights some of the difficulties of predicting binding affinities of small molecular fragments to protein receptors as well as the benefit of using training data.

Aromatic rings in chemical and biological recognition: energetics and structures

This review describes a multidimensional treatment of molecular recognition phenomena involving aromatic rings in chemical and biological systems. It summarizes new results reported since the appearance of an earlier review in 2003 in host-guest chemistry, biological affinity assays and biostructural analysis, data base mining in the Cambridge Structural Database (CSD) and the Protein Data Bank (PDB), and advanced computational studies. Topics addressed are arene-arene, perfluoroarene-arene, S⋅⋅⋅aromatic, cation-π, and anion-π interactions, as well as hydrogen bonding to π systems. The generated knowledge benefits, in particular, structure-based hit-to-lead development and lead optimization both in the pharmaceutical and in the crop protection industry. It equally facilitates the development of new advanced materials and supramolecular systems, and should inspire further utilization of interactions with aromatic rings to control the stereochemical outcome of synthetic transformations.

Structural binding evidence of the trypanocidal drugs berenil and pentacarinate active principles to a serine protease model

Bovine trypsin is a model system for the serine protease class of enzymes, which is an important target for contemporary medicinal chemistry. Some structural and thermodynamic reports are available on its interaction with benzamidine-based compounds but no structural information is available so far on its binding modes to the active principles of the trypanocidal drugs Pentacarinate (pentamidine) and Berenil (diminazene). The crystallographic structures of bovine beta-trypsin in complex with the ligands were determined to a resolution of 1.57 A (diminazene) and 1.70 A (diminazene and pentamidine). The second benzamidine moieties in these inhibitors are bound to the enzyme in different hot spots and only few hydrogen bonds mediate these interactions. Thermodynamic parameters for the association of pentamidine with beta-trypsin reveal that this inhibitor has about 1.3-fold lower affinity than diminazene. Moreover its binding mode resembles other benzamidine-based compounds that assess the aryl binding pocket of the enzyme; however, with almost 2.5-fold higher affinity. This is the first structural evidence of the binding of Berenil and Pentacarinate active principles trypanocidal drugs to serine proteases.

Engineering protein allostery: 1.05 A resolution structure and enzymatic properties of a Na+-activated trypsin

Some trypsin-like proteases are endowed with Na(+)-dependent allosteric enhancement of catalytic activity, but this important mechanism has been difficult to engineer in other members of the family. Replacement of 19 amino acids in Streptomyces griseus trypsin targeting the active site and the Na(+)-binding site were found necessary to generate efficient Na(+) activation. Remarkably, this property was linked to the acquisition of a new substrate selectivity profile similar to that of factor Xa, a Na(+)-activated protease involved in blood coagulation. The X-ray crystal structure of the mutant trypsin solved to 1.05 A resolution defines the engineered Na(+) site and active site loops in unprecedented detail. The results demonstrate that trypsin can be engineered into an efficient allosteric protease, and that Na(+) activation is interwoven with substrate selectivity in the trypsin scaffold.