Computational modelling and optimization studies of electropentamer for molecular imprinting of DJ-1.

Parkinson's disease (PD) is the most prevalent type of incurable movement disorder. Recent research findings propose that the familial PD-associated molecule DJ-1 exists in cerebrospinal fluid (CSF) and that its levels may be altered as Parkinson's disease advances. By using a molecularly imprinted polymer (MIP) as an artificial receptor, it becomes possible to create a functional MIP with predetermined selectivity for various templates, particularly for the DJ-1 biomarker associated with Parkinson's disease. It mostly depends on molecular recognition via interactions between functional monomers and template molecules. So, the computational methods for the appropriate choice of functional monomers for creating molecular imprinting electropolymers (MIEPs) with particular recognition for the detection of DJ-1, a pivotal biomarker involved in PD, are undertaken in this study. Here, molecular docking, molecular dynamics simulations (MD), molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) methods, and quantum mechanical calculation have been applied to investigate the intermolecular interaction between DJ-1 and several functional electropentamers, viz., polypyrrole (PPy), poly(3,4-ethylenedioxythiophene) (PEDOT), poly(o-aminophenol) (POAP), and polythiophene (PTS). In this context, the electropentamers were selected to mimic the imprinted electropolymer system. We analyzed the most stable configurations of the formed complexes involving DJ-1 and electropentamers as a model system for MIEPs. Among these, PEDOT exhibited a more uniform arrangement around DJ-1, engaging in numerous van der Waals, H-bond, electrostatic, and hydrophobic interactions. Hence, it can be regarded as a preferable choice for synthesizing a MIP for DJ-1 recognition. Thus, it will aid in selecting a suitable functional monomer, which is of greater significance in the design and development of selective DJ-1/MIP sensors.

Mechanistic insights into the conformational switch in profilin-1 subject to collective effects of mutation and histidine tautomerism.

Mutations and histidine (His) tautomerism in profilin-1 (PFN1) are associated with amyotrophic lateral sclerosis (ALS). The conformational changes in PFN1 caused by the collective effects of G117V mutation and His tautomeric isomers εε, εδ, δε, and δδ were clarified using molecular dynamics (MD) simulations. The predominant structural variations were seen in α-helices, ß-sheets, turns, and coils and the His tautomer's unique degree of disruption was seen in these conformations. The content of α-helices was 23.2 % in the εε and δδ isomers, but the observed α-helices in the isomers εδ and δε were 20.3 % and 21.7 % respectively. The percentage of ß-sheet was found to be higher (34.1) in the εε isomer than in the εδ, δε, and δδ isomers, and the values were 30.4, 29.7, and 31.9, respectively. Intermolecular water dynamics analysis discloses that His 133 can form an intramolecular H-bond interaction (Nα-H---Nδ), confirming the experimental observations in the simulations of εε, δε, and δδ isomers of G117V PFN1 mutant. It was concluded that these solvent molecules are crucial for aggregation and must be considered in future research on the PFN1 associated with ALS. Overall, the study offers a thorough microscopic understanding of the pathogenic mechanisms behind conformational changes that cause aggregation illnesses like ALS.

Histidine tautomerism dependent conformational transitions driven aggregation of profilin-1: Implications in amyotrophic lateral sclerosis.

Aggregation of profilin-1 (PFN1) causes a fatal neurodegenerative disease, familial amyotrophic lateral sclerosis (fALS). Histidine (His) tautomerism has been linked to the formation of fibril aggregation causing neurodegenerative disease. Characterization of intermediate species that form during aggregation is crucial, however, this has proven very challenging for experimentalists due to their transient nature. Hence, molecular dynamics (MD) simulations have been performed on the His tautomeric isomers εε, εδ, δε, and δδ of PFN1 to explain the structural changes and to correlate them with its aggregation propensity. MD simulations show that His133 presumably plays a major role in the aggregation of PFN1 upon His tautomerism compared to His119. Further, the formation of a new 310-helix is observed in εε and δε but 310-helix is not observed in δδ and εδ isomers. In addition, our findings unveil that ß-sheet dominating conformations are observed in His119(δ)-His133(δ) δδ isomer of PFN1 with significant antiparallel ß-sheets between residues T15-G23, S29-A33, L63-L65, Q68-S76, F83-T89, T97-T105, and K107-K115, suggesting a novel aggregation mechanism possibly occur for the formation of PFN1 aggregates. Overall, these results propose that MD simulations of PFN1 His tautomers can provide a detailed microscopic understanding of the aggregation mechanisms which are hard to probe through experiments.

Biotin and Zn2+ Increase Xylitol Production by Candida tropicalis.

In this study, the medium requirements to increase the production of xylitol by using Candida tropicalis (CT) have been investigated. The technique of single addition or omission of medium components was applied to determine the nutritional requirements. The addition of amino acids such as Asp, Glu, Gln, Asn, Thr, and Gly had no significant effect on CT growth. However, in the absence of other metal ions, there was a higher concentration of cell growth and xylitol production when only Zn2+ was present in the medium. The analysis of various vitamins unveiled that biotin and thiamine were the only vitamins required for the growth of CT. Surprisingly, when only biotin was present in the medium as a vitamin, there was less growth of CT than when the medium was complete, but the amount of xylitol released was significantly higher. Overall, this study will increase the xylitol production using the single omission or addtion technique.

Temperature-dependent conformational dynamics of cytochrome c: Implications in apoptosis.

Heat, electric shock, and burn injuries induce apoptosis by releasing cytochrome c (cyt-c) from mitochondria and by subsequently activating the death protease, caspases-3. During apoptosis, cyt-c undergoes changes in the secondary structure that have been suggested to increase its peroxidase activity. Information about these structural changes will provide better understanding of the apoptotic mechanism. Hence, temperature-dependent conformational dynamics of cyt-c has been investigated through molecular dynamics (MD) simulations to explain the structural changes and to correlate them with its apoptotic behavior. We observe that, at lower temperatures (223, 248, and 300K), the secondary structure of cyt-c, remains stable, while at higher temperatures (323, 373, 423, and 473K), the secondary structural regions change significantly. Further, our MD results indicate that these structural changes are mainly localized on α-helices, turns, ß-sheets, and important loops that were involved in the stabilization of the heme conformation. This conformational transition between specific regions of secondary structure of cyt-c directly affects the electron tunneling properties of the proteins as observed experimentally. We quantify and compare these changes and explain that the temperature plays a vital role in assuring the structural stability of cyt-c and thus its functions. Our findings from this MD study reproduce experimental results at high temperatures and provide evidence for the alteration of the heme through the disruption of the H-bonding interactions between specific regions of cyt-c, thereby enhancing its peroxidase activity which plays a crucial role in the apoptotic process.

Investigation of structural dynamics of Thrombocytopenia Cargeeg mutants of human apoptotic cytochrome c: A molecular dynamics simulation approach.

Naturally occurring mutations to cytochrome c (cyt-c) have been identified recently in patients with mild autosomal dominant thrombocytopenia (low platelet levels), which yield cyt-c mutants with enhanced apoptotic activity. However, the molecular mechanism underlying this low platelet production and enhanced apoptosis remain unclear. Therefore, an attempt is made herein for the first time to investigate the effects of mutations of glycine 41 by serine (G41S) and tyrosine 48 by histidine (Y48H) on the conformational and dynamic changes of apoptotic (Fe3+) cyt-c using all atom molecular dynamics (MD) simulations in explicit water solvent. Our 30ns MD simulations demonstrate considerable structural differences in G41S and Y48H compared to wild type (WT) cyt-c, such as increasing distances between the critical electron transfer residues results in open conformation at the heme active site, large fluctuations in ß-turns and α-helices. Additionally, although the ß-sheets remain mostly unaffected in all the three cyt-c simulations, the α-helices undergo conformational switch to ß-turns in both the mutant simulations. Importantly, this conformational switch of α-helix to ß-turn around heme active site should attributes to the loss of intraprotein H-bonds in the mutant simulations especially between NE2 (His26) and O (Pro44) in agreement with the experimental report. Further, essential dynamics analysis reveals that overall motions of WT cyt-c is mainly involved only in the first eigenvector, but in G41S and Y48H the overall motions are mainly in three and two eigenvectors respectively. Overall, the detailed atomistic level information provide a unifying description for the molecular mechanism of structural destabilization, disregulation of platelet formation and enhanced peroxidase activity of the mutant cyt-c's in the pathology of intrinsic apoptosis.

Effects of different solvents on the conformations of apoptotic cytochrome c: Structural insights from molecular dynamics simulation.

Cytochrome c (cyt-c) upon binding with cardiolipin acquires peroxidase activity and is strictly connected to the pathogenesis of many human diseases including neurodegenerative and cardiovascular diseases. Interaction of cyt-c with cardiolipin mimics partial unfolding/conformational changes of cyt-c in different solvent environments. Dynamic pictures of these conformational changes of cyt-c are crucial in understanding their physiological roles in mitochondrial functions. Therefore, atomistic molecular dynamics (MD) simulations have been carried out to investigate the effect of different solvents (water, urea/water, MeOH and DMSO) on the structure and conformations of apoptotic cyt-c (Fe3+). Our study demonstrates that the structural changes in the protein are solvent dependent. The structural differences are observed majorly on the ß-sheets and α-helical conformations and the degree of their perturbation are specific to the solvent. Although a complete loss of ß-sheets (0%) is observed in MeOH and DMSO, by contrast, well preserved ß-sheets (3.84%) are observed in water and urea/water. A significant decrease in the α-helical contents is observed in MeOH (41.34%) and water (42.46%), a negligible alteration in DMSO (44.25%) and well preserved α-helical (45.19%) contents in urea/water. The distances between the residues critical for electron transfer are decreased considerably for DMSO. Further, the reduction in residue flexibility and the conformational space indicate that the collective motions of cyt-c are reduced when compared to other cosolvents. Essential dynamics analysis implies that the overall motions of cyt-c in water, MeOH and urea/water are involved in three to four eigenvectors and in first eigenvector in DMSO. Overall, we believe that MD simulations of cyt-c in different solvents can provide a detailed microscopic understanding of the physiological roles, electron transport and peroxidase function in the early events of apoptosis which are hard to probe experiments.

Comparative structural and conformational studies on H43R and W32F mutants of copper-zinc superoxide dismutase by molecular dynamics simulation.

Recently, mutations in copper-zinc superoxide dismutase (SOD1) have been linked to familial amyotrophic lateral sclerosis (fALS), a progressive neurodegenerative disease involving motor neuron loss, paralysis and death. It is mainly due to protein misfolding and aggregation resulting from the enhanced peroxidase activity of SOD1 mutants. In this study, we have carried out a 20 ns molecular dynamics simulation for wild type (WT), H43R and W32F mutated SOD1's dimer and compared their structure and conformational properties by extracting several quantitative properties from the trajectory to understand the pathology of fALS disease. Our results show considerable differences in H43R compared to WT and W32F mutated SOD1, such as increasing distances between the critical residues results in open conformation at the active site, strong fluctuations in the important loops (Zinc and electrostatic loops) and weakening of important hydrogen bonds especially between N (His 43/Arg 43) and carbonyl oxygen (His 120) in agreement with the experimental report. The calculated buried surface area of dimer interface for WT, H43R and W32F are 682, 726 and 657 Å(2) respectively, representing the loss of dimerization in H43R. Essential dynamics reveal that overall motions of WT and W32F are mainly involved in three to four eigenvectors, but in H43R the overall motions are mainly in the first eigenvector. These data thus provide a unifying description for the structural destabilization, enhanced peroxidase activity, loss of dismutation activity and increase in aggregation propensity in the pathology of fALS diseases.

Theoretical investigation of quinone metabolites of dopamine interaction with DNA--insights into toxicological effects.

Dopamine (3,4-dihydroxyphenylethylamine, DA), an important neurotransmitter, exists in the cell bodies of the dopaminergic neurons of the substantia nigra. Oxidation of DA to its quinone and subsequent reaction with Adenine and Guanine in DNA result in the formation of depurinating adducts, thus causing DNA damage. In this article, we investigate the interaction of quinone metabolites of dopamine (DMQ) with models representing the structure of DNA using dispersion corrected density functional theory with an aim to evaluate the associated structural changes in DNA upon their interaction. Various binding sites for the DA metabolite on these DNA models have been considered and our computations on the activation barriers allowed us to identify preferential bonding sites for these metabolites analogous to experiments. Analysis of the geometry of these adducts in comparison to free base pairs reveals that the attack of DMQ causes remarkable changes in the structural properties. With our calculations, we propose that these structural alterations induce mutations by favoring the formation of depurinating adducts leading to mutagenic effects such as base mispairing, explaining the toxicological (carcinogenic and neurotoxic) behavior of DMQ.