A view of the three dimensional structure of globular proteins

A view of the three dimensional structure of globular proteins based on continuous networks of hydrogen bonds is proposed. Active sites of enzymes and ion sites are prominent and, within the networks, there are islands of hydrophobic regions giving an overall piebald effect to the appearance of the molecule. This point of view was originally suggested by the results of quantum mechanical computations on the coupling between hydrogen bonds. A formalism for the total energy of a globular protein in water is also suggested.The study of five lines of experimental evidence supports this suggestion. The analysis of the experimental X-ray data for ten globular proteins, using the NETWORK program, revealed the existence of these hydrogen bond networks; X-ray data showed that water molecules tend to occupy fixed positions relative to the protein molecule; a survey has shown that water molecules tend to occupy specific positions relative to the hydrogen bonding side chains; experimental evidence on the bulk properties of lysozyme showed that there exist tightly bound water molecules; graphics studies of the ribonucleaseA molecule demonstrated the networks and the piebald effect. This point of view is pictorially simple and, to illustrate the use of such networks, we discuss the simple ion pairs which occur as substructures within the networks.     

The Ribbon of Hydrogen Bonds in the Three-Dimensional Structure of Globular Proteins

We report quantum mechanical computations and experimental evidence which suggest that the backbone conformation of globular proteins depends generally on the conservation of that part of the hydrogen bond network or ribbon which is joined, in general, directly to the backbone and is largely independent of the remainder of this whole network of hydrogen bonds. The familiar hydrogen bonds of the helix and the sheet form about one-half of this ribbon of hydrogen bonds. Both water molecules and hydrogen bonding side chain groups are involved in the formation of the ribbon.This view of the three-dimensional structure of globular proteins in terms of the `molecule' allows us to deal with the non-secondary structure as well as with the familiar secondary structure. It also suggests that the ribbon contains approximately the same number of hydrogen bonds within all three structures – the helix, the sheet and the coil – and that this is the reason for the ease of interconversion of these three structures.The quantum mechanical computations on hydrogen bonding suggest that delocalised water molecules which have substantial mobility are an essential part of the ribbon. This situation arises because the hydrogen bonding groups of the protein molecule are not free to move to optimise the hydrogen bonding geometries as are the oxygen atoms in the waters and ices. Such delocalised water molecules either have high B values or are invisible in the X-ray data and yet are able to form a structure which is as strong as a normal hydrogen bond.The experimental data on the point mutations of the THRI57 residue of the T4 phage lysome provides an initial test of this model. Both the local backbone conformation and the ribbon of hydrogen bonds are conserved throughout all the mutations of residue 157,providing that the delocalised water molecules are accepted as a genuine part of the structure. These mutations include the introduction of hydrocarbon side chains at position 157 when water molecules or other side chain groups take over the formation of the hydrogen bonds.We suggest that, provided steric effects are not important, many point mutations succeed because they leave the ribbon of hydrogen bonds (and so the backbone conformation) largely unchanged.     

A novel and simple interpretation of the three-dimensional structure of globular proteins based on quantum mechanical computations on small model molecules. II. The clusters of myoglobin

We show that the bridges and clusters of the kind found earlier in rubredoxin also occur in myoglobin. The familiar acid–base(salt) bridges are at the heart of many of the clusters. The data reported by Takano [(1977) J. Mol. Biol., 110 , 537–508] on sperm whale metmyoglobin reveals ten clusters together with the heme structure and some unidentified fragments. The data reported by Hanson and Schoenborn [(1981) J. Mol. Biol. 153 , 117–146] on carbonmonoxymyoglobin shows numerous bridges that generally fit into the cluster picture derived from Takano's data. The data reported by S. E. V. Phillips (personal communication) on oxymyoglobin also fits into the same pattern of cluster. All but one of the myoglobin clusters are formed from two pieces of local structure each of which covers up to fifteen consecutive aminoacid residues. A “long-range” union joins the two pieces of local structure. The clusters are generally involved with the nonhelical structure and with the junctions of the two helices. Two applications of the theory are suggested. First, the clusters offer an explanation for the formation of the α-helices and second the unfolding of the myoglobin molecule at pH 3.8 to 5.3 is readily interpreted as the protonation at this pH of the carboxylate ions that generate the clusters that in turn support the α-helices.     

A theoretical study of cooperative and anticooperative effects on hydrogen-bonded clusters of water and the cyanuric acid

Ab initio and density functional calculations are used to analyse the interaction between a molecule of the cyanuric acid and one, two and three molecules of water at B3LYP/6-311++ G(d,p) and MP2/6-311++ G(d,p) computational levels. Also, the cooperative effect (CE) in terms of the stabilisation energy of clusters is calculated and discussed. Depending on the geometry of clusters under study, the cooperative, non- or anti-CE was found with an increasing cluster size. Red shifts of N–H and C = O stretching frequencies illustrate a good dependence on the CE. The atoms in molecules theory is used to analyse the CE on topological parameters.     

Hydration of Nucleic Acid Bases Studied Using Novel Atom-Atom Potential Functions

Abstract

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contain terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases—uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecule interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interaction play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by “bridges” formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

Aggregation of Kanamycin A: dimer formation with physiological cations

Global cluster geometry optimization has focused so far on clusters of atoms or of compact molecules. We are demonstrating here that present-day techniques also allow to globally optimize clusters of extended, flexible molecules, and that such studies have immediate relevance to experiment. For example, recent experimental findings point to production of larger clusters of an aminoglycoside closely related to Kanamycin A (KA), together with certain preferred physiological cations, by Pseudomonas aeruginosa. The present study provides first theoretical support for KA clustering, with a close examination of the monomer, the bare dimer, and dimers with sodium and potassium cations, employing global cluster structure optimization, in conjunction with force fields, semiempirical methods, DFT and ab-initio approaches. Interestingly, already at this stage the theoretical findings support the experimental observation that sodium cations are preferred over potassium cations in KA clusters, due to fundamentally different cationic embedding. Theoretically predicted NMR and IR spectra for these species indicate that it should be possible to experimentally detect the aggregation state and even the cationic embedding mode in such clusters.     

Hydration of oligonucleotides in crystals

The solvent molecules found around crystallized oligonucleotides after X-ray refinement are analysed in terms of interaction sites to bases, phosphates and sugars in the three main forms of nucleic acid structures, the A-form, the B-form and the Z-form. The average numbers of contacts to nucleic acid atoms made by solvent molecules are identical in the three forms, but it appears that the average number of contacts solvent molecules make with each other depends on the resolution of the structure. The phosphate anionic oxygen atoms are the most hydrated, while the O(3′) and O(5′) backbone atoms and the ring oxygen atom O(4′) are the least hydrated. Among the hydrophilic atoms of the bases, there is a modulation of the relative water affinities with the nucleic acid form. Numerous hydration sites are such that water molecules can bridge hydrophilic atoms of the same residue, of adjacent residues on the same strand, of distant residues on the two strands, or belonging to symmetry-related residues. Through the helical periodicity of the nucleic acid structure, those bridges can lead to regular and striking hydration networks involving several water molecules and characteristic of the nucleic acid form. Solvent dynamics, as seen by temperature factor versus occupancy plots, seems intimately related to nucleic acid structure and dynamics, since they depend on hydration sites around the nucleic acids.     

Correlations between internal mobility and stability of globular proteins.

The recent work is surveyed which leads to the suggestions that the conformation of globular proteins in solution corresponds to a dynamic ensemble of rapidly interconverting spatial structures, that clusters of hydrophobic amino acid side chains have an important role in the architecture of protein molecules, and that mechanistic aspects of protein denaturation can be correlated with internal mobility seen in the native conformation. These conclusions resulted originally from high resolution 1H nuclear magnetic resonance (NMR) studies of aromatic ring mobility, exchange of interior amide protons and thermal denaturation of the basic pancreatic trypsin inhibitor and a group of related proteins. Various new approaches to further characterize proteins in solution have now been taken and preliminary data are presented. These include computer graphics to outline hydrophobic clusters in globular protein structures, high resolution 1H-NMR experiments at variable hydrostatic pressure and 13C-NMR relaxation measurements. At the present early stage of these new investigations it appears that the hydrophobic cluster model for globular proteins is compatible with the data obtained.     

Hydration of nucleic acid bases studied using novel atom-atom potential functions

New simple atom-atom potential functions for simulating behavior of nucleic acids and their fragments in aqueous solutions are suggested. These functions contains terms which are inversely proportional to the first (electrostatics), sixth (or tenth for the atoms, forming hydrogen bonds) and twelfth (repulsion of all the atoms) powers of interatomic distance. For the refinement of the potential function parameters calculations of ice lattice energy, potential energy and configuration of small clusters consisting of water and nucleic acid base molecules as well as Monte Carlo simulation of liquid water were performed. Calculations using new potential functions give rise to more linear hydrogen bonds between water and base molecules than using other potentials. Sites of preferential hydration of five nucleic bases - uracil, thymine, cytosine, guanine and adenine as well as of 6,6,9-trimethyladenine were found. In the most energetically favourable sites water molecular interacts with two adjacent hydrophilic centres of the base. Studies of interaction of the bases with several water molecules showed that water-water interactions play an important role in the arrangement of the nearest to the base water molecules. Hydrophilic centres are connected by "bridges" formed by hydrogen bonded water molecules. The results obtained are consistent with crystallographic and mass-spectrometric data.     

Visualisation of extensive water ribbons and networks in a DNA minor-groove drug complex.

The crystal structure is reported of a complex between an ethyl derivative of the minor-groove drug furamidine and the dodecanucleotide duplex d(CGCGAATTCGCG)2, which has been refined to 1.85 A resolution and an R factor of 16.6% for data collected at -173 degreesC. An exceptionally large number (220) of water molecules have been located. The majority of these occur in the first coordination shell of solvation. There are extensive networks of connected waters, both in the major and minor grooves. In particular, there are 21 water molecules associated with the minor-groove drug, via hydrogen bonds from the four charged nitrogen atoms. One cluster of four waters is situated in the groove itself; the majority are on the outer edge of the groove, and serve to bridge between the outward-directed drug nitrogen atoms and backbone phosphate oxygen atoms. These bridges are both intra- and inter-strand, with the net effect that the outer edge of the drug molecule is covered by ribbons of water molecules.     

Iron (III) can be transferred between ferritin molecules.

The iron-storage molecule ferritin can sequester up to 4500 Fe atoms as the mineral ferrihydrite. The iron-core is gradually built up when FeII is added to apoferritin and allowed to oxidize. Here we present evidence, from M?ssbauer spectroscopic measurements, for the surprising result that iron atoms that are not incorporated into mature ferrihydrite particles, can be transferred between molecules. Experiments were done with both horse spleen ferritin and recombinant human ferritin. M?ssbauer spectroscopy responds only to 57Fe and not to 56Fe and can distinguish chemically different species of iron. In our experiments a small number of 57FeII atoms were added to two equivalent apoferritin solutions and allowed to oxidize (1-5 min or 6 h). Either ferritin containing a small iron-core composed of 56Fe, or an equal volume of NaCl solution, was added and the mixture frozen in liquid nitrogen to stop the reaction at a chosen time. Spectra of the ferritin solution to which only NaCl was added showed a mixture of species including 57FeIII in solitary and dinuclear sites. In the samples to which 150 56FeIII-ferritin had been added the spectra showed that all, or almost all, of the 57FeIII was in large clusters. In these solutions 57FeIII initially present as intermediate species must have migrated to molecules containing large clusters. Such migration must now be taken into account in any model of ferritin iron-core formation.     

Learning about protein folding via potential functions

Over the last few years we have developed an empirical potential function that solves the protein structure recognition problem: given the sequence for an n-residue globular protein and a collection of plausible protein conformations, including the native conformation for that sequence, identify the correct, native conformation. Having determined this potential on the basis of only some 6500 native/nonnative pairs of structures for 58 proteins, we find it recognizes the native conformation for essentially all compact, soluble, globular proteins having known native conformations in comparisons with 10⁴ to 10⁶ reasonable alternative conformations apiece. In this sense, the potential encodes nearly all the essential features of globular protein conformational preference. In addition it “knows” about many additional factors in protein folding, such as the stabilization of multimeric proteins, quaternary structure, the role of disulfide bridges and ligands, proproteins vs. processed proteins, and minimal strand lengths in globular proteins. Comparisons are made with other sorts of protein folding problems, and applications in protein conformational determination and prediction are discussed. © 1994 Wiley-Liss, Inc.     

Hydration of transfer RNA molecules: a crystallographic study

Four crystal structures of transfer RNA molecules were refined at 3 A resolution with the inclusion of the solvent molecules found in the difference maps: yeast tRNA-phe in the orthorhombic form, yeast tRNA-phe in the monoclinic form and yeast tRNA-asp in the A and B forms. Over 100 solvent molecules were located in each tRNA crystal. Several hydration schemes are found repeatedly in the 4 crystals. The tertiary interactions in the corner of the L-shaped molecule attract numerous solvent molecules which bridge the ribose hydroxyl O(2') atoms, base exocyclic atoms and phosphate anionic oxygen atoms. Conservation of bases leads to conservative localized hydration patterns. Several solvent molecules are found stabilizing unusual base pairs like the G-U pairs and those involving the pseudouridine base. Water bridges between the O(2') and the exocyclic atom O2 of pyrimidines or the N3 atom of purines are common. Water bridges occur frequently between successive anionic oxygen atoms of each strand as well as between N7 or other exocyclic atoms of successive bases in the major groove. Magnesium ions or spermine molecules are found to bind in the major groove of tRNA helices without specific interactions.     

The molecular stoichiometric hydration model (SHM) as applied to tendon/collagen, globular proteins and cells

This report describes and documents the presence of multiple water-of-hydration fractions on proteins and in cells. Initial studies of hydration fractions in g of water/g of DM (dry mass) for tendon/collagen led to the development of the molecular SHM (stoichiometric hydration model) and the development of methods for calculating the size of hydration fractions on a number of different proteins of known amino acid composition. The water fractions have differences in molecular motion and other physical properties due to electrostatic interactions of polar water molecules with electric fields generated by covalently bound pairs of opposite partial charge on the protein backbone. The methods allow calculation of the size of four hydration fractions: single water bridges, double water bridges, dielectric water clusters over polar-hydrophilic surfaces and water clusters over hydrophobic surfaces. These four fractions provide monolayer water coverage. The predicted SHM hydration fractions match closely measured hydration fraction values for collagen and for globular proteins. This report also presents water sorption findings that support the SHM. The SHM is applicable for cell systems where it has been studied. In seven cell systems studied, more than half of all of the cell water had properties unlike those of bulk water. The SHM predicts and explains the commonly cited and measured bound water fraction of 0.2-0.4 g of water/g of DM on proteins. The commonly accepted concept that water beyond this bound water fraction can be considered bulk-like water in its physical properties is unwarranted.     

Germ-cell cluster formation in the telotrophic meroistic ovary of <Emphasis Type="Italic">Tribolium castaneum</Emphasis> (Coleoptera,Polyphaga, Tenebrionidae) and its implication on insect phylogeny

Tribolium castaneum has telotrophic meroistic ovarioles of the Polyphaga type. During larval stages, germ cells multiply in a first mitotic cycle forming many small, irregularly branched germ-cell clusters which colonize between the anterior and posterior somatic tissues in each ovariole. Because germ-cell multiplication is accompanied by cluster splitting, we assume a very low number of germ cells per ovariole at the beginning of ovariole development. In the late larval and early pupal stages, we found programmed cell death of germ-cell clusters that are located in anterior and middle regions of the ovarioles. Only those clusters survive that rest on posterior somatic tissue. The germ cells that are in direct contact with posterior somatic cells transform into morphologically distinct pro-oocytes. Intercellular bridges interconnecting pro-oocytes are located posteriorly and are filled with fusomes that regularly fuse to form polyfusomes. Intercellular bridges connecting pro-oocytes to pro-nurse cells are always positioned anteriorly and contain small fusomal plugs. During pupal stages, a second wave of metasynchronous mitoses is initiated by the pro-oocytes, leading to linear subclusters with few bifurcations. We assume that the pro-oocytes together with posterior somatic cells build the center of determination and differentiation of germ cells throughout the larval, pupal, and adult stages. The early developmental pattern of germ-cell multiplication is highly similar to the events known from the telotrophic ovary of the Sialis type. We conclude that among the common ancestors of Neuropterida and Coleoptera, a telotrophic meroistic ovary of the Sialis type evolved, which still exists in Sialidae, Raphidioptera, and a myxophagan Coleoptera family, the Hydroscaphidae. Consequently, the telotrophic ovary of the Polyphaga type evolved from the Sialis type. Electronic supplementary material Supplementary material is available in the online version of this article at and is accessible for authorized users.     

Molecular dynamics simulation of solvated azurin: correlation between surface solvent accessibility and water residence times

A system containing the globular protein azurin and 3,658 water molecules has been simulated to investigate the influence on water dynamics exerted by a protein surface. Evaluation of water mean residence time for elements having different secondary structure did not show any correlation. Identically, comparison of solvent residence time for atoms having different charge and polarity did not show any clear trend. The main factor influencing water residence time in proximity to a specific site was found to be its solvent accessibility. In detail for atoms belonging to lateral chains and having solvent-accessible surface lower than approximately 16 A(2)a relation is found for which charged and polar atoms are surrounded by water molecules characterized by residence times longer than the non polar ones. The involvement of the low accessible protein atom in an intraprotein hydrogen bond further modulates the length of the water residence time. On the other hand for surfaces having high solvent accessibility, all atoms, independently of their character, are surrounded by water molecules which rapidly exchange with the bulk solvent. Proteins 2000;39:56-67.     

Determination of the accessible surface of globular proteins by means of tritium planigraphy

Results are presented for proteins with known three-dimensional structure (lysozyme, myoglobin, ribonuclease), which show that the probability of label incorporation upon bombardment by hot tritium atoms may be quantitatively linked with the surface area of the protein accessible to water molecules. Possible deviations from simple linear dependency caused by particular mechanisms of label introduction are discussed. The data obtained in experiments with model systems were used to determine the accessible surface area of human serum albumin, for which structural data is not sufficiently accurate to allow estimation of accessible surface area. Experimental data correlate reasonably well with estimations based on conventional concepts of the relationship between accessible surface area and molecular weight for globular proteins. Correspondence to: A. V. Volynskaya     

Formation of germ-line cysts with a central cytoplasmic core is accompanied by specific orientation of mitotic spindles and partitioning of existing intercellular bridges

Animal germ cells tend to form clonal groups known as clusters or cysts. Germ cells within the cyst (cystocytes) are interconnected by intercellular bridges and thus constitute a syncytium. Our knowledge of the mechanisms that control the formation of germ-cell clusters comes from extensive studies carried on model organisms (Drosophila, Xenopus). Germ-cell clusters have also been described in worms (annelids, flat worms and nematodes), although their architecture differs significantly from that known in arthropods or vertebrates. Their peculiar feature is the presence of a central anucleate cytoplasmic core (cytophore, rachis) around which the cystocytes are clustered. Each cystocyte in such a cluster always has one intercellular bridge connecting it to the central cytoplasmic core. The way that such clusters are formed has remained a riddle for decades. By means of light, fluorescence and electron microscopy, we have analysed the formation and architecture of cystocyte clusters during early stages of spermatogenesis and oogenesis in a few species belonging to clitellate (oligochaetous) annelids. Our data indicate that the appearance of germ cells connected via a central cytophore is accompanied by a specific orientation of the mitotic spindles during cystocyte divisions. Spindle long axes are always oriented tangentially to the surface of the cytophore. In consequence, cystocytes divide perpendicularly to the plane of the existing intercellular bridge. Towards the final stages of cytokinesis, the contractile ring of the cleavage furrow merges with the rim of the intercellular bridge that connects the dividing cystocyte with the cytophore and forces partition of the existing bridge into two new bridges. This work was supported by the following research grants: 2P04C004 28 from the Ministry of Science and Informatization (to P. Świątek and J. Klag) and DS/1018/IZ/2007 (to J. Kubrakiewicz).     

Monte Carlo studies on equilibrium globular protein folding. I. Homopolymeric lattice models of beta-barrel proteins

Dynamic Monte Carlo studies have been performed on various diamond lattice models of β-proteins. Unlike previous work, no bias toward the native state is introduced; instead, the protein is allowed to freely hunt through all of phase space to find the equilibrium conformation. Thus, these systems may aid in the elucidation of the rules governing protein folding from a given primary sequence; in particular, the interplay of short- vs long-range interaction can be explored. Three distinct models (A? C) were examined. In model A, in addition to the preference for trans (t) over gauche states (g⁺ and g^?) (thereby perhaps favoring β-sheet formation), attractive interactions are allowed between all nonbonded, nearest neighbor pairs of segments. If the molecules possess a relatively large fraction of t states in the denatured form, on cooling spontaneous collapse to a well-defined β-barrel is observed. Unfortunately, in model A the denatured state exhibits too much secondary structure to correctly model the globular protein collapse transition. Thus in models B and C, the local stiffness is reduced. In model B, in the absence of long-range interactions, t and g states are equally weighted, and cooperativity is introduced by favoring formation of adjacent pairs of nonbonded (but not necessarily parallel) t states. While the denatured state of these systems behaves like a random coil, their native globular structure is poorly defined. Model C retains the cooperativity of model B but allows for a slight preference of t over g states in the short-range interactions. Here, the denatured state is indistinguishable from a random coil, and the globular state is a well-defined β-barrel. Over a range of chain lengths, the collapse is well represented by an all-or-none model. Hence, model C possesses the essential qualitative features observed in real globular proteins. These studies strongly suggest that the uniqueness of the globular conformation requires some residual secondary structure to be present in the denatured state.