Solvation studies of DMP323 and A76928 bound to HIV protease: analysis of water sites using grand canonical Monte Carlo simulations.

We examine the water solvation of the complex of the inhibitors DMP323 and A76928 bound to HIV-1 protease through grand canonical Monte Carlo simulations, and demonstrate the ability of this method to reproduce crystal waters and effectively predict water positions not seen in the DMP323 or A76928 structures. The simulation method is useful for identifying structurally important waters that may not be resolved in the crystal structures. It can also be used to identify water positions around a putative drug candidate docked into a binding pocket. Knowledge of these water positions may be useful in designing drugs to utilize them as bridging groups or displace them in the binding pocket. In addition, the method should be useful in finding water sites in homology models of enzymes for which crystal structures are unavailable.

Enzyme-inhibitor association thermodynamics: explicit and continuum solvent studies.

Studying the thermodynamics of biochemical association reactions at the microscopic level requires efficient sampling of the configurations of the reactants and solvent as a function of the reaction pathways. In most cases, the associating ligand and receptor have complementary interlocking shapes. Upon association, loosely connected or disconnected solvent cavities at and around the binding site are formed. Disconnected solvent regions lead to severe statistical sampling problems when simulations are performed with explicit solvent. It was recently proposed that, when such limitations are encountered, they might be overcome by the use of the grand canonical ensemble. Here we investigate one such case and report the association free energy profile (potential of mean force) between trypsin and benzamidine along a chosen reaction coordinate as calculated using the grand canonical Monte Carlo method. The free energy profile is also calculated for a continuum solvent model using the Poisson equation, and the results are compared to the explicit water simulations. The comparison shows that the continuum solvent approach is surprisingly successful in reproducing the explicit solvent simulation results. The Monte Carlo results are analyzed in detail with respect to solvation structure. In the binding site channel there are waters bridging the carbonyl oxygen groups of Asp189 with the NH2 groups of benzamidine, which are displaced upon inhibitor binding. A similar solvent-bridging configuration has been seen in the crystal structure of trypsin complexed with bovine pancreatic trypsin inhibitor. The predicted locations of other internal waters are in very good agreement with the positions found in the crystal structures, which supports the accuracy of the simulations.

Structure-based drug design: computational advances.

Structure-based computational methods continue to enhance progress in the discovery and refinement of therapeutic agents. Several such methods and their applications are described. These include molecular visualization and molecular modeling, docking, fragment methods, 3-D database techniques, and free-energy perturbation. Related issues that are discussed include the use of simplified potential energy functions and the determination of the positions of tightly bound waters. Strengths and weaknesses of the various methods are described.

Theoretical study of inhibition of adenosine deaminase by (8R)-coformycin and (8R)-deoxycoformycin.

Molecular dynamics and free energy simulations were performed to examine the binding of (8R)-deoxycoformycin and (8R)-coformycin to adenosine deaminase. The two inhibitors differ only at the 2' position of the sugar ring; the sugar moiety of conformycin is ribose, while it is deoxyribose for deoxycoformycin. The 100 ps molecular dynamics trajectories reveal that Asp 19 and His 17 interact strongly with the 5' hydroxyl group of the sugar moiety of both inhibitors and appear to play an important role in binding the sugar. The 2' and 3' groups of the sugars are near the protein-water interface and can be stabilized by either protein residues or water. The flexibility of the residues at the opening of the active site helps to explain the modest difference in binding of the two inhibitors and how substrates/inhibitors can enter an otherwise inaccessible binding site.

Computational science new horizons and relevance to pharmaceutical design.

Computer methods are used extensively in the design and refinement of drug leads. A short summary is given for several computational methods followed by a description of how some of these methods have been applied to design drugs targeted to the renin-angiotensin system and to cholinergic synapses. These methods include quantitative structure-activity relationship (QSAR) methods, comparative molecular field analyses (CoMFA), 3D database searching, de novo design of ligands, docking, and computational alchemy [free energy perturbation (FEP) and thermodynamic integration (MCTI)]. Most of these methods can be used whether or not detailed structural information about the binding site is available, although without an x-ray structure, the analyses are more qualitative. All of these methods are used extensively in the commercial design of pharmaceuticals. The main problem with most of these methods is in the scoring (ranking) of interactions or matches. Advances in this area and others (methods development and increases in capabilities of computers) will increase the predictive power of these methods and help to speed the time to market of new pharmaceuticals. (Trends Cardiovasc Med 1996;6:198-203).