Interaction of substrate uridyl 3',5'-adenosine with ribonuclease A:a molecular dynamics study

A wealth of information available from x-ray crystallographic structures of enzyme-ligand complexes makes it possible to study interactions at the molecular level. However, further investigation is needed when i) the binding of the natural substrate must be characterized, because ligands in the stable enzyme-ligand complexes are generally inhibitors or the analogs of substrate and transition state, and when ii) ligand binding is in part poorly characterized. We have investigated these aspects i? the binding of substrate uridyl 3',5'-adenosine (UpA) to ribonuclease A (RNase A). Based on the systematically docked RNase A-UpA complex resulting from our previous study, we have undertaken a molecular dynamics simulation of the complex with solvent molecules. The molecular dynamics trajectories of this complex are analyzed to provide structural explanations for varied experimental observations on the ligand binding at the B2 subsite of ribonuclease A. The present study suggests that B2 subsite stabilization can be effected by different active site groups, depending on the substrate conformation. Thus when adenosine ribose pucker is O4'-endo, Gln69 and Glu111 form hydrogen-bonding contacts with adenine base, and when it is C2'-endo, Asn71 is the only amino acid residue in direct contact with this base. The latter observation is in support of previous mutagenesis and kinetics studies. Possible roles for the solvent molecules in the binding subsites are described. Furthermore, the substrate conformation is also examined along the simulation pathway to see if any conformer has the properties of a transition state. This study has also helped us to recognize that small but concerted changes in the conformation of the substrate can result in substrate geometry favorable for 2',3' cyclization. The identified geometry is suitable for intraligand proton transfer between 2'-hydroxyl and phosphate oxygen atom. The possibility of intraligand proton transfer as suggested previously and the mode of transfer before the formation of cyclic intermediate during transphosphorylation are discussed.

Structure-Based Phylogeny as a Diagnostic for Functional Characterization of Proteins with a Cupin Fold

Background: The members of cupin superfamily exhibit large variations in their sequences, functions, organization of domains, quaternary associations and the nature of bound metal ion, despite having a conserved beta-barrel structural scaffold. Here, an attempt has been made to understand structure-function relationships among the members of this diverse superfamily and identify the principles governing functional diversity. The cupin superfamily also contains proteins for which the structures are available through world-wide structural genomics initiatives but characterized as ``hypothetical''. We have explored the feasibility of obtaining clues to functions of such proteins by means of comparative analysis with cupins of known structure and function. Methodology/Principal Findings: A 3-D structure-based phylogenetic approach was undertaken. Interestingly, a dendrogram generated solely on the basis of structural dissimilarity measure at the level of domain folds was found to cluster functionally similar members. This clustering also reflects an independent evolution of the two domains in bicupins. Close examination of structural superposition of members across various functional clusters reveals structural variations in regions that not only form the active site pocket but are also involved in interaction with another domain in the same polypeptide or in the oligomer. Conclusions/Significance: Structure-based phylogeny of cupins can influence identification of functions of proteins of yet unknown function with cupin fold. This approach can be extended to other proteins with a common fold that show high evolutionary divergence. This approach is expected to have an influence on the function annotation in structural genomics initiatives.

Two-level k-means clustering algorithm for k-tau relationship establishment and linear-time classification

Partitional clustering algorithms, which partition the dataset into a pre-defined number of clusters, can be broadly classified into two types: algorithms which explicitly take the number of clusters as input and algorithms that take the expected size of a cluster as input. In this paper, we propose a variant of the k-means algorithm and prove that it is more efficient than standard k-means algorithms. An important contribution of this paper is the establishment of a relation between the number of clusters and the size of the clusters in a dataset through the analysis of our algorithm. We also demonstrate that the integration of this algorithm as a pre-processing step in classification algorithms reduces their running-time complexity.

[B4O9H2] Cyclic Borate Units as the Building Unit in a Family of Zinc Borate Structures

A solvothermal reaction of ZnO, boric acid (B(OH)(3)), and aliphatic airlines in a water-pyridine mixture gave four zinc borate phases of different dimensionalities: [Zn(B4O8H2)(C3H10N2)], I (one-dimensional); [Zn(B4O8H2)(C3H10N2)] H2O, II (two-dimensional); [Zn(B5O10H3)(C10H24N4)]center dot H2O, III (two-dimensional): and [Zn-2(B8O15H2)(C3H10N2)(2)], IV (three-dimensional). The structures are formed by the connectivity involving polyborate chains and layers with Zn2+ species. In all the compounds, the amine molecules act its file ligand binding either the same or different zn centers. The formation of two different structures, II and IV, from the same amine by varying the reaction time is noteworthy. Transformation studies on II indicate that the formation of IV. from II, is facile and has been investigated for the first time. Two of file compounds, I and III, exhibit activity for second-order nonlinear optical behavior. The UV exposure of the sample indicates the absorption of all the UV radiation suggesting that the zinc borate compounds could be exploited for UV-blocking applications. The compounds have been characterized by powder X-ray diffraction, infrared spectroscopy, thermogravimetric analysis, UV-vis, photoluminescence, and NMR studies.

Mycobacterial Stress Regulation: The Dps ``Twin Sister'' Defense Mechanism and Structure-Function Relationship

In this work, we have tried to emphasize the connection between mycobacterial growth and regulation of gene expression. Utilization of multiple carbon sources and diauxic growth helps bacteria to regulate gene expression at an optimum level so that the inhospitable conditions encountered during nutrient depletion can be circumvented. These aspects will be discussed with respect to mycobacterial growth in subsequent sections. Identification and characterization of genes induced under such conditions is helpful to understand the physiology of the bacterium. Although it is necessary to compare the total expression profile of proteins as they transit from vegetative growth to stationary phase, at times a lot of insights can be deciphered from the expression pattern of one or two proteins. We have compared the protein expression and sigma factor selectivity of two such proteins in M. smegmatis to understand the differential regulation of genes playing diverse function in the same species. Some newer insights on the structure and function of one of the Dps proteins are also explained.

Allostery and conformational free energy changes in human tryptophanyl-tRNA synthetase from essential dynamics and structure networks

The interdependence of the concept of allostery and enzymatic catalysis, and they being guided by conformational mobility is gaining increased prominence. However, to gain a molecular level understanding of llostery and hence of enzymatic catalysis, it is of utter importance that the networks of amino acids participating in allostery be deciphered. Our lab has been exploring the methods of network analysis combined with molecular dynamics simulations to understand allostery at molecular level. Earlier we had outlined methods to obtain communication paths and then to map the rigid/flexible regions of proteins through network parameters like the shortest correlated paths, cliques, and communities. In this article, we advance the methodology to estimate the conformational populations in terms of cliques/communities formed by interactions including the side-chains and then to compute the ligand-induced population shift. Finally, we obtain the free-energy landscape of the protein in equilibrium, characterizing the free-energy minima accessed by the protein complexes. We have chosen human tryptophanyl-tRNA synthetase (hTrpRS), a protein esponsible for charging tryptophan to its cognate tRNA during protein biosynthesis for this investigation. This is a multidomain protein exhibiting excellent allosteric communication. Our approach has provided valuable structural as well as functional insights into the protein. The methodology adopted here is highly generalized to illuminate the linkage between protein structure networks and conformational mobility involved in the allosteric mechanism in any protein with known structure.

Random network behaviour of protein structures

Geometric and structural constraints greatly restrict the selection of folds adapted by protein backbones, and yet, folded proteins show an astounding diversity in functionality. For structure to have any bearing on function, it is thus imperative that, apart from the protein backbone, other tunable degrees of freedom be accountable. Here, we focus on side-chain interactions, which non-covalently link amino acids in folded proteins to form a network structure. At a coarse-grained level, we show that the network conforms remarkably well to realizations of random graphs and displays associated percolation behavior. Thus, within the rigid framework of the protein backbone that restricts the structure space, the side-chain interactions exhibit an element of randomness, which account for the functional flexibility and diversity shown by proteins. However, at a finer level, the network exhibits deviations from these random graphs which, as we demonstrate for a few specific examples, reflect the intrinsic uniqueness in the structure and stability, and perhaps specificity in the functioning of biological proteins.

Hydrothermal synthesis and structural characterization of novel alpha-Keggin unit-supported Cu(II)- and Mn(II)-bipyridine complexes from a tri-lacunary precursor

Blue [{Cu(2,2'-bipy)(2)}(2){alpha-SiW12O40}] (bipy = bipyridyl) (1) and pale yellow [Mn(2,2'-bipy)(3)](2)[alpha-SiW12O40] (2) have been synthesized hydrothermally and characterized by IR spectroscopy and single crystal X-ray structure analysis. In 1, the [alpha-SiW12O40](4-) ion acts as a bridge between the two [{Cu(2,2'-bipy)(2)](2+) moieties via coordination through the terminal oxygen atoms, while in 2, the [Mn(2,2'-bipy)(3)](2+) ion balances the charge on the polyoxo anion without forming any covalent bond. To the best of our knowledge, this is the first example of transition metal-mediated transformation of [alpha-SiW9O34](10-) to [alpha-SiW12O40](4-).

Inter-helical Interactions in Membrane Proteins: Analysis Based on the Local Backbone Geometry and the Side Chain Interactions

The availability of a significant number of the Structures of helical membrane proteins has prompted us to investigate the mode of helix-helix packing. In the present study, we have considered a dataset of alpha-helical membrane proteins representing Structures solved from all the known superfamilies. We have described the geometry of all the helical residues in terms of local coordinate axis at the backbone level. Significant inter-helical interactions have been considered as contacts by weighing the number of atom-atom contacts, including all the side-chain atoms. Such a definition of local axis and the contact criterion has allowed us to investigate the inter-helical interaction in a systematic and quantitative manner. We show that a single parameter (designated as alpha), which is derived from the parameters representing the Mutual orientation of local axes, is able to accurately Capture the details of helix-helix interaction. The analysis has been carried Out by dividing the dataset into parallel, anti-parallel, and perpendicular orientation of helices. The study indicates that a specific range of alpha value is preferred for interactions among the anti-parallel helices. Such a preference is also seen among interacting residues of parallel helices, however to a lesser extent. No such preference is seen in the case of perpendicular helices, the contacts that arise mainly due to the interaction Of Surface helices with the end of the trans-membrane helices. The Study Supports the prevailing view that the anti-parallel helices are well packed. However, the interactions between helices of parallel orientation are non-trivial. The packing in alpha-helical membrane proteins, which is systematically and rigorously investigated in this study, may prove to be useful in modeling of helical membrane proteins.

Correlation between conductivity or diffusivity and activation energy in amorphous solids

There exist many investigations of ionic transport in a variety of glasses. These studies exhibit strong correlation between ionic conductivity and activation energy: Typically, it is found that higher conductivity is associated with lower activation energies and vice versa. Although there are explanations for this at a phenomenological level, there is no consistent physical picture to explain the correlation between conductivity and activation energy. We have carried out molecular dynamics simulation as a function of the size of the impurity atom or diffusant (both neutral and charged) in a host amorphous matrix. We find that there is a maximum in self-diffusivity as a function of the size of the impurity atom suggesting that there is an appropriate size for which the diffusivity is maximum. The activation energy is found to be the lowest for this size of the impurity. A similar maximum has been previously found in other condensed phases, such as confined fluids and dense liquids, and has its origin in the levitation effect. The implications of this result for understanding ionic conductivity in glasses are discussed. Our results suggest that there is a relation between microscopic structure of the amorphous solid, diffusivity or conductivity, and activation energy. The nature of this relationship is discussed in terms of the levitation parameter showing that diffusivity is maximum when the size of the neck or doorway radius is comparable with the size of the diffusant. Our computational results here are in excellent agreement with independent experimental results of Nascimento et al. [Braz. J. Phys. 35, 626 (2005)] that structural features of the glass are important in determining the ionic conductivity.

Quasi-periodic, global oscillations in sea level pressure on intraseasonal timescales

The sea level pressure (SLP) variability in 30-60 day intraseasonal timescales is investigated using 25 years of reanalysis data addressing two issues. The first concerns the non-zero zonal mean component of SLP near the equator and its meridional connections, and the second concerns the fast eastward propagation (EP) speed of SLP compared to that of zonal wind. It is shown that the entire globe resonates with high amplitude wave activity during some periods which may last for few to several months, followed by lull periods of varying duration. SLP variations in the tropical belt are highly coherent from 25A degrees S to 25A degrees N, uncorrelated with variations in mid latitudes and again significantly correlated but with opposite phase around 60A degrees S and 65A degrees N. Near the equator (8A degrees S-8A degrees N), the zonal mean contributes significantly to the total variance in SLP, and after its removal, SLP shows a dominant zonal wavenumber one structure having a periodicity of 40 days and EP speeds comparable to that of zonal winds in the Indian Ocean. SLP from many of the atmospheric and coupled general circulation models show similar behaviour in the meridional direction although their propagation characteristics in the tropical belt differ widely.

Critical test for Altshuler-Aronov theory: Evolution of the density of states singularity in double perovskite Sr2FeMoO6 with controlled disorder

With high-resolution photoemission spectroscopy measurements, the density of states (DOS) near the Fermi level (E-F) of double perovskite Sr2FeMoO6 having different degrees of Fe/Mo antisite disorder has been investigated with varying temperature. The DOS near E-F showed a systematic depletion with increasing degree of disorder, and recovered with increasing temperature. Altshuler-Aronov (AA) theory of disordered metals well explains the dependences of the experimental results. Scaling analysis of the spectra provides experimental indication for the functional form of the AA DOS singularity.

Electrochemical Functionalization of a Gold Electrode with Redox-Active Self-Assembled Monolayer for Electroanalytical Application

The electrochemical functionalization of a Au electrode with a redox-active monolayer and the electroanalytical applications of the functionalized electrode are described. Reaction of the electrochemically derived o-quinone on the self-assembled monolayer (SAM) of 6-mercaptopurine (MPU) on a Au electrode gives a redox-active 4-(6-mercapto-purin-9-yl)benzene-1,2-diol (MPBD) self-assembly under optimized conditions. Electrochemical quartz crystal microbalance technique has been employed to follow the functionalization of the electrode in real time. Electrochemically derived o-quinone reacts at the N(9) position of the self-assembled MPU in neutral pH. Raman spectral measurement confirms the reaction of o-quinone on MPU self-assembly. MPBD shows a well-defined reversible redox response, characteristic of a surface-confined redox mediator at 0.21 V in neutral pH. The anodic peak potential (Epa) of MPBD shifts by −60 mV while changing the solution pH by 1 unit, indicating that the redox reaction involves two electrons and two protons. The surface coverage (Γ) of MPBD was 7.2 ± 0.3 × 10-12 mol/cm2. The apparent heterogeneous rate constant (ksapp) for MPBD was 268 ± 6 s-1. MPBD efficiently mediates the oxidation of nicotinamide adenine dinucleotide (NADH) and ascorbate (AA). A large decrease in the overpotential and significant increase in the peak current with respect to the unmodified electrode has been observed. Surface-confined MPBD has been successfully used for the amperometric sensing of NADH and AA in neutral pH at the nanomolar level.

Investigation of inter-ion interactions in N,N,N',N'-tetramethylethylenediammonium dithiocyanate via experimental and theoretical charge density studies

The crystal structure of the N,N,N',N'-tetramethylethylenediammonium dithiocyanate salt has been examined by experimental charge density studies from high-resolution X-ray diffraction data. The corresponding results are compared with multipole refinements, using theoretical structure factors obtained from a periodic density functional theory calculation at the B3LYP level with a 6-31G** basis set. The salt crystallizes in space group P (1) over bar and contains only a single ion pair with an inversion center in the cation. The salt has thus one unique classical N+-H center dot center dot center dot(NCS)(-) hydrogen bond but also has six other weaker interactions: four C-H center dot center dot center dot S, one C-H center dot center dot center dot N, and one C-H center dot center dot center dot C-pi. The nature of all these interactions has been examined topologically using Bader's quantum theory of "atoms in molecules" and all eight of the Koch-Popelier criteria. The experimental and theoretical approaches agree well and both show that the inter-ion interactions, even in this simplest of systems, play an integrated and complex role in the packing of the ions in the crystal. Electrostatic potential maps are derived from experimental charge densities. This is the first time such a system has been examined in detail by these methods.

Hybrid algorithms for power system unit commitment

Four hybrid algorithms has been developed for the solution of the unit commitment problem. They use simulated annealing as one of the constituent techniques, and produce lower cost schedules; two of them have less overhead than other soft computing techniques. They are also more robust to the choice of parameters. A special technique avoids the generating of infeasible schedules, and thus reduces computation time.