Factor Xa: Simulation studies with an eye to inhibitor design

Factor Xa is a serine protease which activates thrombin and plays a key regulatory role in the blood-coagulation cascade. Factor Xa is at the crossroads of the extrinsic and intrinsic pathways of coagulation and, hence, has become an important target for the design of anti-thrombotics (inhibitors). It is not known to be involved in other processes than hemostasis and its binding site is different to that of other serine proteases, thus facilitating selective inhibition. The design of high-affinity selective inhibitors of factor Xa requires knowledge of the structural and dynamical characteristics of its active site. The three-dimensional structure of factor Xa was resolved by X-ray crystallography and refined at 2.2 Å resolution by Padmanabhan and collaborators. In this article we present results from molecular dynamics simulations of the catalytic domain of factor Xa in aqueous solution. The simulations were performed to characterise the mobility and flexibility of the residues delimiting the unoccupied binding site of the enzyme, and to determine hydrogen bonding propensities (with protein and with solvent atoms) of those residues in the active site that could interact with a substrate or a potential inhibitor. The simulation data is aimed at facilitating the design of high-affinity selective inhibitors of factor Xa.     

Efficiency of connected binary block designs when a single observation is unavailable

In this paper the problem of finding the design efficiency is considered when a single observation is unavailable in a connected binary block design. The explicit expression of efficiency is found for the resulting design when the original design is a balanced incomplete block design or a group divisible, singular or semiregular or regular with ₁>0, design. The efficiency does not depend on the position of the unavailable observation. For a regular group divisible design with ₁>0, the efficiency depends on the position of the unavailable observation. The bounds, both lower and upper, on the efficiency are given in this situation. The efficiencies of designs resulting from a balanced incomplete block design and a group divisible design are in fact high when a single observation is unavailable.The work of the first author is sponsored by the Air Force Office of Scientific Research under Grant AFOSR-90-0092.On leave from Indian Statistical Institute, Calcutta, India. The work of the third author was supported by a grant from the CMDS, Indian Institute of Management, Calcutta.     

Modified AutoDock for accurate docking of protein kinase inhibitors

Protein kinases are an important class of enzymes controlling virtually all cellular signaling pathways. Consequently, selective inhibitors of protein kinases have attracted significant interest as potential new drugs for many diseases. Computational methods, including molecular docking, have increasingly been used in the inhibitor design process [1]. We have considered several docking packages in order to strengthen our kinase inhibitor work with computational capabilities. In our experience, AutoDock offered a reasonable combination of accuracy and speed, as opposed to methods that specialize either in fast database searches or detailed and computationally intensive calculations.However, AutoDock did not perform well in cases where extensive hydrophobic contacts were involved, such as docking of SB203580 to its target protein kinase p38. Another shortcoming was a hydrogen bonding energy function, which underestimated the attraction component and, thus, did not allow for sufficiently accurate modeling of the key hydrogen bonds in the kinase-inhibitor complexes.We have modified the parameter set used to model hydrogen bonds, which increased the accuracy of AutoDock and appeared to be generally applicable to many kinase-inhibitor pairs without customization. Binding to largely hydrophobic sites, such as the active site of p38, was significantly improved by introducing a correction factor selectively affecting only carbon and hydrogen energy grids, thus, providing an effective, although approximate, treatment of solvation.     

Modelling of GC and HPLC separations of simazin and atrazin by experimental design methodology

Summary Experimental Design methodology allows the modelling and optimization of the chromatographic separation of similar pesticides (triazine family) by GC and HPLC. The GC separation of simazin and atrazin is well modelled by a first degree equation, involving injected volume, carrier gas pressure and rising oven temperature. The LC is modelled by a second degree equation, depending on injected volume, eluent flow and composition. These calculated models allow easy optimization of the separations, using isoresponse curves.     

Computational combinatorial ligand design: Application to human α-thrombin

Summary A new method is presented for computer-aided ligand design by combinatorial selection of fragments that bind favorably to a macromolecular target of known three-dimensional structure. Firstly, the multiple-copy simultaneous-search procedure (MCSS) is used to exhaustively search for optimal positions and orientations of functional groups on the surface of the macromolecule (enzyme or receptor fragment). The MCSS minima are then sorted according to an approximated binding free energy, whose solvation component is expressed as a sum of separate electrostatic and nonpolar contributions. The electrostatic solvation energy is calculated by the numerical solution of the linearized Poisson-Boltzmann equation, while the nonpolar contribution to the binding free energy is assumed to be proportional to the loss in solvent-accessible surface area. The program developed for computational combinatorial ligand design (CCLD) allows the fast and automatic generation of a multitude of highly diverse compounds, by connecting in a combinatorial fashion the functional groups in their minimized positions. The fragments are linked as two atoms may be either fused, or connected by a covalent bond or a small linker unit. To avoid the combinatorial explosion problem, pruning of the growing ligand is performed according to the average value of the approximated binding free energy of its fragments. The method is illustrated here by constructing candidate ligands for the active site of human -thrombin. The MCSS minima with favorable binding free energy reproduce the interaction patterns of known inhibitors. Starting from these fragments, CCLD generates a set of compounds that are closely related to high-affinity thrombin inhibitors. In addition, putative ligands with novel binding motifs are suggested. Probable implications of the MCSS-CCLD approach for the evolving scenario of drug discovery are discussed.     

Functionality map analysis of the active site cleft of human thrombin

Summary The Multiple Copy Simultaneous Search methodology has been used to construct functionality maps for an extended region of human thrombin, including the active site. This method allows the determination of energetically favorable positions and orientations for functional groups defined by the user on the three-dimensional surface of a protein. The positions of 10 functional group sites are compared with those of corresponding groups of four thrombin-inhibitor complexes. Many, but not all features, of known thrombin inhibitors are reproduced by the method. The results indicate that certain aspects of the binding modes of these inhibitors are not optimal. In addition, suggestions are made for improving binding by interaction with functional group sites on the thrombin surface that are not used by the thrombin inhibitors. Abbreviations: MCSS, multiple copy simultaneous search; PPACK, d-phenylalanyl-l-propyl-l-arginine chloromethane; NAPAP, N -(2-naphthylsulfonylglycyl)-d-para-amidinophenylalanylpiperidine; argatroban, (2R,4R)-4-methyl-1-[N -(3-methyl-1,2,3,4-tetrahydro-8-quinolinylsulfonyl)-l-arginyl]-2-piperidine carboxylic acid; rms, root mean square. The thrombin residues are numbered according to the chymotrypsin-based numbering by Bode et al. [8]. P1, P2, P3, etc., denote the peptide inhibitor residues on the amino-terminal side of the scissile peptide bond, and S1, S2, S3, etc., the corresponding subsites of thrombin     

Thermodynamic Parameters for the Association of Fluorinated Benzenesulfonamides with Bovine Carbonic Anhydrase II

This paper describes a calorimetric study of the association of a series of seven fluorinated benzenesulfonamide ligands (C₆H_nF_5?nSO₂NH₂) with bovine carbonic anhydrase II (BCA). Quantitative structure–activity relationships between the free energy, enthalpy, and entropy of binding and pK_a and log P of the ligands allowed the evaluation of the thermodynamic parameters in terms of the two independent effects of fluorination on the ligand: its electrostatic potential and its hydrophobicity. The parameters were partitioned to the three different structural interactions between the ligand and BCA: the Zn^II cofactor–sulfonamide bond (≈65 % of the free energy of binding), the hydrogen bonds between the ligand and BCA (≈10 %), and the contacts between the phenyl ring of the ligand and BCA (≈25 %). Calorimetry revealed that all of the ligands studied bind in a 1:1 stoichiometry with BCA; this result was confirmed by ¹⁹F NMR spectroscopy and X‐ray crystallography (for complexes with human carbonic anhydrase II).     

Effects of DNA secondary structure on oligonucleotide probe binding efficiency

Secondary structure motifs in nucleic acid probes generally impair intended hybridization reactions and so efforts to predict and avoid such structures are commonly employed in probe design schemes. Another key facet of probe design that has received much less attention, however, is that secondary structure at targeted probe binding site regions may also impair hybridization. Thus, evaluation of both probe and target site secondary structures together should improve hybridization prediction and design effectiveness. Several challenges confound this goal, including imperfect empirical rules and parameters underlying predictions and the fact that folding algorithms scale poorly with respect to sequence length. Here, we attempt to quantify the consequences of target site structure on predicted hybridization using sequences sampled from the human genome. We also provide a methodology for choosing a reasonable “window size” around target sites that is as small as possible without compromising folding algorithm prediction accuracy.     

Stereoselective allyl amine synthesis through enantioselective addition of diethylzinc and [1,3]-chirality transfer: synthesis of lentiginosine and polyoxamic acid derivative

A new synthetic method for the preparation of allyl amines has been developed. The key steps of this method are enantioselective addition of diethylzinc and [1,3]-chirality transfer through the [3.3] sigmatropic rearrangement of allyl cyanates. Stereocontrolled syntheses of lentiginosine (1) and polyoxamic acid derivative 2 from a common intermediate 7 derived from D-tartaric acid (8), have been accomplished.     

基于Labview的光学空间滤波远程虚拟实验

提出了一种利用Labview实现光学空间滤波模拟的新方法,并对经典的阿贝-波特实验设计了仿真程序.该程序可根据用户设置的参量生成二维网格光栅作为整个物理过程的输入物,允许用户选择多种不同的滤波方式,如低通、高通、方向滤波等,并可依次操纵成像的各个步骤.程序可脱离编程环境独立运行,具有跨平台、界面友好、易于操作等优点,并可架设基于浏览器的远程虚拟实验.