Towards a temperature-guided molecular switch: an unusual reversible low-temperature polymorphic phase transition in a conformational locked environment

A conformationally locked tetraacetate undergoes, quite akin to a temperature-guided molecular switch, a reversible thermal switching between two polymorphic modifications; the room-temperature alpha-form converted at -4 degrees C to a low-temperature denser beta-form, which displayed an unusual kinetic stability till 67 degrees C and transformed back to the alpha-form beyond this temperature.

Unusual structural stability and ligand induced alterations in oligomerization of a galectin

Abstract L-14, a 14-kDa S-type lectin shows the jelly roll tertiary structural fold akin to legume lectins yet, unlike them, it does not dissociate on thermal unfolding. In the absence of ligand L-14 displays denaturation transitions corresponding to tetrameric and octameric entities. The presence of complementary ligand reduces the association of L-14, which is in stark contrast with legume lectins where no alterations in quaternary structures are brought about by saccharides. From the magnitude of the increase in denaturation temperature induced by disaccharides the binding constants calculated from differential scanning calorimetry are comparable with those extrapolated from titration calorimetry indicating that L-14 interacts with ligands essentially in the folded state.

Stability of assembled epitopes on human chorionic gonadotropin to covalent modifications

A measure of stability of a given epitope is an important parameter in the exploration of the utility of a desired MAb. It defines the conditions necessary for using MAbs as an investigative tool in several research methodologies and therapeutic protocols. Despite these obvious interests the lack of simple and rapid assay systems for quantitating MAb-Ag interactions has largely hampered these studies. A single step MAb-Solid Phase Radioimmunoassay (SS-SPRIA), is described which eliminates errors that may arise with multistep sandwich assays. SS-SPRIA has been used to demonstrate the differential stability of the assembled epitopes on gonadotropins. Differential stability towards specific reagents can be exploited to identify aminoacid residues at the epitopic site. Inactivation of an epitopic region is indicative of the presence of the group modified, provided conformational relaxations are not induced due to modifications at distant sites. Here we provide evidence to validate these conclusions.

Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

Conformational study of valinomycin: a molecular dynamics approach

Valinomycin is a highly flexible cyclic dodecadepsipeptide that transports ions across membranes. Such a flexibility in the conformation is required for its biological function since it has to encounter a variety of environments and liganding state. Exploration of conformational space of this molecule is therefore important and is one of the objectives of the present study that has been carried out by means of high temperature Molecular Dynamics. Further, the stability of the known bracelet-like structure of the uncomplexed valinomycin and the inherent flexibility around this structure has been investigated. The uncomplexed form of valinomycin has been simulated at 75–100 K for 1 ns in order to elucidate the average conformational properties. An alanine-analog of valinomycin has been simulated under identical conditions in order to evaluate the effect of sidechain on the conformational properties, The studies confirm the effect of sidechain on conformational equilibrium.

Insight into the early stages of thermal unfolding of peanut agglutinin by molecular dynamics simulations

Peanut agglutinin is a homotetrameric nonglycosylated protein. The protein has a unique open quaternary structure. Molecular dynamics simulations have been employed follow the atomistic details of its unfolding at different temperatures. The early events of the deoligomerization of the protein have been elucidated in the present study. Simulation trajectories of the monomer as well as those of the tetramer have been compared and the tetramer is found to be substantially more stable than its monomeric counterpart. The tetramer shows retention of most of its.. secondary structure but considerable loss of the tertiary structure at high temperature. e generation of a This observation impies the molten globule-like intermediate in the later stages of deoligomerization. The quaternary structure of the protein has weakened to a large extent, but none of the subunits are separated. In addition, the importance of the metal-binding to the stability of the protein structure has also been investigated. Binding of the metal ions not only enhances the local stability of the metal-ion binding loop, but also imparts a global stability to the overall structure. The dynamics of different interfaces vary significantly as probed through interface clusters. The differences are substantially enhanced at higher temperatures. The dynamics and the stability of the interfaces have been captured mainly by cluster analysis, which has provided detailed information on the thermal deoligomerization of the protein.

Structure and bonding of MCB5H7 and its sandwiched dimer CB5H6M–MCB5H6 (M = Si, Ge, Sn): Isomer stability and preference for slip distorted structure

We find sandwiched metal dimers CB5H6M–MCB5H6 (M = Si, Ge, Sn) which are minima in the potential energy surface with a characteristic M–M single bond. The NBO analysis and the M–M distances (Å) (2.3, 2.44 and 2.81 for M = Si, Ge, Sn) indicate substantial M–M bonding. Formal generation of CB5H6M–MCB5H6 has been studied theoretically. Consecutive substitution of two boron atoms in B7H−27 by M (Si, Ge, Sn) and carbon, respectively followed by dehydrogenation may lead to our desired CB5H6M–MCB5H6. We find that the slip distorted geometry is preferred for MCB5H7 and its dehydrogenated dimer CB5H6M–MCB5H6. The slip-distortion of M–M bond in CB5H6M–MCB5H6 is more than the slip distortion of M–H bond in MCB5H7. Molecular orbital analysis has been done to understand the slip distortion. Larger M–M bending (CB5H6M–MCB5H6) in comparison with M–H bending (MCB5H7) is suspected to be encouraged by stabilization of one of the M–M π bonding MO’s. Preference of M to occupy the apex of pentagonal skeleton of MCB5H7 over its icosahedral analogue MCB10H11 has been observed.

Effect of substitution on molecular conformation and packing features in a series of aryl substituted ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4-tetrahydropyrimidine-5-carboxylates

The supramolecular structures of eight aryl protected ethyl-6-methyl-4-phenyl-2-thioxo-1,2,3,4 tetrahydropyrimidine-5-carboxyl ates were analyzed in order to understand the effect of variations in functional groups on molecular geometry, conformation and packing of molecules in the crystalline lattice. It is observed that the existence of a short intra-molecular C-H center dot center dot center dot pi interaction between the aromatic hydrogen of the aryl ring with the isolated double bond of the six-membered tetrahydropyrimidine ring is a key feature which imparts additional stability to the molecular conformation in the solid state. The compounds pack via the cooperative involvement of both N-H center dot center dot center dot S=C and N-H center dot center dot center dot O=C intermolecular dimers forming a sheet like structure. In addition, weak C-H center dot center dot center dot O and C-H center dot center dot center dot pi intermolecular interactions provide additional stability to the crystal packing.

A network representation of protein structures: Implications for protein stability

This study views each protein structure as a network of noncovalent connections between amino acid side chains. Each amino acid in a protein structure is a node, and the strength of the noncovalent interactions between two amino acids is evaluated for edge determination. The protein structure graphs (PSGs) for 232 proteins have been constructed as a function of the cutoff of the amino acid interaction strength at a few carefully chosen values. Analysis of such PSGs constructed on the basis of edge weights has shown the following: 1), The PSGs exhibit a complex topological network behavior, which is dependent on the interaction cutoff chosen for PSG construction. 2), A transition is observed at a critical interaction cutoff, in all the proteins, as monitored by the size of the largest cluster (giant component) in the graph. Amazingly, this transition occurs within a narrow range of interaction cutoff for all the proteins, irrespective of the size or the fold topology. And 3), the amino acid preferences to be highly connected (hub frequency) have been evaluated as a function of the interaction cutoff. We observe that the aromatic residues along with arginine, histidine, and methionine act as strong hubs at high interaction cutoffs, whereas the hydrophobic leucine and isoleucine residues get added to these hubs at low interaction cutoffs, forming weak hubs. The hubs identified are found to play a role in bringing together different secondary structural elements in the tertiary structure of the proteins. They are also found to contribute to the additional stability of the thermophilic proteins when compared to their mesophilic counterparts and hence could be crucial for the folding and stability of the unique three-dimensional structure of proteins. Based on these results, we also predict a few residues in the thermophilic and mesophilic proteins that can be mutated to alter their thermal stability.

Trishanku, a novel regulator of cell-type stability and morphogenesis in Dictyostelium discoideum

We have identified a novel gene, trishanku (triA), by random insertional mutagenesis of Dictyostelium discoideum. TriA is a Broad complex Tramtrack bric-a-brac domain-containing protein that is expressed strongly during the late G2 phase of cell cycle and in presumptive spore (prespore (psp)) cells. Disrupting triA destabilizes cell fate and reduces aggregate size; the fruiting body has a thick stalk, a lowered spore: stalk ratio, a sub-terminal spore mass and small, rounded spores. These changes revert when the wild-type triA gene is re-expressed under a constitutive or a psp-specific promoter. By using short- and long-lived reporter proteins, we show that in triA(-) slugs the prestalk (pst)/psp proportion is normal, but that there is inappropriate transdifferentiation between the two cell types. During culmination, regardless of their current fate, all cells with a history of pst gene expression contribute to the stalk, which could account for the altered cell-type proportion in the mutant.

X-ray and molecular-dynamics studies on Mycobacterium leprae single-stranded DNA-binding protein and comparison with other eubacterial SSB structures

The crystal structures of two forms of Mycobacterium leprae single-stranded DNA-binding protein (SSB) have been determined at 2.05 and 2.8 A resolution. Comparison of these structures with the structures of other eubacterial SSBs indicates considerable variation in their quaternary association, although the DNA-binding domains in all of them exhibit the same OB-fold. This variation has no linear correlation with sequence variation, but could be related to variation in protein stability. Molecular-dynamics simulations have been carried out on tetrameric molecules derived from the two forms and the prototype Escherichia coli SSB and the individual subunits of both proteins. Together, the X-ray studies and molecular-dynamics simulations yield information on the relatively rigid and flexible regions of the molecule and on the effect of oligomerization on flexibility. The simulations provide insight into the changes in subunit structure on oligomerization. They also provide insight into the stability and time evolution of the hydrogen bonds/water bridges that connect the two pairs of monomers in the tetramer.

Molecular dynamics simulations on parallel and antiparallel C.G*G triplexes

Molecular dynamics (MD) studies have been carried out on the Hoogsteen hydrogen bonded parallel and the reverse Hoogsteen hydrogen banded antiparallel C.G*G triplexes. Earlier, the molecular mechanics studies had shown that the parallel structure was energetically more favourable than the antiparallel structure. To characterize the structural stability of the two triplexes and to investigate whether the antiparallel structure can transit to an energetically more favourable structure, due to the local fluctuations in the structure during the MD simulation, the two structures were subjected to 200ps of constant temperature vacuum MD simulations at 300K. Initially no constraints were applied to the structures and it was observed that for the antiparallel tripler, the structure showed a large root mean square deviation from the starting structure within the first 12ps and the N4-H41-O6 hydrogen bond in the WC duplex got distorted due to a high propeller twist and a moderate increase in the opening angle in the basepairs. Starting from an initial value of 30 degrees, helical twist of the average structure from this simulation had a value of 36 degrees, while the parallel structure stabilized at a twist of 33 degrees. In spite of the hydrogen bond distortions in the antiparallel tripler, it was energetically comparable to the parallel tripler. To examine the structural characteristics of an undistorted structure, another MD simulation was performed on the antiparallel tripler by constraining all the hydrogen bonds. This structure stabilized at an average twist of 33 degrees. In the course of the dynamics though the energy of the molecule - compared to the initial structure - improved, it did not become comparable to the parallel structure. Energy minimization studies performed in the presence of explicit water and counterions also showed the two structures to be equally favourable energetically Together these results indicate that the parallel C.G*G tripler with Hoogsteen hydrogen bonds also represents a stereochemically and energetically favourable structure for this class of triplexes.

Sequence Independent Recombination Triple Helices: A Molecular Dynamics Study

Recent experimental studies have shown that the Rec-A mediated homologous recombination reaction involves a triple helical intermediate, in which the third strand base forms hydrogen bonds with both the bases in the major groove of the Watson-Crick duplex. Such 'mixed' hydrogen bonds allow formation of sequence independent triplexes. DNA triple helices involving 'mixed' hydrogen bonds have been studied, using model building, molecular mechanics (MM) and molecular dynamics (MD). Models were built for a tripler comprising all four possible triplets viz., G.C*C, C.G*G, A.T*T and T.A*A. To check the stability of all the 'mixed' hydrogen bonds in such triplexes and the conformational preferences of such tripler structures, MD studies were carried out starting from two structures with 30 degrees and 36 degrees twist between the basepairs. It was observed that though the two triplexes converged towards a similar structure, the various hydrogen bonds between the WC duplex and the third strand showed differential stabilities. An MD simulation with restrained hydrogen bonds showed that the resulting structure was stable and remained close to the starting structure. These studies help us in defining stable hydrogen bond geometries involving the third strand and the WC duplex. It was observed that in the C.G*G triplets the N7 atom of the second strand is always involved in hydrogen bonding. In the G.C*C triplets, either N3 or O2 in the third strand cytosine can interchangeably act as a hydrogen bond acceptor.

Interactions of L-serine at the active site of serine hydroxymethyltransferases: induction of thermal stability

Serine hydroxymethyltransferase (SHMT), EC 2.1.2.1, exhibits broad substrate and reaction specificity. In addition to cleaving many 3-hydroxyamino acids to glycine and an aldehyde, the enzyme also catalyzed the decarboxylation, transamination and racemization of several substrate analogues of amino acids. To elucidate the mechanism of interaction of substrates, especially L-serine with the enzyme, a comparative study of interaction of L-serine with the enzyme from sheep liver and Escherichia coli, was carried out. The heat stability of both the enzymes was enhanced in the presence of serine, although to different extents. Thermal denaturation monitored by spectral changes indicated an alteration in the apparent T, of sheep liver and E. coli SHMTs from 55 +/- 1 degrees C to 72 +/- 3 degrees C at 40 mM serine and from 67 +/- 1 degrees C to 72 +/- 1 degrees C at 20 mM serine, respectively. Using stopped flow spectrophotometry k values of (49 +/- 5)(.)10(-3) s(-1) and (69 +/- 7).10(-3) s(-1) for sheep liver and E. coli enzymes were determined at 50 mM serine. The binding of serine monitored by intrinsic fluorescence and sedimentation velocity measurements indicated that there was no generalized change in the structure of both proteins. However, visible CD measurements indicated a change in the asymmetric environment of pyridoxal 5'-phosphate at the active site upon binding of serine to both the enzymes. The formation of an external aldimine was accompanied by a change in the secondary structure of the enzymes monitored by far UV-CD spectra. Titration microcalorimetric studies in the presence of serine (8 mM) also demonstrated a single class of binding and the conformational changes accompanying the binding of serine to the enzyme resulted in a more compact structure leading to increased thermal stability of the enzyme.

Structure and stability of human telomeric sequence

Telomeric DNA of a variety of vertebrates including humans contains the tandem repeat d(TTAGGG)(n). We have investigated the structural properties of the human telomeric repeat oligonucleotide models d(T(2)AG(3))(4), d(G(3)T(2)A)(3)G(3), and d(G(3)T(2)AG(3)) using CD, gel electrophoresis, and chemical probing techniques. The sequences d(G(3)T(2)A)(3)G(3) and d(T(2)AG(3))(4) assume an antiparallel G quartet structure by intramolecular folding, while the sequence d(G(3)T(2)AG(3)) also adopts an antiparallel G quartet structure but by dimerization of hairpins. In all the above cases, adenines are in the loop. The TTA loops are oriented at the same end of the G tetrad stem in the case of hairpin dimer. Further, the oligonucleotide D(G(3)T(2)AG(3)) forms a higher order structure by the association of two hairpin dimers via stacking of G tetrad planes. Here we show that N-7 of adenine in the hairpin dimer is Hoogsteen hydrogen-bonded. The partial reactivity of loop adenines with DEPC in d(T(2)AG(3))(4) suggests that the intramolecular G quartet structure is highly polymorphic and structures with different loop orientations and topologies are formed in solution. Intra- and interloop hydrogen bonding schemes for the TTA loops are proposed to account for the observed diethyl pyrocarbonate reactivities of adenines. Sodium-induced G quartet structures differ from their potassium-induced counterparts not only in stability but also in loop conformation and interactions. Thus, the overall structure and stability of telomeric sequences are modulated by the cation present, loop sequence, and the number of G tracts, which might be important for the telomere function.