Protein targets for structure-based drug design

Structural basis for isoform-selective inhibition in nitric oxide synthase

Nitric oxide synthase (NOS) converts l-arginine into l-citrulline and releases the important signaling molecule nitric oxide (NO). In the cardiovascular system, NO produced by endothelial NOS (eNOS) relaxes smooth muscle which controls vascular tone and blood pressure. Neuronal NOS (nNOS) produces NO in the brain, where it influences a variety of neural functions such as neural transmitter release. NO can also support the immune system, serving as a cytotoxic agent during infections. Even with all of these important functions, NO is a free radical and, when overproduced, it can cause tissue damage. This mechanism can operate in many neurodegenerative diseases, and as a result the development of drugs targeting nNOS is a desirable therapeutic goal. However, the active sites of all three human isoforms are very similar, and designing inhibitors specific for nNOS is a challenging problem. It is critically important, for example, not to inhibit eNOS owing to its central role in controlling blood pressure. In this Account, we summarize our efforts in collaboration with Rick Silverman at Northwestern University to develop drug candidates that specifically target NOS using crystallography, computational chemistry, and organic synthesis. As a result, we have developed aminopyridine compounds that are 3800-fold more selective for nNOS than eNOS, some of which show excellent neuroprotective effects in animal models. Our group has solved approximately 130 NOS-inhibitor crystal structures which have provided the structural basis for our design efforts. Initial crystal structures of nNOS and eNOS bound to selective dipeptide inhibitors showed that a single amino acid difference (Asp in nNOS and Asn in eNOS) results in much tighter binding to nNOS. The NOS active site is open and rigid, which produces few large structural changes when inhibitors bind. However, we have found that relatively small changes in the active site and inhibitor chirality can account for large differences in isoform-selectivity. For example, we expected that the aminopyridine group on our inhibitors would form a hydrogen bond with a conserved Glu inside the NOS active site. Instead, in one group of inhibitors, the aminopyridine group extends outside of the active site where it interacts with a heme propionate. For this orientation to occur, a conserved Tyr side chain must swing out of the way. This unanticipated observation taught us about the importance of inhibitor chirality and active site dynamics. We also successfully used computational methods to gain insights into the contribution of the state of protonation of the inhibitors to their selectivity. Employing the lessons learned from the aminopyridine inhibitors, the Silverman lab designed and synthesized symmetric double-headed inhibitors with an aminopyridine at each end, taking advantage of their ability to make contacts both inside and outside of the active site. Crystal structures provided yet another unexpected surprise. Two of the double-headed inhibitor molecules bound to each enzyme subunit, and one molecule participated in the generation of a novel Zn(2+) site that required some side chains to adopt alternate conformations. Therefore, in addition to achieving our specific goal, the development of nNOS selective compounds, we have learned how subtle differences in dynamics and structure can control protein-ligand interactions and often in unexpected ways.

Dihydrofolate reductase as a therapeutic target for infectious diseases: opportunities and challenges

Infectious diseases caused by parasites continue to take a massive toll on human health in the poor regions of the world. Filling the anti-infective drug-discovery pipeline has never been as challenging as it is now. The organisms responsible for these diseases have interesting biology with many potential biochemical targets. Inhibition of metabolic enzymes has been established as an attractive strategy for anti-infectious drug development. In this field, dihydrofolate reductase (DHFR) is an important enzyme in nucleic and amino acid synthesis and an extensively studied drug target over the past 50 years. The challenges for novel DHFR inhibition-based chemotherapeutics for the treatment of infectious diseases are now focused on overcoming the resistance problem as well as cost–effectiveness. Each year, the large number of literature citations attest the continued popularity of DHFR. It becomes truly the ‘enzyme of choice for all seasons and almost all reasons’. Herein, we summarize the opportunities and challenges in developing novel lead based on this target.

X-ray crystallography: Assessment and validation of protein-small molecule complexes for drug discovery

INTRODUCTION: Crystallography is the key initial component for structure-based and fragment-based drug design and can often generate leads that can be developed into high potency drugs. Therefore, huge sums of money are committed based on the outcome of crystallography experiments and their interpretation. AREAS COVERED: This review discusses how to evaluate the correctness of an X-ray structure, focusing on the validation of small molecule-protein complexes. Various types of inaccuracies found within the PDB are identified and the ramifications of these errors are discussed. The reader will gain an understanding of the key parameters that need to be inspected before a structure can be used in drug discovery efforts, as well as an appreciation of the difficulties of correctly interpreting electron density for small molecules. The reader will also be introduced to methods for validating small molecules within the context of a macromolecular structure. EXPERT OPINION: One of the reasons that ligand identification and positioning, within a macromolecular crystal structure, is so difficult is that the quality of small molecules widely varies in the PDB. For this reason, the PDB can not always be considered a reliable repository of structural information pertaining to small molecules, and this makes the derivation of general principles that govern small molecule-protein interactions more difficult.

Environmental effects on charge densities of biologically active molecules: do molecule crystal environments indeed approximate protein surroundings?

In the present paper, we investigate whether crystal and enzyme environments influence the electron density (ED) of active compounds in a similar manner. This supposition is essential for high-resolution X-ray studies, which use the EDs obtained from crystals of the pure active compound as approximations for the ED of the active compound in its complex with the target enzyme. The EDs of such complexes determine the molecular recognition process between the targeted enzyme and active compound and are, hence, extremely useful tools for rational drug design. The approximation of such EDs by data obtained from crystals of the pure active compound is needed since high-resolution X-ray experiments of the target-ligand complexes are still extremely demanding. Quantum mechanical/molecular mechanical (QM/MM) and pure QM calculations are employed to determine the EDs of two inhibitors, the reversible trans-4-(aminomethyl)cyclohexane-1-carboxylic acid (AMCHA) and the irreversible E64c in four different environments (the enzyme-inhibitor complex, crystals of the pure compounds, a continuum solvation model, and the gas phase). Our investigation shows that the environment inside of the crystal of the pure active compound generally influences the ED of an active compound in a very similar way as the enzyme surrounding in the complex between the active compound and target enzyme. However, this does not hold any more if the geometrical arrangement of the inhibitor in the enzyme differs significantly from that in the crystal. While EDs computed for gas-phase environments deviate strongly from those in crystal and protein surroundings, polar solvent environments provide rather similar electron distributions. Thus, such continuum solvation models are very well suited to compute density databases which are to be employed for the determination of the ED of macromolecules.

Countering Cooperative Effects in Protease Inhibitors Using Constrained β-Strand-Mimicking Templates in Focused Combinatorial Libraries

A major problem in de novo design of enzyme inhibitors is the unpredictability of the induced fit, with the shape of both ligand and enzyme changing cooperatively and unpredictably in response to subtle structural changes within a ligand. We have investigated the possibility of dampening the induced fit by using a constrained template as a replacement for adjoining segments of a ligand. The template preorganizes the ligand structure, thereby organizing the local enzyme environment. To test this approach, we used templates consisting of constrained cyclic tripeptides, formed through side chain to main chain linkages, as structural mimics of the protease-bound extended beta-strand conformation of three adjoining amino acid residues at the N- or C-terminal sides of the scissile bond of substrates. The macrocyclic templates were derivatized to a range of 30 structurally diverse molecules via focused combinatorial variation of nonpeptidic appendages incorporating a hydroxyethylamine transition-state isostere. Most compounds in the library were potent inhibitors of the test protease (HIV-1 protease). Comparison of crystal structures for five protease-inhibitor complexes containing an N-terminal macrocycle and three protease-inhibitor complexes containing a C-terminal macrocycle establishes that the macrocycles fix their surrounding enzyme environment, thereby permitting independent variation of acyclic inhibitor components with only local disturbances to the protease. In this way, the location in the protease of various acyclic fragments on either side of the macrocyclic template can be accurately predicted. This type of templating strategy minimizes the problem of induced fit, reducing unpredictable cooperative effects in one inhibitor region caused by changes to adjacent enzyme-inhibitor interactions. This idea might be exploited in template-based approaches to inhibitors of other proteases, where a beta-strand mimetic is also required for recognition, and also other protein-binding ligands where different templates may be more appropriate.

A fully integrated protein crystallization platform for small-molecule drug discovery

Structure-based drug discovery in the pharmaceutical industry benefits from cost-efficient methodologies that quickly assess the feasibility of specific, often refractory, protein targets to form well-diffracting crystals. By tightly coupling construct and purification diversity with nanovolume crystallization, the Structural Biology Group at Syrrx has developed such a platform to support its small-molecule drug-discovery program. During the past 18 months of operation at Syrrx, the Structural Biology Group has executed several million crystallization and imaging trials on over 400 unique drug-discovery targets. Here, key components of the platform, as well as an analysis of some experimental results that allowed for platform optimization, will be described.

Beta-strand mimicking macrocyclic amino acids: templates for protease inhibitors with antiviral activity

New amino acids are reported in which component macrocycles are constrained to mimic tripeptides locked in a beta-strand conformation. The novel amino acids involve macrocycles functionalized with both an N- and a C-terminus enabling addition of appendages at either end to modify receptor affinity, selectivity, or membrane permeability. We show that the cycles herein are effective templates within inhibitors of HIV-1 protease. Eleven compounds originating from such bifunctionalized cyclic templates are potent inhibitors of HIV-1 protease (Ki 0.3-50 nM; pH 6.5, I = 0.1 M). Unlike normal peptides comprising amino acids, five of these macrocycle-containing compounds are potent antiviral agents with sub-micromolar potencies (IC(50) 170-900 nM) against HIV-1 replication in human MT2 cells. The most active antiviral agents are the most lipophilic, with calculated values of LogD(6.5) > or = 4. All molecules have a conformationally constrained 17-membered macrocyclic ring that has been shown to structurally mimic a tripeptide segment (Xaa)-(Val/Ile)-(Phe/Tyr) of a peptide substrate in the extended conformation. The presence of two trans amide bonds and a para-substituted aromatic ring prevents intramolecular hydrogen bonds and fixes the macrocycle in the extended conformation. Similarly constrained macrocycles may be useful templates for the creation of inhibitors for the many other proteins and proteases that recognize peptide beta-strands.

Molecular mechanics calculations of proteins. Comparison of different energy minimization strategies

A general strategy for performing energy minimization of proteins using the SYBYL molecular modelling program has been developed. The influence of several variables including energy minimization procedure, solvation, dielectric function and dielectric constant have been investigated in order to develop a general method, which is capable of producing high quality protein structures. Avian pancreatic polypeptide (APP) and bovine pancreatic phospholipase A2 (BP PLA2) were selected for the calculations, because high quality X-ray structures exist and because all classes of secondary structure are represented in the structures. The energy minimized structures were evaluated relative to the corresponding X-ray structures. The overall similarity was checked by calculating RMS distances for all atom positions. Backbone conformation was checked by Ramachandran plots and secondary structure elements evaluated by the length on hydrogen bonds. The dimensions of active site in BP PLA2 is very dependent on electrostatic interactions, due to the presence of the positively charged calcium ion. Thus, the distances between calcium and the calcium-coordinating groups were used as a quality index for this protein. Energy minimized structures of the trimeric PLA2 from Indian cobra (N.n.n. PLA2) were used for assessing the impact of protein-protein interactions. Based on the above mentioned criteria, it could be concluded that using the following conditions: Dielectric constant epsilon = 4 or 20; a distance dependent dielectric function and stepwise energy minimization, it is possible to reproduce X-ray structures very accurately without including explicit solvent molecules.

Active-site-directed 3D database searching: pharmacophore extraction and validation of hits

Two new computational tools, PRO_PHARMEX and PRO_SCOPE, for use in active-site-directed searching of 3D databases are described. PRO_PHARMEX is a flexible, graphics-based program facilitating the extraction of pharmacophores from the active site of a target macromolecule. These pharmacophores can then be used to search a variety of databases for novel lead compounds. Such searches can often generate many 'hits' of varying quality. To aid the user in setting priorities for purchase, synthesis or testing, PRO_SCOPE can be used to dock molecules rapidly back into the active site and to assign them a score using an empirical scoring function correlated to the free energy of binding. To illustrate how these tools can add value to existing 3D database software, their use in the design of novel thrombin inhibitors is described.

N-peptidyl, O-acyl hydroxamates: comparison of the selective inhibition of serine and cysteine proteinases

Two series of N-aminoacyl, O-benzoyl hydroxamates were designed to investigate the influence of the substituted benzoyl residue on the hydrolytic stability and the reactivity of these potential inhibitors towards selected cysteine and serine proteinases. The inactivators react more rapidly with cysteine proteinases than with the serine enzymes tested. While Z-Phe-Gly-NHO-Nbz is the most reactive inhibitor of cathepsin L, inhibiting the target protein by a second order rate constant of 932.000 M-1 s-1, the bacterial serine proteinase thermitase is inhibited best by Z-Gly-Phe-NHO-Nbz, exhibiting a second-order rate constant of 1.170 M-1 s-1. Thiolsubtilisin, having the thiol-group as the reactive nucleophile instead of serine, exhibits specificity constants of the inactivation two orders of magnitude smaller than subtilisin. The degree of selectivity of the inhibitors relative to cathepsin B, cathepsin L, cathepsin S and papain varies up to two orders of magnitude with respect to their second order rate constant of inactivation. The inhibitory reactivity of these compounds varies only up to sixfold depending on the benzoyl substituent. Similarly, the rate constants for the hydrolytic decomposition of the compounds vary by a factor of about 6, suggesting that the structural and mechanistic features of the compounds which are responsible for decomposition as well as for the enzyme inhibition are the same. Comparing both reactions, the data allow the calculation of an acceleration factor of 2.4 x 10(10) for the inhibition of cathepsin L by its most effective inhibitor, clearly characterizing this enzyme inhibition reaction as enzyme-activated.

PRO_LIGAND: an approach to de novo molecular design. 4. Application to the design of peptides

In some instances, peptides can play an important role in the discovery of lead compounds. This paper describes the peptide design facility of the de novo drug design package, PRO_LIGAND. The package provides a unified framework for the design of peptides that are similar or complementary to a specified target. The approach uses single amino acid residues, selected from preconstructed libraries of different residues and conformations, and places them on top of predefined target interaction sites. This approach is a well-tested methodology for the design of organics but has not been used for peptides before. Peptides represent a difficulty because of their great conformational flexibility and a study of the advantages and disadvantages of this simple approach is an important step in the development of design tools. After a description of our general approach, a more detailed discussion of its adaptation to peptides is given. The method is then applied to the design of peptide-based inhibitors to HIV-1 protease and the design of structural mimics of the surface region of lysozyme. The results are encouraging and point the way towards further development of interaction site-based approaches for peptide design.

PRO-LIGAND: an approach to de novo molecular design. 1. Application to the design of organic molecules

An approach to de novo molecular design, PRO-LIGAND, has been developed that, in the environment of a large, integrated molecular design and simulation system, provides a unified framework for the generation of novel molecules which are either similar or complementary to a specified target. The approach is based on a methodology that has proved to be effective in other studies--placing molecular fragments upon target interaction sites-but incorporates many novel features such as the use of a rapid graph-theoretical algorithm for fragment placing, a generalised driver for structure generation which offers a large variety of fragment assembly strategies to the user and the pre-screening of library fragments. After a detailed description of the relevant modules of the package, PRO-LIGAND's efficacy in aiding rational drug design is demonstrated by its ability to design mimics of methotrexate and potential inhibitors for dihydrofolate reductase and HIV-1 protease.

Structure-based drug design: applications in immunopharmacology and immunosuppression

Structure-based drug design (SBDD) combines the power of many scientific disciplines, such as X-ray crystallography, nuclear magnetic resonance, medicinal chemistry, molecular modeling, biology, enzymology and biochemistry, in a functional paradigm of drug development. The current strength of SBDD lies in parlaying enzyme inhibitors into drugs, but a variety of technological advances over the past few years now makes it possible to address complex biological targets, such as those regulating immunosuppression and immunoactivation. Manual Navia and Debra Peattie discuss the SBDD paradigm and consider several of its achievements and challenges in immunopharmacology, particularly as these apply to the design of novel, potent immunosuppressants.

Structure-based drug design: applications in immunopharmacology and immunosuppression

Structure-based drug design (SBDD) combines the power of many scientific disciplines, such as X-ray crystallography, nuclear magnetic resonance, medicinal chemistry, molecular modeling, biology, enzymology and biochemistry, in a functional paradigm of drug development. The current strength of SBDD lies in parlaying enzyme inhibitors into drugs, but a variety of technological advances over the past few years now makes it possible to address complex biological targets, such as those regulating immunosuppression and immunoactivation. Manuel Navia and Debra Peattie discuss the SBDD paradigm and consider several of its achievements and challenges in immunopharmacology, particularly as these apply to the design of novel, potent immunosuppressants.