Network pharmacological study and in vitro studies validation-Molecular dynamics simulation of Cistanche deserticola in promoting periodontitis and bone remodeling

OBJECTIVE

Periodontitis is a chronic infectious disease, characterized by loss of alveolar bone and supporting tissues. Cistanche deserticola(Cd), a local medicinal herb in Xinjiang, possesses favorable biological characteristics and potential applications. Our aim is to investigate the remodeling properties of Cd extract and elucidate the specific mechanisms underlying its therapeutic effects on periodontitis, by employing a combination of basic experimental and network pharmacology approaches.

METHODS

Firstly, UHPLC-QTOF-MS analysis was conducted on Cd extract to identify its main components, with several compounds were identified by standard. Subsequently, in vitro studies were performed using the Cd extract on MC3T3-E1 cells. Cell proliferation viability was assessed using CCK-8 and apoptosis assays, while ALP and ARS staining and quantitative experiments, qRT-PCR, and Western blot assays were employed to evaluate the osteogenic differentiation capability. Network pharmacology analysis was then carried out using the identified compounds to establish a database of Cd components and targets, along with a database of periodontitis. The intersection of these databases revealed the network relationship between Cd components-mapped genes-signaling pathways. KEGG/GO pathway analysis of the targets was performed to filter potential enriched pathways. PPI/CytoHubba protein interaction network analysis was utilized to identify hub genes. Molecular docking and molecular dynamics simulations were employed to analyze the docking and interaction between core gene and Cd components.

RESULTS

We detected 38 major components in the Cd extract, with Echinacoside, Acteoside, Tubuloside A, and Cistanoside A undergoing standard substance verification. In vitro studies indicated that the Cd, at concentrations below 100 μg/ mL, did not affect cell proliferation and inhibited apoptosis. Osteogenesis assays demonstrated that Cd at concentrations of 1 μg/ mL, 10 μg/ mL, and 100 μg/ mL significantly promoted the osteogenic differentiation ability of MC3T3-E1 cells. It also notably upregulated the mRNA and protein levels of Alp, Bmp2, Runx2, and Opn, and the optimal concentration was 10 μg/mL. Network pharmacology results revealed the network relationship between Cd's components, crossed targets and signaling pathways. Combined with KEGG/GO pathway analysis and PPI/CytoHubba protein interaction network analysis. The key pathway and hub genes of Cd regulating periodontitis are both related to hypoxia pathway and HIF-1α. Molecular docking results showed a strong binding affinity between Cd compounds and hub genes, and molecular dynamics simulation results indicated the stability of the complexes formed between HIF-1α and several Cd compounds.

CONCLUSION

Cistanche deserticola exhibits a notable capacity to promote bone regeneration, and its mechanism of action in regulating periodontitis is associated with the hypoxia signaling pathway. HIF-1α may serve as a potential core gene. Future research will focus on exploring the mechanism of Cd in intervene periodontitis and promoting bone remodeling in hypoxic environment.

Structural mechanism underlying variations in DNA binding by the androgen receptor

The androgen receptor (AR) is a steroid activated transcription factor which recognizes DNA motifs resembling inverted repeats of a conserved 5'-AGAACA-3'-like hexanucleotides separated by a three-nucleotide spacer from a similar, but less conserved hexanucleotide. Here, we report the structures of the human AR DNA binding domain (DBD) bound to two natural AREs (C3 and MTV) in head-to-head dimer conformations, diffracting at 2.05 Å and 2.25 Å, respectively. These structures help to explain the impact of androgen insensitivity mutations on the structure integrity, DNA binding and DBD dimerization. The binding affinity of the AR DBD to different DNA motifs were measured by the BioLayer Interferometry (BLI) and further validated by Molecular Dynamics (MD) simulations. This shows that the high binding affinity of the first DBD to the upstream 5'-AGAACA-3' motif induces the cooperative binding of the second DBD to the second hexanucleotide. Our data indicate identical interaction of the DBDs to the upstream hexanucleotides, while forming an induced closer contact of the second DBD on the non-canonical hexanucleotides. The variation in binding between the DBD monomers are the result of differences in DNA occupancy, protein-protein interactions, DNA binding affinity, and DNA binding energy profiles. We propose this has functional consequences.

Role of phosphatidic acid lipids on plasma membrane association of the Ebola virus matrix protein VP40

The Ebola virus matrix protein VP40 is responsible for the formation of the viral matrix by localizing at the inner leaflet of the human plasma membrane (PM). Various lipid types, including PI(4,5)P2 (i.e. PIP2) and phosphatidylserine (PS), play active roles in this process. Specifically, the negatively charged headgroups of both PIP2 and PS interact with the basic residues of VP40 and stabilize it at the membrane surface, allowing for eventual egress. Phosphatidic acid (PA), resulting from the enzyme phospholipase D (PLD), is also known to play an active role in viral development. In this work, we performed a biophysical and computational analysis to investigate the effects of the presence of PA on the membrane localization and association of VP40. We used coarse-grained molecular dynamics simulations to quantify VP40 hexamer interactions with the inner leaflet of the PM. Analysis of the local distribution of lipids shows enhanced lipid clustering when PA is abundant in the membrane. We observed that PA lipids have a similar role to that of PS lipids in VP40 association due to the geometry and charge. Complementary experiments performed in cell culture demonstrate competition between VP40 and a canonical PA-binding protein for the PM. Also, inhibition of PA synthesis reduced the detectable budding of virus-like particles. These computational and experimental results provide new insights into the early stages of Ebola virus budding and the role that PA lipids have on the VP40-PM association.

In-vitro and computational analysis of Urolithin-A for anti-inflammatory activity on Cyclooxygenase 2 (COX-2)

Cyclooxygenase 2 (COX-2) participates in the inflammation process by converting arachidonic acid into prostaglandin G2 which increases inflammation, pain and fever. COX-2 has an active site and a heme pocket and blocking these sites stops the inflammation. Urolithin A is metabolite of ellagitannin produced from humans and animals gut microbes. In the current study, Urolithin A showed good pharmacokinetic properties. Molecular docking of the complex of Urolithin A and COX-2 revealed the ligand affinity of -7.97 kcal/mol with the ligand binding sites at TYR355, PHE518, ILE517 and GLN192 with the 4-H bonds at a distance of 2.8 Å, 2.3 Å, 2.5 Å and 1.9 Å. The RMSD plot for Urolithin A and COX-2 complex was observed to be constant throughout the duration of dynamics. A total of 3 pair of hydrogen bonds was largely observed on average of 3 simulation positions for dynamics duration of 500 ns. The MMPBSA analysis showed that active site amino acids had a binding energy of -22.0368 kJ/mol indicating that throughout the simulation the protein of target was bounded by Urolithin A. In-silico results were validated by biological assays. Urolithin A strongly revealed to exhibit anti-inflammatory effect on COX-2 with an IC50 value of 44.04 µg/mL. The anti-inflammatory capability was also depicted through reduction of protein denaturation that showed 37.6 ± 0.1 % and 43.2 ± 0.07 % reduction of protein denaturation for BSA and egg albumin respectively at 500 µg/mL. The present study, suggests Urolithin A to be an effective anti-inflammatory compound for therapeutic use.

Exome sequence analysis of rare frequency variants in Late-Onset Alzheimer Disease

Alzheimer disease (AD) is a leading cause of dementia in elderly patients who continue to live between 3 and 11 years of diagnosis. A steep rise in AD incidents is observed in the elderly population in East-Asian countries. The disease progresses through several changes, including memory loss, behavioural issues, and cognitive impairment. The etiology of AD is hard to determine because of its complex nature. The whole exome sequences of late-onset AD (LOAD) patients of Korean origin are investigated to identify rare genetic variants that may influence the complex disorder. Computational annotation was performed to assess the function of candidate variants in LOAD. The in silico pathogenicity prediction tools such as SIFT, Polyphen-2, Mutation Taster, CADD, LRT, PROVEAN, DANN, VEST3, fathmm-MKL, GERP +  + , SiPhy, phastCons, and phyloP identified around 17 genes harbouring deleterious variants. The variants in the ALDH3A2 and RAD54B genes were pathogenic, while in 15 other genes were predicted to be variants of unknown significance. These variants can be potential risk candidates contributing to AD. In silico computational techniques such as molecular docking, molecular dynamic simulation and steered molecular dynamics were carried out to understand the structural insights of RAD54B with ATP. The simulation of mutant (T459N) RAD54B with ATP revealed reduced binding strength of ATP at its binding site. In addition, lower binding free energy was observed when compared to the wild-type RAD54B. Our study shows that the identified uncommon variants are linked to AD and could be probable predisposing genetic factors of LOAD.

Application of molecular dynamics simulation in the field of food enzymes: improving the thermal-stability and catalytic ability

Enzymes can produce high-quality food with low pollution, high function, high acceptability, and medical aid. However, most enzymes, in their native form, do not meet the industrial requirements. Sequence-based and structure-based methods are the two main strategies used for enzyme modification. Molecular Dynamics (MD) simulation is a sufficiently comprehensive technology, from a molecular perspective, which has been widely used for structure information analysis and enzyme modification. In this review, we summarize the progress and development of MD simulation, particularly for software, force fields, and a standard procedure. Subsequently, we review the application of MD simulation in various food enzymes for thermostability and catalytic improvement was reviewed in depth. Finally, the limitations and prospects of MD simulation in food enzyme modification research are discussed. This review highlights the significance of MD simulation and its prospects in food enzyme modification.

Micellar Composition Affects Lipid Accretion Kinetics in Molecular Dynamics Simulations: Support for Lipid Network Reproduction

Mixed lipid micelles were proposed to facilitate life through their documented growth dynamics and catalytic properties. Our previous research predicted that micellar self-reproduction involves catalyzed accretion of lipid molecules by the residing lipids, leading to compositional homeostasis. Here, we employ atomistic Molecular Dynamics simulations, beginning with 54 lipid monomers, tracking an entire course of micellar accretion. This was done to examine the self-assembly of variegated lipid clusters, allowing us to measure entry and exit rates of monomeric lipids into pre-micelles with different compositions and sizes. We observe considerable rate-modifications that depend on the assembly composition and scrutinize the underlying mechanisms as well as the energy contributions. Lastly, we describe the measured potential for compositional homeostasis in our simulated mixed micelles. This affirms the basis for micellar self-reproduction, with implications for the study of the origin of life.

Non-enzymatic glycation of human serum albumin modulates its binding efficacy towards bioactive flavonoid chrysin: A detailed study using multi-spectroscopic and computational methods

The non-enzymatic glycation of plasma proteins by reducing sugars have important consequences on the conformational and functional properties of protein. The formation of advanced glycation end products (AGEs) is responsible for cell death and other pathological conditions. We have synthesized the glycated human serum albumin (gHSA) and characterized the same by using differential spectroscopic measurements. The aim of the present study is to determine the effect of glycation on the binding of human serum albumin (HSA) with bioactive flavonoid chrysin, which possesses anti-cancer, anti-inflammatory and anti-oxidant activities. The interaction of chrysin with HSA and gHSA was studied using multi-spectroscopic, molecular docking and molecular dynamics (MD) simulation techniques. Chrysin quenched the intrinsic fluorescence of both HSA and gHSA by static quenching mechanism. The value of the binding constant (K_b) for the interaction of HSA-chrysin complex (4.779 ± 0.623 × 10⁵ M^-1 at 300 K) was found to be higher than that of gHSA-chrysin complex (2.206 ± 0.234 × 10⁵ M^-1 at 300 K). Hence, non-enzymatic glycation of HSA significantly reduced its binding affinity towards chrysin. The % α-helicity of HSA was found to get enhanced upon binding with chrysin, and minimal changes were observed for the gHSA-chrysin complex. Site marker probe studies indicated that chrysin binds to subdomain IIA and IIIA of both HSA and gHSA. The results from molecular docking and MD simulation studies correlated well with the experimental findings. Electrostatic interactions followed by hydrogen bonding and hydrophobic interactions played major roles in the binding process. These observations may have some useful insights into the field of pharmaceutics.

Substitution impact of highly conserved arginine residue at position 75 in GJB1 gene in association with X-linked Charcot-Marie-tooth disease: A computational study

X-linked Charcot-Marie-Tooth type 1 X (CMTX1) disease is a subtype of Charcot-Marie-Tooth (CMT), which is mainly caused by mutations in the GJB1 gene. It is also known as connexin 32 (Cx32) that leads to Schwann cell abnormalities and peripheral neuropathy. CMTX1 is considered as the second most common form of CMT disease. The aim of this study is to computationally predict the potential impact of different single amino acid substitutions at position 75 of Cx32, from arginine (R) to proline (P), glutamine (Q) and tryptophan (W). This position is known to be highly conserved among the family of connexin. To understand the structural and functional changes due to these single amino acid substitutions, we employed a homology-modeling technique to build the three-dimensional structure models for the native and mutant proteins. The protein structures were further embedded into a POPC lipid bilayer, inserted into a water box, and subjected to molecular dynamics simulation for 50 ns. Our results show that the mutants R75P, R75Q and R75W display variable structural conformation and dynamic behavior compared to the native protein. Our data proves useful in predicting the potential pathogenicity of the mutant proteins and is expected to serve as a platform for drug discovery for patients with CMT.

Optimal Magnetic Field for Crossing Super-Para-Magnetic Nanoparticles through the Brain Blood Barrier: A Computational Approach

This paper scrutinizes the magnetic field effect to deliver the superparamagnetic nanoparticles (SPMNs) through the Blood Brain Barrier (BBB). Herein we study the interaction between the nanoparticle (NP) and BBB membrane using Molecular Dynamic (MD) techniques. The MD model is used to enhance our understanding of the dynamic behavior of SPMNs crossing the endothelial cells in the presence of a gradient magnetic field. Actuation of NPs under weak magnetic field offers the great advantage of a non-invasive drug delivery without the risk of causing injury to the brain. Furthermore, a weak magnetic portable stimulator can be developed using low complexity prototyping techniques. Based on MD simulation results in this paper, SPMNs can cross the cell membrane while experiencing very weak mechanical forces in the range of pN. This study also derives guidelines for the design of the SPMNs dedicated to crossing the BBB using external magnetic fields.

An anti-DNA antibody prefers damaged dsDNA over native

DNA-protein interactions, including DNA-antibody complexes, have both fundamental and practical significance. In particular, antibodies against double-stranded DNA play an important role in the pathogenesis of autoimmune diseases. Elucidation of structural mechanisms of an antigen recognition and interaction of anti-DNA antibodies provides a basis for understanding the role of DNA-containing immune complexes in human pathologies and for new treatments. Here we used Molecular Dynamic simulations of bimolecular complexes of a segment of dsDNA with a monoclonal anti-DNA antibody's Fab-fragment to obtain detailed structural and physical characteristics of the dynamic intermolecular interactions. Using a computationally modified crystal structure of a Fab-DNA complex (PDB: 3VW3), we studied in silico equilibrium Molecular Dynamics of the Fab-fragment associated with two homologous dsDNA fragments, containing or not containing dimerized thymine, a product of DNA photodamage. The Fab-fragment interactions with the thymine dimer-containing DNA was thermodynamically more stable than with the native DNA. The amino acid residues constituting a paratope and the complementary nucleotide epitopes for both Fab-DNA constructs were identified. Stacking and electrostatic interactions were shown to play the main role in the antibody-dsDNA contacts, while hydrogen bonds were less significant. The aggregate of data show that the chemically modified dsDNA (containing a covalent thymine dimer) has a higher affinity toward the antibody and forms a stronger immune complex. These findings provide a mechanistic insight into formation and properties of the pathogenic anti-DNA antibodies in autoimmune diseases, such as systemic lupus erythematosus, associated with skin photosensibilization and DNA photodamage.

The role of Ca(2+) ions in the complex assembling of protein Z and Z-dependent protease inhibitor: A structure and dynamics investigation

We investigated the solution structure and dynamics of the human anti-coagulation protein Z (PZ) in the complex with protein Zdependent protease inhibitor (ZPI) to order to understand key structural changes in the presence and absence of Ca(2+). Structural features of the complete complex of PZ-ZPI are poorly understood due to lack of complete atomic model of the PZ-ZPI complex. We have constructed a model of the complete PZ-ZPI complex and molecular dynamics (MD) simulation of the solvated PZ-ZPI complex with and without Ca(2+) was achieved for 100ns. It is consider that the Ω-loop of GLA domains interacts with negatively charged biological membranes in the presence of Ca(2+) ions. The PZ exerts its role as cofactor in a similar way. However, we used solvent-equilibrated dynamics to show structural features of the PZ-ZPI complex in the presence and the absence of Ca(2+)ions. We observed that the distance between the interacting sites of the ZPI with the PZ and the GLA domain decreases in the presence of Ca(2+) ions. Further, we postulated that the calculated distance between the dominant plane of the Ca(2+) ions and Ser196 of the pseudo-catalytic triad of the PZ is similar to the equivalent distance of FXa. This suggests that the central role of the PZ in the blood coagulation may be to align the inhibitory site of the ZPI with the active site of the FXa, which is depends on the interaction of the calcium bound GLA domain of the PZ with the active membrane.