Mesoscope: A Web-based Tool for Mesoscale Data Integration and Curation.

Interest is growing for 3D models of the biological mesoscale, the intermediate scale between the nanometer scale of molecular structure and micrometer scale of cellular biology. However, it is currently difficult to gather, curate and integrate all the data required to define such models. To address this challenge we developed Mesoscope (mesoscope.scripps.edu/beta), a web-based data integration and curation tool. Mesoscope allows users to begin with a listing of molecules (such as data from proteomics), and to use resources at UniProt and the PDB to identify, prepare and validate appropriate structures and representations for each molecule, ultimately producing a portable output file used by CellPACK and other modeling tools for generation of 3D models of the biological mesoscale. The availability of this tool has proven essential in several exploratory applications, given the high complexity of mesoscale models and the heterogeneity of the available data sources.

Cuttlefish: Color Mapping for Dynamic Multi-Scale Visualizations.

Visualizations of hierarchical data can often be explored interactively. For example, in geographic visualization, there are continents, which can be subdivided into countries, states, counties and cities. Similarly, in models of viruses or bacteria at the highest level are the compartments, and below that are macromolecules, secondary structures (such as α-helices), amino-acids, and on the finest level atoms. Distinguishing between items can be assisted through the use of color at all levels. However, currently, there are no hierarchical and adaptive color mapping techniques for very large multi-scale visualizations that can be explored interactively. We present a novel, multi-scale, color-mapping technique for adaptively adjusting the color scheme to the current view and scale. Color is treated as a resource and is smoothly redistributed. The distribution adjusts to the scale of the currently observed detail and maximizes the color range utilization given current viewing requirements. Thus, we ensure that the user is able to distinguish items on any level, even if the color is not constant for a particular feature. The coloring technique is demonstrated for a political map and a mesoscale structural model of HIV. The technique has been tested by users with expertise in structural biology and was overall well received.

Analysis of large-scale sequencing of small RNAs.

The advent of large-scale sequencing has opened up new areas of research, such as the study of Piwi-interacting small RNAs (piRNAs). piRNAs are longer than miRNAs, close to 30 nucleotides in length, involved in various functions, such as the suppression of transposons in germline. Since a large number of them (many tens of thousands) are generated from a wide range of positions in the genome, large-scale sequencing is the only way to study them. The key to understanding their genesis and biological roles is efficient analysis, which is complicated by the large volumes of sequence data. Taking account of the underlying biology is also important. We describe here novel analyses techniques and tools applied to small RNAs from germ cells in D. melanogaster, that allowed us to infer mechanism and biological function.

Recognition templates for predicting adenylate-binding sites in proteins.

Recognition templates encapsulate the structural and energetic features for the specific recognition of a given ligand by a protein active site. These templates identify the major interactions used for specific recognition and may be used to find specific binding sites in proteins of unknown function. We present a grid-based method for deriving recognition templates for adenylate groups from a set of diverse nucleotide-binding proteins. The templates reveal the basis of specific binding of adenylate, including tight shape complementarity, specific hydrogen bonds, and underscoring the importance of a key steric contact for excluding guanylate from adenylate-specific sites. We demonstrate the utility of recognition templates in identifying specific adenylate-binding sites in a diverse set of dinucleotide-binding proteins.

Viral evolution in response to the broad-based retroviral protease inhibitor TL-3.

TL-3 is a protease inhibitor developed using the feline immunodeficiency virus protease as a model. It has been shown to efficiently inhibit replication of human, simian, and feline immunodeficiency viruses and therefore has broad-based activity. We now demonstrate that TL-3 efficiently inhibits the replication of 6 of 12 isolates with confirmed resistance mutations to known protease inhibitors. To dissect the spectrum of molecular changes in protease and viral properties associated with resistance to TL-3, a panel of chronological in vitro escape variants was generated. We have virologically and biochemically characterized mutants with one (V82A), three (M46I/F53L/V82A), or six (L24I/M46I/F53L/L63P/V77I/V82A) changes in the protease and structurally modeled the protease mutant containing six changes. Virus containing six changes was found to be 17-fold more resistant to TL-3 in cell culture than was wild-type virus but maintained similar in vitro replication kinetics compared to the wild-type virus. Analyses of enzyme activity of protease variants with one, three, and six changes indicated that these enzymes, compared to wild-type protease, retained 40, 47, and 61% activity, respectively. These results suggest that deficient protease enzymatic activity is sufficient for function, and the observed protease restoration might imply a selective advantage, at least in vitro, for increased protease activity.

Analysis of a data set of paired uncomplexed protein structures: new metrics for side-chain flexibility and model evaluation.

We compiled and analyzed a data set of paired protein structures containing proteins for which multiple high-quality uncomplexed atomic structures were available in the Protein Data Bank. Side-chain flexibility was quantified, yielding a set of residue- and environment-specific confidence levels describing the range of motion around chi1 and chi2 angles. As expected, buried residues were inflexible, adopting similar conformations in different crystal structure analyses. Ile, Thr, Asn, Asp, and the large aromatics also showed limited flexibility when exposed on the protein surface, whereas exposed Ser, Lys, Arg, Met, Gln, and Glu residues were very flexible. This information is different from and complementary to the information available from rotamer surveys. The confidence levels are useful for assessing the significance of observed side-chain motion and estimating the extent of side-chain motion in protein structure prediction. We compare the performance of a simple 40 degrees threshold with these quantitative confidence levels in a critical evaluation of side-chain prediction with the program SCWRL.

Transmembrane alpha-helices in the gap junction membrane channel: systematic search of packing models based on the pair potential function.

Recent progress in the field of electron cryo-microscopy and image analysis has shown that there is an overwhelming need to interpret medium resolution (5 to 10 A) three-dimensional maps. Traditional methods of fitting amino acid residues into electron density using molecular modeling programs must be supplemented with further analysis. We have used a potential of mean force (PMF) method, derived from Boltzmann statistics in protein structure, to generate models for the packing of alpha-helices, using pairwise potentials between amino acid residues. The approach was tested using the three-dimensional map of a recombinant cardiac gap junction membrane channel provided by electron cryo-crystallography (Unger et al., 1997; 1999a, 1999b) which had a resolution of 7.5 A in the membrane plane and 21 A in the vertical direction. The dodecameric channel was formed by the end-to-end docking of two hexamers, each of which displayed 24 rods of density in the membrane interior, which was consistent with an alpha-helical conformation for the four transmembrane domains of each connexin subunit. Based on the three-dimensional map and the amino acid sequence for the 4 transmembrane domains determined by hydropathy analysis, we used the modeling utility SymServ (Macke et al., 1998) to build hexameric connexons with 24 transmembrane alpha-helices. Canonical alpha-helices were aligned to the axes of the rods of density and translated along the density so that the center of masses coincided. The PMF function was used to evaluate 162,000 conformations for each of the 24 possible alpha-helical packing models. Since the different packing models yielded different energy distributions, the pair potential function appears to be a promising tool for evaluating the packing of alpha-helices in membrane proteins. The analysis will be refined by energy calculations based on the expectations that the outer boundary of the channel will be formed by hydrophobic residues in contact with the lipids.

Identification and analysis of the acyl carrier protein (ACP) docking site on beta-ketoacyl-ACP synthase III.

The molecular details that govern the specific interactions between acyl carrier protein (ACP) and the enzymes of fatty acid biosynthesis are unknown. We investigated the mechanism of ACP-protein interactions using a computational analysis to dock the NMR structure of ACP with the crystal structure of beta-ketoacyl-ACP synthase III (FabH) and experimentally tested the model by the biochemical analysis of FabH mutants. The activities of the mutants were assessed using both an ACP-dependent and an ACP-independent assay. The ACP interaction surface was defined by mutations that compromised FabH activity in the ACP-dependent assay but had no effect in the ACP-independent assay. ACP docked to a positively charged/hydrophobic patch adjacent to the active site tunnel on FabH, which included a conserved arginine (Arg-249) that was required for ACP docking. Kinetic analysis and direct binding studies between FabH and ACP confirmed the identification of Arg-249 as critical for FabH-ACP interaction. Our experiments reveal the significance of the positively charged/hydrophobic patch located adjacent to the active site cavities of the fatty acid biosynthesis enzymes and the high degree of sequence conservation in helix II of ACP across species.

Developmental switch in excitability, Ca(2+) and K(+) currents of retinal ganglion cells and their dendritic structure.

In the retina of teleost fish, continuous neuronal development occurs at the margin, in the peripheral growth zone (PGZ). We prepared tissue slices from the retina of rainbow trout that include the PGZ and that comprise a time line of retinal development, in which cells at progressive stages of differentiation are present side by side. We studied the changes in dendritic structure and voltage-dependent Ca(2+), Na(+), and K(+) currents that occur as ganglion cells mature. The youngest ganglion cells form a distinct bulge. Cells in the bulge have spare and short dendritic trees. Only half express Ca(2+) currents and then only high-voltage-activated currents with slow inactivation (HVAslow). Bulge cells are rarely electrically excitable. They express a mixture of rapidly inactivating and noninactivating K(+) currents (IKA and IKdr). The ganglion cells next organize into a transition zone, consisting of a layered structure two to three nuclei thick, before forming the single layered structure characteristic of the mature retina. In the transition zone, the dendritic arbor is elaborately branched and extends over multiple laminae in the inner plexiform layer, without apparent stratification. The arbor of the mature cells is stratified, and the span of the dendritic arbor is well over five times the cell body's diameter. The electrical properties of cells in the transition and mature zones differ significantly from those in the bulge cells. Correlated with the more elaborate dendritic structures are the expression of both rapidly inactivating HVA (HVAfast) and of low-voltage-activated (LVA) Ca(2+) currents and of a high density of Na(+) currents that renders the cells electrically excitable. The older ganglion cells also express a slowly activating K(+) current (IKsa).

Structural symmetry and protein function.

The majority of soluble and membrane-bound proteins in modern cells are symmetrical oligomeric complexes with two or more subunits. The evolutionary selection of symmetrical oligomeric complexes is driven by functional, genetic, and physicochemical needs. Large proteins are selected for specific morphological functions, such as formation of rings, containers, and filaments, and for cooperative functions, such as allosteric regulation and multivalent binding. Large proteins are also more stable against denaturation and have a reduced surface area exposed to solvent when compared with many individual, smaller proteins. Large proteins are constructed as oligomers for reasons of error control in synthesis, coding efficiency, and regulation of assembly. Symmetrical oligomers are favored because of stability and finite control of assembly. Several functions limit symmetry, such as interaction with DNA or membranes, and directional motion. Symmetry is broken or modified in many forms: quasisymmetry, in which identical subunits adopt similar but different conformations; pleomorphism, in which identical subunits form different complexes; pseudosymmetry, in which different molecules form approximately symmetrical complexes; and symmetry mismatch, in which oligomers of different symmetries interact along their respective symmetry axes. Asymmetry is also observed at several levels. Nearly all complexes show local asymmetry at the level of side chain conformation. Several complexes have reciprocating mechanisms in which the complex is asymmetric, but, over time, all subunits cycle through the same set of conformations. Global asymmetry is only rarely observed. Evolution of oligomeric complexes may favor the formation of dimers over complexes with higher cyclic symmetry, through a mechanism of prepositioned pairs of interacting residues. However, examples have been found for all of the crystallographic point groups, demonstrating that functional need can drive the evolution of any symmetry.

Ionization state and molecular docking studies for the macrophage migration inhibitory factor: the role of lysine 32 in the catalytic mechanism.

The macrophage migration inhibitory factor (MIF) is a cytokine that is structurally similar to certain isomerases and for which multiple immune and catalytic roles have been proposed. Different catalytic activities have been reported for MIF, yet the exact mechanism by which MIF acts is not completely known. As a tautomerase, the enzyme uses a general acid-base mechanism of proton transfer in which the amino-terminal proline has been shown to function as the catalytic base. We report the results of molecular docking simulations of macrophage migration inhibitory factor with three substrates, D-dopachrome, L-dopachrome methyl ester and p-(hydroxyphenyl)pyruvate. Electrostatic pK(a) predictions were also performed for the free and complexed forms of the enzyme. The predicted binding mode of p-(hydroxyphenyl)pyruvate is in agreement with the recently published X-ray structure. A model for the binding mode of D-dopachrome and L-dopachrome methyl ester to MIF is proposed which offers insights into the catalytic mechanism of D-dopachrome tautomerase activity of MIF. The proposed catalytic mechanism is further supported by the pK(a) predictions, which suggest that residue Lys32 acts as the general acid for the enzymatic catalysis of D-dopachrome.

Alteration of substrate and inhibitor specificity of feline immunodeficiency virus protease.

Feline immunodeficiency virus (FIV) protease is structurally very similar to human immunodeficiency virus (HIV) protease but exhibits distinct substrate and inhibitor specificities. We performed mutagenesis of subsite residues of FIV protease in order to define interactions that dictate this specificity. The I37V, N55M, M56I, V59I, and Q99V mutants yielded full activity. The I37V, N55M, V59I, and Q99V mutants showed a significant increase in activity against the HIV-1 reverse transcriptase/integrase and P2/nucleocapsid junction peptides compared with wild-type (wt) FIV protease. The I37V, V59I, and Q99V mutants also showed an increase in activity against two rapidly cleaved peptides selected by cleavage of a phage display library with HIV-1 protease. Mutations at Q54K, I98P, and L101I dramatically reduced activity. Mutants containing a I35D or I57G substitution showed no activity against either FIV or HIV substrates. FIV proteases all failed to cut HIV-1 matrix/capsid, P1/P6, P6/protease, and protease/reverse transcriptase junctions, indicating that none of the substitutions were sufficient to change the specificity completely. The I37V, N55M, M56I, V59I, and Q99V mutants, compared with wt FIV protease, all showed inhibitor specificity more similar to that of HIV-1 protease. The data also suggest that FIV protease prefers a hydrophobic P2/P2' residue like Val over Asn or Glu, which are utilized by HIV-1 protease, and that S2/S2' might play a critical role in distinguishing FIV and HIV-1 protease by specificity. The findings extend our observations regarding the interactions involved in substrate binding and aid in the development of broad-based inhibitors.

Structural studies of FIV and HIV-1 proteases complexed with an efficient inhibitor of FIV protease.

Three forms of feline immunodeficiency virus protease (FIV PR), the wild type (wt) and two single point mutants, V59I and Q99V, as well as human immunodeficiency virus type 1 protease (HIV-1 PR), were cocrystallized with the C2-symmetric inhibitor, TL-3. The mutants of FIV PR were designed to replace residues involved in enzyme-ligand interactions by the corresponding HIV-1 PR residues at the structurally equivalent position. TL-3 shows decreased (improved) inhibition constants with these FIV PR mutants relative to wt FIV PR. Despite similar modes of binding of the inhibitor to all PRs (from P3 to P3'), small differences are evident in the conformation of the Phe side chains of TL-3 at the P1 and P1' positions in the complexes with the mutated FIV PRs. The differences mimick the observed binding of TL-3 in HIV-1 PR and correlate with a significant improvement in the inhibition constants of TL-3 with the two mutant FIV PRs. Large differences between the HIV-1 and FIV PR complexes are evident in the binding modes of the carboxybenzyl groups of TL-3 at P4 and P4'. In HIV-1 PR:TL-3, these groups bind over the flap region, whereas in the FIV PR complexes, the rings are located along the major axis of the active site. A significant difference in the location of the flaps in this region of the HIV-1 and FIV PRs correlates with the observed conformational changes in the binding mode of the peptidomimetic inhibitor at the P4 and P4' positions. These findings provide a structural explanation of the observed Ki values for TL-3 with the different PRs and will further assist in the development of improved inhibitors.

Connexin 26 expression and extensive gap junctional coupling in cultures of GT1-7 cells secreting gonadotropin-releasing hormone.

Gap junctions (GJs) are transmembrane channels that permit rapid intercellular transit of various small molecules including ions, second messengers and metabolites. GJs promote communication and coordinated activity between coupled neurons, and may help facilitate the synchronous release and pulsatile secretion of neurohormones. A previous study using GnRH-secreting GT1-7 cells reported that connexin 26 was the major GJ subunit present, and that about 20% of the cultured cells engaged in GJ coupling as assayed by fluorescence recovery after photobleaching of 5,6-carboxyfluorescein diacetate (MW 460 D). To reassess GJ connectivity with a more permeant probe, we grew GT1-7 cells to 70% confluency on Matrigel-coated glass coverslips and microinjected Neurobiotin(TM) (MW 322 D) into single cells. Dye was allowed to diffuse for 30 min before cultures were fixed, and subsequently immunostained for Neurobiotin with 3,3'-diaminobenzidine HCl and examined by light microscopy. Dye coupling between 2 or more GT1-7 cells was observed after 75% of all microinjections. Connectivity involved the somata and neurites of an average of 6.6 +/- 2.0 adjoining cells, but in one instance was seen in a group of 32 GT1-7 neighbors. Western blotting and immunofluorescence staining confirmed that connexin 26 was the predominant GJ subunit expressed by GT1-7 cultures. Our results using Neurobiotin suggest these GJ channels may be smaller than anticipated. In addition, functional GJ connectivity between subconfluent GT1-7 cells is more extensive than previously reported, occurring with higher frequency and coupling significantly greater numbers of cultured cells. Since cAMP, IP3, and Ca(2+) are able to pass through GJs and can elicit secretion of GnRH by GT1 cell cultures, GJs may play an important role in the coordination and synchronization of GnRH release.

Revisiting catalysis by chymotrypsin family serine proteases using peptide substrates and inhibitors with unnatural main chains.

Chymotrypsin family serine proteases play essential roles in key biological and pathological processes and are frequently targets of drug discovery efforts. This large enzyme family is also among the most advanced model systems for detailed studies of enzyme mechanism and structure/function relationships. Productive interactions between these enzymes and their substrates are widely believed to mimic the "canonical" interactions between serine proteases and "standard" inhibitors observed in numerous protease-inhibitor complexes. To test this central hypothesis we have synthesized and characterized a series of peptide analogs, based on model substrates and inhibitors of trypsin, that contain unnatural main chains. These results call into question a long accepted theory regarding the interaction of chymotrypsin family serine proteases with substrates and suggest that the canonical interactions observed between these enzymes and standard inhibitors may represent nonproductive rather than productive, substrate-like interactions.

Docking of 4-oxalocrotonate tautomerase substrates: implications for the catalytic mechanism.

The enzyme 4-oxalocrotonate tautomerase catalyzes the ketonization of dienols, which after further processing become intermediates in the Krebs cycle. The enzyme uses a general acid-base mechanism for proton transfer: the amino-terminal proline has been shown to function as the catalytic base and Arg39 has been implicated as the catalytic acid. We report the results of molecular docking simulations of 4-oxalocrotonate tautomerase with two substrates, 2-hydroxymuconate and 5-carboxymethyl-2-hydroxymuconate. pKa calculations are also performed for the free enzyme. The predicted binding mode of 2-hydroxymuconate is in agreement with experimental data. A model for the binding mode of 5-carboxymethyl-2-hydroxymuconate is proposed which explains the lower catalytic efficiency of the enzyme toward this substrate. The pKa predictions and docking simulations support residue Arg39 as the general acid for the enzyme catalysis.

Integrating computation and visualization for biomolecular analysis: an example using python and AVS.

One of the challenges in biocomputing is to enable the efficient use of a wide variety of fast-evolving computational methods to simulate, analyze, and understand the complex properties and interactions of molecular systems. Our laboratory investigates several areas including molecular visualization, protein-ligand docking, protein-protein docking, molecular surfaces, and the derivation of phenomenological potentials. In this paper we present an approach based on the Python programming language to achieve a high level of integration between these different computational methods and our primary visualization system AVS. This approach removes many limitations of AVS while increasing dramatically the inter-operability of our computational tools. Several examples are shown to illustrate how this approach enables a high level of integration and inter-operability between different tools, while retaining modularity and avoiding the creation of a large monolithic package that is difficult to extend and maintain.

Coevolution and subsite decomposition for the design of resistance-evading HIV-1 protease inhibitors.

Drug resistance sharply limits the effectiveness of human immunodeficiency virus (HIV) protease inhibitors in acquired immunodeficiency syndrome therapy. In previous work, we presented methods for design of resistance-evading inhibitors using a computational coevolution technique. Here, we report subsite decomposition experiments that examine the relative importance and roles of each subsite in HIV protease, and the constraints on robust inhibitor design that are imposed by possible resistance mutations in each subsite. The results identify several structural features of robust resistance-evading inhibitors for use in drug design, and show their basis in the constraints imposed by the range of allowable mutation in the protease. In particular, the results identify the P3 and P3' sites as being particularly sensitive to protease mutation: inhibitors designed to fill the S3 and S3' sites of the wild-type protease will be susceptible to viral resistance, but inhibitors with side-chains smaller than a phenylalanine residue at P3 and P3', preferably medium-sized amino acids in the range from valine to leucine and isoleucine residues, will be more robust in the face of protease resistance mutation.

Coevolutionary analysis of resistance-evading peptidomimetic inhibitors of HIV-1 protease.

We have developed a coevolutionary method for the computational design of HIV-1 protease inhibitors selected for their ability to retain efficacy in the face of protease mutation. For HIV-1 protease, typical drug design techniques are shown to be ineffective for the design of resistance-evading inhibitors: An inhibitor that is a direct analogue of one of the natural substrates will be susceptible to resistance mutation, as will inhibitors designed to fill the active site of the wild-type or a mutant enzyme. Two design principles are demonstrated: (i) For enzymes with broad substrate specificity, such as HIV-1 protease, resistance-evading inhibitors are best designed against the immutable properties of the active site-the properties that must be conserved in any mutant protease to retain the ability to bind and cleave all of the native substrates. (ii) Robust resistance-evading inhibitors can be designed by optimizing activity simultaneously against a large set of mutant enzymes, incorporating as much of the mutational space as possible.

Modelling of factor Xa-inhibitor complexes: a computational flexible docking approach.

In order to understand the structural basis of Factor Xa (FXa) specificity, structural complexes of FXa with its synthetic inhibitors are determined using a computational docking approach. The AutoDock suite of programs is used to determine the binding modes of the synthetic inhibitors such as 3- and 4-amidinobenzylphenyl ether (ABP), amidinophenyl pyruvic acid (APPA), diamidinobenzofuranyl ethene (DABE), and DX-9065a 2-(5'-amidino-2'-benzofuranyl)-3-(7'amidino-2'-napthyl)-propionic acid (ABAP) to FXa. The synthetic inhibitors docked in the present study are different in size, nature of linkage, and properties. Two sets of simulations were carried out for synthetic inhibitors docking to FXa. In the first set of simulations, no explicit water molecules were included. In the second set of simulations two explicit solvent molecules were considered. In all the computationally predicted synthetic inhibitor complexes of FXa, the specificity pocket residue Asp-189 is involved in hydrogen bonding with the bound inhibitor. The active site water molecule WAT522 is involved in hydrogen bonding with all the bound inhibitors. The computed energies clearly discriminate the high affinity from low affinity binders.