Basis set dependence of the molecular electrostatic potential topography. A case study of substituted benzenes

The basis set dependence of the topographical structure of the molecular electrostatic potential (MESP), as well as the effect of substituents on the MESP distribution, has been investigated with substituted benzenes as test cases. The molecules are studied at HF-SCF 3�21G and 6�31G** levels, with a further MESP topographical investigation at the 3�21G, double-zeta, 6�31G*, 6�31G**, double-zeta polarized and triple-zeta polarized levels. The MESP critical points for a 3�21G optimized/6�31G** basis are similar to the corresponding 6�31G** optimized/6�31G** ones. More generally, the qualitative features of the MESP topography computed at the polarized level are independent of the level at which optimization is carried out. For a proper representation of oxygen lone pairs, however, optimization using a polarized basis set is required. The nature of the substituent drastically changes the MESP distribution over the phenyl ring. The values and positions of MESP minima indicate the most active site for electrophilic attack. This point is strengthened by a study of disubstituted benzenes.

Physics based basis function for vibration analysis of high speed rotating beams

The natural frequencies of continuous systems depend on the governing partial differential equation and can be numerically estimated using the finite element method. The accuracy and convergence of the finite element method depends on the choice of basis functions. A basis function will generally perform better if it is closely linked to the problem physics. The stiffness matrix is the same for either static or dynamic loading, hence the basis function can be chosen such that it satisfies the static part of the governing differential equation. However, in the case of a rotating beam, an exact closed form solution for the static part of the governing differential equation is not known. In this paper, we try to find an approximate solution for the static part of the governing differential equation for an uniform rotating beam. The error resulting from the approximation is minimized to generate relations between the constants assumed in the solution. This new function is used as a basis function which gives rise to shape functions which depend on position of the element in the beam, material, geometric properties and rotational speed of the beam. The results of finite element analysis with the new basis functions are verified with published literature for uniform and tapered rotating beams under different boundary conditions. Numerical results clearly show the advantage of the current approach at high rotation speeds with a reduction of 10 to 33% in the degrees of freedom required for convergence of the first five modes to four decimal places for an uniform rotating cantilever beam.

Structural basis for poor uracil excision from hairpin DNA

Two-dimensional NMR and molecular dynamics simulations have been used to determine the three-dimensional structures of two hairpin DNA structures: d-CTAGAG GATCCUTTTGGATCCT (abbreviated as U1-hairpin) and d-CTAGAGGATCCTTUTGGATCCT (abbreviated as U3-hairpin). The (1) H resonances of both of these hairpin structures have been assigned almost completely. NMR restrained molecular dynamics and energy minimization procedures have been used to describe the three-dimensional structures of these hairpins. This study and concurrent NMR structural studies on two other d-CTAGAGGA TCCTUTTGGATCCT (abbreviated as U2-hairpin) and d-CTAGAGGATCCTTTUGGATCCT (abbreviated as U4-hairpin) have shed light upon various interactions reported between Echerichia coli uracil DNA glycosylase (UDG) and uracil-containing DNA. The backbone torsion angles, which partially influence the local conformation of U12 and U14 in U1 and U3-hairpins, respectively, are probably locked in the trans conformation as in the case of U-13 in the U2-hairpin. Such a stretched-out backbone conformation in the vicinity of U-12 and U-14 is thought to be the reason why the K-m value is poor for U1- and U3-hairpins as it is for the U2-hairpin. Furthermore, the bases U-12 and U-14 in both U1- and U3-hairpins adopt an anti conformation, in contrast with the base conformation of U-13 in the U2-hairpin, which adopts a syn conformation. The clear discrepancy observed in the U-base orientation with respect to the sugar moieties could explain why the V-max value is 10- to 20-fold higher for the U1- and U3-hairpins compared with the U2-hairpin. Taken together, these observations support our interpretation that the unfavourable backbone results in a poor K-m value, whereas the unfavourable nucleotide conformation results in a poor V-max value. These two parameters therefore make the U1- and U3-hairpins better substrates for UDG compared with the U2-hairpin, as reported earlier [Kumar, N. V. & Varshney, U. (1997) Nucleic Acids Res. 25, 2336-2343.].

A Grobner basis for certain affine monomial curves

Let K be a field and let m(0),...,m(e-1) be a sequence of positive integers. Let W be a monomial curve in the affine e-space A(K)(e), defined parametrically by X-0 = T-m0,...,Xe-1 = Tme-1 and let p be the defining ideal of W. In this article, we assume that some e-1 terms of m(0), m(e-1) form an arithmetic sequence and produce a Grobner basis for p.

Conformational basis for substrate recognition and regulation of catalytic activity in Staphylococcus aureus nucleoside di-phosphate kinase

Nucleoside diphosphate kinases (NDK) are characterized by high catalytic turnover rates and diverse substrate specificity. These features make this enzyme an effective activator of a pro-drug an application that has been actively pursued for a variety of therapeutic strategies. The catalytic mechanism of this enzyme is governed by a conserved histidine that coordinates a magnesium ion at the active site. Despite substantial structural and biochemical information on NDK, the mechanistic feature of the phospho-transfer that leads to auto-phosphorylation remains unclear. While the role of the histidine residue is well documented, the other active site residues, in particular the conserved serine remains poorly characterized. Studies on some homologues suggest no role for the serine residue at the active site, while others suggest a crucial role for this serine in the regulation and quaternary association of this enzyme in some species. Here we report the biochemical features of the Staphylococcus aureus NDK and the mutant enzymes. We also describe the crystal structures of the apo-NDK, as a transition state mimic with vanadate and in complex with different nucleotide substrates. These structures formed the basis for molecular dynamics simulations to understand the broad substrate specificity of this enzyme and the role of active site residues in the phospho-transfer mechanism and oligomerization. Put together, these data suggest that concerted changes in the conformation of specific residues facilitate the stabilization of nucleotide complexes thereby enabling the steps involved in the ping-pong reaction mechanism without large changes to the overall structure of this enzyme. (C) 2011 Elsevier B.V. All rights reserved.

Structural basis for the functional and inhibitory mechanisms of beta-hydroxyacyl-acyl carrier protein dehydratase (FabZ) of Plasmodium falciparum

The beta-hydroxyacyl-acyl carrier protein dehydratase of Plasmodium falciparum (PfFabZ) catalyzes the third and important reaction of the fatty acid elongation cycle. The crystal structure of PfFabZ is available in hexameric (active) and dimeric (inactive) forms. However, PfFabZ has not been crystallized with any bound inhibitors until now. We have designed a new condition to crystallize PfFabZ with its inhibitors bound in the active site, and determined the crystal structures of four of these complexes. This is the first report on any FabZ enzyme with active site inhibitors that interact directly with the catalytic residues. Inhibitor binding not only stabilized the substrate binding loop but also revealed that the substrate binding tunnel has an overall shape of ``U''. In the crystal structures, residue Phe169 located in the middle of the tunnel was found to be in two different conformations, open and closed. Thus, Phe169, merely by changing its side chain conformation, appears to be controlling the length of the tunnel to make it suitable for accommodating longer substrates. The volume of the substrate binding tunnel is determined by the sequence as well as by the conformation of the substrate binding loop region and varies between organisms for accommodating fatty acids of different chain lengths. This report on the crystal structures of the complexes of PfFabZ provides the structural basis of the inhibitory mechanism of the enzyme that could be used to improve the potency of inhibitors against an important component of fatty acid synthesis common to many infectious organisms. (C) 2011 Elsevier Inc. All rights reserved.

Conformational Basis for Substrate Recruitment in Protein Tyrosine Phosphatase 10D

The coordinated activity of protein tyrosine phosphatases (PTPs) is crucial for the initiation, modulation, and termination of diverse cellular processes. The catalytic activity of this protein depends on a nucleophilic cysteine at the active site that mediates the hydrolysis of the incoming phosphotyrosine substrate. While the role of conserved residues in the catalytic mechanism of PTPs has been extensively examined, the diversity in the mechanisms of substrate recognition and modulation of catalytic activity suggests that other, less conserved sequence and structural features could contribute to this process. Here we describe the crystal structures of Drosophila melanogaster PTP10D in the apo form as well as in a complex with a substrate peptide and an inhibitor. These studies reveal the role of aromatic ring stacking interactions at the boundary of the active site of PTPs in mediating substrate recruitment. We note that phenylalanine 76, of the so-called KNRY loop, is crucial for orienting the phosphotyrosine residue toward the nucleophilic cysteine. Mutation of phenylalanine 76 to leucine results in a 60-fold decrease in the catalytic efficiency of the enzyme. Fluorescence measurements with a competitive inhibitor, p-nitrocatechol sulfate, suggest that Phe76 also influences the formation of the enzyme-substrate intermediate. The structural and biochemical data for PTP10D thus highlight the role of relatively less conserved residues in PTP domains in both substrate recruitment and modulation of reaction kinetics.

Composite instantaneous comparators basis for design and their transient performance

Principles of design of composite instantaneous comparators (a combination of amplitude- and phase- comparison techniques) are laid out to provide directional, directional-reactance, nonoffset-resistance and conductance characteristices. The respective signals provided by the voltage transformer and the current transformer are directly used as relaying signals without resorting to any form of mixing. Phase shifts required, are obtained by using magnetic ferrite cores in a novel manner. Sampling units employing a combination of ferrite cores and semiconductor devices provide highly reliable designs. Special attention is paid to the choice of relaying signals, to eliminate the need for any synchronisation or modification and to avoid `image¿ characteristics. These factors have resulted in a considerable simplification of the practical circuitry. A thyristor AND circuit is employed in dual comparator units to provide the final tripping, and leads to a circuit which is much less sensitive to extraneous signals than a single-thyristor unit.

An (O)over-bar(m(2)n) randomized algorithm to compute a minimum cycle basis of a directed graph

We consider the problem of computing a minimum cycle basis in a directed graph G. The input to this problem is a directed graph whose arcs have positive weights. In this problem a {- 1, 0, 1} incidence vector is associated with each cycle and the vector space over Q generated by these vectors is the cycle space of G. A set of cycles is called a cycle basis of G if it forms a basis for its cycle space. A cycle basis where the sum of weights of the cycles is minimum is called a minimum cycle basis of G. The current fastest algorithm for computing a minimum cycle basis in a directed graph with m arcs and n vertices runs in O(m(w+1)n) time (where w < 2.376 is the exponent of matrix multiplication). If one allows randomization, then an (O) over tilde (m(3)n) algorithm is known for this problem. In this paper we present a simple (O) over tilde (m(2)n) randomized algorithm for this problem. The problem of computing a minimum cycle basis in an undirected graph has been well-studied. In this problem a {0, 1} incidence vector is associated with each cycle and the vector space over F-2 generated by these vectors is the cycle space of the graph. The fastest known algorithm for computing a minimum cycle basis in an undirected graph runs in O(m(2)n + mn(2) logn) time and our randomized algorithm for directed graphs almost matches this running time.