Polypharmacology-Driven Discovery and Design of Highly Selective, Dual and Multitargeting Inhibitors of Mycobacterium tuberculosis - A Review.

The increasing demand for novel antitubercular agents has been the main 'force' of many TB research efforts due to the uncontrolled growing number of drug-resistant strains of M. tuberculosis in the clinical setting. Many strategies have been employed to address the drug-resistant issue, including a trend that is gaining attention, which is the design and discovery of Mtb inhibitors that are either dual- or multitargeting. The multiple-target design concept is not new in medicinal chemistry. With a growing number of newly discovered Mtb proteins, numerous targets are now available for developing new biochemical/cell-based assays and computer-aided drug design (CADD) protocols. To describe the achievements and overarching picture of this field in anti- infective drug discovery, we provide in this review small molecules that exhibit profound inhibitory activity against the tubercle bacilli and are identified to trace two or more Mtb targets. This review also presents emerging design methodologies for developing new anti-TB agents, particularly tailored to structure-based CADD. Dedicated on the special occasion of the 70th birthday of Prof. Dr. Ma. Alicia Aguinaldo, whose scientific efforts elevated antituberculosis drug discovery in the Philippines.

How to repair medial meniscal ramp lesions: A systematic review of surgical techniques.

Purpose: to provide a comprehensive overview of all the surgical techniques published in the literature for repairing meniscal ramp lesions focusing on the technical aspects and the pros and cons of every procedure. Such lesions can be managed using various approaches, each of this with its specific advantages and disadvantages. Methods: Pubmed Central, Scopus, and EMBASE databases were systematically reviewed according to the preferred reporting items for systematic reviews and meta-analysis (PRISMA) guidelines for studies on surgical techniques for repairing meniscal ramp lesions through May 2023. Overall, 32 articles matched the selection criteria and were included in the study. Results: Debridement alone may be sufficient for small stable meniscal ramp lesions but, for tears in the menisco-capsular junction that affect the stability of the medial meniscus, it seems reasonable to repair it, even though the clinical results available in literature are contrasting. All-inside sutures through anterior portals seems to be an effective solution for meniscal ramp lesions with MTL tears. All-inside sutures through posteromedial portals are particularly useful for large meniscal ramp lesions, in which an inside-out suture can also be performed. Conclusion: Meniscal ramp lesions can be managed using various approaches, each of this with its specific advantages and disadvantages. Further research is required to determine the optimal technique that can be considered as the gold standard and can provide the better results. Level of Evidence: Level III, systematic review.

Social dysfunction relates to shifts within socioaffective brain systems among Schizophrenia and Alzheimer's disease patients.

Social dysfunction represents one of the most common signs of neuropsychiatric disorders, such as Schizophrenia (SZ) and Alzheimer's disease (AD). Perturbed socioaffective neural processing is crucially implicated in SZ/AD and generally linked to social dysfunction. Yet, transdiagnostic properties of social dysfunction and its neurobiological underpinnings remain unknown. As part of the European PRISM project, we examined whether social dysfunction maps onto shifts within socioaffective brain systems across SZ and AD patients. We probed coupling of social dysfunction with socioaffective neural processing, as indexed by an implicit facial emotional processing fMRI task, across SZ (N = 46), AD (N = 40) and two age-matched healthy control (HC) groups (N = 26 HC-younger and N = 27 HC-older). Behavioural (i.e., social withdrawal, interpersonal dysfunction, diminished prosocial or recreational activity) and subjective (i.e., feelings of loneliness) aspects of social dysfunction were assessed using the Social Functioning Scale and De Jong-Gierveld loneliness questionnaire, respectively. Across SZ/AD/HC participants, more severe behavioural social dysfunction related to hyperactivity within fronto-parieto-limbic brain systems in response to sad emotions (P = 0.0078), along with hypoactivity of these brain systems in response to happy emotions (P = 0.0418). Such relationships were not found for subjective experiences of social dysfunction. These effects were independent of diagnosis, and not confounded by clinical and sociodemographic factors. In conclusion, behavioural aspects of social dysfunction across SZ/AD/HC participants are associated with shifts within fronto-parieto-limbic brain systems. These findings pinpoint altered socioaffective neural processing as a putative marker for social dysfunction, and could aid personalized care initiatives grounded in social behaviour.

Surrogate modeling and control of medical digital twins.

The vision of personalized medicine is to identify interventions that maintain or restore a person's health based on their individual biology. Medical digital twins, computational models that integrate a wide range of health-related data about a person and can be dynamically updated, are a key technology that can help guide medical decisions. Such medical digital twin models can be high-dimensional, multi-scale, and stochastic. To be practical for healthcare applications, they often need to be simplified into low-dimensional surrogate models that can be used for optimal design of interventions. This paper introduces surrogate modeling algorithms for the purpose of optimal control applications. As a use case, we focus on agent-based models (ABMs), a common model type in biomedicine for which there are no readily available optimal control algorithms. By deriving surrogate models that are based on systems of ordinary differential equations, we show how optimal control methods can be employed to compute effective interventions, which can then be lifted back to a given ABM. The relevance of the methods introduced here extends beyond medical digital twins to other complex dynamical systems.

Control of Medical Digital Twins with Artificial Neural Networks.

The objective of personalized medicine is to tailor interventions to an individual patient's unique characteristics. A key technology for this purpose involves medical digital twins, computational models of human biology that can be personalized and dynamically updated to incorporate patient-specific data collected over time. Certain aspects of human biology, such as the immune system, are not easily captured with physics-based models, such as differential equations. Instead, they are often multi-scale, stochastic, and hybrid. This poses a challenge to existing model-based control and optimization approaches that cannot be readily applied to such models. Recent advances in automatic differentiation and neural-network control methods hold promise in addressing complex control problems. However, the application of these approaches to biomedical systems is still in its early stages. This work introduces dynamics-informed neural-network controllers as an alternative approach to control of medical digital twins. As a first use case for this method, the focus is on agent-based models, a versatile and increasingly common modeling platform in biomedicine. The effectiveness of the proposed neural-network control method is illustrated and benchmarked against other methods with two widely-used agent-based model types. The relevance of the method introduced here extends beyond medical digital twins to other complex dynamical systems.

Studying Venom Toxin Variation Using Accurate Masses from Liquid Chromatography-Mass Spectrometry Coupled with Bioinformatic Tools.

This study provides a new methodology for the rapid analysis of numerous venom samples in an automated fashion. Here, we use LC-MS (Liquid Chromatography-Mass Spectrometry) for venom separation and toxin analysis at the accurate mass level combined with new in-house written bioinformatic scripts to obtain high-throughput results. This analytical methodology was validated using 31 venoms from all members of a monophyletic clade of Australian elapids: brown snakes (Pseudonaja spp.) and taipans (Oxyuranus spp.). In a previous study, we revealed extensive venom variation within this clade, but the data was manually processed and MS peaks were integrated into a time-consuming and labour-intensive approach. By comparing the manual approach to our new automated approach, we now present a faster and more efficient pipeline for analysing venom variation. Pooled venom separations with post-column toxin fractionations were performed for subsequent high-throughput venomics to obtain toxin IDs correlating to accurate masses for all fractionated toxins. This workflow adds another dimension to the field of venom analysis by providing opportunities to rapidly perform in-depth studies on venom variation. Our pipeline opens new possibilities for studying animal venoms as evolutionary model systems and investigating venom variation to aid in the development of better antivenoms.

Case report: A familial B-acute lymphoblastic leukemia associated with a new germline pathogenic variant in PAX5. The first report in Mexico.

B-cell acute lymphoblastic leukemia (B-ALL) is one of the most common childhood cancers worldwide. Although most cases are sporadic, some familial forms, inherited as autosomal dominant traits with incomplete penetrance, have been described over the last few years. Germline pathogenic variants in transcription factors such as PAX5, IKZF1, and ETV6 have been identified as causal in familial forms. The proband was a 7-year-old Mexican girl diagnosed with high-risk B-ALL at five years and 11 months of age. Family history showed that the proband's mother had high-risk B-ALL at 16 months of age. She received chemotherapy and was discharged at nine years of age without any evidence of recurrence of leukemia. The proband's father was outside the family nucleus, but no history of leukemia or cancer was present up to the last contact with the mother. We performed exome sequencing on the proband and the proband's mother and identified the PAX5 variant NM_016734.3:c.963del: p.(Ala322LeufsTer11), located in the transactivation domain of the PAX5 protein. The variant was classified as probably pathogenic according to the ACMG criteria. To the best of our knowledge, this is the first Mexican family with an inherited increased risk of childhood B-ALL caused by a novel germline pathogenic variant of PAX5. Identifying individuals with a hereditary predisposition to cancer is essential for modern oncological practice. Individuals at high risk of leukemia would benefit from hematopoietic stem cell transplantation, but family members carrying the pathogenic variant should be excluded as hematopoietic stem cell donors.

Control of Medical Digital Twins with Artificial Neural Networks.

The objective of personalized medicine is to tailor interventions to an individual patient's unique characteristics. A key technology for this purpose involves medical digital twins, computational models of human biology that can be personalized and dynamically updated to incorporate patient-specific data collected over time. Certain aspects of human biology, such as the immune system, are not easily captured with physics-based models, such as differential equations. Instead, they are often multi-scale, stochastic, and hybrid. This poses a challenge to existing model-based control and optimization approaches that cannot be readily applied to such models. Recent advances in automatic differentiation and neural-network control methods hold promise in addressing complex control problems. However, the application of these approaches to biomedical systems is still in its early stages. This work introduces dynamics-informed neural-network controllers as an alternative approach to control of medical digital twins. As a first use case for this method, the focus is on agent-based models, a versatile and increasingly common modeling platform in biomedicine. The effectiveness of the proposed neural-network control method is illustrated and benchmarked against other methods with two widely-used agent-based model types. The relevance of the method introduced here extends beyond medical digital twins to other complex dynamical systems.

Toward mechanistic medical digital twins: some use cases in immunology.

A fundamental challenge for personalized medicine is to capture enough of the complexity of an individual patient to determine an optimal way to keep them healthy or restore their health. This will require personalized computational models of sufficient resolution and with enough mechanistic information to provide actionable information to the clinician. Such personalized models are increasingly referred to as medical digital twins. Digital twin technology for health applications is still in its infancy, and extensive research and development is required. This article focuses on several projects in different stages of development that can lead to specific-and practical-medical digital twins or digital twin modeling platforms. It emerged from a two-day forum on problems related to medical digital twins, particularly those involving an immune system component. Open access video recordings of the forum discussions are available.

Polymer Physics: From Theory to Experimental Applications.

The significance of polymer processing techniques cannot be overstated in the production of polymer components [...].

Temporal Stability of Management Zone Patterns: Case Study with Contact and Non-Contact Soil Electrical Conductivity Sensors in Dryland Pastures.

Precision agriculture (PA) intends to validate technological tools that capture soil and crop spatial variability, which constitute the basis for the establishment of differentiated management zones (MZs). Soil apparent electrical conductivity (ECa) sensors are commonly used to survey soil spatial variability. It is essential for surveys to have temporal stability to ensure correct medium- and long-term decisions. The aim of this study was to assess the temporal stability of MZ patterns using different types of ECa sensors, namely an ECa contact-type sensor (Veris 2000 XA, Veris Technologies, Salina, KS, USA) and an electromagnetic induction sensor (EM-38, Geonics Ltd., Mississauga, ON, Canada). These sensors were used in four fields of dryland pastures in the Alentejo region of Portugal. The first survey was carried out in October 2018, and the second was carried out in September 2020. Data processing involved synchronizing the geographic coordinates obtained using the two types of sensors in each location and establishing MZs based on a geostatistical analysis of elevation and ECa data. Although the basic technologies have different principles (contact versus non-contact sensors), the surveys were carried out at different soil moisture conditions and were temporarily separated (about 2 years); the ECa measurements showed statistically significant correlations in all experimental fields (correlation coefficients between 0.449 and 0.618), which were reflected in the spatially stable patterns of the MZ maps (averaging 52% of the total area across the four experimental fields). These results provide perspectives for future developments, which will need to occur in the creation of algorithms that allow the spatial variability and temporal stability of ECa to be validated through smart soil sampling and analysis to generate recommendations for sustained soil amendment or fertilization.

AP-1γ2 is an adaptor protein 1 variant required for endosome-to-Golgi trafficking of the mannose-6-P receptor (CI-MPR) and ATP7B copper transporter.

Selective retrograde transport from endosomes back to the trans-Golgi network (TGN) is important for maintaining protein homeostasis, recycling receptors, and returning molecules that were transported to the wrong compartments. Two important transmembrane proteins directed to this pathway are the Cation-Independent Mannose-6-phosphate receptor (CI-MPR) and the ATP7B copper transporter. Among CI-MPR functions is the delivery of acid hydrolases to lysosomes, while ATP7B facilitates the transport of cytosolic copper ions into organelles or the extracellular space. Precise subcellular localization of CI-MPR and ATP7B is essential for the proper functioning of these proteins. This study shows that both CI-MPR and ATP7B interact with a variant of the clathrin adaptor 1 (AP-1) complex that contains a specific isoform of the γ-adaptin subunit called γ2. Through synchronized anterograde trafficking and cell-surface uptake assays, we demonstrated that AP-1γ2 is dispensable for ATP7B and CI-MPR exit from the TGN while being critically required for ATP7B and CI-MPR retrieval from endosomes to the TGN. Moreover, AP-1γ2 depletion leads to the retention of endocytosed CI-MPR in endosomes enriched in retromer complex subunits. These data underscore the importance of AP-1γ2 as a key component in the sorting and trafficking machinery of CI-MPR and ATP7B, highlighting its essential role in the transport of proteins from endosomes.

Forum on immune digital twins: a meeting report.

Medical digital twins are computational models of human biology relevant to a given medical condition, which are tailored to an individual patient, thereby predicting the course of disease and individualized treatments, an important goal of personalized medicine. The immune system, which has a central role in many diseases, is highly heterogeneous between individuals, and thus poses a major challenge for this technology. In February 2023, an international group of experts convened for two days to discuss these challenges related to immune digital twins. The group consisted of clinicians, immunologists, biologists, and mathematical modelers, representative of the interdisciplinary nature of medical digital twin development. A video recording of the entire event is available. This paper presents a synopsis of the discussions, brief descriptions of ongoing digital twin projects at different stages of progress. It also proposes a 5-year action plan for further developing this technology. The main recommendations are to identify and pursue a small number of promising use cases, to develop stimulation-specific assays of immune function in a clinical setting, and to develop a database of existing computational immune models, as well as advanced modeling technology and infrastructure.

Association Between Years of Education and Amyloid Burden in Patients With Subjective Cognitive Decline, MCI, and Alzheimer Disease.

OBJECTIVES: Higher-educated patients with Alzheimer disease (AD) can harbor greater neuropathologic burden than those with less education despite similar symptom severity. In this study, we assessed whether this observation is also present in potential preclinical AD stages, namely in individuals with subjective cognitive decline and clinical features increasing AD likelihood (SCD+). METHODS: Amyloid-PET information ([18F]Flutemetamol or [18F]Florbetaben) of individuals with SCD+, mild cognitive impairment (MCI), and AD were retrieved from the AMYPAD-DPMS cohort, a multicenter randomized controlled study. Group classification was based on the recommendations by the SCD-I and NIA-AA working groups. Amyloid PET images were acquired within 8 months after initial screening and processed with AMYPYPE. Amyloid load was based on global Centiloid (CL) values. Educational level was indexed by formal schooling and subsequent higher education in years. Using linear regression analysis, the main effect of education on CL values was tested across the entire cohort, followed by the assessment of an education-by-diagnostic-group interaction (covariates: age, sex, and recruiting memory clinic). To account for influences of non-AD pathology and comorbidities concerning the tested amyloid-education association, we compared white matter hyperintensity (WMH) severity, cardiovascular events, depression, and anxiety history between lower-educated and higher-educated groups within each diagnostic category using the Fisher exact test or χ2 test. Education groups were defined using a median split on education (Md = 13 years) in a subsample of the initial cohort, for whom this information was available. RESULTS: Across the cohort of 212 individuals with SCD+ (M(Age) = 69.17 years, F 42.45%), 258 individuals with MCI (M(Age) = 72.93, F 43.80%), and 195 individuals with dementia (M(Age) = 74.07, F 48.72%), no main effect of education (ß = 0.52, 95% CI -0.30 to 1.58), but a significant education-by-group interaction on CL values, was found (p = 0.024) using linear regression modeling. This interaction was driven by a negative association of education and CL values in the SCD+ group (ß = -0.11, 95% CI -4.85 to -0.21) and a positive association in the MCI group (ß = 0.15, 95% CI 0.79-5.22). No education-dependent differences in terms of WMH severity and comorbidities were found in the subsample (100 cases with SCD+, 97 cases with MCI, 72 cases with dementia). DISCUSSION: Education may represent a factor oppositely modulating subjective awareness in preclinical stages and objective severity of ongoing neuropathologic processes in clinical stages.

Community composition of black flies during and after the 2020 vesicular stomatitis virus outbreak in Southern New Mexico, USA.

BACKGROUND: Vesicular stomatitis virus (VSV), a vector-borne pathogen of livestock, emerges periodically in the western US. In New Mexico (NM), US, most cases occur close to the Rio Grande River, implicating black flies (Simulium spp.) as a possible vector. In 2020, VS cases were reported in NM from April to May, although total black fly abundance remained high until September. We investigated the hypothesis that transience of local VSV transmission results from transient abundance of key, competent black fly species. Additionally, we investigated whether irrigation canals in southern NM support a different community of black flies than the main river. Lastly, to gain insight into the source of local black flies, in 2023 we collected black fly larvae prior to the release of water into the Rio Grande River channel. METHODS: We randomly sub-sampled adult black flies collected along the Rio Grande during and after the 2020 VSV outbreak. We also collected black fly adults along the river in 2021 and 2022 and at southern NM farms and irrigation canals in 2022. Black fly larvae were collected from dams in the area in 2023. All collections were counted, and individual specimens were subjected to molecular barcoding for species identification. RESULTS: DNA barcoding of adult black flies detected four species in 2020: Simulium meridionale (N = 158), S. mediovittatum (N = 83), S. robynae (N = 26) and S. griseum/notatum (N = 1). Simulium robynae was only detected during the VSV outbreak period, S. meridionale showed higher relative abundance, but lower absolute abundance, during the outbreak than post-outbreak period, and S. mediovittatum was rare during the outbreak period but predominated later in the summer. In 2022, relative abundance of black fly species did not differ significantly between the Rio Grande sites and farm and irrigation canals. Intriguingly, 63 larval black flies comprised 56% Simulium vittatum, 43% S. argus and 1% S. encisoi species that were either extremely rare or not detected in previous adult collections. CONCLUSIONS: Our results suggest that S. robynae and S. meridionale could be shaping patterns of VSV transmission in southern NM. Thus, field studies of the source of these species as well as vector competence studies are warranted.

The Impact Of Thymectomy In Thymomatous And Nonthymomatous Myasthenia Gravis - The Experience Of A Tertiary Center.

INTRODUCTION: Thymectomy remains a mainstay of treatment in Thymomatous (T) and Nonthymomatous (nT) Myasthenia Gravis (MG), with improved clinical outcomes and reduced need for medical treatment, however, there is little research regarding long-term follow-up. We aim to assess the impact of surgery on the long-term outcome of patients with MG at our center. METHODS: Retrospective analyses of MG patients submitted to thymectomy between 2007 and 2017 at the thoracic surgery department of CHUC. Clinical assessment was performed according to the MG Foundation of America (MGFA) Clinical Classification (cMGFA). The follow-up was categorized according to the MGFA Post-intervention Status (MGFA-PIS) and cMGFA. Statistical analysis was performed with SPSS, to a significance level of 5%. RESULTS: Thirty-seven patients underwent extended thymectomy and 67.6% were female. Median age at diagnosis was 46.68±19.2 years. Most patients (83.8%) had anti-acetylcholine receptor antibodies and 81.1% had generalized forms of MG. Many patients (67.6%) had surgery less than 12 months after the clinical diagnosis. TMG was present in 19 (51.4%) patients. Compared to nTMG, these patients were older (54.06±17.9 vs 40.17±19.4 years) and most were men (52.9% vs 16.7%). We obtained a good outcome in most patients in the first (81.1%), second (86.1%), and fifth (84.8%) year of follow-up. There was a shift towards better prognosis categories in the good outcome group: 9.1% CSR, 3.0% PR, and 66,7% MM in the fifth year. Preoperative medical treatment did not influence the long-term follow-up outcome. A shorter time to surgery (< 12 months) correlated with better outcomes at year 5 (p=0.016). CONCLUSION: Thymectomy led to a sustained clinical improvement in our cohort, allowing for a reduced need for medication. A shorter time to surgery seems to have a positive influence on long-term prognosis. We expect that an extended follow-up would improve our results.

RHO-3 plays a significant role in hyphal extension rate, conidiation, and the integrity of the Spitzenkörper in Neurospora crassa.

The Rho family of monomeric GTPases act as signaling proteins to establish and maintain cell polarity and other essential cellular processes. Rho3 is a GTPase of the Rho family that is exclusive of fungi that regulate cell polarity in yeast. However, studies have yet to explore its function in filamentous fungi. In this work, we investigated the role of RHO-3 in the model organism Neurospora crassa. Confocal microscopy analysis revealed that RHO-3 localizes in the outer region of the Spitzenkörper (Spk), in the plasma membrane from region II to the beginning of region III, and in the septa of mature hyphae. The phenotypic effect of the rho-3 deletion was analyzed. The results revealed that the rho-3 null strain showed severe defects in growth rate, aerial hyphae length, and conidia production. The organization of the Spk is also affected in the absence of RHO-3. Co-expression analysis of GFP-RHO-3 with glucan synthase 1 (GS-1-mChFP) and chitin synthase 1 (CHS-1-mChFP) revealed that RHO-3 localizes in the external region of the Spk in the macrovesicles zone. In summary, our results suggest that RHO-3 is not essential for the polarized growth of hyphae but plays a significant role in hyphal extension rate, conidiation, sexual reproduction and the integrity of the Spk, possibly regulating the delivery of macrovesicles to the apical dome.

Scale-Free Chaos in the 2D Harmonically Confined Vicsek Model.

Animal motion and flocking are ubiquitous nonequilibrium phenomena that are often studied within active matter. In examples such as insect swarms, macroscopic quantities exhibit power laws with measurable critical exponents and ideas from phase transitions and statistical mechanics have been explored to explain them. The widely used Vicsek model with periodic boundary conditions has an ordering phase transition but the corresponding homogeneous ordered or disordered phases are different from observations of natural swarms. If a harmonic potential (instead of a periodic box) is used to confine particles, then the numerical simulations of the Vicsek model display periodic, quasiperiodic, and chaotic attractors. The latter are scale-free on critical curves that produce power laws and critical exponents. Here, we investigate the scale-free chaos phase transition in two space dimensions. We show that the shape of the chaotic swarm on the critical curve reflects the split between the core and the vapor of insects observed in midge swarms and that the dynamic correlation function collapses only for a finite interval of small scaled times. We explain the algorithms used to calculate the largest Lyapunov exponents, the static and dynamic critical exponents, and compare them to those of the three-dimensional model.

HIV-1 Nef Changes the Proteome of T Cells Extracellular Vesicles Depleting IFITMs and Other Antiviral Factors.

Extracellular vesicles (EVs) are biomolecule carriers for intercellular communication in health and disease. Nef is a HIV virulence factor that is released from cells within EVs and is present in plasma EVs of HIV-1 infected individuals. We performed a quantitative proteomic analysis to fully characterize the Nef-induced changes in protein composition of T cell-derived EVs and identify novel host targets of HIV. Several proteins with well-described roles in infection or not previously associated with HIV pathogenesis were specifically modulated by Nef in EVs. Among the downregulated proteins are the interferon-induced transmembrane 1, 2, and 3 (IFITM1-3) proteins, broad-spectrum antiviral factors known to be cell-to-cell transferable by EVs. We demonstrate that Nef depletes IFITM1-3 from EVs by excluding these proteins from the plasma membrane and lipid rafts, which are sites of EVs biogenesis in T cells. Our data establish Nef as a modulator of EVs' global protein content and as an HIV factor that antagonizes IFITMs.

Coarse-graining the vertex model and its response to shear.

Tissue dynamics and collective cell motion are crucial biological processes. Their biological machinery is mostly known, and simulation models such as the active vertex model exist and yield reasonable agreement with experimental observations such as tissue fluidization or fingering. However, a good and well-founded continuum description for tissues remains to be developed. In this work, we derive a macroscopic description for a two-dimensional cell monolayer by coarse-graining the vertex model through the Poisson bracket approach. We obtain equations for cell density, velocity, and the cellular shape tensor. We then study the homogeneous steady states, their stability (which coincides with thermodynamic stability), and especially their behavior under an externally applied shear. Our results contribute to elucidate the interplay between flow and cellular shape. The obtained macroscopic equations present a good starting point for adding cell motion, morphogenetic, and other biologically relevant processes.