Computational method to reduce the search space for directed protein evolution

We introduce a computational method to optimize the in vitro evolution of proteins. Simulating evolution with a simple model that statistically describes the fitness landscape, we find that beneficial mutations tend to occur at amino acid positions that are tolerant to substitutions, in the limit of small libraries and low mutation rates. We transform this observation into a design strategy by applying mean-field theory to a structure-based computational model to calculate each residue's structural tolerance. Thermostabilizing and activity-increasing mutations accumulated during the experimental directed evolution of subtilisin E and T4 lysozyme are strongly directed to sites identified by using this computational approach. This method can be used to predict positions where mutations are likely to lead to improvement of specific protein properties.

A selective peroxisome proliferator-activated receptor δ agonist promotes reverse cholesterol transport

The peroxisome proliferator-activated receptors (PPARs) are dietary lipid sensors that regulate fatty acid and carbohydrate metabolism. The hypolipidemic effects of the fibrate drugs and the antidiabetic effects of the glitazone drugs in humans are due to activation of the α (NR1C1) and γ (NR1C3) subtypes, respectively. By contrast, the therapeutic potential of the δ (NR1C2) subtype is unknown, due in part to the lack of selective ligands. We have used combinatorial chemistry and structure-based drug design to develop a potent and subtype-selective PPARδ agonist, GW501516. In macrophages, fibroblasts, and intestinal cells, GW501516 increases expression of the reverse cholesterol transporter ATP-binding cassette A1 and induces apolipoprotein A1-specific cholesterol efflux. When dosed to insulin-resistant middle-aged obese rhesus monkeys, GW501516 causes a dramatic dose-dependent rise in serum high density lipoprotein cholesterol while lowering the levels of small-dense low density lipoprotein, fasting triglycerides, and fasting insulin. Our results suggest that PPARδ agonists may be effective drugs to increase reverse cholesterol transport and decrease cardiovascular disease associated with the metabolic syndrome X.

Minimizing a binding domain from protein A.

We present a systematic approach to minimizing the Z-domain of protein A, a three-helix bundle (59 residues total) that binds tightly (Kd = 10 nM) to the Fc portion of an immunoglobin IgG1. Despite the fact that all the contacts seen in the x-ray structure of the complex with the IgG are derived from residues in the first two helices, when helix 3 is deleted, binding affinity is reduced > 10(5)-fold (Kd > 1 mM). By using structure-based design and phage display methods, we have iteratively improved the stability and binding affinity for a two-helix derivative, 33 residues in length, such that it binds IgG1, with a Kd of 43 nM. This was accomplished by stepwise selection of random mutations from three regions of the truncated Z-peptide: the 4 hydrophobic residues from helix 1 and helix 2 that contacted helix 3 (the exoface), followed by 5 residues between helix 1 and helix 2 (the intraface), and lastly by 19 residues at or near the interface that interacts with Fc (the interface). As selected mutations from each region were compiled (12 in total), they led to progressive increases in affinity for IgG, and concomitant increases in alpha-helical content reflecting increased stabilization of the two-helix scaffold. Thus, by sequential increases in the stability of the structure and improvements in the quality of the intermolecular contacts, one can reduce larger binding domains to smaller ones. Such mini-protein binding domains are more amenable to synthetic chemistry and thus may be useful starting points for the design of smaller organic mimics. Smaller binding motifs also provide simplified and more tractable models for understanding determinants of protein function and stability.

Considerations on the folding topology and evolutionary origin of cadherin domains.

Cell-cell adhesion in zonula adherens and desmosomal junctions is mediated by cadherins, and recent crystal structures of the first domain from murine N-cadherin provide a plausible molecular basis for this adhesive action. A structure-based sequence analysis of this adhesive domain indicates that its fold is common to all extracellular cadherin domains. The cadherin folding topology is also shown to be similar to immunoglobulin-like domains and to other Greek-key beta-sandwich structures, as diverse as domains from plant cytochromes, bacterial cellulases, and eukaryotic transcription factors. Sequence similarities between cadherins and these other molecules are very low, however, and intron patterns are also different. On balance, independent origins for a favorable folding topology seem more likely than evolutionary divergence from an ancestor common to cadherins and immunoglobulins.

A unified theory of carcinogenesis based on order–disorder transitions in DNA structure as studied in the human ovary and breast

Fourier transform-infrared/statistics models demonstrate that the malignant transformation of morphologically normal human ovarian and breast tissues involves the creation of a high degree of structural modification (disorder) in DNA, before restoration of order in distant metastases. Order–disorder transitions were revealed by methods including principal components analysis of infrared spectra in which DNA samples were represented by points in two-dimensional space. Differences between the geometric sizes of clusters of points and between their locations revealed the magnitude of the order–disorder transitions. Infrared spectra provided evidence for the types of structural changes involved. Normal ovarian DNAs formed a tight cluster comparable to that of normal human blood leukocytes. The DNAs of ovarian primary carcinomas, including those that had given rise to metastases, had a high degree of disorder, whereas the DNAs of distant metastases from ovarian carcinomas were relatively ordered. However, the spectra of the metastases were more diverse than those of normal ovarian DNAs in regions assigned to base vibrations, implying increased genetic changes. DNAs of normal female breasts were substantially disordered (e.g., compared with the human blood leukocytes) as were those of the primary carcinomas, whether or not they had metastasized. The DNAs of distant breast cancer metastases were relatively ordered. These findings evoke a unified theory of carcinogenesis in which the creation of disorder in the DNA structure is an obligatory process followed by the selection of ordered, mutated DNA forms that ultimately give rise to metastases.

Toward optimized carbohydrate-based anticancer vaccines: Epitope clustering, carrier structure, and adjuvant all influence antibody responses to Lewisy conjugates in mice

The feasibility of using carbohydrate-based vaccines for the immunotherapy of cancer is being actively explored at the present time. Although a number of clinical trials have already been conducted with glycoconjugate vaccines, the optimal design and composition of the vaccines has yet to be determined. Among the candidate antigens being examined is Lewisy (Ley), a blood group-related antigen that is overexpressed on the majority of human carcinomas. Using Ley as a model for specificity, we have examined the role of epitope clustering, carrier structure, and adjuvant on the immunogenicity of Ley conjugates in mice. A glycolipopeptide containing a cluster of three contiguous Ley-serine epitopes and the Pam3Cys immunostimulating moiety was found to be superior to a similar construct containing only one Ley-serine epitope in eliciting antitumor cell antibodies. Because only IgM antibodies were produced by this vaccine, the effect on immunogenicity of coupling the glycopeptide to keyhole limpet hemocyanin was examined; although both IgM and IgG antibodies were formed, the antibodies reacted only with the immunizing structure. Reexamination of the clustered Ley-serine Pam3Cys conjugate with the adjuvant QS-21 resulted in the identification of both IgG and IgM antibodies reacting with tumor cells, thus demonstrating the feasibility of an entirely synthetic carbohydrate-based anticancer vaccine in an animal model.

PartsList: a web-based system for dynamically ranking protein folds based on disparate attributes, including whole-genome expression and interaction information

As the number of protein folds is quite limited, a mode of analysis that will be increasingly common in the future, especially with the advent of structural genomics, is to survey and re-survey the finite parts list of folds from an expanding number of perspectives. We have developed a new resource, called PartsList, that lets one dynamically perform these comparative fold surveys. It is available on the web at http://bioinfo.mbb.yale.edu/partslist and http://www.partslist.org. The system is based on the existing fold classifications and functions as a form of companion annotation for them, providing ‘global views’ of many already completed fold surveys. The central idea in the system is that of comparison through ranking; PartsList will rank the approximately 420 folds based on more than 180 attributes. These include: (i) occurrence in a number of completely sequenced genomes (e.g. it will show the most common folds in the worm versus yeast); (ii) occurrence in the structure databank (e.g. most common folds in the PDB); (iii) both absolute and relative gene expression information (e.g. most changing folds in expression over the cell cycle); (iv) protein–protein interactions, based on experimental data in yeast and comprehensive PDB surveys (e.g. most interacting fold); (v) sensitivity to inserted transposons; (vi) the number of functions associated with the fold (e.g. most multi-functional folds); (vii) amino acid composition (e.g. most Cys-rich folds); (viii) protein motions (e.g. most mobile folds); and (ix) the level of similarity based on a comprehensive set of structural alignments (e.g. most structurally variable folds). The integration of whole-genome expression and protein–protein interaction data with structural information is a particularly novel feature of our system. We provide three ways of visualizing the rankings: a profiler emphasizing the progression of high and low ranks across many pre-selected attributes, a dynamic comparer for custom comparisons and a numerical rankings correlator. These allow one to directly compare very different attributes of a fold (e.g. expression level, genome occurrence and maximum motion) in the uniform numerical format of ranks. This uniform framework, in turn, highlights the way that the frequency of many of the attributes falls off with approximate power-law behavior (i.e. according to V–b, for attribute value V and constant exponent b), with a few folds having large values and most having small values.

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.

The crystal structure of an N-terminal two-domain fragment of vascular cell adhesion molecule 1 (VCAM-1): a cyclic peptide based on the domain 1 C-D loop can inhibit VCAM-1-alpha 4 integrin interaction.

Vascular cell adhesion molecule 1 (VCAM-1) represents a structurally and functionally distinct class of immunoglobulin superfamily molecules that bind leukocyte integrins and are involved in inflammatory and immune functions. X-ray crystallography defines the three-dimensional structure of the N-terminal two-domain fragment that participates in ligand binding. Residues in domain 1 important for ligand binding reside in the C-D loop, which projects markedly from one face of the molecule near the contact between domains 1 and 2. A cyclic peptide that mimics this loop inhibits binding of alpha 4 beta 1 integrin-bearing cells to VCAM-1. These data demonstrate how crystallographic structural information can be used to design a small molecule inhibitor of biological function.

Definition of a metal-dependent/Li(+)-inhibited phosphomonoesterase protein family based upon a conserved three-dimensional core structure.

Inositol polyphosphate 1-phosphatase, inositol monophosphate phosphatase, and fructose 1,6-bisphosphatase share a sequence motif, Asp-Pro-(Ile or Leu)-Asp-(Gly or Ser)-(Thr or Ser), that has been shown by crystallographic and mutagenesis studies to bind metal ions and participate in catalysis. We compared the six alpha-carbon coordinates of this motif from the crystal structures of these three phosphatases and found that they are superimposable with rms deviations ranging from 0.27 to 0.60 A. Remarkably, when these proteins were aligned by this motif a common core structure emerged, defined by five alpha-helices and 11 beta-strands comprising 155 residues having rms deviations ranging from 1.48 to 2.66 A. We used the superimposed structures to align the sequences within the common core, and a distant relationship was observed suggesting a common ancestor. The common core was used to align the sequences of several other proteins that share significant similarity to inositol monophosphate phosphatase, including proteins encoded by fungal qa-X and qutG, bacterial suhB and cysQ (identical to amtA), and yeast met22 (identical to hal2). Evolutionary comparison of the core sequences indicate that five distinct branches exist within this family. These proteins share metal-dependent/Li(+)-sensitive phosphomonoesterase activity, and each predicted tree branch exhibits unique substrate specificity. Thus, these proteins define an ancient structurally conserved family involved in diverse metabolic pathways including inositol signaling, gluconeogenesis, sulfate assimilation, and possibly quinone metabolism. Furthermore, we suggest that this protein family identifies candidate enzymes to account for both the therapeutic and toxic actions of Li+ as it is used in patients treated for manic depressive disease.