Tissue Engineering and Regenerative Medicine Volume 7 - 2019 | https://doi.org/10.3389/fbioe.2019.00109 Foam sclerotherapy is clinically employed to treat varicose veins It involves intravenous injection of foamed surfactant agents causing endothelial wall damage and vessel shrinkage Foam production methods used clinically include manual techniques such as the Double Syringe System (DSS) and Tessari (TSS) methods Pre-clinical in-vitro studies are conducted to characterize the performance of sclerosing agents; however the experimental models used often do not replicate physiologically relevant physical and biological conditions physical vein models (PVMs) were developed and employed for the first time to characterize the flow behavior of sclerosing foams PVMs were fabricated in polydimethylsiloxane (PDMS) by replica molding and were designed to mimic qualitative geometrical characteristics of veins Foam behavior was investigated as a function of different physical variables namely (i) geometry of the vein model (i.e. The experimental set-up consisted of a PVM positioned on an inclined platform a syringe pump to control the flow rate of a blood substitute The static pressure of the blood surrogate at the PVM inlet was measured upon foam administration The recorded pressure-time curves were analyzed to quantify metrics of foam behavior with a particular focus on foam expansion and degradation dynamics Results showed that DSS and TSS foams had similar expansion rate in the physiological PVM whilst DSS foam had lower expansion rate in the varicose PVM compared to TSS foam The degradation rate of DSS foam was lower than TSS foam the background flow rate had a significant effect on foam behavior enhancing foam displacement rate in both types of PVM It was however designed for usage under static fluidic conditions and it did not replicate the varicose vein architecture In order to address these limitations of previous test methods the work in this study aims to develop physical models replicating qualitative architectural characteristics of varicose veins and to employ them as a screening platform for comparing the flow behavior of different foam formulation methods The developed biomimetic-inspired vein model (referred to as physical vein model or PVM) allows recapitulating features of physiological and varicose veins and physiologically relevant flow conditions PVMs were employed to compare the flow performance of polidocanol-based PCFs as a function of vessel geometry (straight vs it was demonstrated that models can be coated with endothelial cells enabling future investigations of both mechanical and biological performance of sclerosing agents Assembled PVMs representing simplified physiological (A) and varicose (B) vein models and demonstration of perfusion using a red dye These models were employed to test the flow behavior of foams The main channel was punctured with a needle (16G) in order to mimic the clinical process of injection more closely Inlet/outlet ports were connected with silicone tubing (4 mm outer diameter Schematic of PVMs production via replica molding (A) a CAD design of positive molds was generated and 3D printed (B) two PDMS layers were produced from the molds and permanently bonded via oxygen plasma treatment to obtain a fully circular channel cross-section (C) the main channel was punctured with a 16G needle (D) the foam was produced by passing the polidocanol solution (liquid phase) from one syringe 10 times into and out of the other syringe initially containing room air The pressure transducer was connected to a National Instruments I/O module (NI-DAQ The NI-DAQ system supports analog and digital inputs and communicates with the NI-DAQ software (National Instrument USA) script was employed to store pressure data in an automated fashion Schematic illustrating the set-up for evaluating the flow behavior of foams in the PVM model A steady flow is imposed using a syringe pump (A) and a pressure transducer (B) is positioned in line with the inlet tubing and prior to the PVM (C) The foam is injected through the 16G needle into the main channel (D) The pressure transducer is connected to a National Instruments I/O module A Matlab code is employed to store pressure data in an automated fashion (E) The corresponding inlet Reynolds number was calculated using Equation 1 where ρ is the density of the blood surrogate (kg/m3) v is the mean velocity of the blood surrogate (m/s) D is the hydraulic diameter of the vein model (m) and μ is the dynamic viscosity of the blood surrogate (Pa·s) The Reynolds number in these experiments ranged from 1.65 to 3.34, which is ~100–200 times lower than physiological values (Raju et al., 2004) in order to replicate quasi-static or impaired flow conditions occurring in diseased veins (ii) an almost linear increase in pressure due to the expansion of the foam within the PVM and (iii) an almost linear decrease in pressure caused by foam degradation and “washing out.” The CfAS calculated the slope of phases (ii) and (iii) which were referred to as expansion rate (ER) and degradation rate (DR) expansion time (ET) and degradation time (DT) were also quantified The plots were divided into three phases: (i) an increase of the pressure due to the expansion of the foam inside the channel (ii) a decrease of the pressure caused by the degradation of the foam and (iii) an initial peak due to the injection A representative pressure profile is illustrated in Figure 4 from which four different phases can be identified: This phase was associated with a rapid pressure spike likely due to the insertion of the needle within the PVM or other mechanical perturbations associated with the injection procedure Peak pressure values in this phase ranged from 10 to 50 mmHg The CfAS allowed quantifying the slope of the plug expansion phase which was herein referred to as expansion rate (ER) It is hypothesized that more cohesive foams would fractionate more slowly and dilute less rapidly with the blood substitute thus resulting in lower ER and higher peak pressures The expansion time (ET) was also calculated from the pressure-time curve; higher ET corresponds to a longer contact time between the foam and the inner surface of the PVM The CfAS allowed quantifying the slope of the plug degradation phase which was referred to as degradation rate (DR) It is hypothesized that more cohesive foams would destabilize more slowly The degradation time (DT) was also calculated from the pressure-time curve; higher DT corresponds to a longer contact time between the foam and the inner surface of the PVM The pressure level at the end of the degradation phase could be equal or greater than the initial pressure (i.e. bubbles accumulated at the top surface of the vein model they would remain in place for the duration of the pressure recording (up to 100 s) causing the residual pressure to be greater than the initial pressure Comparisons between foam production methods were performed using an unpaired Student's t-test for two groups analysis or one-way ANOVA in the case of more than two groups Statistical significance was assumed for p < 0.05 All statistical tests were performed with Prism software (GraphPad Software Inc. All data were reported as the mean ± SD of at least six independent repeats of the same experiment HUVECs (human umbilical venous endothelial cells) were seeded in the PVM models HUVECs were extracted from umbilical cords and this was carried out in accordance with the Human Tissue Act (2004) and the recommendations of Southampton & South West Hampshire Research Ethics Committee B with Governance provided by the University of Southampton Research Governance Office Umbilical cords were collected from the Princess Anne Hospital (Southampton UK) from non-complicated natural vaginal births following agreed ethical collection protocols [Local Research Ethical Committee (LREC); Ref: 07/H0502/83] the PDMS devices were washed four times with 70% ethanol for sterilization HBSS (Hanks' Balanced Salt Solution USA) was conveyed through the channels to remove ethanol traces the inner surfaces of the device were coated with 3 mL of different proteins: 50 μl/mL of rat type I collagen (100 μg/mL; Gibco™ UK) or (ii) 100 μl/mL of fibronectin (Sigma The device was subsequently placed at 37°C in a 5 % CO2 incubator for 2 h The HUVECs suspension was injected into the proteins-coated channels (at a concentration of 4–5 × 106 cells/mL) The device was finally incubated overnight to promote cell attachment Bright field images of HUVECs within the PVM models were acquired with an optical microscope (Olympus Images were taken of live samples every 24–48 h The PVM models developed in this study have been employed to characterize the flow behavior of sclerosing foams by measuring static pressure of a blood surrogate at the PVM inlet only the pressure spike corresponding to foam injection was present due to the absence of a background pressure-driven flow Pressure plots obtained from the CfAS while injecting DSS (A,C) and TSS (B,D) foams The upper panels show the recordings obtained using the straight channel geometry whereas the lower panels show the recording obtained using the serpentine-like channel geometry the pressure was measured without injecting foam at a background constant flow rate (72 mL/h) The pressure inside the vein models was also measured while injecting the foam in the absence of a blood-surrogate flow (F) Expansion rate values at the different flow rates investigated and different foam production methods (TSS blue bars Measurements were obtained from the straight (A) and curved (B) channel geometry Data represent the average of 6 measurements ± SD An asterisk (*) indicates that differences between mean values are statistically significant (p < 0.05) and results are reported as mean value ± standard deviation A one-way ANOVA was conducted to evaluate whether ER of both types of foam depended on the inlet flow rate a significant difference was observed between all flow rates investigated; indeed the average ER value increased with increasing the inlet flow rate increasing the background flow caused a reduction of ET likely due to a “washing out” effect of the blood surrogate DSS foam demonstrated greater ability to oppose this effect resulting in higher contact time with the vessel wall Expansion time values at the different flow rates investigated Two asterisks (**) indicate that differences between mean values are very statistically significant (p < 0.01) Degradation rate values at the different flow rates investigated Figure 8 shows DR values determined for both TSS and DSS foams DSS foam had a lower DR (0.26 ± 0.01 mmHg/s) compared to TSS foam (0.32 ± 0.03 mmHg/s) Significant difference in DR between TSS (0.45 ± 0.06 mmHg/s) and DSS (0.37 ± 0.09 mmHg/s) foams was found at 72.0 mL/h At the highest flow rate (125.0 mL/h) both types of foam had similar DR (0.48 ± 0.03 mmHg/s for TSS and 0.47 ± 0.07 mmHg/s for DSS) suggesting that foam degradation performance is dominated by the background flow at these higher flow rates A one-way ANOVA was performed to evaluate the effect of background flow rate on DR a significant difference was observed with increasing the inlet flow rate; indeed the average DR value increased with increasing the flow rate no significant difference was found by varying the flow rate Both foams had comparable degradation performance across the two PVM geometries suggesting that once a foam plug has been established into the vein its degradation dynamics is not significantly affected by the vessel architecture Figure 9 shows DT values determined for both TSS and DSS foams at all flow rates investigated DSS foam had a statistically higher DT (20.73 ± 2.60 s) compared to TSS foam (15.46 ± 2.30 s) the average DT value was still higher compared to TSS (16.38 ± 3.70 s and 10.90 ± 3.40 s At the highest flow rate (125.0 mL/h) both PCFs had comparable DT (9.97 ± 5 s for TSS and 7.8 ± 5.6 s for DSS) foams with lower degradation rate had a longer degradation time Similar observations were made using the varicose PVM model where at the lowest flow rate (62.5 mL/h) DSS foam had statistically higher DT (24.8 ± 6.20 s) compared to TSS Differences between foam types reduced with increasing the inlet flow rate both types of foam presented similar DT (9.1 ± 5.3 s for TSS and 11.07 ± 4.1 s for DSS) A one-way ANOVA was conducted to evaluate the effect of flow rate on DT for both types of foam results show that there is a significant difference between DTs measured at increasing flow rates DT was not significantly influenced by the PVM geometry no significant difference was found by varying the background flow rate Degradation time values at the different flow rates investigated An asterisk (*) indicates that differences between mean values are statistically significant (p < 0.05) and two asterisks (**) indicate that differences between mean values are very statistically significant (p < 0.01) Bright field microscope images of HUVECs cultured within the fully circular PVM channels Images on the left show the lower channel wall coated with collagen (A) and fibronectin (C) Images on the right show the upper channel wall coated with collagen (B) and fibronectin (D) Images (4x magnification) were taken after 48 h from cell seeding a novel experimental method to quantify and compare the flow behavior of sclerosing foams was developed The method provided a quantitative determination of fluid pressure upon foam administration within models of either physiological or varicose veins (referred to as physical vein model When a cohesive foam is injected into the PVM, it forms a plug that displaces the blood substitute. The foam plug however degrades over time, due to its intrinsic instability and the “washing out” action of the background flow. Using our model system, we were able to characterize these phenomenological behaviors for the first time, by measuring the static pressure of a blood surrogate at the PVM inlet (Figure 5) It is well-known that sclerosing foams produced using different techniques differ in their “static” physical properties we evaluated for the first time the dynamic flow behavior of sclerosing foams by analyzing their expansion and degradation within qualitative models of both physiological and varicose veins the behavior of different PCFs was compared at varying volumetric flow rates (in the range 62.5−125.0 mL/h) Results also demonstrated that the flow field within the target vein can significantly influence the expansion dynamics of sclerosing foams DSS foam was slightly less sensitive to changes in the background flow rate suggesting that more cohesive foams may offer higher resistance to the “washing out” effect of the blood flow during expansion Reducing blood flow rate during administration (i.e. via vein compression) may thus be preferable to enhance therapeutic efficacy With respect to the degradation dynamics of PCFs (Figures 8, 9), at the lowest flow rate investigated DSS foam had lower degradation rate compared to TSS foam. This was likely due to the slower coarsening and drainage rate of DSS foams, coherently with previous studies (Carugo et al., 2015) Increasing the inlet flow rate resulted in PCFs having comparable degradation rate suggesting that foam degradation performance is dominated by the background flow in these conditions there was no significant difference in the degradation dynamics between the two PVM geometries investigated; suggesting that once a foam plug has been established into the vein It is important to highlight that expansion and degradation dynamics taken at the lowest flow rate are likely to be more representative of the flow conditions in a diseased (i.e. DSS presented a slightly superior performance compared to TSS Finally, it was demonstrated that PVM models can be lined with endothelial cells in order to recreate the endothelial layer (Figure 10) The degree of endothelial damage upon treatment with foam can be employed as an indicator of therapeutic efficacy we will employ these cell-coated PVM models to investigate the biological effects of sclerosing foams and correlate them with foam mechanical behavior we described the development of physical vein models replicating the qualitative architecture of physiological and varicose veins and their utility as model platforms to screen the flow behavior of sclerosing foams upon different formulation and administration conditions A simple method to manufacture vein models was developed which aimed at generating channels with circular section and with a geometry that recapitulates some characteristics of the varicose vein An experimental protocol was also established to investigate the flow performance of foams at conditions relevant to their clinical administration the experimental set-up replicated some aspects of the clinical process of foam injection and the presence of a background blood flow Fluid pressure at the PVM inlet was measured during foam administration which revealed different phases of the foam expansion and degradation dynamics Particular emphasis was given to expansion and degradation of the foam plug As reported in previous studies (Carugo et al., 2015), the cohesiveness of foams is highly dependent on their rheological properties, which in turn are influenced by the bubble size distribution and foam drainage kinetics. Previous results showed that foam produced using the DSS method were more stable and presented longer dwell time compared to TSS foams (Carugo et al., 2015) Consistently with these previous observations in our dynamic study DSS foam had longer degradation time and slower degradation rate than TSS foam no significant difference between the two foam formulations was found in the physiological vein model whereas DSS had slower expansion in the varicose model Differences in foam behavior across different model geometries could be attributed to the broader bubble size distribution of TSS foam compared to DSS foam; although these aspects merit further investigation the physical vein models and experimental methods developed in this study provide a novel technology platform to measure the behavior of different formulations of sclerosing foams at physical conditions that resemble their clinical administration They could therefore be employed as an additional test method in the pre-clinical pipeline to innovate foam formulation and administration procedures we demonstrated that PVM models are suitable for coating with endothelial cells which enables future investigations to correlate flow performance of sclerosing agents with their biological effects It should be noted that the PVM models reported in this study do not replicate the presence of venous valves or the mechanical properties of the vein wall which may affect the flow behavior of foams Ongoing research is focusing on the incorporation of these additional architectural and functional characteristics a more faithful replication of the physiological boundary conditions (including changes due to clinical practices; i.e. vein compression) will be considered in the future EB conducted all experiments and data analyses DC co-designed all experimental procedures co-wrote the manuscript and supervised the project and VP contributed to the design and implementation of the research and to the analysis of the results All the other authors helped supervise the project This research was supported by Doctoral Training Partnership funded from EPSRC and Biocompatibles UK Ltd The funding was awarded by the Faculty of Engineering and Environment (University of Southampton) EB is in receipt of a Doctoral Training Partnership funded from EPSRC and Biocompatibles UK Ltd and VP are paid employees of Biocompatible UK Ltd The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest Studies on foam decay trend and influence of temperature jump on foam stability in sclerotherapy Hydrophilization and hydrophobic recovery of PDMS by oxygen plasma and chemical treatment—An SEM investigation Sclerosant foam structure and stability is strongly influenced by liquid air fraction The role of clinically-relevant parameters on the cohesiveness of sclerosing foams in a biomimetic vein model Benefits of polidocanol endovenous microfoam (Varithena®) compared with physician-compounded foams Foam sclerotherapy techniques: different gases and methods of preparation Fabrication of a circular PDMS microchannel for constructing a three-dimensional endothelial cell layer Cell migration into scaffolds under co-culture conditions in a microfluidic platform Detergent sclerosants are deactivated and consumed by circulating blood cells Size of sclerosing foams prepared by ultrasound and the handmade tessari method for treatment of varicose veins: sclerosing foams for treatment of varicose veins Polidocanol for endovenous microfoam sclerosant therapy The care of patients with varicose veins and associated chronic venous diseases: clinical practice guidelines of the society for vascular surgery and the American venous forum Sclerotherapy: Treatment of Varicose and Telangiectatic Leg Veins Google Scholar Evaluation of the efficacy of polidocanol in the form of foam compared with liquid form in sclerotherapy of the greater saphenous vein: initial results Comparison of 1% and 3% polidocanol foam in ultrasound guided sclerotherapy of the great saphenous vein: a randomised Influence of surfactant on gas bubble stability Systematic review of foam sclerotherapy for varicose veins perfusable 3D microvascular networks on a chip CrossRef Full Text | Google Scholar Rheology of colloidal gas aphrons (microfoams) Comparative stability of sodium tetradecyl sulphate (std) and polidocanol foam: impact on vein damage in an in-vitro model Properties of polidocanol foam in view of its use in sclerotherapy Foam stability in the presence and absence of hydrocarbons: From bubble- to bulk-scale Efficacy of polidocanol foam versus liquid in sclerotherapy of the great saphenous vein: a multicentre randomised controlled trial with a 2-year follow-up The lytic effects of detergent sclerosants on erythrocytes endothelial cells and microparticles are attenuated by albumin and other plasma components in vitro An investigation into the influence of various gases and concentrations of sclerosants on foam stability An investigation of side-effects and efficacy of foam-based sclerotherapy with carbon dioxide or room air in the treatment of reticular leg veins: a pilot study Blood viscosity in tube flow: dependence on diameter and hematocrit Venous flow restriction: the role of vein wall motion in venous admixture Stability of foam in sclerotherapy: differences between sodium tetradecyl sulfate and polidocanol and the type of connector used in the double-syringe system technique Diameter-reflux relationship in perforating veins of patients with varicose veins Foam and liquid sclerotherapy for varicose veins CrossRef Full Text | Google Scholar Chronic venous insufficiency: a frequently underdiagnosed and undertreated pathology Novel developments in foam sclerotherapy: Focus on Varithena® (polidocanol endovenous microfoam) in the management of varicose veins Preliminary experience with a new sclerosing foam in the treatment of varicose veins Microfluidic approaches for engineering vasculature CrossRef Full Text | Google Scholar Polidocanol concentration and time affect the properties of foam used for sclerotherapy Physiochemical properties and reproducibility of air-based sodium tetradecyl sulphate foam using the Tessari method Basic physiochemical and rheological properties of detergent sclerosants Millar TM and Carugo D (2019) Physical Vein Models to Quantify the Flow Performance of Sclerosing Foams Received: 31 October 2018; Accepted: 01 May 2019; Published: 21 May 2019 Copyright © 2019 Bottaro, Paterson, Zhang, Hill, Patel, Jones, Lewis, Millar and Carugo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY) distribution or reproduction in other forums is permitted provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited in accordance with accepted academic practice distribution or reproduction is permitted which does not comply with these terms *Correspondence: Dario Carugo, ZC5jYXJ1Z29Ac290b24uYWMudWs= Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher 94% of researchers rate our articles as excellent or goodLearn more about the work of our research integrity team to safeguard the quality of each article we publish The dates displayed for an article provide information on when various publication milestones were reached at the journal that has published the article activities on preceding journals at which the article was previously under consideration are not shown (for instance submission All content on this site: Copyright © 2025 Elsevier B.V. Metrics details This paper provides an analysis of microfluidic techniques for the production of nanoscale lipid-based vesicular systems In particular we focus on the key issues associated with the microfluidic production of liposomes production method (including microchannel architecture) and drug loading in determining liposome characteristics we propose microfluidic architectures for the mass production of liposomes with a view to potential industrial translation of this technology They are composed of one or more closed shells consisting of a phospholipid bilayer and enclosing a small volume of aqueous liquid Their dimensions can vary from tens of nanometres up to hundreds of micrometres depending on the preparation protocol and the final use liposomes used in medical applications are unilamellar with an average size of ~100 nm although they are less frequently employed in association with commercial formulations Over the past few decades different techniques have been proposed for the preparation of liposomes and liposome size, dispersity (Đ)9 lamellarity and entrapped liquid volume may all vary depending on the method used The majority of the existing methods can be categorised into one of two groups of bulk methods where the term “bulk” is used throughout this manuscript to describe macroscale or batch techniques More recently, improvements have been made with the development of microfluidic production methods10 in which liposome formation occurs within a confined microenvironment Microfluidic methods have demonstrated potential for achieving higher control over the physical properties of the end product There are however concerns about the viability of this technology on an industrial scale These concerns arise mostly from the lack of experimental evidence demonstrating: (i) a systematic superiority compared to macroscale bulk methods (ii) suitability for producing clinically relevant liposome formulations (iii) efficiency in loading bioactive compounds it should be noted that efforts have been made in the industrial environment to develop controlled precipitation processes for producing liposomes of clinically relevant standards we attempt to address these issues by critically analysing state-of-the art microfluidic methods reported in the literature and by performing additional experiments Many of the aforementioned techniques however result in the production of relatively large vesicular systems (i.e. Giant or Large Unilamellar Vesicles) or microtubules limiting their potential for the production of nanoscale drug delivery vehicles being simply a miniaturised version of their analogous bulk counterparts many of the aforementioned methods cannot be regarded as strictly microfluidic the microfluidic hydrodynamic focusing (MHF) approach presents the typical physical characteristics of microfluidic systems (i.e. including low Reynolds number and diffusion dominated mass transfer) and represents the only viable microfluidic method for producing lipid-based nanoscale vesicular systems with potential for clinical application Only limited attention has been devoted to this technique in previous publications5 however and a critical evaluation of the method has not been carried out which are routinely employed in the industry may not represent a scaled up version of microfluidic flow focusing regimes which may be more suitable for medical or biotechnological applications (e.g as delivery systems for anticancer drugs or as transfection reagents) and which are routinely used in the industry Although this effect could be reduced by significantly diluting the liposome suspension prior to dimensional analysis this may not be applicable to many MHF techniques which suffer from a relatively low liposome concentration in the end-product there has been little or no systematic characterisation of drug encapsulation efficiency or the effect of drug encapsulation on liposome size different device architectures have been considered for liposome production Devices include microscale chips (cross-sectional dimensional range of 100–320 μm) with mixing channel displaying distinct architectural features (i.e. and containing micropillar structures) and scaled-up versions of microscale flow focusing architectures with cross-sectional channel dimension in the millimetre range Considering the spectrum of different architectures described in the literature the following general principles have been considered in this study for designing and constructing chips employed for liposome production The total width and length of the chip should be defined to comply with the dimensional constraints imposed by the typical size of microscope stages that are routinely used to assess the correct functioning of the chips the microchannel architecture should fit within the maximum width and length of conventional glass slides for microscopy (i.e. the thickness of the device should comply with the working distance of microscope objectives It should be noted that whilst in situ microscope observation may not be required in an industrial setting we consider it as a desirable requirement at a research and development stage The inlet channels should be long enough to allow the fluid flows to fully develop before they intersect with each other allowing for a stable and predictable flow focusing regime Channels should also be designed so as to guarantee sufficient spacing between the inlet ports allowing for robust and practical connection with tubing and pumping units A typical inlet channel length for MHF chips is in the range of 5–10 mm; note that longer channels result in higher backpressure which could potentially compromise device usability the mixing channel should be designed to allow for complete mixing of the solvents (ethanol or isopropanol and water) at the selected operating flow rates and to comply with the dimensional requirements imposed by microscope interfacing (iii) Relative orientation of the inlet channels the angle between the side and central inlet channels should be defined so as to minimise fluid dynamic perturbations at the intersection between flows particularly if devices are designed to operate at high throughput regimes An angle in the range 30°–60° was deemed suitable for the applications described in the present study The constitutive materials should satisfy different requirements compatibility with microfabrication and bonding techniques PDMS has been identified as the optimal material choice given its compatibility with pharmaceutical-grade solvents (i.e. ease of surface treatment by exposure to oxygen plasma and potential for bonding with commercially available glass slides Traditional soft lithographic methods based on SU-8 moulding can be employed for producing PDMS microchannel architectures which can be irreversibly coupled to glass substrates using plasma bonding techniques glass is chemically inert and its surfaces smoothness is suitable for optical microscopy With an aim to develop cost-effective and facile microfabrication methods with potential for large-scale chip fabrication we also employed an in-house developed technique which combines micromilling with epoxy-based replica moulding Both these methods allow for obtaining well-defined cross-sectional channel features which may be complex to achieve using wet etching techniques Microfluidic flow focusing was selected as a microfluidic-based regime for nanoscale liposome production and its performance compared with other micromixing geometries MHF operated at laminar flow conditions has been previously demonstrated to allow for precise control over the interfacial boundaries between solvent and co-solvent streams resulting in diffusion dominated mass transfer and leading to liposomes of relatively uniform physical properties liposome size and dispersity in MHF chips has the advantage of being potentially controlled on-demand by finely adjusting the hydraulic boundary conditions MHF devices were compared with other micromixing geometries including those containing serpentine-like microchannels and micropillar arrays Panel A shows a schematic representation of a microfluidic device and the process of liposome (SUV) self-assembly Panel B shows a schematic representation of the ethanol injection procedure Panel C illustrates the geometrical characteristics of the chips employed for the microfluidic experiments #chip1-MHF comprises three inlet microchannels with 30° intersection angle between each other #chip2-YJ comprises two inlet microchannels with 120° intersection angle between each other #chip3-MP comprises three inlet microchannels with 90° intersection angle between each other and a mixing microchannel with three pillar mixing elements TeflonTM tubes were used to connect the microfluidic platform with syringes The volumetric flow rate was controlled using syringe pumps a detailed analysis is provided of the effects of the operating parameters (especially FRR) on liposome dimensions taking into consideration the effect of residual alcohol on the viscosity of the liposomal samples and thus on the determination of liposome size by light scattering measurements Liposomes produced using MHF were compared with those obtained by the bulk ethanol injection method Ethanol injection was selected as a bulk technique since in our opinion it strongly resembles MHF liposome formation (in terms of simplicity and physical conditions) relative to other bulk techniques this method results in the formation of small unilamellar vesicles (SUV); therefore it does not require post-processing homogenisation steps phospholipids are firstly dissolved in ethanol; a small amount of the lipid solution is then injected into water Effect of the variation of lipid concentration (A) ethanol content (B) and simultaneous variation of both parameters (C) on the dimension of liposomes produced by ethanol injection data relative to liposomes prepared by MHF microfluidics are also reported in panel C (dashed line) Liposomes were constituted of PC/cholesterol 4.0−0.4 mM and data represent the mean of three independent samples Effect of the variation of lipid composition on the size (upper part panel B) of liposomes produced by MHF microfluidics Liposomes were constituted of PC/DDAB 9.0−1.0 mM (filled circles dashed line) or PC/cholesterol 9.0−1.0 mM (open circles and are reported as the mean of three independent samples cryo-TEM and macroscopic aspect (insets) of empty PC/cholesterol (A) and PC/DDAB (C) are reported images of the corresponding ivermectin loaded liposomes are also reported (B,D) It can therefore be concluded that microfluidic methods are potentially suitable for the preparation of liposomal suspensions to be included in commercial in order to minimise the content of toxic constituents we demonstrate that ethanol is a suitable substitute for the more toxic isopropanol in microfluidic experiments particularly in view of a potential translation of this technology into the industrial environment Data refer to Z-average of liposomes produced by #chip1-MHF (circles) #chip2-YJ (hexagons) or #chip3-PM (triangles) Liposomes were constituted of PC/DDAB 9.0−1.0 mM produced at the indicated FRR and TFR = 37.50 μl/min Data are reported as the mean of three independent samples the size of empty (water-filled) liposomes is also reported (open circles) Liposomes were produced by #chip1-MHF at TFR = 37.50 μl/min and FRR = 30 (A) Geometrical characteristics of easy-to-build chips for liposome production and potential parallelized network for mass production of liposomes (C) Schematic of #chip5-MHF-LC and related fabrication method developed in house (μMi-REM) (A,B) Dependence of liposome mean size (A) and dispersity index (B) on the flow rate ratio (FRR ranging from 5 to 100) at a fixed TFR of 6 ml/min (C,D) Dependence of liposome mean size (A) and dispersity index (B) on the total flow rate (TFR ranging from 3 to 18 ml/min) at a fixed FRR of 100 Liposome size was observed to reduce with increasing TFR (at a fixed FRR of 100) while increasing FRR resulted in increased liposome size This finding appears to be contrary to previous results obtained using different microscale geometries and therefore merits further investigations (i.e. particularly on the fluid dynamic and mass transport phenomena occurring within these scale-up architectures) these components also present low inter-sample variability which is of particular importance for parallelisation purposes Despite the potential advantages offered by the proposed microfluidic systems a more pervasive understanding of the effect of the microfluidic architecture on liposome characteristics may be beneficial to optimise device performance and their scaling-up for mass production of liposomes This will represent the subject of future investigations The potential for microfluidic-based technology to produce nanoscale lipid vesicles for medical applications has already been demonstrated the translation of this method to the industrial environment has been hindered by several limiting factors we critically review the state-of-the-art microfluidic methods and perform additional experiments in an attempt to address the associated research questions The following concluding remarks can be drawn from our studies: (i) A suitable bulk counterpart to MHF microfluidic methods has been identified (referred to as “controlled ethanol injection”) which makes it possible to carry out a systematic comparison between the two methods we have demonstrated that liposomes produced using microfluidics were smaller and more uniform in size than the ones produced by controlled injection whilst no appreciable differences in liposome dimensional stability were found over time (data not shown) MHF thus has the potential advantage of being able to generate a wider range of mean liposome sizes whilst maintaining lower size dispersity compared to bulk counterparts application-specific tuning of liposome size through changes of the hydrodynamic boundary conditions (ii) Both lipid and ethanol concentration have a significant effect on liposome properties (in both bulk and microfluidic methods) This is of particular importance when results taken at different FRRs are compared and interpreted changing FRR not only causes changes in the fluid dynamic field (i.e. width of the focused stream) but also in the physico-chemical properties of the fluidic environment (iii) The lipid formulation plays an important role in determining liposome properties charged lipids generated smaller vesicles compared to uncharged lipids MHF microfluidic methods were also found to be suitable for producing different liposome formulations (iv) Microfluidic methods are suitable for producing small and uniform liposomes at relatively high concentrations of lipids and using solvents with relatively low toxicity (i.e. This has important implications for the potential utility of this technology in the pharmaceutical industry freezing and thawing) is required with MHF methods to obtain liposomes of desirable characteristics (v) It is possible to efficiently encapsulate biologically active compounds in liposomes produced using microfluidics on-chip liposome loading with bioactive compounds (both hydrophilic and lipophilic) has seen a limited number of advancements in recent years and may represent an exciting avenue of research in the near future (vi) High-throughput and easy-to-build microfluidic architectures can be constructed for mass production of liposomes without significantly influencing the quality of the end-product compared to conventional MHF devices we have demonstrated liposome production at volumetric flow rates of up to 18 ml/min Highly pure phosphatidylcholine (PC) 90% from soybean was purchased from Phospolipon 90G Lipoid Germany; and dimethyldioactdecyl-ammoniumbromide (DDAB) from Sigma-Aldrich 1,2-dimyristoyl-sn-slycero-3-phosphocholine (DMPC) and 1,2-dimyristoyl-sn-glycero-3-phosphoethanolamine-N-[methoxy(polyethylene glycol)-2000] (ammonium salt) (14:0 PEG2000 PE) were purchased from Avanti Polar Lipids; polydimethylsiloxane (PDMS) monomer Sylgard® 184 and curing agent were purchased from Dow Corning (USA); SU-8 photoresist was obtained from Chestech (UK); trichloro (1 H,1 H,2 H,2 H)-perfluorooctylsilane and suramin sodium salt were purchased from Sigma-Aldrich (UK) All the other regents and solvents not included in this section were also purchased from Sigma-Aldrich (UK) The “Cross-flow” #chip1-MHF consists of polydimethylsiloxane (PDMS) and glass layers and was produced with conventional soft lithography techniques an SU-8 mould with the designed microchannel pattern consisting of three inlets and one reaction channel was prepared following the standard lithography protocols the mould was covered with a layer of a 10: 1 (w/w) PDMS precursor and curing agent mixture and heated for 1 h at 80 °C for the polymer to cure The PDMS sheet with the microchannel architecture fabricated on the surface was then removed from the mould and permanently bonded to a glass slide after oxidizing its surface via plasma treatment The “Snake mixer slide” #chip2-YJ was obtained from Thinxxs The mixer was integrated on a microscope slide made of cyclic olefin copolymer (COC) For the fabrication of the “three-inlets pillars” #chip3-PM a photolithography/wet etching procedure was used the channel network was designed using AutoCAD drawing software first A film negative of the desired final size was then prepared by a commercial photo mask producer to form the optical mask Crown white glass plates (thickness of 1.5 mm) coated with a thin layer of chromium metal mask plus an upper layer of positive photoresist were used for channel network fabrication the pattern of interconnecting channels was transferred from the negative film to the photoresist layer on the glass which was then developed and removed to reveal the channel areas of the glass to be etched The glass plate was baked in an oven at 80 °C overnight to dry and harden the mask on glass The channels were then etched using 1% hydrofluoric acid buffered with 5% ammonium fluoride solution at 65 °C under ultrasonic agitation (Ultrasonic Cleaner the etched glass was thermally bonded (595 °C for 3 h) with a top plate of the same material into which outlet and inlet holes had been previously drilled the cross-sectional profile of the etched microchannel was measured using a surface profilometer The ethanol injection method was applied to prepare liposomes The required quantities of lipids were dissolved in ethanol: PC 90G (40–90 mM) The resulting organic solution (0.5 ml) was then injected by a syringe pump (KD scientific New Era pump System) at flow rate of 500 μl/min under magnetic stirring (300 rpm) in an appropriate volume of water (4.5 ml) Spontaneous liposome formation occurred as soon as the ethanolic solution came into contact with the aqueous phase The liposome suspension was then kept under stirring by vortexing for 5 minutes at room temperature The obtained liposomal suspension was stored at 20 °C Liposomes were prepared by injecting a lipid mixture (PC 90G 90 mM DDAB 10 mM) dissolved in ethanol into the central channel of the microfluidic network of #chip1-MHF; water was injected into two oblique side channels intersecting with the central channel The flow rate ratio (FRR) is defined here as the ratio between the water volumetric flow rate and the ethanol volumetric flow rate Liposome formation at different shear forces was investigated by changing the total flow rate (TFR) from 18.75 to 75.00 μl/min The flow focusing was observed and monitored with a Dino-Eye Eyepiece camera (Dino-lite Digital Microscope) Liposomes were prepared via a similar procedure to that applied for the #chip1-MHF apart from the fact that the lipid mixture (PC 90 mM and DDAB 10 mM) dissolved in ethanol was injected into one channel of the Y geometry while water was injected into the other channel TFR was set between 18.75 and 75.00 μl/min The preparation of liposomes by the “Three-inlets pillars” chip (#chip3-PM) was analogous to the procedure described for “cross-flow” (#chip1-MHF) The TFR was varied from 18.75 to 75.00 μl/min and the FRR was varied from 10 to 50 Microchip constituents were obtained from IDEX Health & Science #chip4-OFF3 components were made from PEEK and specifically designed with a 0.006” thru-hole that delivers a low swept volume and includes F-112 and P-416 fittings To connect the “off-the-shelf” devices to the syringes TeflonTM tubes with an inner diameter of 750 μm were employed USA) were used to control the flow rate of liquids pumped through the devices Liposomes were again prepared using a similar procedure to that above for #chip1-MHF apart from the fact that the central ethanol stream contained DMPC Cholesterol and 14:0 PEG2000 PE at a concentration of 10 The TFR was varied from 3 to 18 ml/min and the FRR was varied from 5 to 100 All dimensional analysis of the liposome dispersions were performed using a DLS Zetasizer Nano-ZS (Malvern Instruments UK) with a backscattering detection angle of 173° and a 4.0 mW power source was used to report the intensity mean diameter (Z-average) and the dispersity of the liposome formulations The mean particle size was obtained from the results of three experiments The size distribution was evaluated in terms of the dispersity index (Đ) The measurements of vesicle size and dispersity were carried out at 21 °C in water A 3 mL aliquot of sample solution was pipetted on to plasma-treated (Gatan Solarus Model 950 Advanced Plasma System time 30 s) carbon copper grids (Quantifoil R 3.5/1) in the environmental chamber of a fully automated vitrification device for plunge freezing (FEI Vibrot) having relative air humidity of 100% and a temperature of 22 °C The excess solution was removed by blotting with filter paper for 2 s followed by 1 s draining and plunging of the samples into 1:1 mixture of liquid ethane and liquid propane Vitrified samples were cryo-transferred into a Jeol JEM3200FSC cryo-TEM operating at −194 °C The temperature of the samples was −187 °C during image acquisition The microscope was operated in the bright field mode using a 300 kV acceleration voltage; the in-column energy filter was set to 0–20 eV energy-loss range (zero-loss imaging) Micrographs were recorded with a Gatan Ultrascan 4000 CCD camera The determination of drug entrapment efficiency in liposomes was performed by loading 250 μl of a liposome dispersion into a Sheparose® 4B column (1.0 cm in diameter and 10 cm long; GE Healthcare The void volume peak fractions containing drug loaded liposomes were collected and quantitated for drug content UV-VIS measurement (at 253 nm) was conducted by diluting 300 μl of each fraction with 700 μl of ethanol UV spectra were recorded with a Hewlett-Packard 8452 diode array spectrophotometer Drug encapsulation efficiency (EE) was calculated as follows: where C0 and C1 correspond to the total and liposome-associated amount of drug Liposome production by microfluidics: potential and limiting factors Preparation of vesicles (liposomes) In Encyclopedia of Nanoscience and Nanotechnology (ed Development of liposomal formulations: From concept to clinical investigations Liposomal drug delivery systems: From concept to clinical applications Effectiveness and Safety of Short Course Liposomal Amphotericin B (AmBisome) as First Line Treatment for Visceral Leishmaniasis in Bangladesh Eleven years of Inflexal V-a virosomal adjuvanted influenza vaccine Porcine surfactant (Curosurf) for acute respiratory failure after near-drowning in 12 year old Clinical development of liposome based drugs: formulation Liposome Drug Products-Guidance for Industry Dispersity in polymer science (IUPAC Recommendations 2009) Functional reconstitution of a voltage-gated potassium channel in giant unilamellar vesicles e25529; doi: 10.1371/journal.pone.0025529 (2011) Vesicles of variable sizes produced by a rapid extrusion procedure Unilamellar vesicle formation and encapsulation by microfluidic jetting Controllable monodisperse multiple emulsions Novel method for obtaining homogeneous giant vesicles from a monodisperse water-in-oil emulsion prepared with a microfluidic device Preparation of cell-encapsulation devices in confined microenvironment Octanol-assisted liposome assembly on chip A facile route to the synthesis of monodisperse nanoscale liposomes using 3D microfluidic hydrodynamic focusing in a concentric capillary array Microfluidic directed formation of liposomes of controlled size Microfluidic mixing and the formation of nanoscale lipid vesicles Microfluidic directed self-assembly of liposome-hydrogel hybrid nanoparticles Formation of liposome by microfluidic flow focusing and its application in gene delivery Microfluidic hydrodynamic focusing based synthesis of POPC liposomes for model biological systems Continuous flow production of cationic liposomes at high lipid concentration in microfluidic devices for gene delivery applications Microfluidic synthesis of PEG- and folate-conjugated liposomes for one-step formation of targeted stealth nanocarriers Microfluidic preparation of liposomes to determine particle size influence on cellular uptake mechanisms High-throughput manufacturing of size-tuned liposomes by a new microfluidics method using enhanced statistical tools for characterization Microfluidic-controlled manufacture of liposomes for the solubilisation of a poorly water soluble drug Microfluidic and lab-on-a-chip preparation routes for organic nanoparticles and vesicular systems for nanomedicine applications High-throughput continuous flow production of nanoscale liposomes by microfluidic vertical flow focusing Controlled self-assembly of monodisperse niosomes by microfluidic hydrodynamic focusing characterization and applications of liposomes: State of the Art Download references Authors would like to thank Dr Xunli Zhang (University of Southampton UK) for providing microfluidic devices fabricated using photolithography and wet etching and James Fisk and David Salisbury (Institute of Biomedical Engineering UK) for fabricating PMMA moulds used in the μMI-REM microfabrication method Mechatronics and Bioengineering Science research groups Faculty of Engineering and the Environment Department of Life Science and Biotechnology Elisabetta Bottaro & Claudio Nastruzzi All authors contributed to the writing process of the paper and all reviewed the final manuscript performed the experiments and prepared some of the figures designed the experiments and made the remaining figures The authors declare no competing financial interests Reprints and permissions Download citation Anyone you share the following link with will be able to read this content: a shareable link is not currently available for this article Sign up for the Nature Briefing newsletter — what matters in science 1: Structure of the guanidinase enzyme of the comammox species Nitrospira inopinata The presumed entrance to a tunnel leading to the active site is highlighted in the left image the tunnel is shown as a green line and guanidine as a stick model The structure was elucidated by the team led by Kristina Djinović-Carugo C: Kristina Djinović-Carugo/University of Vienna Research team identifies unconventional energy source for recently discovered "green" nitrifying bacteria An international research team led by the Centre for Microbiology and Environmental Systems Science (CeMESS) at the University of Vienna has discovered that comammox bacteria This unique ability opens new avenues for targeted cultivation of these enigmatic microbes and could also provide a key to reducing agricultural nitrous oxide emissions The research findings were recently published as an article in the prestigious journal Nature the conversion of ammonia via nitrite to nitrate is carried out by specialized microorganisms called nitrifiers This process is extremely important for the global biogeochemical nitrogen cycle in virtually all ecosystems but it plays an ambivalent role in global change nitrification contributes to the emission of the potent greenhouse gas and ozone-depleting substance nitrous oxide and leads to massive fertilizer losses in agriculture resulting in the eutrophication of water bodies nitrification is indispensable as a biological purification step for nutrient removal in wastewater treatment plants thus protecting water bodies from excessive nitrogen input from wastewater The study authors have now found a way that may promote nitrifiers in the environment that emit less nitrous oxide Comammox bacteria are considered "green" nitrifiers because they produce only small amounts of nitrous oxide as a byproduct of their metabolism and efficiently remove nitrogen compounds from wastewater in treatment plants Since the discovery of nitrifiers in the 19th century it was assumed that these microorganisms could only respire ammonia and urea the research groups led by Michael Wagner and Holger Daims demonstrated that some nitrifiers could also use the chemically unstable cyanate for their energy metabolism our team has now shown that comammox bacteria can also grow with the unconventional substrate guanidine," explains Marton Palatinszky "The comammox bacteria use a transporter and an enzyme structurally and functionally characterized in detail by us which allows them to produce ammonium from guanidine in a highly energy-efficient manner within the cell." Guanidine is a metabolic product of microorganisms and plants Little is known about its role in human and animal metabolism It is formed in soils during the degradation of synthetic fertilizer additives and in wastewater during the breakdown of the commonly used drug metformin little is known about the distribution and further processing of guanidine in the environment including microbiologists from the Helmholtz Centre for Environmental Research in Leipzig; Germany and Aalborg University in Denmark demonstrated that guanidine is present not only in human urine but also in livestock excreta and that comammox bacteria utilize guanidine in wastewater treatment plants They also showed that guanidine is metabolized by nitrifiers in agricultural soils New Opportunities for Cultivation and Nitrous Oxide Reduction The Vienna microbiologists are now attempting to enrich and isolate the widespread comammox bacteria from environmental samples using guanidine as only one strain is currently available in pure culture worldwide "This seems particularly promising as none of the other nitrifier strains we tested could grow with guanidine as the sole energy and nitrogen source," explains Katharina Kitzinger The team also wants to investigate whether adding guanidine to agricultural fertilizers could increase the abundance of comammox bacteria in arable soils thereby reducing agricultural nitrous oxide emissions "This work would not have been possible without the close collaboration of many researchers involved in the 'Microbiomes Drive Planetary Health' Cluster of Excellence We extend our sincere thanks to the Austrian Science Fund (FWF) for this special support," says study leader Michael Wagner Palatinszky M , Herbold CW, Sedlacek CJ, Pühringer D, Kitzinger K, Giguere AT, Wasmund K, Nielsen PH, Dueholm MKD, Jehmlich N, Gruseck R, Legin A, Kostan J, Krasnici N, Schreiner C, Palmetzhofer J, Hofmann T, Zumstein M, Djinovic-Carugo K, Daims H, Wagner M. Growth of complete ammonia oxidizers on guanidine. Nature.DOI: 10.1038/s41586-024-07832-z Volume 10 - 2023 | https://doi.org/10.3389/fmolb.2023.1155629 This article is part of the Research TopicWhen Predictions Meet Experiments: The Future of Structure DeterminationView all 6 articles Protein structure prediction and structural biology have entered a new era with an artificial intelligence-based approach encoded in the AlphaFold2 and the analogous RoseTTAfold methods More than 200 million structures have been predicted by AlphaFold2 from their primary sequences and the models as well as the approach itself have naturally been examined from different points of view by experimentalists and bioinformaticians we assessed the degree to which these computational models can provide information on subtle structural details with potential implications for diverse applications in protein engineering and chemical biology and focused the attention on chalcogen bonds formed by disulphide bridges We found that only 43% of the chalcogen bonds observed in the experimental structures are present in the computational models suggesting that the accuracy of the computational models is insufficient to allow the detection of chalcogen bonds according to the usual stereochemical criteria High-resolution experimentally derived structures are therefore still necessary when the structure must be investigated in depth based on fine structural aspects In 2021, sensational progress was made in protein structure prediction with AlphaFold2, the artificial intelligence system developed by DeepMind (Jumper et al., 2021). These predictions became freely available in the AlphaFold Protein Structure Database (AlphaFold DB), created by EMBL-EBI, which presently includes more than 200 million predictions (Tunyasuvunakool et al., 2021) (Varadi et al., 2022) Here, we seek to assess the degree to which these computational models can provide information on subtle details that may be important in various applications in protein engineering, chemical biology, and biotechnology. As an example, we focus on chalcogen bonds (referred to as ChB according to (Aekeroy et al., 2019)) formed by disulphide bridges (Aekeroy et al., 2019) (Vogel et al., 2019) This is an interesting test case because chalcogen bonds are not yet parameterized in any molecular mechanics/dynamics force field their presence cannot be affected by energy minimization protocols these moderate clashes can be tolerated if compensated by a good and native-like packing around them involving attractive interactions like hydrogen bonds (A) Schematic representation of a chalcogen bond; the position of the nucleophilic atom relative to the sulfur atom is monitored with two variables and α = 180 – θ where θ is the angle defined by the nucleophile and the atom covalently bound to the sulfur the differences between the distances d observed in the predicted structures (Alpha Fold DB) and in the experimental structures (PDB) the differences between the angles α observed in the predicted structures (Alpha Fold DB) and in the experimental structures (PDB) (D) Scatter plot of the Delta-α versus the Delta-d values Delta-d and Delta-α values are given in Å and degrees ChBs are an interesting test case because they are not yet parameterized in any molecular mechanics/dynamics force field their presence cannot be favoured by energy minimization protocols they would be considered inter-atomic clashes and consequently disfavoured A ChB can be formed by a nucleophilic atom and a chalcogen atom that is covalently bound to another atom (Figure 1A). The nucleophilic atom must be positioned along the prolongation of the covalent bond, or along the prolongation of one of the two covalent bonds if the chalcogen is divalent. As a consequence, there are two parameters that must be monitored, the distance d between the nucleophilic and the chalcogen atom and the angle α (Figure 1A) According to the IUPAC recommendations, the distance d must be shorter than the sum S of the van der Waals radii of the nucleophilic and chalcogen atoms (Aekeroy et al., 2019), despite the fact that the use of van der Waals radii in determining non-bonding interactions may need to be revised (Politzer et al., 2007) to account for the lower accuracy of macromolecular structures relative to small molecule structures a in a ChB d must not be larger than S + 0.1 Å which is supplementary to the angle θ In chemical crystallography and material science a ChB is usually characterized by a value smaller than 20° to account for the lower accuracy of macromolecular structures we increased this threshold value to 25° Similar settings were previously used (Carugo, 2023) (Carugo, 2023) and can be compared with estimated average positional standard errors of 0.046 Å which imply estimate errors of about 0.06 Å on bond distances and of about 2° on bond angles About one-half (43%) of the ChBs observed in the experimental structures are present in the computational models if the same stereochemical criteria are used This means that about one-half of the ChBs observed experimentally in high resolution crystal structures are not observed in the models deposited in AlphaFold DB Does this indicate that these models are wrong models available in AlphaFold DB are extremely similar to the experimental crystal structures which are very small (0.046 ± 0.002 Å) as expected for high-resolution structures the average absolute value of Delta-d is about three times greater than that calculated on the contacts (shorter than 3.5 Å) between main-chain oxygen and nitrogen atoms involving the residues that form ChBs (0.106 ± 0.008 Å) ChB predictions are much less accurate than main-chain hydrogen bond predictions Figure 2A shows an example in which the ChB is detected in both the experimental structures and in the computational model. The sulfur-oxygen distance d is even shorter in the AlphaFold DB model and the angle α is even closer to 0° in the AlphaFold model. Figure 2B shown an example in which the ChB is detected in the experimental structure and not in the AlphaFold DB model The local stereochemistry is quite well predicted but the sulfur-oxygen distance is slightly too long (the threshold is sum of the van der Waals radii 1.52 Å for oxygen and 1.80 Å for sulfur with and small positive tolerance of 0.1 Å but in terms of chemical interactions is crucial Comparison between the experimental and computation ChBs in three selected examples The experimental structure is represented with ball and sticks and the computational model only with sticks In (A) the ChB is detectable in both structures; in (B) the ChB is detectable only in the experimental structures since the distance d is slightly too long in the computational model; in (C) the ChB is detected in the experimental structures and it cannot be detected in the computational model where the disulphide bond is absent–this is a region partially unstructured of the protein carbon light grey in the experimental structure and azure in the AlphaFold DB model it would be possible to increase the threshold values of d and α that allow to automatically detect ChBs in such a way to increase the number of ChBs in the models of AlphaFold DB the values of these thresholds strictly depend on the laws of chemistry and physics and nothing indicates here that this is justified AlphaFold2 is a powerful tool for predicting protein three-dimensional structures There are only a few cases where the predicted structure is very different from the experimental one. For example, in ten cases with Delta-α > 30°, most of them have large Delta-ds (Figure 1D) the average pLDDT values of the residues bridged by the ChB are <90 indicating that side-chains might not be predicted reliably suggesting that predictions should be treated with caution Only three predicted structures (PDB: 3soj despite Delta-α >30° Figure 2C shows one of the rare examples where AlphaFold DB models seem to be completely inadequate The local stereochemistry–this a partially disordered part of the protein–is wrong the cysteine 4 is misplaced and the disulphide bonds is broken We conclude that computational models produced with AlphaFold2 and stored in AlphaFold DB are accurate–we note that for these proteins a high-resolution crystal structure is available in the PDB they show contacts between sulphur atoms of disulphide bridges and protein nucleophilic atoms that are comparable to the experimental ones the accuracy of the computational models is insufficient to allow the detection of ChBs The original contributions presented in the study are included in the article/Supplementary Material further inquiries can be directed to the corresponding authors OC and KD-C contributed to the conception and design of the study OC performed the statistical analysis and wrote the first draft of the manuscript All authors contributed to the article and approved the submitted version The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1155629/full#supplementary-material Definition of the chalcogen bond (IUPAC Recommendations 2019) CrossRef Full Text | Google Scholar A structural biology community assessment of AlphaFold2 applications PubMed Abstract | CrossRef Full Text | Google Scholar The epigenetic dimension of protein structure PubMed Abstract | CrossRef Full Text | Google Scholar PubMed Abstract | CrossRef Full Text | Google Scholar The protein Data Bank: A computer-based archival file for macromolecular structures PubMed Abstract | CrossRef Full Text | Google Scholar Improved prediction of protein-protein interactions using AlphaFold2 PubMed Abstract | CrossRef Full Text | Google Scholar Can AlphaFold2 predict the impact of missense mutations on structure PubMed Abstract | CrossRef Full Text | Google Scholar Interplay between hydrogen and chalcogen bonds in cysteine PubMed Abstract | CrossRef Full Text | Google Scholar Chalcogen bonds involving selenium in protein structures PubMed Abstract | CrossRef Full Text | Google Scholar Survey of the Intermolecular disulfide bonds observed in protein crystal structures deposited in the protein data bank PubMed Abstract | CrossRef Full Text | Google Scholar Chalcogen bonds formed by protein sulfur atoms in proteins A survey of high-resolution structures deposited in the protein data bank CrossRef Full Text | Google Scholar Online_DPI: A web server to calculate the diffraction precision index for a protein structure CrossRef Full Text | Google Scholar (2022) Protein complex prediction with AlphaFold-Multimer CrossRef Full Text | Google Scholar CD-HIT: Accelerated for clustering the next generation sequencing data PubMed Abstract | CrossRef Full Text | Google Scholar AlphaFill: Enriching AlphaFold models with ligands and cofactors PubMed Abstract | CrossRef Full Text | Google Scholar Evaluation of AlphaFold2 structures as docking targets PubMed Abstract | CrossRef Full Text | Google Scholar Computed structures of core eukaryotic protein complexes PubMed Abstract | CrossRef Full Text | Google Scholar Highly accurate protein structure prediction with AlphaFold PubMed Abstract | CrossRef Full Text | Google Scholar RoseTTAFold and modeller: A case study involving the use of G-protein-coupled receptors CrossRef Full Text | Google Scholar Cd-Hit: A fast program for clustering and comparing large sets of protein or nucleotide sequences PubMed Abstract | CrossRef Full Text | Google Scholar (2022) The protein-folding problem: Not yet solved PubMed Abstract | CrossRef Full Text | Google Scholar AI-based structure prediction empowers integrative structural analysis of human nuclear pores PubMed Abstract | CrossRef Full Text | Google Scholar Pak, M. A., Markhieva, K. A., Novikova, M. S., Petrov, D. S., Vorobyev, I. S., Maksimova, E. S., et al. (2021). Using AlphaFold to predict the impact of single mutations on protein stability and function. Prepr. bioRxiv. Available at: https//www.biorxiv.org/content/101101/20210919460937v1 Google Scholar The origin of chalcogen-bonding interactions PubMed Abstract | CrossRef Full Text | Google Scholar PubMed Abstract | CrossRef Full Text | Google Scholar PubMed Abstract | CrossRef Full Text | Google Scholar Sanchez Rodriguez Using deep-learning predictions of inter-residue distances for model validation CrossRef Full Text | Google Scholar The chalcogen bond in crystalline solids: A world parallel to halogen bond PubMed Abstract | CrossRef Full Text | Google Scholar Assessing PDB macromolecular crystal structure confidence at the individual amino acid residue level PubMed Abstract | CrossRef Full Text | Google Scholar Protein model refinement for cryo-EM maps using AlphaFold2 and the DAQ score CrossRef Full Text | Google Scholar Tunyasuvunakool Highly accurate protein structure prediction for the human proteome PubMed Abstract | CrossRef Full Text | Google Scholar AlphaFold protein structure database: Massively expanding the structural coverage of protein-sequence space with high-accuracy models PubMed Abstract | CrossRef Full Text | Google Scholar PubMed Abstract | CrossRef Full Text | Google Scholar The universal protein resource (UniProt): An expanding universe of protein information PubMed Abstract | CrossRef Full Text | Google Scholar wwPDB Consortium (2019) Protein Data Bank: The single global archive for 3D macromolecular structure data PubMed Abstract | CrossRef Full Text | Google Scholar AlphaPulldown-a python package for protein-protein interaction screens using AlphaFold-Multimer PubMed Abstract | CrossRef Full Text | Google Scholar Citation: Carugo O and Djinović-Carugo K (2023) Automated identification of chalcogen bonds in AlphaFold protein structure database files: is it possible? Received: 31 January 2023; Accepted: 26 June 2023;Published: 06 July 2023 Copyright © 2023 Carugo and Djinović-Carugo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use *Correspondence: Oliviero Carugo, b2xpdmllcm8uY2FydWdvQHVuaXZpZS5hYy5hdA== Kristina Djinović-Carugo, a3Jpc3RpbmEuZGppbm92aWNAdW5pdmllLmFjLmF0 Nearly a third of renal cell carcinoma patients develop metastatic disease oncologists have struggled to understand why.Credit: Rasi Bhadramani/ Getty Images Most patients with renal cell carcinoma (RCC) with slow-growing tumours that typically respond to treatment But nearly a third will develop aggressive metastatic disease from which the likelihood of surviving five years past diagnosis can be as low as 15% Clinicians have struggled to understand why New research from a team led by oncologist Giannicola Genovese and postdoctoral fellow Luigi Perelli at The University of Texas MD Anderson Cancer Center in collaboration with Alessandro Carugo at IRBM in Italy now reveals the series of genetic changes that make relatively indolent tumours malignant (Perelli The Genovese lab develops experimental systems for reconstructing the natural history of cancer “We genetically engineer the key drivers of the disease in experimental model systems — mostly the mouse — and then follow tumour evolution,” says Genovese “We’re trying to understand how tumours become aggressive the team engineered mouse models in which a subset of kidney cells were genetically manipulated to include various mutations with an established link to aggressive RCC One change proved particularly problematic the deletion of a section of mouse chromosome 4 that contains two genes that regulate cell division gave rise to extremely aggressive cancers Genovese and colleagues determined that this disruption creates conditions that lead to the acquisition of further genomic damage cells found in metastatic growths routinely lost a sizable segment of a second mouse chromosome containing genes that allow cells to respond to immune signaling molecules known as interferons Genovese and Perelli believe this deletion may make tumour cells resistant to the immune response allowing the cancer to progress more rapidly the researchers were able to map all the changes they observed in mice back to equivalent chromosomal disruptions in patient-derived RCC tumours “It parallels with the human data,” says Perelli This highlights the relevance of a well-designed model for studying human cancers and Genovese’s team is now looking into opportunities to apply their findings to improve care for RCC patients checkpoint inhibitor drugs that dramatically boost the anti-tumour immune response can be highly effective against RCC but less than half of patients respond meaningfully to this treatment The group is looking into how the disruption of the interferon response observed in this study might influence the success or durability of checkpoint inhibitor treatment They also see opportunities to use insights from this mouse model to guide how clinicians approach early-stage tumours based on the presence or absence of genomic features associated with metastasis “This will be very important for teasing out patients who should be treated less aggressively right now and to identify those who may need additional treatment.” This same approach to modelling tumour evolution could also illuminate distinctive roads to progression followed by other tumour types and Genovese says his team is now looking at the extent to which these routes overlap or diverge in pancreatic cancer as well as rarer forms of kidney cancer “We're working with clinicians to model those tumour types and understand if there are vulnerabilities that can be exploited.” To read the full paper in Nature Cancer, click here. Volume 3 - 2023 | https://doi.org/10.3389/fruro.2023.1335414 This article is part of the Research TopicModelling the Impact of Medical Devices in the Urinary Tract: From Bench to BedsideView all 4 articles Ureteral stents are hollow tubes that are inserted into the ureter to maintain the flow of urine from the kidney to the bladder the use of these indwelling stents is associated with potential complications an organized consortium of bacterial species embedded within a self-producing extracellular matrix can attach to the outer and inner surfaces of ureteral stents encrustation - defined as the buildup of mineral deposits on the stent surface - can occur independently or in parallel with biofilm formation Both phenomena can cause stent obstruction which can lead to obstructive pyelonephritis and make stent removal difficult Understanding the influence of flow on the development of biofilm and encrustation and the impact of small mechanical environmental changes (e.g. wall shear stress distribution) is key to improve the long-term performance of stents Identifying the optimal stent properties to prevent early bacterial attachment and/or crystal deposition and their growth would represent a breakthrough in reducing biofilm-/encrustation-associated complications This review identifies the most prevalent bacterial strains and crystal types associated with ureteral stents and the process of their association with the stent surface which often depends on patient comorbidities we focus on the often-overlooked role of fluid dynamics on biofilm and encrustation development in ureteral stents from micro- to macro-scale) with the aim of providing a knowledge base to inform the development of safer and more effective ureteral stents We refer to encrustation as the build-up of mineral deposits/crystals on the surface of ureteral stents or within the matrix of a biofilm translating computational and experimental findings into actual stent improvements is a challenging endeavor that is often ignored or underestimated This comprehensive review aims to introduce and discuss the methods used to identify and manage biofilm and encrustation in ureteral stents It first describes the bacterial species involved in biofilm formation and the interplay between biofilm and encrustation It subsequently discusses ways to prevent or minimize the growth of biofilms on ureteral stents that are specifically optimized based on fluid dynamics These include methods based on computer simulations and experimental tests to spatially resolve relevant flow metrics particularly WSS at different dimensional scales it identifies a range of WSS values acting on the stent surface and elucidates the impact of WSS on biofilm formation and crystal deposition This finding underlies the necessity to explore new techniques to prevent the conditioning film from forming in the first place in addition to preventing bacterial attachment to the stent material While numerous factors contribute to biofilm development, bacteria remain the primary cause of its formation (a summary on biofilm formation in ureteral stents is provided in Box 1). Biofilms often comprise a combination of bacterial species, which compete and/or cooperate to form multi-species biofilms (51) To effectively combat biofilm formation in ureteral stents it is crucial to identify and understand the specific bacterial species that thrive in this environment Summary of biofilm formation in ureteral stents • Ureteral stents are susceptible to biofilm development • Conditioning film forms from proteins and polysaccharides on stent surfaces • Bacteria adhere to stents and the conditioning film • Bacteria use flagella for strong catch bonds to resist shear forces • EPS matrix provides protection and nutrients to bacteria • Biofilms are usually flat under high-shear stress and thick under low-shear stress Influence of bacterial presence on crystal types in ureteral stents • Crystal type on ureteral stents can vary based on the bacterial presence calcium oxalate monohydrate crystals are the main crystals in ureteral stents struvite and hydroxyapatite crystals predominate • Certain bacteria accelerate stent encrustation these studies highlight the significance of the interplay between bacteria and crystal deposition in the development of encrustation Bacteria-induced encrustation forms more rapidly and is often thicker Since bacteria survive within the encrusting material effective removal of encrusting fragments is key to preventing recurrence Furthermore, while some antimicrobial coatings can reduce biofilm formation on polymeric stents, thus enhancing their durability, the effectiveness of these coatings can be compromised by the conditioning film and subsequent crystal deposition (57). Since crystals have a higher affinity to various components of the conditioning film than to the stent material itself (39) it is crucial to address the problem of crystal deposition in parallel with biofilm reduction Detection and analysis of bacterial colonization and encrustation in ureteral stents • Sterile urine cultures may underestimate the bacterial colonization on stents • Stent culture after vortexing and sonication is more effective than standard urine culture • Techniques like 16S rRNA sequencing and SEM improve understanding of stent microbiomes • SEM with EDX is commonly used to identify crystal types in stents and µCT are used to quantify stent encrustation Besides identifying the crystal type, some studies also used SEM as a qualitative method to estimate the extent of encrustation on different sections of ureteral stents. An in vitro study performed by Cauda et al. (76) used this technique to compare different stent materials qualitatively identifying those developing less encrustation and the percentage of the analyzed ureteral stents associated with the presence of bacteria (colonization rate) and research challenges in ureteral stent studies • Escherichia coli and Enterococcus faecalis and Pseudomonas aeruginosa are also frequently found on ureteral stents • Longer indwelling time increases microbial colonization of stents • Bacterial load and encrustation are normally low under six weeks • The high variability of the research methods used in the field complicates the Interpretation of the results and the drawing of clear conclusions Table 1 Characteristics of ureteral stents and associated bacterial colonization it emerges that Escherichia coli is consistently identified as the predominant bacterial species found in ureteral stents retrieved from patients are also identified in most of these studies Staphylococcus epidermis is the most prevalent one more than five bacterial species were found in most studies the identification of a ‘safe indwelling time’ is challenging due to the independent nature of most studies and variations in patients’ conditions and comorbidities the integration of obstruction monitoring methods with quantification techniques holds the potential to determine the appropriate time for stent removal and potentially identify regions of the stent that are most problematic Ureteral stents can exhibit different encrustation and biofilm patterns depending on the specific segment of the stent (a summary of the distribution patterns is provided in Box 5) A stent can be divided into its two pigtails (located in the renal pelvis and the bladder respectively) and a straight central segment This straight part can be subdivided into three regions: i) proximal section (near the ureteropelvic junction or UPJ) ii) middle section (within the middle ureter) and iii) distal section (near the ureterovesical junction or UVJ) Location matters in ureteral stent bacterial colonization and encrustation • Microbial colonization is more pronounced in the distal and proximal sections of ureteral stents • Encrustations are typically more visible in the proximal section • Colonization and encrustation patterns can vary based on the patient type • Future consideration should be given to developing patient-specific stents which is in disagreement with most other studies in this area the possibility of developing patient-specific stents should be considered in the future future scientific efforts should focus on combining information on stent sections that are more prone to develop encrustation (also based on patient-specific conditions) with the analysis of local fluid dynamics in these regions and WSS ranges observed on both the internal and external walls of the stent and side holes Table 2 Summary of study parameters and wall shear stress value ranges in ureteral stents Firstly, microfluidic Stent-on-Chip (SOC) models (32) were employed to quantify the WSS distribution in a ureteral stent and the effect of side hole shape and stent wall thickness on the WSS field these values should be interpreted cautiously as the intrinsic simplifications of microfluidic models prevent them from reproducing the full complexity of the stented ureter (e.g. In an attempt to provide more physiological WSS values, Mosayyebi et al. (33) used CFD to compute WSS in a tapered stented ureter that mimicked the ureter of a pig, and identified variations in WSS levels depending on the longitudinal position along the ureter. They also evaluated WSS acting over the inner wall of side holes, as shown in Table 2 the ureteropelvic junction (UPJ) had a larger diameter than the ureterovesical junction (UVJ) This configuration resulted in elevated WSS levels within the mid and distal sections A comprehensive CFD and particle image velocimetry study by Zheng et al. (2) employed a full-scale model to assess WSS distribution over the inner walls of side holes the ureter was modeled as an unobstructed tube with a constant diameter The mean WSS was higher in the first and last side holes of the straight part of the stent while the side holes in the middle section exhibited lower WSS These differences may be attributed to the passive nature of side holes in the middle section since there is no significant flow exchange through these holes (‘passive’ side holes) whereas the first and last side holes promote such flow exchange (‘active’ side holes) Their findings also indicated that smaller side holes experienced significantly greater WSS levels providing a potential means to mitigate encrustation and biofilm formation in these regions of the stent due to the correlation between these phenomena and WSS Vogt et al. (31) conducted a small-scale CFD simulations on DJ ureteral stents emphasizing the importance of utilizing smooth stent designs since grooves and imperfections on the stent surface lead to stagnation zones (with low WSS levels) which can act as traps for crystals and promote encrustation corresponding to a wide range of fluid shear forces Wall shear stress values across various ureteral stent locations • Stent design affects shear stress and encrustation patterns  o Proximal side holes show very low wall shear stress (WSS) below 0.001 Pa  o Distal side holes have higher WSS  o WSS varies from 10-6 Pa to 0.19 Pa on the stent’s internal wall accurate WSS measurements in stented ureters depend on various factors and the study design none of the studies reviewed have considered the impact of reflux which substantially influence WSS values and add complexities to a model To better understand the true impact of WSS on ureteral stents future research must diligently consider and incorporate these factors into their models (whether in vitro and/or in silico) At the micro-scale, the relationship between WSS and bacterial interactions with surfaces has been extensively explored in various domains (115) (see Section 4). However, this critical association remains largely unexplored in the context of ureteral stents. De Grazia et al. (34) conducted experimental investigations using a microfluidic SOC approach were the primary sites of bacterial attachment followed by side holes and the intra-luminal surface This study marked the first attempt to examine bacterial attachment on ureteral stent architectures The study revealed a correlation between bacteria coverage area Pseudomonas fluorescens was used as a bacterial model which is not a highly prevalent bacterium in the urinary tract and may limit the clinical relevance of the study and its generalization streamlined) side holes could reduce particle deposits in stents as a result of increased WSS levels Besides the SOC study, Mosayyebi et al. (33) also investigated the relationship between WSS and crystal deposition at the macroscale level Their study stands out as the only macroscale model investigating this relationship It was observed that the accumulation of particles was more pronounced in the side holes situated in the proximal region of the stent compared to those in the distal region as expected by the higher shear stress levels in the former their full-scale CFD study showed that during the voiding stage the average WSS at side holes located in the distal section is larger than in the middle and proximal sections could be perfused with both artificial urine and urinary bacteria Significance of wall shear stress in bacterial attachment and encrustation in ureteral stents • Experimental studies showed that areas with low WSS (below 0.04 Pa) are primary sites for bacterial attachment • Experimental studies showed that particle deposition is more common in areas with lower WSS thinner walls and streamlined side holes) influences shear stress and encrustation patterns • Particle accumulation is higher in proximal side holes (WSS < 0.001 Pa) compared to distal ones (WSS > 0.01 Pa) • Experimental studies showed that vesicoureteral reflux (VUR) reduces encrustation • There is a lack of comprehensive studies integrating flow dynamics with biofilm formation and encrustation Valuable insights into bacterial attachment and wall shear stress from related research • Despite limited experimental studies directly linking bacterial attachment and WSS in the urinary tract valuable insights are derived from other research areas • Microfluidics is used to explore how shear stress affects biofilm formation • Experimental studies showed that bacterial colonization decreases with WSS > 0.02 Pa • Different bacteria can exhibit different attachment rates depending on the level of WSS the development of fluid-mechanical-based stent designs maximizing WSS to reduce bacterial attachment Table 3 provides an overview of microfluidic studies that have established a correlation between WSS levels (within the range of values observed in ureteral stents) and the behavior of bacteria that are prevalent in the urinary tract The studies were identified and selected using the Scopus database with specific keywords (“*” AND (“shear stress” OR “biofilm formation” OR “biofilm growth”) AND (“microfluidics” OR “microfluidic” OR “micro-fluidic” OR “flow chamber”) where * represent all the bacteria selected above) Table 3 The impact of wall shear stress on biofilm formation in microfluidic experiments coli to bladder T24 transitional cells and type IV collagen reaches its highest level under minimal shear stress (0.01 Pa) conditions even when subjected to higher shear stress (0.59 Pa) coli cells remain attached after binding to host cells or collagen though detachment rates escalate with the magnitude of shear stress the reported WSS in the distal section of the ureteral stent spans from 10-2 to 10-1 Pa those in the proximal section range from 10-5 to 10-4 Pa Considering several reports suggesting the complete suppression or substantial reduction of biofilm thickness beyond 0.02 Pa this becomes a possible reason for the observed phenomenon a flawless stent design has not yet been achieved and none of the existing stents can completely prevent the formation of biofilms or encrustation on their surface Advanced materials and coatings: Continued exploration of novel materials and coatings could provide solutions to reduce encrustation and biofilm formation Researchers might focus on materials that minimize friction during placement and removal while maximizing long-term performance Biofilm and encrustation control: Exploring innovative strategies such as targeted antimicrobial/encrustation therapies or biofilm/encrustation-disrupting techniques could offer new avenues to prevent stent-related complications Patient-specific and practice-driven design: For short-term stenting (less than two weeks) clinicians may prefer ureteral stents with better cost-effectiveness due to the small likelihood of encrustation and biofilm formation a special focus on the current stent designs is crucial since small changes already proved beneficial in reducing attachment Special care should be taken in the pigtail sections Micro- and macro-scale research areas are also presented followed by a description of the WSS ranges within ureteral stents Innovative stents: New strategies and technologies are emerging such as smart stents that incorporate sensors enabling real-time monitoring of urinary tract conditions This could allow for early detection of complications and timely interventions Multidisciplinary approach: Encouraging collaboration between diverse disciplines will be essential to tackle the challenges associated with stent complications Regulatory Guidelines: Collaboration between researchers and regulatory bodies could lead to the development of standardized guidelines for stent testing Writing – review & editing The author(s) declare financial support was received for the research The work is supported by the Swiss National Science Foundation (SNSF DC is co-inventor on patent application PCT/GB2018/052522 “Stent with streamlined side holes” The author(s) declared that they were an editorial board member of Frontiers This had no impact on the peer review process and the final decision Clinical use of long-term indwelling silicone rubber ureteral splints inserted cystoscopically Fluid mechanical performance of ureteral stents: The role of side hole and lumen size New ureteral stent design does not improve patient quality of life: a randomized Latest advancements in ureteral stent technology Transl Androl Urol (2019) 8(Suppl 4):S436–s441 Quantitative evaluation of encrustations in double-j ureteral stents with micro-computed tomography and semantic segmentation Patient experiences and preferences with ureteral stent removal The change of antibiotic susceptibility in febrile urinary tract infection in childhood and adolescence during the last decade Investig Clin Urol (2022) 63(1):99–106 Ureteral stents: morbidity and impact on quality of life PubMed Abstract | Google Scholar Bilateral ureteral stent removal after 15 years: A case report J Minim Invasive Gynecol (2018) 25(5):920–2 Google Scholar Advances in ureteral stent design and materials PubMed Abstract | CrossRef Full Text | Google Scholar Experience with new double j ureteral catheter stent PubMed Abstract | CrossRef Full Text | Google Scholar A randomized clinical trial evaluating the short-term results of ureteral stent encrustation in urolithiasis patients undergoing ureteroscopy: micro-computed tomography evaluation Reducing deposition of encrustation in ureteric stents by changing the stent architecture: A microfluidic-based investigation Investigating the flow dynamics in the obstructed and stented ureter by means of a biomimetic artificial model Comparative assessment of ureteral stent biomaterial encrustation Comparison of silicone versus polyurethane ureteral stents: a prospective controlled study Heparin coating on ureteral double j stents prevents encrustations: an in vivo case study Diamond-like carbon coatings on ureteral stents–a new strategy for decreasing the formation of crystalline bacterial biofilms Use of triclosan-eluting ureteral stents in patients with long-term stents Numerical analysis of the urine flow in a stented ureter with no peristalsis BioMed Mater Eng (2015) 26(Suppl 1):S215–23 Influence of single stent size and tandem stents subject to extrinsic ureteral obstruction and stent occlusion on stent failure Analysis of urine flow in three different ureter models Comput Math Methods Med (2017) 2017:5172641 Numerical simulation of the urine flow in a stented ureter PubMed Abstract | CrossRef Full Text | Google Scholar Arrangement of side holes in a double j stent for high urine flow in a stented ureter J Mechanical Sci Technol (2020) 34(2):949–54 Ureteric stents: investigating flow and encrustation Proc Inst Mech Eng H (2008) 222(4):551–61 PubMed Abstract | CrossRef Full Text | Google Scholar The effect of ureteric stents on urine flow: reflux Fluid structural analysis of urine flow in a stented ureter Comput Math Methods Med (2016) 2016:5710798 Failure of ureteral stents subject to extrinsic ureteral obstruction and stent occlusions Int Urol Nephrol (2021) 53(8):1535–41 Investigating the encrustation of reinforced ureteral stents by computational flow dynamic simulations PubMed Abstract | CrossRef Full Text | Google Scholar Particle accumulation in ureteral stents is governed by fluid dynamics: In vitro study using a "Stent-on-Chip" model The accumulation of particles in ureteric stents is mediated by flow dynamics: Full-scale computational and experimental modeling of the occluded and unoccluded ureter A microfluidic-based investigation of bacterial attachment in ureteral stents The binding of urinary components and uropathogens to a silicone latex urethral catheter Cells Materials (1992) 2(3):253-–60 Google Scholar Understanding the fundamental basis for biofilm formation on plastic surfaces: Role of conditioning films Observations of bacterial biofilm on ureteral stent and studies on the distribution of pathogenic bacteria and drug resistance Understanding urinary conditioning film components on ureteral stents: profiling protein components and evaluating their role in bacterial colonization Uropathogens preferrentially interact with conditioning film components on the surface of indwelling ureteral stents rather than stent material Presence of five conditioning film proteins are highly associated with early stent encrustation Adhesins involved in attachment to abiotic surfaces by gram-negative bacteria doi: 10.1128/microbiolspec.MB-0018-2015 CrossRef Full Text | Google Scholar Emerg Infect Dis (2002) 8(9):881–90 PubMed Abstract | CrossRef Full Text | Google Scholar Binding of the bacterial adhesin FimH to its natural multivalent high-mannose type glycan targets Catch-bond mechanism of force-enhanced adhesion: counterintuitive Cell Host Microbe (2008) 4(4):314–23 Not just going with the flow: The effects of fluid flow on bacteria and plankton Annu Rev Cell Dev Biol (2019) 35:213–37 doi: 10.1146/annurev-cellbio-100818-125119 Effect of shear stress and growth conditions on detachment and physical properties of biofilms The EPS matrix: the "house of biofilm cells" PubMed Abstract | CrossRef Full Text | Google Scholar Biofilms: an emergent form of bacterial life Nat Rev Microbiol (2016) 14(9):563–75 and hydrodynamic conditions on bacterial surface sensing and their initial adhesion Front Bioengineering Biotechnol (2021) 9:643722 Mechanics of biofilms formed of bacteria with fimbriae appendages PubMed Abstract | CrossRef Full Text | Google Scholar Current understanding of multi-species biofilms Double-j ureteric stent encrustations: clinical study on crystal formation on polyurethane stents Study on concretions developed around urinary catheters and mechanisms of renal calculi development Uric acid as inducer of calcium oxalate crystal development Scand J Urol Nephrol (2007) 41(1):26–31 Risk of symptomatic kidney stones during and after pregnancy Am J Kidney Dis (2021) 78(3):409–17 Update on biofilm infections in the urinary tract The encrustation and blockage of long-term indwelling bladder catheters: a way forward in prevention and control Extraction and quantification of biofilm bacteria: Method optimized for urinary catheters Colonization of internal ureteral stent and bacteriuria PubMed Abstract | CrossRef Full Text | Google Scholar Improved detection of microbial ureteral stent colonisation by sonication Bacterial colonization of double j stents and bacteriuria frequency Kaohsiung J Med Sci (2013) 29(12):658–61 Microbial adhesion and biofilm formation on ureteral stents in vitro and in vivo Clinical significance of DJ stent culture in patients with indwelling ureteral stents prior to endourological intervention Microbial ureteral stent colonization in renal transplant recipients: frequency and influence on the short-time functional outcome Transpl Infect Dis (2012) 14(1):57–63 doi: 10.1111/j.1399-3062.2011.00671.x Ureteral stent colonization and urinary tract infection in children undergoing minimally invasive pyeloplasty Eur J Pediatr Surg (2023) 33(01):047–52 Characteristics of bacterial colonization after indwelling ureteral stents in urinary stone patients with diabetes mellitus and chronic kidney disease Evaluation of bacterial colonization and bacteriuria secondary to internal ureteral stent Mymensingh Med J (2010) 19(3):366–71 PubMed Abstract | Google Scholar Predictive value of urine cultures in evaluation of bacterial colonization of ureteral stents Mymensingh Med J (2012) 21(2):300–5 PubMed Abstract | Google Scholar Predictive value of urinary cultures in assessment of microbial colonization of ureteral stents Ureteral stent microbiota is associated with patient comorbidities but not antibiotic exposure Incidence of bacterial colonisation after indwelling of double-j ureteral stent Arch Ital Urol Androl (2016) 87(4):291–4 Encrustations on ureteral stents from patients without urinary tract infection reveal distinct urotypes and a low bacterial load Comparison of four DNA extraction methods for comprehensive assessment of 16S rRNA bacterial diversity in marine biofilms using high-throughput sequencing Detection of microbial colonization of the urinary tract of patients prior to secondary ureterorenoscopy is highly variable between different types of assessment: results of a prospective observational study Surface analysis of ureteral stent before and after implantation in the bodies of child patients Ureteral double-j stents performances toward encrustation after long-term indwelling in a dynamic in vitro model J BioMed Mater Res B Appl Biomater (2017) 105(8):2244–2253 Development of a model for assessment of biomaterial encrustation in the upper urinary tract Assessment of encrustations on polyurethane ureteral stents Surface- treated pellethanes: Comparative quantification of encrustation in artificial urine solution Severely encrusted polyurethane ureteral stents: management and analysis of potential risk factors In vitro encrustation of a semi-permanent polymer-covered nitinol ureter stent: an artificial urine model Chemical and mineralogical analysis of ureteral stent encrustation and associated risk factors Encrustation of urologic double pigtail catheters-an ex vivo optical coherence tomography (OCT) study An in vitro bladder model with physiological dynamics: Vesicoureteral reflux alters stent encrustation pattern Front Bioengineering Biotechnol (2022) 10:1028325 Fluid mechanical modeling of the upper urinary tract Characteristics of bacterial colonization after indwelling double-j ureteral stents for different time duration Characterization of biofilm and encrustation on ureteric stents in vivo doi: 10.1111/j.1464-410x.1994.tb07557.x Characteristics of bacterial colonization and urinary tract infection after indwelling of double-j ureteral stent Evaluation of ureteral stent colonization in live-donor renal transplant recipients doi: 10.1016/j.transproceed.2017.01.004 Bacterial colonization of ureteral double-j stents in patients with negative urine culture Google Scholar Risk factors for lower urinary tract infection and bacterial stent colonization in patients with a double j ureteral stent Bacteriology of urinary tract infection associated with indwelling j ureteral stents Study of risk factors for urinary colonization in patients with the double j catheter CrossRef Full Text | Google Scholar Relationship between urinary tract infection and self-retaining double-j catheter colonization Does bacterial colonization influence ureteral stent-associated morbidity Analysis of bacterial biofilm formation in patients with malignancy undergoing double j stent indwelling and its influencing factors Infection Drug Resistance (2021) 14:5209–17 Int Urol Nephrol (2002) 34(2):199–202 Is there a role for double j stent culture in contemporary urology Urologia Internationalis (2018) 100(2):203–8 The risk of bacteriuria and ureteric stent colonization in immune-compromised patients with double j stent insertion Int Urol Nephrol (2010) 42(2):343–7 Epidemiology and risk factors for ureteral stent-associated urinary tract infections in non-transplanted renal patients: a systematic review of the literature Conditioning film deposition on ureteral stents after implantation PubMed Abstract | CrossRef Full Text | Google Scholar PubMed Abstract | CrossRef Full Text | Google Scholar Biofilm formation on ureteral stents – incidence and coloring: morbidity related to indwelling times Polyurethane internal ureteral stents in treatment of stone patients: morbidity related to indwelling times Ureteral stenting and urinary stone management: a systematic review Real-time monitoring of urinary encrustation using a quartz crystal microbalance Analytical Chem (2018) 90(3):1531–5 CrossRef Full Text | Google Scholar Electromechanically functionalized ureteral stents for wireless obstruction monitoring ACS Biomater Sci Eng (2023) 9(7):4392–403 Ureteral stents harbor complex biofilms with rich microbiome-metabolite interactions Management strategies and results for severely encrusted retained ureteral stents PubMed Abstract | CrossRef Full Text | Google Scholar Effects of rhamnolipids and shear on initial attachment of pseudomonas aeruginosa PAO1 in glass flow chambers Environ Sci pollut Res (2010) 17(9):1529–38 Quorum sensing is crucial to escherichia coli O157:H7 biofilm formation under static or very slow laminar flow conditions Effect of shear stress on the formation of bacterial biofilm in a microfluidic channel CrossRef Full Text | Google Scholar A shear stress micromodel of urinary tract infection by the escherichia coli producing dr adhesin CFD study on vesicoureteral reflux in the urinary tract with double j stent doi: 10.1016/j.compbiomed.2022.105456 Formation and characteristics of laminar vortices in microscale environments within an obstructed and stented ureter: a computational study in: Proceedings of the 18th International Conference on Miniaturized Systems for Chemistry and Life Sciences Google Scholar Bacterial aggregation and biofilm formation in a vortical flow PubMed Abstract | CrossRef Full Text | Google Scholar Structural and population characterization of MrkD A microfluidic approach for studying microcolonization of escherichia coli O157:H7 on leaf trichome-mimicking surfaces under fluid shear stress Biotechnol Bioengineering (2022) 119(6):1556–66 Bacterial biofilm under flow: First a physical struggle to stay Effect of shear stress on the reduction of bacterial adhesion to antifouling polymers Bioinspiration Biomimetics (2018) 13(6):065001 Combining chemistry and topography to fight biofilm formation: Fabrication of micropatterned surfaces with a peptide-based coating Colloids Surfaces B: Biointerfaces (2020) 196:111365 The effects of surface properties on escherichia coli adhesion are modulated by shear stress Colloids Surf B Biointerfaces (2014) 123:1–7 Influence of shear on microbial adhesion to PEO-brushes and glass by convective-diffusion and sedimentation in a parallel plate flow chamber Colloids Surfaces B: Biointerfaces (2005) 46(1):1–6 Microfluidic devices for studying growth and detachment of staphylococcus epidermidis biofilms Biomed Microdevices (2008) 10(4):489–98 Bacterial adhesion and growth on a polymer brush-coating doi: 10.1016/j.biomaterials.2008.07.014 Real-time monitoring of mono- and dual-species biofilm formation and eradication using microfluidic platform A novel planar flow cell for studies of biofilm heterogeneity and flow-biofilm interactions Biotechnol Bioeng (2011) 108(11):2571–82 Modeling growth and quorum sensing in biofilms grown in microfluidic chambers Ann BioMed Eng (2009) 37(6):1206–16 Characterizing the adhesion of motile and nonmotile escherichia coli to a glass surface using a parallel-plate flow chamber Biotechnol Bioengineering (2002) 78(2):179–89 CrossRef Full Text | Google Scholar Physico-chemical characterization of caesium and strontium using fluorescent intensity of bacteria in a microfluidic platform The effects of fluid composition and shear conditions on bacterial adhesion to an antifouling peptide-coated surface CrossRef Full Text | Google Scholar Micro- and macro-flow systems to study escherichia coli adhesion to biomedical materials CrossRef Full Text | Google Scholar A selection of platforms to evaluate surface adhesion and biofilm formation in controlled hydrodynamic conditions A novel medical device coating prevents staphylococcus aureus biofilm formation on medical device surfaces Submicron-textured biomaterial surface reduces staphylococcal bacterial adhesion and biofilm formation Acta Biomaterialia (2012) 8(1):72–81 Shear stress effects on bacterial adhesion and leukocyte oxidative capacity on a polyetherurethane J Biomed Materials Res (1999) 46(4):511–519 doi: 10.1002/(SICI)1097-4636(19990915)46:4<511::AID-JBM9>3.0.CO;2-M Quantitative comparison of shear-dependent staphylococcus aureus adhesion to three polyurethane ionomer analogs with distinct surface properties J Biomed Materials Res (1997) 36(2):152–62 doi: 10.1002/(SICI)1097-4636(199708)36:2<152::AID-JBM3>3.0.CO;2-J Bacterial adhesion to polyurethane surfaces in the presence of pre-adsorbed high molecular weight kininogen J Biomaterials Science Polymer Edition (1996) 7(9):769–80 CrossRef Full Text | Google Scholar Pathogenic distribution and drug resistance of bacterial biofilm on ureteral stent Chin J Tissue Eng Res (2020) 24(16):2556–60 CrossRef Full Text | Google Scholar Burkhard F and Clavica F (2024) The interplay between bacterial biofilms and wall shear stress in ureteral stents: a review across scales Received: 09 November 2023; Accepted: 19 December 2023;Published: 16 January 2024 Copyright © 2024 Amado, Zheng, Lange, Carugo, Waters, Obrist, Burkhard and Clavica. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY) *Correspondence: Pedro Amado, cGVkcm8ucGVyZWlyYUB1bmliZS5jaA== Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher. 94% of researchers rate our articles as excellent or goodLearn more about the work of our research integrity team to safeguard the quality of each article we publish. Takotsubo syndrome (TTS) is a transient heart disease that has been historically related to the occurrence of psychological (emotional) factors (“broken heart” syndrome). We aimed to conduct a systematic review analyzing the role of psychological factors in TTS. All studies on TTS and psychological factors from January 1991 through April 2019 were scrutinized according to the Cochrane Collaboration and the PRISMA statements. Selected studies were additionally evaluated for the Risk of Bias according to the Newcastle-Ottawa Scale (NOS). Fifteen case-control studies (by Mayo Clinic criteria) were finally selected. Most studies analyzed stressful life-events or trauma, although with conflicting findings, while a likely role of long-lasting psychological distress seemed to be a homogenous result. Among life-time psychopathology, only anxiety appeared to have a significant role. Some studies outlined a likely role of personality, but findings are conflicting. Our findings do not lead to any definitive assumption on the specific role of psychological factors in TTS, also for scant strong methodology of the most part of the studies. More studies with stronger research methodology are needed to better characterize psychological elements in TTS. Volume 10 - 2019 | https://doi.org/10.3389/fpsyg.2019.02742 Objective: Takotsubo syndrome (TTS) is a transient heart disease that has been historically related to the occurrence of psychological (emotional) factors (“broken heart” syndrome) We aimed to conduct a systematic review analyzing the role of psychological factors in TTS Methods: All studies on TTS and psychological factors from January 1991 through April 2019 were scrutinized according to the Cochrane Collaboration and the PRISMA statements Selected studies were additionally evaluated for the Risk of Bias according to the Newcastle-Ottawa Scale (NOS) Results: Fifteen case-control studies (by Mayo Clinic criteria) were finally selected Most studies analyzed stressful life-events or trauma while a likely role of long-lasting psychological distress seemed to be a homogenous result only anxiety appeared to have a significant role Some studies outlined a likely role of personality Conclusion: Our findings do not lead to any definitive assumption on the specific role of psychological factors in TTS also for scant strong methodology of the most part of the studies More studies with stronger research methodology are needed to better characterize psychological elements in TTS Takotsubo syndrome (TTS) is a form of transient heart failure syndrome often mimicking acute myocardial infarction. It is also known as stress cardiomyopathy, broken heart syndrome, or apical ballooning syndrome, and was firstly described in 1990 (Dote et al., 1991), even if some Authors dated it earlier (Wittstein, 2008) no systematic reviews have been realized analyzing the likely causative link between psychological factors and TTS with the existing reviews detecting other factors (e.g. the role of drugs or pathophysiologic mechanisms) The objective of the present study was to systematically review all the studies on psychological factors (as antecedents) (traumatic/stressful events psychopathology and personality) in patients with TTS diagnosis in order to understand their likely role in this syndrome We conducted a systematic review of the literature on psychological factors (psychopathologic disorders stressful life-events/psychological trauma All observational studies were included in the review by the ascertainment of a case-control study design The search was limited to English-written publications and to the period from January 1991 to April 2019 An additional analysis of the reference list in each selected paper was also performed • Studies with an analytical study design as defined by Grimes and Schulz (2002) (i.e. an observational study with a comparison or control group Both retrospective and perspective studies have been included to consider the highest number of studies • Diagnosis of TTS by the Mayo Clinic criteria or by the new TTS criteria (Parodi et al., 2014; Lyon et al., 2016) • Studies adopting standardized and validated tests • Studies written in English language • Studies with intra-group control (e.g. • Pharmacological and behavioral intervention trials • Number of subjects per group ≤ 5 Study selection was performed by two independent reviewers with research expertise in clinical psychology and cardiology (FG and FB) who assessed the relevance of the study for the objectives of this review This first round of selection was based on the title If the reviewers did not reach a consensus or the abstract did not contain sufficient information the full text was reviewed In the second phase (screening), full-text reports were evaluated to detect whether the studies met the inclusion criteria (Figure 1) PRISMA flow diagram of literature search and selection of publications all full-texts were retrieved and a final check was made to exclude papers not responding to inclusion/exclusion criteria and reaching the final consensus to decide the final number of studies to be selected A standardized data extraction form was prepared; data were independently extracted by two of the authors (FG and FB) and inserted in a study database. A process of discussion/consensus moderated by a third reviewer (SC) (Furlan et al., 2009) resolved discrepancies between reviewers For variables expressed as median and 25th–75th percentile it was not possible to calculate Cohen’s d Quality assessment of each of the included studies was evaluated following the Newcastle-Ottawa Scale (NOS) for case-control studies on a 9-star model (Wells et al., 2014) Studies scoring above the median NOS value were considered as high quality (low risk of bias) and those scoring below the median value were considered as low quality (high risk of bias) two reviewers (FB and FG) independently extracted relevant information and data from all eligible reports that met the above inclusion criteria We found 15 studies meeting inclusion criteria (Figure 1), for a total of 2581 subjects (1152 TTS patients; 1069 other heart diseases; 360 healthy controls). The description of the samples, psychological variables, study design, psychometric scales (tests, interview or retrospective medical records analysis), key findings, Cohen’s d or Odd ratios and 95% confidence intervals of selected studies are reported in Table 1 Of interest, the time elapsing between the cardiac event and the psychological assessment was extremely variable, ranging from hours (Delmas et al., 2013) to years (Goh et al., 2016) History of long-lasting psychological distress not temporally related to the cardiac events was evidenced in four studies (Delmas et al., 2013; Kastaun et al., 2014; Lacey et al., 2014; Rosman et al., 2017) Some studies differentiated emotional (27–38%) from physical (36–50%) triggers (Compare et al., 2014; Templin et al., 2015; Goh et al., 2016; Rosman et al., 2017) with 12–28% of patients with indeterminable reasons acute myocardial infarction with emotion triggers vs TTS without emotion trigger and found a significant prevalence only in the two groups with emotion triggers Other studies (Compare et al., 2014; Kastaun et al., 2014) did not find any role for psychopathology in the history of TTS patients Half of the studies reflected the median value (μ = 5), four were above it and three below (Table 2). Four studies were quoted as high quality (low risk of bias) by NOS (see Table 2) Ottawa-Newcastle risk of bias for case-control studies Although the role of psychological factors has been extensively studied in TTS only fifteen studies fulfilled the criteria to perform a systematic review because of the small number of selected studies and the heterogeneous methodology used for the psychological assessment (no studies shared the same psychological assessment tools) Most studies attempted to understand if stressful events (or trauma) could have a role (trigger) in TTS, but findings are conflicting. As suggested by the recent Expert Consensus Document on TTS (Jelena-Rima et al., 2018) one of the key questions to answer is which role triggering factors have in the stress response of the heart the etymology of “stress” cardiomyopathy requires specific attention for the role of psychological stressor as possible etiological factor our review does not allow any conclusion by this side The first point that warrants attention is the difference between studies drawing data from standardized psychometric tools or from retrospective assessment of medical records none allowed any conclusion toward a role for psychological trigger events in comparison to other cardiac events (control group) all studies that draw data from medical records (usually based on clinical interview at admission) evidenced a role for psychological trigger events compared to controls this opens both to methodological and clinical considerations the use of standardized measures of assessment would lead to strongest conclusions but in a direction making questionable the evidence of a role of psychological factors closely involved in the etiology of TTS homogenous clinical observations by medical records suggest implementing further case-control studies to support the role of psychological triggers in TTS The unanswered question is why a person develops TTS and another one other disorders Not finding any differences between TTS and myocardial infarction as regards to psychological triggers may delineate at some level the involvement of similar mechanisms Furthermore, homogeneous (albeit still limited) findings are driven from the analysis of studies evidencing a role for long-lasting psychological distress. A role for early traumatic psychological experience has been evidenced as predisposing factor for patients with cardiovascular diseases (Thurston et al., 2014, 2017; Bomhof-Roordink et al., 2015; Winning et al., 2016) and deserves more attention in future studies This finding may help to explain why only some patients experience a stressful event as trigger of TTS the comorbid occurrence of anxiety and depression should be considered as non-specific of TTS A final note on the unique study (Saffari et al., 2017) evidencing worsening of the quality of life and sexuality in TTS there is the impossibility to make a meta-analysis the small number of studies and the differences in the psychometric tools timing of observation and different design among different studies did not allow to pursuit the initial aim A or D personality) not totally supported by strong evidence (not included in DSM 5) so that any sound conclusion on this topic is not allowed and further research need to be addressed on this topic the rigor of methodology we relied on (Cochrane Collaboration and the PRISMA statements) allowed to get strong results and conclusion on the basis of our systematic review we cannot evidence a clear-cut role for psychological trauma preceding TTS onset but a possible role of long-lasting emotional distress we can suggest a role for life-time anxiety disorders (more than depression) but studies are needed to clarify if differences exist with other cardiac events we cannot conclude in the direction of specific patterns differentiating TTS from other cardiac disorders We need studies with stronger methodology addressing the involvement of emotional events by structured interviews conducted shortly after the onset of TTS The timing of interviewing patients should be carefully delineated (no more than 6 months) to avoid recall bias multicentre studies are warranted to recruit a large number of patients and increase sample size for this relatively rare entity The choice of an adequate control groups needs attention because one of the main questions is whether TTS actually differs from other cardiac disorders as regards to personality and comorbid psychopathologic disorders Finally, we stress the importance of a multidisciplinary approach to TTS; such an approach should involve a collaborative process between cardiologists and clinical psychologists from the diagnosis to treatment. Evidence are accumulating on the efficacy of psychological interventions for cardiac diseases (Richards et al., 2018) and SC contributed to the conception and design of the study FG organized the database and wrote the first draft of the manuscript FG and FB made the bibliographic research and selected papers for the systematic review (in case of doubt made confirmation with SC) FB and SC read and approved the submitted version of the manuscript All authors contributed to the manuscript revision ESC working group position paper on myocardial infarction with non-obstructive coronary arteries Google Scholar American Psychiatric Association [APA] Google Scholar Emotion risk-factor in patients with cardiac diseases: the role of cognitive emotion regulation strategies positive affect and negative affect (A Case-Control Study) Associations between life stress and subclinical cardiovascular disease are partly mediated by depressive and anxiety symptoms Google Scholar Myocardial lesions in victims of homicidal assaults without internal injuries depression and anxiety in takotsubo cardiomyopathy Statistical Power Analysis for the Behavioral Sciences Google Scholar Type D personality is associated with the developmentof stress cardiomyopathy following emotional triggers Google Scholar The role of emotional competence in takotsubo cardiomyopathy Stress-Induced cardiomyopathy and psychological wellbeing 1 year after an acute event Anxiety trait in patients with stress-induced cardiomyopathy: a case–control study Anxiodepressive disorders and chronic psychological stress are associated with tako-tsubo cardiomyopathy Type D personality-A potential risk factor refined PubMed Abstract | Google Scholar Myocardial stunning due to simultaneous multivessel coronary spasms: a review of 5 cases Posttraumatic stress disorder prevalence and risk of recurrence in acute coronary syndrome patients: a meta-analytic review PubMed Abstract | Google Scholar “Principles of psychosomatic assessment,” in The Psychosomatic Assessment Google Scholar Association of specific overt behaviour pattern with blood and cardiovascular findings: blood cholesterol level Google Scholar Updated method guidelines for systematic reviews in the Cochrane back review group Google Scholar Atrial fibrillation and Psychological factors: a systematic review International expert consensus document on takotsubo syndrome (Part I): clinical characteristics Comparing anxiety and depression in patients with takotsubo stress cardiomyopathy to those with acute coronary syndrome An overview of clinical research: the lay of the land PubMed Abstract | Google Scholar Symptoms of anxiety and depression and risk of acute myocardial infarction: the HUNT 2 study Higgins, J. P. T., and Green, S. (2011). Cochrane Handbook for Systematic Reviews of Interventions Version 5.1.0. The Cochrane Collaboration. Available at: www.handbook.cochrane.org (accessed March 2019) Google Scholar Hospitalization for depression is associated with an increased risk for myocardial infarction not explained by lifestyle Google Scholar International expert consensus document on Takotsubo syndrome (part II): diagnostic workup personality and psychiatric profiles in patientswith takotsubo cardiomyopathy Major depression and coronary artery disease in the Swedish twin registry: phenotypic Broken heart syndrome – is it a psychosomatic disorder Anxiety disorders and stressful events in Takotsubo syndrome Type A personality and mortality: competitiveness but not speed is associated with increased risk doi: 10.1016/j.atherosclerosis.2017.04.016 Current state of knowledge on takotsubo syndrome: a position statement from the taskforce on takotsubo syndrome of the heart failure association of the European society of cardiology Attempted suicide as a trigger of takotsubo syndrome: a mini review of available case reports (Letter) Google Scholar preferred reporting items for systematic reviews and meta-analyses: the PRISMA statement Google Scholar Pre-existing psychiatric illness is associated with increased risk of recurrent takotsubo cardiomyopathy Psychiatric Illness inTakotsubo (Stress) cardiomyopathy: a review Psychological distress and mortality in stable coronary heart disease: Persistence of high distress means increased risk Google Scholar Revised clinical diagnostic criteria for Tako-tsubo syndrome: the Tako-tsubo Italian Network proposal Google Scholar Human stress cardiomyopathy mimicking acute myocardial infarction Google Scholar Apical ballooning syndrome (Tako-Tsubo or stress cardiomyopathy): a mimic of acute myocardial infarction Google Scholar Psychological interventions for coronary heart disease: cochrane systematic review and meta-analysis Cumulative Impact of stressful life events on the developmentof takotsubo cardiomyopathy psychological distress and quality of life in female patients with Takotsubo cardiomyopathy: a prospective case-control study Salmoirago-Blotcher and personality characteristics among incident female cases of takotsubo cardiomyopathy: a case-control study Takotsubo syndrome is associated with mood disorders and antidepressants use not with anxiety and impairment of quality of life due to the psychiatric disorder Psychological distress and personality factors in takotsubo cardiomyopathy Clinical features and outcomes of Takotsubo (Stress) cardiomyopathy Child abuse and neglect and subclinical cardiovascular disease among midlife women Abuse and subclinical cardiovascular disease among midlife women: the study of women’s health across the nation A new insight into sudden cardiac death in young people The Newcastle-Ottawa scale (NOS) for Assessing the Quality of Non-Randomised Studies in Meta-Analyses Ottawa: Ottawa Hospital Research Institute Google Scholar Childhood psychological distress as a mediator in the relationship between early-life social disadvantage and adult cardiometabolic risk PubMed Abstract | Google Scholar Bursi F and Carugo S (2019) Traumatic Events Personality and Psychopathology in Takotsubo Syndrome: A Systematic Review Copyright © 2019 Galli, Bursi and Carugo. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY) *Correspondence: Federica Galli, ZmVkZXJpY2EuZ2FsbGkxQHVuaW1pLml0 Metrics details is a cell surface protein involved in homotypic cell–cell adhesion via intercellular oligomerization and proliferative signalling via proteolytic cleavage Despite its use as a diagnostic marker and being a drug target structural details of this conserved vertebrate-exclusive protein remain unknown Here we present the crystal structure of a heart-shaped dimer of the extracellular part of human EpCAM The structure represents a cis-dimer that would form at cell surfaces and may provide the necessary structural foundation for the proposed EpCAM intercellular trans-tetramerization mediated by a membrane-distal region we show how proteolytic processing at various sites could influence structural integrity oligomeric state and associated functionality of the molecule We also describe the epitopes of this therapeutically important protein and explain the antigenicity of its regions the outlined aspects of EpCAM biology are far from being well-defined especially from the structural point of view since to date no representative structure was available to interpret the experimental observations we present the structure of a non-glycosylated mutant form of EpEX (EpEXΔ) at 1.86-Å resolution a first and exemplary structure of the unique GA733 protein family EpEXΔ forms a cis-dimer corresponding to a half of the proposed trans-tetrameric intercellular unit We reveal the cis-dimerization interface between the two subunits involving the thyroglobulin type-1A (TY) domain and thereby describe a novel role of the versatile TY protein module EpCAM cis-dimer could be additionally stabilized by interaction of the TM helices of the two cis-oriented subunits as we explore by molecular dynamics (MD) simulations We show that the lateral protein surfaces of the wt EpCAM cis-dimer are partially covered with glycan chains as inferred from the position of known N-glycosylation sites The structure explains the high antigenicity of the surface-exposed small amino-terminal domain that harbours conformational epitopes of the majority of anti-EpCAM antibodies We provide a model of the trans-tetrameric intercellular unit where the trans-interactions involve the glycan-free membrane-distal platform-like surface of EpEXΔ we comment on the design of efficient EpCAM-targeted binder molecules and explain the detrimental effect of CTE-causing mutations on EpCAM structure and function It lacks the TM and intracellular part (EpIC) of full-length EpCAM it contains three N-glycosylation-abolishing mutations (indicated by black lines) (b) Cartoon representation of the complete EpEXΔ chain contained in the asymmetric unit shown in two different orientations Secondary structure elements are labelled as they appear along the polypeptide chain from N to C terminus (for example The three glutamine residues (N-glycosylation abolishing mutations) are shown as black sticks and all disulphide bridges are shown as sticks The N-terminal region tethered to the βA sheet is coloured in pale cyan Zoomed-in part shows contact between ND (green) and its copy from the adjacent asymmetric unit (grey) Detailed view (top) of TY-loop residues (dark purple) interacting with residues of βC sheet (grey side chains within 4 Å of TY loop shown in dark pink) The Arg80 and Arg81 (dark teal sticks) interact with Asp194 Gln74 and Gln111 (black sticks) introduced by glycosylation-abolishing mutations within TY correspond to glycosylated Asn residues in wt EpEX The three disulphide bonds of TY are shown as sticks the 6 and 28 kDa fragments (ND plus α1 and the rest of TY plus CD respectively) remain connected by the Cys66-Cys99 bond the C-terminal region of EpEXΔ (from Met261 to Lys265 plus His266) is not shown The top membrane-distal part of the molecule formed by RCD and ND is free of glycan chains the same fragment pattern could be generated using cathepsin K which also cleaves at the Gly79-Arg80 site such downward orientation would bring the C-terminal regions of the two subunits into close proximity it is plausible to assume that also the TM part (single helix per chain) forms a dimer with potential to additionally stabilize the EpCAM cis-dimer To investigate the dimerization propensity of the EpCAM TM helix we performed coarse-grained (CG) MD simulations of two modelled EpCAM TM helices embedded in a lipid bilayer The two EpCAM TM helices were allowed to diffuse freely within the bilayer patch from their initial positions separated by 50 Å (a) Interhelical distance calculated from six CG MD simulations started from the same initial state (b) Interhelical crossing angle distribution in combined simulation runs 1 and 3–6 (red) For combined runs only those parts where the two helices formed a dimer were used for calculation (c) Representative CG dimer model TM helices of EpCAM calculated by MD simulation Backbones of the two helices are shown as light and dark gray sticks isoleucine and valine residues are shown as blue Labels follow residue numbering in full-length EpCAM which are known modifiers of adhesion and signalling functions of TM molecules which is formed by the ND and the α2 helix plus RCD of the CD may represent the environment-sensing part of EpCAM while the lateral parts may be involved in interactions with other (trans)membrane proteins on the surface of the same cell the RCD and possibly the nearby α2 are proposed to be directly involved in EpCAM trans-tetramerization as described in our paper The role of all these cleavages and the exact mechanism by which they occur remains to be elucidated This latter mutation could impair formation of EpCAM cis-dimers or even prevent normal folding of EpCAM and adoption of the compact overall structure The identified epitope regions and the most antigenic N-domain are depicted as intensively shaded surface regions Listed below are other anti-EpEX antibodies for which the exact epitope is unknown Antibody names are color-coded by the domains they recognize (ND 17-1A and MT201 are targeted by Catumaxomab no adhesive function has been reported to date It is important that the similarities/differences are taken into account when developing diagnostic/prognostic and therapeutic approaches and great care is needed when analysing the results with regard to both EpCAM and Trop2 We prepared two EpEX variants—wt EpEX (residues 24-265 of EpCAM) and the non-glycosylated mutant variant (EpEXΔ) both with an uncleavable His6-tag at the C terminus the region of EpCAM cDNA (clone ID IRAKp961G0321Q; Source BioScience imaGenes) coding for Gln24-Lys265 was amplified by PCR In case of EpEXΔ three mutations were introduced: Asn74Gln 5′- and 3′-ends of the DNA fragment were extended by PCR to introduce the N-terminal melittin signal peptide (MKFLVNVALVFMVVYISYIYA) and a C-terminal His6-tag The resulting fragment was cloned into pFastBac1 which was in turn used to prepare recombinant baculoviruses using the Bac-to-Bac system (Invitrogen) Expression was performed using suspension cultures of Spodoptera frugiperda Sf9 cells at multiplicity of infection of 5 The recombinant protein was collected from the medium 36 h post infection by dialysis and subsequent immobilized metal affinity chromatography using Ni2+-charged HisTrap columns (GE Healthcare) The peak eluate fractions were pooled and EDTA (pH 8.0) was added to a final concentration of 5 mM Following dialysis and ion-exchange chromatography on Mono Q column (GE Healthcare) the final buffer exchange was done by size-exclusion chromatography on a Superdex 200 column (GE Healthcare) equilibrated in 20 mM Na-HEPES pH 8.0 The peak fraction was pooled and concentrated using Amicon centrifugal filter unit with 10 kDa cutoff membrane (Millipore) EpEX was prepared using the same procedure; however, the ion-exchange chromatography step was omitted. This recombinant EpEX variant is heterogeneously glycosylated as revealed by the presence of multiple bands on SDS–polyacrylamide gel electrophoresis (SDS–PAGE) gel, all at higher Mr than the EpEXΔ variant (Supplementary Fig. 1a) The interhelical distance was defined as the distance between centres of the two helices The interhelical crossing angle was defined as the angle between vectors describing the two helices; vector of each helix was drawn between centres of two groups of backbone atoms (Ala270-Val273 and Ala281-Val284) EpEX and EpEXΔ were incubated with human recombinant cathepsins L a HEPES buffer with pH 7.2 and an additional low-pH buffer were used (sodium acetate pH 5.5 or Bis-Tris pH 6.0) All buffers were 0.1 M and contained 2 mM EDTA and 2 mM DTT Cathepsin/EpEX(Δ) molar ratio was 10−2 or 10−3 Reaction mixtures were incubated at 37 °C for 1 h Reactions were stopped by adding reducing SDS–PAGE sample buffer followed by boiling Fragments were visualized on 15% SDS–PAGE by Coomassie staining Cleavage site was determined by N-terminal amino-acid sequencing (Edman degradation) of electrophoretically separated fragments transferred to a polyvinylidene difluoride membrane Crosslinking of intact or trypsin-cleaved EpEXΔ was performed in PBS by treating the protein with 1.3 mM BS3 (Pierce) at 21 °C for 30 min The reaction was quenched by adding Tris buffer (pH 7.4) to a final concentration of 50 mM followed by 30 min incubation at 21 °C Reaction mixtures were analysed on 15% SDS–PAGE under non-reducing and reducing conditions Protein bands were visualized by Coomassie staining The N-terminal pyroGlu residue was changed to Gln because modified residues are not generally accepted in docking simulations Histidine protonation states and flexible segments were defined automatically and centre of mass restraint was used to enforce contact between the molecules C2 point symmetry restraint was imposed between the two dimers (AB–CD) An additional calculation was performed where C2 point symmetry restrains between subunits in each subunit pair (A–C B–D) were used to limit calculation results to tetrameric units with perpendicular C2 rotation axes Docking solutions were first ranked according to their HADDOCK and Z-scores the latter indicating separation of a cluster of solutions from the average in terms of score (better solutions have more negative Z-scores) and further ranked according to the buried surface area and contacts between ND In the models from the best docking cluster (Z=−1.2) the relative angle between the two cis-dimers was 30° while in the models from second (Z=−0.9) and third best clusters (Z=0.1) the relative angle was 150° (ND’s already in contact) and 75° the best model from additional symmetry-restrained docking run was used Accession codes: Coordinates and structural factors have been deposited in the Protein Data Bank under accession code 4MZV Crystal structure and its bearing towards an understanding of key biological functions of EpCAM Efficient selection of human tumor growth-inhibiting monoclonal antibodies Frequent high-level expression of the immunotherapeutic target Ep-CAM in colon CD marker expression profiles of human embryonic stem cells and their neural derivatives reveal a novel CD marker for exclusion of pluripotent stem cells Characterization of epithelial cell adhesion molecule as a surface marker on undifferentiated human embryonic stem cells Epithelial cell adhesion molecule-targeted drug delivery for cancer therapy immunogenicity and bioactivity of the therapeutic antibody catumaxomab intraperitoneally administered to cancer patients Absence of cell-surface EpCAM in congenital tufting enteropathy Ep-CAM: a human epithelial antigen is a homophilic cell-cell adhesion molecule Evidence for a role of the epithelial glycoprotein 40 (Ep-CAM) in epithelial cell-cell adhesion The structural analysis of adhesions mediated by Ep-CAM Oligomeric state of the colon carcinoma-associated glycoprotein GA733-2 (Ep-CAM/EGP40) and its role in GA733-mediated homotypic cell-cell adhesion Epidermal growth factor-like repeats mediate lateral and reciprocal interactions of Ep-CAM molecules in homophilic adhesions Cytoplasmic tail regulates the intercellular adhesion function of the epithelial cell adhesion molecule Epithelial cell adhesion molecule (Ep-CAM) modulates cell-cell interactions mediated by classic cadherins Expression of Ep-CAM shifts the state of cadherin-mediated adhesions from strong to weak The cell-cell adhesion molecule EpCAM interacts directly with the tight junction protein claudin-7 EpCAM contributes to formation of functional tight junction in the intestinal epithelium by recruiting claudin proteins Epithelial cell adhesion molecule (EpCAM) regulates claudin dynamics and tight junctions Nuclear signalling by tumour-associated antigen EpCAM Initial activation of EpCAM cleavage via cell-to-cell contact Determination of disulphide bond assignments and N-glycosylation sites of the human gastrointestinal carcinoma antigen GA733-2 (CO17-1A Solution structure of mouse Cripto CFC domain and its inactive variant Trp107Ala a new tool for fast protein structure alignment in three dimensions WW domain-containing proteins: retrospectives and the future Characterization of the type-1 repeat from thyroglobulin a cysteine-rich module found in proteins from different families Structure of the RNA binding domain of a DEAD-box helicase bound to its ribosomal RNA target reveals a novel mode of recognition by an RNA recognition motif Structure of the mature ectodomain of the human receptor-type protein-tyrosine phosphatase IA-2 crystallization and preliminary X-ray characterization of the human epithelial cell-adhesion molecule ectodomain Inference of macromolecular assemblies from crystalline state Glycosylation is crucial for stability of tumour and cancer stem cell antigen EpCAM Biosynthesis and glycosylation of the carcinoma-associated antigen recognized by monoclonal antibody KS1/4 Biochemical analysis of a human epithelial surface antigen: differential cell expression and processing Isolation and characterization of a cDNA encoding the KS1/4 epithelial carcinoma marker Biochemical and immunological characterization of the human carcinoma-associated antigen MH 99/KS 1/4 Helix-helix interactions in membrane proteins: coarse-grained simulations of glycophorin a helix dimerization Major histocompatibility complex class II-associated p41 invariant chain fragment is a strong inhibitor of lysosomal cathepsin L Dual concentration-dependent activity of thyroglobulin type-1 domain of testican: specific inhibitor and substrate of cathepsin L Structural basis for the inhibition of insulin-like growth factors by insulin-like growth factor-binding proteins Diversity and evolution of the thyroglobulin type-1 domain superfamily and tetraspanins promotes colorectal cancer progression In vivo imaging of extracellular pH using 1H MRSI Cysteine cathepsins and the cutting edge of cancer invasion Tumor-specific glycosylation of the carcinoma-associated epithelial cell adhesion molecule EpCAM in head and neck carcinomas Crystal structure of ICAM-2 reveals a distinctive integrin recognition surface C-cadherin ectodomain structure and implications for cell adhesion mechanisms Urinary EpCAM in urothelial bladder cancer patients: characterisation and evaluation of biomarker potential Regulated intramembrane proteolysis and degradation of murine epithelial cell adhesion molecule mEpCAM Crystal structure of an active form of BACE1 an enzyme responsible for amyloid beta protein production EpCAM proteolysis: new fragments with distinct functions Identification of EpCAM as the gene for congenital tufting enteropathy Trop2 identifies a subpopulation of murine and human prostate basal cells with stem cell characteristics Trop2: a possible therapeutic target for late stage epithelial carcinomas Regulated proteolysis of Trop2 drives epithelial hyperplasia and stem cell self-renewal via β-catenin signaling PHENIX: a comprehensive Python-based system for macromolecular structure solution Substructure search procedures for macromolecular structures structure refinement and density modification with the PHENIX AutoBuild wizard Automated protein model building combined with iterative structure refinement ‘Conditional Restraints’: restraining the free atoms in ARP/wARP Refinement of macromolecular structures by the maximum-likelihood method Coot: model-building tools for molecular graphics MolProbity: all-atom structure validation for macromolecular crystallography GlyProt: in silico glycosylation of proteins Knowledge-based protein secondary structure assignment MultiSeq: unifying sequence and structure data for evolutionary analysis The MARTINI force field: coarse grained model for biomolecular simulations GROMACS 4: algorithms for highly efficient The HADDOCK web server for data-driven biomolecular docking Download references supported by the national GRID Initiatives of Belgium the Netherlands (via the Dutch BiG Grid project) Taiwan and the Latin America GRID infrastructure via the Gisela project is acknowledged for the use of web portals This work was supported by the Slovenian Research Agency (grant J1-2017) Faculty of Chemistry and Chemical Technology Kristina Djinović-Carugo & Brigita Lenarčič Department of Structural and Computational Biology analysed the results and wrote the manuscript Download citation Metrics details We compared the topologies of protein and small molecule crystals which have many common features – both are molecular crystals with intermolecular interactions much weaker than intramolecular interactions They also have different features – a considerably large fraction of the volume of protein crystals is occupied by liquid water while no room is available to other molecules in small molecule crystals We analyzed the overall and local topology and performed multilevel topological analyses (with the software package ToposPro) of carefully selected high quality sets of protein and small molecule crystal structures Given the suboptimal packing of protein crystals which is due the special shape and size of proteins it would be reasonable to expect that the topology of protein crystals is different from the topology of small molecule crystals we discovered that these two types of crystalline compounds have strikingly similar topologies This might suggest that molecular crystal formations share symmetry rules independent of molecular dimension This was explained by the model of deformable molecules in contrast to the model of rigid ones used by Kitaigorodskii We discarded structures containing nucleic acids and structures with an average B-factor smaller than 10 Å2 or larger than 40 Å2 We considered only structures with 50–500 amino acids in the asymmetric unit We discarded structures deposited without their experimental diffraction data Redundancy was reduced to 30% sequence identity Structures with missing protein atoms or with protein atoms deposited with zero occupancy were discarded as well as the structures where non-protein and non-water atoms are more than 5% of the total number of atoms resolution and R-factor thresholds ensure that unreliable structures are excluded from the analyzed data sets; moreover we considered a homogeneous set containing only low temperature crystal structures since packing might depend on temperature; furthermore limitations on the average B-factor and on protein dimension ensure that anomalous structures – much larger or smaller or much more or less flexible than customary proteins – are removed from the data sets; and eventually redundancy reduction to 30% sequence identity It is also important to observe that the exclusion of structures containing too many (more than 5%) hetero-atoms ensures that crystal contacts due to the presence of co-crystallized small molecules (the so called packing bridges) are minimized where some of the protein atoms/residues were undetected ensures that no protein-protein crystal packing contacts are neglected We prepared three ensembles of protein crystal structures we collected only monomeric proteins that crystallized with only one molecule per asymmetric unit (monomer set; 394 structures) we grouped dimeric proteins that crystallized with only one dimer per asymmetric unit (dimer set; 207 structures) we pulled together structures of monomeric proteins with two molecules in the asymmetric unit (double set; 164 structures) The identification codes of all these crystal structures are listed in the Supplementary Information (Table S1) Crystallographic data for organic molecular crystals consisting of chemically equivalent molecules were taken from the CSD (release 5.38) We have not analyzed molecular crystals of coordination or organometallic compounds which have different chemical nature compared to proteins which contain one independent molecule and 10,240 structures with two independent chemically equivalent molecules in the asymmetric unit or disordered structures as well as those with Rf > 10% were excluded; no other restrictions were applied The crystal contacts were identified as described previously13 The asymmetric unit was transformed according to all symmetry operations and translated up to three times along the axes a and c in both negative and positive directions is symmetry related to the first – were considered to be in contact if at least one atom of one of them is closer than 4.5 Å from an atom of the other one We refer to this method below as ‘Distance’ method The ‘Domains’ method is designed to process large samples of crystal structures in an automated mode this method was never used for proteins; thus it was important to check if this method gives similar results for proteins as the ‘Distance’ method combined from the faces of their Voronoi polyhedra which correspond to the intermolecular contacts This method was also used to compute coordination numbers of some protein molecules for comparison with the simple geometrical approach described above its center of mass and centers of mass of four neighboring molecules (pink balls) which represent a fragment of the underlying net where N is the series of coordination numbers for all independent nodes indicating the dimensionality of the net (C respectively) and n is the number of the net with a given ND sequence in the TTD collection the 16T4 symbol means a three-periodic net with 16-coordinated nodes (i.e any molecule has 16 neighbors in the packing) and it is the fourth 16-coordinated net in the TTD collection Since the strength of intermolecular contacts is quite different the topology of a molecular packing strongly depends on which contacts are taken into account Ignoring interactions of a particular level of strength we generate different underlying nets for the same crystal structure Each underlying net specifies the packing topology at a given level of interaction and corresponds to a structure representation hcp topology; (b) 10 coordinated underlying net bct topology; (c) 8 coordinated underlying net hex topology; (d) 6 coordinated underlying net hxl topology; (e) 4 coordinated underlying net connected at different levels of interactions are in yellow strengths of crystal packing contacts were estimated based on their dissociation free energies as described below which were generated by applying all the symmetry operations to the asymmetric unit and by applying three translations along all the axes in the negative and in the positive directions Intermolecular contacts were defined as pairs of atoms one in the reference molecule and the other in a satellite molecule which are determined by the ‘Domains’ method Distances between the molecular centers of mass are given in Å Thus protein crystals follow the close packing topologies but with some gaps which are obviously filled by the solvent molecules and decrease coordination numbers of protein molecules It is interesting to observe that while few topologies are extremely more frequent than others in small molecule crystals 30% of the small molecule crystals are associated with the most frequent topology (bcu-x) and 65% of them have one of the four most common topologies (bcu-x only 6% of the protein crystals are associated with the most frequent topology (fcu) and only 19% of them have one of the four most common topologies (fcu This suggests that the suboptimal packing of the protein molecules in the crystal state allows a wider number of topologies and none of them can be much more frequent than the others Distribution of the coordination numbers (%) for all the structures of proteins examined (see full data in Table S2) compared to the CN for small molecules we see that the most frequent coordination number of small molecules (14) is observed in a large fraction of crystals (52%) while the most frequent coordination number for proteins (7 or 8) is observed only in a smaller fraction of crystals (about 20%) small molecules can adopt 25 different coordination number values (from 4 to 32) while proteins can have only 14 different coordination number values (ranging from 1 to 16 – coordination numbers 2 and 15 are never observed) All these observations point out that proteins which are suboptimally packed in their crystals can hardly by surrounded by numerous other proteins and that the palette of their coordination numbers is considerably more limited This means that from 5 to 26% of protein molecules essentially differ from small molecules by their ability to be packed In order to reach a better understanding of the differences between the topologies of protein and small molecule crystals we performed the multilevel analysis for the molecular packings in 394 monomeric protein structures and 105,549 structures of small molecules 1,799 and 1,144,539 structure representations were generated Distribution of the underlying nets observed in the multilevel topological analysis of the proteins and small molecule crystals ranked by CN and scaled to 100 on the most frequent in each set (dia for proteins and layers sql for small molecules) Most typical underlying motifs in packings of protein and small organic molecules protein crystals nucleation occurs at very high levels of supersaturation often two or three orders of magnitude greater than that required to sustain crystal growth proteins may assume several distinctive solid states that include amorphous precipitates Further studies are necessary and additional data must be considered to find out the rationale why topology seems to be independent of packing efficiency and crystallization Die Chemismus in der thierischen Organisation Criteria to extract high quality Protein Data Bank subsets for structure users Download references Krissinel for his advices on the use of PISA and on the interpretation of the results thanks Kristina Djinović for her kind hospitality in Vienna are grateful to the Russian Government (Grant 14.B25.31.0005) for support Samara Center for Theoretical Materials Science (SCTMS) School of Materials Science and Engineering The authors declare that they have no competing interests Download citation DOI: https://doi.org/10.1038/s41598-017-12699-4 Sorry, a shareable link is not currently available for this article. Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily. What would a Kim Dolce & Gabbana show be like prepping for the Kim Dolce & Gabbana show. We're doing fittings and everything's kind of crazy I've been spending a lot of time over the years I have a really big Dolce & Gabbana archive. So when they saw a lot of the looks that I had, And we just started having these conversations with this like retrospective of the last 25 years And then when Kourtney got married at their home, it kind of reignited a whole new love and friendship. We all have such different style and personalities. So to see everyone wear Dolce & Gabbana There's something just extra sensual about Italian style by the '90s supermodels and early 2000s. And to get every single one of those looks reimagined I didn't really understand that each look to express what we're trying to express. that there's a science behind it all. I have, I think every single piece in the show. Just on a regular day in Milan wearing this? This whole process has been so eye opening I wore these bottoms to Kourtney and Travis's wedding. you cannot beat Dolce & Gabbana early '90s, Some of my favorites is this little bolero And we had it remade with all of the crystal, which I feel like was a little bit more me. And then of course we have all of the original corsets. it came in, I think three different colors of crystals. Who doesn't remember Gisele wearing this bra. and all I wanted to do was save up for this little bra. I think crystals are definitely making a comeback Seeing their relationship together is so inspiring. and that's always what I've been about. I signed your praises earlier so I was like... So thank you Vogue for coming on this journey with us. And I can't wait for you to see the show. Kristina Djinovic-Carugo (bottom left) and her team at the Center for Optimized Structural Studies (COSS) who have successfully extended their funding as a "Laura Bassi Centre of Excellence" boosting structural research in Vienna with 1 million euro until 2016 colored by its electrostatic potential (red representing negative charges and blue representing positive charges) Utl.: "Laura Bassi Centre of Excellence" of Kristina Djinovic-Carugo extended for another three years