Calculation of Chemical Potentials and Occupancies in Clathrate Hydrates through Monte Carlo Molecular Simulations

The flexibility of the water lattice in clathrate hydrates and guest-guest interactions has been shown in previous studies to significantly affect the values of the thermodynamic properties, such as chemical potentials and free energies. Here we describe methods for computing occupancies, chemical potentials, and free energies that account for the flexibility of water lattice and guest-guest interactions in the hydrate phase. The methods are validated for a wide variety of guest molecules, such as methane, ethane, carbon dioxide, and tetrahydrodfuran by comparing the predicted occupancy values of guest molecules with those obtained from isothermal isobaric semigrand Monte Carlo simulations. The proposed methods extend the van der Waals and Platteuw theory for clathrate hydrates, and the Langmuir constant is calculated based on the structure of the empty hydrate lattice. These methods in combination with development of advanced molecular models for water and guest molecules should lead to a more thermodynamically consistent theory for clathrate hydrates.

Advances in the molecular design of potential anticancer agents via targeting of human telomeric DNA

Telomerases are an attractive drug target to develop new generation drugs against cancer. A telomere appears from the chromosomal termini and protects it from double-stranded DNA degradation. A short telomere promotes genomic instability, like end-to-end fusion and regulates the over-expression of the telomere repairing enzyme, telomerase. The telomerase maintains the telomere length, which may lead to genetically abnormal situations, leading to cancer. Thus, the design and synthesis of an efficient telomerase inhibitor is a viable strategy toward anticancer drugs development. Accordingly, small molecule induced stabilization of the G-quadruplex structure, formed by the human telomeric DNA, is an area of contemporary scientific art. Several such compounds efficiently stabilize the G-quadruplex forms of nucleic acids, which often leads to telomerase inhibition. This Feature article presents the discovery and development of the telomere structure, function and evolution in telomere targeted anticancer drug design and incorporates the recent advances in this area, in addition to discussing the advantages and disadvantages in the methods, and prospects for the future.

Calculation of three-phase methane-ethane binary clathrate hydrate phase equilibrium from Monte Carlo molecular simulations

Methane and ethane are the simplest hydrocarbon molecules that can form clathrate hydrates. Previous studies have reported methods for calculating the three-phase equilibrium using Monte Carlo simulation methods in systems with a single component in the gas phase. Here we extend those methods to a binary gas mixture of methane and ethane. Methane-ethane system is an interesting one in that the pure components form sII clathrate hydrate whereas a binary mixture of the two can form the sII clathrate. The phase equilibria computed from Monte Carlo simulations show a good agreement with experimental data and are also able to predict the sI-sII structural transition in the clathrate hydrate. This is attributed to the quality of the TIP4P/Ice and TRaPPE models used in the simulations. (C) 2014 Elsevier B.V. All rights reserved.

Fine-tuning solid-state luminescence in NPIs (1,8-naphthalimides): impact of the molecular environment and cumulative interactions

An investigation of a series of seven angular ``V'' shaped NPIs (1-7) is presented. The effect of substitution of these structurally similar NPIs on their photophysical properties in the solution-state and the solid-state is presented and discussed in light of experimental and computational findings. Compounds 1-7 show negligible to intensely strong emission yields in their solid-state depending on the nature of substituents appended to the oxoaryl moiety. The solution and solid-state properties of the compounds can be directly correlated with their structural rigidity, nature of substituents and intermolecular interactions. The versatile solid-state structures of the NPI siblings are deeply affected by the pendant substituents. All of the NPIs (1-7) show antiparallel dimeric pi-pi stacking interactions in their solid-state which can further extend in a parallel, alternate, orthogonal or lateral fashion depending on the steric and electronic nature of the C-4' substituents. Structural investigations including Hirshfeld surface analysis methods reveal that where strongly interacting systems show weak to moderate emission in their condensed states, weakly interacting systems show strong emission yields under the same conditions. The nature of packing and extended structures also affects the emission colors of the NPIs in their solid-states. Furthermore, DFT computational studies were utilized to understand the molecular and cumulative electronic behaviors of the NPIs. The comprehensive studies provide insight into the condensed-state luminescence of aggregationprone small molecules like NPIs and help to correlate the structure-property relationships.

Quorum Sensing and Biofilm Formation in Mycobacteria: Role of c-di-GMP and Methods to Study This Second Messenger

Bacteria have evolved to survive the ever-changing environment using intriguing mechanisms of quorum sensing (QS). Very often, QS facilitates formation of biofilm to help bacteria to persist longer and the formation of such biofilms is regulated by c-di-GMP. It is a well-known second messenger also found in mycobacteria. Several methods have been developed to study c-di-GMP signaling pathways in a variety of bacteria. In this review, we have attempted to highlight a connection between c-di-GMP and biofilm formation and QS in mycobacteria and several methods that have helped in better understanding of c-di-GMP signaling. (c) 2014 IUBMB Life, 66(12):823-834, 2014

CHEXVIS: a tool for molecular channel extraction and visualization

Background: Understanding channel structures that lead to active sites or traverse the molecule is important in the study of molecular functions such as ion, ligand, and small molecule transport. Efficient methods for extracting, storing, and analyzing protein channels are required to support such studies. Further, there is a need for an integrated framework that supports computation of the channels, interactive exploration of their structure, and detailed visual analysis of their properties. Results: We describe a method for molecular channel extraction based on the alpha complex representation. The method computes geometrically feasible channels, stores both the volume occupied by the channel and its centerline in a unified representation, and reports significant channels. The representation also supports efficient computation of channel profiles that help understand channel properties. We describe methods for effective visualization of the channels and their profiles. These methods and the visual analysis framework are implemented in a software tool, CHEXVIS. We apply the method on a number of known channel containing proteins to extract pore features. Results from these experiments on several proteins show that CHEXVIS performance is comparable to, and in some cases, better than existing channel extraction techniques. Using several case studies, we demonstrate how CHEXVIS can be used to study channels, extract their properties and gain insights into molecular function. Conclusion: CHEXVIS supports the visual exploration of multiple channels together with their geometric and physico-chemical properties thereby enabling the understanding of the basic biology of transport through protein channels. The CHEXVIS web-server is freely available at http://vgl.serc.iisc.ernet.in/chexvis/. The web-server is supported on all modern browsers with latest Java plug-in.

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.

MS3ALIGN: an efficient molecular surface aligner using the topology of surface curvature

Background: Aligning similar molecular structures is an important step in the process of bio-molecular structure and function analysis. Molecular surfaces are simple representations of molecular structure that are easily constructed from various forms of molecular data such as 3D atomic coordinates (PDB) and Electron Microscopy (EM) data. Methods: We present a Multi-Scale Morse-Smale Molecular-Surface Alignment tool, MS3ALIGN, which aligns molecular surfaces based on significant protrusions on the molecular surface. The input is a pair of molecular surfaces represented as triangle meshes. A key advantage of MS3ALIGN is computational efficiency that is achieved because it processes only a few carefully chosen protrusions on the molecular surface. Furthermore, the alignments are partial in nature and therefore allows for inexact surfaces to be aligned. Results: The method is evaluated in four settings. First, we establish performance using known alignments with varying overlap and noise values. Second, we compare the method with SurfComp, an existing surface alignment method. We show that we are able to determine alignments reported by SurfComp, as well as report relevant alignments not found by SurfComp. Third, we validate the ability of MS3ALIGN to determine alignments in the case of structurally dissimilar binding sites. Fourth, we demonstrate the ability of MS3ALIGN to align iso-surfaces derived from cryo-electron microscopy scans. Conclusions: We have presented an algorithm that aligns Molecular Surfaces based on the topology of surface curvature. Awebserver and standalone software implementation of the algorithm available at http://vgl.serc.iisc.ernet. in/ms3align.