Quantum clustering and network analysis of MD simulation trajectories to probe the conformational ensembles of protein-ligand interactions

In this article, we present a novel application of a quantum clustering (QC) technique to objectively cluster the conformations, sampled by molecular dynamics simulations performed on different ligand bound structures of the protein. We further portray each conformational population in terms of dynamically stable network parameters which beautifully capture the ligand induced variations in the ensemble in atomistic detail. The conformational populations thus identified by the QC method and verified by network parameters are evaluated for different ligand bound states of the protein pyrrolysyl-tRNA synthetase (DhPylRS) from D. hafniense. The ligand/environment induced re-distribution of protein conformational ensembles forms the basis for understanding several important biological phenomena such as allostery and enzyme catalysis. The atomistic level characterization of each population in the conformational ensemble in terms of the re-orchestrated networks of amino acids is a challenging problem, especially when the changes are minimal at the backbone level. Here we demonstrate that the QC method is sensitive to such subtle changes and is able to cluster MD snapshots which are similar at the side-chain interaction level. Although we have applied these methods on simulation trajectories of a modest time scale (20 ns each), we emphasize that our methodology provides a general approach towards an objective clustering of large-scale MD simulation data and may be applied to probe multistate equilibria at higher time scales, and to problems related to protein folding for any protein or protein-protein/RNA/DNA complex of interest with a known structure.

Water-mediated ionic interactions in protein structures

It is well known that water molecules play an indispensable role in the structure and function of biological macromolecules. The water-mediated ionic interactions between the charged residues provide stability and plasticity and in turn address the function of the protein structures. Thus, this study specifically addresses the number of possible water-mediated ionic interactions, their occurrence, distribution and nature found in 90% non-redundant protein chains. Further, it provides a statistical report of different charged residue pairs that are mediated by surface or buried water molecules to form the interactions. Also, it discusses its contributions in stabilizing various secondary structural elements of the protein. Thus, the present study shows the ubiquitous nature of the interactions that imparts plasticity and flexibility to a protein molecule.

Statistical thermodynamics of lattice models in zeolites: Implications of local versus global mean field interactions

The statistical thermodynamics of adsorption in caged zeolites is developed by treating the zeolite as an ensemble of M identical cages or subsystems. Within each cage adsorption is assumed to occur onto a lattice of n identical sites. Expressions for the average occupancy per cage are obtained by minimizing the Helmholtz free energy in the canonical ensemble subject to the constraints of constant M and constant number of adsorbates N. Adsorbate-adsorbate interactions in the Brag-Williams or mean field approximation are treated in two ways. The local mean field approximation (LMFA) is based on the local cage occupancy and the global mean field approximation (GMFA) is based on the average coverage of the ensemble. The GMFA is shown to be equivalent in formulation to treating the zeolite as a collection of interacting single site subsystems. In contrast, the treatment in the LMFA retains the description of the zeolite as an ensemble of identical cages, whose thermodynamic properties are conveniently derived in the grand canonical ensemble. For a z coordinated lattice within the zeolite cage, with epsilon(aa) as the adsorbate-adsorbate interaction parameter, the comparisons for different values of epsilon(aa)(*)=epsilon(aa)z/2kT, and number of sites per cage, n, illustrate that for -1 0. We compare the isotherms predicted with the LMFA with previous GMFA predictions [K. G. Ayappa, C. R. Kamala, and T. A. Abinandanan, J. Chem. Phys. 110, 8714 (1999)] (which incorporates both the site volume reduction and a coverage-dependent epsilon(aa)) for xenon and methane in zeolite NaA. In all cases the predicted isotherms are very similar, with the exception of a small steplike feature present in the LMFA for xenon at higher coverages. (C) 1999 American Institute of Physics. [S0021-9606(99)70333-8].

A new class of three dimensional D-pi-A trigonal cryptand derivatives for second-order nonlinear optics

Synthesis, crystal structures, linear and nonlinear optical properties of tris D-pi-A cryptand derivatives with C-3 symmetry are reported. Three fold symmetry inherent in the cryptand molecules has been utilized for designing these molecules. Molecular nonlinearities have been measured by hyper-Rayleigh scattering (HRS) experiments. Among the compounds studied, L-1 adopts non-centrosymmetric crystal structure. Compounds L-1, L-2, L-3 and L-4 show a measurable SHG powder signal. These molecules are more isotropic and have significantly higher melting points than the classical p-nitroaniline based dipolar NLO compounds, making them useful for further device applications. Besides, different acceptor groups can be attached to the cryptand molecules to modulate their NLO properties.

First-principles electronic structure of the delafossites ABO(2) (A=Cu, Ag, Au; B=Al, Ga, Sc, In, Y): evolution of d(10)-d(10) interactions

While bonding between d(10) atoms and ions in molecular systems has been well studied, less attention has been paid to interactions between such seemingly closed shell species in extended inorganic solids. In this contribution, we present visualizations of the electronic structures of the delafossites ABO(2) (A = Cu, Ag, Au) with particular emphasis on the nature of d(10)-d(10) interactions in the close packed plane of the coinage metal ion. We find that on going from Cu to Ag to Au, the extent of bonding between A and A increases. However, the structures (in terms of distances) of these compounds are largely determined by the strongly ionic 13,11 0 interaction and for the larger B ions Sc, In and Y, the A atoms are sufficiently well-separated that A-A bonding is almost negligible. We also analyze some interesting differences between Ag and Au, including the larger A-O covalency of the Au. The trends in electronic structure suggest that the Ag and Au compounds are not good candidate transparent conducting oxides. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

Interactions of interfacial arc cracks

The interaction of two interfacial arc cracks around a circular elastic inclusion embedded in an elastic matrix is examined. New results for stress intensity factors for a pair of interacting cracks are derived for a concentrated force acting in the matrix. For verifying the point load solutions, stress intensity factors under uniform loading are obtained by superposing point force results. For achieving this objective, a general method for generating desired stress fields inside a test region using point loads is described. The energetics of two interacting interfacial arc cracks is discussed in order to shed more light on the debonding of hard or soft inclusions from the matrix. The analysis based on complex variables is developed in a general way to handle the interactions of multiple interfacial arc cracks/straight cracks.

Supramolecular hydrogen-bonded structures in organic amine squarates

Structures of monohydrogen squarates of methylamine, ethylenediamine, 1,3-diaminopropane, 1,4-diaminobutane, 1,5-diaminopentane, N,N'-diemethylpiperazine and N,N,N,N-tetramethylguanidine have been studied in detail. The supramolecular hydrogen-bonded molecular networks are formed by the monoanion of squaric acid by itself or in association with the parent acid. Three types of hydrogen-bonded motifs are observed in these compounds, namely a liner chain, a cyclic dimer and a cyclic tetramer. These hydrogen-bonded motifs formed by the squaric acid species interact with the amine through N-H...O hydrogen-bonding and give rise to predominantly layered structures, while some of them also exhibit three-dimensional structures. Two of the monohydrogen squarate structures also exhibit pi-pi interactions between two squarate rings. The various hydrogen-bonding parameters in the amine squarates are discussed at length. (C) 2002 Elsevier Science B.V. All rights reserved.

Synthesis and characterization of mono- and bis-D-pi-A cryptand derivatives for second-order nonlinear optics and its modulation with different metal ion inputs

Two new classes of mono- and bis-D-pi-A cryptand derivatives with a flexible and a rigid cryptand core have been synthesized. The linear and nonlinear optical properties of these molecules are probed. The three dimensional cavity of the cryptand moiety has been utilized to modulate the SHG intensity to different extents in solution with metal ion inputs such as Ni-II,Cu-II,Zn-II, and Cd-II. We also report that decomplexation events can be used to reversibly modulate their NLO responses.

Twisted aromatics, 9-anthryl and 1-pyrenyl terpyridines organize into novel multi-directional 'ladder-like' motifs in the solid state

9-Anthryl and 1-pyrenyl terpyridines (1 and 2, respectively), key precursors for the design of novel fluorescent sensors have been synthesized and characterized by H-1 NMR, mass spectroscopy and X-ray crystallography. Twisted molecular conformations for each 1 and 2 were observed in their single crystal structures. Energy minimization calculations for the 1 and 2 using the semi-empirical AM1 method show that the 'twisted' conformation is intrinsic to these systems. We observe interconnected networks of edge-to-face CH...pi interactions, which appear to be cooperative in nature, in each of the crystal structures. The two twisted molecules, although having differently shaped polyaromatic hydrocarbon substituents, show similar patterns of edge-to-face CH...pi interactions.The presently described systems comprise of two aromatic surfaces that are almost orthogonal to each other. This twisted or orthogonal nature of the molecules leads to the formation of interesting multi-directional ladder like supramolecular organizations. A combination of edge-to-face and face-to-face packing modes helps to stabilize these motifs. The ladder like architecture in 1 is helical in nature. (C) 2002 Published by Elsevier Science B.V.

Nature and strength of C-H center dot center dot center dot O interactions involving formyl hydrogen atoms: computational and experimental studies of small aldehydes

Structural and electronic properties of C-H center dot center dot center dot O contacts in compounds containing a formyl group are investigated from the perspective of both hydrogen bonding and dipole-dipole interactions, in a systematic and graded approach. The effects of a-substitution and self-association on the nature of the formyl H-atom are studied with the NBO and AIM methodologies. The relative dipole-dipole contributions in formyl C-H center dot center dot center dot O interactions are obtained for aldehyde dimers. The stabilities and energies of aldehyde clusters (dimer through octamer) have been examined computationally. Such studies have an implication in crystallization mechanisms. Experimental X-ray crystal structures of formaldehyde, acrolein and N-methylformamide have been determined in order to ascertain the role of C-H center dot center dot center dot O interactions in the crystal packing of formyl compounds.

Effect of screening of intermicellar interactions on the linear and nonlinear rheology of a viscoelastic gel

We report our studies of the linear and nonlinear rheology of aqueous solutions of the surfactant cetyl trimethylammonium tosylate (CTAT) with varying amounts of sodium chloride (NaCl). The CTAT concentration is fixed at 42 mM, and the salt concentration is varied between 0 and 120 mM. On increasing the salt (NaCl) concentration, we see three distinct regimes in the zero-shear viscosity and the high-frequency plateau modulus data. In regime 1, the zero-shear viscosity shows a weak increase with salt concentration due to enhanced micellar growth. The decrease in the zero-shear viscosities with salt concentration in regimes II and III can be explained in terms of intermicellar branching. The most intriguing feature of our data, however, is the anomalous behavior of the high-frequency plateau modulus in regime II (0.12 less than or equal to [NaCl]/[CTAT] less than or equal to 1.42). In this regime, the plateau modulus increases with an increase in NaCl concentration. This is highly interesting, since the correlation length of concentration fluctuations and hence the plateau modulus G(0) are not expected to change appreciably in the semidilute regime. We propose to explain the changes in regime II in terms of a possible unbinding of the organic counterions (tosylate) from the CTA(+) surfaces on the addition of NaCl. In the nonlinear flow curves of the samples with high salt content, significant deviations from the predictions of the Giesekus model for entangled micelles are observed.

Synthesis and characterization of novel cationic lipid and cholesterol-coated gold nanoparticles and their interactions with dipalmitoylphosphatidylcholine membranes

Novel gold nanoparticles bearing cationic single-chain, double-chain, and cholesterol based amphiphilic units have been synthesized. These nanoparticles represent size-stable entities in which various cationic lipids have been immobilized through their thiol group onto the gold nanoparticle core. The resulting colloids have been characterized by UV-vis, (1)H NMR, FT-IR spectroscopy, and transmission electron microscopy. The average size of the resultant nanoparticles could be controlled by the relative bulkiness of the capping agent. Thus, the average diameters of the nanoparticles formed from the cationic single-chain, double-chain, and cholesterol based thiolate-coated materials were 5.9,2.9, and 2.04 nm, respectively. We also examined the interaction of these cationic gold nanoparticles with vesicular membranes generated from dipalmitoylphosphatidylcholine (DPPC) lipid suspensions. Nanoparticle doped DPPC vesicular suspensions displayed a characteristic surface plasmon band in their UV-vis spectra. Inclusion of nanoparticles in vesicular suspensions led to increases in the aggregate diameters, as evidenced from dynamic light scattering. Differential scanning calorimetric examination indicated that incorporation of single-chain, double-chain, and cholesteryl-linked cationic nanoparticles exert variable effects on the DPPC melting transitions. While increased doping of single-chain nanoparticles in DPPC resulted in the phases that melt at higher temperatures, inclusion of an incremental amount of double-chain nanoparticles caused the lowering of the melting temperature of DPPC. On the other hand, the cationic cholesteryl nanoparticle interacted with DPPC in membranes in a manner somewhat analogous to that of cholesterol itself and caused broadening of the DPPC melting transition.

Weak interactions involving organic fluorine: analysis of structural motifs in Flunazirine and Haloperidol

The crystal structure of Flunazirine, an anticonvulsant drug, is analyzed in terms of intermolecular interactions involving fluorine. The structure displays motifs formed by only weak interactions C–H⋯F and C–H⋯π. The motifs thus generated show cavities, which could serve as hosts for complexation. The structure of Flunazirine displays cavities formed by C–H⋯F and C–H⋯π interactions. Haloperidol, an antipsychotic drug, shows F⋯F interactions in the crystalline lattice in lieu of Cl⋯Cl interactions. However, strong O–H⋯N interactions dominate packing. The salient features of the two structures in terms of intermolecular interactions reveal, even though organic fluorine has lower tendency to engage in hydrogen bonding and F⋯F interactions, these interactions could play a significant role in the design of molecular assemblies via crystal engineering.

Effects of polarization on study of anomalous VVH interactions at a Linear Collider

We investigate the use of beam polarisation as well as nal state polarisation efects in probing the interaction of the Higgs boson with a pair of heavy vector bosons in the process e+e! ffH, where f is any light fermion. The sensitivity of the International Linear Collider (ILC) operating at ps = 500 GeV, to such V V H(V = W=Z) couplings is examined in a model independent way. The efects of ISR and beamstrahlung are discussed.