Biomedical subjects

A Caflisch

Publications and source records attributed to A Caflisch.

25 records · Page 2Linked to original sources

PubMed1999

Computational ligand design.

A variety of computational tools that are used to assist drug design are reviewed. Particular emphasis is given to the limitations and merits of different methodologies. Recently, a number of general methods have been proposed for clustering compounds in classes of drug-like and non-drug-like molecules. The usefulness of this classification for drug design is discussed. The estimation of (relative) binding affinities is from a theoretical point of view the most challenging part of ligand design. We review three methods for the estimation of binding energies. Firstly, quantitative structure-activity relationships (QSAR) are presented. These have gained significantly from recent developments of experimental techniques for combinatorial synthesis and high-throughput screening as well as the use of powerful computational procedures like genetic algorithms and neural networks for the derivation of models. Secondly, empirical energy functions are shown to lead to more general models than standard QSAR, since they are fitted to a variety of complexes. They have been used recently with considerable success. Thirdly, we briefly outline free energy calculations based on molecular dynamics simulations, the method with the most sound theoretical foundation. Recent developments are reestablishing the interest in this approach. In the last part of this review structure-based ligand design programs are described. These are closely related to docking, with the difference that in design, unlike in most docking procedures, ligands are built on a fragment-by-fragment basis. Finally, a short description of our approach to computational combinatorial ligand design is given.

Drug Design↗

PubMed1996

Computational combinatorial ligand design: application to human alpha-thrombin.

A new method is presented for computer-aided ligand design by combinatorial selection of fragments that bind favorably to a macromolecular target of known three-dimensional structure. Firstly, the multiple-copy simultaneous-search procedure (MCSS) is used to exhaustively search for optimal positions and orientations of functional groups on the surface of the macromolecule (enzyme or receptor fragment). The MCSS minima are then sorted according to an approximated binding free energy, whose solvation component is expressed as a sum of separate electrostatic and nonpolar contributions. The electrostatic solvation energy is calculated by the numerical solution of the linearized Poisson-Boltzmann equation, while the nonpolar contribution to the binding free energy is assumed to be proportional to the loss in solvent-accessible surface area. The program developed for computational combinatorial ligand design (CCLD) allows the fast and automatic generation of a multitude of highly diverse compounds, by connecting in a combinatorial fashion the functional groups in their minimized positions. The fragments are linked as two atoms may be either fused, or connected by a covalent bond or a small linker unit. To avoid the combinatorial explosion problem, pruning of the growing ligand is performed according to the average value of the approximated binding free energy of its fragments. The method is illustrated here by constructing candidate ligands for the active site of human alpha-thrombin. The MCSS minima with favorable binding free energy reproduce the interaction patterns of known inhibitors. Starting from these fragments, CCLD generates a set of compounds that are closely related to high-affinity thrombin inhibitors. In addition, putative ligands with novel binding motifs are suggested. Probable implications of the MCSS-CCLD approach for the evolving scenario of drug discovery are discussed.

Binding Sites↗

PubMed1995

Acid and thermal denaturation of barnase investigated by molecular dynamics simulations.

The transition in barnase from the native state to a partially unfolded conformation has been studied by molecular dynamics simulations with explicit water molecules at 360 K and low pH(450 ps), and at 600 K and neutral pH (three simulations of 120, 250 and 200 ps each). The use of several simulations provides evidence that the results are not sensitive to initial conditions. To mimic low pH conditions, the acidic sidechains in barnase were neutralized and the two histidine residues were doubly protonated. Runs at 300 K showed that the solvated structures at low pH (300 ps) and neutral pH (310 ps) are very similar. The main structural differences involved the acidic residues, histidine residues, and the beta-turn connecting strands 4 and 5. When the temperature is raised to 360 K at low pH and to 600 K at neutral pH the barnase molecule begins to unfold. The molecule rapidly expands (Rg changes from 13.9 A to 15.3 A in 450 ps at 360 K and from 13.7 A to between 15.1 and 15.5 A in 120 ps at 600 K). However, the expansion is not uniform. In all the simulations, the chain termini, loops and the N-terminal parts of the main alpha-helix (helix 1) show a continuous and progressive unfolding. An essential step in the denaturation process is that the major alpha-helix (helix 1) separates from the beta-sheet; this is coupled to the exposure of the principal hydrophobic core, many of whose non-polar side chains become solvated by hydrogen-bonded water molecules. The barnase-water interaction energy improves during unfolding at the expense of the barnase self-energy. The deterioration of the intramolecular van der Waals energy suggests that the rupture of the tight packing during the initial unfolding phase contributes to the energy barrier of the denaturation process. The mutationally well-analyzed Asp8-Arg110-Asp12 double salt-bridge on the barnase surface is found to be marginally stable in the folded form in the simulations. A Poisson-Boltzmann calculation indicates that the salt-bridge is unstable; this is probably due to an overestimate of the solvation energy. A detailed analysis of the main hydrophobic core reveals that increase in solvent-accessible surface area and penetration of water molecules are simultaneous in the high-temperature simulation; at lower temperatures there is significant cavity formation and the entrance of the water molecules is somewhat delayed. The cavities occur in the neighborhood of the hydrophobic sidechains; the region formed by the sidechains of Val10, Leu14, Leu20, Tyr24, Ala74, Ile76 and Tyr90 is involved. The loosening of the core packing is coupled to an increase in the number of dihedral transitions.(ABSTRACT TRUNCATED AT 400 WORDS)

Bacterial Proteins↗

PubMed1994

Molecular dynamics simulation of protein denaturation: solvation of the hydrophobic cores and secondary structure of barnase.

The transition in barnase from the native state to a compact globule has been studied with high-temperature molecular dynamics simulations. A partial destruction of the alpha-helices and the outer strands of the beta-sheet is observed with water molecules replacing the hydrogen bonds of the secondary structural elements. Simultaneously, the main alpha-helix moves away from the beta-sheet and exposes the principal hydrophobic core, many of whose nonpolar side chains, beginning with the ones near the surface, become solvated by hydrogen-bonded water molecules. This step involves a significant increase in the solvent-exposed surface area; the resulting loss of stability due to the hydrophobic effect may be the major source of the activation barrier in the unfolding reaction. The detailed mechanism described here for the first stage of the denaturation of barnase, including the essential role of water molecules, is likely to be representative of protein denaturation, in general.

Bacillus↗

PubMed1993

Multiple copy simultaneous search and construction of ligands in binding sites: application to inhibitors of HIV-1 aspartic proteinase.

Rational ligand design is a complex problem that can be divided into three parts: the search for optimal positions and orientations of functional groups in the binding site, the connection of such positions to form candidate ligands, and the estimation of their binding constants. Approaches for addressing the first two parts of the problem are described in the present work. They are applied to the construction of peptide ligands in the binding site of the human immunodeficiency virus 1 (HIV-1) proteinase. The primary objective is to test the method by comparison of the results with the MVT-101 complex structure for which coordinates are available; the results obtained with the liganded and unliganded proteinase structure are used to examine the utility of the latter for binding studies. A secondary objective is to show how to find new inhibitor candidates. The multiple copy simultaneous search (MCSS) method is utilized to search for optimal positions and orientations of a set of functional groups. For peptide ligands, functional groups corresponding to the protein main chain (N-methylacetamide) and to protein side chains (e.g., methanol, ethyl guanidinium) are used. The resulting N-methylacetamide minima are connected to form hexapeptide main chains with a simple pseudoenergy function that permits a complete search of all possible ways of connecting the minima. Side chains are added to the main-chain candidates by application of the same pseudoenergy function to the appropriate functional group minima. A set of 15 hexapeptides with the sequence of MVT-101 is then minimized by a Monte Carlo scheme, which allows for escape from local minima. Comparison of the MCSS results with the structure of MVT-101 in the HIV-1 binding site showed that all of its functional group positions correspond (within 2.4 A) to some (usually more than one) MCSS minima. There were also many other low-energy MCSS minima which do not appear in any known inhibitors, e.g., methyl ammonium minima in the neighborhood of the catalytic aspartates. Among the 15 lowest minima are seven hexapeptides with the same main-chain orientation as the one found by X-ray crystallography for the inhibitor MVT-101 in the binding site and eight with the main chain oriented in the opposite direction; the latter tend to be more stable. [Addendum: These results are in agreement with recent high-resolution crystallographic data provided after the study was completed.(ABSTRACT TRUNCATED AT 400 WORDS)

Aspartic Acid Endopeptidases↗

PubMed1992

Monte Carlo docking of oligopeptides to proteins.

A new two-step procedure has been developed for the docking of flexible oligopeptide chains of unknown conformation to static proteins of known structure. In the first step positions and conformations are sampled and the association energy minimized starting from an approximate preselected docking position. The resulting conformations are further optimized in the second step by a Metropolis Monte Carlo minimization, which optimizes each of these structures. The method has been tested on the HIV-1 aspartic proteinase complex with an inhibitor, whose crystallographic structure is known at 2.3 A resolution. Furthermore, the application of this method to the docking of the hendecapeptide 58-68 of the influenza A virus matrix protein to the HLA-A2 molecule produced results which are in agreement with experimental observations in identifying side chains critical for T cell recognition and residues responsible of MHC protein binding.

Algorithms↗

PubMed1992

Monte Carlo minimization with thermalization for global optimization of polypeptide conformations in cartesian coordinate space.

A new minimization procedure for the global optimization in cartesian coordinate space of the conformational energy of a polypeptide chain is presented. The Metropolis Monte Carlo minimization is thereby supplemented by a thermalization process, which is initiated whenever a structure becomes trapped in an area containing closely located local minima in the conformational space. The method has been applied to the endogenous opioid pentapeptide methionine enkephalin. Five among 13 different starting conformations led to the same apparent global minimum of an in-house developed energy function, a type II' reverse turn, the central residues of which are Gly-3-Phe-4. A comparison between the ECEPP/2 global minimum conformation of methionine enkephalin and the apparent one achieved by the present method shows that minimum-energy conformations having a certain similarity can be generated by relatively different force fields.

Amino Acid Sequence↗