Advanced Protein Modeling Method: Benchmarking and Applications in Computer-Aided Drug Discovery

Development of homology modeling methods will remain an area of active research. These methods aim to develop and model increasingly accurate three-dimensional structures of yet uncrystallized therapeutically relevant proteins e.g. Class A G-Protein Coupled Receptors. Incorporating protein flexibility is one way to achieve this goal. Here, I will discuss the enhancement and validation of the ligand-steered modeling, originally developed by Dr. Claudio Cavasotto, via cross modeling of the newly crystallized GPCR structures. This method uses known ligands and known experimental information to optimize relevant protein binding sites by incorporating protein flexibility. The ligand-steered models were able to model, reasonably reproduce binding sites and the co-crystallized native ligand poses of the β2 adrenergic and Adenosine 2A receptors using a single template structure. They also performed better than the choice of template, and crude models in a small scale high-throughput docking experiments and compound selectivity studies. Next, the application of this method to develop high-quality homology models of Cannabinoid Receptor 2, an emerging non-psychotic pain management target, is discussed. These models were validated by their ability to rationalize structure activity relationship data of two, inverse agonist and agonist, series of compounds. The method was also applied to improve the virtual screening performance of the β2 adrenergic crystal structure by optimizing the binding site using β2 specific compounds. These results show the feasibility of optimizing only the pharmacologically relevant protein binding sites and applicability to structure-based drug design projects.

Structural homologies with ATP- and folate-binding enzymes in the crystal structure of folylpolyglutamate synthetase

Folylpolyglutamate synthetase, which is responsible for the addition of a polyglutamate tail to folate and folate derivatives, is an ATP-dependent enzyme isolated from eukaryotic and bacterial sources, where it plays a key role in the retention of the intracellular folate pool. Here, we report the 2.4-Å resolution crystal structure of the MgATP complex of the enzyme from Lactobacillus casei. The structural analysis reveals that folylpolyglutamate synthetase is a modular protein consisting of two domains, one with a typical mononucleotide-binding fold and the other strikingly similar to the folate-binding enzyme dihydrofolate reductase. We have located the active site of the enzyme in a large interdomain cleft adjacent to an ATP-binding P-loop motif. Opposite this site, in the C domain, a cavity likely to be the folate binding site has been identified, and inspection of this cavity and the surrounding protein structure suggests that the glutamate tail of the substrate may project into the active site. A further feature of the structure is a well defined Ω loop, which contributes both to the active site and to interdomain interactions. The determination of the structure of this enzyme represents the first step toward the elucidation of the molecular mechanism of polyglutamylation of folates and antifolates.

Crystal structure of chemically synthesized [N33A] stromal cell-derived factor 1α, a potent ligand for the HIV-1 “fusin” coreceptor

Stromal cell-derived factor-1α (SDF-1α ) is a member of the chemokine superfamily and functions as a growth factor and chemoattractant through activation of CXCR4/LESTR/Fusin, a G protein-coupled receptor. This receptor also functions as a coreceptor for T-tropic syncytium-inducing strains of HIV-1. SDF-1α antagonizes infectivity of these strains by competing with gp120 for binding to the receptor. The crystal structure of a variant SDF-1α ([N33A]SDF-1α ) prepared by total chemical synthesis has been refined to 2.2-Å resolution. Although SDF-1α adopts a typical chemokine β-β-β-α topology, the packing of the α-helix against the β-sheet is strikingly different. Comparison of SDF-1α with other chemokine structures confirms the hypothesis that SDF-1α may be either an ancestral protein from which all other chemokines evolved or the chemokine that is the least divergent from a primordial chemokine. The structure of SDF-1α reveals a positively charged surface ideal for binding to the negatively charged extracellular loops of the CXCR4 HIV-1 coreceptor. This ionic complementarity is likely to promote the interaction of the mobile N-terminal segment of SDF-1α with interhelical sites of the receptor, resulting in a biological response.

A peptide that inhibits hydroxyapatite growth is in an extended conformation on the crystal surface

Proteins play an important role in the biological mechanisms controlling hard tissue development, but the details of molecular recognition at inorganic crystal interfaces remain poorly characterized. We have applied a recently developed homonuclear dipolar recoupling solid-state NMR technique, dipolar recoupling with a windowless sequence (DRAWS), to directly probe the conformation of an acidic peptide adsorbed to hydroxyapatite (HAP) crystals. The phosphorylated hexapeptide, DpSpSEEK (N6, where pS denotes phosphorylated serine), was derived from the N terminus of the salivary protein statherin. Constant-composition kinetic characterization demonstrated that, like the native statherin, this peptide inhibits the growth of HAP seed crystals when preadsorbed to the crystal surface. The DRAWS technique was used to measure the internuclear distance between two 13C labels at the carbonyl positions of the adjacent phosphoserine residues. Dipolar dephasing measured at short mixing times yielded a mean separation distance of 3.2 ± 0.1 Å. Data obtained by using longer mixing times suggest a broad distribution of conformations about this average distance. Using a more complex model with discrete α-helical and extended conformations did not yield a better fit to the data and was not consistent with chemical shift analysis. These results suggest that the peptide is predominantly in an extended conformation rather than an α-helical state on the HAP surface. Solid-state NMR approaches can thus be used to determine directly the conformation of biologically relevant peptides on HAP surfaces. A better understanding of peptide and protein conformation on biomineral surfaces may provide design principles useful for the modification of orthopedic and dental implants with coatings and biological growth factors that are designed to enhance biocompatibility with surrounding tissue.

Crystal structure of the BTB domain from PLZF

The BTB domain (also known as the POZ domain) is an evolutionarily conserved protein–protein interaction motif found at the N terminus of 5–10% of C2H2-type zinc-finger transcription factors, as well as in some actin-associated proteins bearing the kelch motif. Many BTB proteins are transcriptional regulators that mediate gene expression through the control of chromatin conformation. In the human promyelocytic leukemia zinc finger (PLZF) protein, the BTB domain has transcriptional repression activity, directs the protein to a nuclear punctate pattern, and interacts with components of the histone deacetylase complex. The association of the PLZF BTB domain with the histone deacetylase complex provides a mechanism of linking the transcription factor with enzymatic activities that regulate chromatin conformation. The crystal structure of the BTB domain of PLZF was determined at 1.9 Å resolution and reveals a tightly intertwined dimer with an extensive hydrophobic interface. Approximately one-quarter of the monomer surface area is involved in the dimer intermolecular contact. These features are typical of obligate homodimers, and we expect the full-length PLZF protein to exist as a branched transcription factor with two C-terminal DNA-binding regions. A surface-exposed groove lined with conserved amino acids is formed at the dimer interface, suggestive of a peptide-binding site. This groove may represent the site of interaction of the PLZF BTB domain with nuclear corepressors or other nuclear proteins.

Crystal structures of the copper and nickel complexes of RNase A: Metal-induced interprotein interactions and identification of a novel copper binding motif

We report the crystal structures of the copper and nickel complexes of RNase A. The overall topology of these two complexes is similar to that of other RNase A structures. However, there are significant differences in the mode of binding of copper and nickel. There are two copper ions per molecule of the protein, but there is only one nickel ion per molecule of the protein. Significant changes occur in the interprotein interactions as a result of differences in the coordinating groups at the common binding site around His-105. Consequently, the copper- and nickel-ion-bound dimers of RNase A act as nucleation sites for generating different crystal lattices for the two complexes. A second copper ion is present at an active site residue His-119 for which all the ligands are from one molecule of the protein. At this second site, His-119 adopts an inactive conformation (B) induced by the copper. We have identified a novel copper binding motif involving the α-amino group and the N-terminal residues.

A link between protein structure and enzyme catalyzed hydrogen tunneling

We present evidence that the size of an active site side chain may modulate the degree of hydrogen tunneling in an enzyme-catalyzed reaction. Primary and secondary kH/kT and kD/kT kinetic isotope effects have been measured for the oxidation of benzyl alcohol catalyzed by horse liver alcohol dehydrogenase at 25°C. As reported in earlier studies, the relationship between secondary kH/kT and kD/kT isotope effects provides a sensitive probe for deviations from classical behavior. In the present work, catalytic efficiency and the extent of hydrogen tunneling have been correlated for the alcohol dehydrogenase-catalyzed hydride transfer among a group of site-directed mutants at position 203. Val-203 interacts with the opposite face of the cofactor NAD+ from the alcohol substrate. The reduction in size of this residue is correlated with diminished tunneling and a two orders of magnitude decrease in catalytic efficiency. Comparison of the x-ray crystal structures of a ternary complex of a high-tunneling (Phe-93 → Trp) and a low-tunneling (Val-203 → Ala) mutant provides a structural basis for the observed effects, demonstrating an increase in the hydrogen transfer distance for the low-tunneling mutant. The Val-203 → Ala ternary complex crystal structure also shows a hyperclosed interdomain geometry relative to the wild-type and the Phe-93 → Trp mutant ternary complex structures. This demonstrates a flexibility in interdomain movement that could potentially narrow the distance between the donor and acceptor carbons in the native enzyme and may enhance the role of tunneling in the hydride transfer reaction.

Inhibition of regulator of G protein signaling function by two mutant RGS4 proteins

Regulators of G protein signaling (RGS) proteins limit the lifetime of activated (GTP-bound) heterotrimeric G protein α subunits by acting as GTPase-activating proteins (GAPs). Mutation of two residues in RGS4, which, based on the crystal structure of RGS4 complexed with Giα1-GDP-AlF4−, directly contact Giα1 (N88 and L159), essentially abolished RGS4 binding and GAP activity. Mutation of another contact residue (S164) partially inhibited both binding and GAP activity. Two other mutations, one of a contact residue (R167M/A) and the other an adjacent residue (F168A), also significantly reduced RGS4 binding to Giα1-GDP-AlF4−, but in addition redirected RGS4 binding toward the GTPγS-bound form. These two mutant proteins had severely impaired GAP activity, but in contrast to the others behaved as RGS antagonists in GAP and in vivo signaling assays. Overall, these results are consistent with the hypothesis that the predominant role of RGS proteins is to stabilize the transition state for GTP hydrolysis. In addition, mutant RGS proteins can be created with an altered binding preference for the Giα-GTP conformation, suggesting that efficient RGS antagonists can be developed.

Structural details of an interaction between cardiolipin and an integral membrane protein

Anionic lipids play a variety of key roles in biomembrane function, including providing the immediate environment for the integral membrane proteins that catalyze photosynthetic and respiratory energy transduction. Little is known about the molecular basis of these lipid–protein interactions. In this study, x-ray crystallography has been used to examine the structural details of an interaction between cardiolipin and the photoreaction center, a key light-driven electron transfer protein complex found in the cytoplasmic membrane of photosynthetic bacteria. X-ray diffraction data collected over the resolution range 30.0–2.1 Å show that binding of the lipid to the protein involves a combination of ionic interactions between the protein and the lipid headgroup and van der Waals interactions between the lipid tails and the electroneutral intramembrane surface of the protein. In the headgroup region, ionic interactions involve polar groups of a number of residues, the protein backbone, and bound water molecules. The lipid tails sit along largely hydrophobic grooves in the irregular surface of the protein. In addition to providing new information on the immediate lipid environment of a key integral membrane protein, this study provides the first, to our knowledge, high-resolution x-ray crystal structure for cardiolipin. The possible significance of this interaction between an integral membrane protein and cardiolipin is considered.

Crystal structure of the Sec18p N-terminal domain

Yeast Sec18p and its mammalian orthologue N-ethylmaleimide-sensitive fusion protein (NSF) are hexameric ATPases with a central role in vesicle trafficking. Aided by soluble adapter factors (SNAPs), Sec18p/NSF induces ATP-dependent disassembly of a complex of integral membrane proteins from the vesicle and target membranes (SNAP receptors). During the ATP hydrolysis cycle, the Sec18p/NSF homohexamer undergoes a large-scale conformational change involving repositioning of the most N terminal of the three domains of each protomer, a domain that is required for SNAP-mediated interaction with SNAP receptors. Whether an internal conformational change in the N-terminal domains accompanies their reorientation with respect to the rest of the hexamer remains to be addressed. We have determined the structure of the N-terminal domain from Sec18p by x-ray crystallography. The Sec18p N-terminal domain consists of two β-sheet-rich subdomains connected by a short linker. A conserved basic cleft opposite the linker may constitute a SNAP-binding site. Despite structural variability in the linker region and in an adjacent loop, all three independent molecules in the crystal asymmetric unit have the identical subdomain interface, supporting the notion that this interface is a preferred packing arrangement. However, the linker flexibility allows for the possibility that other subdomain orientations may be sampled.

Crystal structure of truncated human apolipoprotein A-I suggests a lipid-bound conformation

The structure of truncated human apolipoprotein A-I (apo A-I), the major protein component of high density lipoprotein, has been determined at 4-Å resolution. The crystals comprise residues 44–243 (exon 4) of apo A-I, a fragment that binds to lipid similarly to intact apo A-I and that retains the lipid-bound conformation even in the absence of lipid. The molecule consists almost entirely of a pseudo-continuous, amphipathic α-helix that is punctuated by kinks at regularly spaced proline residues; it adopts a shape similar to a horseshoe of dimensions 125 × 80 × 40 Å. Four molecules in the asymmetric unit associate via their hydrophobic faces to form an antiparallel four-helix bundle with an elliptical ring shape. Based on this structure, we propose a model for the structure of apo A-I bound to high density lipoprotein.

The crystal structure of the Rev binding element of HIV-1 reveals novel base pairing and conformational variability

The crystal and molecular structure of an RNA duplex corresponding to the high affinity Rev protein binding element (RBE) has been determined at 2.1-Å resolution. Four unique duplexes are present in the crystal, comprising two structural variants. In each duplex, the RNA double helix consists of an annealed 12-mer and 14-mer that form an asymmetric internal loop consisting of G-G and G-A noncanonical base pairs and a flipped-out uridine. The 12-mer strand has an A-form conformation, whereas the 14-mer strand is distorted to accommodate the bulges and noncanonical base pairing. In contrast to the NMR model of the unbound RBE, an asymmetric G-G pair with N2-N7 and N1-O6 hydrogen bonding, is formed in each helix. The G-A base pairing agrees with the NMR structure in one structural variant, but forms a novel water-mediated pair in the other. A backbone flip and reorientation of the G-G base pair is required to assume the RBE conformation present in the NMR model of the complex between the RBE and the Rev peptide.

Crystal structure of the breakpoint cluster region-homology domain from phosphoinositide 3-kinase p85α subunit

Proteins such as the product of the breakpoint cluster region, chimaerin, and the Src homology 3-binding protein 3BP1, are GTPase activating proteins (GAPs) for members of the Rho subfamily of small GTP-binding proteins (G proteins or GTPases). A 200-residue region, named the breakpoint cluster region-homology (BH) domain, is responsible for the GAP activity. We describe here the crystal structure of the BH domain from the p85 subunit of phosphatidylinositol 3-kinase at 2.0 Å resolution. The domain is composed of seven helices, having a previously unobserved arrangement. A core of four helices contains most residues that are conserved in the BH family. Their packing suggests the location of a G-protein binding site. This structure of a GAP-like domain for small GTP-binding proteins provides a framework for analyzing the function of this class of molecules.

Crystal structure of the Holliday junction DNA in complex with a single RuvA tetramer

In the major pathway of homologous DNA recombination in prokaryotic cells, the Holliday junction intermediate is processed through its association with RuvA, RuvB, and RuvC proteins. Specific binding of the RuvA tetramer to the Holliday junction is required for the RuvB motor protein to be loaded onto the junction DNA, and the RuvAB complex drives the ATP-dependent branch migration. We solved the crystal structure of the Holliday junction bound to a single Escherichia coli RuvA tetramer at 3.1-Å resolution. In this complex, one side of DNA is accessible for cleavage by RuvC resolvase at the junction center. The refined junction DNA structure revealed an open concave architecture with a four-fold symmetry. Each arm, with B-form DNA, in the Holliday junction is predominantly recognized in the minor groove through hydrogen bonds with two repeated helix-hairpin-helix motifs of each RuvA subunit. The local conformation near the crossover point, where two base pairs are disrupted, suggests a possible scheme for successive base pair rearrangements, which may account for smooth Holliday junction movement without segmental unwinding.

Evolution of an enzyme active site: The structure of a new crystal form of muconate lactonizing enzyme compared with mandelate racemase and enolase

Muconate lactonizing enzyme (MLE), a component of the β-ketoadipate pathway of Pseudomonas putida, is a member of a family of related enzymes (the “enolase superfamily”) that catalyze the abstraction of the α-proton of a carboxylic acid in the context of different overall reactions. New untwinned crystal forms of MLE were obtained, one of which diffracts to better than 2.0-Å resolution. The packing of the octameric enzyme in this crystal form is unusual, because the asymmetric unit contains three subunits. The structure of MLE presented here contains no bound metal ion, but is very similar to a recently determined Mn2+-bound structure. Thus, absence of the metal ion does not perturb the structure of the active site. The structures of enolase, mandelate racemase, and MLE were superimposed. A comparison of metal ligands suggests that enolase may retain some characteristics of the ancestor of this enzyme family. Comparison of other residues involved in catalysis indicates two unusual patterns of conservation: (i) that the position of catalytic atoms remains constant, although the residues that contain them are located at different points in the protein fold; and (ii) that the positions of catalytic residues in the protein scaffold are conserved, whereas their identities and roles in catalysis vary.