Lysis dynamics and membrane oligomerization pathways for Cytolysin A (ClyA) pore-forming toxin

Pore-forming toxins are known for their ability to efficiently form transmembrane pores which eventually leads to cell lysis. The dynamics of lysis and underlying self-assembly or oligomerization pathways leading to pore formation are incompletely understood. In this manuscript the pore-forming kinetics and lysis dynamics of Cytolysin-A (ClyA) toxins on red blood cells (RBCs) are quantified and compared with experimental lysis data. Lysis experiments are carried out on a fixed mass of RBCs, under isotonic conditions in phosphate-buffered saline, for different initial toxin concentrations ranging from 2.94-14.7 nM. Kinetic models which account for monomer binding, conformation and oligomerization to form the dodecameric ClyA pore complex are developed and lysis is assumed to occur when the number of pores per RBC (n(p)) exceeds a critical number, n(pc). By analysing the model in a sublytic regime (n(p) < n(pc)) the number of pores per RBC to initiate lysis is found to lie between 392 and 768 for the sequential oligomerization mechanism and between 5300 and 6300 for the non-sequential mechanism. Rupture rates which are first order in the number of RBCs are seen to provide the best agreement with the lysis experiments. The time constants for pore formation are estimated to lie between 1 and 20 s and monomer conformation time scales were found to be 2-4 times greater than the oligomerization times. Cell rupture takes places in 100s of seconds, and occurs predominantly with a steady number of pores ranging from 515 to 11 000 on the RBC surface for the sequential mechanism. Both the sequential irreversible and non-sequential kinetics provide similar predictions of the hemoglobin release dynamics, however the hemoglobin released as a function of the toxin concentration was accurately captured only with the sequential model. Each mechanism develops a distinct distribution of mers on the surface, providing a unique experimentally observable fingerprint to identify the underlying oligomerization pathways. Our study offers a method to quantify the extent and dynamics of lysis which is an important aspect of developing novel drug and gene delivery strategies based on pore-forming toxins.

Conserved water-mediated recognition and dynamics of NAD(+) (carboxamide group) to hIMPDH enzyme: water mimic approach toward the design of isoform-selective inhibitor

Inosine monophosphate dehydrogenase (IMPDH) enzyme involves in GMP biosynthesis pathway. Type I hIMPDH is expressed at lower levels in all cells, whereas type II is especially observed in acute myelogenous leukemia, chronic myelogenous leukemia cancer cells, and 10 ns simulation of the IMP-NAD(+) complex structures (PDB ID. 1B3O and 1JCN) have revealed the presence of a few conserved hydrophilic centers near carboxamide group of NAD(+). Three conserved water molecules (W1, W, and W1 `) in di-nucleotide binding pocket of enzyme have played a significant role in the recognition of carboxamide group (of NAD(+)) to D274 and H93 residues. Based on H-bonding interaction of conserved hydrophilic (water molecular) centers within IMP-NAD(+)-enzyme complexes and their recognition to NAD(+), some covalent modification at carboxamide group of di-nucleotide (NAD(+)) has been made by substituting the -CONH(2)group by -CONHNH2 (carboxyl hydrazide group) using water mimic inhibitor design protocol. The modeled structure of modified ligand may, though, be useful for the development of antileukemic agent or it could be act as better inhibitor for hIMPDH-II.

Dynamical facilitation governs glassy dynamics in suspensions of colloidal ellipsoids

One of the greatest challenges in contemporary condensed matter physics is to ascertain whether the formation of glasses from liquids is fundamentally thermodynamic or dynamic in origin. Although the thermodynamic paradigm has dominated theoretical research for decades, the purely kinetic perspective of the dynamical facilitation (DF) theory has attained prominence in recent times. In particular, recent experiments and simulations have highlighted the importance of facilitation using simple model systems composed of spherical particles. However, an overwhelming majority of liquids possess anisotropy in particle shape and interactions, and it is therefore imperative to examine facilitation in complex glass formers. Here, we apply the DF theory to systems with orientational degrees of freedom as well as anisotropic attractive interactions. By analyzing data from experiments on colloidal ellipsoids, we show that facilitation plays a pivotal role in translational as well as orientational relaxation. Furthermore, we demonstrate that the introduction of attractive interactions leads to spatial decoupling of translational and rotational facilitation, which subsequently results in the decoupling of dynamical heterogeneities. Most strikingly, the DF theory can predict the existence of reentrant glass transitions based on the statistics of localized dynamical events, called excitations, whose duration is substantially smaller than the structural relaxation time. Our findings pave the way for systematically testing the DF approach in complex glass formers and also establish the significance of facilitation in governing structural relaxation in supercooled liquids.

A transmission electron microscopy and X-ray photoelectron spectroscopy study of annealing induced gamma-phase nucleation, clustering, and interfacial dynamics in reactively sputtered amorphous alumina thin films

Pure alpha-Al2O3 exhibits a very high degree of thermodynamical stability among all metal oxides and forms an inert oxide scale in a range of structural alloys at high temperatures. We report that amorphous Al2O3 thin films sputter deposited over crystalline Si instead show a surprisingly active interface. On annealing, crystallization begins with nuclei of a phase closely resembling gamma-Alumina forming almost randomly in an amorphous matrix, and with increasing frequency near the substrate/film interface. This nucleation is marked by the signature appearance of sharp (400) and (440) reflections and the formation of a diffuse diffraction halo with an outer maximal radius of approximate to 0.23 nm enveloping the direct beam. The microstructure then evolves by a cluster-coalescence growth mechanism suggestive of swift nucleation and sluggish diffusional kinetics, while locally the Al ions redistribute slowly from chemisorbed and tetrahedral sites to higher anion coordinated sites. Chemical state plots constructed from XPS data and simple calculations of the diffraction patterns from hypothetically distorted lattices suggest that the true origins of the diffuse diffraction halo are probably related to a complex change in the electronic structure spurred by the a-gamma transformation rather than pure structural disorder. Concurrent to crystallization within the film, a substantially thick interfacial reaction zone also builds up at the film/substrate interface with the excess Al acting as a cationic source. (C) 2015 AIP Publishing LLC.

Soft-spring wall based non-periodic boundary conditions for non-equilibrium molecular dynamics of dense fluids

Non-equilibrium molecular dynamics (MD) simulations require imposition of non-periodic boundary conditions (NPBCs) that seamlessly account for the effect of the truncated bulk region on the simulated MD region. Standard implementation of specular boundary conditions in such simulations results in spurious density and force fluctuations near the domain boundary and is therefore inappropriate for coupled atomistic-continuum calculations. In this work, we present a novel NPBC model that relies on boundary atoms attached to a simple cubic lattice with soft springs to account for interactions from particles which would have been present in an untruncated full domain treatment. We show that the proposed model suppresses the unphysical fluctuations in the density to less than 1% of the mean while simultaneously eliminating spurious oscillations in both mean and boundary forces. The model allows for an effective coupling of atomistic and continuum solvers as demonstrated through multiscale simulation of boundary driven singular flow in a cavity. The geometric flexibility of the model enables straightforward extension to nonplanar complex domains without any adverse effects on dynamic properties such as the diffusion coefficient. (c) 2015 AIP Publishing LLC.

Dynamics of the chemical bond: inter- and intra-molecular hydrogen bond

In this discussion, we show that a static definition of a `bond' is not viable by looking at a few examples for both inter-and intra-molecular hydrogen bonding. This follows from our earlier work (Goswami and Arunan, Phys. Chem. Chem. Phys. 2009, 11, 8974) which showed a practical way to differentiate `hydrogen bonding' from `van der Waals interaction'. We report results from ab initio and atoms in molecules theoretical calculations for a series of Rg center dot center dot center dot HX complexes (Rg = He/Ne/Ar and X = F/Cl/Br) and ethane-1,2-diol. Results for the Rg center dot center dot center dot HX/DX complexes show that Rg center dot center dot center dot DX could have a `deuterium bond' even when Rg center dot center dot center dot HX is not `hydrogen bonded', according to the practical criterion given by Goswami and Arunan. Results for ethane-1,2-diol show that an `intra-molecular hydrogen bond' can appear during a normal mode vibration which is dominated by the O center dot center dot center dot O stretching, though a `bond' is not found in the equilibrium structure. This dynamical `bond' formation may nevertheless be important in ensuring the continuity of electron density across a molecule. In the former case, a vibration `breaks' an existing bond and in the later case, a vibration leads to `bond' formation. In both cases, the molecule/complex stays bound irrespective of what happens to this `hydrogen bond'. Both these cases push the borders on the recent IUPAC recommendation on hydrogen bonding (Arunan et al. Pure. Appl. Chem. 2011, 83 1637) and justify the inclusive nature of the definition.

Catalytic pathway, substrate binding and stability in SAICAR synthetase: A structure and molecular dynamics study

The de novo purine biosynthesis is one of the highly conserved pathways among all organisms and is essential for the cell viability. A clear understanding of the enzymes in this pathway would pave way for the development of antimicrobial and anticancer drugs. Phosphoribosylaminoimidazole-succinocar boxamide (SAICAR) synthetase is one of the enzymes in this pathway that catalyzes ATP dependent ligation of carboxyaminoimidazole ribotide (CAIR) with L-aspartate (ASP). Here, we describe eight crystal structures of this enzyme, in C222(1) and H3 space groups, bound to various substrates and substrate mimics from a hyperthermophilic archaea Pyrococcus horikoshii along with molecular dynamics simulations of the structures with substrates. Complexes exhibit minimal deviation from its apo structure. The CAIR binding site displays a preference for pyrimidine nucleotides. In the ADP.TMP-ASP complex, the ASP binds at a position equivalent to that found in Saccharomyces cerevisiae structure (PDB: 2CNU) and thus, clears the ambiguity regarding ASP's position. A possible mode for the inhibition of the enzyme by CTP and UTP, observed earlier in the yeast enzyme, is clearly illustrated in the structures bound to CMP and UMP. The ADP.Mg2+.PO4.CD/MP complex having a phosphate ion between the ATP and CAIR sites strengthens one of the two probable pathways (proposed in Escherichia coli study) of catalytic mechanism and suggests the possibility of a phosphorylation taking place before the ASP's attack on CAIR. Molecular dynamic simulations of this enzyme along with its substrates at 90 degrees C reveal the relative strengths of substrate binding, possible antagonism and the role of Mg2+ ions. (C) 2015 Elsevier Inc. All rights reserved.

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.

Microwave spectroscopic and theoretical investigations of the strongly hydrogen bonded hexafluoroisopropanol...water complex

This paper reports microwave spectroscopic and theoretical investigations on the interaction of water with hexafluoroisopropanol (HFIP). The HFIP monomer can exist in two conformations, antiperiplanar (AP) and synclinical (SC). The former is about 5 kJ mol(-1) more stable than the latter. Theoretical calculations predicted three potential minima for the complex, two having AP and one having SC conformations. Though, the binding energy for the HFIP(SC)...H2O turned out to be larger than that for the other two conformers having HFIP in the AP form, the global minimum for the complex in the potential energy hypersurface had HFIP in the AP form. Experimental rotational constants for four isotopologues measured using a pulsed nozzle Fourier transform microwave spectrometer, correspond to the global minimum in the potential energy hypersurface. The structural parameters and the internal dynamics of the complex could be determined from the rotational spectra of the four isotopologues. The global minimum has the HFIP(AP) as a hydrogen bond donor forming a strong hydrogen bond with H2O. To characterize the strength of the bonding and to probe the other interactions within the complex, atoms in molecules, non-covalent interaction index and natural bond orbital theoretical analyses have been performed.

The SPL7013 dendrimer destabilizes the HIV-1 gp120-CD4 complex

The poly (l-lysine)-based SPL7013 dendrimer with naphthalene disulphonate surface groups blocks the entry of HIV-1 into target cells and is in clinical trials for development as a topical microbicide. Its mechanism of action against R5 HIV-1, the HIV-1 variant implicated in transmission across individuals, remains poorly understood. Using docking and fully atomistic MD simulations, we find that SPL7013 binds tightly to R5 gp120 in the gp120-CD4 complex but weakly to gp120 alone. Further, the binding, although to multiple regions of gp120, does not occlude the CD4 binding site on gp120, suggesting that SPL7013 does not prevent the binding of R5 gp120 to CD4. Using MD simulations to compute binding energies of several docked structures, we find that SPL7013 binding to gp120 significantly weakens the gp120-CD4 complex. Finally, we use steered molecular dynamics (SMD) to study the kinetics of the dissociation of the gp120-CD4 complex in the absence of the dendrimer and with the dendrimer bound in each of the several stable configurations to gp120. We find that SPL7013 significantly lowers the force required to rupture the gp120-CD4 complex and accelerates its dissociation. Taken together, our findings suggest that SPL7013 compromises the stability of the R5 gp120-CD4 complex, potentially preventing the accrual of the requisite number of gp120-CD4 complexes across the virus-cell interface, thereby blocking virus entry.

Docking, molecular dynamics and QM/MM studies to delineate the mode of binding of CucurbitacinE to F-actin

CucurbitacinE (CurE) has been known to bind covalently to F-actin and inhibit depolymerization. However, the mode of binding of CurE to F-actin and the consequent changes in the F-actin dynamics have not been studied. Through quantum mechanical/molecular mechanical (QM/MM) and density function theory (DFT) simulations after the molecular dynamics (MD) simulations of the docked complex of F-actin and CurE, a detailed transition state (TS) model for the Michael reaction is proposed. The TS model shows nucleophilic attack of the sulphur of Cys257 at the beta-carbon of Michael Acceptor of CurE producing an enol intermediate that forms a covalent bond with CurE. The MD results show a clear difference between the structure of the F-actin in free form and F-actin complexed with CurE. CurE affects the conformation of the nucleotide binding pocket increasing the binding affinity between F-actin and ADP, which in turn could affect the nucleotide exchange. CurE binding also limits the correlated displacement of the relatively flexible domain 1 of F-actin causing the protein to retain a flat structure and to transform into a stable ``tense'' state. This structural transition could inhibit depolymerization of F-actin. In conclusion, CurE allosterically modulates ADP and stabilizes F-actin structure, thereby affecting nucleotide exchange and depolymerization of F-actin. (C) 2015 Elsevier Inc. All rights reserved.

Complex patterns arise through spontaneous symmetry breaking in dense homogeneous networks of neural oscillators

There has been much interest in understanding collective dynamics in networks of brain regions due to their role in behavior and cognitive function. Here we show that a simple, homogeneous system of densely connected oscillators, representing the aggregate activity of local brain regions, can exhibit a rich variety of dynamical patterns emerging via spontaneous breaking of permutation or translational symmetries. Upon removing just a few connections, we observe a striking departure from the mean-field limit in terms of the collective dynamics, which implies that the sparsity of these networks may have very important consequences. Our results suggest that the origins of some of the complicated activity patterns seen in the brain may be understood even with simple connection topologies.

Modelling within-host spatiotemporal dynamics of invasive bacterial disease

Mechanistic determinants of bacterial growth, death, and spread within mammalian hosts cannot be fully resolved studying a single bacterial population. They are also currently poorly understood. Here, we report on the application of sophisticated experimental approaches to map spatiotemporal population dynamics of bacteria during an infection. We analyzed heterogeneous traits of simultaneous infections with tagged Salmonella enterica populations (wild-type isogenic tagged strains [WITS]) in wild-type and gene-targeted mice. WITS are phenotypically identical but can be distinguished and enumerated by quantitative PCR, making it possible, using probabilistic models, to estimate bacterial death rate based on the disappearance of strains through time. This multidisciplinary approach allowed us to establish the timing, relative occurrence, and immune control of key infection parameters in a true host-pathogen combination. Our analyses support a model in which shortly after infection, concomitant death and rapid bacterial replication lead to the establishment of independent bacterial subpopulations in different organs, a process controlled by host antimicrobial mechanisms. Later, decreased microbial mortality leads to an exponential increase in the number of bacteria that spread locally, with subsequent mixing of bacteria between organs via bacteraemia and further stochastic selection. This approach provides us with an unprecedented outlook on the pathogenesis of S. enterica infections, illustrating the complex spatial and stochastic effects that drive an infectious disease. The application of the novel method that we present in appropriate and diverse host-pathogen combinations, together with modelling of the data that result, will facilitate a comprehensive view of the spatial and stochastic nature of within-host dynamics. © 2008 Grant et al.

Forced Dissociation of Selectin-ligand Complexes Using Steered Molecular Dynamics Simulation

Selectin-ligand interactions are crucial to such biological processes as inflammatory cascade or tumor metastasis. How transient formation and dissociation of selectin-ligand bonds in blood flow are coupled to molecular conformation at atomic level, however, has not been well understood. In this study, steered molecular dynamics (SMD) simulations were used to elucidate the intramolecular and intermolecular conformational evolutions involved in forced dissociation of three selectin-ligand systems: the construct consisting of P-selectin lectin (Lec) and epidermal growth factor (EGF)-like domains (P-LE) interacting with synthesized sulfoglycopeptide or SGP-3, P-LE with sialyl Lewis X (sLeX), and E-LE with sLeX. SMD simulations were based on newly built-up force field parameters including carbohydrate units and sulfated tyrosine(s) using an analogy approach. The simulations demonstrated that the complex dissociation was coupled to the molecular extension. While the intramolecular unraveling in P-LESGP-3 system mainly resulted from the destroy of the two anti-parallel sheets of EGF domain and the breakage of hydrogen-bond cluster at the Lec-EGF interface, the intermolecular dissociation was mainly determined by separation of fucose (FUC) from Ca2+ ion in all three systems. Conformational changes during forced dissociations depended on pulling velocities and forces, as well as on how the force was applied. This work provides an insight into better understanding of conformational changes and adhesive functionality of selectin-ligand interactions under external forces.

A complex systems modelling approach for decentralized simulation of electrical microgrids

ICECCS 2010