Analyzing protein functions in four dimensions

Time-resolved structural studies on biomolecular function are coming of age. Focus has shifted from studies on 'systems of opportunities' to a more problem-oriented approach, addressing significant questions in biology and chemistry. An important step in this direction has been the use of physical and chemical trapping methods to capture and then freeze reaction intermediates in crystals. Subsequent monochromatic data collection at cryogenic temperatures can produce high resolution structures of otherwise elusive intermediates. The combination of diffraction methods with spectroscopic techniques provides a means to directly correlate electronic transitions with structural transitions in the sample, eliminating much of the guesswork from experiments. Studies on cytochrome P450, isopenicillin N synthase, cytochrome cd1 nitrite reductase, copper amine oxidase and bacteriorhodopsin were selected as examples, and the results are discussed.     

Can local energy minimization refine the PDB structures of different resolution universally?

Local energy minimization was statistically tested as the refinement strategy for PDB structure pairs of different resolution. The 13 pairs of structures with the only difference being the resolution were extracted from PDB and represented structures of 11 identical proteins obtained with different x-ray diffraction techniques. The rmsd distribution was calculated for these pairs before and after local energy minimization of each structure. MMFF94 was used for energy calculations and the quasi-Newtonian method was used for local energy minimization. By comparison of these two rmsd distributions, the local energy minimization was proved to statistically increase the structural differences in pairs, so it cannot be used for refinement purposes. To explore the prospects of complex refinement strategies based on energy minimization, randomized structures were obtained by moving the initial PDB structures as far as the minimized structures had been moved in the multidimensional space of atomic coordinates. For these randomized structures the rmsd distribution was calculated and compared with the one for minimized structures. The significant differences in their mean values proved the energy surface of the protein to have only few minima near the conformations of different resolution obtained by x-ray analysis for PDB. Some other results we obtained exploring the energy surface near these conformations are also presented. These results are expected to be useful for the development of new protein refinement strategies based on energy minimization.     

Diving Into the Lipid Bilayer to Investigate the Transmembrane Organization and Conformational State Transitions of P-type Ion ATPases

Although membrane proteins constitute more than 20% of the total proteins, the structures of only a few are known in detail. An important group of integral membrane proteins are ion-transporting ATPases of the P-type family, which share the formation of an acid-stable phosphorylated intermediate as part of their reaction cycle. There are several crystal structures of the sarcoplasmic reticulum Ca(2+) pump (SERCA) revealing different conformations, and recently, crystal structures of the H(+)-ATPase and the Na(+)/K(+)-ATPase were reported as well. However, there are no atomic resolution structures for other P-type ATPases including the plasma membrane calcium pump (PMCA), which is integral to cellular Ca(2+) signaling. Crystallization of these proteins is challenging because there is often no natural source from which the protein can be obtained in large quantities, and the presence of multiple isoforms in the same tissue further complicates efforts to obtain homogeneous samples suitable for crystallization. Alternative techniques to study structural aspects and conformational transitions in the PMCAs (and other P-type ATPases) have therefore been developed. Specifically, information about the structure and assembly of the transmembrane domain of an integral membrane protein can be obtained from an analysis of the lipid-protein interactions. Here, we review recent efforts using different hydrophobic photo-labeling methods to study the non-covalent interactions between the PMCA and surrounding phospholipids under different experimental conditions, and discuss how the use of these lipid probes can reveal valuable information on the membrane organization and conformational state transitions in the PMCA, Na(+)/K(+)-ATPase, and other P-type ATPases.     

A method for biomolecular structural recognition and docking allowing conformational flexibility.

In this work, we present an algorithm developed to handle biomolecular structural recognition problems, as part of an interdisciplinary research endeavor of the Computer Vision and Molecular Biology fields. A key problem in rational drug design and in biomolecular structural recognition is the generation of binding modes between two molecules, also known as molecular docking. Geometrical fitness is a necessary condition for molecular interaction. Hence, docking a ligand (e.g., a drug molecule or a protein molecule), to a protein receptor (e.g., enzyme), involves recognition of molecular surfaces. Conformational transitions by "hinge-bending" involves rotational movements of relatively rigid parts with respect to each other. The generation of docked binding modes between two associating molecules depends on their three dimensional structures (3-D) and their conformational flexibility. In comparison to the particular case of rigid-body docking, the computational difficulty grows considerably when taking into account the additional degrees of freedom intrinsic to the flexible molecular docking problem. Previous docking techniques have enabled hinge movements only within small ligands. Partial flexibility in the receptor molecule is enabled by a few techniques. Hinge-bending motions of protein receptors domains are not addressed by these methods, although these types of transitions are significant, e.g., in enzymes activity. Our approach allows hinge induced motions to exist in either the receptor or the ligand molecules of diverse sizes. We allow domains/subdomains/group of atoms movements in either of the associating molecules. We achieve this by adapting a technique developed in Computer Vision and Robotics for the efficient recognition of partially occluded articulated objects. These types of objects consist of rigid parts which are connected by rotary joints (hinges). Our method is based on an extension and generalization of the Hough transform and the Geometric Hashing paradigms for rigid object recognition. We show experimental results obtained by the successful application of the algorithm to cases of bound and unbound molecular complexes, yielding fast matching times. While the "correct" molecular conformations of the known complexes are obtained with small RMS distances, additional, predictive good-fitting binding modes are generated as well. We conclude by discussing the algorithm's implications and extensions, as well as its application to investigations of protein structures in Molecular Biology and recognition problems in Computer Vision.     

DynaPred: a structure and sequence based method for the prediction of MHC class I binding peptide sequences and conformations

MOTIVATION: The binding of endogenous antigenic peptides to MHC class I molecules is an important step during the immunologic response of a host against a pathogen. Thus, various sequence- and structure-based prediction methods have been proposed for this purpose. The sequence-based methods are computationally efficient, but are hampered by the need of sufficient experimental data and do not provide a structural interpretation of their results. The structural methods are data-independent, but are quite time-consuming and thus not suited for screening of whole genomes. Here, we present a new method, which performs sequence-based prediction by incorporating information obtained from molecular modeling. This allows us to perform large databases screening and to provide structural information of the results. RESULTS: We developed a SVM-trained, quantitative matrix-based method for the prediction of MHC class I binding peptides, in which the features of the scoring matrix are energy terms retrieved from molecular dynamics simulations. At the same time we used the equilibrated structures obtained from the same simulations in a simple and efficient docking procedure. Our method consists of two steps: First, we predict potential binders from sequence data alone and second, we construct protein-peptide complexes for the predicted binders. So far, we tested our approach on the HLA-A0201 allele. We constructed two prediction models, using local, position-dependent (DynaPred(POS)) and global, position-independent (DynaPred) features. The former model outperformed the two sequence-based methods used in our evaluation; the latter shows a much higher generalizability towards other alleles than the position-dependent models. The constructed peptide structures can be refined within seconds to structures with an average backbone RMSD of 1.53 A from the corresponding experimental structures.     

An efficient automated computer vision based technique for detection of three dimensional structural motifs in proteins.

As the number of available three dimensional coordinates of proteins increases, it is now recognized that proteins from different families and topologies are constructed from independent motifs. Detection of specific structural motifs within proteins aids in understanding their role and the mechanism of their operation. To aid in identification and use of these motifs it has become necessary to develop efficient methods for systematic scanning of structural databases. To date, methods of structural protein comparison suffer from at least one of the following limitations: (1) are not fully automated (require human intervention), (2) are limited to relatively similar structures, (3) are constrained to linear alignments of the structures, (4) are sensitive to insertions, deletions or gaps in the sequences or (5) are very time consuming. We present a method to overcome the above limitations. The method discovers and ranks every piece of structural similarity between the structures compared, thus allowing the simultaneous detection of real 3-D motifs in different domains, between domains, in active sites, surfaces etc. The method uses the Geometric Hashing Paradigm which is an efficient technique originally developed for Computer Vision. The algorithm exploits the geometrical constraints of rigid objects, it is especially geared towards recognition of partial structures in rigid objects belonging to large data bases and is straightforwardly parallelizable. Computer Vision techniques are for the first time applied to molecular structure comparison, resulting in an efficient, fully automated tool. The method has been tested in a number of cases, including comparisons of the haemoglobins, immunoglobulins, serine proteinases, calcium binding proteins, DNA binding proteins and others. In all examples our results were equivalent to the published results from previous methods and in some cases additional structural information was obtained by our method.     

Study of polynucleotide conformation by resolution-enhanced ultraviolet spectroscopy poly(rC) and poly(dC).

Self-deconvolution and the fourth derivative of ultraviolet absorption spectra have been used to study stacked single-stranded and double-helix structures of different cytosine-containing polynucleotides for the first time. These compounds were studied under different solution conditions (pH and organic solvents) and at low temperatures. The red shift of the lower band (B2u band plus possibly some n-->pi* transition) of the absorption spectra in the cytosine-containing polynucleotides and the appearance of new peaks in the deconvoluted and derivative spectra in the 280-310 nm region are attributed mainly to cytosine-cytosine stacking interactions. In particular, the fourth-derivative peaks at wavelengths higher than 290 nm can be associated to coupling of electronic transitions of cytosine bases. The nature of the electronic transitions producing the absorption bands which are resolved in the aforementioned fourth-derivative peaks is discussed. It is concluded that the resolution-enhancement techniques used in this work, i.e. self-deconvolution and fourth derivative, complement each other and are useful methods to study structural changes of single-stranded and double-stranded polynucleotides allowing, at the same time, more information to be obtained about specific stacking interactions than classical absorption spectrophotometry.     

Structural genomics studies of human caries pathogen Streptococcus mutans

Gram-positive bacterium Streptococcus mutans is the primary causative agent of human dental caries. To better understand this pathogen at the atomic structure level and to establish potential drug and vaccine targets, we have carried out structural genomics research since 2005. To achieve the goal, we have developed various in-house automation systems including novel high-throughput crystallization equipment and methods, based on which a large-scale, high-efficiency and low-cost platform has been establish in our laboratory. From a total of 1,963 annotated open reading frames, 1,391 non-membrane targets were selected prioritized by protein sequence similarities to unknown structures, and clustered by restriction sites to allow for cost-effective high-throughput conventional cloning. Selected proteins were over-expressed in different strains of Escherichia coli. Clones expressed soluble proteins were selected, expanded, and expressed proteins were purified and subjected to crystallization trials. Finally, protein crystals were subjected to X-ray analysis and structures were determined by crystallographic methods. Using the previously established procedures, we have so far obtained more than 200 kinds of protein crystals and 100 kinds of crystal structures involved in different biological pathways. In this paper we demonstrate and review a possibility of performing structural genomics studies at moderate laboratory scale. Furthermore, the techniques and methods developed in our study can be widely applied to conventional structural biology research practice.     

Small-scale batch crystallization of proteins revisited: an underutilized way to grow large protein crystals

Growth of high-quality crystals is a major obstacle in many structural investigations. In recent years, the techniques for screening crystals have improved dramatically, whereas the methods for obtaining large crystals have progressed more slowly. This is an important issue since, although many structures can be solved from small crystals with synchrotron radiation, it is far easier to solve and refine structures when strong data is recorded from large crystals. In an effort to improve the size of crystals, a strategy for a small-scale batch method has been developed that in many cases yields far larger crystals than attainable by vapor diffusion.     

Biased coarse-grained molecular dynamics simulation approach for flexible fitting of X-ray structure into cryo electron microscopy maps

Several approaches have been introduced to interpret, in terms of high-resolution structure, low-resolution structural data as obtained from cryo-EM. As conformational changes are often observed in biological molecules, these techniques need to take into account the flexibility of proteins. Flexibility has been described in terms of movement between rigid domains and between rigid secondary structure elements, which present some limitations for studying dynamical properties. Normal mode analysis has also been used, but is limited to medium resolution data. All-atom molecular dynamics fitting techniques are more appropriate to fit structures into higher-resolution data as full protein flexibility is considered, but are cumbersome in terms of computational time. Here, we introduce a coarse-grained approach; a Go-model was used to represent biological molecules, combined with biased molecular dynamics to reproduce accurately conformational transitions. Illustrative examples on simulated data are shown. Accurate fittings can be obtained for resolution ranging from 5 to 20 Å. The approach was also tested on experimental data of Elongation Factor G and Escherichia coli RNA polymerase, where its validity is compared to previous models obtained from different techniques. This comparison demonstrates that quantitative flexible techniques, as opposed to manual docking, need to be considered to interpret low-resolution data.     

Combination of NMR spectroscopy and X-ray crystallography offers unique advantages for elucidation of the structural basis of protein complex assembly

NMR spectroscopy and X-ray crystallography are two premium methods for determining the atomic structures of macro-biomolecular complexes. Each method has unique strengths and weaknesses. While the two techniques are highly complementary, they have generally been used separately to address the structure and functions of biomolecular complexes. In this review, we emphasize that the combination of NMR spectroscopy and X-ray crystallography offers unique power for elucidating the structures of complicated protein assemblies. We demonstrate, using several recent examples from our own laboratory, that the exquisite sensitivity of NMR spectroscopy in detecting the conformational properties of individual atoms in proteins and their complexes, without any prior knowledge of conformation, is highly valuable for obtaining the high quality crystals necessary for structure determination by X-ray crystallography. Thus NMR spectroscopy, in addition to answering many unique structural biology questions that can be addressed specifically by that technique, can be exceedingly powerful in modern structural biology when combined with other techniques including X-ray crystallography and cryo-electron microscopy.     

Protein O-glycosylation analysis

Abstract This review provides an overview on the methods available for analysis of O-glycosylation. Three major themes are addressed: analysis of released O-glycans including different O-glycan liberation, derivatization, and detection methods; analysis of formerly O-glycosylated peptides yielding information on O-glycan attachment sites; analysis of O-glycopeptides, representing by far the most informative but also most challenging approach for O-glycan analysis. Although there are various techniques available for the identification of O-linked oligosaccharides, the focus here is on MS fragmentation techniques such as collision-induced fragmentation, electron capture dissociation, and electron transfer dissociation. Finally, the O-glycan analytical challenges that need to be met will be discussed.     

An algebraic representation of RNA secondary structures.

This paper develops mathematical methods for describing and analyzing RNA secondary structures. It was motivated by the need to develop rigorous yet efficient methods to treat transitions from one secondary structure to another, which we propose here may occur as motions of loops within RNAs having appropriate sequences. In this approach a molecular sequence is described as a vector of the appropriate length. The concept of symmetries between nucleic acid sequences is developed, and the 48 possible different types of symmetries are described. Each secondary structure possible for a particular nucleotide sequence determines a symmetric, signed permutation matrix. The collection of all possible secondary structures is comprised of all matrices of this type whose left multiplication with the sequence vector leaves that vector unchanged. A transition between two secondary structures is given by the product of the two corresponding structure matrices. This formalism provides an efficient method for describing nucleic acid sequences that allows questions relating to secondary structures and transitions to be addressed using the powerful methods of abstract algebra. In particular, it facilitates the determination of possible secondary structures, including those containing pseudoknots. Although this paper concentrates on RNA structure, this formalism also can be applied to DNA.     

An Algebraic Representation of RNA Secondary Structures

Abstract

This paper develops mathematical methods for describing and analyzing RNA secondary structures. It was motivated by the need to develop rigorous yet efficient methods to treat transitions from one secondary structure to another, which we propose here may occur as motions of loops within RNAs having appropriate sequences. In this approach a molecular sequence is described as a vector of the appropriate length. The concept of symmetries between nucleic acid sequences is developed, and the 48 possible different types of symmetries are described. Each secondary structure possible for a particular nucleotide sequence determines a symmetric, signed permutation matrix. The collection of all possible secondary structures is comprised of all matrices of this type whose left multiplication with the sequence vector leaves that vector unchanged. A transition between two secondary structures is given by the product of the two corresponding structure matrices. This formalism provides an efficient method for describing nucleic acid sequences that allows questions relating to secondary structures and transitions to be addressed using the powerful methods of abstract algebra. In particular, it facilitates the determination of possible secondary structures, including those containing pseudoknots. Although this paper concentrates on RNA structure, this formalism also can be applied to DNA     

Hydrogenases: active site puzzles and progress

Recent research on the hydrogenase reactions has sought to probe beyond the information that is provided by X-ray diffraction structures. The major challenge of locating 'transient' hydrogen atoms in species that are potential catalytic intermediates is being addressed, using advanced electron paramagnetic resonance (EPR) techniques and theoretical methods. This article discusses recent progress towards a consensus on the structures of different states of the active site of hydrogenases, the mechanisms of activation and hydrogen cycling.     

Characterization of novel proteins based on known protein structures

The genome sciences face the challenge to characterize structure and function of a vast number of novel genes. Sequence search techniques are used to infer functional and structural information from similarities to experimentally characterized genes or proteins. The persistent goal is to refine these techniques and to develop alternative and complementary methods to increase the range of reliable inference.Here, we focus on the structural and functional assignments that can be inferred from the known three-dimensional structures of proteins. The study uses all structures in the Protein Data Bank that were known by the end of 1997. The protein structures released in 1998 were then characterized in terms of functional and structural similarity to the previously known structures, yielding an estimate of the maximum amount of information on novel protein sequences that can be obtained from inference techniques.The 147 globular proteins corresponding to 196 domains released in 1998 have no clear sequence similarity to previously known structures. However, 75 % of the domains have extensive structure similarity to previously known folds, and most importantly, in two out of three cases similarity in structure coincides with related function. In view of this analysis, full utilization of existing structure data bases would provide information for many new targets even if the relationship is not accessible from sequence information alone. Currently, the most sophisticated techniques detect of the order of one-third of these relationships.     

Role of heme in structural organization of cytochrome c probed by semisynthesis

The heme prosthetic group of cytochrome c is covalently attached to the protein through thioether bonds to two cysteine side chains. The role of covalent heme attachment to cytochrome c is not understood, and most heme proteins bind the prosthetic group by iron ion ligation and tertiary interactions only. A two-armed attachment seems redundant if the role of covalent connection is to limit heme group orientation or to decouple heme affinity from redox potential. These considerations suggested that one role for covalent attachment of the rigid planar heme might be in organizing the cytochrome c protein structure. Indeed, porphyrin cytochrome c (in which the heme iron ion has been removed) is substantially more ordered than apocytochrome c, having characteristics consistent with a molten globule state. To assess the importance of planar rigidity in ordering this protein, semisynthesis was used to substitute porphyrin by two hydrophobic surrogates, one based on biphenyl and the other on phenanthrene, which have different degrees of planarity and rigidity. The expected two-armed covalent attachment of each surrogate was confirmed in the protein products by a variety of methods including mass spectrometry and NMR. Despite being only about half the size of the porphyrin macrocycle, and lacking any possibility for ligation or polar group interactions with the surrounding protein, the two surrogates confer helix contents that are comparable to that of the molten globule formed by porphyrin cytochrome c under similar solution conditions. The pH titrations of the derivatives monitored by circular dichroism exhibit reversible, bell-shaped folding and unfolding transitions, implying that charge group interactions in the protein are involved in stabilizing the helical structures formed. The thermal transitions of the two derivatives at neutral pH are cooperative, with similar midpoints. The similarity of helical content and structural stability in the two derivatives indicates that the increase in conformational freedom by the biphenyl surrogate does not substantially reduce protein structural stability. The similarity of the two derivatives to porphyrin cytochrome c suggests that the common feature among the three covalently attached groups-their hydrophobicity-is by far the dominant factor in organizing stable structures in the protein.     

Quantitative organization of the known protein x-ray structures. I. Methods and short-length-scale results

We address herein the problem of delineating the relationships between the known protein structures. In order to study this problem, methods have been developed to represent arbitrarily sized fragments of biopolymer backbone, and to compare distributions of such fragments. These methods are applied to a classification of 123 structures representing the entire set of known x-ray structures. The resulting data are analyzed (on the four-C alpha length scale) to determine both the large-scale organization of the set of known structures (i.e., the relationships between large groups of structures, each comprised of proteins that are structurally related) and its local structure (i.e., the quantitative degree of similarity between any two specific structures). It is shown that the set of structures forms a continuum of structural types, ranging from all-helical to all-sheet/barrel proteins. It is further demonstrated that the density of protein structures is not uniform across this continuum, but rather that structures cluster in certain regions, separated by regions of lower population. The properties of the various regions of the structural space are determined. The existence is demonstrated of strong quantitative correlations between the contents of different types of four-C alpha fragments within protein structures, which imply significant constraints on the types of architecture that can occur in proteins. Analysis of the distribution of structures demonstrates some hitherto unsuspected similarities and suggests that, in some circumstances, neither structural similarity nor sequence homology may be necessary conditions for evolutionary relationship between proteins. It is also suggested that these unsuspected similarities may imply similar folding mechanisms for structures of apparently different global architecture. Cases are also noted in which apparently similar structures may fold by different mechanisms. The connection between structure and dynamic properties is discussed, and a possible role of dynamics in the evolution of protein structures is suggested. The sensitivity of the methods presented herein to anomalies of structure refinement is demonstrated. It is suggested that the present results provide a framework for analyzing experimental results on structural similarity obtained using vibrational circular dichroism spectra, which are sensitive to local backbone structure.     

Structural heterogeneity of the various forms of apomyoglobin: implications for protein folding.

Temperature-induced denaturation transitions of different structural forms of apomyoglobin were studied monitoring intrinsic tryptophan fluorescence. It was found that the tryptophans are effectively screened from solvent both in native and acid forms throughout most of the temperature range tested. Thus, the tryptophans' surrounding do not show a considerable change in structure where major protein conformational transitions have been found in apomyoglobin using other techniques. At high temperatures and under strong destabilizing conditions, the tryptophans' fluorescence parameters show sigmoidal thermal denaturation. These results, combined with previous studies, show that the structure of this protein is heterogeneous, including native-like (tightly packed) and molten globule-like substructures that exhibit conformation (denaturation) transitions under different conditions of pH and temperature (and denaturants). The results suggest that the folding of this protein proceeds via two "nucleation" events whereby native-like contacts are formed. One of these events, which involves AGH "core" formation, appears to occur very early in the folding process, even before significant hydrophobic collapse in the rest of the protein molecule. From the current studies and other results, a rather detailed picture of the folding of myoglobin is presented, on the level of specific structures and their thermodynamical properties as well as formation kinetics.     

Structural analysis of amyloid beta peptide fragment (25-35) in different microenvironments

Amyloid beta (Abeta) peptides are one of the classes of amphiphilic molecules that on dissolution in aqueous solvents undergo interesting conformational transitions. These conformational changes are known to be associated with their neuronal toxicity. The mechanism of structural transition involved in the monomeric Abeta to toxic assemblage is yet to be understood at the molecular level. Early results indicate that oriented molecular crowding has a profound effect on their assemblage formation. In this work, we have studied how different microenvironments affect the conformational transitions of one of the active amyloid beta-peptide fragments (Abeta(25-35)). Spectroscopic techniques such as CD and Fourier transform infrared spectroscopy were used. It was observed that a stored peptide concentrates on dissolution in methanol adopts a minor alpha-helical conformation along with unordered structures. On changing the methanol concentration in the solvated film form, the conformation switches to the antiparallel beta-sheet structure on the hydrophilic surface, whereas the peptide shows transition from a mixture of helix and unordered structure into predominantly a beta-sheet with minor contribution of helix structure on the hydrophobic surface. Our present investigations indicate that the conformations induced by the different surfaces dictate the gross conformational preference of the peptide concentrate.