Effect of amino acid substitutions at the subunit interface on the stabiity and aggregation properties of a dimeric protein: role of Arg 178 and Arg 218 at the dimer interface of thymidylate synthase

The significance of two interface arginine residues on the structural integrity of an obligatory dimeric enzyme thymidylate synthase (TS) from Lactobacillus casei was investigated by thermal and chemical denaturation. While the R178F mutant showed apparent stability to thermal denaturation by its decreased tendency to aggregate, the Tm of the R218K mutant was lowered by 5 degrees C. Equilibrium denaturation studies in guanidinium chloride (GdmCl) and urea indicate that in both the mutants, replacement of Arg residues results in more labile quaternary and tertiary interactions. Circular dichroism studies in aqueous buffer suggest that the protein interior in R218K may be less well-packed as compared to the wild type protein. The results emphasize that quaternary interactions may influence the stability of the tertiary fold of TS. The amino acid replacements also lead to notable alteration in the ability of the unfolding intermediate of TS to aggregate. The aggregated state of partially unfolded intermediate in the R178F mutant is stable over a narrower range of denaturant concentrations. In contrast, there is an exaggerated tendency on the part of R218K to aggregate in intermediate concentrations of the denaturant. The 3 A crystal structure of the R178F mutant reveals no major structural change as a consequence of amino acid substitution. The results may be rationalized in terms of mutational effects on both the folded and unfolded state of the protein. Site specific amino acid substitutions are useful in identifying specific regions of TS involved in association of non-native protein structures.

1-Anilino-8-naphthalene sulfonate (ANS) binding to proteins revisited: Investigation of dye binding to molten globule states by electrospray ionization mass spectrometry

The binding of 1-anilino-8-naphthalene-sulfonic acid to globular proteins at acidic pH has been investigated by electrospray ionization mass spectrometry (ESIMS). Mass spectra of apomyoglobin recorded in the pH range 2−7 establish that maximal ANS binding is observed at pH 4.0. As many as seven distinct species may be observed in the gas phase which correspond to protein molecules containing one to six molecules of bound ANS. At neutral pH only a single molecule of ANS is bound. In the case of cytochrome c, maximal binding is observed at pH 4.0, with five molecules being bound. Binding is suppressed at neutral pH. In both cases ESIMS demonstrates maximal ANS binding at pH values where the proteins have been reported to exist in molten globule states. ANS binding is not observed for lysozyme, which has a tightly folded structure over the entire pH range. Reduction of disulfide bonds in lysozyme leads to the detection of ANS-bound species at neutral pH. Binding is suppressed at low pH due to complete unfolding of the reduced protein. The results suggest that ESIMS may provide a convenient method of probing the stoichiometry and distribution of dye complexes with molten protein globules

Modelling, Design and Analysis of Low Frequency Platform for Attenuating Micro-Vibration in Spacecraft

One of the most important factors that affect the pointing of precision payloads and devices in space platforms is the vibration generated due to static and dynamic unbalanced forces of rotary equipments placed in the neighborhood of payload. Generally, such disturbances are of low amplitude, less than 1 kHz, and are termed as ‘micro-vibrations’. Due to low damping in the space structure, these vibrations have long decay time and they degrade the performance of payload. This paper addresses the design, modeling and analysis of a low frequency space frame platform for passive and active attenuation of micro-vibrations. This flexible platform has been designed to act as a mount for devices like reaction wheels, and consists of four folded continuous beams arranged in three dimensions. Frequency and response analysis have been carried out by varying the number of folds, and thickness of vertical beam. Results show that lower frequencies can be achieved by increasing the number of folds and by decreasing the thickness of the blade. In addition, active vibration control is studied by incorporating piezoelectric actuators and sensors in the dynamic model. It is shown using simulation that a control strategy using optimal control is effective for vibration suppression under a wide variety of loading conditions.

Nucleotide sequence of the 3terminal region of belladonna mottle virus (renamed Physalis mottle virus) RNA and analysis of the evolutionary relationships among tymoviral coat proteins

The 3' terminal 1255 nt sequence of Physalis mottle virus (PhMV) genomic RNA has been determined from a set of overlapping cDNA clones. The open reading frame (ORF) at the 3' terminus corresponds to the amino acid sequence of the coat protein (CP) determined earlier except for the absence of the dipeptide, Lys-Leu, at position 110-111. In addiition, the sequence upstream of the CP gene contains the message coding for 178 amino acid residues of the C-terminus of the putative replicase protein (RP). The sequence downstream of the CP gene contains an untranslated region whose terminal 80 nucleotides can be folded into a characteristic tRNA-like structure. A phylogenetic tree constructed after aligning separately the sequence of the CP, the replicase protein (RP) and the tRNA-like structure determined in this study with the corresponding sequences of other tymoviruses shows that PhMV wrongly named belladonna mottle virus [BDMV(I)] is a separate tymovirus and not another strain of BDMV(E) as originally envisaged. The phylogenetic tree in all the three cases is identical showing that any subset of genomic sequence of sufficient length can be used for establishing evolutionary relationships among tymoviruses.

Playing Games on Sets and Models

The most prominent objective of the thesis is the development of the generalized descriptive set theory, as we call it. There, we study the space of all functions from a fixed uncountable cardinal to itself, or to a finite set of size two. These correspond to generalized notions of the universal Baire space (functions from natural numbers to themselves with the product topology) and the Cantor space (functions from natural numbers to the {0,1}-set) respectively. We generalize the notion of Borel sets in three different ways and study the corresponding Borel structures with the aims of generalizing classical theorems of descriptive set theory or providing counter examples. In particular we are interested in equivalence relations on these spaces and their Borel reducibility to each other. The last chapter shows, using game-theoretic techniques, that the order of Borel equivalence relations under Borel reduciblity has very high complexity. The techniques in the above described set theoretical side of the thesis include forcing, general topological notions such as meager sets and combinatorial games of infinite length. By coding uncountable models to functions, we are able to apply the understanding of the generalized descriptive set theory to the model theory of uncountable models. The links between the theorems of model theory (including Shelah's classification theory) and the theorems in pure set theory are provided using game theoretic techniques from Ehrenfeucht-Fraïssé games in model theory to cub-games in set theory. The bottom line of the research declairs that the descriptive (set theoretic) complexity of an isomorphism relation of a first-order definable model class goes in synch with the stability theoretical complexity of the corresponding first-order theory. The first chapter of the thesis has slightly different focus and is purely concerned with a certain modification of the well known Ehrenfeucht-Fraïssé games. There we (me and my supervisor Tapani Hyttinen) answer some natural questions about that game mainly concerning determinacy and its relation to the standard EF-game

A Deterministic Optimization Approach to Protein Sequence Design Using Continuous Models

Determining the sequence of amino acid residues in a heteropolymer chain of a protein with a given conformation is a discrete combinatorial problem that is not generally amenable for gradient-based continuous optimization algorithms. In this paper we present a new approach to this problem using continuous models. In this modeling, continuous "state functions" are proposed to designate the type of each residue in the chain. Such a continuous model helps define a continuous sequence space in which a chosen criterion is optimized to find the most appropriate sequence. Searching a continuous sequence space using a deterministic optimization algorithm makes it possible to find the optimal sequences with much less computation than many other approaches. The computational efficiency of this method is further improved by combining it with a graph spectral method, which explicitly takes into account the topology of the desired conformation and also helps make the combined method more robust. The continuous modeling used here appears to have additional advantages in mimicking the folding pathways and in creating the energy landscapes that help find sequences with high stability and kinetic accessibility. To illustrate the new approach, a widely used simplifying assumption is made by considering only two types of residues: hydrophobic (H) and polar (P). Self-avoiding compact lattice models are used to validate the method with known results in the literature and data that can be practically obtained by exhaustive enumeration on a desktop computer. We also present examples of sequence design for the HP models of some real proteins, which are solved in less than five minutes on a single-processor desktop computer Some open issues and future extensions are noted.

The crystal structure of the dipotassium salt of uridine 5'-diphosphate

C9H12N2Ot2P22-. 2K + .3H20 is orthorhombic, P2~2~2p with a = 18.977 (5), b - 22.597 (6), c = 8.995 (2) A, Z = 8. The structure was refined to R = 0.059 for 2587 observed reflexions. The two molecules of the asymmetric unit have very similar conformations with a 2'- endo sugar pucker and a folded pyrophosphate chain. They form a dimer, coordinated by the K + ions but without direct bridging between the base and the pyrophosphate within each individual molecule. One uracil base has the keto-enol and the other the diketo form. The extended structure shows alternating hydrophobic and hydrophilic regions.

Interaction Energy Based Protein Structure Networks

The three dimensional structure of a protein is formed and maintained by the noncovalent interactions among the amino acid residues of the polypeptide chain These interactions can be represented collectively in the form of a network So far such networks have been investigated by considering the connections based on distances between the amino acid residues Here we present a method of constructing the structure network based on interaction energies among the amino acid residues in the protein We have investigated the properties of such protein energy based networks (PENs) and have shown correlations to protein structural features such as the clusters of residues involved in stability formation of secondary and super secondary structural units Further we demonstrate that the analysis of PENs in terms of parameters such as hubs and shortest paths can provide a variety of biologically important information such as the residues crucial for stabilizing the folded units and the paths of communication between distal residues in the protein Finally the energy regimes for different levels of stabilization in the protein structure have clearly emerged from the PEN analysis

Temperature-Sensitive Src-Transformed Mdck Cells nn 3d Cultures as a Model of Malignant Transformation of Epithelial Cells

The equilibrium between cell proliferation, differentiation, and apoptosis is crucial for maintaining homeostasis in epithelial tissues. In order for the epithelium to function properly, individual cells must gain normal structural and functional polarity. The junctional proteins have an important role both in binding the cells together and in taking part in cell signaling. Cadherins form adherens junctions. Cadherins initiate the polarization process by first recognizing and binding the neighboring cells together, and then guiding the formation of tight junctions. Tight junctions form a barrier in dividing the plasma membranes to apical and basolateral membrane domains. In glandular tissues, single layered and polarized epithelium is folded into tubes or spheres, in which the basal side of the epithelial layer faces the outer basal membrane, and the apical side the lumen. In carcinogenesis, the differentiated architecture of an epithelial layer is disrupted. Filling of the luminal space is a hallmark of early epithelial tumors in tubular and glandular structures. In order for the transformed tumor cells to populate the lumen, enhanced proliferation as well as inhibition of apoptosis is required. Most advances in cancer biology have been achieved by using two-dimensional (2D) cell culture models, in which the cells are cultured on flat surfaces as monolayers. However, the 2D cultures are limited in their capacity to recapitulate the structural and functional features of tubular structures and to represent cell growth and differentiation in vivo. The development of three-dimensional (3D) cell culture methods enables the cells to grow and to be studied in a more natural environment. Despite the wide use of 2D cell culture models and the development of novel 3D culture methods, it is not clear how the change of the dimensionality of culture conditions alters the polarization and transformation process and the molecular mechanisms behind them. Src is a well-known oncogene. It is found in focal and adherens junctions of cultured cells. Active src disrupts cell-cell junctions and interferes with cell-matrix binding. It promotes cell motility and survival. Src transformation in 2D disrupts adherens junctions and the fibroblastic phenotype of the cells. In 3D, the adherens junctions are weakened, and in glandular structures, the lumen is filled with nonpolarized vital cells. Madin-Darby canine kidney (MDCK) cells are an epithelial cell type commonly used as a model for cell polarization. Its-src-transformed variants are useful model systems for analyzing the changes in cell morphology, and they play a role in src-induced malignant transformation. This study investigates src-transformed cells in 3D cell cultures as a model for malignant transformation. The following questions were posed. Firstly: What is the role of the composition and stiffness of the extracellular matrix (ECM) on the polarization and transformation of ts v-src MDCK cells in 3D cell cultures? Secondly: How do the culture conditions affect gene expression? What is the effect of v-src transformation in 2D and in 3D cell models? How does the shift from 2D to 3D affect cell polarity and gene expression? Thirdly: What is the role of survivin and its regulator phosphatase and tensin homolog protein (PTEN) in cell polarization and transformation, and in determining cell fate? How does their expression correlate with impaired mitochondrial function in transformed cells? In order to answer the above questions, novel methods of culturing and monitoring cells had to be created: novel 3D methods of culturing epithelial cells were engineered, enabling real time monitoring of a polarization and transformation process, and functional testing of 3D cell cultures. Novel 3D cell culture models and imaging techniques were created for the study. Attention was focused especially on confocal microscopy and live-cell imaging. Src-transformation disturbed the polarization of the epithelium by disrupting cell adhesion, and sensitized the cells to their environment. With active src, the morphology of the cell cluster depended on the composition and stiffness of the matrix. Gene expression studies revealed a broader impact of src transformation than mere continuous activity of src-kinase. In 2D cultures, src transformation altered the expression of immunological, actin cytoskeleton and extracellular matrix (ECM). In 3D, the genes regulating cell division, inhibition of apoptosis, cell metabolism, mitochondrial function, actin cytoskeleton and mechano-sensing proteins were altered. Surprisingly, changing the culture conditions from 2D to 3D affected also gene expression considerably. The microarray hit survivin, an inhibitor of apoptosis, played a crucial role in the survival and proliferation of src-transformed cells.

Interactions of tris buffer with nucleotides: the crystal structure of tris(hydroxymethyl)methylammonium adenosine 5'-diphosphate dihydrate

The crystal and molecular structures of the Tris salt of adenosine 5'-diphosphate were determined from X-ray diffraction data. The crystals are monoclinic, space P21, and Z = 2 with a=9.198 (2) A, b=6.894 (1) A, c=18.440 (4) A, and beta = 92.55 (2) degrees. Intensity data were collected on an automated diffractometer. The structure was solved by the heavy-atom technique and refined by least squares to R = 0.047. The ADP molecule adopts a folded conformation. The conformation about the glycosidic bond is anti. The conformation of the ribose ring is close to a perfect C(2')-endo-C-(3')-exo puckering. The conformation about C(4')-C(5') is gauche-gauche, similar to other nucleotide structures. The pyrophosphate chain displays a nearly eclipsed geometry when viewed down the P-P vector, unlike the staggered conformation observed in crystal structures of other pyrophosphates. The less favorable eclipsed conformation probably results from the observed association of Tris molecules with the polar diphosphate chain through electrostatic interactions and hydrogen bonds. Such interactions may play an important role in Tris-buffered aqueous solutions of nucleotides and metal ions.

Temperature dependent free energy surface of polymer folding from equilibrium and quench studies

Langevin dynamics simulation studies have been employed to calculate the temperature dependent free energy surface and folding characteristics of a 500 monomer long linear alkane (polyethylene) chain with a realistic interaction potential. Both equilibrium and temperature quench simulation studies have been carried out. Using the shape anisotropy parameter (S) of the folded molecule as the order parameter, we find a weakly first order phase transition between the high-temperature molten globule and low-temperature rodlike crystalline states separated by a small barrier of the order of k(B)T. Near the melting temperature (580 K), we observe an intriguing intermittent fluctuation with pronounced ``1/f noise characteristics'' between these two states with large difference in shape and structure. We have also studied the possibilities of different pathways of folding to states much below the melting point. At 300 K starting from the all-trans linear configuration, the chain folds stepwise into a very regular fourfold crystallite with very high shape anisotropy. Whereas, when quenched from a high temperature (900 K) random coil regime, we identify a two step transition from the random coiled state to a molten globulelike state and, further, to a anisotropic rodlike state. The trajectory reveals an interesting coupling between the two order parameters, namely, radius of gyration (R-g) and the shape anisotropy parameter (S). The rodlike final state of the quench trajectory is characterized by lower shape anisotropy parameter and significantly larger number of gauche defects as compared to the final state obtained through equilibrium simulation starting from all-trans linear chain. The quench study shows indication of a nucleationlike pathway from the molten globule to the rodlike state involving an underlying rugged energy landscape. (C) 2010 American Institute of Physics. doi:10.1063/1.3509398]

An optimal parallel algorithm for sparse matrix bandwidth reduction on a linear array

A simple and efficient algorithm for the bandwidth reduction of sparse symmetric matrices is proposed. It involves column-row permutations and is well-suited to map onto the linear array topology of the SIMD architectures. The efficiency of the algorithm is compared with the other existing algorithms. The interconnectivity and the memory requirement of the linear array are discussed and the complexity of its layout area is derived. The parallel version of the algorithm mapped onto the linear array is then introduced and is explained with the help of an example. The optimality of the parallel algorithm is proved by deriving the time complexities of the algorithm on a single processor and the linear array.

Conformational Polymorphism in Telomeric Structures: Loop Orientation and Interloop Pairing in d(G4TnG4),

Sequence repeats constituting the telomeric regions of chromosomes are known to adopt a variety of unusual structures, consisting of a G tetraplex stem and short stretches of thymines or thymines and adenines forming loops over the stem. Detailed model building and molecular mechanics studies have been carried out for these telomeric sequences to elucidate different types of loop orientations and possible conformations of thymines in the loop. The model building studies indicate that a minimum of two thymines have to be interspersed between guanine stretches to form folded-back structures with loops across adjacent strands in a G tetraplex (both over the small as well as large groove), while the minimum number of thymines required to build a loop across the diagonal strands in a G tetraplex is three. For two repeat sequences, these hairpins, resulting from different types of folding, can dimerize in three distinct ways-i.e., with loops across adjacent strands and on same side, with loops across adjacent strands and on opposite sides, and with loops across diagonal strands and on opposite sides-to form hairpin dimer structures. Energy minimization studies indicate that all possible hairpin dimers have very similar total energy values, though different structures are stabilized by different types of interactions. When the two loops are on the same side, in the hairpin dimer structures of d(G(4)T(n)G(4)), the thymines form favorably stacked tetrads in the loop region and there is interloop hydrogen bonding involving two hydrogen bonds for each thymine-thymine pair. Our molecular mechanics calculations on various folded-back as well as parallel tetraplex structures of these telomeric sequences provide a theoretical rationale for the experimentally observed feature that the presence of intervening thymine stretches stabilizes folded-back structures, while isolated stretches of guanines adopt a parallel tetraplex structure

Conformational Polymorphism in G-tetraplex Structures: Strand Reversal by Base Flipover or Sugar Flipover

Guanine rich sequences adopt a variety of four stranded structures, which differ in strand orientation and conformation about the glycosidic bond even though they are all stabilised by Hoogsteen hydrogen bonded guanine tetrads. Detailed model building and molecular mechanics calculations have been carried out to investigate various possible conformations of guanines along a strand and different possible orientations of guanine strands in a G-tetraplex structure. It is found that for an oligo G stretch per se, a parallel four stranded structure with all guanines in anti conformation is favoured over other possible tetraplex structures. Hence an alternating syn-anti arrangement of guanines along a strand is likely to occur only in folded back tetraplex structures with antiparallel G strands. Our study provides a theoretical rationale for the observed alternation of glycosidic conformation and the inverted stacking arrangement arising from base flipover, in antiparallel G-tetraplex structures and also highlights the various structural features arising due to different types of strand orientations. The molecular mechanics calculations help in elucidating the various interactions which stabilize different G-tetraplex structures and indicate that screening of phosphate charge by counterions could have a dramatic effect on groove width in these four stranded structures.

Stability of synchronization in a multi-cellular system

Networks of biochemical reactions regulated by positive-and negative-feedback processes underlie functional dynamics in single cells. Synchronization of dynamics in the constituent cells is a hallmark of collective behavior in multi-cellular biological systems. Stability of the synchronized state is required for robust functioning of the multi-cell system in the face of noise and perturbation. Yet, the ability to respond to signals and change functional dynamics are also important features during development, disease, and evolution in living systems. In this paper, using a coupled multi-cell system model, we investigate the role of system size, coupling strength and its topology on the synchronization of the collective dynamics and its stability. Even though different coupling topologies lead to synchronization of collective dynamics, diffusive coupling through the end product of the pathway does not confer stability to the synchronized state. The results are discussed with a view to their prevalence in biological systems. Copyright (C) EPLA, 2010