Discovery of dual ZAP70 and Syk kinases inhibitors by docking into a rare C-helix-out conformation of Syk

The non-receptor tyrosine kinase Syk (spleen tyrosine kinase) is a pharmaceutical relevant target because its over-activation is observed in several autoimmune diseases, allergy, and asthma. Here we report the identification of two novel inhibitors of Syk by high-throughput docking into a rare C-helix-out conformation published recently. Interestingly, both compounds are slightly more active on ZAP70 (Zeta-chain-associated protein kinase 70), which is the kinase closest to Syk in the phylogenetic tree of human kinases. Taken together, the docking pose and experimental results suggest that the higher affinity of the inhibitors for ZAP70 than Syk originates from a more populated C-helix-out conformation in ZAP70. The latter observation is congruent with the 100-fold lower intrinsic activity of ZAP70 than Syk, as the C-helix-out conformation is inactive. The pharmacophore features of DFG-in, C-helix-out compounds are analyzed in relation to DFG-out inhibitors.     

Foundations of vectorial metabolism and osmochemistry

Chemical transformations, like osmotic translocations, are transport processes when looked at in detail. In chemiosmotic systems, the pathways of specific ligand conduction are spatially orientated through osmoenzymes and porters in which the actions of chemical group, electron and solute transfer occur as vectorial (or higher tensorial order) diffusion processes down gradients of total potential energy that represent real spatially-directed fields of force. Thus, it has been possible to describe classical bag-of-enzymes biochemistry as well as membrane biochemistry in terms of transport. But it would not have been possible to explain biological transport in terms of classical transformational biochemistry or chemistry. The recognition of this conceptual asymmetry in favour of transport has seemed to be upsetting to some biochemists and chemists; and they have resisted the shift towards thinking primarily in terms of the vectorial forces and co-linear displacements of ligands in place of their much less informative scalar products that correspond to the conventional scalar energies. Nevertheless, considerable progress has been made in establishing vectorial metabolism and osmochemistry as acceptable biochemical disciplines embracing transport and metabolism, and bioenergetics has been fundamentally transformed as a result.     

Use of restrained molecular dynamics in water to determine three-dimensional protein structure: prediction of the three-dimensional structure of Ecballium elaterium trypsin inhibitor II

Refinement of distance geometry (DG) structures of EETI-II (Heitz et al.: Biochemistry 28:2392-2398, 1989), a member of the squash family trypsin inhibitor, have been carried out by restrained molecular dynamics (RMD) in water. The resulting models show better side chain apolar/polar surface ratio and estimated solvation free energy than structures refined "in vacuo." The consistent lower values of residual NMR constraint violations, apolar/polar surface ratio, and solvation free energy for one of these refined structures allowed prediction of the 3D folding and disulfide connectivity of EETI-II. Except for the few first residues for which no NMR constraints were available, this computer model fully agreed with X-ray structures of CMTI-I (Bode et al.: FEBS Lett. 242:285-292, 1989) and EETI-II complexed with trypsin that appeared after the RMD simulation was completed. Restrained molecular dynamics in water is thus proved to be highly valuable for refinement of DG structures. Also, the successful use of apolar/polar surface ratio and of solvation free energy reinforce the analysis of Novotny et al. (Proteins 4:19-30, 1988) and shows that these criteria are useful indicators of correct versus misfolded models.     

Contribution of unusual Arginine-Arginine short-range interactions to stabilization and recognition in proteins

Although the majority of the ion pairs found in proteins consists of two charges of opposite sign, the observation of some unusual arrangements of two arginines led us to a search of such occurrences in the Brookhaven Protein Data Bank. We have found 41 Arginine-Arginine interactions with a C...C distance less than 5 å. Computer graphics analysis of these structures shows that most of the Arg-Arg pairs are found in the vicinity of the surface of the proteins, in an easily hydrated region. In order to determine which factors could stabilize such arrangements of species of similar charge, we have carried out AM1 semi-empirical calculations on a model of two guanidinium ions surrounded by several water molecules. The results show the existence of stable clusters with six or more water molecules, with distances between C atoms around 3 å. The bridging role of the water molecules is an important structural and energetic feature and we find bridges of two and three molecules between the guanidinium ions. These results are in good agreement with the structures found in our search of the experimental data. Enhancement of the electrostatic potential around these clusters, when compared to one of the guanidinium ions alone, is also demonstrated.     

Intrinsic pKas of ionizable residues in proteins: An explicit solvent calculation for lysozyme

Molecular dynamics simulations of triclinic hen egg white lysozyme in aqueous solution were performed to calculate the intrinsic pK_as of 14 ionizable residues. An all-atom model was used for both solvent and solute, and a single 180 ps simulation in conjunction with a Gaussian fluctuation analysis method was used. An advantage of the Gaussian fluctuation method is that it only requires a single simulation of the system in a reference state to calculate all the pK_as in the protein, in contrast to multiple simulations for the free energy perturbation method. pK_int shifts with respect to reference titratable residues were evaluated and compared to results obtained using a finite difference Poisson-Boltzmann (FDPB) method with a continuum solvent model; overall agreement with the direction of the shifts was generally observed, though the magnitude of the shifts was typically larger with the explicit solvent model. The contribution of the first solvation shell to the total charging free energies of the titratable groups was explicitly evaluated and found to be significant. Dielectric shielding between pairs of titratable groups was examined and found to be smaller than expected. The effect of the approximations used to treat the long-range interactions on the pK_int shifts is discussed. © 1994 Wiley-Liss, Inc.     

Macromolecular solvation energies derived from small molecule crystal morphology.

The morphology of small molecule crystals provides a model for evaluating surface solvation energies in a system with similar packing density to that observed for amino acid residues in proteins. The solvation energies associated with the transfer of methylene and carboxyl groups between vacuum and aqueous phases are estimated to be approx. $40 and -260 cal/A2, respectively, from an analysis of the morphology of succinic acid crystals. These solvation energies predict values for contact angles in reasonable agreement with measurements determined from macroscopic monolayer surfaces. Transfer free energies between vapor and water phases for a series of carboxylic acids are also predicted reasonably well by these solvation energies, provided the surface exposure of different groups is quantitated with the molecular surface area rather than the more traditional accessible surface area. In general, molecular surfaces and molecular surface areas are seen to have important advantages for characterizing the structure and energetics of macromolecular surfaces. Crystal faces of succinic acid with the lowest surface energies in aqueous solution are characteristically smooth. Increasing surface roughness and apolarity are associated with higher surface energies, which suggests an approach for modifying the surface properties of proteins and other macromolecules.     

The use of NMR chemical shifts to analyse the MD trajectories: simulation of bovine pancreatic trypsin inhibitor dynamics in water as a test case for solvent influences.

In this paper the NMR secondary chemical shifts, that are estimated from a set of 3D-structures, are compared with the observed ones to appraise the behaviour of a known x-ray diffraction structure (of the bovine pancreatic trypsin inhibitor protein) when various molecular dynamics are applied. The results of a 200 ps molecular dynamics under various conditions are analysed and different ways to modify the molecular dynamics are considered. With the purpose of avoiding the time-consuming explicit representation of the solvent (water) molecules, an attempt was made to understand the role of the solvent and to develop an implicit representation, which may be refined. A simulation of hydrophobic effects in an aqueous environment is also proposed which seems to provide a better approximation of the observed solution structure of the protein.     

Atomic and residue hydrophilicity in the context of folded protein structures

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes. © 1995 Wiley-Liss, Inc.     

Solvation of platinum anti-cancer drugs in methanol-water mixtures

Solubilities and transfer chemical potentials of carboplatin, cisplatin, iproplatin, and several related platinum complexes have been determined in methanol-water mixtures. the range of solvation behaviour is discussed in relation to possible oral administration of complexes of this type.     

Visualizing biomolecular electrostatics in virtual reality with UnityMol‐APBS

Virtual reality is a powerful tool with the ability to immerse a user within a completely external environment. This immersion is particularly useful when visualizing and analyzing interactions between small organic molecules, molecular inorganic complexes, and biomolecular systems such as redox proteins and enzymes. A common tool used in the biomedical community to analyze such interactions is the Adaptive Poisson‐Boltzmann Solver (APBS) software, which was developed to solve the equations of continuum electrostatics for large biomolecular assemblages. Numerous applications exist for using APBS in the biomedical community including analysis of protein ligand interactions and APBS has enjoyed widespread adoption throughout the biomedical community. Currently, typical use of the full APBS toolset is completed via the command line followed by visualization using a variety of two‐dimensional external molecular visualization software. This process has inherent limitations: visualization of three‐dimensional objects using a two‐dimensional interface masks important information within the depth component. Herein, we have developed a single application, UnityMol‐APBS, that provides a dual experience where users can utilize the full range of the APBS toolset, without the use of a command line interface, by use of a simple graphical user interface (GUI) for either a standard desktop or immersive virtual reality experience.