Essential Dynamics Research Articles

Stable Communications in Green Unmanned Aerial Relaying Systems

Green unmanned aerial relaying (UAR) systems, featuring solar-powered unmanned aerial vehicles (UAVs) to provide agile and flexible communication services, have recently gained significant attention for their wide applications. Nonetheless, due to the essential dynamics of wireless channel conditions and solar energy supply, maintaining system stability is a critical concern in practical green UAR systems, as fluctuation in these two factors can significantly affect communication performance. To address this issue, this paper proposes a scheme design for stable communications in a UAR system with dynamic solar energy supply and wireless channel conditions. Specifically, we first formulate an optimization problem, called OFL, to control the power allocation and flow data rate assignment to minimize long-term time-averaged energy consumption while ensuring the demanded data rate and system stability. Considering that solving OFL needs complete prior information about the solar power supply and wireless channel conditions, which is hardly acquired, we then explore the Lyapunov theory to transform OFL as an online optimization problem, called ONL. To find the solution to ONL, we subsequently design a distributed algorithm that enables each UAV to make its own decision locally, thus reducing the computational overhead of each UAV. Particularly, we prove that stability can be guaranteed by employing the designed algorithm. Extensive simulations are conducted to show the performance achieved by the proposed scheme.

The kinetic Ising model encapsulates essential dynamics of land pattern change.

A land pattern change represents a globally significant trend with implications for the environment, climate and societal well-being. While various methods have been developed to predict land change, our understanding of the underlying change processes remains inadequate. To address this issue, we investigate the suitability of the two-dimensional kinetic Ising model (IM), an idealized model from statistical mechanics, for simulating land change dynamics. We test the IM on a variety of patterns, each with different focus land type. Specifically, we investigate four sites characterized by distinct patterns, presumably driven by different physical processes. Each site is observed on eight occasions between 2001 and 2019. Given the observed pattern at the time ti we find two parameters of the IM such that the model-evolved land pattern at ti+1 resembles the observed land pattern at that time. The data support simulating seven such transitions per site. Our findings indicate that the IM produces approximate matches to the observed patterns in terms of layout, composition, texture and patch size distributions. Notably, the IM simulations even achieve a high degree of cell-scale pattern accuracy in two of the sites. Nevertheless, the IM has certain limitations, including its inability to model linear features, account for the formation of new large patches and handle pattern shifts.

An exploration of binding of Hesperidin, Rutin, and Thymoquinone to acetylcholinesterase enzyme using multi-level computational approaches

Alzheimer’s disease, an intricate neurological disorder, is impacting an ever-increasing number of individuals globally, particularly among the aging population. For several decades phytochemicals were used as Ayurveda to treat both communicable and non-communicable diseases. Acetylcholinesterase (AChE) is a widely chosen therapeutic target for the development of early prevention and effective management of neurodegenerative diseases. The primary objective of the present study was to investigate the binding potential between Rutin Thymoquinone, Hesperidin and the FDA-approved drug Donepezil with AChE. Additionally, a comparative analysis was conducted. These phytochemicals were docked with the binding site of the AChE experimental complex. The molecular dockings demonstrated that the Hesperidinh showed a better binding affinity of −22.0631 kcal/mol. The ADME/T investigations revealed that the selected phytochemicals are non-toxic and drug-like candidates. Molecular dynamics simulations were implemented to determine the conformational changes of Rutin, hesperidin, Thymoquinone, and Donepezil complexed with AChE. Hesperidin and Donepezil were more stable than Rutin, Thymoquinone complexed with AChE. Next, essential dynamics and defining the secondary structure of protein were to determine the conformational changes in AChE complexed with selected phytochemicals during simulations. Overall, the MD Simulations demonstrated that all complexes in this study achieved stability until 100 ns of the simulation period was performed thrice. The structural analysis of AChE was done using multiple search engines to explore the molecular functions, biological processes, and pathways in which AChE proteins are involved and to identify potential drug targets for various diseases. This present study concludes that Hesperidin was found to be a more potent AChE inhibitors than Rutin, and further experiments are required to determine the effectivity of Hesperidin against neurodegenerative diseases. Communicated by Ramaswamy H. Sarma

Effect of cysteine mutation at Ca2+ coordinating residues to the autolysis, folding and hydrophobicity of full length and mature Rand protease: molecular dynamics simulation and essential dynamics

Rand protease is a serine protease that shared common characteristics with members of the MEROPS S8 subtilisin family. It is thermostable, highly stable in organic solvent and broad in specificity. Many structures of homologous protein solved by X-ray crystallography and NMR have been deposited to Protein Data Bank (PDB) which allowed this study to rely on structure prediction by deep learning to build three-dimensional (3D) structure of full length and mature Rand protease (flRP and mRP). In silico cysteine mutation to 7 predicted high affinity Ca2+ coordinating residues were introduced, and the mutants were subjected to molecular dynamics simulation to study its effect on flRP and mRP. MD simulation showed a marked increase in flexibility of the pro-peptide segment indicating the impact of single cysteine substitution at high affinity Ca2+ coordinating residues to autolysis of flRP. MD simulation for mRP reaffirmed the role of Ca2+ coordinating sites in providing stability to Rand protease. In addition, these residues also affect the autolysis, folding and hydrophobicity of RP. Essential dynamics observed large contribution of the first few eigenvectors of flRP, mRP and their high affinity Ca2+ coordinating residues mutants to the TMSF values which indicates that these values account for a large portion of the overall atomic fluctuations. These results have given a more comprehensive understanding on the role of cysteine substituted Ca2+ coordinating surface loop to the structure of flRP and mRP which are important in contributing to the structural stability of subtilisin. Communicated by Ramaswamy H. Sarma

A multifactorial framework of psychobehavioral determinants of coping behaviors: an online survey at the early stage of the COVID-19 pandemic.

Coronavirus disease 2019 dramatically changed people's behavior because of the need to adhere to infection prevention and to overcome general adversity resulting from the implementation of infection prevention measures. However, coping behavior has not been fully distinguished from risk perception, and a comprehensive picture of demographic, risk-perception, and psychobehavioral factors that influence the major coping-behavior factors remain to be elucidated. In this study, we recruited 2,885 Japanese participants. Major coping-behavior and risk-perception factors were identified via exploratory factor analysis of 50 candidate items. Then, we conducted a hierarchical multiple regression analysis to investigate factors associated with each coping-behavior factor. We identified four types of coping behavior [CB1 (mask-wearing), CB2 (information-seeking), CB3 (resistance to social stagnation), and CB4 (infection-prevention)] and three risk-perception factors [RP1 (shortages of daily necessities), RP2 (medical concerns), and RP3 (socioeconomic concerns)]. CB1 was positively associated with female sex and etiquette. CB2 was positively related to RP1 and RP3. CB3 was positively related to RP1 and leadership, and negatively associated with etiquette. CB4 was positively associated with female sex, etiquette, and active well-being. This parsimonious model may help to elucidate essential social dynamics and provide a theoretical framework for coping behavior during a pandemic.

Histologic Examination of Vocal Fold Mucosal Wave and Vibration.

Despite gross anatomic and histologic differences between human and canine vocal folds, similar wave patterns have been described yet not fully characterized. We reconstructed vocal fold (VF) vibration in a canine hemilarynx and performed histologic examination of the same vocal fold. We demonstrate comparable wave patterns while exploring the importance of certain anatomic architectures. An in vivo canine hemilarynx was phonated against a glass prism at low and high muscle activation conditions. Vibration was captured using high-speed video, and trajectories of VF medial surface tattooed landmarks were 3D-reconstructed. The method of empirical eigenfunctions was used to capture the essential dynamics of vibratory movement. Histologic examination of the hemilarynx was performed. Oscillation patterns were highly similar between the in vivo canine and previous reports of ex vivo human models. The two most dominant eigenfunctions comprised over 90% of total variance of movement, representing opening/closing and convergent/divergent movement patterns, respectively. We demonstrate a vertical phase difference during the glottal cycle. The time delay between the inferior and superior VF was greater during opening than closing for both activation conditions. Histological examination of canine VF showed not only a thicker lamina propria layer superiorly but also a distinct pattern of thyroarytenoid muscle fibers and fascicles as described in human studies. Histologic and vibratory examination of the canine vocal fold demonstrated human vocal fold vibratory patterns despite certain microstructural differences. This study suggests that the multilayered lamina propria may not be fundamental to vibratory patterns necessary for human-like voice production. NA (Basic science study) Laryngoscope, 134:264-271, 2024.

Differential gene expression analysis combined with molecular dynamics simulation study to elucidate the novel potential biomarker involved in pulmonary TB

Tuberculosis (TB) is a lethal multisystem disease that attacks the lungs' first line of defense. A substantial threat to public health and a primary cause of death is pulmonary TB. This study aimed to identify and investigate the probable differentially expressed genes (DEGs) primarily involved in Pulmonary TB. Accordingly, three independent gene expression data sets, numbered GSE139825, GSE139871, and GSE54992, were utilized for this purpose. The identified DEGs were used for bioinformatics-based analysis, including physical gene interaction, Gene Ontology (GO), network analysis and pathway studies using the Kyoto Encyclopedia of Genes and Genomes pathway (KEGG). The computational analysis predicted that TNFAIP6 is the significant DEG in the gene expression profiling of TB datasets. According to gene ontology analysis, TNFAIP6 is also essential in injury and inflammation. Further, TNFA1P6 is strongly linked to arsenic poisoning, evident from the results of NetworkAnalyst, a comprehensive and interactive platform for gene expression profiling via network visual analytics. As a result, the TNFAIP6 gene was ultimately chosen as a candidate DEG and subsequently employed for in silico structural characterization studies. The tertiary structure of TNFAIP6 was modelled using the ROBETTA server, followed by validation with SAVES and ProSA webserver. Additionally, structural dynamic studies, including molecular dynamics simulation (MDS) and essential dynamics analysis, including principal component (PC) based free energy landscape (FEL) analysis, was used for checking the stability of TNFAIP6 models. The dynamics result established the structural rigidity of modelled TNFAIP6 through RMSD, RMSF and RoG results. The FEL analysis revealed the restricted conformational flexibility of TNFAIP6 by displaying a single minimum energy basin in the contour plot. The comprehensive computational analysis established that TNFAIP6 could serve as a viable biomarker to assess the severity of pulmonary TB.

A computational perspective towards the identification of promising lead molecules against 6-hydroxy-methyl dihydropterin pyrophosphokinase (HPPK) from Acinetobacter baumannii

Acinetobacter baumannii is an ESKAPE pathogen that causes endocarditis, pneumonia, blood infections, urinary tract infections, and several other illnesses. In addition, it is mainly responsible for nosocomial infection-related mortality. Gram-negative A. baumannii bacterium (AYE Strain) has high MDR and XDR levels. Due to its function in synthesizing purines and amino acids, folic acid is a significant molecule necessary for the growth of bacteria. The metabolic pathway of folate production is therefore a potential therapeutic target to inhibit bacterial growth. In the current study, the three-dimensional model of 6-Hydroxy-methyl dihydropterinpyrophosphokinase (HPPK) was predicted and subsequently processed through a virtual high throughput screening (vHTS) against compounds from Enamine HTSC library, that could bind to its active site. Three lead candidates (Z73322064, Z354558542, and Z906123504) and a control molecule (7,8 dihydro-7,7-dimethyl-6-hydroxymethlypterin; Accession Number: DB02278) were identified using several screening criteria namely estimated binding affinity, estimated inhibition constant, drug-like properties, ADME properties, mode of binding, and interaction patterns of the screened compounds. The physiological behavior of ligand binding on the HPPK enzyme was then studied using molecular dynamics simulations of apo and ligand bound complexes. This study proposed the following three molecules: Z73322064, Z354338542, and Z906123504 as promising lead candidates against the substrate-binding site of the HPPK enzyme from A. baumannii using global, essential dynamics studies along with MM/PBSA based binding free energy analysis. Communicated by Ramaswamy H. Sarma

Schooling of light reflecting fish.

One of the hallmarks of the collective movement of large schools of pelagic fish are waves of shimmering flashes that propagate across the school, usually following an attack by a predator. Such flashes arise when sunlight is reflected off the specular (mirror-like) skin that characterizes many pelagic fishes, where it is otherwise thought to offer a means for camouflage in open waters. While it has been suggested that these 'shimmering waves' are a visual manifestation of the synchronized escape response of the fish, the phenomenon has been regarded only as an artifact of esthetic curiosity. In this study we apply agent-based simulations and deep learning techniques to show that, in fact, shimmering waves contain information on the behavioral dynamics of the school. Our analyses are based on a model that combines basic rules of collective motion and the propagation of light beams in the ocean, as they hit and reflect off the moving fish. We use the resulting reflection patterns to infer the essential dynamics and inter-individual interactions which are necessary to generate shimmering waves. Moreover, we show that light flashes observed by the school members themselves may extend the range at which information can be communicated across the school. Assuming that fish pay heed to this information, for example by entering an apprehensive state of reduced response-time, our analysis suggests that it can speed up the propagation of information across the school. Further still, we use an artificial neural network to show that light flashes are, on their own, indicative of the state and dynamics of the school, and are sufficient to infer the direction of attack and the shape of the school with high accuracy.

'Drug decriminalization' in Canada: a plea for a nuanced, evidence-informed, and realistic approach towards improved health outcomes.

Canada has been home to an exceptional public health crisis of illicit (mainly synthetic opioid but increasingly other) drug deaths for numerous years. Despite extensive health-oriented interventions, this crisis remains acute and devastating. Earlier this year, the federal and British Columbia governments implemented the 'drug decriminalization initiative' (DDI) which exempts drug users from criminal enforcement of small amounts of personal drug possession; this has been promoted as a 'solution' for the drug-death crisis. This commentary examines and discusses that while decriminalization of drug use is fundamentally required to align drug policy with principles of health, evidence-based consideration of the drug death crisis' essential dynamics as well as enforcement realities substantially curtail the DDI's plausible prospects to significantly reduce the drug-fatality toll. Decriminalization efforts elsewhere ('Portugal model') are set in distinct ecologies and parameters, and not readily transferable to Canadian contexts. While the DDI has symbolic value for policy reform, for decriminalization to provide tangible benefits for drug users and public health, its scope and design require fundamental reconsideration and improvement. Meanwhile, given the principally causal role of toxic drug supply for the current drug-death epidemic, more effective health and life-saving measures-specifically involving 'safer supply'-need to be urgently implemented to alleviate this crisis.

Computational study of SENP1 in cancer by novel natural compounds and ZINC database screening.

Introduction: Sentrin-specific protease 1 (SENP1) is a protein whose main function is deSUMOylation. SENP1 inhibits apoptosis, and increases angiogenesis, estrogen and androgen receptor transcription and c-Jun transcription factor, proliferation, growth, cell migration, and invasion of cancer. The in vivo and in vitro studies also demonstrated which natural compounds, especially phytochemicals, minerals, and vitamins, prevent cancer. More than 3,000 plant species have been reported in modern medicine. Natural compounds have many anti-cancerous andanti-turmeric properties such as antioxidative, antiangiogenic, antiproliferative, and pro-apoptotic properties. Methods: In this study, we investigated the interaction of some natural compounds with SENP1 to inhibit its activity. We also screened the ZINC database including natural compounds. Molecular docking was performed, and toxicity of compounds was determined; then, molecular dynamics simulation (MDS) and essential dynamics (ED) were performed on natural compounds with higher free binding energies and minimal side effects. By searching in a large library, virtual screening of the ZINC database was performed using LibDock and CDOCKER, and the final top 20 compounds were allowed for docking against SENP1. According to the docking study, the top three leading molecules were selected and further analyzed by MDS and ED. Results: The results suggest that resveratrol (from the selected compounds) and ZINC33916875 (from the ZINC database) could be more promising SENP1 inhibitory ligands. Discussion: Because these compounds can inhibit SENP1 activity, then they can be novel candidates for cancer treatment. However, wet laboratory experiments are needed to validate their efficacy as SENP1 inhibitors.

EMBED: Essential MicroBiomE Dynamics, a dimensionality reduction approach for longitudinal microbiome studies

Dimensionality reduction offers unique insights into high-dimensional microbiome dynamics by leveraging collective abundance fluctuations of multiple bacteria driven by similar ecological perturbations. However, methods providing lower-dimensional representations of microbiome dynamics both at the community and individual taxa levels are not currently available. To that end, we present EMBED: Essential MicroBiomE Dynamics, a probabilistic nonlinear tensor factorization approach. Like normal mode analysis in structural biophysics, EMBED infers ecological normal modes (ECNs), which represent the unique orthogonal modes capturing the collective behavior of microbial communities. Using multiple real and synthetic datasets, we show that a very small number of ECNs can accurately approximate microbiome dynamics. Inferred ECNs reflect specific ecological behaviors, providing natural templates along which the dynamics of individual bacteria may be partitioned. Moreover, the multi-subject treatment in EMBED systematically identifies subject-specific and universal abundance dynamics that are not detected by traditional approaches. Collectively, these results highlight the utility of EMBED as a versatile dimensionality reduction tool for studies of microbiome dynamics.

Integrative dynamic structural biology unveils conformers essential for the oligomerization of a large GTPase.

Guanylate binding proteins (GBPs) are soluble dynamin-like proteins that undergo a conformational transition for GTP-controlled oligomerization and disrupt membranes of intracellular parasites to exert their function as part of the innate immune system of mammalian cells. We apply neutron spin echo, X-ray scattering, fluorescence, and EPR spectroscopy as techniques for integrative dynamic structural biology to study the structural basis and mechanism of conformational transitions in the human GBP1 (hGBP1). We mapped hGBP1's essential dynamics from nanoseconds to milliseconds by motional spectra of sub-domains. We find a GTP-independent flexibility of the C-terminal effector domain in the µs-regime and resolve structures of two distinct conformers essential for an opening of hGBP1 like a pocket knife and for oligomerization. Our results on hGBP1's conformational heterogeneity and dynamics (intrinsic flexibility) deepen our molecular understanding relevant for its reversible oligomerization, GTP-triggered association of the GTPase-domains and assembly-dependent GTP-hydrolysis.

What Contributes to the Measured Chiral Optical Response of the Glutathione-Protected Au25 Nanocluster?

The water-soluble glutathione-protected [Au25(GSH)18]-1 nanocluster was investigated by integrating several methodologies such as molecular dynamics simulations, essential dynamics analysis, and state-of-the-art time-dependent density functional theory calculations. Fundamental aspects such as conformational, weak interactions and solvent effects, especially hydrogen-bonds, were included and found to play a fundamental role in assessing the optical response of this system. Our analysis demonstrated not only that the electronic circular dichroism is extremely sensitive to the solvent presence but also that the solvent itself plays an active role in the optical activity of such system, forming a chiral solvation shell around the cluster. Our work demonstrates a successful strategy to investigate in detail chiral interfaces between metal nanoclusters and their environments, applicable, e.g., to chiral electronic interactions between clusters and biomolecules.

MDiGest: A Python package for describing allostery from molecular dynamics simulations.

Many biological processes are regulated by allosteric mechanisms that communicate with distant sites in the protein responsible for functionality. The binding of a small molecule at an allosteric site typically induces conformational changes that propagate through the protein along allosteric pathways regulating enzymatic activity. Elucidating those communication pathways from allosteric sites to orthosteric sites is, therefore, essential to gain insights into biochemical processes. Targeting the allosteric pathways by mutagenesis can allow the engineering of proteins with desired functions. Furthermore, binding small molecule modulators along the allosteric pathways is a viable approach to target reactions using allosteric inhibitors/activators with temporal and spatial selectivity. Methods based on network theory can elucidate protein communication networks through the analysis of pairwise correlations observed in molecular dynamics (MD) simulations using molecular descriptors that serve as proxies for allosteric information. Typically, single atomic descriptors such as α-carbon displacements are used as proxies for allosteric information. Therefore, allosteric networks are based on correlations revealed by that descriptor. Here, we introduce a Python software package that provides a comprehensive toolkit for studying allostery from MD simulations of biochemical systems. MDiGest offers the ability to describe protein dynamics by combining different approaches, such as correlations of atomic displacements or dihedral angles, as well as a novel approach based on the correlation of Kabsch-Sander electrostatic couplings. MDiGest allows for comparisons of networks and community structures that capture physical information relevant to allostery. Multiple complementary tools for studying essential dynamics include principal component analysis, root mean square fluctuation, as well as secondary structure-based analyses.

Laboratory experiments on coastal current separation and vortex formation in front of Cabo Corrientes and Bahía de Banderas, México

A series of rotating tank experiments are carried out to study the behaviour of a coastal current when encountering a cape and an adjacent bay. The experiments are designed to represent some essential dynamics of the interaction of the Mexican Coastal Current with Cabo Corrientes and Bahía de Banderas on the western coast of Mexico. The experiments show that the current surrounds the cape and accumulates negative relative vorticity in front of the downstream bay. For sufficiently strong currents, the flow separates from the coast and the accumulated vorticity develops an anticyclonic eddy that drifts away from the generation region. Afterwards, the current is re-established and the process is repeated. The periodicity is about 16 periods of the rotating system. The current separation and vortex formation depend on the ratio between the inertial length of the current (defined as the characteristic speed over the Coriolis parameter) and the coastal radius of curvature associated with the cape. Another mechanism for vortex formation is verified when the coastal current is switched off. The current separation and eventual vortex shedding are further studied in additional experiments without the cape or the bay, thus increasing the effective coastal radius of curvature. It is verified that the absence of any of these two features inhibits the current separation. The results illustrate the dynamical conditions that may help to explain the emergence of mesoscale oceanic eddies observed on the lee side of Cabo Corrientes, especially during the summer. This assertion is explored by examining historical data from the GLORYS12V1 reanalysis, in which current separation and eddy-shedding events are found and compared with the experimental observations.

Structure-guided identification of potent inhibitors of ROS1 kinase for therapeutic development against non-small cell lung cancer

Proto-oncogene tyrosine-protein kinase ROS (ROS1) is a member of the sevenless receptor, which affects epithelial cell differentiation and is highly expressed in a variety of tumor cells. The elevated expression and dysfunction of ROS1 have been involved in various malignancies, such as non-small cell lung cancer (NSCLC), stomach cancer, ovarian, breast cancer, cholangiocarcinoma, colorectal cancer, adenosarcoma, oesophageal cancer, etc. ROS1 has been postulated as a potential drug target in anticancer therapeutics. In this study, we carried out a virtual screening of phytochemicals against ROS1 to identify its potential inhibitors. The virtual screening process was performed on the ROS1 structure, where two phytochemicals, Helioscopinolide C and Taiwanin C, were identified. These compounds resulted from filters like Lipinski rule of five, PAINS filter, binding affinities values, and all-atom molecular dynamics (MD) simulations followed by principal component analysis (PCA) and essential dynamics. The findings of this study highlight the role of ROS1 in multiple physiological candidates and its therapeutic targeting using phytochemicals. This study suggests Helioscopinolide C and Taiwanin C as potential compounds for therapeutic development targeting ROS1-associated non-small cell lung cancer for clinical applications. Further in vitro and in vivo experiments are required to validate these findings. Communicated by Ramaswamy H. Sarma

Effect of osmolytes on the EcoRI endonuclease: Insights into hydration and protein dynamics from molecular dynamics simulations

Osmolytes play an important role in cellular physiology by modulating the properties of proteins, including their molecular specificity. EcoRI is a model restriction enzyme whose specificity to DNA is altered in the presence of osmolytes. Here, we investigate the effect of two different osmolytes, glycerol and DMSO, on the dynamics and hydration of the EcoRI enzyme using molecular dynamics simulations. Our results show that the osmolytes, alter the essential dynamics of EcoRI. Particularly, we observe that the dynamics of the arm region of EcoRI which is involved in DNA binding is significantly altered. In addition, conformational free energy analyses reveals that the osmolytes bring about a change in the landscape similar to that of EcoRI bound to cognate DNA. We further observe that the hydration of the enzyme for each of the osmolyte is different, indicating that the mechanism of action of each of these osmolytes could be different. Further analyses of interfacial water dynamics using rotational autocorrelation function reveals that while the protein surface contributes to a slower tumbling motion of water, osmolytes, additionally contribute to the slowing of the angular motion of the water molecules. Entropy analysis also corroborates with this finding. We also find that the slowed rotational motion of interfacial waters in the presence of osmolytes contributes to a slowed relaxation of the hydrogen bonds between the interfacial waters and the functionally important residues in the protein. Taken together, our results show that osmolytes alter the dynamics of the protein by altering the dynamics of water. This altered dynamics, mediated by the changes in the water dynamics and hydrogen bonds with functionally important residues, may contribute to the altered specificity of EcoRI in the presence of osmolytes.

The deleterious variants of N-acetylgalactosamine-6-sulfatase (GalN6S) enzyme trigger Morquio a syndrome by disrupting protein foldings

Lysosomal enzymes degrade cellular macromolecules, while their inactivation causes human hereditary metabolic disorders. Mucopolysaccharidosis IVA (MPS IVA; Moquio A syndrome) is one of the lysosomal storage disorders caused by a defective Galactosamine-6-sulfatase (GalN6S) enzyme. In several populations, disease incidence is elevated due to missense mutations brought on by non-synonymous allelic variation in the GalN6S enzyme. Here, we studied the effect of non-synonymous single nucleotide polymorphism (nsSNPs) on the structural dynamics of the GalN6S enzyme and its binding with N-acetylgalactosamine (GalNAc) using all-atom molecular dynamics simulation and an essential dynamics approach. Consequently, in this study, we have identified three functionally disruptive mutations in domain-I and domain-II, that is, S80L, R90W, and S162F, which presumably contribute to post-translational modifications. The study delineated that both domains work cooperatively, and alteration in domain II (S80L, R90W) leads to conformational changes in the catalytic site in domain-I, while mutation S162F mainly provokes higher residual flexibility of domain II. These results show that these mutations impair the hydrophobic core, implying that Morquio A syndrome is caused by misfolding of the GalN6S enzyme. The results also show the instability of the GalN6S-GalNAc complex upon substitution. Overall, the structural dynamics resulting from point mutations give the molecular rationale for Moquio A syndrome and, more importantly, the Mucopolysaccharidoses (MPS) family of diseases, re-establishing MPS IVA as a protein-folding disease. Communicated by Ramaswamy H. Sarma

Identification of potential inhibitors of tropomyosin receptor kinase B targeting CNS-related disorders and cancers

Tropomyosin receptor kinase B (TrkB), also known as neurotrophic tyrosine kinase receptor type 2 (NTRK2), is a protein that belongs to the family of receptor tyrosine kinases (RTKs). NTRK2 plays a crucial role in regulating the development and maturation of the central nervous system (CNS) and peripheral nervous system (PNS). Elevated TrkB expression levels observed in different pathological conditions make it a potential target for therapeutic interventions against neurological disorders, including depression, anxiety, Alzheimer’s disease, Parkinson’s disease, and certain types of cancer. Targeting TrkB using small molecule inhibitors is a promising strategy for the treatment of a variety of neurological disorders. In this research, a systematic virtual screening was carried out on phytoconstituents found in the IMPPAT library to identify compounds potentially inhibiting TrkB. The retrieved compounds from the IMPPAT library were first filtered using Lipinski’s rule of five. The compounds were then sorted based on their docking score and ligand efficiency. In addition, PAINS, ADMET, and PASS evaluations were carried out for selecting drug-like compounds. Finally, in interaction analysis, we found two phytoconstituents, Wedelolactone and 3,8-dihydroxy-1-methylanthraquinone-2-carboxylic acid (DMCA), which possessed considerable docking scores and specificity on the TrkB ATP-binding pocket. The selected compounds were further assessed employing molecular dynamics (MD) simulations and essential dynamics. The results revealed that the elucidated compounds bind well with the TrkB binding pocket and lead to fewer conformations fluctuations. This study highlighted using phytoconstituents, Wedelolactone and DMCA as starting leads in developing novel TrkB inhibitors. Communicated by Ramaswamy H. Sarma