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.

Soluble proteins: size, shape and function.

Why are proteins so big? Why do cells build oligomeric proteins? A visual survey of the protein structures available in the Protein Data Bank sheds new light on these questions.

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.

Morphology of protein-protein interfaces.

BACKGROUND: Most soluble proteins are active as low-number oligomers. Statistical surveys of oligomeric proteins have defined the roles of hydrophobicity and complementarity in the stability of protein interfaces, but tend to average structural features over a diverse set of protein-protein interfaces, blurring information on how individual interfaces are stabilized. RESULTS: We report a visual survey of 136 homodimeric proteins from the Brookhaven Protein Data Bank, with images that highlight the major structural features of each protein-protein interaction surface. Nearly all of these proteins have interfaces formed between two globular subunits. Surprisingly, the pattern of hydrophilicity over the surface of these interfaces is quite variable. Approximately one-third of the interfaces show a recognizable hydrophobic core, with a single large, contiguous, hydrophobic patch surrounded by a ring of intersubunit polar interactions. The remaining two-thirds of the proteins show a varied mixture of small hydrophobic patches, polar interactions and water molecules scattered over the entire interfacial area. Ten proteins in the survey have intertwined interfaces formed by extensive interdigitation of the two subunit chains. These interfaces are very hydrophobic and are associated with proteins that require both stability and internal symmetry. CONCLUSIONS: The archetypal protein interface, with a defined hydrophobic core, is present in only a minority of the surveyed homodimeric proteins. Most homodimeric proteins are stabilized by a combination of small hydrophobic patches, polar interactions and a considerable number of bridging water molecules. The presence or absence of a hydrophobic core within these interfaces does not correlate with specific protein functions.

Structure of tomato busy stunt virus IV. The virus particle at 2.9 A resolution.

The structure of tomato bushy stunt virus has been determined crystallographically to 2.9 A resolution. Details are presented of both the molecular structure and the methods by which it has been solved. The icosahedrally symmetric viral shell is composed of 180 protein subunits (Mr 43,000), with three similar but distinct modes of subunit bonding. This capacity for alternative packing is due to localized flexibility in the folded polypeptide (hinges between domains) and to multiple conformations for surface side-chains. The polypeptide backbone has an essentially invariant fold within a compact domain. A mechanism for correct positioning of the different modes of subunit interaction is evident from the structure of the TBSV particle. Thirty-five residues of the polypeptide chain fold in an ordered way on 60 of the 180 subunits, forming an internal framework. Interaction of folded domains with this framework permits accuracy of long-range geometry (correct curvature and closure) to be determined by unambiguous switching between alternative local contact angles. RNA packs tightly into the particle interior. Protein-RNA interactions occur through parts of the subunit that are flexibly linked to the well-ordered domains of the shell. This variable interaction imposes minimum restrictions on the folding of the RNA chain.

Structure, function and properties of antibody binding sites.

Do antibody combining sites possess general properties that enable them to bind different antigens with varying affinities and to bind novel antigens? Here, we address this question by examining the physical and chemical characteristics most favourable for residues involved in antigen accommodation and binding. Amphipathic amino acids could readily tolerate the change of environment from hydrophilic to hydrophobic that occurs upon antibody-antigen complex formation. Residues that are large and can participate in a wide variety of van der Waals' and electrostatic interactions would permit binding to a range of antigens. Amino acids with flexible side-chains could generate a structurally plastic region, i.e. a binding site possessing the ability to mould itself around the antigen to improve complementarity of the interacting surfaces. Hence, antibodies could bind to an array of novel antigens using a limited set of residues interspersed with more unique residues to which greater binding specificity can be attributed. An individual antibody molecule could thus be cross-reactive and have the capacity to bind structurally similar ligands. The accommodation of variations in antigenic structure by modest combining site flexibility could make an important contribution to immune defence by allowing antibody binding to distinct but closely related pathogens. Tyr and Trp most readily fulfil these catholic physicochemical requirements and thus would be expected to be common in combining sites on theoretical grounds. Experimental support for this comes from three sources, (1) the high frequency of participation by these amino acids in the antigen binding observed in six crystallographically determined antibody-antigen complexes, (2) their frequent occurrence in the putative binding regions of antibodies as determined from structural and sequence data and (3) the potential for movement of their side-chains in known antibody binding sites and model systems. The six bound antigens comprise two small different haptens, non-overlapping regions of the same large protein and a 19 amino acid residue peptide. Out of a total of 85 complementarity determining region positions, only 37 locations (plus 3 framework) are directly involved in antigen interaction. Of these, light chain residue 91 is utilized by all the complexes examined, whilst light chain 32, light chain 96 and heavy chain 33 are employed by five out of the six. The binding sites in known antibody-antigen complexes as well as the postulated combining sites in free Fab fragments show similar characteristics with regard to the types of amino acids present. The possible role of other amino acids is also assessed.(ABSTRACT TRUNCATED AT 400 WORDS)

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.

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.

Constructing lattice models of protein chains with side groups.

An algorithm to construct lattice models of polymers with side chains is presented. A search for the global minimum of the error function for a given lattice-to-chain orientation is done by dynamic programming, making the search both fast and complete. Application of the algorithm is illustrated by constructing lattice models for 12 proteins of different sizes and structural types.

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.

A developmental time line in a retinal slice from rainbow trout.

The retina in teleost fish continues to grow throughout much of the life of the animal, in part by the continuing differentiation of new tissue at the retinal margin, an area termed the peripheral growth zone (PGZ) (Lyall, Q J Micros Sci, 1957:98:101-110). We have developed a retinal slice preparation--including the PGZ--from juvenile rainbow trout (Onchorynchus mykiss), a species in which retinal growth is rapid and the PGZ is correspondingly pronounced. The PGZ slice preparation contains a time line of retinal development, with cells at different stages of maturation present side by side. We present evidence that the birth sequence of the various retinal cell types in the PGZ recapitulates the sequence during embryonic development. We also report data on the rate of growth of the PGZ in juvenile trout in vivo. Finally, we have used the PGZ slice preparation to make whole-cell voltage clamp recordings from individual retinal GCs at both early and late stages of maturation. We report that the amplitude of delayed rectifier and A-type potassium currents increases during GC maturation.

Conditioned media from the injured lower vertebrate CNS promote neurite outgrowth from mammalian brain neurons in vitro.

Unlike most pathways of the mature mammalian central nervous system (CNS), the CNS of lower vertebrates can regenerate after jury, a capacity that may be due to the secretion of neurite-promoting factors from the injured CNS. We report that conditioned media (CM) from the injured optic nerve of the mature goldfish promoted marked neurite outgrowth from dissociated embryonic rat cortical and hindbrain neurons in serum-free, neuron-enriched culture. This property was not shared by CM from intact goldfish optic nerve, or from intact or injured optic nerve of mature rats. Neurite-promoting activity was obtained at concentrations as low as 100 ng total protein/ml of CM from injured goldfish optic nerve, and was associated with a distinctive morphology of neurite outgrowth. Due its properties of non-dialyzability, heat lability, and trypsin sensitivity, the neurite-promoting factor(s) appeared to be one or more protein species of MW greater than 12,000. Factors secreted by the regenerating CNS of lower vertebrates can directly promote outgrowth of mammalian CNS neurons.

Interactive modeling of supramolecular assemblies.

The modeling of supramolecular structure presents two major challenges: (1) managing the large amount of sequence, structural and biochemical data, and (2) presenting the data to the user in a flexible and comprehensible manner that addresses these problems. We describe a visualization environment for the creation and analysis of supramolecular models. A set of modular symmetry tools, collectively called SymGen, has been created, providing a flexible platform for the creation of complex assemblies, with interactive control of all symmetry elements and their parameters. A second tool, SymSearch, allows a range of parameters defined within SymGen to be sampled and the resulting conformations to be evaluated. The environment avoids information overload, caused by the large number of atoms in supramolecular complexes, by using a multiresolution spherical harmonic representation that allows the user to display only essential features. Spherical harmonics also enables control of the triangulation level, allowing the user to reduce the complexity of the geometric description to retain interactive speed. The visual fidelity of the surface data is retained by using texture maps that are independent of the resolution of the underlying triangulation. We describe the design and implementation of this environment, and three illustrative examples of its utility.

Probes for double helical DNA sequence information: molecular mechanics study of a proposed model.

The biological activity of many molecules which bind to double helical DNA is related to the sequence specificity of their binding. The development of new substances of this type is challenging; a general solution to the problem does not exist. A new mechanism for small molecule-duplex DNA interaction termed intercalative displacement is proposed. The approach is promising for the design of new substances which will recognize sequence information on DNA. Molecular mechanics calculations in the absence of solvent and counterions predict that the proposed model is structurally and energetically closely related to the well known process of standard intercalation. The implications of these calculations for experimental studies are discussed.

Visualizing the future of molecular graphics.

The field of computer graphics has played an important role in the advancement of structural molecular biology and in the development of structure-based drug design. This article will provide a brief background on the development of this technology, and then focus on the current trends and future directions in molecular graphics and how they will impact the practice of molecular modeling and design. Specific areas that will be covered include: 1) the development of surface and volume based representations of molecular properties and interactions; 2) new approaches to modeling flexible and multi-component structures, and 3) the impact of object-oriented graphics-based programming and the rapidly growing use of network based computing.

Automated docking and the search for HIV protease inhibitors.

This article will discuss the motivations, technologies, and future directions of computational automated docking in the context of the structure-based rational design of HIV-1 protease inhibitors. Docking simulations are widely used for screening of compound libraries to identify new drug leads, employing a simple model for rapid testing of thousands of compounds. Docking simulations are also useful for lead enhancement, using more detailed models to analyze the atomic interactions between inhibitors and target macromolecules. Major advances have been reported in the development of empirical force fields, which now allow assessment of relative binding strength and drug specificity, and extensions of automated docking techniques allow de novo drug design.

Serials standards work: the next frontier.

Serials, one of the more complicated areas of library technical endeavors, has lacked the benefit of standards for a long time. Even now, with standards beginning to be available, the majority of institutions are not working within standard serials formats. A survey to determine the use of serials standards in libraries was conducted in 1988 by the American Library Association, Resources and Technical Services Division, Serials Section, Committee to Study Serials Standards. In the spring of 1988 a survey was sent to a group encompassing the Association of Research Libraries members, CONSER participants, United States Newspaper Program participants, Microform Project libraries, and some vendors and librarians who attended the Committee meetings on a regular basis. The survey questionnaire assessed the current level of serials standards awareness of librarians and vendors. Topics included the type of serials systems used, standards relevant to serials control and union listing and whether or not they are implemented, types and levels of training staff received in the application of standards, benefits of the standards, and areas where standards are most needed.

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.

Real time surface reconstruction for moving molecular fragments.

Recently we introduced the Reduced Surface as an efficient tool to built molecular surfaces. We describe here how this geometric construct can be used to efficiently reconstruct the solvent excluded surface of a protein for which the coordinates of a subset of atoms are changing. We show that, the complexity of that operation is not dependent upon the size of the molecule and is in O[tlog(t)] where t is the maximum of the number of probes and atoms involved in the reconstruction of the surface. The algorithms described here have been implemented and tested on several proteins. The triangulation of the solvent excluded surface of proteins in which a side chain was changing conformation could be updated at rates ranging from 7 to 22 frames per second. We also applied this method to compute the surface area fluctuation of the FIV protease undergoing a constrained molecular dynamics simulation (16 mobile residues). Rate of 6 frames per second were obtained in this case.