PubMed Health⌕ Search

Biomedical subjects

Jordi Villà-Freixa

Publications and source records attributed to Jordi Villà-Freixa.

4 recordsLinked to original sources

The role of residue stability in transient protein-protein interactions involved in enzymatic phosphate hydrolysis. A computational study.

Finding why protein-protein interactions (PPIs) are so specific can provide a valuable tool in a variety of fields. Statistical surveys of so-called transient complexes (like those relevant for signal transduction mechanisms) have shown a tendency of polar residues to participate in the interaction region. Following this scheme, residues in the unbound partners have to compete between interacting with water or interacting with other residues of the protein. On the other hand, several works have shown that the notion of active site electrostatic preorganization can be used to interpret the high efficiency in enzyme reactions. This preorganization can be related to the instability of the residues important for catalysis. In some enzymes, in addition, conformational changes upon binding to other proteins lead to an increase in the activity of the enzymatic partner. In this article the linear response approximation version of the semimacroscopic protein dipoles Langevin dipoles (PDLD/S-LRA) model is used to evaluate the stability of several residues in two phosphate hydrolysis enzymes upon complexation with their activating partners. In particular, the residues relevant for PPI and for phosphate hydrolysis in the CDK2/Cyclin A and Ras/GAP complexes are analyzed. We find that the evaluation of the stability of residues in these systems can be used to identify not only active site regions but it can also be used as a guide to locate "hot spots" for PPIs. We also show that conformational changes play a major role in positioning interfacing residues in a proper "energetic" orientation, ready to interact with the residues in the partner protein surface. Thus, we extend the preorganization theory to PPIs, extrapolating the results we obtained from the above-mentioned complexes to a more general case. We conclude that the correlation between stability of a residue in the surface and the likelihood that it participates in the interaction can be a general fact for transient PPIs.

Adenosine Triphosphate↗

Framework-based design of a new all-purpose molecular simulation application: the Adun simulator.

Here we present Adun, a new molecular simulator that represents a paradigm shift in the way scientific programs are developed. The traditional algorithm centric methods of scientific programming can lead to major maintainability and productivity problems when developing large complex programs. These problems have long been recognized by computer scientists; however, the ideas and techniques developed to deal with them have not achieved widespread adoption in the scientific community. Adun is the result of the application of these ideas, including pervasive polymorphism, evolutionary frameworks, and refactoring, to the molecular simulation domain. The simulator itself is underpinned by the Adun Framework, which separates the structure of the program from any underlying algorithms, thus giving a completely reusable design. The aims are twofold. The first is to provide a platform for rapid development and implementation of different simulation types and algorithms. The second is to decrease the learning barrier for new developers by providing a rigorous and well-defined structure. We present some examples on the use of Adun by performing simple free-energy simulations for the adiabatic charging of a single ion, using both free-energy perturbation and the Bennett's method. We also illustrate the power of the design by detailing the ease with which ASEP/MD, an elaborated mean field QM/MM method originally written in FORTRAN 90, was implemented into Adun.

Algorithms↗

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↗

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↗