Identification of phosphorylase kinase as a novel therapeutic target through high-throughput screening for anti-angiogenesis compounds in zebrafish.

Angiogenesis is essential for development and tumor progression. With the aim of identifying new compound inhibitors of the angiogenesis process, we used an established enhanced green fluorescent protein-transgenic zebrafish line to develop an automated assay that enables high-throughput screening of compound libraries in a whole-organism setting. Using this system, we have identified novel kinase inhibitor compounds that show anti-angiogenic properties in both zebrafish in-vivo system and in human endothelial cell in-vitro angiogenesis models. Furthermore, we have determined the kinase target of these compounds and have identified and validated a previously uncharacterized involvement of phosphorylase kinase subunit G1 (PhKG1) in angiogenesis in vivo. In addition, we have found that PhKG1 is upregulated in human tumor samples and that aberrations in gene copy number of PhK subunits are a common feature of human tumors. Our results provide a novel insight into the angiogenesis process, as well as identify new potential targets for anti-angiogenic therapies.

Design, synthesis, and biological evaluation of potent and selective amidino bicyclic factor Xa inhibitors.

Thrombotic diseases are a major cause of death and morbidity. Factor Xa (fXa) plays a vital role in the regulation of normal homeostasis and abnormal intravascular thrombus development in the blood coagulation cascade. A novel series of fXa inhibitors incorporating an amidino 6,5-fused bicyclic moiety at the P1 position has been designed and synthesized based on molecular modeling studies. Structure-activity relationship (SAR) studies have led to selective subnanomolar fXa inhibitors. The most potent fXa inhibitor in this series (72, SE170) has a potent inhibition constant (K(i) = 0.3 nM), is 350-fold selective for fXa over trypsin, and also shows good in vivo efficacy in a rabbit arterio-venous thrombosis model (ID(50) = 0.14 micromol/kg/h). An X-ray crystal structure of 72 complexed to bovine trypsin was completed, and a binding mode of 72 with fXa has been proposed based on modeling with human des-Gla-fXa.

Synthesis and activity studies of conformationally restricted alpha-ketoamide factor Xa inhibitors.

Conformationally restricted borolysine compounds containing a 2-(2-cyanophenylthio) benzoyl in the P3 position unexpectedly led to enhanced factor Xa inhibition. In an effort to improve both the potency and selectivity of this series by extending into the S' domain, we have replaced the boronic acid with alpha-ketoamides, utilizing a novel process that was developed in our labs.

Fast prediction and visualization of protein binding pockets with PASS.

PASS (Putative Active Sites with Spheres) is a simple computational tool that uses geometry to characterize regions of buried volume in proteins and to identify positions likely to represent binding sites based upon the size, shape, and burial extent of these volumes. Its utility as a predictive tool for binding site identification is tested by predicting known binding sites of proteins in the PDB using both complexed macromolecules and their corresponding apoprotein structures. The results indicate that PASS can serve as a front-end to fast docking. The main utility of PASS lies in the fact that it can analyze a moderate-size protein (approximately 30 kDa) in under 20 s, which makes it suitable for interactive molecular modeling, protein database analysis, and aggressive virtual screening efforts. As a modeling tool, PASS (i) rapidly identifies favorable regions of the protein surface, (ii) simplifies visualization of residues modulating binding in these regions, and (iii) provides a means of directly visualizing buried volume, which is often inferred indirectly from curvature in a surface representation. PASS produces output in the form of standard PDB files, which are suitable for any modeling package, and provides script files to simplify visualization in Cerius2, InsightII, MOE, Quanta, RasMol, and Sybyl. PASS is freely available to all.

The de novo design and synthesis of cyclic urea inhibitors of factor Xa: optimization of the S4 ligand.

In this report refinements to the S4 ligand group leads to compound 19, an inhibitor of fXa with good potency in vitro and an improved pharmacokinetic profile in rabbit. The X-ray crystallographic study of a representative analogue confirms our binding model for this series.

Preparation of pyrrolidine and isoxazolidine benzamidines as potent inhibitors of coagulation factor Xa.

The serine protease factor Xa is a critical enzyme in the blood coagulation cascade. Recently, the inhibition of factor Xa has begun to emerge as an attractive strategy for the discovery of novel antithrombotic agents. Here we describe pyrrolidine and isoxazolidine benzamidines as novel and potent inhibitors of factor Xa.

Rational design and synthesis of novel, potent bis-phenylamidine carboxylate factor Xa inhibitors.

The molecular modeling studies, rational design, and synthesis of a novel series of bisphenylamidine carboxylate compounds which are inhibitors of factor Xa in the blood coagulation cascade are described. Inhibition of blood coagulation has been proposed to have several potential therapeutic utilities (Kaiser and Hauptmann, Cardiovasc. Drug Rev. 1994, 12, 225-236). Factor Xa (fXa) holds a central position in the coagulation cascade (Coleman et al. in Hemostasis and Thrombosis: Basic Principles and Clinical Practice, 1994, pp 3-18). Its major role is the generation of thrombin by the proteolytic cleavage of prothrombin. Inhibition of fXa would serve to reduce the formation of platelet clots. The fXa dimer crystal structure (Tulinsky et al., J. Mol. Biol. 1993, 232, 947-966) was used in our molecular modeling studies to design a novel series of fXa inhibitors. We initially docked and minimized isolated small molecule fragments in the S1 and S4 aryl-binding subsites. Subsequently, these fragments were connected with a tether, so as not to disturb the orientation of the fragments in their respective pockets. These modeling studies led to the initial compound (1) which was found to have significant inhibitory potency for fXa (Ki = 34 nM). The synthesis of the core structure, structure-activity relationships (SAR), and proposed binding orientation based on molecular modeling for this novel bis-phenylamidine series of fXa inhibitors are described.

The de novo design and synthesis of cyclic urea inhibitors of factor Xa: initial SAR studies.

In this report we discuss the design, synthesis, and validation of a novel series of cyclic urea inhibitors of the blood coagulation protein Factor Xa. This work culminates in compound 11, a monoamidine inhibitor of fXa employing a new S4 ligand that reduces the cationic character of these analogs. Compound 11 represents a lead for a series of more potent and selective inhibitors.

A comparison of structural and dynamic properties of different simulation methods applied to SH3.

The dynamic and static properties of molecular dynamics simulations using various methods for treating solvent were compared. The SH3 protein domain was chosen as a test case because of its small size and high surface-to-volume ratio. The simulations were analyzed in structural terms by examining crystal packing, distribution of polar residues, and conservation of secondary structure. In addition, the "essential dynamics" method was applied to compare each of the molecular dynamics trajectories with a full solvent simulation. This method proved to be a powerful tool for the comparison of large concerted atomic motions in SH3. It identified methods of simulation that yielded significantly different dynamic properties compared to the full solvent simulation. Simulating SH3 using the stochastic dynamics algorithm with a vacuum (reduced charge) force field produced properties close to those of the full solvent simulation. The application of a recently described solvation term did not improve the dynamic properties. The large concerted atomic motions in the full solvent simulation as revealed by the essential dynamics method were analyzed for possible biological implications. Two loops, which have been shown to be involved in ligand binding, were seen to move in concert to open and close the ligand-binding site.

A novel search method for protein sequence--structure relations using property profiles.

In protein engineering and design it is very important that residues can be inspected in their specific environment. A standard relational database system cannot serve this purpose adequately because it cannot handle relations between individual residues. With SCAN3D we introduce a new database system for integrated sequence and structure analysis of proteins. It uses the relational paradigm wherever possible. Its main power, however, stems from the ability to retrieve stretches of consecutive residues with certain properties by comparing a property profile with all stretches of residues in the database, exploiting the ordered character of proteins. In doing so, it bypasses the large number of join operations that would be required by relational database systems. An additional advantage of using property profile matching is that searches can be carried out allowing a pre-set number of mismatches. Also, as the database is read-only, SCAN3D does not need interactive data update mechanisms. Queries typical of a molecular engineering environment are demonstrated with specific examples: analysis of peptides that induce local structure, analysis of site-dependent rotamers and residue--residue contact analysis.

How does the switch II region of G-domains work?

The transition of guanine nucleotide binding proteins between the 'on' (GTP-bound) and 'off' (GDP-bound) states has become a paradigm of molecular switching after a chemical reaction. The mechanism by which the switch signal is transmitted to the downstream recipients in the intracellular signal pathway has been extensively studied by biochemical, biophysical and genetic methods, but a clear picture of this process has yet to emerge. Based on the similarities of ras-p21 and elongation factor Tu we propose here a model of the GDP state of ras-p21 that is in agreement with all relevant experimental evidence. The model provides important clues about: (1) a possible molecular mechanism for signal transmission from the site of GTP hydrolysis to downstream effectors; (2) a major conformational change during signal generation and a key residue involved in this process (Tyr-64); and (3) regions in ras-p21 that can be differentially recognized by binding to external partners in a GTP/GDP state dependent fashion, most notably residues D69, Q70, R73, T74, R102, K104, D105 at the end of the alpha-helices 2 and 3.

New triple-helical model for the shaft of the adenovirus fibre.

The adenovirus fibre is a trimeric protein with a globular head on a long thin shaft that projects from the twelve fivefold vertices of the virion. The shaft region of the fibre primary sequence has a unique pseudo-repeating motif of 15 residues. Using constraints derived from sequence analysis, the trimeric nature of the fibre, the experimental determination of the shaft length and general knowledge about protein structure, an atomic model of the fibre shaft has been constructed by computer modelling techniques. In the final model the three monomers form a left-handed triple-helical structure with threefold symmetry and with successive 15-residue repeats on the same chain related by an axial rise of 13.1 A and a left-handed azimuthal rotation of close to 300 degrees. Three threefold related beta-sheets with short strands are formed by inter-monomer main-chain hydrogen bonds giving rise to superhelical ribbons covering the surface of the shaft. The model satisfies criteria of extensive hydrogen bonding, reasonable backbone torsion angles, burial of most hydrophobic residues and good packing of the hydrophobic core. Furthermore, the model is consistent with the observed shaft length of about 290 A and its calculated X-ray fibre diffraction patterns shows the characteristic features found in the diffraction pattern of crystals of fibre, notably layer lines with a spacing of about 1/26 A-1 and strong meridional intensity at 1/4.4 A-1.

Interchain cysteine bridges control entry of progesterone to the central cavity of the uteroglobin dimer.

The progesterone-binding protein uteroglobin has been expressed in Escherichia coli in an unfused, soluble form. Like mature uteroglobin from rabbit endometrium (UG), the E.coli produced uteroglobin (UG1) dimerizes in vitro, forms an antiparallel dimer with Cys3-Cys69' and Cys69-Cys3' disulfide bonds and binds progesterone under reducing conditions. In order to analyze the dimerization and the reduction dependence of progesterone binding in more detail, we separately replaced cysteine 3 and cysteine 69 by serines. Under reducing conditions, both uteroglobin variants (UG1-3Ser and UG1-69Ser) bind progesterone with the same affinity as the wild-type suggesting that both cysteine residues are not directly involved in progesterone binding. In contrast to the wild-type protein, both cysteine variants also bind progesterone with high affinity in the absence of reducing agents. In addition, UG1-3Ser and UG1-69Ser both form covalently linked homodimers. Thus, unnatural Cys69-69' and Cys3-3' disulfide bonds exist in UG1-3Ser and UG1-69Ser, respectively. These data together with computer models based on X-ray diffraction data strongly support the idea that progesterone reaches its binding site located in an internal hydrophobic cavity via a hydrophobic tunnel along helices 1 and 4. Under non-reducing conditions the tunnel is closed by two disulfide bridges (Cys3-Cys69' and Cys69-Cys3') that lie in the most flexible region of the dimer. Reduction or replacement of a cysteine residue enables conformational changes that open the channel allowing progesterone to enter.

Structure-function relationship for the highly toxic crotoxin from Crotalus durissus terrificus.

The three-dimensional structure of the highly toxic crotoxin from Crotalus durissus terrificus was modelled based on sequence analysis and the refined structure of calcium-free phospholipase of Crotalus atrox venom. Small-angle x-ray scattering experiments were performed on aqueous solutions of crotoxin. The radial distribution function derived from these scattering experiments and the one calculated from the model structure are in good agreement. Crotoxin consists of a basic and an acidic subunit. The model strongly suggests that the overall folding motif of phospholipases has been preserved in both subunits. The basic domain has an intact active site. The residues that are expected to contact the lipid tails of the phospholipid are different from other phospholipases, but they are all hydrophobic. The acidic domain consists of three independent chains interconnected by disulfide bonds. Compared to other phospholipases the active site for the greater part has been preserved in this domain, but it is not very well shielded from solvent. Most residues normally in contact with the lipid tails of the phospholipid are missing, which might explain the acidic subunit's lack of phospholipase activity. A homology between the third chain of the acidic domain and neurophysins suggests that the acidic domain may act as a chaperone for the basic domain.