Biomedical subjects

Hugo Gutiérrez-de-Terán

Publications and source records attributed to Hugo Gutiérrez-de-Terán.

8 recordsLinked to original sources

PubMed2006

Computational analysis of plasmepsin IV bound to an allophenylnorstatine inhibitor.

The plasmepsin proteases from the malaria parasite Plasmodium falciparum are attracting attention as putative drug targets. A recently published crystal structure of Plasmodium malariae plasmepsin IV bound to an allophenylnorstatine inhibitor [Clemente, J.C. et al. (2006) Acta Crystallogr. D 62, 246-252] provides the first structural insights regarding interactions of this family of inhibitors with plasmepsins. The compounds in this class are potent inhibitors of HIV-1 protease, but also show nM binding affinities towards plasmepsin IV. Here, we utilize automated docking, molecular dynamics and binding free energy calculations with the linear interaction energy LIE method to investigate the binding of allophenylnorstatine inhibitors to plasmepsin IV from two different species. The calculations yield excellent agreement with experimental binding data and provide new information regarding protonation states of active site residues as well as conformational properties of the inhibitor complexes.

Animals↗

PubMed2006

Inhibitor binding to the plasmepsin IV aspartic protease from Plasmodium falciparum.

Plasmepsin IV (Plm IV) is one of the aspartic proteases present in the food vacuole of the malaria parasite Plasmodium falciparum involved in host hemoglobin degradation by the parasite. Using a series of previously synthesized plasmepsin inhibitors [Ersmark, K., et al. (2005) J. Med. Chem. 48, 6090-106], we report here experimental data and theoretical analysis of their inhibitory activity toward Plm IV. All compounds share a 1,2-dihydroxyethylene unit as the transition state mimic. They possess symmetric P1 and P1' side chains and either a diacylhydrazine, a five-membered oxadiazole ring, or a retroamide at the P2 and P2' positions. Experimental binding affinities are compared to those predicted by the linear interaction energy (LIE) method and an empirical scoring function, using both a crystal structure and a homology model for the enzyme. Molecular dynamics (MD) simulations of the modeled complexes allow a rational interpretation of the structural determinants for inhibitor binding. A ligand bearing a P2 and P2' symmetric oxadiazole which is devoid of amide bonds is identified both experimentally and theoretically as the most potent inhibitor of Plm IV. For the P2 and P2' asymmetric compounds, the results are consistent with earlier predictions regarding the mode of binding of this class of inhibitors to Plm II. Theoretical estimation of selectivity for some compounds is also reported. Significant features of the Plm IV binding pocket are discussed in comparison to related enzymes, and the results obtained here should be helpful for further optimization of inhibitors.

Animals↗

PubMed2005

Macrocyclic inhibitors of the malarial aspartic proteases plasmepsin I, II, and IV.

The first macrocyclic inhibitor of the Plasmodium falciparum aspartic proteases plasmepsin I, II, and IV with considerable selectivity over the human aspartic protease cathepsin D has been identified. A series of macrocyclic compounds were designed and synthesized. Cyclizations were accomplished using ring-closing metathesis with the second generation Grubbs catalyst. These compounds contain either a 13-membered or a 16-membered macrocycle and incorporate a 1,2-dihydroxyethylene as transition state mimicking unit. The binding mode of this new class of compounds was predicted with automated docking and molecular dynamics simulations, with an estimation of the binding affinities through the linear interaction energy (LIE) method.

Animals↗

PubMed2004

Towards a MIP-based alignment and docking in computer-aided drug design.

Structural alignment of ligands in their biological conformation is a crucial step in the building of pharmacophoric models in structure-based drug design. In addition, docking algorithms are limited in some cases by the quality of the scoring functions and the limited flexibility of the environment that the different programs allow. On the other hand, GRID molecular interaction potentials (MIPs) have been used for a long time in 3D-QSAR studies. However, in most of these studies the alignment of the molecules is performed on the basis of geometrical or physico-chemical criteria that differ from the MIPs used in the partial least squares statistical analysis. We have previously developed a method to use the same scoring function for the molecular alignment and for 3D-QSAR studies. This methodology, based on the use of GRID potentials, consists in the weighted averaging of similarities of the relevant MIPs of the molecules to be aligned. Here we present a method to obtain the weights for the different GRID probes in the average based on the structural information on protein-ligand complexes for relevant systems. The method, implemented in MIPSIM, is shown to yield good accuracy in the prediction of the alignments for two systems: a set of three inhibitors of dihydrofolate reductase and a set of fifteen non-nucleoside HIV-1 reverse transcriptase inhibitors (NNRTIs). The smooth GRID potentials are shown to capture the flexible character of the active site, as opposed to traditional docking scoring energy functions.

Binding Sites↗

PubMed2004

Comparative analysis of putative agonist-binding modes in the human A1 adenosine receptor.

A recent study reported a model of the human A(1) adenosine receptor and its agonist binding site, proposing two putative binding modes in the same binding site for the natural agonist, adenosine. The present work investigates the flexibility of this binding site by exhaustive exploration with the natural agonist and with three other adenosine derivatives: N6-cyclopentyladenosine (CPA), 2-chloro-N6-cyclopentyladenosine (CCPA), and 5'-N-ethylcarboxamidoadenosine (NECA). Our aim was to find a common binding mode for agonists that would explain the role in the binding process of the different substitutions allowed at the 2, N6, and 5' positions of adenosine. This problem was addressed through docking simulations, molecular dynamics studies, and estimations of the ligand-binding free energy with both the AUTODOCK scoring function and the linear interaction energy (LIE) approach. The results point to a single receptor-binding position that explains the effects of the different chemical modifications on the adenosine derivatives considered here.

Adenosine↗

PubMed2004

Novel approaches for modeling of the A1 adenosine receptor and its agonist binding site.

The present work describes the building of a human A(1) adenosine receptor (hA(1)AR) model, based on the X-ray crystal structure of bovine rhodopsin, and its use as a basis for the investigation of some important structural characteristics of the receptor. One of the issues investigated was the protonation position of two histidine residues known to influence ligand binding, with protonation of His251 (6.52) in epsilon position and His278 (7.43) in delta position showing the best agreement with experimental evidence. The model was also used to study the position and structural role of water molecules present in the helical bundle. Finally, the binding site location and the ligand docking were investigated using an objective strategy. A suitable site for the binding of the ribose moiety of adenosine was first postulated and further confirmed by means of a novel chemometric strategy based on GRIND descriptors. Using this position as an anchor point, the binding of adenosine was studied by docking and molecular dynamics simulations obtaining two putative binding positions in good agreement with experimental data.

Adenosine↗

PubMed2003

On the generation of catalytic antibodies by transition state analogues.

The effective design of catalytic antibodies represents a major conceptual and practical challenge. It is implicitly assumed that a proper transition state analogue (TSA) can elicit a catalytic antibody (CA) that will catalyze the given reaction in a similar way to an enzyme that would evolve (or was evolved) to catalyze this reaction. However, in most cases it was found that the TSA used produced CAs with relatively low rate enhancement as compared to the corresponding enzymes, when these exist. The present work explores the origin of this problem, by developing two approaches that examine the similarity of the TSA and the corresponding transition state (TS). These analyses are used to assess the proficiency of the CA generated by the given TSA. Both approaches focus on electrostatic effects that have been found to play a major role in enzymatic reactions. The first method uses molecular interaction potentials to look for the similarity between the TSA and the TS and, in principle, to help in designing new haptens by using 3D quantitative structure-activity relationships. The second and more quantitative approach generates a grid of Langevin dipoles, which are polarized by the TSA, and then uses the grid to bind the TS. Comparison of the resulting binding energy with the binding energy of the TS to the grid that was polarized by the TS provides an estimate of the proficiency of the given CA. Our methods are used in examining the origin of the difference between the catalytic power of the 1F7 CA and chorismate mutase. It is demonstrated that the relatively small changes in charge and structure between the TS and TSA are sufficient to account for the difference in proficiency between the CA and the enzyme. Apparently the environment that was preorganized to stabilize the TSA charge distribution does not provide a sufficient stabilization to the TS. The general implications of our findings and the difficulties in designing a perfect TSA are discussed. Finally, the possible use of our approach in screening for an optimal TSA is pointed out.

Antibodies, Catalytic↗

PubMed2003

Use of alignment-free molecular descriptors in diversity analysis and optimal sampling of molecular libraries.

The selection of a sample of diverse compounds is a common strategy for exploring large molecular libraries. However, the success of such approach depends on the selection of relevant molecular descriptors and the use of appropriate sampling methods. In the context of pharmaceutical research, the molecular descriptors should be based on physicochemical properties related with the pharmacological behaviour of the compounds. In this sense, the alignment-free GRIND and VolSurf molecular descriptors are promising candidates since they have been successfully used in the modelling of both pharmacodynamic and pharmacokinetic properties of drugs. This work describes the use of such descriptors in the diversity sampling of a library of primary amines and compares the results with those obtained in a previous study that used quantum-mechanical descriptors. As in the previous work, principal component (PC) analysis was applied to reduce the dimensionality and remove redundant information of the original descriptors, and the compounds were sampled on the basis of k-means clustering on the space of the selected PCs. The results of the present study show that VolSurf and GRIND provide similar quality sampling regarding global features of the molecules such as hydrophilicity, however the topology of the compounds is considered differently. The similarity between particular compounds strongly depends on the original descriptors used. However all the sample selections done in the PC space after k-means clustering provide the same apparent diversity in comparison to the whole dataset. The results indicate that there is no best set of descriptors on a diversity basis. The selection of descriptors must be based on the drug features to be investigated.

Combinatorial Chemistry Techniques↗