by Keyword: Dynamics

El Hauadi K, Resina L, Zanuy D, Esteves T, Ferreira FC, Pérez-Madrigal MM, Alemán C, (2022). Dendritic Self-assembled Structures from Therapeutic Charged Pentapeptides Langmuir 38, 12905-12914

CRENKA [Cys-Arg-(NMe)Glu-Lys-Ala, where (NMe)Glu refers to N-methyl-Glu], an anti-cancer pentapeptide that induces prostate tumor necrosis and significant reduction in tumor growth, was engineered to increase the resistance to endogenous proteases of its parent peptide, CREKA (Cys-Arg-Glu-Lys-Ala). Considering their high tendency to aggregate, the self-assembly of CRENKA and CREKA into well-defined and ordered structures has been examined as a function of peptide concentration and pH. Spectroscopic studies and atomistic molecular dynamics simulations reveal significant differences between the secondary structures of CREKA and CRENKA. Thus, the restrictions imposed by the (NMe)Glu residue reduce the conformational variability of CRENKA with respect to CREKA, which significantly affects the formation of well-defined and ordered self-assembly morphologies. Aggregates with poorly defined morphology are obtained from solutions with low and moderate CREKA concentrations at pH 4, whereas well-defined dendritic microstructures with fractal geometry are obtained from CRENKA solutions with similar peptide concentrations at pH 4 and 7. The formation of dendritic structures is proposed to follow a two-step mechanism: (1) pseudo-spherical particles are pre-nucleated through a diffusion-limited aggregation process, pre-defining the dendritic geometry, and (2) such pre-nucleated structures coalesce by incorporating conformationally restrained CRENKA molecules from the solution to their surfaces, forming a continuous dendritic structure. Instead, no regular assembly is obtained from solutions with high peptide concentrations, as their dynamics is dominated by strong repulsive peptide-peptide electrostatic interactions, and from solutions at pH 10, in which the total peptide charge is zero. Overall, results demonstrate that dendritic structures are only obtained when the molecular charge of CRENKA, which is controlled through the pH, favors kinetics over thermodynamics during the self-assembly process.

JTD Keywords: aggregation, amphiphilic peptides, breast-cancer, cells, design, oxidative stress, resistance, strategy, Molecular-dynamics

Cable, J, Arlotta, P, Parker, KK, Hughes, AJ, Goodwin, K, Mummery, CL, Kamm, RD, Engle, SJ, Tagle, DA, Boj, SF, Stanton, AE, Morishita, Y, Kemp, ML, Norfleet, DA, May, EE, Lu, A, Bashir, R, Feinberg, AW, Hull, SM, Gonzalez, AL, Blatchley, MR, Pulido, NM, Morizane, R, McDevitt, TC, Mishra, D, Mulero-Russe, A, (2022). Engineering multicellular living systems-A Keystone Symposia report Annals Of The New York Academy Of Sciences , 1-13

The ability to engineer complex multicellular systems has enormous potential to inform our understanding of biological processes and disease and alter the drug development process. Engineering living systems to emulate natural processes or to incorporate new functions relies on a detailed understanding of the biochemical, mechanical, and other cues between cells and between cells and their environment that result in the coordinated action of multicellular systems. On April 3-6, 2022, experts in the field met at the Keystone symposium "Engineering Multicellular Living Systems" to discuss recent advances in understanding how cells cooperate within a multicellular system, as well as recent efforts to engineer systems like organ-on-a-chip models, biological robots, and organoids. Given the similarities and common themes, this meeting was held in conjunction with the symposium "Organoids as Tools for Fundamental Discovery and Translation".

JTD Keywords: Brain organoids, Cell diversity, Computational, Dynamics, Engineered living, Engineered organs, Heart, Maturation, Model, Multicellular, Mycobacterium-tuberculosis, Quantitative-analysis, Systems, Tissue deformation

De Corato, M, Arroyo, M, (2022). A theory for the flow of chemically responsive polymer solutions: Equilibrium and shear-induced phase separation Journal Of Rheology 66, 813-835

Chemically responsive polymers are macromolecules that respond to local variations of the chemical composition of the solution by changing their conformation, with notable examples including polyelectrolytes, proteins, and DNA. The polymer conformation changes can occur in response to changes in the pH, the ionic strength, or the concentration of a generic solute that interacts with the polymer. These chemical stimuli can lead to drastic variations of the polymer flexibility and even trigger a transition from a coil to a globule polymer conformation. In many situations, the spatial distribution of the chemical stimuli can be highly inhomogeneous, which can lead to large spatial variations of polymer conformation and of the rheological properties of the mixture. In this paper, we develop a theory for the flow of a mixture of solute and chemically responsive polymers. The approach is valid for generic flows and inhomogeneous distributions of polymers and solutes. To model the polymer conformation changes introduced by the interactions with the solute, we consider the polymers as linear elastic dumbbells whose spring stiffness depends on the solute concentration. We use Onsager's variational formalism to derive the equations governing the evolution of the variables, which unveils novel couplings between the distribution of dumbbells and that of the solute. Finally, we use a linear stability analysis to show that the governing equations predict an equilibrium phase separation and a distinct shear-induced phase separation whereby a homogeneous distribution of solute and dumbbells spontaneously demix. Similar phase transitions have been observed in previous experiments using stimuli-responsive polymers and may play an important role in living systems. (C) 2022 The Society of Rheology.

JTD Keywords: Coil-globule transition, Constitutive equation, Dilute-solutions, Dumbbell model, Dynamics, Macromolecules, Nonequilibrium thermodynamics, Polyelectrolytes, Polymer migration, Polymer phase separation, Polymers, Predictions, Rheology, Shear-induced phase separation, Solute-polymer interactions, Stress, Viscoelasticity

Mesquida-Veny, F, Martinez-Torres, S, Del Rio, JA, Hervera, A, (2022). Nociception-Dependent CCL21 Induces Dorsal Root Ganglia Axonal Growth via CCR7-ERK Activation Frontiers In Immunology 13, 880647

While chemokines were originally described for their ability to induce cell migration, many studies show how these proteins also take part in many other cell functions, acting as adaptable messengers in the communication between a diversity of cell types. In the nervous system, chemokines participate both in physiological and pathological processes, and while their expression is often described on glial and immune cells, growing evidence describes the expression of chemokines and their receptors in neurons, highlighting their potential in auto- and paracrine signalling. In this study we analysed the role of nociception in the neuronal chemokinome, and in turn their role in axonal growth. We found that stimulating TRPV1(+) nociceptors induces a transient increase in CCL21. Interestingly we also found that CCL21 enhances neurite growth of large diameter proprioceptors in vitro. Consistent with this, we show that proprioceptors express the CCL21 receptor CCR7, and a CCR7 neutralizing antibody dose-dependently attenuates CCL21-induced neurite outgrowth. Mechanistically, we found that CCL21 binds locally to its receptor CCR7 at the growth cone, activating the downstream MEK-ERK pathway, that in turn activates N-WASP, triggering actin filament ramification in the growth cone, resulting in increased axonal growth.

JTD Keywords: Actin dynamics, Axonal growth, Ccl21, Ccr7, Cell-migration, Central-nervous-system, Chemokine, Ligands, Mek-erk, Microglia, Neurons, Neuropathic pain, Nociception, Phosphorylation, Regeneration

Wagner, AM, Eto, H, Joseph, A, Kohyama, S, Haraszti, T, Zamora, RA, Vorobii, M, Giannotti, MI, Schwille, P, Rodriguez-Emmenegger, C, (2022). Dendrimersome Synthetic Cells Harbor Cell Division Machinery of Bacteria Advanced Materials 34, 2202364

The integration of active cell machinery with synthetic building blocks is the bridge toward developing synthetic cells with biological functions and beyond. Self-replication is one of the most important tasks of living systems, and various complex machineries exist to execute it. In Escherichia coli, a contractile division ring is positioned to mid-cell by concentration oscillations of self-organizing proteins (MinCDE), where it severs membrane and cell wall. So far, the reconstitution of any cell division machinery has exclusively been tied to liposomes. Here, the reconstitution of a rudimentary bacterial divisome in fully synthetic bicomponent dendrimersomes is shown. By tuning the membrane composition, the interaction of biological machinery with synthetic membranes can be tailored to reproduce its dynamic behavior. This constitutes an important breakthrough in the assembly of synthetic cells with biological elements, as tuning of membrane-divisome interactions is the key to engineering emergent biological behavior from the bottom-up.

JTD Keywords: Bacterial cell division, Bottom-up synthetic biology, Dendrimersomes, Dynamic min patterns, Dynamics, Ftsz assembly, Ftsz filaments, Mind, Organization, Pole oscillation, Polymersome membranes, Proteins, Rapid pole, Synthetic cells, Vesicles

Tuveri, GM, Ceccarelli, M, Pira, A, Bodrenko, IV, (2022). The Optimal Permeation of Cyclic Boronates to Cross the Outer Membrane via the Porin Pathway Antibiotics 11, 840

We investigated the diffusion of three cyclic boronates formulated as beta-lactamase inhibitors through the porin OmpF to evaluate their potential to cross OM via the porin pathway. The three nonbeta-lactam molecules diffuse through the porin eyelet region with the same mechanism observed for beta-lactam molecules and diazobicyclooctan derivatives, with the electric dipole moment aligned with the transversal electric field. In particular, the BOH group can interact with both the basic ladder and the acidic loop L3, which is characteristic of the size-constricted region of this class of porins. On one hand, we confirm that the transport of small molecules through enterobacter porins has a common general mechanism; on the other, the class of cyclic boronate molecules does not seem to have particular difficulties in diffusing through enterobacter porins, thus representing a good scaffold for new anti-infectives targeting Gram-negative bacteria research.

JTD Keywords: Antibiotics, Beta-lactamase inhibitors, Cyclic boronates, Diffusion, Diffusion current, Discovery, Electric-field, Metadynamics, Molecular dynamics simulations, Molecular-dynamics simulations, Nanopores, Permeability, Permeation, Porins, Rules, Translocation

Casanellas, I, Lagunas, A, Vida, Y, Perez-Inestrosa, E, Rodriguez-Pereira, C, Magalhaes, J, Andrades, JA, Becerra, J, Samitier, J, (2022). Nanoscale ligand density modulates gap junction intercellular communication of cell condensates during chondrogenesis Nanomedicine 17, 775-791

Aim: To unveil the influence of cell-matrix adhesions in the establishment of gap junction intercellular communication (GJIC) during cell condensation in chondrogenesis. Materials & methods: Previously developed nanopatterns of the cell adhesive ligand arginine-glycine-aspartic acid were used as cell culture substrates to control cell adhesion at the nanoscale. In vitro chondrogenesis of mesenchymal stem cells was conducted on the nanopatterns. Cohesion and GJIC were evaluated in cell condensates. Results: Mechanical stability and GJIC are enhanced by a nanopattern configuration in which 90% of the surface area presents adhesion sites separated less than 70 nm, thus providing an onset for cell signaling. Conclusion: Cell-matrix adhesions regulate GJIC of mesenchymal cell condensates during in vitro chondrogenesis from a threshold configuration at the nanoscale.

JTD Keywords: Actin, Adhesion, Arginine-glycine-aspartic acid, Cell adhesion, Collagen, Condensation, Connexin-43, Dendrimer-based nanopatterning, Dynamics, Extracellular-matrix, Fibronectin, Gap junction intercellular communication, Mesenchymal stem cells, Permeability, Phenotype, Vinculin

Rosales-Rojas, R, Zuniga-Bustos, M, Salas-Sepulveda, F, Galaz-Araya, C, Zamora, RA, Poblete, H, (2022). Self-Organization Dynamics of Collagen-like Peptides Crosslinking Is Driven by Rose-Bengal-Mediated Electrostatic Bridges Pharmaceutics 14, 1148

The present work focuses on the computational study of the structural micro-organization of hydrogels based on collagen-like peptides (CLPs) in complex with Rose Bengal (RB). In previous studies, these hydrogels computationally and experimentally demonstrated that when RB was activated by green light, it could generate forms of stable crosslinked structures capable of regenerating biological tissues such as the skin and cornea. Here, we focus on the structural and atomic interactions of two collagen-like peptides (collagen-like peptide I (CLPI), and collagen-like peptide II, (CLPII)) in the presence and absence of RB, highlighting the acquired three-dimensional organization and going deep into the stabilization effect caused by the dye. Our results suggest that the dye could generate a ternary ground-state complex between collagen-like peptide fibers, specifically with positively charged amino acids (Lys in CLPI and Arg in CLPII), thus stabilizing ordered three-dimensional structures. The discoveries generated in this study provide the structural and atomic bases for the subsequent rational development of new synthetic peptides with improved characteristics for applications in the regeneration of biological tissues during photochemical tissue bonding therapies.

JTD Keywords: Anastomosis, Collagen-like peptide, Crosslinking, Green light, Mm simulations, Molecular dynamics, Molecular-dynamics, Photochemical tissue bonding therapies, Qm, Rose bengal

Ebisuya, M, Trepat, X, (2022). Tension hones body segmentation around the clock Nature 605, 432-433

JTD Keywords: Dynamics

Martí D, Martín-Martínez E, Torras J, Betran O, Turon P, Alemán C, (2022). In silico study of substrate chemistry effect on the tethering of engineered antibodies for SARS-CoV-2 detection: Amorphous silica vs gold Colloids And Surfaces B-Biointerfaces 213, 112400

The influence of the properties of different solid substrates on the tethering of two antibodies, IgG1-CR3022 and IgG1-S309, which were specifically engineered for the detection of SARS-CoV-2, has been examined at the molecular level using conventional and accelerated Molecular Dynamics (cMD and aMD, respectively). Two surfaces with very different properties and widely used in immunosensors for diagnosis, amorphous silica and the most stable facet of the face-centered cubic gold structure, have been considered. The effects of such surfaces on the structure and orientation of the immobilized antibodies have been determined by quantifying the tilt and hinge angles that describe the orientation and shape of the antibody, respectively, and the dihedrals that measure the relative position of the antibody arms with respect to the surface. Results show that the interactions with amorphous silica, which are mainly electrostatic due to the charged nature of the surface, help to preserve the orientation and structure of the antibodies, especially of the IgG1-CR3022, indicating that the primary sequence of those antibodies also plays some role. Instead, short-range van der Waals interactions with the inert gold surface cause a higher degree tilting and fraying of the antibodies with respect to amorphous silica. The interactions between the antibodies and the surface also affect the correlation among the different angles and dihedrals, which increases with their strength. Overall, results explain why amorphous silica substrates are frequently used to immobilize antibodies in immunosensors. © 2022 The Authors

JTD Keywords: amorphous silica, enzyme, gol d, immobilization, immunosensor, molecu l a r dynamics, protein adsorption, sars-cov-2 immunosensor, simulations, spike protein, surface interactions, target, vaccine, Amorphous silica, Antibodies, Antibody engineering, Antibody immobilization, Antibody structure, Article, Chemical detection, Computer model, Controlled study, Dihedral angle, Gold, In-silico, Molecular dynamics, Molecular levels, Molecular-dynamics, Nonhuman, Property, Sars, Sars-cov-2 immunosensor, Severe acute respiratory syndrome coronavirus 2, Silica, Silico studies, Silicon dioxide, Solid substrates, Structure analysis, Substrate chemistry, Substrates, Van der waals forces, Virus detection

Espinoso A, Andrzejak RG, (2022). Phase irregularity: A conceptually simple and efficient approach to characterize electroencephalographic recordings from epilepsy patients Physical Review e 105, 34212

The severe neurological disorder epilepsy affects almost 1% of the world population. For patients who suffer from pharmacoresistant focal-onset epilepsy, electroencephalographic (EEG) recordings are essential for the localization of the brain area where seizures start. Apart from the visual inspection of the recordings, quantitative EEG signal analysis techniques proved to be useful for this purpose. Among other features, regularity versus irregularity and phase coherence versus phase independence allowed characterizing brain dynamics from the measured EEG signals. Can phase irregularities also characterize brain dynamics? To address this question, we use the univariate coefficient of phase velocity variation, defined as the ratio of phase velocity standard deviation and the mean phase velocity. Beyond that, as a bivariate measure we use the classical mean phase coherence to quantify the degree of phase locking. All phase-based measures are combined with surrogates to test null hypotheses about the dynamics underlying the signals. In the first part of our analysis, we use the Rössler model system to study our approach under controlled conditions. In the second part, we use the Bern-Barcelona EEG database which consists of focal and nonfocal signals extracted from seizure-free recordings. Focal signals are recorded from brain areas where the first seizure EEG signal changes can be detected, and nonfocal signals are recorded from areas that are not involved in the seizure at its onset. Our results show that focal signals have less phase variability and more phase coherence than nonfocal signals. Once combined with surrogates, the mean phase velocity proved to have the highest discriminative power between focal and nonfocal signals. In conclusion, conceptually simple and easy to compute phase-based measures can help to detect features induced by epilepsy from EEG signals. This holds not only for the classical mean phase coherence but even more so for univariate measures of phase irregularity. © 2022 American Physical Society.

JTD Keywords: brain, entropy, epileptogenic networks, functional connectivity, hilbert transform, seizure onset, surrogate data, synchronization, time-series, Biomedical signal processing, Brain areas, Brain dynamics, Dynamics, Electroencephalographic signals, Electroencephalography, Electrophysiology, Intracranial eeg signals, Localisation, Neurological disorders, Neurology, Phase based, Phase coherence, Signal detection, Simple++, Univariate, Velocity, World population

Valenti, Sofia, del Valle, Luis Javier, Romanini, Michela, Mitjana, Meritxell, Puiggalí, Jordi, Tamarit, Josep Lluís, Macovez, Roberto, (2022). Drug-Biopolymer Dispersions: Morphology- and Temperature- Dependent (Anti)Plasticizer Effect of the Drug and Component-Specific Johari–Goldstein Relaxations International Journal Of Molecular Sciences 23, 2456

Amorphous molecule-macromolecule mixtures are ubiquitous in polymer technology and are one of the most studied routes for the development of amorphous drug formulations. For these applications it is crucial to understand how the preparation method affects the properties of the mixtures. Here, we employ differential scanning calorimetry and broadband dielectric spectroscopy to investigate dispersions of a small-molecule drug (the Nordazepam anxiolytic) in biodegradable polylactide, both in the form of solvent-cast films and electrospun microfibres. We show that the dispersion of the same small-molecule compound can have opposite (plasticizing or antiplasticizing) effects on the segmental mobility of a biopolymer depending on preparation method, temperature, and polymer enantiomerism. We compare two different chiral forms of the polymer, namely, the enantiomeric pure, semicrystalline L-polymer (PLLA), and a random, fully amorphous copolymer containing both L and D monomers (PDLLA), both of which have lower glass transition temperature (Tg) than the drug. While the drug has a weak antiplasticizing effect on the films, consistent with its higher Tg, we find that it actually acts as a plasticizer for the PLLA microfibres, reducing their Tg by as much as 14 K at 30%-weight drug loading, namely, to a value that is lower than the Tg of fully amorphous films. The structural relaxation time of the samples similarly depends on chemical composition and morphology. Most mixtures displayed a single structural relaxation, as expected for homogeneous samples. In the PLLA microfibres, the presence of crystalline domains increases the structural relaxation time of the amorphous fraction, while the presence of the drug lowers the structural relaxation time of the (partially stretched) chains in the microfibres, increasing chain mobility well above that of the fully amorphous polymer matrix. Even fully amorphous homogeneous mixtures exhibit two distinct Johari–Goldstein relaxation processes, one for each chemical component. Our findings have important implications for the interpretation of the Johari–Goldstein process as well as for the physical stability and mechanical properties of microfibres with small-molecule additives.

JTD Keywords: amorphous pharmaceuticals, beta-relaxation, constant loss, crystallization, dielectric spectroscopy, dynamics, formulation morphology, glass transition, molecular mobility, nanofibers, polylactide, polymer enantiomerism, secondary relaxations, valium metabolite, viscous-liquids, Amorphous pharmaceuticals, Glass-transition, Secondary relaxations

Dhiman, Shikha, Andrian, Teodora, Gonzalez, Beatriz Santiago, Tholen, Marrit ME., Wang, Yuyang, Albertazzi, Lorenzo, (2022). Can super-resolution microscopy become a standard characterization technique for materials chemistry? Chemical Science 13, 2152-2166

The characterization of newly synthesized materials is a cornerstone of all chemistry and nanotechnology laboratories. For this purpose, a wide array of analytical techniques have been standardized and are used routinely by laboratories across the globe. With these methods we can understand the structure, dynamics and function of novel molecular architectures and their relations with the desired performance, guiding the development of the next generation of materials. Moreover, one of the challenges in materials chemistry is the lack of reproducibility due to improper publishing of the sample preparation protocol. In this context, the recent adoption of the reporting standard MIRIBEL (Minimum Information Reporting in Bio–Nano Experimental Literature) for material characterization and details of experimental protocols aims to provide complete, reproducible and reliable sample preparation for the scientific community. Thus, MIRIBEL should be immediately adopted in publications by scientific journals to overcome this challenge. Besides current standard spectroscopy and microscopy techniques, there is a constant development of novel technologies that aim to help chemists unveil the structure of complex materials. Among them super-resolution microscopy (SRM), an optical technique that bypasses the diffraction limit of light, has facilitated the study of synthetic materials with multicolor ability and minimal invasiveness at nanometric resolution. Although still in its infancy, the potential of SRM to unveil the structure, dynamics and function of complex synthetic architectures has been highlighted in pioneering reports during the last few years. Currently, SRM is a sophisticated technique with many challenges in sample preparation, data analysis, environmental control and automation, and moreover the instrumentation is still expensive. Therefore, SRM is currently limited to expert users and is not implemented in characterization routines. This perspective discusses the potential of SRM to transition from a niche technique to a standard routine method for material characterization. We propose a roadmap for the necessary developments required for this purpose based on a collaborative effort from scientists and engineers across disciplines.

JTD Keywords: blinking, fluorophore, intramolecular spirocyclization, localization, nanoparticles, resolution limit, reveals, single-molecule fluorescence, stimulated-emission, Characterization techniques, Diffraction, Distributed computer systems, Environmental management, Information reporting, Material chemistry, Materials characterization, Minimum information, Optical reconstruction microscopy, Optical resolving power, Sample preparation, Structure dynamics, Structure functions, Super-resolution microscopy, Synthesized materials

Bertran, Oscar, Martí, Didac, Torras, Juan, Turon, Pau, Alemán, Carlos, (2022). Computer simulations on oxidative stress-induced reactions in SARS-CoV-2 spike glycoprotein: a multi-scale approach Molecular Diversity ,

Abstract Oxidative stress, which occurs when an organism is exposed to an adverse stimulus that results in a misbalance of antioxidant and pro-oxidants species, is the common denominator of diseases considered as a risk factor for SARS-CoV-2 lethality. Indeed, reactive oxygen species caused by oxidative stress have been related to many virus pathogenicity. In this work, simulations have been performed on the receptor binding domain of SARS-CoV-2 spike glycoprotein to study what residues are more susceptible to be attacked by ·OH, which is one of the most reactive radicals associated to oxidative stress. The results indicate that isoleucine (ILE) probably plays a crucial role in modification processes driven by radicals. Accordingly, QM/MM-MD simulations have been conducted to study both the ·OH-mediated hydrogen abstraction of ILE residues and the induced modification of the resulting ILE radical through hydroxylation or nitrosylation reactions. All in all, in silico studies show the importance of the chemical environment triggered by oxidative stress on the modifications of the virus, which is expected to help for foreseeing the identification or development of antioxidants as therapeutic drugs. Graphic abstract

JTD Keywords: atom abstraction, damage, density functionals, hydrogen abstraction, isoleucine, molecular dynamics, pathogenesis, protein, reactive oxygen species, receptor binding domain, residues, spike protein, Amino-acids, Hydrogen abstraction, Isoleucine, Molecular dynamics, Reactive oxygen species, Receptor binding domain, Spike protein

Garreta E, Nauryzgaliyeva Z, Marco A, Safi W, Montserrat N, (2022). Dissecting nephron morphogenesis using kidney organoids from human pluripotent stem cells Current Opinion In Genetics & Development 72, 22-29

During kidney development the emergence of complex multicellular shapes such as the nephron (the functional unit of the kidney) rely on spatiotemporally coordinated developmental programs. These involve gene regulatory networks, signaling pathways and mechanical forces, that work in concert to shape and form the nephron(s). The generation of kidney organoids from human pluripotent stem cells now represent an unprecedented experimental set up to study these processes. Here we discuss the potential applications of kidney organoids to advance our knowledge of how mechanical forces and cell fate specification spatiotemporally interact to orchestrate nephron patterning and morphogenesis in humans. Progress in innovative techniques for quantifying and perturbing these processes in a controlled manner will be crucial. A mechanistic understanding of the multicellular dynamic processes occurring during nephrogenesis will pave the way to unveil new mechanisms of human kidney disease. © 2021

JTD Keywords: differentiation, dynamics, induction, lumen formation, models, mouse, organogenesis, reveals, tubules, Divergent features

Martí, Didac, Alemán, Carlos, Ainsley, Jon, Ahumada, Oscar, Torras, Juan, (2022). IgG1-b12–HIV-gp120 Interface in Solution: A Computational Study Journal Of Chemical Information And Modeling 62, 359-371

The use of broadly neutralizing antibodies against human immunodeficiency virus type 1 (HIV-1) has been shown to be a promising therapeutic modality in the prevention of HIV infection. Understanding the b12-gp120 binding mechanism under physiological conditions may assist the development of more broadly effective antibodies. In this work, the main conformations and interactions between the receptor-binding domain (RBD) of spike glycoprotein gp120 of HIV-1 and the IgG1-b12 mAb are studied. Accelerated molecular dynamics (aMD) and ab initio hybrid molecular dynamics have been combined to determine the most persistent interactions between the most populated conformations of the antibody-antigen complex under physiological conditions. The results show the most persistent receptor-binding mapping in the conformations of the antibody-antigen interface in solution. The binding-free-energy decomposition reveals a small enhancement in the contribution played by the CDR-H3 region to the b12-gp120 interface compared to the crystal structure.

JTD Keywords: antibody, complex, functionals, gp120 envelope glycoprotein, hiv, immunodeficiency-virus, noncovalent interactions, simulations, software integration, Ab initio, Accelerated molecular dynamics, Accelerated molecular-dynamics, Antibodies, Antigens, Binding energy, Binding mechanisms, Computational studies, Crystal structure, Diseases, Free energy, Hiv infection, Human immunodeficiency virus, Molecular dynamics, Neutralizing antibodies, Physiological condition, Physiology, Receptor-binding domains, Therapeutic modality, Viruses

Amil, Adrián Fernández, Verschure, Paul F.M.J., (2021). Supercritical dynamics at the edge-of-chaos underlies optimal decision-making Journal Of Physics-Complexity 2,

Abstract Critical dynamics, characterized by scale-free neuronal avalanches, is thought to underlie optimal function in the sensory cortices by maximizing information transmission, capacity, and dynamic range. In contrast, deviations from criticality have not yet been considered to support any cognitive processes. Nonetheless, neocortical areas related to working memory and decision-making seem to rely on long-lasting periods of ignition-like persistent firing. Such firing patterns are reminiscent of supercritical states where runaway excitation dominates the circuit dynamics. In addition, a macroscopic gradient of the relative density of Somatostatin (SST+) and Parvalbumin (PV+) inhibitory interneurons throughout the cortical hierarchy has been suggested to determine the functional specialization of low- versus high-order cortex. These observations thus raise the question of whether persistent activity in high-order areas results from the intrinsic features of the neocortical circuitry. We used an attractor model of the canonical cortical circuit performing a perceptual decision-making task to address this question. Our model reproduces the known saddle-node bifurcation where persistent activity emerges, merely by increasing the SST+/PV+ ratio while keeping the input and recurrent excitation constant. The regime beyond such a phase transition renders the circuit increasingly sensitive to random fluctuations of the inputs -i.e., chaotic-, defining an optimal SST+/PV+ ratio around the edge-of-chaos. Further, we show that both the optimal SST+/PV+ ratio and the region of the phase transition decrease monotonically with increasing input noise. This suggests that cortical circuits regulate their intrinsic dynamics via inhibitory interneurons to attain optimal sensitivity in the face of varying uncertainty. Hence, on the one hand, we link the emergence of supercritical dynamics at the edge-of-chaos to the gradient of the SST+/PV+ ratio along the cortical hierarchy, and, on the other hand, explain the behavioral effects of the differential regulation of SST+ and PV+ interneurons by neuromodulators like acetylcholine in the presence of input uncertainty.

JTD Keywords: attractor model, cortex, cortical networks, edge-of-chaos, model, nmda receptors, Attractor model, Cortical hierarchies, Decision making, Dynamics, Edge of chaos, Edge-of-chaos, High-order, Higher-order, Inhibitory interneurons, Neurons, Optimal decision making, Persistent activities, Persistent activity, Supercritical, Supercriticality

Casellas NM, Albertazzi L, Pujals S, Torres T, García-Iglesias M, (2021). Unveiling Polymerization Mechanism in pH-regulated Supramolecular Fibers in Aqueous Media Chemistry-A European Journal 27, 11056-11060

An amine functionalized C3-symmetric benzotrithiophene (BTT) monomer has been designed and synthetized in order to form pH responsive one-dimensional supramolecular polymers in aqueous media. While most of the reported studies looked at the effect of pH on the size of the aggregates, herein, a detailed mechanistic study is reported, carried out upon modifying the pH to trigger the formation of positively charged ammonium groups. A dramatic and reversible change in the polymerization mechanism and size of the supramolecular fibers is observed and ascribed to the combination of Coulombic repulsive forces and higher monomer solubility. Furthermore, the induced frustrated growth of the fibers is further employed to finely control the one-dimensional supramolecular polymerisation and copolymerization processes.

JTD Keywords: dynamics, ph responsivity, polymerization mechanism, self-assembly, supramolecular chemistry, supramolecular polymers, Ph responsivity, Polymerization mechanism, Polymers, Self-assembly, Supramolecular chemistry, Supramolecular polymers

Estefan DP, Zucca R, Arsiwalla X, Principe A, Zhang H, Rocamora R, Axmacher N, Verschure PFMJ, (2021). Volitional learning promotes theta phase coding in the human hippocampus Proceedings Of The National Academy Of Sciences Of The United States Of America 118, e2021238118

© 2021 National Academy of Sciences. All rights reserved. Electrophysiological studies in rodents show that active navigation enhances hippocampal theta oscillations (4–12 Hz), providing a temporal framework for stimulus-related neural codes. Here we show that active learning promotes a similar phase coding regime in humans, although in a lower frequency range (3–8 Hz). We analyzed intracranial electroencephalography (iEEG) from epilepsy patients who studied images under either volitional or passive learning conditions. Active learning increased memory performance and hippocampal theta oscillations and promoted a more accurate reactivation of stimulus-specific information during memory retrieval. Representational signals were clustered to opposite phases of the theta cycle during encoding and retrieval. Critically, during active but not passive learning, the temporal structure of intracycle reactivations in theta reflected the semantic similarity of stimuli, segregating conceptually similar items into more distant theta phases. Taken together, these results demonstrate a multilayered mechanism by which active learning improves memory via a phylogenetically old phase coding scheme.

JTD Keywords: active learning, dynamics, gamma-power, hippocampus, intracranial eeg, movement, navigation, neural phase coding, oscillations, representations, retrieval, rhythm, theta oscillations, toolbox, Active learning, Theta oscillations, Working-memory

Marti, D, Martin-Martinez, E, Torras, J, Bertran, O, Turon, P, Aleman, C, (2021). In silico antibody engineering for SARS-CoV-2 detection Computational And Structural Biotechnology Journal 19, 5525-5534

Engineered immunoglobulin-G molecules (IgGs) are of wide interest for the development of detection elements in protein-based biosensors with clinical applications. The strategy usually employed for the de novo design of such engineered IgGs consists on merging fragments of the three-dimensional structure of a native IgG, which is immobilized on the biosensor surface, and of an antibody with an exquisite target specificity and affinity. In this work conventional and accelerated classical molecular dynamics (cMD and aMD, respectively) simulations have been used to propose two IgG-like antibodies for COVID-19 detection. More specifically, the crystal structure of the IgG1 B12 antibody, which inactivates the human immunodeficiency virus-1, has been merged with the structure of the antibody CR3022 Fab tightly bounded to SARS-CoV-2 receptor-binding domain (RBD) and the structure of the 5309 antibody Fab fragment complexed with SARS-CoV-2 RBD. The two constructed antibodies, named IgG1-CR3022 and IgG1-S309, respectively, have been immobilized on a stable gold surface through a linker. Analyses of the influence of both the merging strategy and the substrate on the stability of the two constructs indicate that the IgG1-S309 antibody better preserves the neutralizing structure than the IgG1-CR3022 one. Overall, results indicate that the IgG1-S309 is appropriated for the generation of antibody based sensors for COVID-19 diagnosis. (C) 2021 The Author(s). Published by Elsevier B.V. on behalf of Research Network of Computational and Structural Biotechnology.

JTD Keywords: cr3022, igg1, molecular engineering, s309, Antibodies, Antibody engineering, Biosensors, Chemical detection, Clinical application, Cov, Cr3022, Crystal structure, Design, Diseases, Gold nanoparticles, Igg1, Igg1 antibody, Immobilization, Immunoglobulin g, Immunosensor, In-silico, Merging, Molecular dynamics, Molecular engineering, Orientation, Protein-based biosensors, Receptor-binding domains, S309, Sars, Sensor, Spike protein, Target, Vaccine, Viruses

Bertran, O., Saldías, C., Díaz, D. D., Alemán, C., (2020). Molecular dynamics simulations on self-healing behavior of ionene polymer-based nanostructured hydrogels Polymer 211, 123072

The microscopic mechanism accounting for the self-healing attribute of aromatic ionene-forming hydrogels derived from 1,4-diazabicyclo [2.2.2]octane (DABCO) and N,N’-(x-phenylene)dibenzamide (x = ortho-/meta-/para-) is unknown. Interestingly, the self-healing property of such DABCO-containing hydrogels is largely dependent on the polymer topology, the ortho ionene being the only self-healable without adding oppositely charged species. In this work, Molecular Dynamics (MD) simulations have been conducted to evaluate the influence of the topology on ionene···ionene and ionene··water interactions, as well as their effect on the self-healing behavior. For this purpose, destabilized and structurally damaged models were produced for ionene hydrogels with ortho, meta and para topologies and used as starting geometries for simulations. These models were allowed to evolve without any restriction during MD production runs and, subsequently, the temporal evolution of ionene···ionene and water···ionene interactions was examined. Analysis of the results indicated that the ortho-isomer rapidly forms unique interactions that are not detected for other two isomers. Thus, in addition to the interactions also identified for the meta-and para-ionenes, the ortho-isomer exhibits the formation of strong intermolecular three-centered (N–)H⋯O (=C)⋯H (–N) hydrogen bonds, intramolecular planar sandwich π-π stacking interactions and Cl−···N+ electrostatic interactions. Furthermore, the amount of intermolecular π-π stacking interactions and the strength of water···polymer interaction are also influenced by the topology, favoring the stabilization of the ortho-ionene reconstituted hydrogels. Overall, the arrangement of the functional groups in the ortho topology favors the formation of more types of ionene···ionene interactions, as well as stronger interactions, than in the meta and para topologies.

JTD Keywords: DABCO, Econstituted hydrogels, Molecular dynamics, Polyelectrolyte hydrogels, Self-healing mechanism

Raote, Ishier, Chabanon, Morgan, Walani, Nikhil, Arroyo, Marino, Garcia-Parajo, Maria F., Malhotra, Vivek, Campelo, Felix, (2020). A physical mechanism of TANGO1-mediated bulky cargo export eLife 9, e59426

The endoplasmic reticulum (ER)-resident protein TANGO1 assembles into a ring around ER exit sites (ERES), and links procollagens in the ER lumen to COPII machinery, tethers, and ER-Golgi intermediate compartment (ERGIC) in the cytoplasm (Raote et al., 2018). Here, we present a theoretical approach to investigate the physical mechanisms of TANGO1 ring assembly and how COPII polymerization, membrane tension, and force facilitate the formation of a transport intermediate for procollagen export. Our results indicate that a TANGO1 ring, by acting as a linactant, stabilizes the open neck of a nascent COPII bud. Elongation of such a bud into a transport intermediate commensurate with bulky procollagens is then facilitated by two complementary mechanisms: (i) by relieving membrane tension, possibly by TANGO1-mediated fusion of retrograde ERGIC membranes and (ii) by force application. Altogether, our theoretical approach identifies key biophysical events in TANGO1-driven procollagen export.

JTD Keywords: Membrane tension, Procollagen export, Secretory pathway, Membrane curvature, Membrane dynamics, Budding

Redondo-Morata, Lorena, Losada-Pérez, Patricia, Giannotti, Marina Inés, (2020). Lipid bilayers: Phase behavior and nanomechanics Current Topics in Membranes (ed. Levitan, Irena, Trache, Andreea), Academic Press (Berlin, Germany) 86, 1-55

Lipid membranes are involved in many physiological processes like recognition, signaling, fusion or remodeling of the cell membrane or some of its internal compartments. Within the cell, they are the ultimate barrier, while maintaining the fluidity or flexibility required for a myriad of processes, including membrane protein assembly. The physical properties of in vitro model membranes as model cell membranes have been extensively studied with a variety of techniques, from classical thermodynamics to advanced modern microscopies. Here we review the nanomechanics of solid-supported lipid membranes with a focus in their phase behavior. Relevant information obtained by quartz crystal microbalance with dissipation monitoring (QCM-D) and atomic force microscopy (AFM) as complementary techniques in the nano/mesoscale interface is presented. Membrane morphological and mechanical characterization will be discussed in the framework of its phase behavior, phase transitions and coexistence, in simple and complex models, and upon the presence of cholesterol.

JTD Keywords: Lipid phase behavior, Phase transition, Phase coexistence, Nanomechanics, Thermodynamics, Atomic force microscopy (AFM), Quartz crystal microbalance with dissipation monitoring (QCM-D)

Lozano-García, M., Estrada-Petrocelli, L., Moxham, J., Rafferty, G. F., Torres, A., Jolley, C. J., Jané, R. , (2019). Noninvasive assessment of inspiratory muscle neuromechanical coupling during inspiratory threshold loading IEEE Access 7, 183634-183646

Diaphragm neuromechanical coupling (NMC), which reflects the efficiency of conversion of neural activation to transdiaphragmatic pressure (Pdi), is increasingly recognized to be a useful clinical index of diaphragm function and respiratory mechanics in neuromuscular weakness and cardiorespiratory disease. However, the current gold standard assessment of diaphragm NMC requires invasive measurements of Pdi and crural diaphragm electromyography (oesEMGdi), which complicates the measurement of diaphragm NMC in clinical practice. This is the first study to compare invasive measurements of diaphragm NMC (iNMC) using the relationship between Pdi and oesEMGdi, with noninvasive assessment of NMC (nNMC) using surface mechanomyography (sMMGlic) and electromyography (sEMGlic) of lower chest wall inspiratory muscles. Both invasive and noninvasive measurements were recorded in twelve healthy adult subjects during an inspiratory threshold loading protocol. A linear relationship between noninvasive sMMGlic and sEMGlic measurements was found, resulting in little change in nNMC with increasing inspiratory load. By contrast, a curvilinear relationship between invasive Pdi and oesEMGdi measurements was observed, such that there was a progressive increase in iNMC with increasing inspiratory threshold load. Progressive recruitment of lower ribcage muscles, serving to enhance the mechanical advantage of the diaphragm, may explain the more linear relationship between sMMGlic and sEMGlic (both representing lower intercostal plus costal diaphragm activity) than between Pdi and crural oesEMGdi. Noninvasive indices of NMC derived from sEMGlic and sMMGlic may prove to be useful indices of lower chest wall inspiratory muscle NMC, particularly in settings that do not have access to invasive measures of diaphragm function.

JTD Keywords: Cardiovascular system, Diaphragms, Diseases, Electromyography, Medical signal processing, Neurophysiology, Patient monitoring, Pneumodynamics, Inspiratory muscle neuromechanical coupling, Diaphragm neuromechanical coupling, Neural activation, Transdiaphragmatic pressure, Diaphragm function, Respiratory mechanics, Diaphragm NMC, Invasive measurements, Crural diaphragm electromyography, iNMC, Noninvasive assessment, nNMC, Lower chest wall inspiratory muscles, Inspiratory threshold loading protocol, Noninvasive sMMGlic measurements, sEMGlic measurements, oesEMGdi measurements, Inspiratory threshold load, Lower ribcage muscles, Lower intercostal plus costal diaphragm activity, Crural oesEMGdi, Noninvasive indices, sEMGlic sMMGlic, Lower chest wall inspiratory muscle NMC, Surface mechanomyography, Electromyography, Inspiratory threshold loading, Mechanomyography, Neuromechanical coupling, Respiratory muscles

Torres-Sanchez, A., Millan, D., Arroyo, M., (2019). Modelling fluid deformable surfaces with an emphasis on biological interfaces Journal of Fluid Mechanics 872, 218-271

Fluid deformable surfaces are ubiquitous in cell and tissue biology, including lipid bilayers, the actomyosin cortex or epithelial cell sheets. These interfaces exhibit a complex interplay between elasticity, low Reynolds number interfacial hydrodynamics, chemistry and geometry, and govern important biological processes such as cellular traffic, division, migration or tissue morphogenesis. To address the modelling challenges posed by this class of problems, in which interfacial phenomena tightly interact with the shape and dynamics of the surface, we develop a general continuum mechanics and computational framework for fluid deformable surfaces. The dual solid–fluid nature of fluid deformable surfaces challenges classical Lagrangian or Eulerian descriptions of deforming bodies. Here, we extend the notion of arbitrarily Lagrangian–Eulerian (ALE) formulations, well-established for bulk media, to deforming surfaces. To systematically develop models for fluid deformable surfaces, which consistently treat all couplings between fields and geometry, we follow a nonlinear Onsager formalism according to which the dynamics minimizes a Rayleighian functional where dissipation, power input and energy release rate compete. Finally, we propose new computational methods, which build on Onsager’s formalism and our ALE formulation, to deal with the resulting stiff system of higher-order partial differential equations. We apply our theoretical and computational methodology to classical models for lipid bilayers and the cell cortex. The methods developed here allow us to formulate/simulate these models in their full three-dimensional generality, accounting for finite curvatures and finite shape changes.

JTD Keywords: Capsule/cell dynamics, Computational methods, Membranes

Ruiz, A. D., Mejía, J. S., López, J. M., Giraldo, B. F., (2019). Characterization of cardiac and respiratory system of healthy subjects in supine and sitting position Pattern Recognition and Image Analysis ibPRIA 2019: Iberian Conference on Pattern Recognition and Image Analysis (Lecture Notes in Computer Science) , Springer, Cham (Madrid, Spain) 11867, 367-377

Studies based on the cardiac and respiratory system have allowed a better knowledge of their behavior to contribute with the diagnosis and treatment of diseases associated with them. The main goal of this project was to analyze the behavior of the cardiorespiratory system in healthy subjects, depending on the body position. The electrocardiography and respiratory flow signals were recorded in two positions, supine and sitting. Each signal was analyzed considering sliding windows of 30 s, with and overlapping of 50%. Temporal and spectral features were extracted from each signal. A total of 187 features were extracted for each window. According to statistical analysis, 148 features showed significant differences when comparing the position of the subject. Afterwards, the classifications methods based on decision trees, k-nearest neighbor and support vector machines were applied to identify the best classification model. The most advantageous performance model was obtained with a linear support vector machine method, with an accuracy of 99.5%, a sensitivity of 99.2% and a specificity of 99.6%. In conclusion, we have observed that the position of the body (supine or sitting) could modulate the cardiac and respiratory system response. New statistical models might provide new tools to analyze the behavior of these systems and the cardiorespiratory interaction complexity.

JTD Keywords: Cardiac dynamics, Respiratory dynamics, Statistical models, Supine and sitting posture

Hernández-Vega, Amayra, Marsal, María, Pouille, Philippe-Alexandre, Tosi, Sébastien, Colombelli, Julien, Luque, Tomás, Navajas, Daniel, Pagonabarraga, Ignacio, Martín-Blanco, Enrique, (2017). Polarized cortical tension drives zebrafish epiboly movements EMBO Journal 36, (1), 25-41

The principles underlying the biomechanics of morphogenesis are largely unknown. Epiboly is an essential embryonic event in which three tissues coordinate to direct the expansion of the blastoderm. How and where forces are generated during epiboly, and how these are globally coupled remains elusive. Here we developed a method, hydrodynamic regression (HR), to infer 3D pressure fields, mechanical power, and cortical surface tension profiles. HR is based on velocity measurements retrieved from 2D+T microscopy and their hydrodynamic modeling. We applied HR to identify biomechanically active structures and changes in cortex local tension during epiboly in zebrafish. Based on our results, we propose a novel physical description for epiboly, where tissue movements are directed by a polarized gradient of cortical tension. We found that this gradient relies on local contractile forces at the cortex, differences in elastic properties between cortex components and the passive transmission of forces within the yolk cell. All in all, our work identifies a novel way to physically regulate concerted cellular movements that might be instrumental for the mechanical control of many morphogenetic processes.

JTD Keywords: Epiboly, Hydrodynamics, Mechanics, Morphogenesis, Zebrafish

A. R. Dalton, J., Lans, I., Rovira, X., Malhaire, F., Gómez-Santacana, X., Pittolo, S., Gorostiza, P., Llebaria, A., Goudet, C., Pin, J-P., Giraldo, J., (2016). Shining light on an mGlu5 photoswitchable NAM: A theoretical perspective Current Neuropharmacology , 14, (5), 441-454

Metabotropic glutamate receptors (mGluRs) are important drug targets because of their involvement in several neurological diseases. Among mGluRs, mGlu5 is a particularly high-profile target because its positive or negative allosteric modulation can potentially treat schizophrenia or anxiety and chronic pain, respectively. Here, we computationally and experimentally probe the functional binding of a novel photoswitchable mGlu5 NAM, termed alloswitch-1, which loses its NAM functionality under violet light. We show alloswitch-1 binds deep in the allosteric pocket in a similar fashion to mavoglurant, the co-crystallized NAM in the mGlu5 transmembrane domain crystal structure. Alloswitch-1, like NAM 2-Methyl-6-(phenylethynyl)pyridine (MPEP), is significantly affected by P655M mutation deep in the allosteric pocket, eradicating its functionality. In MD simulations, we show alloswitch-1 and MPEP stabilize the co-crystallized water molecule located at the bottom of the allosteric site that is seemingly characteristic of the inactive receptor state. Furthermore, both NAMs form H-bonds with S809 on helix 7, which may constitute an important stabilizing interaction for NAM-induced mGlu5 inactivation. Alloswitch-1, through isomerization of its amide group from trans to cis is able to form an additional interaction with N747 on helix 5. This may be an important interaction for amide-containing mGlu5 NAMs, helping to stabilize their binding in a potentially unusual cis-amide state. Simulated conformational switching of alloswitch-1 in silico suggests photoisomerization of its azo group from trans to cis may be possible within the allosteric pocket. However, photoexcited alloswitch-1 binds in an unstable fashion, breaking H-bonds with the protein and destabilizing the co-crystallized water molecule. This suggests photoswitching may have destabilizing effects on mGlu5 binding and functionality.

JTD Keywords: Allosteric modulation, Docking, Metabotropic glutamate receptor, Molecular dynamics, Mutation, Protein structure, Transmembrane domain

Garde, Ainara, Voss, Andreas, Caminal, Pere, Benito, Salvador, Giraldo, Beatriz F., (2013). SVM-based feature selection to optimize sensitivity-specificity balance applied to weaning Computers in Biology and Medicine , 43, (5), 533-540

Classification algorithms with unbalanced datasets tend to produce high predictive accuracy over the majority class, but poor predictive accuracy over the minority class. This problem is very common in biomedical data mining. This paper introduces a Support Vector Machine (SVM)-based optimized feature selection method, to select the most relevant features and maintain an accurate and well-balanced sensitivity–specificity result between unbalanced groups. A new metric called the balance index (B) is defined to implement this optimization. The balance index measures the difference between the misclassified data within each class. The proposed optimized feature selection is applied to the classification of patients' weaning trials from mechanical ventilation: patients with successful trials who were able to maintain spontaneous breathing after 48 h and patients who failed to maintain spontaneous breathing and were reconnected to mechanical ventilation after 30 min. Patients are characterized through cardiac and respiratory signals, applying joint symbolic dynamic (JSD) analysis to cardiac interbeat and breath durations. First, the most suitable parameters (C+,C−,σ) are selected to define the appropriate SVM. Then, the feature selection process is carried out with this SVM, to maintain B lower than 40%. The best result is obtained using 6 features with an accuracy of 80%, a B of 18.64%, a sensitivity of 74.36% and a specificity of 82.42%.

JTD Keywords: Support vector machines, Balance index, Sensitivity-specificity balance, Cardiorespiratory interaction, Joint symbolic dynamics, Feature selection, Weaning procedure

Giraldo, B. F., Tellez, J. P., Herrera, S., Benito, S., (2013). Analysis of heart rate variability in elderly patients with chronic heart failure during periodic breathing CinC 2013 Computing in Cardiology Conference (CinC) , IEEE (Zaragoza, Spain) , 991-994

Assessment of the dynamic interactions between cardiovascular signals can provide valuable information that improves the understanding of cardiovascular control. Heart rate variability (HRV) analysis is known to provide information about the autonomic heart rate modulation mechanism. Using the HRV signal, we aimed to obtain parameters for classifying patients with and without chronic heart failure (CHF), and with periodic breathing (PB), non-periodic breathing (nPB), and Cheyne-Stokes respiration (CSR) patterns. An electrocardiogram (ECG) and a respiratory flow signal were recorded in 36 elderly patients: 18 patients with CHF and 18 patients without CHF. According to the clinical criteria, the patients were classified into the follow groups: 19 patients with nPB pattern, 7 with PB pattern, 4 with Cheyne-Stokes respiration (CSR), and 6 non-classified patients (problems with respiratory signal). From the HRV signal, parameters in the time and frequency domain were calculated. Frequency domain parameters were the most discriminant in comparisons of patients with and without CHF: PTot (p = 0.02), PLF (p = 0.022) and fpHF (p = 0.021). For the comparison of the nPB vs. CSR patients groups, the best parameters were RMSSD (p = 0.028) and SDSD (p = 0.028). Therefore, the parameters appear to be suitable for enhanced diagnosis of decompensated CHF patients and the possibility of developed periodic breathing and a CSR pattern.

JTD Keywords: cardiovascular system, diseases, electrocardiography, frequency-domain analysis, geriatrics, medical signal processing, patient diagnosis, pneumodynamics, signal classification, Cheyne-Stokes respiration patterns, ECG, autonomic heart rate modulation mechanism, cardiovascular control, cardiovascular signals, chronic heart failure, decompensated CHF patients, dynamic interaction assessment, elderly patients, electrocardiogram, enhanced diagnosis, frequency domain parameters, heart rate variability analysis, patient classification, periodic breathing, respiratory flow signal recording, Electrocardiography, Frequency modulation, Frequency-domain analysis, Heart rate variability, Senior citizens, Standards

Giraldo, B. F., Chaparro, J. A., Caminal, P., Benito, S., (2013). Characterization of the respiratory pattern variability of patients with different pressure support levels Engineering in Medicine and Biology Society (EMBC) 35th Annual International Conference of the IEEE , IEEE (Osaka, Japan) , 3849-3852

One of the most challenging problems in intensive care is still the process of discontinuing mechanical ventilation, called weaning process. Both an unnecessary delay in the discontinuation process and a weaning trial that is undertaken too early are undesirable. In this study, we analyzed respiratory pattern variability using the respiratory volume signal of patients submitted to two different levels of pressure support ventilation (PSV), prior to withdrawal of the mechanical ventilation. In order to characterize the respiratory pattern, we analyzed the following time series: inspiratory time, expiratory time, breath duration, tidal volume, fractional inspiratory time, mean inspiratory flow and rapid shallow breathing. Several autoregressive modeling techniques were considered: autoregressive models (AR), autoregressive moving average models (ARMA), and autoregressive models with exogenous input (ARX). The following classification methods were used: logistic regression (LR), linear discriminant analysis (LDA) and support vector machines (SVM). 20 patients on weaning trials from mechanical ventilation were analyzed. The patients, submitted to two different levels of PSV, were classified as low PSV and high PSV. The variability of the respiratory patterns of these patients were analyzed. The most relevant parameters were extracted using the classifiers methods. The best results were obtained with the interquartile range and the final prediction errors of AR, ARMA and ARX models. An accuracy of 95% (93% sensitivity and 90% specificity) was obtained when the interquartile range of the expiratory time and the breath duration time series were used a LDA model. All classifiers showed a good compromise between sensitivity and specificity.

JTD Keywords: autoregressive moving average processes, feature extraction, medical signal processing, patient care, pneumodynamics, signal classification, support vector machines, time series, ARX, autoregressive modeling techniques, autoregressive models with exogenous input, autoregressive moving average model, breath duration time series, classification method, classifier method, discontinuing mechanical ventilation, expiratory time, feature extraction, final prediction errors, fractional inspiratory time, intensive care, interquartile range, linear discriminant analysis, logistic regression analysis, mean inspiratory flow, patient respiratory volume signal, pressure support level, pressure support ventilation, rapid shallow breathing, respiratory pattern variability characterization, support vector machines, tidal volume, weaning trial, Analytical models, Autoregressive processes, Biological system modeling, Estimation, Support vector machines, Time series analysis, Ventilation

Hernando, D., Alcaine, A., Pueyo, E., Laguna, P., Orini, M., Arcentales, A., Giraldo, B., Voss, A., Bayes-Genis, A., Bailon, R., (2013). Influence of respiration in the very low frequency modulation of QRS slopes and heart rate variability in cardiomyopathy patients CinC 2013 Computing in Cardiology Conference (CinC) , IEEE (Zaragoza, Spain) , 117-120

This work investigates the very low frequency (VLF) modulation of QRS slopes and heart rate variability (HRV). Electrocardiogram (ECG) and respiratory flow signal were acquired from patients with dilated cardiomyopathy and ischemic cardiomyopathy. HRV as well as the upward QRS slope (IUS) and downward QRS slope (IDS) were extracted from the ECG. The relation between HRV and QRS slopes in the VLF band was measured using ordinary coherence in 5-minute segments. Partial coherence was then used to remove the influence that respiration simultaneously exerts on HRV and QRS slopes. A statistical threshold was determined, below which coherence values were considered not to represent a linear relation. 7 out of 276 segments belonging to 5 out of 29 patients for IUS and 10 segments belonging to 5 patients for IDS presented a VLF modulation in QRS slopes, HRV and respiration. In these segments spectral coherence was statistically significant, while partial coherence decreased, indicating that the coupling HRV and QRS slopes was related to respiration. 4 segments had a partial coherence value below the threshold for IUS, 3 segments for IDS. The rest of the segments also presented a notable decrease in partial coherence, but still above the threshold, which means that other non-linearly effects may also affect this modulation.

JTD Keywords: diseases, electrocardiography, feature extraction, medical signal processing, pneumodynamics, statistical analysis, ECG, QRS slopes, cardiomyopathy patients, dilated cardiomyopathy, electrocardiogram, feature extraction, heart rate variability, ischemic cardiomyopathy, ordinary coherence, partial coherence value, respiration, respiratory flow signal acquisition, spectral coherence, statistical threshold, time 5 min, very low frequency modulation, Coherence, Educational institutions, Electrocardiography, Frequency modulation, Heart rate variability

Gonzalez, H., Acevedo, H., Arizmendi, C., Giraldo, B. F., (2013). Methodology for determine the moment of disconnection of patients of the mechanical ventilation using discrete wavelet transform Complex Medical Engineering (CME) 2013 ICME International Conference , IEEE (Beijing, China) , 483-486

The process of weaning from mechanical ventilation is one of the challenges in intensive care units. 66 patients under extubation process (T-tube test) were studied: 33 patients with successful trials and 33 patients who failed to maintain spontaneous breathing and were reconnected. Each patient was characterized using 7 time series from respiratory signals, and for each serie was evaluated the discrete wavelet transform. It trains a neural network for discriminating between patients from the two groups.

JTD Keywords: discrete wavelet transforms, neural nets, patient treatment, pneumodynamics, time series, ventilation, T-tube test, discrete wavelet transform, extubation process, intensive care units, mechanical ventilation, moment of disconnection, neural network, patients, respiratory signals, spontaneous breathing, time series, weaning, Mechanical Ventilation, Neural Networks, Time series from respiratory signals, Wavelet Transform

Jané, R., Lazaro, J., Ruiz, P., Gil, E., Navajas, D., Farre, R., Laguna, P., (2013). Obstructive Sleep Apnea in a rat model: Effects of anesthesia on autonomic evaluation from heart rate variability measures CinC 2013 Computing in Cardiology Conference (CinC) , IEEE (Zaragoza, Spain) , 1011-1014

Rat model of Obstructive Sleep Apnea (OSA) is a realistic approach for studying physiological mechanisms involved in sleep. Rats are usually anesthetized and autonomic nervous system (ANS) could be blocked. This study aimed to assess the effect of anesthesia on ANS activity during OSA episodes. Seven male Sprague-Dawley rats were anesthetized intraperitoneally with urethane (1g/kg). The experiments were conducted applying airway obstructions, simulating 15s-apnea episodes for 15 minutes. Five signals were acquired: respiratory pressure and flow, SaO2, ECG and photoplethysmography (PPG). In total, 210 apnea episodes were studied. Normalized power spectrum of Pulse Rate Variability (PRV) was analyzed in the Low Frequency (LF) and High Frequency (HF) bands, for each episode in consecutive 15s intervals (before, during and after the apnea). All episodes showed changes in respiratory flow and SaO2 signal. Conversely, decreases in the amplitude fluctuations of PPG (DAP) were not observed. Normalized LF presented extremely low values during breathing (median=7,67%), suggesting inhibition of sympathetic system due to anesthetic effect. Subtle increases of LF were observed during apnea. HRV and PPG analysis during apnea could be an indirect tool to assess the effect and deep of anesthesia.

JTD Keywords: electrocardiography, fluctuations, medical disorders, medical signal detection, medical signal processing, neurophysiology, photoplethysmography, pneumodynamics, sleep, ECG, SaO2 flow, SaO2 signal, airway obstructions, amplitude fluctuations, anesthesia effects, anesthetized nervous system, autonomic evaluation, autonomic nervous system, breathing, heart rate variability, high-frequency bands, low-frequency bands, male Sprague-Dawley rats, normalized power spectrum, obstructive sleep apnea, photoplethysmography, physiological mechanisms, pulse rate variability, rat model, respiratory flow, respiratory pressure, signal acquisition, sympathetic system inhibition, time 15 min, time 15 s, Abstracts, Atmospheric modeling, Computational modeling, Electrocardiography, Rats, Resonant frequency

Giraldo, B. F., Tellez, J. P., Herrera, S., Benito, S., (2013). Study of the oscillatory breathing pattern in elderly patients Engineering in Medicine and Biology Society (EMBC) 35th Annual International Conference of the IEEE , IEEE (Osaka, Japan) , 5228-5231

Some of the most common clinical problems in elderly patients are related to diseases of the cardiac and respiratory systems. Elderly patients often have altered breathing patterns, such as periodic breathing (PB) and Cheyne-Stokes respiration (CSR), which may coincide with chronic heart failure. In this study, we used the envelope of the respiratory flow signal to characterize respiratory patterns in elderly patients. To study different breathing patterns in the same patient, the signals were segmented into windows of 5 min. In oscillatory breathing patterns, frequency and time-frequency parameters that characterize the discriminant band were evaluated to identify periodic and non-periodic breathing (PB and nPB). In order to evaluate the accuracy of this characterization, we used a feature selection process, followed by linear discriminant analysis. 22 elderly patients (7 patients with PB and 15 with nPB pattern) were studied. The following classification problems were analyzed: patients with either PB (with and without apnea) or nPB patterns, and patients with CSR versus PB, CSR versus nPB and PB versus nPB patterns. The results showed 81.8% accuracy in the comparisons of nPB and PB patients, using the power of the modulation peak. For the segmented signal, the power of the modulation peak, the frequency variability and the interquartile ranges provided the best results with 84.8% accuracy, for classifying nPB and PB patients.

JTD Keywords: cardiovascular system, diseases, feature extraction, geriatrics, medical signal processing, oscillations, pneumodynamics, signal classification, time-frequency analysis, Cheyne-Stokes respiration, apnea, cardiac systems, chronic heart failure, classification problems, discriminant band, diseases, elderly patients, feature selection process, frequency variability, interquartile ranges, linear discriminant analysis, nonperiodic breathing, oscillatory breathing pattern, periodic breathing, respiratory How signal, respiratory systems, signal segmentation, time 5 min, time-frequency parameters, Accuracy, Aging, Frequency modulation, Heart, Senior citizens, Time-frequency analysis

Arimon, M., Sanz, F., Giralt, E., Carulla, N., (2012). Template-assisted lateral growth of amyloid-β42 fibrils studied by differential labeling with gold nanoparticles Bioconjugate Chemistry , 23, (1), 27-32

Amyloid-β protein (Aβ) aggregation into amyloid fibrils is central to the origin and development of Alzheimer’s disease (AD), yet this highly complex process is poorly understood at the molecular level. Extensive studies have shown that Aβ fibril growth occurs through fibril elongation, whereby soluble molecules add to the fibril ends. Nevertheless, fibril morphology strongly depends on aggregation conditions. For example, at high ionic strength, Aβ fibrils laterally associate into bundles. To further study the mechanisms leading to fibril growth, we developed a single-fibril growth assay based on differential labeling of two Aβ42 variants with gold nanoparticles. We used this assay to study Aβ42 fibril growth under different conditions and observed that bundle formation is preceded by lateral interaction of soluble Aβ42 molecules with pre-existing fibrils. Based on this data, we propose template-assisted lateral fibril growth as an additional mechanism to elongation for Aβ42 fibril growth.

JTD Keywords: AFM, Beta-Amyloid Fibrils, Polymorphism, Association, Elongation, Dynamics, State

Valle-Delgado, J. J., Liepina, I., Lapidus, D., Sabaté, R., Ventura, S., Samitier, J., Fernàndez-Busquets, X., (2012). Self-assembly of human amylin-derived peptides studied by atomic force microscopy and single molecule force spectroscopy Soft Matter 8, (4), 1234-1242

The self-assembly of peptides and proteins into amyloid fibrils of nanometric thickness and up to several micrometres in length, a phenomenon widely observed in biological systems, has recently aroused a growing interest in nanotechnology and nanomedicine. Here we have applied atomic force microscopy and single molecule force spectroscopy to study the amyloidogenesis of a peptide derived from human amylin and of its reverse sequence. The spontaneous formation of protofibrils and their orientation along well-defined directions on graphite and DMSO-coated graphite substrates make the studied peptides interesting candidates for nanotechnological applications. The measured binding forces between peptides correlate with the number of hydrogen bonds between individual peptides inside the fibril structure according to molecular dynamics simulations.

JTD Keywords: Amyloid fibril, Amyloidogenesis, Binding forces, Fibril structure, Graphite substrate, Molecular dynamics simulations, Nanometrics, Protofibrils, Single molecule force spectroscopy, Spontaneous formation, Atomic force microscopy, Atomic spectroscopy, Graphite, Hydrogen bonds, Medical nanotechnology, Molecular dynamics, Molecular physics, Self assembly, Thickness measurement, Peptides

Redondo-Morata, Lorena, Oncins, Gerard, Sanz, Fausto, (2012). Force spectroscopy reveals the effect of different ions in the nanomechanical behavior of phospholipid model membranes: The case of potassium cation Biophysical Journal , 102, (1), 66-74

How do metal cations affect the stability and structure of phospholipid bilayers? What role does ion binding play in the insertion of proteins and the overall mechanical stability of biological membranes? Investigators have used different theoretical and microscopic approaches to study the mechanical properties of lipid bilayers. Although they are crucial for such studies, molecular-dynamics simulations cannot yet span the complexity of biological membranes. In addition, there are still some experimental difficulties when it comes to testing the ion binding to lipid bilayers in an accurate way. Hence, there is a need to establish a new approach from the perspective of the nanometric scale, where most of the specific molecular phenomena take place. Atomic force microscopy has become an essential tool for examining the structure and behavior of lipid bilayers. In this work, we used force spectroscopy to quantitatively characterize nanomechanical resistance as a function of the electrolyte composition by means of a reliable molecular fingerprint that reveals itself as a repetitive jump in the approaching force curve. By systematically probing a set of bilayers of different composition immersed in electrolytes composed of a variety of monovalent and divalent metal cations, we were able to obtain a wealth of information showing that each ion makes an independent and important contribution to the gross mechanical resistance and its plastic properties. This work addresses the need to assess the effects of different ions on the structure of phospholipid membranes, and opens new avenues for characterizing the (nano)mechanical stability of membranes.

JTD Keywords: Molecular-dynamics simulation, Liquid expanded monolayers, Lipid-bilayers, Hofmeister series, Monovalent salt, Phosphatidylcholine, Microscopy, Binding, Surfaces, NaCl

Giraldo, B.F., Gaspar, B.W., Caminal, P., Benito, S., (2012). Analysis of roots in ARMA model for the classification of patients on weaning trials Engineering in Medicine and Biology Society (EMBC) 34th Annual International Conference of the IEEE , IEEE (San Diego, USA) , 698-701

One objective of mechanical ventilation is the recovery of spontaneous breathing as soon as possible. Remove the mechanical ventilation is sometimes more difficult that maintain it. This paper proposes the study of respiratory flow signal of patients on weaning trials process by autoregressive moving average model (ARMA), through the location of poles and zeros of the model. A total of 151 patients under extubation process (T-tube test) were analyzed: 91 patients with successful weaning (GS), 39 patients that failed to maintain spontaneous breathing and were reconnected (GF), and 21 patients extubated after the test but before 48 hours were reintubated (GR). The optimal model was obtained with order 8, and statistical significant differences were obtained considering the values of angles of the first four poles and the first zero. The best classification was obtained between GF and GR, with an accuracy of 75.3% on the mean value of the angle of the first pole.

JTD Keywords: Analytical models, Biological system modeling, Computational modeling, Estimation, Hospitals, Poles and zeros, Ventilation, Autoregressive moving average processes, Patient care, Patient monitoring, Pneumodynamics, Poles and zeros, Ventilation, ARMA model, T-tube test, Autoregressive moving average model, Extubation process, Mechanical ventilation, Optimal model, Patient classification, Respiratory flow signal, Roots, Spontaneous breathing, Weaning trials

Chaparro, J.A., Giraldo, B.F., Caminal, P., Benito, S., (2012). Performance of respiratory pattern parameters in classifiers for predict weaning process Engineering in Medicine and Biology Society (EMBC) 34th Annual International Conference of the IEEE , IEEE (San Diego, USA) , 4349-4352

Weaning trials process of patients in intensive care units is a complex clinical procedure. 153 patients under extubation process (T-tube test) were studied: 94 patients with successful trials (group S), 38 patients who failed to maintain spontaneous breathing and were reconnected (group F), and 21 patients with successful test but that had to be reintubated before 48 hours (group R). The respiratory pattern of each patient was characterized through the following time series: inspiratory time (TI), expiratory time (TE), breathing cycle duration (TTot), tidal volume (VT), inspiratory fraction (TI/TTot), half inspired flow (VT/TI), and rapid shallow index (f/VT), where f is respiratory rate. Using techniques as autoregressive models (AR), autoregressive moving average models (ARMA) and autoregressive models with exogenous input (ARX), the most relevant parameters of the respiratory pattern were obtained. We proposed the evaluation of these parameters using classifiers as logistic regression (LR), linear discriminant analysis (LDA), support vector machines (SVM) and classification and regression tree (CART) to discriminate between patients from groups S, F and R. An accuracy of 93% (98% sensitivity and 82% specificity) has been obtained using CART classification.

JTD Keywords: Accuracy, Indexes, Logistics, Regression tree analysis, Support vector machines, Time series analysis, Autoregressive moving average processes, Medical signal processing, Pattern classification, Pneumodynamics, Regression analysis, Sensitivity, Signal classification, Support vector machines, Time series, SVM, T-tube testing, Autoregressive models-with-exogenous input, Autoregressive moving average models, Breathing cycle duration, Classification-and-regression tree, Expiratory time, Extubation process, Half inspired flow, Inspiratory fraction, Inspiratory time, Intensive care units, Linear discriminant analysis, Logistic regression, Rapid shallow index, Respiratory pattern parameter performance, Sensitivity, Spontaneous breathing, Support vector machines, Tidal volume, Time 48 hr, Time series, Weaning process classifiers

Garde, A., Giraldo, B.F., Jané, R., Latshang, T.D., Turk, A.J., Hess, T., Bosch, M-.M., Barthelmes, D., Hefti, J.P., Maggiorini, M., Hefti, U., Merz, T.M., Schoch, O.D., Bloch, K.E., (2012). Periodic breathing during ascent to extreme altitude quantified by spectral analysis of the respiratory volume signal Engineering in Medicine and Biology Society (EMBC) 34th Annual International Conference of the IEEE , IEEE (San Diego, USA) , 707-710

High altitude periodic breathing (PB) shares some common pathophysiologic aspects with sleep apnea, Cheyne-Stokes respiration and PB in heart failure patients. Methods that allow quantifying instabilities of respiratory control provide valuable insights in physiologic mechanisms and help to identify therapeutic targets. Under the hypothesis that high altitude PB appears even during physical activity and can be identified in comparison to visual analysis in conditions of low SNR, this study aims to identify PB by characterizing the respiratory pattern through the respiratory volume signal. A number of spectral parameters are extracted from the power spectral density (PSD) of the volume signal, derived from respiratory inductive plethysmography and evaluated through a linear discriminant analysis. A dataset of 34 healthy mountaineers ascending to Mt. Muztagh Ata, China (7,546 m) visually labeled as PB and non periodic breathing (nPB) is analyzed. All climbing periods within all the ascents are considered (total climbing periods: 371 nPB and 40 PB). The best crossvalidated result classifying PB and nPB is obtained with Pm (power of the modulation frequency band) and R (ratio between modulation and respiration power) with an accuracy of 80.3% and area under the receiver operating characteristic curve of 84.5%. Comparing the subjects from 1st and 2nd ascents (at the same altitudes but the latter more acclimatized) the effect of acclimatization is evaluated. SaO2 and periodic breathing cycles significantly increased with acclimatization (p-value <; 0.05). Higher Pm and higher respiratory frequencies are observed at lower SaO2, through a significant negative correlation (p-value <; 0.01). Higher Pm is observed at climbing periods visually labeled as PB with >; 5 periodic breathing cycles through a significant positive correlation (p-value <; 0.01). Our data demonstrate that quantification of the respiratory volum- signal using spectral analysis is suitable to identify effects of hypobaric hypoxia on control of breathing.

JTD Keywords: Frequency domain analysis, Frequency modulation, Heart, Sleep apnea, Ventilation, Visualization, Cardiology, Medical disorders, Medical signal processing, Plethysmography, Pneumodynamics, Sensitivity analysis, Sleep, Spectral analysis, Cheyne-Stokes respiration, Climbing periods, Dataset, Heart failure patients, High altitude PB, High altitude periodic breathing, Hypobaric hypoxia, Linear discriminant analysis, Pathophysiologic aspects, Physical activity, Physiologic mechanisms, Power spectral density, Receiver operating characteristic curve, Respiratory control, Respiratory frequency, Respiratory inductive plethysmography, Respiratory pattern, Respiratory volume signal, Sleep apnea, Spectral analysis, Spectral parameters

Mesquita, J., Poree, F., Carrault, G., Fiz, J. A., Abad, J., Jané, R., (2012). Respiratory and spontaneous arousals in patients with Sleep Apnea Hypopnea Syndrome Engineering in Medicine and Biology Society (EMBC) 34th Annual International Conference of the IEEE , IEEE (San Diego, USA) , 6337-6340

Sleep in patients with Sleep Apnea-Hypopnea Syndrome (SAHS) is frequently interrupted with arousals. Increased amounts of arousals result in shortening total sleep time and repeated sleep-arousal change can result in sleep fragmentation. According to the American Sleep Disorders Association (ASDA) an arousal is a marker of sleep disruption representing a detrimental and harmful feature for sleep. The nature of arousals and its role on the regulation of the sleep process raises controversy and has sparked the debate in the last years. In this work, we analyzed and compared the EEG spectral content of respiratory and spontaneous arousals on a database of 45 SAHS subjects. A total of 3980 arousals (1996 respiratory and 1984 spontaneous) were analyzed. The results showed no differences between the spectral content of the two kinds of arousals. Our findings raise doubt as to whether these two kinds of arousals are truly triggered by different organic mechanisms. Furthermore, they may also challenge the current beliefs regarding the underestimation of the importance of spontaneous arousals and their contribution to sleep fragmentation in patients suffering from SAHS.

JTD Keywords: Adaptive filters, Correlation, Databases, Electroencephalography, Hospitals, Sleep apnea, Electroencephalography, Medical signal processing, Pneumodynamics, Sleep, EEG spectral content, Organic mechanism, Respiratory, Sleep apnea hypopnea syndrome, Sleep fragmentation, Spectral content, Spontaneous arousal

Angelini, Thomas E., Hannezo, Edouard, Trepat, Xavier, Marquez, Manuel, Fredberg, Jeffrey J., Weitz, David A., (2011). Glass-like dynamics of collective cell migration Proceedings of the National Academy of Sciences of the United States of America 108, (12), 4714-4719

Collective cell migration in tissues occurs throughout embryonic development, during wound healing, and in cancerous tumor invasion, yet most detailed knowledge of cell migration comes from single-cell studies. As single cells migrate, the shape of the cell body fluctuates dramatically through cyclic processes of extension, adhesion, and retraction, accompanied by erratic changes in migration direction. Within confluent cell layers, such subcellular motions must be coupled between neighbors, yet the influence of these subcellular motions on collective migration is not known. Here we study motion within a confluent epithelial cell sheet, simultaneously measuring collective migration and subcellular motions, covering a broad range of length scales, time scales, and cell densities. At large length scales and time scales collective migration slows as cell density rises, yet the fastest cells move in large, multicell groups whose scale grows with increasing cell density. This behavior has an intriguing analogy to dynamic heterogeneities found in particulate systems as they become more crowded and approach a glass transition. In addition we find a diminishing self-diffusivity of short-wavelength motions within the cell layer, and growing peaks in the vibrational density of states associated with cooperative cell-shape fluctuations. Both of these observations are also intriguingly reminiscent of a glass transition. Thus, these results provide a broad and suggestive analogy between cell motion within a confluent layer and the dynamics of supercooled colloidal and molecular fluids approaching a glass transition.

JTD Keywords: Active matter, Cell mechanics, Jamming, Collective cell dynamics, Nonequilibrium

Melchels, Ferry P. W., Tonnarelli, Beatrice, Olivares, Andy L., Martin, Ivan, Lacroix, Damien, Feijen, Jan, Wendt, David J., Grijpma, Dirk W., (2011). The influence of the scaffold design on the distribution of adhering cells after perfusion cell seeding Biomaterials 32, (11), 2878-2884

In natural tissues, the extracellular matrix composition, cell density and physiological properties are often non-homogeneous. Here we describe a model system, in which the distribution of cells throughout tissue engineering scaffolds after perfusion seeding can be influenced by the pore architecture of the scaffold. Two scaffold types, both with gyroid pore architectures, were designed and built by stereolithography: one with isotropic pore size (412 ± 13 [mu]m) and porosity (62 ± 1%), and another with a gradient in pore size (250-500 [mu]m) and porosity (35%-85%). Computational fluid flow modelling showed a uniform distribution of flow velocities and wall shear rates (15-24 s-1) for the isotropic architecture, and a gradient in the distribution of flow velocities and wall shear rates (12-38 s-1) for the other architecture. The distribution of cells throughout perfusion-seeded scaffolds was visualised by confocal microscopy. The highest densities of cells correlated with regions of the scaffolds where the pores were larger, and the fluid velocities and wall shear rates were the highest. Under the applied perfusion conditions, cell deposition is mainly determined by local wall shear stress, which, in turn, is strongly influenced by the architecture of the pore network of the scaffold.

JTD Keywords: Scaffolds, Microstructure, Cell adhesion, Confocal microscopy, Image analysis, Computational fluid dynamics

Llorens, Franc, Hummel, Manuela, Pastor, Xavier, Ferrer, Anna, Pluvinet, Raquel, Vivancos, Ana, Castillo, Ester, Iraola, Susana, Mosquera, Ana M., Gonzalez, Eva, Lozano, Juanjo, Ingham, Matthew, Dohm, Juliane C., Noguera, Marc, Kofler, Robert, Antonio del Rio, Jose, Bayes, Monica, Himmelbauer, Heinz, Sumoy, Lauro, (2011). Multiple platform assessment of the EGF dependent transcriptome by microarray and deep tag sequencing analysis BMC Genomics 12, 326

Background: Epidermal Growth Factor (EGF) is a key regulatory growth factor activating many processes relevant to normal development and disease, affecting cell proliferation and survival. Here we use a combined approach to study the EGF dependent transcriptome of HeLa cells by using multiple long oligonucleotide based microarray platforms (from Agilent, Operon, and Illumina) in combination with digital gene expression profiling (DGE) with the Illumina Genome Analyzer. Results: By applying a procedure for cross-platform data meta-analysis based on RankProd and GlobalAncova tests, we establish a well validated gene set with transcript levels altered after EGF treatment. We use this robust gene list to build higher order networks of gene interaction by interconnecting associated networks, supporting and extending the important role of the EGF signaling pathway in cancer. In addition, we find an entirely new set of genes previously unrelated to the currently accepted EGF associated cellular functions. Conclusions: We propose that the use of global genomic cross-validation derived from high content technologies (microarrays or deep sequencing) can be used to generate more reliable datasets. This approach should help to improve the confidence of downstream in silico functional inference analyses based on high content data.

JTD Keywords: Gene-expression measurements, Quality-control maqc, Cancer-cell-lines, Real-time pcr, Oligonucleotide microarrays, Phosphorylation dynamics, In-vivo, Networks, Signal, Technologies

Garcia-Manyes, S., Redondo-Morata, L., Oncins, G., Sanz, F., (2010). Nanomechanics of lipid bilayers: Heads or tails? Journal of the American Chemical Society American Chemical Society 132, (37), 12874-12886

Understanding the effect of mechanical stress on membranes is of primary importance in biophysics. Here we use force spectroscopy AFM to quantitatively characterize the nanomechanical stability of supported lipid bilayers as a function of their chemical composition. The onset of plastic deformation reveals itself as a repetitive jump in the approaching force curve, which represents a molecular fingerprint for the bilayer mechanical stability. By systematically probing a set of chemically distinct supported lipid bilayers (SLBs), we first show that both the headgroup and tail have a decisive effect on their mechanical properties. While the mechanical stability of the probed SLBs linearly increases by 3.3 nN upon the introduction of each additional -CH2- in the chain, it exhibits a significant dependence on the phospholipid headgroup, ranging from 3 nN for DPPA to 66 nN for DPPG. Furthermore, we also quantify the reduction of the membrane mechanical stability as a function of the number of unsaturations and molecular branching in the chemical structure of the apolar tails. Finally, we demonstrate that, upon introduction of cholesterol and ergosterol, contrary to previous belief the mechanical stability of membranes not only increases linearly in the liquid phase (DLPC) but also for phospholipids present in the gel phase (DPPC). Our results are discussed in the framework of the continuum nucleation model. This work highlights the compelling effect of subtle variations in the chemical structure of phospholipid molecules on the membrane response when exposed to mechanical forces, a mechanism of common occurrence in nature.

JTD Keywords: Atomic-force microscopy, Molecular-dynamics simulation, Aqueous-electrolyte solutions, Supported planar membranes, Phospholipid-bilayers, Biological-membranes, Physical-properties, Fluid membranes, Model membranes, Chain-length

Santoro, R., Olivares, A. L., Brans, G., Wirz, D., Longinotti, C., Lacroix, D., Martin, I., Wendt, D., (2010). Bioreactor based engineering of large-scale human cartilage grafts for joint resurfacing Biomaterials 31, (34), 8946-8952

Apart from partial or total joint replacement, no surgical procedure is currently available to treat large and deep cartilage defects associated with advanced diseases such as osteoarthritis. In this work, we developed a perfusion bioreactor system to engineer human cartilage grafts in a size with clinical relevance for unicompartmental resurfacing of human knee joints (50 mm diameter x 3 mm thick). Computational fluid dynamics models were developed to optimize the flow profile when designing the perfusion chamber. Using the developed system, human chondrocytes could be seeded throughout large 50 mm diameter scaffolds with a uniform distribution. Following two weeks culture, tissues grown in the bioreactor were viable and homogeneously cartilaginous, with biomechanical properties approaching those of native cartilage. In contrast, tissues generated by conventional manual production procedures were highly inhomogeneous and contained large necrotic regions. The unprecedented engineering of human cartilage tissues in this large-scale opens the practical perspective of grafting functional biological substitutes for the clinical treatment for extensive cartilage defects, possibly in combination with surgical or pharmacological therapies to support durability of the implant. Ongoing efforts are aimed at integrating the up-scaled bioreactor based processes within a fully automated and closed manufacturing system for safe, standardized, and GMP compliant production of large-scale cartilage grafts.

JTD Keywords: Bioreactor, Cartilage repair, Computational fluid dynamics, Scale-up, Regenerative medicine, Tissue engineering

Estevez, M., Fernandez-Ulibarri, I., Martinez, E., Egea, G., Samitier, J., (2010). Changes in the internal organization of the cell by microstructured substrates Soft Matter 6, (3), 582-590

Surface features at the micro and nanometre scale have been shown to influence and even determine cell behaviour and cytoskeleton organization through direct mechanotransductive pathways. Much less is known about the function and internal distribution of organelles of cells grown on topographically modified surfaces. In this study, the nanoimprint lithography technique was used to manufacture poly(methyl methacrylate) (PMMA) sheets with a variety of features in the micrometre size range. Normal rat kidney (NRK) fibroblasts were cultured on these substrates and immunofluorescence staining assays were performed to visualize cell adhesion, the organization of the cytoskeleton and the morphology and subcellular positioning of the Golgi complex. The results show that different topographic features at the micrometric scale induce different rearrangements of the cell cytoskeleton, which in turn alter the positioning and morphology of the Golgi complex. Microposts and microholes alter the mechanical stability of the Golgi complex by modifying the actin cytoskeleton organization leading to the compaction of the organelle. These findings prove that physically modified surfaces are a valuable tool with which to study the dynamics of cell cytoskeleton organization and its subsequent repercussion on internal cell organization and associated function.

JTD Keywords: Actin stress fibers, Golgi-complex, Focal adhesions, Cytoskeletal organization, Osteoblast adhesion, Mammalian-cells, Micron-scale, Nanoscale, Dynamics, Rho

Garde, A., Schroeder, R., Voss, A., Caminal, P., Benito, S., Giraldo, B., (2010). Patients on weaning trials classified with support vector machines Physiological Measurement , 31, (7), 979-993

The process of discontinuing mechanical ventilation is called weaning and is one of the most challenging problems in intensive care. An unnecessary delay in the discontinuation process and an early weaning trial are undesirable. This study aims to characterize the respiratory pattern through features that permit the identification of patients' conditions in weaning trials. Three groups of patients have been considered: 94 patients with successful weaning trials, who could maintain spontaneous breathing after 48 h ( GSucc ); 39 patients who failed the weaning trial ( GFail ) and 21 patients who had successful weaning trials, but required reintubation in less than 48 h ( GRein ). Patients are characterized by their cardiorespiratory interactions, which are described by joint symbolic dynamics (JSD) applied to the cardiac interbeat and breath durations. The most discriminating features in the classification of the different groups of patients ( GSucc , GFail and GRein ) are identified by support vector machines (SVMs). The SVM-based feature selection algorithm has an accuracy of 81% in classifying GSucc versus the rest of the patients, 83% in classifying GRein versus GSucc patients and 81% in classifying GRein versus the rest of the patients. Moreover, a good balance between sensitivity and specificity is achieved in all classifications.

JTD Keywords: Mechanical ventilation, Weaning, Support vector machines, Joint symbolic dynamics

Morgenstern, C., Schwaibold, M., Randerath, W., Bolz, A., Jané, R., (2010). Automatic non-invasive differentiation of obstructive and central hypopneas with nasal airflow compared to esophageal pressure Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) , 6142-6145

The differentiation of obstructive and central respiratory events is a major challenge in the diagnosis of sleep disordered breathing. Esophageal pressure (Pes) measurement is the gold-standard method to identify these events but its invasiveness deters its usage in clinical routine. Flattening patterns appear in the airflow signal during episodes of inspiratory flow limitation (IFL) and have been shown with invasive techniques to be useful to differentiate between central and obstructive hypopneas. In this study we present a new method for the automatic non-invasive differentiation of obstructive and central hypopneas solely with nasal airflow. An overall of 36 patients underwent full night polysomnography with systematic Pes recording and a total of 1069 hypopneas were manually scored by human experts to create a gold-standard annotation set. Features were automatically extracted from the nasal airflow signal to train and test our automatic classifier (Discriminant Analysis). Flattening patterns were non-invasively assessed in the airflow signal using spectral and time analysis. The automatic non-invasive classifier obtained a sensitivity of 0.71 and an accuracy of 0.69, similar to the results obtained with a manual non-invasive classification algorithm. Hence, flattening airflow patterns seem promising for the non-invasive differentiation of obstructive and central hypopneas.

JTD Keywords: Practical, Experimental/ biomedical measurement, Feature extraction, Flow measurement, Medical disorders, Medical signal processing, Patient diagnosis, Pneumodynamics, Pressure measurement, Signal classification, Sleep, Spectral analysis/ automatic noninvasive differentiation, Obstructive hypopnea, Central hypopnea, Inspiratory flow limitation, Nasal airflow, Esophageal pressure, Polysomnography, Feature extraction, Discriminant analysis, Spectral analysis

Garde, A., Sörnmo, L., Jané, R., Giraldo, B. F., (2010). Correntropy-based nonlinearity test applied to patients with chronic heart failure Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) , 2399-2402

In this study we propose the correntropy function as a discriminative measure for detecting nonlinearities in the respiratory pattern of chronic heart failure (CHF) patients with periodic or nonperiodic breathing pattern (PB or nPB, respectively). The complexity seems to be reduced in CHF patients with higher risk level. Correntropy reflects information on both, statistical distribution and temporal structure of the underlying dataset. It is a suitable measure due to its capability to preserve nonlinear information. The null hypothesis considered is that the analyzed data is generated by a Gaussian linear stochastic process. Correntropy is used in a statistical test to reject the null hypothesis through surrogate data methods. Various parameters, derived from the correntropy and correntropy spectral density (CSD) to characterize the respiratory pattern, presented no significant differences when extracted from the iteratively refined amplitude adjusted Fourier transform (IAAFT) surrogate data. The ratio between the powers in the modulation and respiratory frequency bands R was significantly different in nPB patients, but not in PB patients, which reflects a higher presence of nonlinearities in nPB patients than in PB patients.

JTD Keywords: Practical, Theoretical or Mathematical, Experimental/cardiology diseases, Fourier transforms, Medical signal processing, Pattern classification, Pneumodynamics, Spectral analysis, Statistical analysis, Stochastic processes/ correntropy based nonlinearity test, Chronic heart failure, Correntropy function, Respiratory pattern nonlinearities, CHF patients, Nonperiodic breathing pattern, Dataset statistical distribution, Dataset temporal structure, Nonlinear information, Null hypothesis, Gaussian linear stochastic process, Statistical test, Correntropy spectral density, Iteratively refined amplitude adjusted Fourier transform, Surrogate data, Periodic breathing pattern

Correa, L. S., Laciar, E., Mut, V., Giraldo, B. F., Torres, A., (2010). Multi-parameter analysis of ECG and Respiratory Flow signals to identify success of patients on weaning trials Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) -----, 6070-6073

Statistical analysis, power spectral density, and Lempel Ziv complexity, are used in a multi-parameter approach to analyze four temporal series obtained from the Electrocardiographic and Respiratory Flow signals of 126 patients on weaning trials. In which, 88 patients belong to successful group (SG), and 38 patients belong to failure group (FG), i.e. failed to maintain spontaneous breathing during trial. It was found that mean values of cardiac inter-beat and breath durations give higher values for SG than for FG; Kurtosis coefficient of the spectrum of the rapid shallow breathing index is higher for FG; also Lempel Ziv complexity mean values associated with the respiratory flow signal are bigger for FG. Patients were then classified with a pattern recognition neural network, obtaining 80% of correct classifications (81.6% for FG and 79.5% for SG).

JTD Keywords: Electrocardiography, Medical signal processing, Neural nets, Pattern recognition, Pneumodynamics, Signal classification, Statistical analysis, ECG, Kurtosis coefficient, Lempel Ziv complexity, Breath durations, Cardiac interbeat durations, Electrocardiography, Multiparameter analysis, Pattern recognition neural network, Power spectral density, Respiratory flow signals, Signal classification, Spontaneous breathing, Statistical analysis, Weaning trials

Leder, R. S., Schlotthauer, G., Penzel, T., Jané, R., (2010). The natural history of the sleep and respiratory engineering track at EMBC 1988 to 2010 Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) , 288-291

Sleep science and respiratory engineering as medical subspecialties and research areas grew up side-by-side with biomedical engineering. The formation of EMBS in the 1950's and the discovery of REM sleep in the 1950's led to parallel development and interaction of sleep and biomedical engineering in diagnostics and therapeutics.

JTD Keywords: Practical/ biomedical equipment, Biomedical measurement, Patient diagnosis, Patient monitoring, Patient treatment, Pneumodynamics, Sleep/ sleep engineering, Respiratory engineering, Automatic sleep analysis, Automatic sleep interpretation systems, Breathing, Biomedical, Engineering, Diagnostics, Therapeutics, REM sleep, Portable, Measurement, Ambulatory measurement, Monitoring

Mesquita, J., Fiz, J. A., Solà, J., Morera, J., Jané, R., (2010). Regular and non regular snore features as markers of SAHS Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) , 6138-6141

Sleep Apnea-Hypopnea Syndrome (SAHS) diagnosis is still done with an overnight multi-channel polysomnography. Several efforts are being made to study profoundly the snore mechanism and discover how it can provide an opportunity to diagnose the disease. This work introduces the concept of regular snores, defined as the ones produced in consecutive respiratory cycles, since they are produced in a regular way, without interruptions. We applied 2 thresholds (TH/sub adaptive/ and TH/sub median/) to the time interval between successive snores of 34 subjects in order to select regular snores from the whole all-night snore sequence. Afterwards, we studied the effectiveness that parameters, such as time interval between successive snores and the mean intensity of snores, have on distinguishing between different levels of SAHS severity (AHI (Apnea-Hypopnea Index)<5h/sup -1/, AHI<10 h/sup -1/, AHI<15h/sup -1/, AHI<30h/sup -1/). Results showed that TH/sub adaptive/ outperformed TH/sub median/ on selecting regular snores. Moreover, the outcome achieved with non-regular snores intensity features suggests that these carry key information on SAHS severity.

JTD Keywords: Practical, Experimental/ acoustic signal processing, Bioacoustics, Biomedical measurement, Diseases, Feature extraction, Medical signal processing, Patient diagnosis, Pneumodynamics, Sleep/ nonregular snore features, SAHS markers, Sleep apnea hypopnea syndrome, Overnight multichannel polysomnography, Snore mechanism

Arcentales, A., Giraldo, B. F., Caminal, P., Diaz, I., Benito, S., (2010). Spectral analysis of the RR series and the respiratory flow signal on patients in weaning process Engineering in Medicine and Biology Society (EMBC) 32nd Annual International Conference of the IEEE , IEEE (Buenos Aires, Argentina) , 2485-2488

A considerable number of patients in weaning process have problems to keep spontaneous breathing during the trial and after it. This study proposes to extract characteristic parameters of the RR series and respiratory flow signal according to the patients' condition in weaning test. Three groups of patients have been considered: 93 patients with successful trials (group S), 40 patients that failed to maintain spontaneous breathing (group F), and 21 patients who had successful weaning trials, but that had to be reintubated before 48 hours (group R). The characterization was performed using spectral analysis of the signals, through the power spectral density, cross power spectral density and Coherence method. The parameters were extracted on the three frequency bands (VLF, LF and HF), and the principal statistical differences between groups were obtained in bands of VLF and HF. The results show an accuracy of 76.9% in the classification of the groups S and F.

JTD Keywords: Biomedical measurement, Electrocardiography, Medical signal processing, Pneumodynamics, Spectral analysis, RR series, Coherence method, Cross power spectral density, Electrocardiography, Principal statistical differences, Respiratory flow signal, Spectral analysis, Spontaneous breathing, Weaning test

Fumagalli, L., Ferrari, G., Sampietro, M., Gomila, G., (2009). Quantitative nanoscale dielectric microscopy of single-layer supported biomembranes Nano Letters 9, (4), 1604-1608

We present the experimental demonstration of low-frequency dielectric constant imaging of single-layer supported biomembranes at the nanoscale. The dielectric constant image has been quantitatively reconstructed by combining the thickness and local capacitance obtained using a scanning force microscope equipped with a sub-attofarad low-frequency capacitance detector. This work opens new possibilities for studying bioelectric phenomena and the dielectric properties of biological membranes at the nanoscale.

JTD Keywords: Atomic-force microscopy, Nnear-field microscopy, Purple membrane, Scanning capacitance, Biological-systems, Fluid, Spectroscopy, Resolution, Proteins, Dynamics

Roca-Cusachs, P., Gauthier, N. C., del Rio, A., Sheetz, M. P., (2009). Clustering of alpha(5)beta(1) integrins determines adhesion strength whereas alpha(v)beta(3) and talin enable mechanotransduction Proceedings of the National Academy of Sciences of the United States of America 106, (38), 16245-16250

A key molecular link between cells and the extracellular matrix is the binding between fibronectin and integrins alpha(5)beta(1) and alpha(v)beta(3). However, the roles of these different integrins in establishing adhesion remain unclear. We tested the adhesion strength of fibronectin-integrin-cytoskeleton linkages by applying physiological nanonewton forces to fibronectin-coated magnetic beads bound to cells. We report that the clustering of fibronectin domains within 40 nm led to integrin alpha(5)beta(1) recruitment, and increased the ability to sustain force by over six-fold. This force was supported by alpha(5)beta(1) integrin clusters. Importantly, we did not detect a role of either integrin alpha(v)beta(3) or talin 1 or 2 in maintaining adhesion strength. Instead, these molecules enabled the connection to the cytoskeleton and reinforcement in response to an applied force. Thus, high matrix forces are primarily supported by clustered alpha(5)beta(1) integrins, while less stable links to alpha(v)beta(3) integrins initiate mechanotransduction, resulting in reinforcement of integrin-cytoskeleton linkages through talin-dependent bonds.

JTD Keywords: Cell-adhesion, Mechanical force, Vinculin-binding, Fibronectin, Activation, Dynamics, Domain, Alpha-v-beta-3, Translocation, Bonds

Milan, J. L., Planell, J. A., Lacroix, D., (2009). Computational modelling of the mechanical environment of osteogenesis within a polylactic acid-calcium phosphate glass scaffold Biomaterials 30, (25), 4219-4226

A computational model based on finite element method (FEM) and computational fluid dynamics (CFD) is developed to analyse the mechanical stimuli in a composite scaffold made of polylactic acid (PLA) matrix with calcium phosphate glass (Glass) particles. Different bioreactor loading conditions were simulated within the scaffold. In vitro perfusion conditions were reproduced in the model. Dynamic compression was also reproduced in an uncoupled fluid-structure scheme: deformation level was studied analyzing the mechanical response of scaffold alone under static compression while strain rate was studied considering the fluid flow induced by compression through fixed scaffold. Results of the model show that during perfusion test an inlet velocity of 25mum/s generates on scaffold surface a fluid flow shear stress which may stimulate osteogenesis. Dynamic compression of 5% applied on the PLA-Glass scaffold with a strain rate of 0.005s(-1) has the benefit to generate mechanical stimuli based on both solid shear strain and fluid flow shear stress on large scaffold surface area. Values of perfusion inlet velocity or compression strain rate one order of magnitude lower may promote cell proliferation while values one order of magnitude higher may be detrimental for cells. FEM-CFD scaffold models may help to determine loading conditions promoting bone formation and to interpret experimental results from a mechanical point of view.

JTD Keywords: Bone tissue engineering, Scaffold, Finite element analysis, Computational fluid dynamics, Mechanical stimuli

Olivares, A. L., Marshal, E., Planell, J. A., Lacroix, D., (2009). Finite element study of scaffold architecture design and culture conditions for tissue engineering Biomaterials 30, (30), 6142-6149

Tissue engineering scaffolds provide temporary mechanical support for tissue regeneration and transfer global mechanical load to mechanical stimuli to cells through its architecture. In this study the interactions between scaffold pore morphology, mechanical stimuli developed at the cell microscopic level, and culture conditions applied at the macroscopic scale are studied on two regular scaffold structures. Gyroid and hexagonal scaffolds of 55% and 70% porosity were modeled in a finite element analysis and were submitted to an inlet fluid flow or compressive strain. A mechanoregulation theory based on scaffold shear strain and fluid shear stress was applied for determining the influence of each structures on the mechanical stimuli on initial conditions. Results indicate that the distribution of shear stress induced by fluid perfusion is very dependent on pore distribution within the scaffold. Gyroid architectures provide a better accessibility of the fluid than hexagonal structures. Based on the mechanoregulation theory, the differentiation process in these structures was more sensitive to inlet fluid flow than axial strain of the scaffold. This study provides a computational approach to determine the mechanical stimuli at the cellular level when cells are cultured in a bioreactor and to relate mechanical stimuli with cell differentiation.

JTD Keywords: Tissue engineering, Scaffold, Rapid prototyping, Computational fluid dynamics, Finite element

Gimenez-Oya, V., Villacanas, O., Fernàndez-Busquets, X., Rubio-Martinez, J., Imperial, S., (2009). Mimicking direct protein-protein and solvent-mediated interactions in the CDP-methylerythritol kinase homodimer: a pharmacophore-directed virtual screening approach Journal of Molecular Modeling , 15, (8), 997-1007

The 2C-methylerythritol 4-phosphate (MEP) pathway for the biosynthesis of isopentenyl pyrophosphate and its isomer dimethylallyl pyrophosphate, which are the precursors of isoprenoids, is present in plants, in the malaria parasite Plasmodium falciparum and in most eubacteria, including pathogenic agents. However, the MEP pathway is absent from fungi and animals, which have exclusively the mevalonic acid pathway. Given the characteristics of the MEP pathway, its enzymes represent potential targets for the generation of selective antibacterial, antimalarial and herbicidal molecules. We have focussed on the enzyme 4-(cytidine 5'-diphospho)-2-C-methyl-D: -erythritol kinase (CMK), which catalyses the fourth reaction step of the MEP pathway. A molecular dynamics simulation was carried out on the CMK dimer complex, and protein-protein interactions analysed, considering also water-mediated interactions between monomers. In order to find small molecules that bind to CMK and disrupt dimer formation, interactions observed in the dynamics trajectory were used to model a pharmacophore used in database searches. Using an intensity-fading matrix-assisted laser desorption/ionisation time-of-flight mass spectrometry approach, one compound was found to interact with CMK. The data presented here indicate that a virtual screening approach can be used to identify candidate molecules that disrupt the CMK-CMK complex. This strategy can contribute to speeding up the discovery of new antimalarial, antibacterial, and herbicidal compounds.

JTD Keywords: Solvent-mediated interactions, Protein-protein interactions, Molecular dynamics, Drug design, Intensisty-fading MALDI-TOF mass spectrometry