Ole Hvilsted Olsen

Structural and Evolutionary Relationships among Protein Tyrosine Phosphatase Domains

With the current access to the whole genomes of various organisms and the completion of the first draft of the human genome, there is a strong need for a structure-function classification of protein families as an initial step in moving from DNA databases to a comprehensive understanding of human biology. As a result of the explosion in nucleic acid sequence information and the concurrent development of methods for high-throughput functional characterization of gene products, the genomic revolution also promises to provide a new paradigm for drug discovery, enabling the identification of molecular drug targets in a significant number of human diseases. This molecular view of diseases has contributed to the importance of combining primary sequence data with three-dimensional structure and has increased the awareness of computational homology modeling and its potential to elucidate protein function. In particular, when important proteins or novel therapeutic targets are identified—like the family of protein tyrosine phosphatases (PTPs) (reviewed in reference 53)—a structure-function classification of such protein families becomes an invaluable framework for further advances in biomedical science. Here, we present a comparative analysis of the structural relationships among vertebrate PTP domains and provide a comprehensive resource for sequence analysis of phosphotyrosine-specific PTPs.

Protein secondary structure and homology by neural networks The α-helices in rhodopsin

Neural networks provide a basis for semiempirical studies of pattern matching between the primary and secondary structures of proteins. Networks of the perceptron class have been trained to classify the amino‐acid residues into two categories for each of three types of secondary feature: α‐helix or not, β‐sheet or not, and random coil or not. The explicit prediction for the helices in rhodopsin is compared with both electron microscopy results and those of the Chou‐Fasman method. A new measure of homology between proteins is provided by the network approach, which thereby leads to quantification of the differences between the primary structures of proteins.

Rational design of coagulation factor VIIa variants with substantially increased intrinsic activity

A trace amount of coagulation factor VII (FVII) circulates in the blood in the activated form, FVIIa (EC 3.4.21.21), formed by internal proteolysis. To avoid disseminated thrombus formation, FVIIa remains in a conformation with zymogen-like properties. Association with tissue factor (TF), locally exposed upon vascular injury, is necessary to render FVIIa biologically active and initiate blood clotting. We have designed potent mutants of FVIIa by replacing residues believed to function as determinants for the inherent zymogenicity. The TF-independent rate of factor X activation was dramatically improved, up to about 100-fold faster than that obtained with the wild-type enzyme and close to that of the FVIIa-soluble TF complex. The mutants appear to retain the substrate specificity of the parent enzyme and can be further stimulated by TF. Insights into the mechanism behind the increased activity of the mutants, presumably also pertinent to the TF-induced, allosteric stimulation of FVIIa activity, were obtained by studying their calcium dependence and the accessibility of the N terminus of the protease domain to chemical modification. The FVIIa analogues promise to offer a more efficacious treatment of bleeding episodes especially in hemophiliacs with inhibitory antibodies precluding conventional replacement therapy.

Ligand-Induced Conformational Changes: Improved Predictions of Ligand Binding Conformations and Affinities

A computational docking strategy using multiple conformations of the target protein is discussed and evaluated. A series of low molecular weight, competitive, nonpeptide protein tyrosine phosphatase inhibitors are considered for which the x-ray crystallographic structures in complex with protein tyrosine phosphatase 1B (PTP1B) are known. To obtain a quantitative measure of the impact of conformational changes induced by the inhibitors, these were docked to the active site region of various structures of PTP1B using the docking program FlexX. Firstly, the inhibitors were docked to a PTP1B crystal structure cocrystallized with a hexapeptide. The estimated binding energies for various docking modes as well as the RMS differences between the docked compounds and the crystallographic structure were calculated. In this scenario the estimated binding energies were not predictive inasmuch as docking modes with low estimated binding energies corresponded to relatively large RMS differences when aligned with the corresponding crystal structure. Secondly, the inhibitors were docked to their parent protein structures in which they were cocrystallized. In this case, there was a good correlation between low predicted binding energy and a correct docking mode. Thirdly, to improve the predictability of the docking procedure in the general case, where only a single target protein structure is known, we evaluate an approach which takes possible protein side-chain conformational changes into account. Here, side chains exposed to the active site were considered in their allowed rotamer conformations and protein models containing all possible combinations of side-chain rotamers were generated. To evaluate which of these modeled active sites is the most likely binding site conformation for a certain inhibitor, the inhibitors were docked against all active site models. The receptor rotamer model corresponding to the lowest estimated binding energy is taken as the top candidate. Using this protocol, correct inhibitor binding modes could successfully be discriminated from proposed incorrect binding modes. Moreover, the ranking of the estimated ligand binding energies was in good agreement with experimentally observed binding affinities.

PROPERTIES OF SPIN AND FLUORESCENT LABELS AT A RECEPTOR-LIGAND INTERFACE

Site-directed labeling was used to obtain local information on the binding interface in a receptor-ligand complex. As a model we have chosen the specific association of the extracellular part of tissue factor (sTF) and factor VIIa (FVIIa), the primary initiator of the blood coagulation cascade. Different spectroscopic labels were covalently attached to an engineered cysteine in position 140 in sTF, a position normally occupied by a Phe residue previously characterized as an important contributor to the sTF:FVIIa interaction. Two spin labels, IPSL [N-(1-oxyl-2,2,5, 5-tetramethyl-3-pyrrolidinyl)iodoacetamide] and MTSSL [(1-oxyl-2,2,5, 5-tetramethylpyrroline-3-methyl)methanethiosulfonate], and two fluorescent labels, IAEDANS [5-((((2-iodoacetyl)amino) ethyl)amino)naphthalene-1-sulfonic acid] and BADAN [6-bromoacetyl-2-dimethylaminonaphthalene], were used. Spectral data from electron paramagnetic resonance (EPR) and fluorescence spectroscopy showed a substantial change in the local environment of all labels when the sTF:FVIIa complex was formed. However, the interaction was probed differently by each label and these differences in spectral appearance could be attributed to differences in label properties such as size, polarity, and/or flexibility. Accordingly, molecular modeling data suggest that the most favorable orientations are unique for each label. Furthermore, line-shape simulations of EPR spectra and calculations based on fluorescence depolarization measurements provided additional details of the local environment of the labels, thereby confirming a tight protein-protein interaction between FVIIa and sTF when the complex is formed. The tightness of this local interaction is similar to that seen in the interior of globular proteins.

Augmented intrinsic activity of Factor VIIa by replacement of residues 305, 314, 337 and 374: evidence of two unique mutational mechanisms of activity enhancement.

Coagulation Factor VIIa (FVIIa) lacks the ability to spontaneously complete the conversion to a fully active enzyme after specific cleavage of an internal peptide bond (Arg152-Ile153) in the zymogen. Recently, several variants of FVIIa with enhanced intrinsic activity have been constructed. The in vitro characterization of these variants has shed light on molecular determinants that put restrictions on FVIIa in favour of a zymogen-like conformation and warrants continued efforts. Here we describe a new FVIIa variant with high intrinsic activity containing the mutations Leu305-->Val, Ser314-->Glu, Lys337-->Ala, and Phe374-->Tyr. The variant, called FVIIa(VEAY), processes a tripeptidyl substrate very efficiently because of an unprecedented, 5.5-fold lowering of the K(m) value. Together with a 4-fold higher substrate turnover rate this gives the variant a catalytic efficiency 22 times that of wild-type FVIIa, which is reflected in a considerably enhanced susceptibility to inhibition by antithrombin and other inhibitors. For instance, the affinity of FVIIa(VEAY) for the S1 probe and inhibitor p -aminobenzamidine is represented by an 8-fold lower K(i) value compared with that of FVIIa. Activation of Factor X in solution occurs about 10 times faster with FVIIa(VEAY) than with FVIIa, due virtually exclusively to an increased kcat value. The high activity of FVIIa(VEAY) is not accompanied by an increased burial of the N-terminus of the protease domain. A comparison of the kinetic parameters and molecular properties of FVIIa(VEAY) with those of the previously described mutant V158D/E296V/M298Q-FVIIa (FVIIa(IIa)), and the locations of the substitutions in the two variants, reveals what appear to be two profoundly different structural mechanisms dictating improvements in enzymic performance.

Ring-Shaped Quasi-Soliton Solutions to the Two- and Three-Dimensional Sine-Gordon Equation

Ring-shaped solitary wave solutions to the Sine-Gordon equation in two and three spatial dimensions are investigated by numerical computation. Each expanding wave exhibits a return effect. The reflection of the shrinking wave at the singularity at the center of the wave is investigated in a particular case. Collision experiments in numero for expanding and shrinking concentric ring waves show that the solutions possess quasisoliton properties. A Backlund transformation for the non-symmetric three-dimensional case is given.

Allosteric Activation of Coagulation Factor VIIa Visualized by Hydrogen Exchange

Coagulation factor VIIa (FVIIa) is a serine protease that, after binding to tissue factor (TF), plays a pivotal role in the initiation of blood coagulation. We used hydrogen exchange monitored by mass spectrometry to visualize the details of FVIIa activation by comparing the exchange kinetics of distinct molecular states, namely zymogen FVII, endoproteolytically cleaved FVIIa, TF-bound zymogen FVII, TF-bound FVIIa, and FVIIa in complex with an active site inhibitor. The hydrogen exchange kinetics of zymogen FVII and FVIIa are identical indicating highly similar solution structures. However, upon tissue factor binding, FVIIa undergoes dramatic structural stabilization as indicated by decreased exchange rates localized throughout the protease domain and in distant parts of the light chain, spanning across 50Å and revealing a concerted interplay between functional sites in FVIIa. The results provide novel insights into the cofactor-induced activation of this important protease and reveal the potential for allosteric regulation in the trypsin family of proteases.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. I. ligand-induced changes in the protein motions.

Activity of enzymes, such as protein tyrosine phosphatases (PTPs), is often associated with structural changes in the enzyme, resulting in selective and stereospecific reactions with the substrate. To investigate the effect of a substrate on the motions occurring in PTPs, we have performed molecular dynamics simulations of PTP1B and PTP1B complexed with a high-affinity peptide DADEpYL, where pY stands for phosphorylated tyrosine. The peptide sequence is derived from the epidermal growth factor receptor (EGFR988-993). Simulations were performed in water for 1 ns, and the concerted motions in the protein were analyzed using the essential dynamics technique. Our results indicate that the predominately internal motions in PTP1B occur in a subspace of only a few degrees of freedom. Upon substrate binding, the flexibility of the protein is reduced by approximately 10%. The largest effect is found in the protein region, where the N-terminal of the substrate is located, and in the loop region Val198-Gly209. Displacements in the latter loop are associated with the motions in the WPD loop, which contains a catalytically important aspartic acid. Estimation of the pKa of the active-site cysteine along the trajectory indicates that structural inhomogeneity causes the pKa to vary by approximately +/-1 pKa unit. In agreement with experimental observations, the active-site cysteine is negatively charged at physiological pH.

Molecular dynamics simulations of protein-tyrosine phosphatase 1B. II. substrate-enzyme interactions and dynamics.

Molecular dynamics simulations of protein tyrosine phosphatase 1B (PTP1B) complexed with the phosphorylated peptide substrate DADEpYL and the free substrate have been conducted to investigate 1) the physical forces involved in substrate-protein interactions, 2) the importance of enzyme and substrate flexibility for binding, 3) the electrostatic properties of the enzyme, and 4) the contribution from solvation. The simulations were performed for 1 ns, using explicit water molecules. The last 700 ps of the trajectories was used for analysis determining enthalpic and entropic contributions to substrate binding. Based on essential dynamics analysis of the PTP1B/DADEpYL trajectory, it is shown that internal motions in the binding pocket occur in a subspace of only a few degrees of freedom. In particular, relatively large flexibilities are observed along several eigenvectors in the segments: Arg(24)-Ser(28), Pro(38)-Arg(47), and Glu(115)-Gly(117). These motions are correlated to the C- and N-terminal motions of the substrate. Relatively small fluctuations are observed in the region of the consensus active site motif (H/V)CX(5)R(S/T) and in the region of the WPD loop, which contains the general acid for catalysis. Analysis of the individual enzyme-substrate interaction energies revealed that mainly electrostatic forces contribute to binding. Indeed, calculation of the electrostatic field of the enzyme reveals that only the field surrounding the binding pocket is positive, while the remaining protein surface is characterized by a predominantly negative electrostatic field. This positive electrostatic field attracts negatively charged substrates and could explain the experimentally observed preference of PTP1B for negatively charged substrates like the DADEpYL peptide.