Polarizable Atomic Multipole Water Model for Molecular Mechanics Simulation

A new classical empirical potential is proposed for water. The model uses a polarizable atomic multipole description of electrostatic interactions. Multipoles through the quadrupole are assigned to each atomic center based on a distributed multipole analysis (DMA) derived from large basis set molecular orbital calculations on the water monomer. Polarization is treated via self-consistent induced atomic dipoles. A modified version of Thole’s interaction model is used to damp induction at short range. Repulsion-dispersion (vdW) effects are computed from a buffered 14-7 potential. In a departure from most current water potentials, we find that significant vdW parameters are necessary on hydrogen as well as oxygen. The new potential is fully flexible and has been tested versus a variety of experimental data and quantum calculations for small clusters, liquid water, and ice. Overall, excellent agreement with experimental and high level ab initio results is obtained for numerous properties, including cluster structures and energetics and bulk thermodynamic and structural measures. The parametrization scheme described here is easily extended to other molecular systems, and the resulting water potential should provide a useful explicit solvent model for organic solutes and biopolymer modeling

A Simple Crowdfunding Model

Honors (Bachelor's)EconomicsUniversity of Michiganhttps://deepblue.lib.umich.edu/bitstream/2027.42/139649/1/pren.pd

Large Domain Motions in Ago Protein Controlled by the Guide DNA-Strand Seed Region Determine the Ago-DNA-mRNA Complex Recognition Process

The recognition mechanism and cleavage activity of argonaute (Ago), miRNA, and mRNA complexes are the core processes to the small non-coding RNA world. The 5′ nucleation at the ‘seed’ region (position 2–8) of miRNA was believed to play a significant role in guiding the recognition of target mRNAs to the given miRNA family. In this paper, we have performed all-atom molecular dynamics simulations of the related and recently revealed Ago-DNA:mRNA ternary complexes to study the dynamics of the guide-target recognition and the effect of mutations by introducing “damaging” C·C mismatches at different positions in the seed region of the DNA-RNA duplex. Our simulations show that the A-form-like helix duplex gradually distorts as the number of seed mismatches increases and the complex can survive no more than two such mismatches. Severe distortions of the guide-target heteroduplex are observed in the ruinous 4-sites mismatch mutant, which give rise to a bending motion of the PAZ domain along the L1/L2 “hinge-like” connection segment, resulting in the opening of the nucleic-acid-binding channel. These long-range interactions between the seed region and PAZ domain, moderated by the L1/L2 segments, reveal the central role of the seed region in the guide-target strands recognition: it not only determines the guide-target heteroduplex’s nucleation and propagation, but also regulates the dynamic motions of Ago domains around the nucleic-acid-binding channel

PerturbScore: Connecting Discrete and Continuous Perturbations in NLP

With the rapid development of neural network applications in NLP, model robustness problem is gaining more attention. Different from computer vision, the discrete nature of texts makes it more challenging to explore robustness in NLP. Therefore, in this paper, we aim to connect discrete perturbations with continuous perturbations, therefore we can use such connections as a bridge to help understand discrete perturbations in NLP models. Specifically, we first explore how to connect and measure the correlation between discrete perturbations and continuous perturbations. Then we design a regression task as a PerturbScore to learn the correlation automatically. Through experimental results, we find that we can build a connection between discrete and continuous perturbations and use the proposed PerturbScore to learn such correlation, surpassing previous methods used in discrete perturbation measuring. Further, the proposed PerturbScore can be well generalized to different datasets, perturbation methods, indicating that we can use it as a powerful tool to study model robustness in NLP.Comment: Accepted by Findings of EMNLP202

Constant-pH simulations with the polarizable atomic multipole AMOEBA force field

Accurately predicting protein behavior across diverse pH environments remains a significant challenge in biomolecular simulations. Existing constant-pH molecular dynamics (CpHMD) algorithms are limited to fixed-charge force fields, hindering their application to biomolecular systems described by permanent atomic multipoles or induced dipoles. This work overcomes these limitations by introducing the first polarizable CpHMD algorithm in the context of the Atomic Multipole Optimized Energetics for Biomolecular Applications (AMOEBA) force field. Additionally, our implementation in the open-source Force Field X (FFX) software has the unique ability to handle titration state changes for crystalline systems including flexible support for all 230 space groups. The evaluation of constant-pH molecular dynamics (CpHMD) with the AMOEBA force field was performed on 11 crystalline peptide systems that span the titrating amino acids (Asp, Glu, His, Lys, and Cys). Titration states were correctly predicted for 15 out of the 16 amino acids present in the 11 systems, including for the coordination of Z

Trypsin-ligand binding free energies from explicit and implicit solvent simulations with polarizable potential

We have calculated the binding free energies of a series of benzamidine-like inhibitors to trypsin with a polarizable force field using both explicit and implicit solvent approaches. Free energy perturbation has been performed for the ligands in bulk water and in protein complex with molecular dynamics simulations. The calculated binding free energies are well within the accuracy of experimental measurement and the direction of change is predicted correctly in call cases. We analyzed the molecular dipole moments of the ligands in gas, water and protein environments. Neither binding affinity nor ligand solvation free energy in bulk water shows much dependence on the molecular dipole moments of the ligands. Substitution of the aromatic or the charged group in the ligand results in considerable change in the solvation energy in bulk water and protein whereas the binding affinity varies insignificantly due to cancellation. The effect of chemical modification on ligand charge distribution is mostly local. Replacing benzene with diazine has minimal impact on the atomic multipoles at the amidinium group. We have also utilized an implicit solvent based end-state approach to evaluate the binding free energies of these inhibitors. In this approach, the polarizable multipole model combined with Poisson-Boltzmann/surface area (PMPB/SA) provides the electrostatic interaction energy and the polar solvation free energy. Overall the relative binding free energies obtained from the PMPB/SA model are in good agreement with the experimental data

Biocontrol potential and action mechanism of Bacillus amyloliquefaciens DB2 on Bipolaris sorokiniana

IntroductionBipolaris sorokiniana is the popular pathogenic fungi fungus which lead to common root rot and leaf spot on wheat. Generally, chemical fungicides are used to control diseases. However, the environmental pollution resulting from fungicides should not be ignored. It is important to study the mode of antagonistic action between biocontrol microbes and plant pathogens to design efficient biocontrol strategies.ResultsAn antagonistic bacterium DB2 was isolated and identified as Bacillus amyloliquefaciens. The inhibition rate of cell-free culture filtrate (CF, 20%, v/v) of DB2 against B. sorokiniana reached 92.67%. Light microscopy and scanning electron microscopy (SEM) showed that the CF significantly altered the mycelial morphology of B. sorokiniana and disrupted cellular integrity. Fluorescence microscopy showed that culture filtrate destroyed mycelial cell membrane integrity, decreased the mitochondrial transmembrane potential, induced reactive oxygen species (ROS) accumulation, and nuclear damage which caused cell death in B. sorokiniana. Moreover, the strain exhibited considerable production of protease and amylase, and showed a significant siderophore and indole-3-acetic acid (IAA) production. In the detached leaves and potted plants control assay, B. amyloliquefacien DB2 had remarkable inhibition activity against B. sorokiniana and the pot control efficacy was 75.22%. Furthermore, DB2 suspension had a significant promotion for wheat seedlings growth.ConclusionB. amyloliquefaciens DB2 can be taken as a potential biocontrol agent to inhibit B. sorokiniana on wheat and promote wheat growth

A Model of a MAPK•Substrate Complex in an Active Conformation: A Computational and Experimental Approach

The mechanisms by which MAP kinases recognize and phosphorylate substrates are not completely understood. Efforts to understand the mechanisms have been compromised by the lack of MAPK-substrate structures. While MAPK-substrate docking is well established as a viable mechanism for bringing MAPKs and substrates into close proximity the molecular details of how such docking promotes phosphorylation is an unresolved issue. In the present study computer modeling approaches, with restraints derived from experimentally known interactions, were used to predict how the N-terminus of Ets-1 associates with ERK2. Interestingly, the N-terminus does not contain a consensus-docking site ((R/K)2-3-X2-6-ΦA-X-ΦB, where Φ is aliphatic hydrophobic) for ERK2. The modeling predicts that the N-terminus of Ets-1 makes important contributions to the stabilization of the complex, but remains largely disordered. The computer-generated model was used to guide mutagenesis experiments, which support the notion that Leu-11 and possibly Ile-13 and Ile-14 of Ets-1 1-138 (Ets) make contributions through binding to the hydrophobic groove of the ERK2 D-recruiting site (DRS). Based on the modeling, a consensus-docking site was introduced through the introduction of an arginine at residue 7, to give the consensus 7RK-X2-ΦA-X-ΦB13. This results in a 2-fold increase in kcat/Km for the phosphorylation of Ets by ERK2. Similarly, the substitution of the N-terminus for two different consensus docking sites derived from Elk-1 and MKK1 also improves kcat/Km by two-fold compared to Ets. Disruption of the N-terminal docking through deletion of residues 1-23 of Ets results in a 14-fold decrease in kcat/Km, with little apparent change in kcat. A peptide that binds to the DRS of ERK2 affects Km, but not kcat. Our kinetic analysis suggests that the unstructured N-terminus provides 10-fold uniform stabilization of the ground state ERK2•Ets•MgATP complex and intermediates of the enzymatic reaction