Dissociation Constants of Weak Acids from ab Initio Molecular Dynamics Using Metadynamics: Influence of the Inductive Effect and Hydrogen Bonding on pK(a) Values

The theoretical estimation of the dissociation constant, or pK(a), of weak acids continues to be a challenging field. Here, we show that ab initio CarParrinello molecular dynamics simulations in conjunction with metadynamics calculations of the free-energy profile of the dissociation reaction provide reasonable estimates of the pK(a) value. Water molecules, sufficient to complete the three hydration shells surrounding the acid molecule, were included explicitly in the computation procedure. The free-energy profiles exhibit two distinct minima corresponding to the dissociated and neutral states of the acid, and the difference in their values provides the estimate for pK(a). We show for a series of organic acids that CPMD simulations in conjunction with metadynamics can provide reasonable estimates of pK(a) values. The acids investigated were aliphatic carboxylic acids, chlorine-substituted carboxylic acids, cis- and trans-butenedioic acid, and the isomers of hydroxybenzoic acid. These systems were chosen to highlight that the procedure could correctly account for the influence of the inductive effect as well as hydrogen bonding on pK(a) values of weak organic acids. In both situations, the CPMD metadynamics procedure faithfully reproduces the experimentally observed trend and the magnitudes of the pK(a) values.

Molecular Dissection of Mycobacterium tuberculosis Integration Host Factor Reveals Novel Insights into the Mode of DNA Binding and Nucleoid Compaction

Background: mIHF belongs to a subfamily of proteins, distinct from E. coli IHF. Results: Functionally important amino acids of mIHF and the mechanism(s) underlying DNA binding, DNA bending, and site-specific recombination are distinct from that of E. coli IHF. Conclusion: mIHF functions could contribute beyond nucleoid compaction. Significance: Because mIHF is essential for growth, the molecular mechanisms identified here can be exploited in drug screening efforts. The annotated whole-genome sequence of Mycobacterium tuberculosis revealed that Rv1388 (Mtihf) is likely to encode for a putative 20-kDa integration host factor (mIHF). However, very little is known about the functional properties of mIHF or the organization of the mycobacterial nucleoid. Molecular modeling of the mIHF three-dimensional structure, based on the cocrystal structure of Streptomyces coelicolor IHF duplex DNA, a bona fide relative of mIHF, revealed the presence of Arg-170, Arg-171, and Arg-173, which might be involved in DNA binding, and a conserved proline (Pro-150) in the tight turn. The phenotypic sensitivity of Escherichia coli ihfA and ihfB strains to UV and methyl methanesulfonate could be complemented with the wild-type Mtihf but not its alleles bearing mutations in the DNA-binding residues. Protein-DNA interaction assays revealed that wild-type mIHF, but not its DNA-binding variants, binds with high affinity to fragments containing attB and attP sites and curved DNA. Strikingly, the functionally important amino acid residues of mIHF and the mechanism(s) underlying its binding to DNA, DNA bending, and site-specific recombination are fundamentally different from that of E. coli IHF. Furthermore, we reveal novel insights into IHF-mediated DNA compaction depending on the placement of its preferred binding sites; mIHF promotes DNA compaction into nucleoid-like or higher order filamentous structures. We therefore propose that mIHF is a distinct member of a subfamily of proteins that serve as essential cofactors in site-specific recombination and nucleoid organization and that these findings represent a significant advance in our understanding of the role(s) of nucleoid-associated proteins.

Sunlight-Induced Covalent Marriage of Two Triply Interlocked Pd-6 Cages and Their Facile Thermal Separation

A template-free triply interlocked Pd-6 cage (2) was synthesized by two-component self-assembly of cis-blocked 90 degrees acceptor cis-(tmen)Pd(NO3)(2) (M) and 1,3,5-tris((E)-2-(pyridin-3-yl)vinyl)benzene (L). Assembly 2 was characterized by H-1 NMR and ESI-MS, and the structure was confirmed by X-ray crystallography, which revealed a parallel conformation of the olefin double bonds belonging to the adjacent cages in the solid state at a distance of 3.656 angstrom, thereby indicating the feasibility of 2+2] photochemical reaction. Two adjacent interlocked cages were covalently married together by intermolecular 2+2] cycloaddition in a single crystal-to-single crystal fashion upon exposure to sunlight/UV irradiation. Most surprisingly, the covalently married pair was easily separated thermally in aqueous medium under mild reaction conditions.

Poly(alkylene itaconate)s - an interesting class of polyesters with periodically located exo-chain double bonds susceptible to Michael addition

Itaconic acid is a bio-sourced dicarboxylic acid that carries a double bond; although several reports have dealt with the radical-initiated chain polymerization of dialkyl itaconates, only a few studies have utilized it as a di-acid monomer to prepare polyesters. In this study, we demonstrate that dibutyl itaconate can be melt-condensed with aliphatic diols to generate unsaturated polyesters; importantly, we show that the double bonds remain unaffected during the melt polymerization. A particularly useful attribute of these polyesters is that the exo-chain double bonds are conjugated to the ester carbonyl and, therefore, can serve as excellent Michael acceptors. A variety of organic thiols, such as alkane thiols, MPEG thiol, thioglycerol, derivatized cysteine etc., were shown to quantitatively Michael-add to the exo-chain double bonds and generate interesting functionalized polyesters. Similarly, organic amines, such as N-methyl-benzylamine, diallyl amine and proline, also add across the double bond; thus, these poly(alkylene itaconate)s could serve as potentially bio-benign polyesters that could be quantitatively transformed into a variety of interesting and potentially useful functionalized polymers.

NrichD database: sequence databases enriched with computationally designed protein-like sequences aid in remote homology detection

NrichD ( ext-link-type=''uri'' xlink:href=''http://proline.biochem.iisc.ernet.in/NRICHD/'' xlink:type=''simple''>http://proline.biochem.iisc.ernet.in/NRICHD/)< /named-content> is a database of computationally designed protein-like sequences, augmented into natural sequence databases that can perform hops in protein sequence space to assist in the detection of remote relationships. Establishing protein relationships in the absence of structural evidence or natural `intermediately related sequences' is a challenging task. Recently, we have demonstrated that the computational design of artificial intermediary sequences/linkers is an effective approach to fill naturally occurring voids in protein sequence space. Through a large-scale assessment we have demonstrated that such sequences can be plugged into commonly employed search databases to improve the performance of routinely used sequence search methods in detecting remote relationships. Since it is anticipated that such data sets will be employed to establish protein relationships, two databases that have already captured these relationships at the structural and functional domain level, namely, the SCOP database and the Pfam database, have been `enriched' with these artificial intermediary sequences. NrichD database currently contains 3 611 010 artificial sequences that have been generated between 27 882 pairs of families from 374 SCOP folds. The data sets are freely available for download. Additional features include the design of artificial sequences between any two protein families of interest to the user.

Al3+-Ion-Triggered Conformational Isomerization of a Rhodamine B Derivative Evidenced by a Fluorescence Signal - A Crystallographic Proof

A newly designed rhodamine B anisaldehyde hydrazone exhibits Al3+-ion-induced cis (L) to trans (L) conformational isomerization with respect to the xanthene moiety through a rotation about a N-N bond; the isomerization is indicated by a detectable naked-eye color change and a turn-on red fluorescence in 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES) buffer (EtOH/Water 1:9 v/v; pH 7.4) at 25 degrees C. In support of this observation, detailed spectroscopic and physicochemical studies along with density function theory (DFT) calculations have been performed. This cis-to-trans conformational isomerization is due to Al3+ ion coordination, which induces this visual color change and the turn-on fluorescence response. To strengthen our knowledge of the conformational isomerization, detailed structural characterizations of the cis and trans isomers in the solid state were performed by single-crystal X-ray diffraction. To the best of our knowledge, this is the first structural report of both cis and trans conformational isomers for this family of compounds. Moreover, this noncytotoxic probe could be used to image the accumulation of Al3+ ions in HeLa and MCF-7 cell lines.

Implication of a sigma-Methane Complex en Route to Elimination of Methane from a Ruthenium Complex: An Experimental and Theoretical Investigation

The five-coordinated 16-electron complex Ru(Me)(dppe)(2)]OTf] (3) undergoes methane elimination at room temperature to afford the ortho-metalated species (dppe){(C6H5)(C6H4)PCH2CH2P(C6H5)(2)}Ru]OTf] (7). Methane elimination, monitored using NMR spectroscopy, revealed no intermediate throughout the reaction. The NOE between Ru-Me protons and ortho phenyl protons and an agostic interaction trans to the methyl group were found in complex 3 by NMR spectroscopy, which form the basis for three plausible pathways for methane elimination and ortho metalation: pathway I (through spatial interaction), pathway II (through oxidative addition and reductive elimination), and pathway III (through agostic interaction). Methane elimination from complex 3 via pathway I was discounted, since it involves interactions through space and not through bonds. Moreover, the calculated energy barrier for the pathway I transition state was quite high (71.3 kcal/mol), which also indicates that this pathway is very unlikely. Furthermore, no spectroscopic evidence for oxidatively added seven-coordinated Ru(IV) species was found and the computed energy barrier of the transition state for pathway II was moderately high (41.1 kcal/mol), which suggests that this cannot be the right pathway for methane elimination and ortho-metalation of complex 3. On the other hand, indirect evidence in the form of chemical reactions point to the most plausible pathway for methane elimination, pathway III, via the intermediacy of a sigma-CH4 complex that could not be found spectroscopically. DFT calculations at several levels on this pathway showed an initial low-barrier rearrangement through TS1 to a square-pyramidal intermediate wherein methyl and agostic C-H are cis to each other. Migration of hydrogen from agostic C-H and elimination of methane proceed through the transition state TS2, which retains a weak metal-H bonding through most parts of the reaction coordinate. Upon comparison of all three pathways, pathway III was found to be the most likely for methane elimination and ortho-metalation of complex 3.

HIF-1 alpha (1772C > T) polymorphism as marker for breast cancer development

Hypoxia-inducible factor 1 alpha (HIF-1 alpha) is an important transcription factor that regulates different cellular responses to hypoxia. HIF-1 alpha is rapidly degraded by von Hippel-Lindau (VHL) protein under normoxic conditions and stabilized under hypoxia. A common variant of HIF-1 alpha (1772C > T) (rs 11549465) polymorphism, corresponding to an amino acid change from proline to serine at 582 position within the oxygen-dependent degradation domain, results in increased stability of the protein and altered transactivation of its target genes. The present study was aimed to find the association between HIF-1 alpha (1772C > T) (rs 11549465) polymorphism and breast cancer development. For this purpose, 348 primary breast cancer patients and 320 healthy and age-matched controls were genotyped through PCR-RFLP method. The genotype frequencies were compared between patients and controls, and their influence on clinical characteristics of breast cancer patients was analyzed. Our study revealed a significant increase of TT genotype in breast cancer patients compared to controls (p = 0.038). Further, TT genotype and T allele were found to be associated with progesterone receptor (PR)-negative status (p < 0.09). None of the clinical variables revealed significant association with HIF-1 alpha (1772C > T) (rs 11549465) polymorphism.

Formal synthesis of degraded sterol (+)-aplykurodinone-1

The formal synthesis of aplykurodinone-1 is accomplished starting from a suitably functionalized bicyclic lactone having the requisite cis-fused ring junction with a quaternary chiral center that was assembled following a Cp2TiCl-mediated radical cyclization protocol. Our synthetic route further elaborates implementation of Grubbs ring closing metathesis (RCM), Eschenmoser-Claisen rearrangement and iodo-lactonization reactions for the synthesis of the final tricyclic precursor of the target molecule. (C) 2015 Elsevier Ltd. All rights reserved.

Molecular Dissection of Mycobacterium tuberculosis Integration Host Factor Reveals Novel Insights into the Mode of DNA Binding and Nucleoid Compaction

The annotated whole-genome sequence of Mycobacterium tuberculosis indicated that Rv1388 (Mtihf) likely encodes a putative 20 kDa integration host factor (mIHF). However, very little is known about the functional properties of mIHF or organization of mycobacterial nucleoid. Molecular modeling of the mIHF three-dimensional structure, based on the cocrystal structure of Streptomyces coelicolor IHF-duplex DNA, a bona fide relative of mIHF, revealed the presence of Arg170, Arg171, and Arg173, which might be involved in DNA binding, and a conserved proline (P150) in the tight turn. The phenotypic sensitivity of Escherichia coli Delta ihfA and Delta ihfB strains to UV and methylmethanesulfonate could be complemented with the wild-type Mtihf, but not its alleles bearing mutations in the DNA-binding residues. Protein DNA interaction assays revealed that wild-type mIHF, but not its DNA-binding variants, bind with high affinity to fragments containing attB and attP sites and curved DNA. Strikingly, the functionally important amino acid residues of mIHF and the mechanism(s) underlying its binding to DNA, DNA bending, and site-specific recombination are fundamentally different from that of E. coli IHF alpha beta. Furthermore, we reveal novel insights into IHF-mediated DNA compaction depending on the placement of its preferred binding sites; mIHF promotes compaction of DNA into nucleoid-like or higher-order filamentous structures. We hence propose that mIHF is a distinct member of a subfamily of proteins that serve as essential cofactors in site-specific recombination and nucleoid organization and that these findings represent a significant advance in our understanding of the role(s) of nucleoid-associated proteins.

Splicing Functions and Global Dependency on Fission Yeast Slu7 Reveal Diversity in Spliceosome Assembly

The multiple short introns in Schizosaccharomyces pombe genes with degenerate cis sequences and atypically positioned polypyrimidine tracts make an interesting model to investigate canonical and alternative roles for conserved splicing factors. Here we report functions and interactions of the S. pombe slu7(+) (spslu7(+)) gene product, known from Saccharomyces cerevisiae and human in vitro reactions to assemble into spliceosomes after the first catalytic reaction and to dictate 3' splice site choice during the second reaction. By using a missense mutant of this essential S. pombe factor, we detected a range of global splicing derangements that were validated in assays for the splicing status of diverse candidate introns. We ascribe widespread, intron-specific SpSlu7 functions and have deduced several features, including the branch nucleotide-to-3' splice site distance, intron length, and the impact of its A/U content at the 5' end on the intron's dependence on SpSlu7. The data imply dynamic substrate-splicing factor relationships in multiintron transcripts. Interestingly, the unexpected early splicing arrest in spslu7-2 revealed a role before catalysis. We detected a salt-stable association with U5 snRNP and observed genetic interactions with spprp1(+), a homolog of human U5-102k factor. These observations together point to an altered recruitment and dependence on SpSlu7, suggesting its role in facilitating transitions that promote catalysis, and highlight the diversity in spliceosome assembly.

De-DUFing the DUFs: Deciphering distant evolutionary relationships of Domains of Unknown Function using sensitive homology detection methods

Background: In the post-genomic era where sequences are being determined at a rapid rate, we are highly reliant on computational methods for their tentative biochemical characterization. The Pfam database currently contains 3,786 families corresponding to ``Domains of Unknown Function'' (DUF) or ``Uncharacterized Protein Family'' (UPF), of which 3,087 families have no reported three-dimensional structure, constituting almost one-fourth of the known protein families in search for both structure and function. Results: We applied a `computational structural genomics' approach using five state-of-the-art remote similarity detection methods to detect the relationship between uncharacterized DUFs and domain families of known structures. The association with a structural domain family could serve as a start point in elucidating the function of a DUF. Amongst these five methods, searches in SCOP-NrichD database have been applied for the first time. Predictions were classified into high, medium and low-confidence based on the consensus of results from various approaches and also annotated with enzyme and Gene ontology terms. 614 uncharacterized DUFs could be associated with a known structural domain, of which high confidence predictions, involving at least four methods, were made for 54 families. These structure-function relationships for the 614 DUF families can be accessed on-line at http://proline.biochem.iisc.ernet.in/RHD_DUFS/. For potential enzymes in this set, we assessed their compatibility with the associated fold and performed detailed structural and functional annotation by examining alignments and extent of conservation of functional residues. Detailed discussion is provided for interesting assignments for DUF3050, DUF1636, DUF1572, DUF2092 and DUF659. Conclusions: This study provides insights into the structure and potential function for nearly 20 % of the DUFs. Use of different computational approaches enables us to reliably recognize distant relationships, especially when they converge to a common assignment because the methods are often complementary. We observe that while pointers to the structural domain can offer the right clues to the function of a protein, recognition of its precise functional role is still `non-trivial' with many DUF domains conserving only some of the critical residues. It is not clear whether these are functional vestiges or instances involving alternate substrates and interacting partners. Reviewers: This article was reviewed by Drs Eugene Koonin, Frank Eisenhaber and Srikrishna Subramanian.

Impact of glycosylation on stability, structure and unfolding of soybean agglutinin (SBA): an insight from thermal perturbation molecular dynamics simulations

Glycosylation has been recognized as one of the most prevalent and complex post-translational modifications of proteins involving numerous enzymes and substrates. Its effect on the protein conformational transitions is not clearly understood yet. In this study, we have examined the effect of glycosylation on protein stability using molecular dynamics simulation of legume lectin soybean agglutinin (SBA). Its glycosylated moiety consists of high mannose type N-linked glycan (Man(9)GlcNAc(2)). To unveil the structural perturbations during thermal unfolding of these two forms, we have studied and compared them to the experimental results. From the perspective of dynamics, our simulations revealed that the nonglycosylated monomeric form is less stable than corresponding glycosylated form at normal and elevated temperatures. Moreover, at elevated temperature thermal destabilization is more prominent in solvent exposed loops, turns and ends of distinct beta sheets. SBA maintains it folded structure due to some important saltbridges, hydrogen bonds and hydrophobic interactions within the protein. The reducing terminal GlcNAc residues interact with the protein residues VAL161, PRO182 and SER225 via hydrophobic and via hydrogen bonding with ASN 9 and ASN 75. Our simulations also revealed that single glycosylation (ASN75) has no significant effect on corresponding cis peptide angle orientation. This atomistic description might have important implications for understanding the functionality and stability of Soybean agglutinin.

A Pd-8 Tetrafacial Molecular Barrel as Carrier for Water Insoluble Fluorophore

A new carbazole-based tetraimidazole ligand 1,3,6,8-tetra(1H-imidazol-1-yl)-9-methyl-9H-carbazole (L) has been synthesized. The unsymmetrical nature of L as well as the rotational freedom of imidazole donor moieties around C-N bond make it a special building unit, which upon treatment with cis-(tmeda)Pd(NO3)(2) produced an unprecedented single linkage-isomeric Pd-8 tetrafacial molecular nanobarrel (PSMBR-1) tmeda N,N,N',N'-tetramethylethane-1,2-diamine]. Unlike closed architectures, open barrel architecture of water-soluble PSMBR-1 makes it an ideal host for some water insoluble polyaromatic hydrocarbons in aqueous medium; one such inclusion complex coroneneCPSMBR-1 was characterized by X-ray diffraction study. Moreover, the potential application of PSMER-1 as carrier in aqueous medium for the transportation of water insoluble fluorophore (perylene) for live cell imaging is explored.

Peptide delta-Turn: Literature Survey and Recent Progress

Among the various types of a-peptide folding motifs, delta-turn, which requires a central cis-amide disposition, has been one of the least extensively investigated. In particular, this main-chain reversal topology has been studied in-depth neither in linear/cyclic peptides nor in proteins. This Minireview article assembles and critically analyzes relevant data from a literature survey on the d-turn conformation in those compounds. Unpublished results from recent conformational energy calculations and a preliminary solution-state analysis on a small model peptide, currently ongoing in our laboratories, are also briefly outlined.