Mechanism of Substrate Translocation in an Alternating Access Transporter.

Transporters shuttle molecules across cell membranes by alternating among distinct conformational states. Fundamental questions remain about how transporters transition between states and how such structural rearrangements regulate substrate translocation. Here, we capture the translocation process by crystallography and unguided molecular dynamics simulations, providing an atomic-level description of alternating access transport. Simulations of a SWEET-family transporter initiated from an outward-open, glucose-bound structure reported here spontaneously adopt occluded and inward-open conformations. Strikingly, these conformations match crystal structures, including our inward-open structure. Mutagenesis experiments further validate simulation predictions. Our results reveal that state transitions are driven by favorable interactions formed upon closure of extracellular and intracellular "gates" and by an unfavorable transmembrane helix configuration when both gates are closed. This mechanism leads to tight allosteric coupling between gates, preventing them from opening simultaneously. Interestingly, the substrate appears to take a "free ride" across the membrane without causing major structural rearrangements in the transporter.

A Zn-dependent structural transition of SOD1 modulates its ability to undergo phase separation.

The misfolding and mutation of Cu/Zn superoxide dismutase (SOD1) is commonly associated with amyotrophic lateral sclerosis (ALS). SOD1 can accumulate within stress granules (SGs), a type of membraneless organelle, which is believed to form via liquid-liquid phase separation (LLPS). Using wild-type, metal-deficient, and different ALS disease mutants of SOD1 and computer simulations, we report here that the absence of Zn leads to structural disorder within two loop regions of SOD1, triggering SOD1 LLPS and amyloid formation. The addition of exogenous Zn to either metal-free SOD1 or to the severe ALS mutation I113T leads to the stabilization of the loops and impairs SOD1 LLPS and aggregation. Moreover, partial Zn-mediated inhibition of LLPS was observed for another severe ALS mutant, G85R, which shows perturbed Zn-binding. By contrast, the ALS mutant G37R, which shows reduced Cu-binding, does not undergo LLPS. In addition, SOD1 condensates induced by Zn-depletion exhibit greater cellular toxicity than aggregates formed by prolonged incubation under aggregating conditions. Overall, our work establishes a role for Zn-dependent modulation of SOD1 conformation and LLPS properties that may contribute to amyloid formation.

Atomistic simulations reveal impacts of missense mutations on the structure and function of SynGAP1.

De novo mutations in the synaptic GTPase activating protein (SynGAP) are associated with neurological disorders like intellectual disability, epilepsy, and autism. SynGAP is also implicated in Alzheimer's disease and cancer. Although pathogenic variants are highly penetrant in neurodevelopmental conditions, a substantial number of them are caused by missense mutations that are difficult to diagnose. Hence, in silico mutagenesis was performed for probing the missense effects within the N-terminal region of SynGAP structure. Through extensive molecular dynamics simulations, encompassing three 150-ns replicates for 211 variants, the impact of missense mutations on the protein fold was assessed. The effect of the mutations on the folding stability was also quantitatively assessed using free energy calculations. The mutations were categorized as potentially pathogenic or benign based on their structural impacts. Finally, the study introduces wild-type-SynGAP in complex with RasGTPase at the inner membrane, while considering the potential effects of mutations on these key interactions. This study provides structural perspective to the clinical assessment of SynGAP missense variants and lays the foundation for future structure-based drug discovery.

A novel mRNA multi-epitope vaccine of Acinetobacter baumannii based on multi-target protein design in immunoinformatic approach.

Acinetobacter baumannii is a gram-negative bacillus prevalent in nature, capable of thriving under various environmental conditions. As an opportunistic pathogen, it frequently causes nosocomial infections such as urinary tract infections, bacteremia, and pneumonia, contributing to increased morbidity and mortality in clinical settings. Consequently, developing novel vaccines against Acinetobacter baumannii is of utmost importance. In our study, we identified 10 highly conserved antigenic proteins from the NCBI and UniProt databases for epitope mapping. We subsequently screened and selected 8 CTL, HTL, and LBL epitopes, integrating them into three distinct vaccines constructed with adjuvants. Following comprehensive evaluations of immunological and physicochemical parameters, we conducted molecular docking and molecular dynamics simulations to assess the efficacy and stability of these vaccines. Our findings indicate that all three multi-epitope mRNA vaccines designed against Acinetobacter baumannii are promising; however, further animal studies are required to confirm their reliability and effectiveness.

Design of Cryptococcus neoformans multi-epitope vaccine based on immunoinformatics method.

Cryptococcus neoformans is a widely distributed opportunistic pathogenic fungus. While C. neoformans commonly infects immunocompromised individuals, it can also affect those who are immunocompetent. Transmission of C. neoformans primarily occurs through the respiratory tract, leading to the development of meningitis. The mortality rate of Cryptococcal meningitis is high, and treatment options are limited. Cryptococcus neoformans infections pose a significant public health threat and currently lack targeted and effective response strategies. This study aimed to screen T lymphocyte (cytotoxic T lymphocyte and helper T lymphocyte) and B lymphocyte epitopes derived from four C. neoformans antigens and develop two multi-epitope vaccines by combining them with various adjuvants. Molecular docking results demonstrated that the vaccines bind stably to Toll-like receptor 4 ( and induce innate immunity. The credibility of the molecular docking results was validated through subsequent molecular dynamics simulations. Furthermore, the results of immune simulation analyses underscored the multi-epitope vaccine's capability to effectively induce robust humoral and cellular immune responses within the host organism. These two vaccines have demonstrated theoretical efficacy against C. neoformans infection as indicated by computer analysis. Nevertheless, additional experimental validation is essential to substantiate the protective efficacy of the vaccines.

A multi-epitope Cryptococcus neoformans vaccine covering the most common A and D phenotypes was designed using bioinformatics methods.

Free Energy Evaluation of Cavity Formation in Metastable Liquid Based on Stochastic Thermodynamics.

Nucleation is a fundamental and general process at the initial stage of first-order phase transition. Although various models based on the classical nucleation theory (CNT) have been proposed to explain the energetics and kinetics of nucleation, detailed understanding at nanoscale is still required. Here, in view of the homogeneous bubble nucleation, we focus on cavity formation, in which evaluation of the size dependence of free energy change is the key issue. We propose the application of a formula in stochastic thermodynamics, the Jarzynski equality, for data analysis of molecular dynamics (MD) simulation to evaluate the free energy of cavity formation. As a test case, we performed a series of MD simulations with a Lennard-Jones (LJ) fluid system. By applying an external spherical force field to equilibrated LJ liquid, we evaluated the free energy change during cavity growth as the Jarzynski's ensemble average of required works. A fairly smooth free energy curve was obtained as a function of bubble radius in metastable liquid of mildly negative pressure conditions.

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective.

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein-protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.

Bioinformatics insights into CENP-T and CENP-W protein-protein interaction disruptive amino acid substitution in the CENP-T-W complex.

Kinetochores are multi-protein assemblies present at the centromere of the human chromosome and play a crucial role in cellular mitosis. The CENP-T and CENP-W chains form a heterodimer, which is an integral part of the inner kinetochore, interacting with the linker DNA on one side and the outer kinetochore on the other. Additionally, the CENP-T-W dimer interacts with other regulatory proteins involved in forming inner kinetochores. The specific roles of different amino acids in the CENP-W at the protein-protein interaction (PPI) interface during the CENP-T-W dimer formation remain incompletely understood. Since cell division goes awry in diseases like cancer, this CENP-T-W partnership is a potential target for new drugs that could restore healthy cell division. We employed molecular docking, binding free energy calculations, and molecular dynamics (MD) simulations to investigate the disruptive effects of amino acids substitutions in the CENP-W chain on CENP-T-W dimer formation. By conducting a molecular docking study and analysing hydrogen bonding interactions, we identified key residues in CENP-W (ASN-46, ARG-53, LEU-83, SER-86, ARG-87, and GLY-88) for further investigation. Through site-directed mutagenesis and subsequent binding free energy calculations, we refined the selection of mutant. We chose four mutants (N46K, R53K, L83K, and R87E) of CENP-W to assess their comparative potential in forming CENP-T-W dimer. Our analysis from 250 ns long revealed that the substitution of LEU83 and ARG53 residues in CENP-W with the LYS significantly disrupts the formation of CENP-T-W dimer. In conclusion, LEU83 and ARG53 play a critical role in CENP-T and CENP-W dimerization which is ultimately required for cellular mitosis. Our findings not only deepen our understanding of cell division but also hint at exciting drug-target possibilities.

Sodium-Dependent Neutral Amino Acid Transporter 2 Can Serve as a Tertiary Carrier for l-Type Amino Acid Transporter 1-Utilizing Prodrugs.

Membrane transporters are the key determinants of the homeostasis of endogenous compounds in the cells and their exposure to drugs. However, the substrate specificities of distinct transporters can overlap. In the present study, the interactions of l-type amino acid transporter 1 (LAT1)-utilizing prodrugs with sodium-coupled neutral amino acid transporter 2 (SNAT2) were explored. The results showed that the cellular uptake of LAT1-utilizing prodrugs into a human breast cancer cell line, MCF-7 cells, was mediated via SNATs as the uptake was increased at higher pH (8.5), decreased in the absence of sodium, and inhibited in the presence of unselective SNAT-inhibitor, (α-(methylamino)isobutyric acid, MeAIB). Moreover, docking the compounds to a SNAT2 homology model (inward-open conformation) and further molecular dynamics simulations and the subsequent trajectory and principal component analyses confirmed the chemical features supporting the interactions of the studied compounds with SNAT2, which was found to be the main SNAT expressed in MCF-7 cells.

Novel and Potential Small Molecule Scaffolds as DYRK1A Inhibitors by Integrated Molecular Docking-Based Virtual Screening and Dynamics Simulation Study.

The dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) is a novel, promising and emerging biological target for therapeutic intervention in neurodegenerative diseases, especially in Alzheimer's disease (AD). The molMall database, comprising rare, diverse and unique compounds, was explored for molecular docking-based virtual screening against the DYRK1A protein, in order to find out potential inhibitors. Ligands exhibiting hydrogen bond interactions with key amino acid residues such as Ile165, Lys188 (catalytic), Glu239 (gk+1), Leu241 (gk+3), Ser242, Asn244, and Asp307, of the target protein, were considered potential ligands. Hydrogen bond interactions with Leu241 (gk+3) were considered key determinants for the selection. High scoring structures were also docked by Glide XP docking in the active sites of twelve DYRK1A related protein kinases, viz. DYRK1B, DYRK2, CDK5/p25, CK1, CLK1, CLK3, GSK3ß, MAPK2, MAPK10, PIM1, PKA, and PKCα, in order to find selective DYRK1A inhibitors. MM/GBSA binding free energies of selected ligand-protein complexes were also calculated in order to remove false positive hits. Physicochemical and pharmacokinetic properties of the selected six hit ligands were also computed and related with the proposed limits for orally active CNS drugs. The computational toxicity webserver ProTox-II was used to predict the toxicity profile of selected six hits (molmall IDs 9539, 11352, 15938, 19037, 21830 and 21878). The selected six docked ligand-protein systems were exposed to 100 ns molecular dynamics (MD) simulations to validate their mechanism of interactions and stability in the ATP pocket of human DYRK1A kinase. All six ligands were found to be stable in the ATP binding pocket of DYRK1A kinase.

Structure-based virtual screening and computational study towards identification of novel inhibitors of hypoxanthine-guanine phosphoribosyltransferase of Trypanosoma cruzi.

Hypoxanthine-guanine phosphoribosyltransferase (HGPRT) is the key regulatory enzyme of the purine salvage pathway present in the members of trypanosomatids. The parasite solely depends on this pathway for the synthesis of nucleotides due to the absence of the de novo pathway. This study intends to identify putative inhibitors towards Trypanosoma cruzi HGPRT (TcHGPRT). Initial virtual screening was performed with substructures of phosphoribosyl pyrophosphate (PRPP), an original substrate of HGPRT. Twenty compounds that had greater binding energy than the substrate was treated as hits and was further screened and narrowed down through induced fit docking which resulted in top five compounds which was distinguished into two groups based on the ligand occupancy within the PRPP binding site of TcHGPRT. Group-I compounds (PubChem CID 130316561 and 134978234) are analogous to PRPP structure with greater occupancy, were preferred over Group-II compounds which had lesser occupancy than the substrate. However, one compound (22404820) among Group II was chosen for further analysis considering its significant electrostatic interactions. Molecular docking studies revealed the requirement of an electronegative moiety like phosphate group to be present in the ligand due to the presence of metal ions in the substrate binding site. The three chosen compounds along with PRPP were subjected to molecular dynamics analysis, which indicated a strong presence of electrostatic interaction. Considering the dynamic stability of interactions as well as pharmacological properties of ligands based on absorption, distribution, metabolism, excretion prediction, Group-I compounds were selected as lead compounds and were subjected to molecular electrostatic potential analysis to determine the charge distribution of the compound. The overall analysis thus suggests both 130316561 and 134978234 can be used as TcHGPRT inhibitors. Furthermore, these computational results emphasize the requirement of phosphorylated ligands which are essential in mediating electrostatic interactions and to compete with the binding affinity of the original substrate.

A computational study of structural differences of binding of NADP+ and G6P substrates to G6PD Mediterraneanc.563T, G6PD A-c.202A/c.376G, G6PD Cairoc.404C and G6PD Gazac.536A mutations.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency is a common X-linked inherited enzymopathic disorder that may lead to transfusion-requiring acute hemolytic anemia (AHA) triggered by fava beans ingestion, infection or some drugs. The gene encoding for G6PD carries a large number of genetic variants that have varying pathogenicity. We reported on three G6PD variants in the Gaza Strip Palestinian population with differing clinical impacts and frequencies: G6PD Mediterraneanc.563T, African G6PD A-c.202A/c.376G, and G6PD Cairoc.404C. We also identified a novel G6PD missense (Ser179Asn) mutation c.536G > A "G6PD Gaza". In this work we explore the effect of these four genetic variants on the structural and substrate (NADP+ and G6P) binding characteristics of the G6PD enzyme using the Monte Carlo (MC) flexible docking and molecular dynamics (MD) simulation approaches. We report that G6PD A-c.202A/c.376G, G6PD Mediterraneanc.563T, G6PD Cairoc.404C and G6PD Gazac.536A mutations cause significant structural changes in G6PD enzyme to induce conformational instability leading to the loss of binding of one or both substrates and are causative of G6PD deficiency.

In silico analyses on the comparative sensing of SARS-CoV-2 mRNA by the intracellular TLRs of humans.

The coronavirus disease-2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has already resulted in a huge setback to mankind in terms of millions of deaths, while the unavailability of an appropriate therapeutic strategy has made the scenario much more severe. Toll-like receptors (TLRs) are crucial mediators and regulators of host immunity and the role of human cell surface TLRs in SARS-CoV-2 induced inflammatory pathogenesis has been demonstrated recently. However, the functional significance of the human intracellular TLRs including TLR3, 7, 8, and 9 is yet unclear. Hitherto, the involvement of these intracellular TLRs in inducing pro-inflammatory responses in COVID-19 has been reported but the identity of the interacting viral RNA molecule(s) and the corresponding TLRs have not been explored. This study hopes to rationalize the comparative binding of the major SARS-CoV-2 mRNAs to the intracellular TLRs, considering the solvent-based force-fields operational in the cytosolic aqueous microenvironment that predominantly drives these interactions. Our in silico study on the binding of all mRNAs with the intracellular TLRs depicts that the mRNA of NSP10, S2, and E proteins of SARS-CoV-2 are possible virus-associated molecular patterns that bind to TLR3, TLR9, and TLR7, respectively, and trigger downstream cascade reactions. Intriguingly, binding of the viral mRNAs resulted in variable degrees of conformational changes in the ligand-binding domain of the TLRs ratifying the activation of the downstream inflammatory signaling cascade. Taken together, the current study is the maiden report to describe the role of TLR3, 7, and 9 in COVID-19 immunobiology and these could serve as useful targets for the conception of a therapeutic strategy against the pandemic.

Prioritization of potential vaccine candidates and designing a multiepitope-based subunit vaccine against multidrug-resistant Salmonella Typhi str. CT18: A subtractive proteomics and immunoinformatics approach.

Salmonella enterica serovar Typhi (S. Typhi), a causative agent of typhoid fever, is a Gram-negative, human-restricted pathogen that causes significant morbidity and mortality, particularly in developing countries. The currently available typhoid vaccines are not recommended to children below six years of age and have poor long-term efficacy. Due to these limitations and the emerging threat of multidrug-resistance (MDR) strains, the development of a new vaccine is urgently needed. The present study aims to design a multiepitope-based subunit vaccine (MESV) against MDR S. Typhi str. CT18 using a computational-based approach comprising subtractive proteomics and immunoinformatics. Firstly, we investigated the proteome of S. Typhi str. CT18 using subtractive proteomics and identified twelve essential, virulent, host non-homologous, and antigenic outer membrane proteins (OMPs) as potential vaccine candidates with low transmembrane helices (≤1) and molecular weight (≤110 kDa). The OMPs were mapped for cytotoxic T lymphocyte(CTL) epitopes, helper T lymphocyte (HTL) epitopes, and linear B lymphocyte (LBL) epitopes using various immunoinformatics tools and servers. A total of 6, 12, and 11 CTL, HTL, and LBL epitopes were shortlisted, respectively, based on their immunogenicity, antigenicity, allergenicity, toxicity, and hydropathicity potential. Four MESV constructs (MESVCs), MESVC-1, MESVC-2, MESVC-3, and MESVC-4, were designed by linking the CTL, HTL, and LBL epitopes with immune-modulating adjuvants, linkers, and PADRE (Pan HLA DR-binding epitope) sequences. The MESVCs were evaluated for their physicochemical properties, allergenicity, antigenicity, toxicity, and solubility potential to ensure their safety and immunogenic behavior. Secondary and tertiary structures of shortlisted MESVCs (MESVC-1, MESVC-3, and MESVC-4) were predicted, modeled, refined, validated, and then docked with various MHC I, MHC II, and TLR4/MD2 complex. Molecular dynamics (MD) simulation of the final selected MESVC-4 with TLR4/MD2 complex confirms its binding affinity and stability. Codon optimization and in silico cloning verified the translation efficiency and successful expression of MESVC-4 in E. coli str. K12. Finally, the efficiency of MESVC-4 to trigger an effective immune response was assessed by an in silico immune simulation. In conclusion, our findings show that the designed MESVC-4 can elicit humoral and cellular immune responses, implying that it may be used for prophylactic or therapeutic purposes. Therefore, it should be subjected to further experimental validations.

Investigating the Mechanistic Inhibitory Discrepancies of Novel Halogen and Alkyl Di-Substituted Oxadiazole-Based Dibenzo-Azepine-Dione Derivatives on Poly (ADP-Ribose) Polymerase-1.

Numerous studies have established the involvement of Poly (ADP-ribose) Polymerase-1 (PARP-1) in cancer presenting it as an important therapeutic target over recent years. Although homology among the PARP protein family makes selective targeting difficult, two compounds [d11 (0.939 µM) and d21 (0.047 µM)] with disparate inhibitory potencies against PARP-1 were recently identified. In this study, free energy calculations and molecular simulations were used to decipher underlying mechanisms of differential PARP-1 inhibition exhibited by the two compounds. The thermodynamics calculation revealed that compound d21 had a relatively higher ΔGbind than d11. High involvement of van der Waal and electrostatic effects potentiated the affinity of d21 at PARP-1 active site. More so, incorporated methyl moiety in d11 accounted for steric hindrance which, in turn, prevented complementary interactions of key site residues such as TYR889, MET890, TYR896, TYR907. Conformational studies also revealed that d21 is more stabilized for interactions in the active site compared to d11. We believe that findings from this study would provide an important avenue for the development of selective PARP-1 inhibitors.

The Interplay of Cholesterol and Ligand Binding in hTSPO from Classical Molecular Dynamics Simulations.

The translocator protein (TSPO) is a 18kDa transmembrane protein, ubiquitously present in human mitochondria. It is overexpressed in tumor cells and at the sites of neuroinflammation, thus representing an important biomarker, as well as a promising drug target. In mammalian TSPO, there are cholesterol-binding motifs, as well as a binding cavity able to accommodate different chemical compounds. Given the lack of structural information for the human protein, we built a model of human (h) TSPO in the apo state and in complex with PK11195, a molecule routinely used in positron emission tomography (PET) for imaging of neuroinflammatory sites. To better understand the interactions of PK11195 and cholesterol with this pharmacologically relevant protein, we ran molecular dynamics simulations of the apo and holo proteins embedded in a model membrane. We found that: (i) PK11195 stabilizes hTSPO structural fold; (ii) PK11195 might enter in the binding site through transmembrane helices I and II of hTSPO; (iii) PK11195 reduces the frequency of cholesterol binding to the lower, N-terminal part of hTSPO in the inner membrane leaflet, while this impact is less pronounced for the upper, C-terminal part in the outer membrane leaflet, where the ligand binding site is located; (iv) very interestingly, cholesterol most frequently binds simultaneously to the so-called CRAC and CARC regions in TM V in the free form (residues L150-X-Y152-X(3)-R156 and R135-X(2)-Y138-X(2)-L141, respectively). However, when the protein is in complex with PK11195, cholesterol binds equally frequently to the CRAC-resembling motif that we observed in TM I (residues L17-X(2)-F20-X(3)-R24) and to CRAC in TM V. We expect that the CRAC-like motif in TM I will be of interest in future experimental investigations. Thus, our MD simulations provide insight into the structural features of hTSPO and the previously unknown interplay between PK11195 and cholesterol interactions with this pharmacologically relevant protein.

Computational Simulation of HIV Protease Inhibitors to the Main Protease (Mpro) of SARS-CoV-2: Implications for COVID-19 Drugs Design.

SARS-CoV-2 is highly homologous to SARS-CoV. To date, the main protease (Mpro) of SARS-CoV-2 is regarded as an important drug target for the treatment of Coronavirus Disease 2019 (COVID-19). Some experiments confirmed that several HIV protease inhibitors present the inhibitory effects on the replication of SARS-CoV-2 by inhibiting Mpro. However, the mechanism of action has still not been studied very clearly. In this work, the interaction mechanism of four HIV protease inhibitors Darunavir (DRV), Lopinavir (LPV), Nelfinavir (NFV), and Ritonavire (RTV) targeting SARS-CoV-2 Mpro was explored by applying docking, molecular dynamics (MD) simulations, and MM-GBSA methods using the broad-spectrum antiviral drug Ribavirin (RBV) as the negative and nonspecific control. Our results revealed that LPV, RTV, and NFV have higher binding affinities with Mpro, and they all interact with catalytic residues His41 and the other two key amino acids Met49 and Met165. Pharmacophore model analysis further revealed that the aromatic ring, hydrogen bond donor, and hydrophobic group are the essential infrastructure of Mpro inhibitors. Overall, this study applied computational simulation methods to study the interaction mechanism of HIV-1 protease inhibitors with SARS-CoV-2 Mpro, and the findings provide useful insights for the development of novel anti-SARS-CoV-2 agents for the treatment of COVID-19.

Molecular Dynamics of Cobalt Protoporphyrin Antagonism of the Cancer Suppressor REV-ERBß.

Nuclear receptor REV-ERBß is an overexpressed oncoprotein that has been used as a target for cancer treatment. The metal-complex nature of its ligand, iron protoporphyrin IX (Heme), enables the REV-ERBß to be used for multiple therapeutic modalities as a photonuclease, a photosensitizer, or a fluorescence imaging agent. The replacement of iron with cobalt as the metal center of protoporphyrin IX changes the ligand from an agonist to an antagonist of REV-ERBß. The mechanism behind that phenomenon is still unclear, despite the availability of crystal structures of REV-ERBß in complex with Heme and cobalt protoporphyrin IX (CoPP). This study used molecular dynamic simulations to compare the effects of REV-ERBß binding to Heme and CoPP, respectively. The initial poses of Heme and CoPP in complex with agonist and antagonist forms of REV-ERBß were predicted using molecular docking. The binding energies of each ligand were calculated using the MM/PBSA method. The computed binding affinity of Heme to REV-ERBß was stronger than that of CoPP, in agreement with experimental results. CoPP altered the conformation of the ligand-binding site of REV-ERBß, disrupting the binding site for nuclear receptor corepressor, which is required for REV-ERBß to regulate the transcription of downstream target genes. Those results suggest that a subtle change in the metal center of porphyrin can change the behavior of porphyrin in cancer cell signaling. Therefore, modification of porphyrin-based agents for cancer therapy should be conducted carefully to avoid triggering unfavorable effects.

Molecular Dynamics Simulations for Effects of Fluoropolymer Binder Content in CL-20/TNT Based Polymer-Bonded Explosives.

In order to better understand the role of binder content, molecular dynamics (MD) simulations were performed to study the interfacial interactions, sensitivity and mechanical properties of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane/2,4,6-trinitrotoluene (CL-20/TNT) based polymer-bonded explosives (PBXs) with fluorine rubber F2311. The binding energy between CL-20/TNT co-crystal (1 0 0) surface and F2311, pair correlation function, the maximum bond length of the N-NO2 trigger bond, and the mechanical properties of the PBXs were reported. From the calculated binding energy, it was found that binding energy increases with increasing F2311 content. Additionally, according to the results of pair correlation function, it turns out that H-O hydrogen bonds and H-F hydrogen bonds exist between F2311 molecules and the molecules in CL-20/TNT. The length of trigger bond in CL-20/TNT were adopted as theoretical criterion of sensitivity. The maximum bond length of the N-NO2 trigger bond decreased very significantly when the F2311 content increased from 0 to 9.2%. This indicated increasing F2311 content can reduce sensitivity and improve thermal stability. However, the maximum bond length of the N-NO2 trigger bond remained essentially unchanged when the F2311 content was further increased. Additionally, the calculated mechanical data indicated that with the increase in F2311 content, the rigidity of CL-20/TNT based PBXs was decrease, the toughness was improved.

NMR- and MD simulation-based structural characterization of the membrane-associating FATC domain of ataxia telangiectasia mutated.

The Ser/Thr protein kinase ataxia telangiectasia mutated (ATM) plays an important role in the DNA damage response, signaling in response to redox signals, the control of metabolic processes, and mitochondrial homeostasis. ATM localizes to the nucleus and at the plasma membrane, mitochondria, peroxisomes, and other cytoplasmic vesicular structures. It has been shown that the C-terminal FATC domain of human ATM (hATMfatc) can interact with a range of membrane mimetics and may thereby act as a membrane-anchoring unit. Here, NMR structural and 15N relaxation data, NMR data using spin-labeled micelles, and MD simulations of micelle-associated hATMfatc revealed that it binds the micelle by a dynamic assembly of three helices with many residues of hATMfatc located in the headgroup region. We observed that none of the three helices penetrates the micelle deeply or makes significant tertiary contacts to the other helices. NMR-monitored interaction experiments with hATMfatc variants in which two conserved aromatic residues (Phe3049 and Trp3052) were either individually or both replaced by alanine disclosed that the double substitution does not abrogate the interaction with micelles and bicelles at the high concentrations at which these aggregates are typically used, but impairs interactions with small unilamellar vesicles, usually used at much lower lipid concentrations and considered a better mimetic for natural membranes. We conclude that the observed dynamic structure of micelle-associated hATMfatc may enable it to interact with differently composed membranes or membrane-associated interaction partners and thereby regulate ATM's kinase activity. Moreover, the FATC domain of ATM may function as a membrane-anchoring unit for other biomolecules.