key: cord-0884773-msbplvg9 authors: Ebrahim-Habibi, Azadeh; Kashani-Amin, Elaheh; Larijani, Bagher title: Modeling and simulation in medical sciences: an overview of specific applications based on research experience in EMRI (Endocrinology and Metabolism Research Institute of Tehran University of Medical Sciences) date: 2021-01-22 journal: J Diabetes Metab Disord DOI: 10.1007/s40200-020-00706-x sha: 7a8602a9fdd808bc013dffe0c69c125be0b2b109 doc_id: 884773 cord_uid: msbplvg9 The concomitant use of various types of models (in silico, in vitro, and in vivo) has been exemplified here within the context of biomedical researches performed in the Endocrinology and Metabolism Research Institute (EMRI) of Tehran University of Medical Sciences. Two main research aeras have been discussed: the search for new small molecules as therapeutics for diabetes and related metabolic conditions, and diseases related to protein aggregation. Due to their multidisciplinary nature, the majority of these studies have needed the collaboration of different specialties. In both cases, a brief overview of the subject is provided through literature examples, and sequential use of these methods is described. "Model" is a generic word, present in the majority (if not all) of scientific disciplines technical vocabulary, and could be represented by a simple diagram that is easily understood by a glance. It could also be a complex entity incorporating hundreds of thousands of interacting partners, that should be analyzed by an information processing system. Either simple or complex, a model is always an approximation, contains errors, and is usually a tool to be used alongside with others. Even within a particular scientific branch, such as medical sciences, a "model" takes different forms. As an example, in the present context of COVID- 19 , there has been a worldwide effort to predict the spread and dynamics of the disease [1] [2] [3] [4] , and a similar international endeavor toward finding potential drugs that would bind the virus proteins [5] [6] [7] . In both cases, sophisticated "models" were used, which have fundamentally different characteristics: one includes mathematical equations and the other is an atomistic representation of proteins. The first modeling works that have been performed in EMRI had a mathematical nature and included modeling of angiogenesis in tumors (aiming at improving diagnostic) [8] , and proposing a suitable time to test bone mass density in post-menopausal women receiving levothyroxine [9] . In the following years, different kinds of mathematical models have been used in numerous studies related with public health (e.g. [10] ) and other clinically-oriented subjects such as osteoporosis (e.g. adjusting a predictive tool to assess fracture risk in osteoporotic patients for the Iranian population) [11] . In this condensed overview, examples of specific research subjects related to medical sciences are discussed with an emphasis on the use of in vivo, in vitro, and in silico models. These subjects have been investigated by students and researchers over the past twelve years in the Modeling and Simulation in Medical Sciences (MSMS) research group of EMRI with the help and contribution of numerous collaborators inside and outside EMRI. We believe that models are important tools that upon extensive development, could reduce the overall cost of research, and ideally replace the need for extended testing. In a recent study concerning the global burden of diseases, increase of the mortality rate and total years of life lost (YLL) due to diabetes mellitus has been observed (around 31 and 25% respectively) between 2006 and 2016, while YLL itself had decreased (around 2%) [12] . A worldwide rise in the prevalence of diabetes is likely, estimated to reach 10.4% by 2040 [13] , and justifies the search for novel therapies to control the disease and prevent its complications. Various molecular targets have been proposed and explored in this regard, with the aim of improving insulin's effect (e.g. by influencing beta-cells function by modulating related receptors), as well as controlling oxidative stress and inflammation to counteract beta-cells apoptosis [14, 15] . A different approach consists of interfering with key carbohydrate digestive enzymes (glycosidases) that are responsible for the postprandial rise of blood glucose level, especially alphaglucosidase [16] , and alpha-amylase [17] . Inhibitors of these enzymes could also be potentially effective in obesity. Many studies are now reporting the glycosidase inhibitory potential of plants and their active components [18] ; earlier studies had also demonstrated the efficacy of flavonoids in this regard [19, 20] , and analyzed putative binding modes of these compounds to the enzyme by in silico techniques [21] . During in vitro experiments, with the use of pancreatic alpha-amylase, and based on the fact that the natural compounds flavonoids were known as amylase inhibitors [22] , we found trans-chalcone (flavonoids precursor) to be an inhibitor of pancreatic amylase [23] . Further experiments showed it to be also effective in vivo, on a rodent diabetic model [24] and in a mice model of fatty-liver, where the compound was able to modulate lipid profile, leptin, and glucose levels, as well as improving liver steatosis [25] ; these results were comparable to what we had seen in a mice model of obesity [26] . Trans-chalcone is actually considered as a potent scaffold, whose numerous derivatives span a wide range of therapeutic properties [27, 28] . The benzothiazole thioflavin T (ThT), with slighter inhibitory activity in vitro, was also found to be potentially antidiabetic in vivo [29] , as well as effective in weight reduction of a mice model of obesity [30] . ThT is a well-known probe for detection of amyloid structures [31] , but interestingly, had also been found to extend the lifespan of C. elegans (as a model for studying ageing and lifespan) [32, 33] , and seems to be an interesting therapeutic scaffold. Synthetic aurones ((Z)-2-benzylidenebenzofuran-3-one derivatives), which bear structural similarity to flavonoids were also tested in vitro against pancreatic amylase as inhibitors, and subsequent in silico tests (docking) suggested their putative binding modes into the enzyme [34] . Docking was done with the use of U-Dock 1.6 program [35] , with the MMFF94s force field and on MOE (Molecular Operating Environment 2012.10) and the target was PDB entry 1OSE [36] In this method, a series of ligands conformers were docked onto the target and poses with best energy were retrieved. This data could be of use for synthesizing more potent inhibitors. Recently, many natural extracts and synthetic new molecules have been investigated with regard to their alpha-amylase inhibitory properties [37] [38] [39] [40] [41] [42] [43] [44] [45] [46] [47] [48] and this line of study seems to be promising with regard to finding effective anti-diabetic aids. On the other hand, activators of the enzyme may potentially show an adverse effect in this regard, although hypothetically beneficial for digestive problems. Interestingly, the commercially available sweetener neohesperidin dihydochalcone (from sweet orange), which has a chalcone structural component was found to activate alpha-amylase from different sources [49, 50] . Xanthine derivatives, pentoxyfilline, theobromine and caffeine were also found to slightly (20-30%) increase enzyme activity [51] . To our knowledge, these should be the only small molecules that have been reported as activator of this enzyme. In another context, sweet taste receptor (STR), a class C Gprotein coupled receptor (C GPCR), has been suggested as a drug target for designing either new low-calorie sweeteners or drugs to control metabolic condition disorders such as type II diabetes mellitus [52, 53] . In order to design appropriate ligands for this receptor, a three-dimensional structure is needed, but since the human receptor structure is still not experimentally elucidated, in silico (molecular modeling) methods were used in this case. First, a new model of STR structure was developed [54] . Extensive docking experiments on the model revealed presence of carbohydrate binding modules as structural and functional motifs [55] , which is applicable in further characterization of the binding site. Docking was performed on an ensemble of STR structures obtained by running a molecular dynamics simulation (MDS) experiment. MDS was run for 50 ns in YASARA program [56] using the YASARA forcefield [57] . The initial structure atoms were protonated at physiological pH, Particle Mesh Ewald (PME) was used, 8.0 A˚cutoff for non-bonded interactions and 4 fs timestep were applied, hydrogen atoms were constrained and the system was at constant pressure and temperature (NPT). Based on carbon alpha RMSDs of the total MDS run, clusters were obtained based on minimum, median and maximum values. Docking was performed on these clusters with AutoDock Vina [58] on the interface present in YASARA which enabled computing H-bonds and HP-interactions, as well as minimization of the best poses to take water molecules into account. Finally, by studying interaction patterns of 316 sweet molecules as well as sweet proteins and antagonists of the receptor, a few compounds were suggested as putative sweet molecules or agonists for the receptor [unpublished data]. Increasing our knowledge about this receptor structure/ function and its agonists/antagonists would surely be of help in the current debate about the beneficial/detrimental effects of sweeteners [59-62]. The integrity of proteins three-dimensional structures is a requirement for their normal functioning. Over disruption of their intramolecular interactions, "misfolded" proteins structures could be formed, which, depending on environmental conditions, may "aggregate" together and form stable clusters. Misfolding may happen upon mutations or other external causes, and lead to "protein diseases". Aggregates may be structured (amyloid) or unstructured (amorphous), although various intermediate forms can also be observed [63] . Amyloid formation is suggested to have a pathogenic role in a range of disorders, including Azheimer's and Parkinson disease, metabolic diseases, and diabetes [64] . In recent years, formation of structured amyloid forms has been suggested to be a generic property of proteins, which means that any protein could be potentially driven toward forming amyloids [65] , and as so, be used as a model to study the process in details. Inhibiting amyloid formation or disruption of amyloid structures by small ligands has been the subject of numerous studies in the past two decades, and based on these reports, there have been attempts to classify ligands and identify the structural prerequisites that make suitable anti-amyloid compounds [66] [67] [68] [69] [70] [71] [72] . One of the research subjects that could provide useful information for drug design is the study of structural changes that occur in proteins and result onto amyloid formation. Disruption of the native and functional architecture of the protein may start from specific locations, and aggregationprone sites in proteins have an important role in triggering aggregates formation. Predicting those regions has been (and still is) the subject of multiple computational studies which have resulted onto design of various tools [73] [74] [75] [76] [77] [78] [79] [80] [81] [82] [83] . An alternative to these tools is the application of molecular dynamics simulation methods to proteins and monitor the changes that happen in their structures under various deleterious conditions. We have applied this method to insulin [84] and two forms of myoglobin (native and glycated) [85] and could spot specific regions that seem to be initiators in the structural change process. In both cases, high temperatures were employed in order to reach the unfolded state more easily [86] . For insulin, different conditions (presence and absence of KCL and NaCL) were tested, at high temperature and acidic pH (obtained by changing the protonation state of the protein). For the actual run, MOE (Molecular Operating Environment MOE.2010.10), and the implemented MMFF94x force field were used for 15 ns runs of each condition. To check the conformational change over time, secondary structure content was analyzed by using various methods including DSSP [87] , STRIDE [88] , PALSSE [89] , P-SEA [90] , STICK [91] , and XTLSSTR [92] . For myoglobin too, the same program (MOE.2010.10), and the implemented MMFF94x force field were used. In this case, in addition to secondary structure content (computed by DSSP), structures' RMSDs were also compared in order to find the most stable regions. Our subsequent ongoing research projects focus on shorter proteins with the aim of pinpointing specific residues in well-defined secondary structure segments. These in silico methods could also provide a mean to screen potential stabilizing ligands that could counteract these structural changes to some extent [unpublished results]. In vitro methods provide the possibility to generate amyloid forms of various proteins, and to test potential antiamyloid molecules. Phenolic, polyphenolic, and other cyclic compounds have been tested in our past (collaborative) studies on various proteins including polygalacturonase, albumin [93] , myoglobin [94] , and insulin [95] . Molecules that were found to inhibit proteins fibrillation, were then tested in vivo, on a rodent model of AD. In order to simulate part of the AD signs, abeta amyloid is injected to the animals hippocampus, which results into plaque formation and troubles in learning. Candidate ligands that were administered to the AD animal models, and could attenuate their symptoms to some extent included indole and trans-chalcone [96] , sylimarin [97] , metformin [98] , thymol [99] , and eugenol [100] . Based on our experience in generating amyloid structures in insulin, we have also proposed a modification in the generation of rodent AD-like model : instead of injecting abeta amyloids to the animals brains, we used insulin amyloids which are easier formed at a lower cost [101] . Insulin amyloid formation is a problem in patients who have to inject the protein subcutaneously to control their blood glucose levels. We found that chemical modification of insulin by anhydrides at its lysine residue could attenuate its amyloid formation, while insulin still remained functional with one of the modifiers (succinic anhydride) and with another, a delayed onset of insulin action was observed [102] . In another set of experiments, we simulated the condition that is observed in diabetic patients by injecting insulin amyloid to rodents. This resulted onto formation of masses composed of adipose cells in which amyloid deposits could be observed [103] . A subsequent test showed turmeric to be effective in reducing the size of those masses [104] and further studies are now conducted in order to further develop this model and check the effect of other compounds. Recently, our method was used by other researchers to assess the effect of serine proteases on amyloid fibrils [105] . In the two main research subjects that have been discussed here, interdisciplinary work was needed, which is strongly dependent on extensive collaborations; whenever the projects included a higher diversity of specialty, more interesting viewpoints were found, and overall better results were obtained. In another set of researches (not discussed here), we had the opportunity to collaborate with clinicians and geneticists, where we could provide a "basic science" insight onto genetic diseases that were manifesting themselves at protein level. After a case was identified, the faulty gene was investigated, and we had a look at the (computed) protein structure of the mutated gene. This experience led us to define a new research section where we try to see whether our used in silico methods can provide more precise information on the proteins' mutation effects on its structure. While in silico and in vitro methods have serious limitations, and in vivo disease models are still far from the reality of pathogenic conditions in human bodies, it is still possible to relate (to some extent) the data obtained from these different methods. We are still a long way from the ideal situation where everything could be computed and there would be no need to perform other experiments, but in order to get closer to that aim, we need to gather more data, and continue to interact between disciplines. Modeling and forecasting trend of COVID-19 epidemic in Iran until May 13, 2020 Optimization method for forecasting confirmed cases of COVID-19 in China Forecasting the novel coronavirus COVID-19 Modeling and forecasting the COVID-19 pandemic in India Drug repositioning an alternative for the treatment of coronavirus COVID-19 Current status of COVID-19 therapies and drug repositioning applications Drug repurposing using computational methods to identify therapeutic options for COVID-19 A stochastic model of tumor angiogenesis Modeling the effect of levothyroxine therapy on bone mass density in postmenopausal women: a different approach leads to new inference Improved population health in Iran from 1979 to 2019; Decreasing mortality rates and increasing life expectancy The risk of osteoporotic fractures and its associating risk factors according to the FRAX model in the Iranian patients: a follow-up cohort Global, regional, and national age-sex specific mortality for 264 causes of death, 1980-2016: a systematic analysis for the Global Burden of Disease Study IDF Diabetes Atlas: Global estimates for the prevalence of diabetes for 2015 and 2040 Novel approaches to drug discovery for the treatment of type 2 diabetes The hot and potential targets of type 2 diabetes mellitus treatment in recent decade Acarbose: safe and effective for lowering postprandial hyperglycaemia and improving cardiovascular outcomes Inhibitors of alpha-amylase from plants-a possibility to treat diabetes mellitus type II by phytotherapy Natural products for the treatment of type 2 diabetes mellitus Inhibition of alpha-glucosidase and amylase by luteolin, a flavonoid Comparative evaluation of quercetin, isoquercetin and rutin as inhibitors of α-glucosidase Flavonoids for controlling starch digestion: structural requirements for inhibiting human α-amylase Inhibition of αglucosidase and α-amylase by flavonoids Trans-chalcone: a novel small molecule inhibitor of mammalian alpha-amylase Core structure of flavonoids precursor as an antihyperglycemic and antihyperlipidemic agent: an in vivo study in rats Effect of trans-chalcone on atheroma plaque formation, liver fibrosis and adiponectin gene expression in cholesterol-fed NMRI mice Acarbose versus trans-chalcone: comparing the effect of two glycosidase inhibitors on obese mice Chalcone: a privileged structure in medicinal chemistry Chalcone derivatives: promising starting points for drug design Thioflavin T effect in diabetic Wistar rats: reporting the antihyperglycemic property of an amyloid probing dye Benzothiazole Thioflavin T improves obesity-related symptoms in mice Binding mode of Thioflavin T and other molecular probes in the context of amyloid fibrils-current status Effects of thioflavin T and GSK-3 inhibition on lifespan and motility in a caenorhabditis elegans model of tauopathy Amyloid-binding compounds maintain protein homeostasis during ageing and extend lifespan Aurones as new porcine pancreatic α-Amylase inhibitors Design, synthesis, and biological evaluation of substituted 2-alkylthio-1, 5-diarylimidazoles as selective COX-2 inhibitors Crystal structure of pig pancreatic α-amylase isoenzyme II, in complex with the carbohydrate inhibitor acarbose Novel curcumin-based pyrano [2, 3-d] pyrimidine anti-oxidant inhibitors for α-amylase and α-glucosidase: Implications for their pleiotropic effects against diabetes complications Curcumin-based antioxidant and glycohydrolase inhibitor compounds: Synthesis and in vitro appraisal of the dual activity against diabetes Bin Asad MHH. Enzymes inhibitors from natural sources with antidiabetic activity: A review α-Amylase and α-glucosidase inhibitors from Zanthoxylum chalybeum Engl. root bark Antimicrobial and alpha-amylase inhibitory activities of organic extracts of selected Sri Lankan bryophytes An insight of alpha-amylase inhibitors as a valuable tool in the management of type 2 diabetes mellitus Alpha-glucosidase and alpha-amylase inhibitory activities of novel abietane diterpenes from Salvia africana-lutea Anti-oxidant activity, alpha-amylase inhibition and toxicity of leaf extracts of cultivated Rapanea melanophloeos (L.) Mez (cape beech) Potential inhibitors of î±-glucosidase and î±-amylase enzymes from locally avaialable fruit wastes by solid state fermentation Alphaglucosidase and amylase inhibitory effects of Eruca vesicaria subsp. longirostris essential oils: synthesis of new 1, 2, 4-triazole-thiol derivatives and 1, 3, 4-thiadiazole with potential inhibitory activity Inhibitory activity of phenolic-rich pistachio green hull extract-enriched pasta on key type 2 diabetes relevant enzymes and glycemic index Inhibitory effects of solvent-partitioned fractions of two nigerian herbs (spondias mombin linn. and mangifera indica L.) on αamylase and α-glucosidase Effect of neohesperidin dihydrochalcone on the activity and stability of alpha-amylase: a comparative study on bacterial, fungal, and mammalian enzymes Neohesperidin dihydrochalcone: presentation of a small molecule activator of mammalian alpha-amylase as an allosteric effector Xanthine derivatives as activators of alpha-amylase: hypothesis on a link with the hyperglycemia induced by caffeine Functional roles of the sweet taste receptor in oral and extraoral tissues Taste receptors in the gastrointestinal system Introducing a new model of sweet taste receptor, a class C G-protein Coupled Receptor (C GPCR) Presence of carbohydrate binding modules in extracellular region of class C G-protein coupled receptors (C GPCR): An in silico investigation on sweet taste receptor YASARA View-molecular graphics for all devices-from smartphones to workstations New ways to boost molecular dynamics simulations AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading Current progress in understanding the structure and function of sweet taste receptor Non-caloric sweeteners, sweetness modulators, and sweetener enhancers The association between artificial sweeteners and obesity Artificial sweeteners: safe or unsafe? Protein aggregation: folding aggregates, inclusion bodies and amyloid Protein misfolding, amyloid formation, and human disease: a summary of progress over the last decade The amyloid state and its association with protein misfolding diseases Inhibition of protein aggregation and amyloid formation by small molecules Anti-amyloid aggregation activity of natural compounds: implications for Alzheimer's drug discovery Design of novel amyloid β aggregation inhibitors using QSAR, pharmacophore modeling, molecular docking and ADME prediction Small molecule probes of protein aggregation Natural product-based amyloid inhibitors Inhibition of amyloid-β aggregation in Alzheimer's disease Inhibition of amyloid fibrillation by small molecules and nanomaterials: strategic development of pharmaceuticals against amyloidosis Prediction of aggregation-prone regions in structured proteins Prediction of aggregation prone regions of therapeutic proteins Aggregation in protein-based biotherapeutics: computational studies and tools to identify aggregation-prone regions Aggregation prone regions and gatekeeping residues in protein sequences Predicting aggregation-prone sequences in proteins Solubis: a webserver to reduce protein aggregation through mutation Driving forces for nonnative protein aggregation and approaches to predict aggregation-prone regions Prediction of aggregation-prone regions in proteins based on the distribution of surface patches AGGRESCAN3D: Toward the prediction of the aggregation propensities of protein structures. Computational drug discovery and design ANuPP: A versatile tool to predict aggregation nucleating regions in peptides and proteins A molecular modeling study on full-length insulin: insight into initial events of amyloid formation Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study Characterization of the transition state of protein unfolding by use of molecular dynamics: chymotrypsin inhibitor 2 DSSP: definition of secondary structure of proteins given a set of 3D coordinates STRIDE: a web server for secondary structure assignment from known atomic coordinates of proteins PALSSE: A program to delineate linear secondary structural elements from protein structures SEA: a new efficient assignment of secondary structure from Cα trace of proteins Defining linear segments in protein structure Assigning secondary structure from protein coordinate data Amyloid formation and inhibition of an all-beta protein: a study on fungal polygalacturonase Heme binding site in apomyoglobin may be effectively targeted with small molecules to control aggregation Inhibition by aromatic compounds Indole and trans-chalcone attenuate amyloid β plaque accumulation in male Wistar rat: in vivo effectiveness of two anti-amyloid scaffolds Silymarin effect on amyloid-β plaque accumulation and gene expression of APP in an Alzheimer's disease rat model Neuroprotective effects of metformin against Aβ-mediated J Diabetes Metab Disord inhibition of long-term potentiation in rats fed a high-fat diet Investigation of thymol effect on learning and memory impairment induced by intrahippocampal injection of amyloid beta peptide in high fat diet-fed rats Effects of eugenol on alzheimer's disease-like manifestations in insulin-and Aβ-induced rat models Injection of insulin amyloid fibrils in the hippocampus of male Wistar rats: Report on memory impairment and formation of amyloid plaques Structure and function of anhydride-modified forms of human insulin: In silico, in vitro and in vivo studies Generating local amyloidosis in mice by the subcutaneous injection of human insulin amyloid fibrils Turmeric effect on subcutaneous insulininduced amyloid mass: an in vivo study Attenuation of subcutaneous insulin induced amyloid mass in vivo using Lumbrokinase and Serratiopeptidase Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgements The authors wish to thank Joulia Alizadeh-Rahrovi, Farnaz Banakar, and Homa Faraji for their valuable help. Conflict of interest The authors declare to have no conflict of interests.