1
|
Biehn SE, Goncalves LM, Lehmann J, Marty JD, Mueller C, Ramirez SA, Tillier F, Sage CR. BioPrint meets the AI age: development of artificial intelligence-based ADMET models for the drug-discovery platform SAFIRE. Future Med Chem 2024; 16:587-599. [PMID: 38372202 DOI: 10.4155/fmc-2024-0007] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/20/2024] Open
Abstract
Background: To prioritize compounds with a higher likelihood of success, artificial intelligence models can be used to predict absorption, distribution, metabolism, excretion and toxicity (ADMET) properties of molecules quickly and efficiently. Methods: Models were trained with BioPrint database proprietary data along with public datasets to predict various ADMET end points for the SAFIRE platform. Results: SAFIRE models performed at or above 75% accuracy and 0.4 Matthew's correlation coefficient with validation sets. Training with both proprietary and public data improved model performance and expanded the chemical space on which the models were trained. The platform features scoring functionality to guide user decision-making. Conclusion: High-quality datasets along with chemical space considerations yielded ADMET models performing favorably with utility in the drug discovery process.
Collapse
Affiliation(s)
- Sarah E Biehn
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | | | - Juerg Lehmann
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | - Jessica D Marty
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | - Christoph Mueller
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | - Samuel A Ramirez
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | - Fabien Tillier
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| | - Carleton R Sage
- Eurofins DiscoveryAI, Eurofins Panlabs, Inc., Saint Charles, MO 63304, USA
| |
Collapse
|
2
|
Ramirez SA, Pablo M, Burk S, Lew DJ, Elston TC. A novel stochastic simulation approach enables exploration of mechanisms for regulating polarity site movement. PLoS Comput Biol 2021; 17:e1008525. [PMID: 34264926 PMCID: PMC8315557 DOI: 10.1371/journal.pcbi.1008525] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2020] [Revised: 07/27/2021] [Accepted: 06/24/2021] [Indexed: 12/23/2022] Open
Abstract
Cells polarize their movement or growth toward external directional cues in many different contexts. For example, budding yeast cells grow toward potential mating partners in response to pheromone gradients. Directed growth is controlled by polarity factors that assemble into clusters at the cell membrane. The clusters assemble, disassemble, and move between different regions of the membrane before eventually forming a stable polarity site directed toward the pheromone source. Pathways that regulate clustering have been identified but the molecular mechanisms that regulate cluster mobility are not well understood. To gain insight into the contribution of chemical noise to cluster behavior we simulated clustering using the reaction-diffusion master equation (RDME) framework to account for molecular-level fluctuations. RDME simulations are a computationally efficient approximation, but their results can diverge from the underlying microscopic dynamics. We implemented novel concentration-dependent rate constants that improved the accuracy of RDME-based simulations, allowing us to efficiently investigate how cluster dynamics might be regulated. Molecular noise was effective in relocating clusters when the clusters contained low numbers of limiting polarity factors, and when Cdc42, the central polarity regulator, exhibited short dwell times at the polarity site. Cluster stabilization occurred when abundances or binding rates were altered to either lengthen dwell times or increase the number of polarity molecules in the cluster. We validated key results using full 3D particle-based simulations. Understanding the mechanisms cells use to regulate the dynamics of polarity clusters should provide insights into how cells dynamically track external directional cues. Cells localize polarity molecules in a small region of the plasma membrane forming a polarity cluster that directs functions such as migration, reproduction, and growth. Guided by external signals, these clusters move across the membrane allowing cells to reorient growth or motion. The polarity molecules continuously and randomly shuttle between the cluster and the cell cytosol and, as a result, the number and distribution of molecules at the cluster constantly changes. Here we present an improved stochastic simulation algorithm to investigate how such molecular-scale fluctuations induce cluster movement across the cell membrane. Unexpectedly, cluster mobility does not correlate with variations in total molecule abundance within the cluster, but rather with changes in the spatial distribution of molecules that form the cluster. Cluster motion is faster when polarity molecules are scarce and when they shuttle rapidly between the cluster and the cytosol. Our results suggest that cells control cluster mobility by regulating the abundance of polarity molecules and biochemical reactions that affect the time molecules spend at the cluster. We provide insights into how cells harness random molecular behavior to perform functions important for survival, such as detecting the direction of external signals.
Collapse
Affiliation(s)
- Samuel A. Ramirez
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- * E-mail: (SAR); (TCE)
| | - Michael Pablo
- Department of Chemistry, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- Program in Molecular and Cellular Biophysics, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Sean Burk
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Daniel J. Lew
- Department of Pharmacology and Cancer Biology, Duke University, Durham, North Carolina, United States of America
| | - Timothy C. Elston
- Department of Pharmacology and Computational Medicine Program, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
- * E-mail: (SAR); (TCE)
| |
Collapse
|
3
|
Ramirez SA, Pablo M, Burk S, Lew DJ, Elston TC. A Novel Stochastic Simulation Approach Enables Exploration of Mechanisms to Regulate Polarization Dynamics. Biophys J 2020. [DOI: 10.1016/j.bpj.2019.11.865] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/01/2022] Open
|
4
|
Abstract
After photoactivation, rhodopsin (R), a G-protein-coupled receptor, rapidly activates multiple transducin G-proteins (G) in an initial amplification step of phototransduction. G-protein activation requires diffusion-mediated association with an active rhodopsin (R*) at the rod disk membrane. Different organizations of R within the membrane have been revealded by several microscopy studies, including static and freely diffusing situations. However, it is unclear how such different scenarios influence the activation rate of G proteins. Through Monte Carlo simulations, we study the association reaction between a photoactivated rhodopsin and transducin under different reported receptor organizations including (a) R monomers diffusing freely, (b) R forming static dispersed crystalline domains made of rows of dimers, and (c) R arranged in static tracks formed by two adjacent rows of dimers. A key parameter in our simulations is the probability of binding following a collision ( p). For high p, the association rate between R* and G is higher in the freely diffusive system than in the static organizations, but for low collision efficiencies, the static organizations can result in faster association rates than the mobile system. We also observe that with low p, increasing the concentration of R increases the association rate significantly in the dispersed crystals configuration and just slightly in the free diffusive system. In summary, the lateral organization of rhodopsin influences the association rate between R* and G in a manner strongly dependent on the collision efficiency.
Collapse
Affiliation(s)
- Samuel A Ramirez
- Department of Pharmacology , University of North Carolina at Chapel Hill , Chapel Hill , North Carolina 27599 ; United States
| | - Chad Leidy
- Department of Physics , Universidad de los Andes , Bogotá 111711 , Colombia
| |
Collapse
|
5
|
Chiou JG, Ramirez SA, Elston TC, Witelski TP, Schaeffer DG, Lew DJ. Principles that govern competition or co-existence in Rho-GTPase driven polarization. PLoS Comput Biol 2018; 14:e1006095. [PMID: 29649212 PMCID: PMC5916526 DOI: 10.1371/journal.pcbi.1006095] [Citation(s) in RCA: 45] [Impact Index Per Article: 7.5] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/29/2018] [Revised: 04/24/2018] [Accepted: 03/16/2018] [Indexed: 11/19/2022] Open
Abstract
Rho-GTPases are master regulators of polarity establishment and cell morphology. Positive feedback enables concentration of Rho-GTPases into clusters at the cell cortex, from where they regulate the cytoskeleton. Different cell types reproducibly generate either one (e.g. the front of a migrating cell) or several clusters (e.g. the multiple dendrites of a neuron), but the mechanistic basis for unipolar or multipolar outcomes is unclear. The design principles of Rho-GTPase circuits are captured by two-component reaction-diffusion models based on conserved aspects of Rho-GTPase biochemistry. Some such models display rapid winner-takes-all competition between clusters, yielding a unipolar outcome. Other models allow prolonged co-existence of clusters. We investigate the behavior of a simple class of models and show that while the timescale of competition varies enormously depending on model parameters, a single factor explains a large majority of this variation. The dominant factor concerns the degree to which the maximal active GTPase concentration in a cluster approaches a "saturation point" determined by model parameters. We suggest that both saturation and the effect of saturation on competition reflect fundamental properties of the Rho-GTPase polarity machinery, regardless of the specific feedback mechanism, which predict whether the system will generate unipolar or multipolar outcomes.
Collapse
Affiliation(s)
- Jian-Geng Chiou
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, North Carolina, United States of America
| | - Samuel A. Ramirez
- Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Timothy C. Elston
- Department of Pharmacology, University of North Carolina at Chapel Hill, Chapel Hill, North Carolina, United States of America
| | - Thomas P. Witelski
- Department of Mathematics, Duke University, Durham, North Carolina, United States of America
| | - David G. Schaeffer
- Department of Mathematics, Duke University, Durham, North Carolina, United States of America
| | - Daniel J. Lew
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, North Carolina, United States of America
- * E-mail:
| |
Collapse
|
6
|
Abstract
Polarity establishment, the spontaneous generation of asymmetric molecular distributions, is a crucial component of many cellular functions. Saccharomyces cerevisiae (yeast) undergoes directed growth during budding and mating, and is an ideal model organism for studying polarization. In yeast and many other cell types, the Rho GTPase Cdc42 is the key molecular player in polarity establishment. During yeast polarization, multiple patches of Cdc42 initially form, then resolve into a single front. Because polarization relies on strong positive feedback, it is likely that the amplification of molecular-level fluctuations underlies the generation of multiple nascent patches. In the absence of spatial cues, these fluctuations may be key to driving polarization. Here we used particle-based simulations to investigate the role of stochastic effects in a Turing-type model of yeast polarity establishment. In the model, reactions take place either between two molecules on the membrane, or between a cytosolic and a membrane-bound molecule. Thus, we developed a computational platform that explicitly simulates molecules at and near the cell membrane, and implicitly handles molecules away from the membrane. To evaluate stochastic effects, we compared particle simulations to deterministic reaction-diffusion equation simulations. Defining macroscopic rate constants that are consistent with the microscopic parameters for this system is challenging, because diffusion occurs in two dimensions and particles exchange between the membrane and cytoplasm. We address this problem by empirically estimating macroscopic rate constants from appropriately designed particle-based simulations. Ultimately, we find that stochastic fluctuations speed polarity establishment and permit polarization in parameter regions predicted to be Turing stable. These effects can operate at Cdc42 abundances expected of yeast cells, and promote polarization on timescales consistent with experimental results. To our knowledge, our work represents the first particle-based simulations of a model for yeast polarization that is based on a Turing mechanism. Many cells need to generate and maintain biochemical signals in specific subcellular regions. This phenomenon is broadly called polarity establishment, and is important in fundamental processes such as cell migration and differentiation. A key polarity factor found in diverse organisms, including yeast and humans, is the protein Cdc42. In yeast, Cdc42-dependent polarization occurs through a self-reinforcing biochemical signaling loop. Directional cues can guide polarity establishment, but interestingly, yeast can polarize in the absence of such a cue. The mechanism thought to underlie this symmetry breaking involves the amplification of inhomogeneities in molecular distributions that arise from molecular-level fluctuations. We investigated the effects of random fluctuations on polarization by performing particle-based simulations of the Cdc42 signaling network. We found that fluctuations can facilitate polarization, allowing faster polarization, and polarization over a broader range of concentrations. Our observations may help understand how polarity works in other systems.
Collapse
Affiliation(s)
- Michael Pablo
- Department of Chemistry, The University of North Carolina, Chapel Hill, NC, United States of America
- Program in Molecular and Cellular Biophysics, The University of North Carolina, Chapel Hill, NC, United States of America
| | - Samuel A. Ramirez
- Department of Pharmacology, The University of North Carolina, Chapel Hill, NC, United States of America
| | - Timothy C. Elston
- Department of Pharmacology, The University of North Carolina, Chapel Hill, NC, United States of America
- * E-mail:
| |
Collapse
|
7
|
Ramirez SA, Raghavachari S, Lew DJ. Dendritic spine geometry can localize GTPase signaling in neurons. Mol Biol Cell 2015; 26:4171-81. [PMID: 26337387 PMCID: PMC4710246 DOI: 10.1091/mbc.e15-06-0405] [Citation(s) in RCA: 22] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2015] [Revised: 08/13/2015] [Accepted: 08/25/2015] [Indexed: 12/02/2022] Open
Abstract
Dendritic spines are the postsynaptic terminals of most excitatory synapses in the mammalian brain. Learning and memory are associated with long-lasting structural remodeling of dendritic spines through an actin-mediated process regulated by the Rho-family GTPases RhoA, Rac, and Cdc42. These GTPases undergo sustained activation after synaptic stimulation, but whereas Rho activity can spread from the stimulated spine, Cdc42 activity remains localized to the stimulated spine. Because Cdc42 itself diffuses rapidly in and out of the spine, the basis for the retention of Cdc42 activity in the stimulated spine long after synaptic stimulation has ceased is unclear. Here we model the spread of Cdc42 activation at dendritic spines by means of reaction-diffusion equations solved on spine-like geometries. Excitable behavior arising from positive feedback in Cdc42 activation leads to spreading waves of Cdc42 activity. However, because of the very narrow neck of the dendritic spine, wave propagation is halted through a phenomenon we term geometrical wave-pinning. We show that this can account for the localization of Cdc42 activity in the stimulated spine, and, of interest, retention is enhanced by high diffusivity of Cdc42. Our findings are broadly applicable to other instances of signaling in extreme geometries, including filopodia and primary cilia.
Collapse
Affiliation(s)
- Samuel A Ramirez
- Program in Computational Biology and Bioinformatics, Duke University Medical Center, Durham, NC 27710 Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, NC 27710
| | | | - Daniel J Lew
- Department of Pharmacology and Cancer Biology, Duke University Medical Center, Durham, NC 27710
| |
Collapse
|
8
|
Ramirez SA, Leidy C. Monte Carlo Study of the Association Rate Between Transducin and Photoactivated Rhodopsin at Disc Membranes. Biophys J 2014. [DOI: 10.1016/j.bpj.2013.11.1774] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/15/2022] Open
|
9
|
Chan C, Billard M, Ramirez SA, Schmidl H, Monson E, Kepler TB. A model for migratory B cell oscillations from receptor down-regulation induced by external chemokine fields. Bull Math Biol 2013; 75:185-205. [PMID: 23296998 PMCID: PMC3547247 DOI: 10.1007/s11538-012-9799-9] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.1] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2012] [Accepted: 11/15/2012] [Indexed: 01/28/2023]
Abstract
A long-standing paradigm in B cell immunology is that effective somatic hypermutation and affinity maturation require cycling between the dark zone and light zone of the germinal center. The cyclic re-entry hypothesis was first proposed based on considerations of the efficiency of affinity maturation using an ordinary differential equations model for B cell population dynamics. More recently, two-photon microscopy studies of B cell motility within lymph nodes in situ have revealed the complex migration patterns of B lymphocytes both in the preactivation follicle and post-activation germinal center. There is strong evidence that chemokines secreted by stromal cells and the regulation of cognate G-protein coupled receptors by these chemokines are necessary for the observed spatial cell distributions. For example, the distribution of B cells within the light and dark zones of the germinal center appears to be determined by the reciprocal interaction between the level of the CXCR4 and CXCR5 receptors and the spatial distribution of their respective chemokines CXCL12 and CXCL13. Computer simulations of individual-based models have been used to study the complex biophysical and mechanistic processes at the individual cell level, but such simulations can be challenging to parameterize and analyze. In contrast, ordinary differential equations are more tractable, but traditional compartment model formalizations ignore the spatial chemokine distribution that drives B cell redistribution. Motivated by the desire to understand the motility patterns observed in an individual-based simulation of B cell migration in the lymph node, we propose and analyze the dynamics of an ordinary differential equation model incorporating explicit chemokine spatial distributions. While there is experimental evidence that B cell migration patterns in the germinal center are driven by extrinsically regulated differentiation programs, the model shows, perhaps surprisingly, that feedback from receptor down-regulation induced by external chemokine fields can give rise to spontaneous interzonal and intrazonal oscillations in the absence of any extrinsic regulation. While the extent to which such simple feedback mechanisms contributes to B cell migration patterns in the germinal center is unknown, the model provides an alternative hypothesis for how complex B cell migration patterns might arise from very simple mechanisms.
Collapse
Affiliation(s)
- Cliburn Chan
- Department of Biostatistics and Bioinformatics, Duke University Medical Center, Durham, NC 27705, USA.
| | | | | | | | | | | |
Collapse
|
10
|
Goel VK, Ramirez SA, Kong W, Gilbertson LG. Cancellous bone Young's modulus variation within the vertebral body of a ligamentous lumbar spine--application of bone adaptive remodeling concepts. J Biomech Eng 1995; 117:266-71. [PMID: 8618378 DOI: 10.1115/1.2794180] [Citation(s) in RCA: 83] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [What about the content of this article? (0)] [Affiliation(s)] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023]
Abstract
Bone remodeling theory based on strain energy density (SED) as the feedback control variable was used in conjunction with the finite element method to analyze the shape of the vertebral bodies within the ligamentous motion segment. The remodeling theory was once again applied to the altered two motion segments model to predict the Young's modulus distribution of the cancellous bone within the vertebral bodies. A three-dimensional finite element model of the two motion segments ligamentous lumbar spine (L3-5) was developed. Bone remodeling response (external as well as internal) of the motion segments to a uniaxial compressive load of 424.7 N was studied. The external shape of the converged model matched the normal shape of a vertebral body. The internal remodeling resulted in regional cancellous bone Young's moduli (or bone density) distributions similar to those reported in the literature; posterocentral regions of the vertebrae were predicted to have greater values of the elastic modulus than that of the outer regions. The results of the present study suggest that vertebral body assumes an adequate/optimum structure in terms of both its shape and its elastic moduli distribution within the cancellous region in response to the applied load. Extensions of the present model and its clinically relevant applications are discussed.
Collapse
Affiliation(s)
- V K Goel
- Iowa Spine Research Center, Department of Biomedical Engineering, College of Engineering, University of Iowa, Iowa City 52242, USA
| | | | | | | |
Collapse
|
11
|
Daskal I, Ramirez SA, Ballal RN, Spohn WH, Wu B, Busch H. Detergent lysis for isolation of intact polysomes of Nivikoff hepatoma ascites cells. Cancer Res 1976; 36:1026-34. [PMID: 1253165] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [Abstract] [MESH Headings] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/26/2022]
Abstract
An efficient, rapid procedure for the isolation of polysomes from Novikoff hepatoma ascites cells that avoids mechanical shear was developed, and Ivory was used as the detergent to lyse the cells. This method yields classes of large polysomes and 3- to 5-fold increased yields over that of routine homogenization procedures. Up to 70% of the polysome population is larger than tetramers. The polysomes and ribosomes are morphologically intact and contain undegraded ribosomal RNA and biologically active messenger RNA. This procedure also yields nuclei free of cytoplasmic contamination as a by-product of polysome isolation. These nuclei do not show evidence of protein degradation and are enzymatically active with respect to RNA polymerases I and II.
Collapse
|
12
|
Abstract
ABSTRACT
Cytogenetic and molecular hybridization techniques were used to examine the arrangement of rRNA genes at the nucleolar organizer region (NOR) of maize. TB-6a stocks involving a reciprocal translocation between a B chromosome and the NOR of chromosome 6 were used. The amount of rDNA in the different stocks used varied from 1.5 to 4 NOR equivalents depending on the number of B6 choromosomes present. Cytological measurements show that the medial break through the NOR results in an equal partitioning of the heterochromatic knob (the NOR) between chromosome 6 and the B chromosome to which it was translocated. DNA-rRNA hybridization experiments show a linear relationship between the amount of rRNA capable of hybridizing with DNA and the number of NOR equivalents present. The data confirm McClintock's conclusion that the heterochromatic knob, rather than the constricted portion, is the true NOR region. Further, they show that the number of ribosomal genes is correlated with the amount of cytologically visible NOR, suggesting a uniform distribution of genes throughout the locus.
Collapse
Affiliation(s)
- S A Ramirez
- Zoology Department, Indiana University, Bloomington, Indiana 47401
| | | |
Collapse
|
13
|
Abstract
ABSTRACT
Ribosomal genes in eukaryotes are highly redundant. Considerable variation in the level of redundancy among species, especially in higher plants, has been reported; but except for deletion and duplication mutants, it is generally accepted that intraspecific variability in redundancy level is small. We have examined the level of redundancy in several lines of maize by DNA-rRNA saturation hybridization. The amount of nuclear DNA which hybridizes with rRNA in the ten lines examined varied from 0.24% to 0.50%. The number of rRNA genes per diploid genome thus ranges from 1.12 × 104 to 2.32 × 104. Results also indicate that the level of redundancy is genetically transmitted.
Collapse
Affiliation(s)
- S A Ramirez
- Zoology Department, Indiana University, Bloomington, Indiana 47401
| | | |
Collapse
|
14
|
Ramirez AA, Ramirez SA, Malbas G, De Jesus G, Yazon F, Mallonga A. Cod liver oil vs. burns. J Philipp Med Assoc 1966; 42:521-8. [PMID: 5976613] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [What about the content of this article? (0)] [MESH Headings] [Subscribe] [Scholar Register] [Indexed: 01/17/2023]
|