Inhibition mechanism of testis-expressed gene 14 (TEX14) in cytokinetic abscission: Well-tempered metadynamics simulation studies.

Cytokinesis requires a apoptosis-linked gene 2 interacting protein X (ALIX) and a 55 kDa midbody centrosomal protein (CEP55) to activate the cell abscission in somatic cells. However, in germ cells, CEP55 forms intercellular bridges with testis-expressed gene 14 (TEX14), which blocks the cell abscission. These intercellular bridges play important roles in the synchronization of the germ cells and facilitate the coordinated passage of organelles and molecules between germ cells. If TEX14 is intentionally removed, intercellular bridges are disrupted, leading to sterility. Hence, a deeper understanding regarding the roles of TEX14 can provide significant insights into the inactivation of abscission and the inhibition of proliferation in cancer cells. Previous experimental studies have shown that the high affinity and low dissociation rate of TEX14 for CEP55 prevent ALIX from binding CEP55 and inactivate the germ cell abscission. However, detailed information about how TEX14 interacts with CEP55 to prevent the cell abscission is still lacking. To gain more specific insights into the interactions between CEP55 and TEX14 and the difference in reactivity between TEX14 and ALIX, we performed well-tempered metadynamics simulations of these protein complexes using atomistic models of CEP55, TEX14, and ALIX. We identified the major binding residues of TEX14 and ALIX with CEP55 by using 2D Gibbs free energy evaluations, the results of which are consistent with previous experimental studies. Our results may help design synthetic TEX14 mimicking peptides, which can bind CEP55 and facilitate the inactivation of abscission in abnormal cells, including cancer cells.

Attitudinal analysis of vaccination effects to lead endemic phases.

To achieve endemic phases, repeated vaccinations are necessary. However, individuals may grapple with whether to get vaccinated due to potential side effects. When an individual is already immune due to previous infections or vaccinations, the perceived risk from vaccination is often less than the risk of infection. Yet, repeated rounds of vaccination can lead to avoidance, impeding the establishment of endemic phases. We explore this phenomenon using an individual-based Monte Carlo simulation, validating our findings with game theory. The Nash equilibrium encapsulates individuals' non-cooperative behavior, while the system's optimal value represents the societal benefits of altruistic cooperation. We define the difference between these as the price of anarchy. Our simulations reveal that the price of anarchy must fall below a threshold of 12.47 for endemic phases to be achieved in a steady state. This suggests that for a basic reproduction number of 10, a consistent vaccination rate greater than 89% is required. These findings offer new insights into vaccination-related decision-making and can inform effective strategies to tackle infectious diseases.

Structure and stability of polydiacetylene membrane systems: Molecular dynamics simulation studies.

We have performed full atomistic molecular dynamics (MD) simulations to investigate structure and stability of bilayer membrane systems consisting of monomeric or polymeric 10,12-pentacosadiynoic acid (PCDA) units connected with lysine groups by amide bonds. The PCDA monomer molecules show a twisted three-rod-domain structure with two kinks but upon polymerization, they possess more elongated conformation. The resulting polydiacetylene (PDA) membrane systems have stable membrane structures at room temperature, which is similar to biological lipid bilayer membranes and maintain their gel-like membrane integrity even up to as high as 370 K. Structural properties such as area per monomer, membrane thickness, density profile, 2D pair distribution function, and orientational correlation function are also calculated to understand the membrane structure and check its stability upon thermal fluctuation with atomistic resolution. This study is expected to provide the understanding about PDA membrane systems in atomistic details as well as significant insights into designing new novel PDA sensors.

Clioquinol as an inhibitor of JmjC-histone demethylase exhibits common and unique histone methylome and transcriptome between clioquinol and hypoxia.

Clioquinol (CQ) is a hypoxic mimicker to activate hypoxia-inducible factor-1α (HIF-1α) by inhibiting HIF-1α specific asparaginyl hypoxylase (FIH-1). The structural similarity of the Jumonji C (JmjC) domain between FIH-1 and JmjC domain-containing histone lysine demethylases (JmjC-KDMs) led us to investigate whether CQ could inhibit the catalytic activities of JmjC-KDMs. Herein, we showed that CQ inhibits KDM4A/C, KDM5A/B, and KDM6B and affects H3K4me3, H3K9me3, and H3K27me3 marks, respectively. An integrative analysis of the histone methylome and transcriptome data revealed that CQ-mediated JmjC-KDM inhibition altered the transcription of target genes through differential combinations of KDMs and transcription factors. Notably, functional enrichment of target genes showed that CQ and hypoxia commonly affected the response to hypoxia, VEGF signaling, and glycolysis, whereas CQ uniquely altered apoptosis/autophagy and cytoskeleton/extracellular matrix organization. Our results suggest that CQ can be used as a JmjC-KDM inhibitor, HIF-α activator, and an alternative therapeutic agent in hypoxia-based diseases.

Chlorosulfolipid (Danicalipin A) Membrane Structure: Hybrid Molecular Dynamics Simulation Studies.

Chlorosulfolipids (CSLs) are major components of flagellar membranes in sea algae. Unlike typical biological lipids, CSLs contain hydrophilic sulfate and chloride groups in the hydrocarbon tail; this has deterred the prediction of the CSL membrane structure since 1960. In this study, we combine coarse-grained (CG) and atomistic molecular dynamics (MD) simulations to gain significant insights into the membrane structure of Danicalipin A, which is one of the typical CSLs. It is observed from the CG MD that Danicalipin A lipids form a stable monolayer membrane structure wherein the hydrocarbon moieties are sandwiched by hydrophilic sulfate and chloride groups in both the head and tail regions. On the basis of the mesoscopic structure, we built the corresponding atomistic model to investigate the integrity of the CSL monolayer membrane structure. The monolayer membrane comprising bent lipids shows high thermal stability up to 313 K. The gel-liquid crystalline phase transition is observed around 300 K.

Interstitially Mixed Self-Assembled Monolayers Enhance Electrical Stability of Molecular Junctions.

Electrical breakdown is a critical problem in electronics. In molecular electronics, it becomes more problematic because ultrathin molecular monolayers have delicate and defective structures and exhibit intrinsically low breakdown voltages, which limit device performances. Here, we show that interstitially mixed self-assembled monolayers (imSAMs) remarkably enhance electrical stability of molecular-scale electronic devices without deteriorating function and reliability. The SAM of the sterically bulky matrix (SC11BIPY rectifier) molecule is diluted with a skinny reinforcement (SCn) molecule via the new approach, so-called repeated surface exchange of molecules (ReSEM). Combined experiments and simulations reveal that the ReSEM yields imSAMs wherein interstices between the matrix molecules are filled with the reinforcement molecules and leads to significantly enhanced breakdown voltage inaccessible by traditional pure or mixed SAMs. Thanks to this, bias-driven disappearance and inversion of rectification is unprecedentedly observed. Our work may help to overcome the shortcoming of SAM's instability and expand the functionalities.

Application of molecular dynamics simulation to improve the theoretical prediction for collisional cross section of aromatic compounds with long alkyl chains in crude oils.

RATIONALE: Molecular dynamics (MD) simulations with finite temperature were performed to improve the theoretical prediction of collisional cross section (CCS) values, especially for aromatic compounds containing long alkyl chains. METHODS: In this study, the CCS values of 11 aromatic compounds with long alkyl chains were calculated by MD simulations while considering internal energy at 300, 500, and 700 K, and the results were compared with experimentally determined values. RESULTS: The CCS values calculated at higher energies showed better agreement with the experimental values. Polycyclic aromatic hydrocarbons (PAHs) such as pentacene and benz[b]anthracene were also investigated, and better agreement between the theoretical and experimental results was observed when higher temperature (or higher internal energy) was considered. CONCLUSIONS: The data presented in this study show that the internal degrees of freedom of ions must be considered to accurately predict the CCS values of aromatic compounds with a flexible structure measured by ion mobility mass spectrometry.

A motor neuron strategy to save time and energy in neurodegeneration: adaptive protein stoichiometry.

Neurofilament proteins (Nf) are a biomarker of disease progression in amyotrophic lateral sclerosis (ALS). This study investigated whether there are major differences in expression from in vivo measurements of neurofilament isoforms, from the light chain, NfL (68 kDa), compared with larger proteins, the medium chain (NfM, 150 kDa) and the heavy (NfH, 200-210 kDa) chains in ALS patients and healthy controls. New immunological methods were combined with Nf subunit stoichiometry calculations and Monte Carlo simulations of a coarse-grained Nf brush model. Based on a physiological Nf subunit stoichiometry of 7 : 3 : 2 (NfL:NfM:NfH), we found an 'adaptive' Nf subunit stoichiometry of 24 : 2.4 : 1.6 in ALS. Adaptive Nf stoichiometry preserved NfL gyration radius in the Nf brush model. The energy and time requirements for Nf translation were 56 ± 27k ATP (5.6 h) in control subjects compared to 123 ± 102k (12.3 h) in ALS with 'adaptive' (24:2.4:1.6) Nf stoichiometry (not significant) and increased significantly to 355 ± 330k (35.5 h) with 'luxury' (7:3:2) Nf subunit stoichiometry (p < 0.0001 for each comparison). Longitudinal disease progression-related energy consumption was highest with a 'luxury' (7:3:2) Nf stoichiometry. Therefore, an energy and time-saving option for motor neurons is to shift protein expression from larger to smaller (cheaper) subunits, at little or no costs on a protein structural level, to compensate for increased energy demands.

On-Demand Modulation of Bacterial Cell Fates on Multifunctional Dynamic Substrates.

This paper reports unprecedented dynamic surfaces based on zwitterionic low-density self-assembled monolayers (LDSAMs) of alkanethiolates on gold, which integrate three interconvertible states-bacteria-adherable, bactericidal, and nonfouling states-through electrical modulations. The conformations of alkanethiolates were electrically modulated to generate zwitterionic, anionic, and cationic surfaces, which responded differently to bacteria and determined the fate of bacteria. Furthermore, the reversible switching of multifunctions of the surface was realized for killing bacteria and subsequently releasing dead bacteria from the surface. For practical application of our strategy, we examined the selective antibacterial effect of our surface for eradication of mycoplasma contaminants in contaminated mammalian cell cultures.

Nanochannel-Confined TAMRA-Polypyrrole Stained DNA Stretching by Varying the Ionic Strength from Micromolar to Millimolar Concentrations.

Large DNA molecules have been utilized as a model system to investigate polymer physics. However, DNA visualization via intercalating dyes has generated equivocal results due to dye-induced structural deformation, particularly unwanted unwinding of the double helix. Thus, the contour length increases and the persistence length changes so unpredictably that there has been a controversy. In this paper, we used TAMRA-polypyrrole to stain single DNA molecules. Since this staining did not change the contour length of B-form DNA, we utilized TAMRA-polypyrrole stained DNA as a tool to measure the persistence length by changing the ionic strength. Then, we investigated DNA stretching in nanochannels by varying the ionic strength from 0.06 mM to 47 mM to evaluate several polymer physics theories proposed by Odijk, de Gennes and recent papers to deal with these regimes.

Structural and biochemical insights into the role of testis-expressed gene 14 (TEX14) in forming the stable intercellular bridges of germ cells.

Intercellular bridges are a conserved feature of spermatogenesis in mammalian germ cells and derive from arresting cell abscission at the final stage of cytokinesis. However, it remains to be fully understood how germ cell abscission is arrested in the presence of general cytokinesis components. The TEX14 (testis-expressed gene 14) protein is recruited to the midbody and plays a key role in the inactivation of germ cell abscission. To gain insights into the structural organization of TEX14 at the midbody, we have determined the crystal structures of the EABR [endosomal sorting complex required for transport (ESCRT) and ALIX-binding region] of CEP55 bound to the TEX14 peptide (or its chimeric peptides) and performed functional characterization of the CEP55-TEX14 interaction by multiexperiment analyses. We show that TEX14 interacts with CEP55-EABR via its AxGPPx3Y (Ala793, Gly795, Pro796, Pro797, and Tyr801) and PP (Pro803 and Pro804) sequences, which together form the AxGPPx3YxPP motif. TEX14 competitively binds to CEP55-EABR to prevent the recruitment of ALIX, which is a component of the ESCRT machinery with the AxGPPx3Y motif. We also demonstrate that a high affinity and a low dissociation rate of TEX14 to CEP55, and an increase in the local concentration of TEX14, cooperatively prevent ALIX from recruiting ESCRT complexes to the midbody. The action mechanism of TEX14 suggests a scheme of how to inactivate the abscission of abnormal cells, including cancer cells.

Rational design of paraoxonase 1 (PON1) for the efficient hydrolysis of organophosphates.

A rational design of paraoxonase 1 based on molecular docking discovered H115W/T332S and I74F/H115W/T332S mutants exhibited a 40-fold increase in catalytic efficiency (kcat/Km) toward the hydrolysis of two toxic and popular organophosphates (diethyl-paraoxon and dimethyl-paraoxon). Moreover, the conversion of the paraoxons (741.3-825.6 mg L(-1)) by the evolved mutants was 42-60-fold faster than that by the wild type.

Insights into the Lactonase Mechanism of Serum Paraoxonase 1 (PON1): Experimental and Quantum Mechanics/Molecular Mechanics (QM/MM) Studies.

Serum paraoxonase 1 (PON1) is a versatile enzyme for the hydrolysis of various substrates (e.g., lactones, phosphotriesters) and for the formation of a promising chemical platform γ-valerolactone. Elucidation of the PON1-catalyzed lactonase reaction mechanism is very important for understanding the enzyme function and for engineering this enzyme for specific applications. Kinetic study and hybrid quantum mechanics/molecular mechanics (QM/MM) method were used to investigate the PON1-catalyzed lactonase reaction of γ-butyrolactone (GBL) and (R)-γ-valerolactone (GVL). The activation energies obtained from the QM/MM calculations were in good agreement with the experiments. Interestingly, the QM/MM energy barriers at MP2/3-21G(d,p) level for the lactonase of GVL and GBL were respectively 14.3-16.2 and 11.5-13.1 kcal/mol, consistent with the experimental values (15.57 and 14.73 kcal/mol derived from respective kcat values of 36.62 and 147.21 s(-1)). The QM/MM energy barriers at MP2/6-31G(d) and MP2/6-31G(d,p) levels were also in relatively good agreements with the experiments. Importantly, the difference in the QM/MM energy barriers at MP2 level with all investigated basis sets for the lactonase of GVL and GBL were in excellent agreement with the experiments (0.9-3.1 and 0.8 kcal/mol, respectively). A detailed mechanism for the PON1-catalyzed lactonase reaction was also proposed in this study.

Structural insights into the efficient CO2-reducing activity of an NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA.

CO2 fixation is thought to be one of the key factors in mitigating global warming. Of the various methods for removing CO2, the NAD-dependent formate dehydrogenase from Candida boidinii (CbFDH) has been widely used in various biological CO2-reduction systems; however, practical applications of CbFDH have often been impeded owing to its low CO2-reducing activity. It has recently been demonstrated that the NAD-dependent formate dehydrogenase from Thiobacillus sp. KNK65MA (TsFDH) has a higher CO2-reducing activity compared with CbFDH. The crystal structure of TsFDH revealed that the biological unit in the asymmetric unit has two conformations, i.e. open (NAD(+)-unbound) and closed (NAD(+)-bound) forms. Three major differences are observed in the crystal structures of TsFDH and CbFDH. Firstly, hole 2 in TsFDH is blocked by helix α20, whereas it is not blocked in CbFDH. Secondly, the sizes of holes 1 and 2 are larger in TsFDH than in CbFDH. Thirdly, Lys287 in TsFDH, which is crucial for the capture of formate and its subsequent delivery to the active site, is an alanine in CbFDH. A computational simulation suggested that the higher CO2-reducing activity of TsFDH is owing to its lower free-energy barrier to CO2 reduction than in CbFDH.

Determination of self-exchange rate of alkanethiolates in self-assembled monolayers on gold using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry.

In this paper, we describe a new method for determining the exchange rates of alkanethiolates in self-assembled monolayers (SAMs) on gold using matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) to analyze the compositions of the alkanethiolate in SAMs rapidly and directly. In particular, to investigate the self-exchange of alkanethiols, we prepared a deuterated alkanethiol that has the same molecular properties as the non-deuterated alkanethiol but a different molecular weight. SAMs consisting of deuterated alkanethiolates were immersed in a solution of the non-deuterated alkanethiol, and the influences of the immersion time, temperature, concentration, and solvent on the self-exchange rates were investigated. Furthermore, we assessed the exchange rates among alkanethiols with different carbon chain lengths and different size of ethylene glycol units. In addition, we performed molecular dynamics simulations using a model SAM system in order to understand the molecular mechanism of the exchange process.

Free energy of PAMAM dendrimer adsorption onto model biological membranes.

We investigated the thermodynamic, structural, and dynamics changes in dendrimer-membrane systems during dendrimer adsorption to biological membrane systems by combining atomistic molecular dynamics simulations with umbrella sampling techniques to understand the atomistic interactions between the dendrimer and biological membranes. An ethylenediamine core polyamidoamine dendrimer (generation 3) with amine terminal groups and both zwitterionic dipalmitoyl-phosphatidyl-choline (DPPC) and anionic palmitoyl-oleoyl-phosphatidyl glycerol (POPG) lipid bilayer membranes were used as the model dendrimer and biological membranes, respectively, in this study. The free energy of the dendrimer adsorption onto two model membranes with different charge states was quantitatively determined. For the zwitterionic DPPC membrane, the dendrimer has a minimum free energy of approximately 50 kcal/mol, which is 15 kcal/mol higher than that observed in previous studies. The dominant contribution to the adsorption potential energy is the van der Waals attraction between the dendrimer and the DPPC membrane. However, the anionic POPG membrane pulls the positively charged dendrimer with an attractive mean force of about 200 pN, finally positioning the dendrimer in the membrane headgroup region. As a result of these strong attractive dendrimer and membrane interactions, the dendrimer structurally undergoes the transition from spherical to a pancake conformation, which slows its lateral mobility, especially in the presence of the POPG membrane. The bilayer lipid membranes are also perturbed by the dendrimer adsorption.

Nanoslit Confined DNA at Low Ionic Strengths.

We present nanoslit confined DNA conformations at very low ionic strengths and a theory to explain most measurements for single DNA molecule size under strong nanoslit confinement. Very low ionic strength conditions not only increase the DNA persistence length dramatically, but also cause DNA molecules to swell to the extent that the effective diameter of DNA becomes larger than the nanoslit height. By accounting for these effects, our results and theory provide a reasonable clue for a current controversy regarding the dependence of the DNA conformation on slit height (h), persistence length (p), and effective diameter (w).

Systems-level modeling with molecular resolution elucidates the rate-limiting mechanisms of cellulose decomposition by cellobiohydrolases.

Interprotein and enzyme-substrate couplings in interfacial biocatalysis induce spatial correlations beyond the capabilities of classical mass-action principles in modeling reaction kinetics. To understand the impact of spatial constraints on enzyme kinetics, we developed a computational scheme to simulate the reaction network of enzymes with the structures of individual proteins and substrate molecules explicitly resolved in the three-dimensional space. This methodology was applied to elucidate the rate-limiting mechanisms of crystalline cellulose decomposition by cellobiohydrolases. We illustrate that the primary bottlenecks are slow complexation of glucan chains into the enzyme active site and excessive enzyme jamming along the crowded substrate. Jamming could be alleviated by increasing the decomplexation rate constant but at the expense of reduced processivity. We demonstrate that enhancing the apparent reaction rate required a subtle balance between accelerating the complexation driving force and simultaneously avoiding enzyme jamming. Via a spatiotemporal systems analysis, we developed a unified mechanistic framework that delineates the experimental conditions under which different sets of rate-limiting behaviors emerge. We found that optimization of the complexation-exchange kinetics is critical for overcoming the barriers imposed by interfacial confinement and accelerating the apparent rate of enzymatic cellulose decomposition.

Conformational properties of interacting neurofilaments: Monte Carlo simulations of cylindrically grafted apposing neurofilament brushes.

Neurofilaments are essential cytoskeletal filaments that impart mechanical stability to axons. They are mostly assembled from three neurofilament proteins that form the core of the filament and its sidearms. Adjacent neurofilaments interact with each other through their apposing sidearms and attain unique conformations depending on the ionic condition, phosphorylation state, and interfilament separations. To understand the conformational properties of apposing sidearms under various conditions and gain insight into interfilament interactions, we performed Monte Carlo simulations of neurofilament pairs. We employed a sequence-based coarse-grained model of apposing NF sidearms that are end-tethered to cylindrical geometries according to the stoichiometry of the three neurofilament subunits. Monte Carlo simulations were conducted under different conditions such as phosphorylation state, ionic condition, and interfilament separations. Under salt-free conditions, apposing sidearms are found to adopt mutually excluding stretched but bent away conformations that are reminiscent of a repulsive type of interaction. Under physiological conditions, apposing sidearms are found to be in a coiled conformation, suggesting a short-range steric repulsive type of interaction. Increased sidearm mutual interpenetration and a simultaneous decrease in the individual brush heights were observed as the interfilament separation was reduced from 60 to 40 nm. The observed conformations suggest entropic interaction as a likely mechanism for sidearm-mediated interfilament interactions under physiological conditions.

Lipase-catalyzed enantioselective synthesis of (R,R)-lactide from alkyl lactate to produce PDLA (poly D-lactic acid) and stereocomplex PLA (poly lactic acid).

R-lactide, a pivotal monomer for the production of poly (D-lactic acid) (PDLA) or stereocomplex poly (lactic acid) (PLA) was synthesized from alkyl (R)-lactate through a lipase-catalyzed reaction without racemization. From among several types of lipase, only lipase B from Candida antarctica (Novozym 435; CAL-B) was effective in the reaction that synthesized (R,R)-lactide. Enantiopure (R,R)-lactide, which consisted of over 99% enantiomeric excess, was synthesized from methyl (R)-lactate through CAL-B catalysis. Removal of the methanol by-product was critical to obtain a high level of lactide conversion. The (R,R)-lactide yield was 56% in a reaction containing 100 mg of Novozym 435, 10 mM methyl (R)-lactate and 1500 mg of molecular sieve 5A in methyl tert-butyl ether (MTBE). The important monomer (R,R)-lactide that is required for the production of the widely recognized bio-plastic PDLA and the PLA stereocomplex can be obtained using this novel synthetic method.