Charge density distribution and electrostatic interactions of ethionamide: an inhibitor of the enoyl acyl carrier protein reductase (inhA) enzyme of Mycobacterium tuberculosis

An experimental charge density analysis of an anti-TB drug ethionamide was carried out from high resolution X-ray diffraction at 100 K to understand its charge density distribution and electrostatic properties. The experimental results were validated from periodic theoretical charge density calculations performed using CRYSTAL09 at the B3LYP/6-31G** level of theory. The electron density rho(bcp)(r) and the Laplacian of electron density del(2)(rho bcp)(r) of the molecule calculated from both the methods display the charge density distribution of the ethionamide molecule in the crystal field. The electrostatic potential map shows a large electropositive region around the pyridine ring and a large electronegative region at the vicinity of the thiol atom. The calculated experimental dipole moment is 10.6D, which is higher than the value calculated from theory (8.2D). The topological properties of C-H center dot center dot center dot S, N-H center dot center dot center dot N and N-H center dot center dot center dot S hydrogen bonds were calculated, revealing their strength. The charge density analysis of the ethionamide molecule determined from both the experiment and theory gives the topological and electrostatic properties of the molecule, which allows to precisely understand the nature of intra and intermolecular interactions.

Sensitivity of polarization fluctuations to the nature of protein-water interactions: Study of biological water in four different protein-water systems

Since the time of Kirkwood, observed deviations in magnitude of the dielectric constant of aqueous protein solution from that of neat water (similar to 80) and slower decay of polarization have been subjects of enormous interest, controversy, and debate. Most of the common proteins have large permanent dipole moments (often more than 100 D) that can influence structure and dynamics of even distant water molecules, thereby affecting collective polarization fluctuation of the solution, which in turn can significantly alter solution's dielectric constant. Therefore, distance dependence of polarization fluctuation can provide important insight into the nature of biological water. We explore these aspects by studying aqueous solutions of four different proteins of different characteristics and varying sizes, chicken villin headpiece subdomain (HP-36), immunoglobulin binding domain protein G (GB1), hen-egg white lysozyme (LYS), and Myoglobin (MYO). We simulate fairly large systems consisting of single protein molecule and 20000-30000 water molecules (varied according to the protein size), providing a concentration in the range of similar to 2-3 mM. We find that the calculated dielectric constant of the system shows a noticeable increment in all the cases compared to that of neat water. Total dipole moment auto time correlation function of water < dM(W) (0)delta M-W (t) > is found to be sensitive to the nature of the protein. Surprisingly, dipole moment of the protein and total dipole moment of the water molecules are found to be only weakly coupled. Shellwise decomposition of water molecules around protein reveals higher density of first layer compared to the succeeding ones. We also calculate heuristic effective dielectric constant of successive layers and find that the layer adjacent to protein has much lower value (similar to 50). However, progressive layers exhibit successive increment of dielectric constant, finally reaching a value close to that of bulk 4-5 layers away. We also calculate shellwise orientational correlation function and tetrahedral order parameter to understand the local dynamics and structural re-arrangement of water. Theoretical analysis providing simple method for calculation of shellwise local dielectric constant and implication of these findings are elaborately discussed in the present work. (C) 2014 AIP Publishing LLC.

Stacking Interactions in RNA and DNA: Roll-Slide Energy Hyperspace for Ten Unique Dinucleotide Steps

Understanding dinucleotide sequence directed structures of nuleic acids and their variability from experimental observation remained ineffective due to unavailability of statistically meaningful data. We have attempted to understand this from energy scan along twist, roll, and slide degrees of freedom which are mostly dependent on dinucleotide sequence using ab initio density functional theory. We have carried out stacking energy analysis in these dinucleotide parameter phase space for all ten unique dinucleotide steps in DNA and RNA using DFT-D by B97X-D/6-31G(2d,2p), which appears to satisfactorily explain conformational preferences for AU/AU step in our recent study. We show that values of roll, slide, and twist of most of the dinucleotide sequences in crystal structures fall in the low energy region. The minimum energy regions with large twist values are associated with the roll and slide values of B-DNA, whereas, smaller twist values correspond to higher stability to RNA and A-DNA like conformations. Incorporation of solvent effect by CPCM method could explain the preference shown by some sequences to occur in B-DNA or A-DNA conformations. Conformational preference of BII sub-state in B-DNA is preferentially displayed mainly by pyrimidine-purine steps and partly by purine-purine steps. The purine-pyrimidine steps show largest effect of 5-methyl group of thymine in stacking energy and the introduction of solvent reduces this effect significantly. These predicted structures and variabilities can explain the effect of sequence on DNA and RNA functionality. (c) 2014 Wiley Periodicals, Inc. Biopolymers 103: 134-147, 2015.

An overview of recent advances in structural bioinformatics of protein-protein interactions and a guide to their principles

Rich data bearing on the structural and evolutionary principles of protein protein interactions are paving the way to a better understanding of the regulation of function in the cell. This is particularly the case when these interactions are considered in the framework of key pathways. Knowledge of the interactions may provide insights into the mechanisms of crucial `driver' mutations in oncogenesis. They also provide the foundation toward the design of protein protein interfaces and inhibitors that can abrogate their formation or enhance them. The main features to learn from known 3-D structures of protein protein complexes and the extensive literature which analyzes them computationally and experimentally include the interaction details which permit undertaking structure-based drug discovery, the evolution of complexes and their interactions, the consequences of alterations such as post-translational modifications, ligand binding, disease causing mutations, host pathogen interactions, oligomerization, aggregation and the roles of disorder, dynamics, allostery and more to the protein and the cell. This review highlights some of the recent advances in these areas, including design, inhibition and prediction of protein protein complexes. The field is broad, and much work has been carried out in these areas, making it challenging to cover it in its entirety. Much of this is due to the fast increase in the number of molecules whose structures have been determined experimentally and the vast increase in computational power. Here we provide a concise overview. (C) 2014 Elsevier Ltd. All rights reserved.

Heterogeneity in Dye-TiO2 Interactions Dictate Charge Transfer Efficiencies for Diketopyrrolopyrrole-Based Polymer Sensitized Solar Cells

Power conversion efficiency of a solar cell is a complex parameter which usually hides the molecular details of the charge generation process. For rationally tailoring the overall device efficiency of the dye-sensitized solar cell, detailed molecular understanding of photoinduced reactions at the dye-TiO2 interface has to be achieved. Recently, near-IR absorbing diketopyrrolopyrrole-based (DPP) low bandgap polymeric dyes with enhanced photostabilities have been used for TiO2 sensitization with moderate efficiencies. To improve the reported device performances, a critical analysis of the polymerTiO(2) interaction and electron transfer dynamics is imperative. Employing a combination of time-resolved optical measurements complemented by low temperature EPR and steady-state Raman spectroscopy on polymerTiO(2) conjugates, we provide direct evidence for photoinduced electron injection from the TDPP-BBT polymer singlet state into TiO2 through the C-O group of the DPP-core. A detailed excited state description of the electron transfer process in films reveals instrument response function (IRF) limited (<110 fs) charge injection from a minor polymer fraction followed by a picosecond recombination. The major fraction of photoexcited polymers, however, does not show injection indicating pronounced ground state heterogeneity induced due to nonspecific polymerTiO(2) interactions. Our work therefore underscores the importance of gathering molecular-level insight into the competitive pathways of ultrafast charge generation along with probing the chemical heterogeneity at the nanoscale within the polymerTiO2 films for optimizing photovoltaic device efficiencies.

From Workstations to Workbenches: Towards Predicting Physicochemically Viable Protein-Protein Interactions Across a Host and a Pathogen

The understanding of protein-protein interactions is indispensable in comprehending most of the biological processes in a cell. Small-scale experiments as well as large-scale high-throughput techniques over the past few decades have facilitated identification and analysis of protein-protein interactions which form the basis of much of our knowledge on functional and regulatory aspects of proteins. However, such rich catalog of interaction data should be used with caution when establishing protein-protein interactions in silico, as the high-throughput datasets are prone to false positives. Numerous computational means developed to pursue genome-wide studies on protein-protein interactions at times overlook the mechanistic and molecular details, thus questioning the reliability of predicted protein-protein interactions. We review the development, advantages, and shortcomings of varied approaches and demonstrate that by providing a structural viewpoint in terms of shape complementarity and interaction energies at protein-protein interfaces coupled with information on expression and localization of proteins homologous to an interacting pair, it is possible to assess the credibility of predicted interactions in biological context. With a focus on human pathogen Mycobacterium tuberculosis H37Rv, we show that such scrupulous use of details at the molecular level can predict physicochemically viable protein-protein interactions across host and pathogen. Such predicted interactions have the potential to provide molecular basis of probable mechanisms of pathogenesis and hence open up ways to explore their usefulness as targets in the light of drug discovery. (c) 2014 IUBMB Life, 66(11):759-774, 2014

Salts of hexamethylenetetramine with organic acids: Enhanced anomeric interactions with a lowering of molecular symmetry revealed by crystal structures

The hexamethylenetetramine (HMT) framework displays interesting stereoelectronic interactions of the anomeric type. In the highly symmetrical parent system, the nitrogen centres act as both donors and acceptors. Protonation lowers symmetry and also leads to an enhancement of the anomeric interaction around the protonated centre. X-ray diffraction crystal structures of four derivatives of HMT - with succinic, (DL)-malic, phthalic and 4-hydroxybenzoic acids - reveal significant trends. (The first three form well-defined salts, 4-hydroxybenzoic acid forming a co-crystalline compound.) Each molecular structure is essentially characterised by a major anomeric interaction involving the protonated centre as acceptor. In two cases (succinic and 4-hydroxybenzoic), secondary protonation leads to a weaker anomeric interaction site that apparently competes with the dominant one. Bond length changes indicate that the anomeric interaction decreases as malic > phthalic > succinic > 4-hydroxybenzoic, which correlates with the degree of proton transfer to the nitrogen centre. Along with other bond length and angle changes, the results offer insight into the applicability of the antiperiplanar lone pair hypothesis (ALPH) in a rigid system. (C) 2014 Elsevier B.V. All rights reserved.

Salts of hexamethylenetetramine with organic acids: Enhanced anomeric interactions with a lowering of molecular symmetry revealed by crystal structures

The hexamethylenetetramine (HMT) framework displays interesting stereoelectronic interactions of the anomeric type. In the highly symmetrical parent system, the nitrogen centres act as both donors and acceptors. Protonation lowers symmetry and also leads to an enhancement of the anomeric interaction around the protonated centre. X-ray diffraction crystal structures of four derivatives of HMT - with succinic, (DL)-malic, phthalic and 4-hydroxybenzoic acids - reveal significant trends. (The first three form well-defined salts, 4-hydroxybenzoic acid forming a co-crystalline compound.) Each molecular structure is essentially characterised by a major anomeric interaction involving the protonated centre as acceptor. In two cases (succinic and 4-hydroxybenzoic), secondary protonation leads to a weaker anomeric interaction site that apparently competes with the dominant one. Bond length changes indicate that the anomeric interaction decreases as malic > phthalic > succinic > 4-hydroxybenzoic, which correlates with the degree of proton transfer to the nitrogen centre. Along with other bond length and angle changes, the results offer insight into the applicability of the antiperiplanar lone pair hypothesis (ALPH) in a rigid system. (C) 2014 Elsevier B.V. All rights reserved.

Jacalin-carbohydrate interactions: distortion of the ligand molecule as a determinant of affinity

Jacalin is among the most thoroughly studied lectins. Its carbohydrate-binding site has also been well characterized. It has been postulated that the lower affinity of beta-galactosides for jacalin compared with beta-galactosides is caused by steric interactions of the substituents in the former with the protein. This issue has been explored energetically and structurally using different appropriate carbohydrate complexes of jacalin. It turns out that the earlier postulation is not correct. The interactions of the substituent with the binding site remain essentially the same irrespective of the anomeric nature of the substitution. This is achieved through a distortion of the sugar ring in beta-galactosides. The difference in energy, and therefore in affinity, is caused by a distortion of the sugar ring in beta-galactosides. The elucidation of this unprecedented distortion of the ligand as a strategy for modulating affinity is of general interest. The crystal structures also provide a rationale for the relative affinities of the different carbohydrate ligands for jacalin.

DEEP WASHINGTON PHOTOMETRY OF INCONSPICUOUS STAR CLUSTER CANDIDATES IN THE LARGE MAGELLANIC CLOUD

We present deep Washington photometry of 45 poorly populated star cluster candidates in the Large Magellanic Cloud (LMC). We have performed a systematic study to estimate the parameters of the cluster candidates by matching theoretical isochrones to the cleaned and dereddened cluster color-magnitude diagrams. We were able to estimate the basic parameters for 33 clusters, out of which 23 are identified as single clusters and 10 are found to be members of double clusters. The other 12 cluster candidates have been classified as possible clusters/asterisms. About 50% of the true clusters are in the 100-300 Myr age range, whereas some are older or younger. We have discussed the distribution of age, location, and reddening with respect to field, as well as the size of true clusters. The sizes and masses of the studied sample are found to be similar to that of open clusters in the Milky Way. Our study adds to the lower end of cluster mass distribution in the LMC, suggesting that the LMC, apart from hosting rich clusters, also has formed small, less massive open clusters in the 100-300 Myr age range.

Admissible fundamental operators(star)

Let F and G be two bounded operators on two Hilbert spaces. Let their numerical radii be no greater than one. This note investigates when there is a Gamma-contraction (S, P) such that F is the fundamental operator of (S, P) and G is the fundamental operator of (S*, P*). Theorem 1 puts a necessary condition on F and G for them to be the fundamental operators of (S, P) and (S*, P*) respectively. Theorem 2 shows that this necessary condition is also sufficient provided we restrict our attention to a certain special case. The general case is investigated in Theorem 3. Some of the results obtained for Gamma-contractions are then applied to tetrablock contractions to figure out when two pairs (F1, F2) and (G(1), G(2)) acting on two Hilbert spaces can be fundamental operators of a tetrablock contraction (A, B, P) and its adjoint (A*, B*, P*) respectively. This is the content of Theorem 3. (C) 2015 Elsevier Inc. All rights reserved.

Percolation on networks with antagonistic and dependent interactions

Drawing inspiration from real world interacting systems, we study a system consisting of two networks that exhibit antagonistic and dependent interactions. By antagonistic and dependent interactions we mean that a proportion of functional nodes in a network cause failure of nodes in the other, while failure of nodes in the other results in failure of links in the first. In contrast to interdependent networks, which can exhibit first-order phase transitions, we find that the phase transitions in such networks are continuous. Our analysis shows that, compared to an isolated network, the system is more robust against random attacks. Surprisingly, we observe a region in the parameter space where the giant connected components of both networks start oscillating. Furthermore, we find that for Erdos-Renyi and scale-free networks the system oscillates only when the dependence and antagonism between the two networks are very high. We believe that this study can further our understanding of real world interacting systems.

Hierarchical sampling for metastable conformers determines biomolecular recognition: the case of malectin and diglucosylated N-glycan interactions

Structural information over the entire course of binding interactions based on the analyses of energy landscapes is described, which provides a framework to understand the events involved during biomolecular recognition. Conformational dynamics of malectin's exquisite selectivity for diglucosylated N-glycan (Dig-N-glycan), a highly flexible oligosaccharide comprising of numerous dihedral torsion angles, are described as an example. For this purpose, a novel approach based on hierarchical sampling for acquiring metastable molecular conformations constituting low-energy minima for understanding the structural features involved in a biologic recognition is proposed. For this purpose, four variants of principal component analysis were employed recursively in both Cartesian space and dihedral angles space that are characterized by free energy landscapes to select the most stable conformational substates. Subsequently, k-means clustering algorithm was implemented for geometric separation of the major native state to acquire a final ensemble of metastable conformers. A comparison of malectin complexes was then performed to characterize their conformational properties. Analyses of stereochemical metrics and other concerted binding events revealed surface complementarity, cooperative and bidentate hydrogen bonds, water-mediated hydrogen bonds, carbohydrate-aromatic interactions including CH-pi and stacking interactions involved in this recognition. Additionally, a striking structural transition from loop to beta-strands in malectin CRD upon specific binding to Dig-N-glycan is observed. The interplay of the above-mentioned binding events in malectin and Dig-N-glycan supports an extended conformational selection model as the underlying binding mechanism.

Synthesis, X-ray structure and in vitro cytotoxicity studies of Cu(I/II) complexes of thiosemicarbazone: special emphasis on their interactions with DNA

4-(p-X-phenyl)thiosemicarbazone of napthaldehyde {where X = Cl (HL1) and X = Br (HL2)}, thiosemicarbazone of quinoline-2-carbaldehyde (HL3) and 4-(p-fluorophenyl) thiosemicarbazone of salicylaldehyde (H2L4) and their copper(I) {Cu(HL1)(PPh3)(2)Br]center dot CH3CN (1) and Cu(HL2)(PPh3)(2)Cl]center dot DMSO (2)} and copper(II) {((Cu2L2Cl)-Cl-3)(2)(mu-Cl)(2)]center dot 2H(2)O (3) and Cu(L-4)(Py)] (4)} complexes are reported herein. The synthesized ligands and their copper complexes were successfully characterized by elemental analysis, cyclic voltammetry, NMR, ESI-MS, IR and UV-Vis spectroscopy. Molecular structures of all the Cu(I) and Cu(II) complexes have been determined by X-ray crystallography. All the complexes (1-4) were tested for their ability to exhibit DNA-binding and - cleavage activity. The complexes effectively interact with CT-DNA possibly by groove binding mode, with binding constants ranging from 10(4) to 10(5) M-1. Among the complexes, 3 shows the highest chemical (60%) as well as photo-induced (80%) DNA cleavage activity against pUC19 DNA. Finally, the in vitro antiproliferative activity of all the complexes was assayed against the HeLa cell line. Some of the complexes have proved to be as active as the clinical referred drugs, and the greater potency of 3 may be correlated with its aqueous solubility and the presence of the quinonoidal group in the thiosemicarbazone ligand coordinated to the metal.

Hybrid Elements for Modelling Squeeze Film Effects Coupled with Structural Interactions in Vibratory MEMS Devices

We present a hybrid finite element based methodology to solve the coupled fluid structure problem of squeeze film effects in vibratory MEMS devices, such as gyroscopes, RF switches, and 2D resonators. The aforementioned devices often have a thin plate like structure vibrating normally to a fixed substrate, and are generally not perfectly vacuum packed. This results in a thin air film being trapped between the vibrating plate and the fixed substrate which behaves like a squeeze film offering both stiffness and damping. For accurate modelling of such devices the squeeze film effects must be incorporated. Extensive literature is available on squeeze film modelling, however only a few studies address the coupled fluid elasticity problem. The majority of the studies that account for the plate elasticity coupled with the fluid equation, either use approximate mode shapes for the plate or use iterative solution strategies. In an earlier work we presented a single step coupled methodology using only one type of displacement based element to solve the coupled problem. The displacement based finite element models suffer from locking issues when it comes to modelling very thin structures with the lateral dimensions much larger than the plate thickness as is typical in MEMS devices with squeeze film effects. In this work we present another coupled formulation where we have used hybrid elements to model the structural domain. The numerical results show a huge improvement in convergence and accuracy with coarse hybrid mesh as compared to displacement based formulations. We further compare our numerical results with experimental data from literature and find them to be in good accordance.