Probing the role of packing specificity in protein design

By using a protein-design algorithm that quantitatively considers side-chain packing, the effect of specific steric constraints on protein design was assessed in the core of the streptococcal protein G β1 domain. The strength of packing constraints used in the design was varied, resulting in core sequences that reflected differing amounts of packing specificity. The structural flexibility and stability of several of the designed proteins were experimentally determined and showed a trend from well-ordered to highly mobile structures as the degree of packing specificity in the design decreased. This trend both demonstrates that the inclusion of specific packing interactions is necessary for the design of native-like proteins and defines a useful range of packing specificity for the design algorithm. In addition, an analysis of the modeled protein structures suggested that penalizing for exposed hydrophobic surface area can improve design performance.

Design and characterization of α-melanotropin peptide analogs cyclized through rhenium and technetium metal coordination

α-Melanocyte stimulating hormone (α-MSH) analogs, cyclized through site-specific rhenium (Re) and technetium (Tc) metal coordination, were structurally characterized and analyzed for their abilities to bind α-MSH receptors present on melanoma cells and in tumor-bearing mice. Results from receptor-binding assays conducted with B16 F1 murine melanoma cells indicated that receptor-binding affinity was reduced to approximately 1% of its original levels after Re incorporation into the cyclic Cys4,10, d-Phe7–α-MSH4-13 analog. Structural analysis of the Re–peptide complex showed that the disulfide bond of the original peptide was replaced by thiolate–metal–thiolate cyclization. A comparison of the metal-bound and metal-free structures indicated that metal complexation dramatically altered the structure of the receptor-binding core sequence. Redesign of the metal binding site resulted in a second-generation Re–peptide complex (ReCCMSH) that displayed a receptor-binding affinity of 2.9 nM, 25-fold higher than the initial Re–α-MSH analog. Characterization of the second-generation Re–peptide complex indicated that the peptide was still cyclized through Re coordination, but the structure of the receptor-binding sequence was no longer constrained. The corresponding 99mTc- and 188ReCCMSH complexes were synthesized and shown to be stable in phosphate-buffered saline and to challenges from diethylenetriaminepentaacetic acid (DTPA) and free cysteine. In vivo, the 99mTcCCMSH complex exhibited significant tumor uptake and retention and was effective in imaging melanoma in a murine-tumor model system. Cyclization of α-MSH analogs via 99mTc and 188Re yields chemically stable and biologically active molecules with potential melanoma-imaging and therapeutic properties.

Up-regulation of specific tyrosinase mRNAs in mouse melanomas with the c2j gene substituted for the wild-type tyrosinase allele: Utilization in design of syngeneic immunotherapy models

The expression of cell-specialization genes is likely to be changing in tumor cells as their differentiation declines. Functional changes in these genes might yield unusual peptide epitopes with anti-tumor potential and could occur without modification in the DNA sequence of the gene. Melanomas undergo a characteristic decline in melanization that may reflect altered contributions of key melanocytic genes such as tyrosinase. Quantitative reverse transcriptase–PCR of the wild-type (C) tyrosinase gene in transgenic (C57BL/6 strain) mouse melanomas has revealed a shift toward alternative splicing of the pre-mRNA that generated increased levels of the Δ1b and Δ1d mRNA splice variants. The spontaneous c2j albino mutation of tyrosinase (in the C57BL/6 strain) changes the pre-mRNA splicing pattern. In c2j/c2j melanomas, alternative splicing was again increased. However, while some mRNAs (notably Δ1b) present in C/C were obligatorily absent, others (Δ3 and Δ1d) were elevated. In c2j/c2j melanomas, the percentage of total tyrosinase transcripts attributable to Δ3 reached approximately 2-fold the incidence in c2j/c2j or C/C skin melanocytes. The percentage attributable to Δ1d rose to approximately 2-fold the incidence in c2j/c2j skin, and to 10-fold that in C/C skin. These differences provide a basis for unique mouse models in which the melanoma arises in skin grafted from a C/C or c2j/c2j transgenic donor to a transgenic host of the same or opposite tyrosinase genotype. Immunotherapy designs then could be based on augmenting those antigenic peptides that are novel or overrepresented in a tumor relative to the syngeneic host.

Network formation of lipid membranes: Triggering structural transitions by chain melting

Phospholipids when dispersed in excess water generally form vesicular membrane structures. Cryo-transmission and freeze-fracture electron microscopy are combined here with calorimetry and viscometry to demonstrate the reversible conversion of phosphatidylglycerol aqueous vesicle suspensions to a three-dimensional structure that consists of extended bilayer networks. Thermodynamic analysis indicates that the structural transitions arise from two effects: (i) the enhanced membrane elasticity accompanying the lipid state fluctuations on chain melting and (ii) solvent-associated interactions (including electrostatics) that favor a change in membrane curvature. The material properties of the hydrogels and their reversible formation offer the possibility of future applications, for example in drug delivery, the design of structural switches, or for understanding vesicle fusion or fission processes.

Synthetic genes for glycoprotein design and the elucidation of hydroxyproline-O-glycosylation codes

Design of hydroxyproline (Hyp)-rich glycoproteins (HRGPs) offers an approach for the structural and functional analysis of these wall components, which are broadly implicated in plant growth and development. HRGPs consist of multiple small repetitive “glycomodules” extensively O-glycosylated through the Hyp residues. The patterns of Hyp-O-glycosylation are putatively coded by the primary sequence as described by the Hyp contiguity hypothesis, which predicts contiguous Hyp residues to be attachment sites of small arabinooligosaccharides (1–5 Ara residues/Hyp); while clustered, noncontiguous Hyp residues are sites of arabinogalactan polysaccharide attachment. As a test, we designed two simple HRGPs as fusion proteins with green fluorescent protein. The first was a repetitive Ser-Hyp motif that encoded only clustered noncontiguous Hyp residues, predicted polysaccharide addition sites. The resulting glycoprotein had arabinogalactan polysaccharide O-linked to all Hyp residues. The second construct, based on the consensus sequence of a gum arabic HRGP, contained both arabinogalactan and arabinooligosaccharide addition sites and, as predicted, gave a product that contained both saccharide types. These results identify an O-glycosylation code of plants.

Design of fast enzymes by optimizing interaction potential in active site

The diffusional encounter between substrate and enzyme, and hence catalytic efficiency, can be enhanced by mutating charged residues on the surface of the enzyme. In this paper we present a simple method for screening such mutations. This is based on our earlier result that electrostatic enhancement of the enzyme-substrate binding rate constant can be accounted for just by the interaction potential within the active site. Assuming that catalytic and structural integrity is maintained, the catalytic efficiency can be optimized by surface charge mutations which lead to stronger interaction potential within the active site. Application of the screening method on superoxide dismutase shows that only charge mutations close to the active site will have practical effect on the catalytic efficiency. This rationalizes a large number of findings obtained in previous simulation and experimental studies.

The design, synthesis, and evaluation of molecules that enable or enhance cellular uptake: Peptoid molecular transporters

Certain proteins contain subunits that enable their active translocation across the plasma membrane into cells. In the specific case of HIV-1, this subunit is the basic domain Tat49–57 (RKKRRQRRR). To establish the optimal structural requirements for this translocation process, and thereby to develop improved molecular transporters that could deliver agents into cells, a series of analogues of Tat49–57 were prepared and their cellular uptake into Jurkat cells was determined by flow cytometry. All truncated and alanine-substituted analogues exhibited diminished cellular uptake, suggesting that the cationic residues of Tat49–57 play a principal role in its uptake. Charge alone, however, is insufficient for transport as oligomers of several cationic amino acids (histidine, lysine, and ornithine) are less effective than Tat49–57 in cellular uptake. In contrast, a 9-mer of l-arginine (R9) was 20-fold more efficient than Tat49–57 at cellular uptake as determined by Michaelis–Menton kinetic analysis. The d-arginine oligomer (r9) exhibited an even greater uptake rate enhancement (>100-fold). Collectively, these studies suggest that the guanidinium groups of Tat49–57 play a greater role in facilitating cellular uptake than either charge or backbone structure. Based on this analysis, we designed and synthesized a class of polyguanidine peptoid derivatives. Remarkably, the subset of peptoid analogues containing a six-methylene spacer between the guanidine head group and backbone (N-hxg), exhibited significantly enhanced cellular uptake compared to Tat49–57 and even to r9. Overall, a transporter has been developed that is superior to Tat49–57, protease resistent, and more readily and economically prepared.

Design of three-dimensional domain-swapped dimers and fibrous oligomers

Three-dimensional (3D) domain-swapped proteins are intermolecularly folded analogs of monomeric proteins; both are stabilized by the identical interactions, but the individual domains interact intramolecularly in monomeric proteins, whereas they form intermolecular interactions in 3D domain-swapped structures. The structures and conditions of formation of several domain-swapped dimers and trimers are known, but the formation of higher order 3D domain-swapped oligomers has been less thoroughly studied. Here we contrast the structural consequences of domain swapping from two designed three-helix bundles: one with an up-down-up topology, and the other with an up-down-down topology. The up-down-up topology gives rise to a domain-swapped dimer whose structure has been determined to 1.5 Å resolution by x-ray crystallography. In contrast, the domain-swapped protein with an up-down-down topology forms fibrils as shown by electron microscopy and dynamic light scattering. This demonstrates that design principles can predict the oligomeric state of 3D domain-swapped molecules, which should aid in the design of domain-swapped proteins and biomaterials.

Structure-based design of selective and potent G quadruplex-mediated telomerase inhibitors

The telomerase enzyme is a potential therapeutic target in many human cancers. A series of potent inhibitors has been designed by computer modeling, which exploit the unique structural features of quadruplex DNA. These 3,6,9-trisubstituted acridine inhibitors are predicted to interact selectively with the human DNA quadruplex structure, as a means of specifically inhibiting the action of human telomerase in extending the length of single-stranded telomeric DNA. The anilino substituent at the 9-position of the acridine chromophore is predicted to lie in a third groove of the quadruplex. Calculated relative binding energies predict enhanced selectivity compared with earlier 3,6-disubstituted compounds, as a result of this substituent. The ranking order of energies is in accord with equilibrium binding constants for quadruplex measured by surface plasmon resonance techniques, which also show reduced duplex binding compared with the disubstituted compounds. The 3,6,9-trisubstututed acridines have potent in vitro inhibitory activity against human telomerase, with EC50 values of up to 60 nM.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having α,β-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme’s catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

Structural basis for specificity and potency of a flavonoid inhibitor of human CDK2, a cell cycle kinase.

The central role of cyclin-dependent kinases (CDKs) in cell cycle regulation makes them a promising target for studying inhibitory molecules that can modify the degree of cell proliferation. The discovery of specific inhibitors of CDKs such as polyhydroxylated flavones has opened the way to investigation and design of antimitotic compounds. A novel flavone, (-)-cis-5,7-dihydroxyphenyl-8-[4-(3-hydroxy-1-methyl)piperidinyl] -4H-1-benzopyran-4-one hydrochloride hemihydrate (L868276), is a potent inhibitor of CDKs. A chlorinated form, flavopiridol, is currently in phase I clinical trials as a drug against breast tumors. We determined the crystal structure of a complex between CDK2 and L868276 at 2.33 angstroms resolution and refined to an Rfactor 20.3%. The aromatic portion of the inhibitor binds to the adenine-binding pocket of CDK2, and the position of the phenyl group of the inhibitor enables the inhibitor to make contacts with the enzyme not observed in the ATP complex structure. The analysis of the position of this phenyl ring not only explains the great differences of kinase inhibition among the flavonoid inhibitors but also explains the specificity of L868276 to inhibit CDK2 and CDC2.

Role of hydrophobic interactions and desolvation in determining the structural properties of a model alpha beta peptide.

Model AB, a 20-amino acid peptide that was designed to adopt an alpha beta tertiary structure stabilized by hydrophobic interactions between residues in adjacent helical and extended segments, exhibited large pKa shifts of several ionizable groups and slow hydrogen/deuterium exchange rates of nearly all the peptide amide groups [Butcher, D. J., Bruch, M. D. & Moe, G. T. (1995) Biopolymers 36, 109-120]. These properties, which depend on structure and hydration, are commonly observed in larger proteins but are quite unusual for small peptides. To identify which of several possible features of the peptide design are most important in determining these properties, several closely related analogs of Model AB were characterized by CD and NMR spectroscopy. The results show that hydrophobic interactions between adjacent helical and extended segments are structure-determining and have the additional effect of altering water-peptide interactions over much of the peptide surface. These results may have important implications for understanding mechanisms of protein folding and for the design of independently folding peptides.

Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution: implications for catalytic proton transfer and inhibitor design.

The three-dimensional structure of murine mitochondrial carbonic anhydrase V has been determined and refined at 2.45-A resolution (crystallographic R factor = 0.187). Significant structural differences unique to the active site of carbonic anhydrase V are responsible for differences in the mechanism of catalytic proton transfer as compared with other carbonic anhydrase isozymes. In the prototypical isozyme, carbonic anhydrase II, catalytic proton transfer occurs via the shuttle group His-64; carbonic anhydrase V has Tyr-64, which is not an efficient proton shuttle due in part to the bulky adjacent side chain of Phe-65. Based on analysis of the structure of carbonic anhydrase V, we speculate that Tyr-131 may participate in proton transfer due to its proximity to zinc-bound solvent, its solvent accessibility, and its electrostatic environment in the protein structure. Finally, the design of isozyme-specific inhibitors is discussed in view of the complex between carbonic anhydrase V and acetazolamide, a transition-state analogue. Such inhibitors may be physiologically important in the regulation of blood glucose levels.

Maintenance of pre-mRNA secondary structure by epistatic selection.

Linkage disequilibrium between polymorphisms in a natural population may result from various evolutionary forces, including random genetic drift due to sampling of gametes during reproduction, restricted migration between subpopulations in a subdivided population, or epistatic selection. In this report, we present evidence that the majority of significant linkage disequilibria observed in introns of the alcohol dehydrogenase locus (Adh) of Drosophila pseudoobscura are due to epistatic selection maintaining secondary structure of precursor mRNA (pre-mRNA). Based on phylogenetic-comparative analysis and a likelihood approach, we propose secondary structure models of Adh pre-mRNA for the regions of the adult intron and intron 2 where clustering of linkage disequilibria has been observed. Furthermore, we applied the likelihood ratio test to the phylogenetically predicted secondary structure in intron 1. In contrast to the other two structures, polymorphisms associated with the more conserved stem-loop structure of intron 1 are in low frequency, and linkage disequilibria have not been observed. These findings are qualitatively consistent with a model of compensatory fitness interactions. This model assumes that mutations disrupting pairing in a secondary structural element are individually deleterious if they destabilize a functionally important structure; a second "compensatory" mutation, however, may restabilize the structure and restore fitness.

Design and biological activities of L-163,191 (MK-0677): a potent, orally active growth hormone secretagogue.

A potent, orally active growth hormone (GH) secretagogue L-163,191 belonging to a recently synthesized structural class has been characterized. L-163,191 releases GH from rat pituitary cells in culture with EC50 = 1.3 +/- 0.09 nM and is mechanistically indistinguishable from the GH-releasing peptide GHRP-6 and the prototypical nonpeptide GH secretagogue L-692,429 but clearly distinguishable from the natural GH secretagogue, GH-releasing hormone. L-163,191 elevates GH in dogs after oral doses as low as 0.125 mg/kg and was shown to be specific in its release of GH without significant effect on plasma levels of aldosterone, luteinizing hormone, thyroxine, and prolactin after oral administration of 1 mg/kg. Only modest increases in cortisol were observed. Based on these properties, L-163,191 has been selected for clinical studies.