An automated approach to network features of protein structure ensembles

Network theory applied to protein structures provides insights into numerous problems of biological relevance. The explosion in structural data available from PDB and simulations establishes a need to introduce a standalone-efficient program that assembles network concepts/parameters under one hood in an automated manner. Herein, we discuss the development/application of an exhaustive, user-friendly, standalone program package named PSN-Ensemble, which can handle structural ensembles generated through molecular dynamics (MD) simulation/NMR studies or from multiple X-ray structures. The novelty in network construction lies in the explicit consideration of side-chain interactions among amino acids. The program evaluates network parameters dealing with topological organization and long-range allosteric communication. The introduction of a flexible weighing scheme in terms of residue pairwise cross-correlation/interaction energy in PSN-Ensemble brings in dynamical/chemical knowledge into the network representation. Also, the results are mapped on a graphical display of the structure, allowing an easy access of network analysis to a general biological community. The potential of PSN-Ensemble toward examining structural ensemble is exemplified using MD trajectories of an ubiquitin-conjugating enzyme (UbcH5b). Furthermore, insights derived from network parameters evaluated using PSN-Ensemble for single-static structures of active/inactive states of 2-adrenergic receptor and the ternary tRNA complexes of tyrosyl tRNA synthetases (from organisms across kingdoms) are discussed. PSN-Ensemble is freely available from http://vishgraph.mbu.iisc.ernet.in/PSN-Ensemble/psn_index.html.

Detecting symmetry in scalar fields using augmented extremum graphs

Visualizing symmetric patterns in the data often helps the domain scientists make important observations and gain insights about the underlying experiment. Detecting symmetry in scalar fields is a nascent area of research and existing methods that detect symmetry are either not robust in the presence of noise or computationally costly. We propose a data structure called the augmented extremum graph and use it to design a novel symmetry detection method based on robust estimation of distances. The augmented extremum graph captures both topological and geometric information of the scalar field and enables robust and computationally efficient detection of symmetry. We apply the proposed method to detect symmetries in cryo-electron microscopy datasets and the experiments demonstrate that the algorithm is capable of detecting symmetry even in the presence of significant noise. We describe novel applications that use the detected symmetry to enhance visualization of scalar field data and facilitate their exploration.

Versatility of azide in serendipitous assembly of Copper(II) Magnetic Polyclusters

Engineering at the molecular level is one of the most exciting new developments for the generation of functional materials. However, the concept of designing polynuclear extended structures from bottom up is still not mature. Although progress has been made with secondary building units (SBUs) in metal organic frameworks (MOFs), the control seems to be just an illusion when it comes to bridging ligands such as the azide ion. When we say that the azido ligand is versatile in its bridging capabilities, what we mean is that it would be difficult to predict or control its bridging properties. However, this kind of serendipity is not always bad news. For example, scientists have shown that the azido ligand can mediate magnetic exchanges between paramagnetic metals in a predictable fashion (usually depending upon the bonding geometries). Therefore, it is a well-respected ligand in polynuclear assemblies. Serendipitous assemblies offer new magnetic structures that we may not otherwise even think about synthesizing. The azido ligand forms a variety of complexes with copper(II) using different blocking amines or pyridine based ligands. Its structural nature changes upon changing the substitution on amine, as well as the amount of blocking ligand. In principle, if we take any of these complexes and provide more coordination sites to the bridging azido ligands by removing a fraction of the blocking ligands, we can get new complexes with intricate structural networks and therefore different magnetic properties with the same components as used for the parent complex. In this Account, we mainly discuss the development of a number of new topological and magnetic exchange systems synthesized using this concept. Not all of these new complexes can be grouped according to their basic building structures or even by the ratio of the metal to blocking ligand. Therefore, we divided the discussion by the nuclearity of the basic building structures. Some of the complexes with the same nuclearities have very similar or even almost identical basic structures. However, the way these building units are joined together (by the azido bridges) to form the overall extended structures differ almost in every case. The complexes having the Cu-6 core are particularly interesting from a structural point of view. Although they have almost identical basic structures, some of them are extended in three dimensions, but two of them are extended in two dimensions by two different bridging networks. In the complexes having linear Cu-4 basic units, we find that using similar ligands does not always give the same bridging networks even within the basic building structures. These complexes have also enriched the field of molecular magnetism. One of the complexes with a Cu-3 building unit has provided us with the opportunity to study the competing behavior of two different kinds of magnetic exchange mechanism (ferromagnetic and antiferromagnetic) acting simultaneously between two metal ions. Through density functional theory calculations, we showed how they work independently and their additive nature to produce the overall effect. The exciting methodology for the generation of copper(II) polyclusters presented in this Account will provide the opportunity to explore analogous serendipitous assembly of diverse structures with interesting magnetic behavior using other transition metal ions having more than one unpaired electrons.

Constraining complex aquifer geometry with geophysics (2-D ERT and MRS measurements) for stochastic modelling of groundwater flow

Stochastic modelling is a useful way of simulating complex hard-rock aquifers as hydrological properties (permeability, porosity etc.) can be described using random variables with known statistics. However, very few studies have assessed the influence of topological uncertainty (i.e. the variability of thickness of conductive zones in the aquifer), probably because it is not easy to retrieve accurate statistics of the aquifer geometry, especially in hard rock context. In this paper, we assessed the potential of using geophysical surveys to describe the geometry of a hard rock-aquifer in a stochastic modelling framework. The study site was a small experimental watershed in South India, where the aquifer consisted of a clayey to loamy-sandy zone (regolith) underlain by a conductive fissured rock layer (protolith) and the unweathered gneiss (bedrock) at the bottom. The spatial variability of the thickness of the regolith and fissured layers was estimated by electrical resistivity tomography (ERT) profiles, which were performed along a few cross sections in the watershed. For stochastic analysis using Monte Carlo simulation, the generated random layer thickness was made conditional to the available data from the geophysics. In order to simulate steady state flow in the irregular domain with variable geometry, we used an isoparametric finite element method to discretize the flow equation over an unstructured grid with irregular hexahedral elements. The results indicated that the spatial variability of the layer thickness had a significant effect on reducing the simulated effective steady seepage flux and that using the conditional simulations reduced the uncertainty of the simulated seepage flux. As a conclusion, combining information on the aquifer geometry obtained from geophysical surveys with stochastic modelling is a promising methodology to improve the simulation of groundwater flow in complex hard-rock aquifers. (C) 2013 Elsevier B.V. All rights reserved.

Low formation energy and kinetic barrier of Stone-Wales defect in infinite and finite silicene

Stone-Wales (SW) defects in materials having hexagonal lattice are the most common topological defects that affect the electronic and mechanical properties. Using first principles density functional theory based calculations, we study the formation energy and kinetic barrier of SW-defect in infinite and finite sheets of silicene. The formation energies as well as the barriers in both the cases are significantly lower than those of graphene. Furthermore, compared with the infinite sheets, the energy barriers and formation energies are lower for finite sheets. However, due to low barriers these defects are expected to heal out of the finite sheets. (C) 2013 Elsevier B.V. All rights reserved.

Combinatorial Exploration of the Structural Landscape of Acid-Pyridine Cocrystals

The structural landscape of acid-pyridine cocrystals is explored by adopting a combinatorial matrix method with 4-substituted benzoic acids and 4-substituted pyridines. The choice of the system restricts the primary synthon to the robust acid-pyridine entity. This methodology accordingly provides hints toward the formation of secondary synthons. The pK(a) rule is validated in the landscape by taking all components of the matrix together and exploring it as a whole. Along with the global features, the exploration of landscapes reveals some local features. Apart from the identification of secondary synthons, it also sheds light on the propensity of hydration in cocrystals, synthon competition, and certain topological similarities. The method described here combines two approaches, namely, database analysis and high throughput crystallography, to extract more information with minimal extra experimental effort.

Enhanced Gas Adsorption on Graphitic Substrates via Defects and Local Curvature: A Density Functional Theory Study

Using van-der-Waals-corrected density functional theory calculations, we explore the possibility of engineering the local structure and morphology of high-surface-area graphene-derived materials to improve the uptake of methane and carbon dioxide for gas storage and sensing. We test the sensitivity of the gas adsorption energy to the introduction of native point defects, curvature, and the application of strain. The binding energy at topological point defect sites is inversely correlated with the number of missing carbon atoms, causing Stone-Wales defects to show the largest enhancement with respect to pristine graphene (similar to 20%). Improvements of similar magnitude are observed at concavely curved surfaces in buckled graphene sheets under compressive strain, whereas tensile strain tends to weaken gas binding. Trends for CO2 and CH4 are, similar, although CO2 binding is generally stronger by similar to 4 to 5 kJ mol(-1). However, the differential between the adsorption of CO2 and CH4 is much higher on folded graphene sheets and at concave curvatures; this could possibly be leveraged for CH4/CO2 flow separation and gasselective sensors.

Precision control of drying using rhythmic dancing of sessile nanoparticle laden droplets

This work analyses the unique spatio-temporal alteration of the deposition pattern of evaporating nanoparticle laden droplets resting on a hydrophobic surface through targeted low frequency substrate vibrations. External excitation near the lowest resonant mode (n = 2) of the droplet initially de-pins and then subsequently re-pins the droplet edge creating pseudo-hydrophilicity (low contact angle). Vibration subsequently induces droplet shape oscillations (cyclic elongation and flattening) resulting in strong flow recirculation. This strong radially outward liquid flow augments nanoparticle transport, vaporization, and agglomeration near the pinned edge resulting in much reduced drying time under certain characteristic frequency of oscillations. The resultant deposit exhibits a much flatter structure with sharp, defined peripheral wedge topology as compared to natural drying. Such controlled manipulation of transport enables tailoring of structural and topological morphology of the deposits and offers possible routes towards controlling the formation and drying timescales which are crucial for applications ranging from pharmaceutics to surface patterning. (C) 2014 AIP Publishing LLC.

Two new atom centered fragment descriptors and scoring function enhance classification of antibacterial activity

Classification of pharmacologic activity of a chemical compound is an essential step in any drug discovery process. We develop two new atom-centered fragment descriptors (vertex indices) - one based solely on topological considerations without discriminating atomor bond types, and another based on topological and electronic features. We also assess their usefulness by devising a method to rank and classify molecules with regard to their antibacterial activity. Classification performances of our method are found to be superior compared to two previous studies on large heterogeneous data sets for hit finding and hit-to-lead studies even though we use much fewer parameters. It is found that for hit finding studies topological features (simple graph) alone provide significant discriminating power, and for hit-to-lead process small but consistent improvement can be made by additionally including electronic features (colored graph). Our approach is simple, interpretable, and suitable for design of molecules as we do not use any physicochemical properties. The singular use of vertex index as descriptor, novel range based feature extraction, and rigorous statistical validation are the key elements of this study.

Intermolecular interactions, charge-density distribution and the electrostatic properties of pyrazinamide anti-TB drug molecule: an experimental and theoretical charge-density study

An experimental charge-density analysis of pyrazinamide (a first line antitubercular drug) was performed using high-resolution X-ray diffraction data (sin theta/lambda)(max) = 1.1 angstrom(-1)] measured at 100 (2) K. The structure was solved by direct methods using SHELXS97 and refined by SHELXL97. The total electron density of the pyrazinamide molecule was modeled using the Hansen-Coppens multipole formalism implemented in the XD software. The topological properties of electron density determined from the experiment were compared with the theoretical results obtained from CRYSTAL09 at the B3LYP/6-31G** level of theory. The crystal structure was stabilized by N-H center dot center dot center dot N and N-H center dot center dot center dot O hydrogen bonds, in which the N3-H3B center dot center dot center dot N1 and N3-H3A center dot center dot center dot O1 interactions form two types of dimers in the crystal. Hirshfeld surface analysis was carried out to analyze the intermolecular interactions. The fingerprint plot reveals that the N center dot center dot center dot H and O center dot center dot center dot H hydrogen-bonding interactions contribute 26.1 and 18.4%, respectively, of the total Hirshfeld surface. The lattice energy of the molecule was calculated using density functional theory (B3LYP) methods with the 6-31G** basis set. The molecular electrostatic potential of the pyrazinamide molecule exhibits extended electronegative regions around O1, N1 and N2. The existence of a negative electrostatic potential (ESP) region just above the upper and lower surfaces of the pyrazine ring confirm the pi-electron cloud.

Shape Memory Alloy based Active Chiral Composite Cells

Wing morphing is one of the emerging methodology towards improving aerodynamic efficiency of flight vehicle structures. In this paper a morphing structural element is designed and studied which has its origin in the well known chiral structures. The new aspect of design and functionality explored in this paper is that the chiral cell is actuated using thermal Shape Memory Alloy (SMA) actuator wires to provide directional motion. Such structure utilizes the potential of different actuations concepts based on actuator embedded in the chiral structure skin. This paper describes a new class of chiral cell structure with integrated SMA wire for actuation. Chiral topological constructs are obtained by considering passive and active load path decoupling and sub-optimal shape changes. Single cell of chiral honeycomb with actuators are analyzed using finite element simulation results and experiments. To this end, a multi-cell plan-form is characterized showing interesting possibilities in structural morphing applications. The applicability of the developed chiral cell to flexible wing skin, variable stiffness based design and controlling longitudinal-to-transverse stiffness ratio are discussed.

Evaluating FuncSION: A software for automated synthesis of design solutions for stimulating ideation during mechanical conceptual design

The goal of the work reported in this paper is to use automated, combinatorial synthesis to generate alternative solutions to be used as stimuli by designers for ideation. FuncSION, a computational synthesis tool that can automatically synthesize solution concepts for mechanical devices by combining building blocks from a library, is used for this purpose. The objectives of FuncSION are to help generate a variety of functional requirements for a given problem and a variety of concepts to fulfill these functions. A distinctive feature of FuncSION is its focus on automated generation of spatial configurations, an aspect rarely addressed by other computational synthesis programs. This paper provides an overview of FuncSION in terms of representation of design problems, representation of building blocks, and rules with which building blocks are combined to generate concepts at three levels of abstraction: topological, spatial, and physical. The paper then provides a detailed account of evaluating FuncSION for its effectiveness in providing stimuli for enhanced ideation.

Molecular topology of three ring nematogens from C-13-H-1 dipolar couplings

Molecules in their liquid crystalline phase undergo rotational motion about the long axis of the molecule and the shape adopted by the rotating molecule plays an important role in influencing the mesophase morphology. In this context, obtaining the topology and the relative orientation of the different sub-units are important steps. For studying the liquid crystalline phase, C-13 NMR spectroscopy is a convenient method and for certain specifically designed nematogens, 2-dimensional separated local field (2D-SLF) NMR spectroscopy provides a particularly simple and straightforward means of arriving at the molecular topology. We demonstrate this approach on two three ring based nematogens designed with a phenyl or a thiophene ring at one of the termini. From the C-13-H-1 dipolar couplings of the terminal carbon obtained using the 2D-SLF NMR technique, the order parameter of the local symmetry axis of the terminal phenyl ring as well as of the long molecular axis could be easily estimated. For the thiophene nematogen, the lack of symmetry of the thiophene moiety necessitates some additional computational steps. The results indicate that the thiophene unit has its local ordering axis oriented away from the long molecular axis by a small angle, consistent with a bent structure expected in view of the thiophene geometry. The experiment also demonstrates the ability of 2D-SLF NMR to provide high resolution spectra by separation of several overlapped resonances in terms of their C-13-H-1 dipolar couplings. The results are consistent with a rod-like topology of the core of the investigated mesogens. The investigation demonstrates the potential of 2D-SLF NMR C-13 spectroscopy for obtaining atomistic level information and its utility for topological studies of different mesogens.

On entanglement entropy functionals in higher-derivative gravity theories

In arXiv:1310.5713 1] and arXiv:1310.6659 2] a formula was proposed as the entanglement entropy functional for a general higher-derivative theory of gravity, whose lagrangian consists of terms containing contractions of the Riemann tensor. In this paper, we carry out some tests of this proposal. First, we find the surface equation of motion for general four-derivative gravity theory by minimizing the holographic entanglement entropy functional resulting from this proposed formula. Then we calculate the surface equation for the same theory using the generalized gravitational entropy method of arXiv:1304.4926 3]. We find that the two do not match in their entirety. We also construct the holographic entropy functional for quasi-topological gravity, which is a six-derivative gravity theory. We find that this functional gives the correct universal terms. However, as in the R-2 case, the generalized gravitational entropy method applied to this theory does not give exactly the surface equation of motion coming from minimizing the entropy functional.

On the Utility of Microscopy and C-Scan Techniques in the Study of Defects Arising From Resin Infusion Process and a Study on the Influence of Voids on the Impact Behavior of CFRP Composites

Carbon Fiber Reinforced Plastic composites were fabricated through vacuum resin infusion technology by adopting two different processing conditions, viz., vacuum only in the first and vacuum plus external pressure in the next, in order to generate two levels of void-bearing samples. They were relatively graded as higher and lower void-bearing ones, respectively. Microscopy and C-scan techniques were utilized to describe the presence of voids arising from the two different processing parameters. Further, to determine the influence of voids on impact behavior, the fabricated +45 degrees/90 degrees/-45 degrees composite samples were subjected to low velocity impacts. The tests show impact properties like peak load and energy to peak load registering higher values for the lower void-bearing case where as the total energy, energy for propagation and ductility indexes were higher for the higher void-bearing ones. Fractographic analysis showed that higher void-bearing samples display lower number of separation of layers in the laminate. These and other results are described and discussed in this report.