Spiral-Wave Dynamics in a Mathematical Model of Human Ventricular Tissue with Myocytes and Fibroblasts

Cardiac fibroblasts, when coupled functionally with myocytes, can modulate the electrophysiological properties of cardiac tissue. We present systematic numerical studies of such modulation of electrophysiological properties in mathematical models for (a) single myocyte-fibroblast (MF) units and (b) two-dimensional (2D) arrays of such units; our models build on earlier ones and allow for zero-, one-, and two-sided MF couplings. Our studies of MF units elucidate the dependence of the action-potential (AP) morphology on parameters such as E-f, the fibroblast resting-membrane potential, the fibroblast conductance G(f), and the MF gap-junctional coupling G(gap). Furthermore, we find that our MF composite can show autorhythmic and oscillatory behaviors in addition to an excitable response. Our 2D studies use (a) both homogeneous and inhomogeneous distributions of fibroblasts, (b) various ranges for parameters such as G(gap), G(f), and E-f, and (c) intercellular couplings that can be zero-sided, one-sided, and two-sided connections of fibroblasts with myocytes. We show, in particular, that the plane-wave conduction velocity CV decreases as a function of G(gap), for zero-sided and one-sided couplings; however, for two-sided coupling, CV decreases initially and then increases as a function of G(gap), and, eventually, we observe that conduction failure occurs for low values of G(gap). In our homogeneous studies, we find that the rotation speed and stability of a spiral wave can be controlled either by controlling G(gap) or E-f. Our studies with fibroblast inhomogeneities show that a spiral wave can get anchored to a local fibroblast inhomogeneity. We also study the efficacy of a low-amplitude control scheme, which has been suggested for the control of spiral-wave turbulence in mathematical models for cardiac tissue, in our MF model both with and without heterogeneities.

Dynamics of edge dislocations in a sheared lamellar mesophase

The dynamics and interactions of edge dislocations in a nearly aligned sheared lamellar mesophase is analysed to provide insights into the relationship between disorder and rheology. First, the mesoscale permeation and momentum equations for the displacement field in the presence of external forces are derived from the model H equations for the concentration and momentum field. The secondary flow generated due to the mean shear around an isolated defect is calculated, and the excess viscosity due to the presence of the defect is determined from the excess energy dissipation due to the secondary flow. The excess viscosity for an isolated defect is found to increase with system size in the cross-stream direction as L-3/2 for an isolated defect, though this divergence is cut-off due to interactions in a defect suspension. As the defects are sheared past each other due to the mean flow, the Peach-Koehler force due to elastic interaction between pairs of defects is found to cause no net displacement relative to each other as they approach from large separation to the distance of closest approach. The equivalent force due to viscous interactions is found to increase the separation for defects of opposite sign, and decrease the separation for defects of same sign. During defect interactions, we find that there is no buckling instability due to dilation of layers for systems of realistic size. However, there is another mechanism, which is the velocity difference generated across a slightly deformed bilayer due to the mean shear, which could result in the creation of new defects. (C) 2013 AIP Publishing LLC.

C-H center dot center dot center dot O-nitrate synthon assisted molecular assembly of hydrogen bonded Ni(II) and Cu(II) complexes

Two Schiff base metal complexes Cu-SPETNNO3 (1) and Ni-SPETNNO3 (2) SPETN=2,2-propane,1,3-diylbis(nitrilomethyldyne)pyridyl,phenolate] ] with hydrogen bonding groups have been synthesized and characterized by single-crystal X-ray diffraction. In both of the compounds nitrates occupy a crystallographic general position. In 1 the lattice nitrates are on the 2(1) screw axis while in 2 they are at the crystallographic inversion center. C-HOnitrate synthons (formed by the nitrate anions and peripheral hydrogen bonding groups of the metal complexes) are non-covalent building blocks in molecular-assembly and packing of the cationic Schiff base metal complexes (M=Ni2+, Cu2+), resulting in 2-D hydrogen bonded networks. The CuCu non-bonding contact in 1 is 3.268 angstrom while the Ni-Ni bonding distance in 2 is 3.437 angstrom.

Confinement induced stochastic sensing of charged coronene and perylene aggregates in alpha-hemolysin nanochannels

Biological nanopores provide optimum dimensions and an optimal environment to study early aggregation kinetics of charged polyaromatic molecules in the nano-confined regime. It is expected that probing early stages of nucleation will enable us to design a strategy for supramolecular assembly and biocrystallization processes. Specifically, we have studied translocation dynamics of coronene and perylene based salts, through the alpha-hemolysin (alpha-HL) protein nanopore. The characteristic blocking events in the time-series signal are a function of concentration and bias voltage. We argue that different blocking events arise due to different aggregation processes as captured by all atomistic molecular dynamics (MD) simulations. These confinement induced aggregations of polyaromatic chromophores during the different stages of translocation are correlated with the spatial symmetry and charge distribution of the molecules.

Layerwise decomposition of water dynamics in reverse micelles: A simulation study of two-dimensional infrared spectrum

We present computer simulation study of two-dimensional infrared spectroscopy (2D-IR) of water confined in reverse micelles (RMs) of various sizes. The present study is motivated by the need to understand the altered dynamics of confined water by performing layerwise decomposition of water, with an aim to quantify the relative contributions of different layers water molecules to the calculated 2D-IR spectrum. The 0-1 transition spectra clearly show substantial elongation, due to in-homogeneous broadening and incomplete spectral diffusion, along the diagonal in the surface water layer of different sized RMs. Fitting of the frequency fluctuation correlation functions reveal that the motion of the surface water molecules is sub-diffusive and indicate the constrained nature of their dynamics. This is further supported by two peak nature of the angular analogue of van Hove correlation function. With increasing system size, the water molecules become more diffusive in nature and spectral diffusion almost completes in the central layer of the larger size RMs. Comparisons between experiments and simulations establish the correspondence between the spectral decomposition available in experiments with the spatial decomposition available in simulations. Simulations also allow a quantitative exploration of the relative role of water, sodium ions, and sulfonate head groups in vibrational dephasing. Interestingly, the negative cross correlation between force on oxygen and hydrogen of O-H bond in bulk water significantly decreases in the surface layer of each RM. This negative cross correlation gradually increases in the central water pool with increasing RMs size and this is found to be partly responsible for the faster relaxation rate of water in the central pool. (C) 2013 AIP Publishing LLC.

In-situ impedance spectroscopy of layer-by-layer self-assembly of polyelectrolytes

In-situ impedance spectroscopy of layer-by-layer self-assembly of weak polyelectrolytes is presented. Interdigitated capacitors with active area of 1×1 mm2 and electrode spacing of 5 μm are fabricated and used for this purpose. Measurement results indicate that the impedance decreases with increase in number of polyelectrolyte layers. About 2.5% of relative change in magnitude of impedance at 104.7 KHz is seen for four bi-layers of Poly(Allylamine Hydrochloride) (PAH)/Poly(Acrylic acid) (PAA). An electrical equivalent for polyelectrolyte binding is obtained.

Enhancing Surface Coverage and Growth in Layer-by-Layer Assembly of Protein Nanoparticles

Thin films of bovine serum albumin (BSA) nanoparticles are fabricated via layer-by-layer assembly. The surface of BSA nanoparticles have two oppositely acting functional groups on the surface: amine (NH2) and carboxylate (COO-). The protonation and deprotonation of these functional groups at different pH vary the charge density on the particle surface, and entirely different growth can be observed by varying the nature of the complementary polymer and the pH of the particles. The complementary polymers used in this study are poly(dimethyldiallylammonium chloride) (PDDAC) and poly(acrylic acid) (PAA). The assembly of BSA nanoparticles based on electrostatic interaction with PDDAC suffers from the poor loading of the nanoparticles. The assembly with PAA aided by a hydrogen bonding interaction shows tremendous improvement in the growth of the assembly over PDDAC. Moreover, the pH of the BSA nanoparticles was observed to affect the loading of nanoparticles in the LbL assembly with PAA significantly.

Switching and release dynamics of an electrostatically actuated MEMS switch under the influence of squeeze-film damping

This paper reports analytical modeling, simulation and experimental validation for switching and release times of an electrostatically actuated micromachined switch. Presented work is an extension of our earlier work [1] that analytically argued, and numerically and experimentally demonstrated, why pull-in time is larger that pull-up time when the actuation voltage is less than twice of the pull-in voltage. In this paper, switching dynamics is investigated under the influence of squeeze-film damping. Tests were performed on SOI (silicon-on-insulator) based parallel beams structures.Typical voltage requirement for actuation is in the range of 10-30 V. All the experiments were performed in normal atmospheric pressure. Measurement results confirm that the quality factor Q has appreciable effect on the release time compared to the switching time. The quality factor Q is extracted from the response measurement and compared with the ANSYS simulation result. In addition, the dynamic pull-in effect has also been studied and reported in this paper. A contribution of this work includes the effect of various phenomena such as squeeze-film damping, dynamic pull-in, and frequency pull-in effects on the switching dynamics of a MEMS switch.

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.

Reaction dynamics under confinement: an exact path integral treatment of a two-stage model of stochastic gene expression

Gene expression in living systems is inherently stochastic, and tends to produce varying numbers of proteins over repeated cycles of transcription and translation. In this paper, an expression is derived for the steady-state protein number distribution starting from a two-stage kinetic model of the gene expression process involving p proteins and r mRNAs. The derivation is based on an exact path integral evaluation of the joint distribution, P(p, r, t), of p and r at time t, which can be expressed in terms of the coupled Langevin equations for p and r that represent the two-stage model in continuum form. The steady-state distribution of p alone, P(p), is obtained from P(p, r, t) (a bivariate Gaussian) by integrating out the r degrees of freedom and taking the limit t -> infinity. P(p) is found to be proportional to the product of a Gaussian and a complementary error function. It provides a generally satisfactory fit to simulation data on the same two-stage process when the translational efficiency (a measure of intrinsic noise levels in the system) is relatively low; it is less successful as a model of the data when the translational efficiency (and noise levels) are high.

Capturability of augmented proportional navigation (APN) guidance with nonlinear engagement dynamics

Proportional Navigation (PN) and its variants are widely used guidance philosophies. However, in the presence of target maneuver, PN guidance law is effective only for a restrictive set of initial geometries. To account for target maneuvers, the concept of Augmented Proportional Navigation (APN) guidance law was introduced and analyzed in a linearized interceptor-target engagement framework presented in literature. However, there is no work in the literature, that addresses the capturability performance of the APN guidance law in a nonlinear engagement framework. This paper presents such an analysis and obtains the conditions for capturability. It also shows that a shorter time of interception is obtained when APN is formulated in the nonlinear framework as proposed in this paper. Simulation results are given to support the theoretical findings.

Concentration Fluctuations and Segmental Dynamics in Weakly Dynamic Asymmetric Blends in the Presence of Surface-Functionalized Multiwall Carbon Nanotubes

Surface-functionalized multiwall carbon nanotubes (MWCNTs) are incorporated in poly(methyl methacrylate)/styrene acrylonitrile (PMMA/SAN) blends and the pretransitional regime is monitored in situ by melt rheology and dielectric spectroscopy. As the blends exhibit weak dynamic asymmetry, the obvious transitions in the melt rheology due to thermal concentration fluctuations are weak. This is further supported by the weak temperature dependence of the correlation length ( approximate to 10-12 angstrom) in the vicinity of demixing. Hence, various rheological techniques in both the temperature and frequency domains are adopted to evaluate the demixing temperature. The spinodal decomposition temperature is manifested in an increase in the miscibility gap in the presence of MWCNTs. Furthermore, MWCNTs lead to a significant slowdown of the segmental dynamics in the blends. Thermally induced phase separation in the PMMA/SAN blends lead to selective localization of MWCNTs in the PMMA phase. This further manifests itself in a significant increase in the melt conductivity.

Efficacious Anticancer Drug Delivery Mediated by a pH-Sensitive Self-Assembly of a Conserved Tripeptide Derived from Tyrosine Kinase NGF Receptor

We present herein a short tripeptide sequence (Lys-Phe-Gly or KFG) that is situated in the juxtamembrane region of the tyrosine kinase nerve growth factor (Trk NGF) receptors. KFG self-assembles in water and shows a reversible and concentration-dependent switching of nanostructures from nanospheres (vesicles) to nanotubes, as evidenced by dynamic light scattering, transmission electron microscopy, and atomic force microscopy. The morphology change was associated with a transition in the secondary structure. The tripeptide vesicles have inner aqueous compartments and are stable at pH7.4 but rupture rapidly at pH approximate to 6. The pH-sensitive response of the vesicles was exploited for the delivery of a chemotherapeutic anticancer drug, doxorubicin, which resulted in enhanced cytotoxicity for both drug-sensitive and drug-resistant cells. Efficient intracellular release of the drug was confirmed by fluorescence-activated cell sorting analysis, fluorescence microscopy, and confocal microscopy.

Molecular Dynamics Perspective on the Protein Thermal Stability: A Case Study Using SAICAR Synthetase

The enzyme SAICAR synthetase ligates aspartate with CAIR (5'-phosphoribosyl-4-carboxy-5-aminoimidazole) forming SAICAR (5-amino-4-imidazole-N-succinocarboxamide ribonucleotide) in the presence of ATP. In continuation with our previous study on the thermostability of this enzyme in hyper-/thermophiles based on the structural aspects, here, we present the dynamic aspects that differentiate the mesophilic (E. coli, E. chaffeensis), thermophilic (G. kaustophilus), and hyperthermophilic (M. jannaschii, P. horikoshii) SAICAR synthetases by carrying out a total of 11 simulations. The five functional dimers from the above organisms were simulated using molecular dynamics for a period of 50 ns each at 300 K, 363 K, and an additional simulation at 333 K for the thermophilic protein. The basic features like root-mean-square deviations, root-mean-square fluctuations, surface accessibility, and radius of gyration revealed the instability of mesophiles at 363 K. Mean square displacements establish the reduced flexibility of hyper-/thermophiles at all temperatures. At the simulations time scale considered here, the long-distance networks are considerably affected in mesophilic structures at 363 K. In mesophiles, a comparatively higher number of short-lived (having less percent existence time) C alpha, hydrogen bonds, hydrophobic interactions are formed, and long-lived (with higher percentage existence time) contacts are lost. The number of time-averaged salt-bridges is at least 2-fold higher in hyperthermophiles at 363 K. The change in surface accessibility of salt-bridges at 363 K from 300 K is nearly doubled in mesophilic protein compared to proteins from other temperature classes.