Li-Dong Gong

Development of a Polarizable Force Field Using Multiple Fluctuating Charges per Atom

A polarizable force field (PFF) using multiple fluctuating charges per atom, ABEEMσπ PFF, is presented in this work. The fluctuating partial charges are obtained from the electronegativity equalization principle applied to the decomposition scheme of atom-bond regions into multiple charge sites: atomic, lone-pair electron, and σ and π bond regions. These multiple partial charges per atom should better account for the polarization effect than single atomic charge in other PFFs. To evaluate the PFF, structural and energetic properties for some organic and biochemical systems, including rotational barriers; binding energies of base pairs; a base-base interaction in a B-DNA decamer; and interaction energies of ten stationary conformers of a water dimer, peptides, and bases with water molecules, have been calculated and compared with the experimental data or ab initio MP2 results. Molecular dynamics simulations using the PFF have been performed for crambin and BPTI protein systems. Better performances in modeling root-mean-square deviations of backbone bond lengths, bond angles, key dihedral angles, the coordinate root-mean-square shift of atoms, and the distribution of hydrogen bonds have been observed in comparison with other PFFs. These results indicate that the fluctuating charge force field, ABEEMσπ/MM, is accurate and reliable and can be applied to wide ranges of organic and biomolecular systems.

Development and applications of the ABEEM fluctuating charge molecular force field in the ion-containing systems

The microscopic mechanisms of ion hydration and ion selectivity in biomolecular systems are long-standing research topics, in which the difficulty is how to reasonably and accurately describe the ion-water and ion-biomolecule interactions. This paper summarizes the development and applications of the atom-bond electronegativity equalization fluctuating charge force field model, ABEEM/MM, in the investigations of ion hydration, metalloproteins and ion-DNA bases systems. Based on high-level quantum chemistry calculations, the parameters were optimized and the molecular potential functions were constructed and applied to studies of structures, activities, energetics, and thermodynamic and kinetic properties of these ion-containing systems. The results show that the performance of ABEEM/MM is generally better than that of the common force fields, and its accuracy can reach or approach that of the high-level ab initio MP2 method. These studies provide a solid basis for further investigations of ion selectivity in biomolecular systems, the structures and properties of metalloproteins and other related ion-containing systems.

Calculating solvation energies by means of a fluctuating charge model combined with continuum solvent model.

Continuum solvent models have shown to be very efficient for calculating solvation energy of biomolecules in solution. However, in order to produce accurate results, besides atomic radii or volumes, an appropriate set of partial charges of the molecule is needed. Here, a set of partial charges produced by a fluctuating charge model-the atom-bond electronegativity equalization method model (ABEEMσπ) fused into molecular mechanics is used to fit for the analytical continuum electrostatics model of generalized-Born calculations. Because the partial atomic charges provided by the ABEEMσπ model can well reflect the polarization effect of the solute induced by the continuum solvent in solution, accurate and rapid calculations of the solvation energies have been performed for series of compounds involving 105 small neutral molecules, twenty kinds of dipeptides and several protein fragments. The solvation energies of small neutral molecules computed with the combination of the GB model with the fluctuating charge protocol (ABEEMσπ∕GB) show remarkable agreement with the experimental results, with a correlation coefficient of 0.97, a slope of 0.95, and a bias of 0.34 kcal∕mol. Furthermore, for twenty kinds of dipeptides and several protein fragments, the results obtained from the analytical ABEEMσπ∕GB model calculations correlate well with those from ab initio and Poisson-Boltzmann calculations. The remarkable agreement between the solvation energies computed with the ABEEMσπ∕GB model and PB model provides strong motivation for the use of ABEEMσπ∕GB solvent model in the simulation of biochemical systems.

Investigation of base pairs containing oxidized guanine using ab initio method and ABEEMσπ polarizable force field.

The integrity of the genetic information is constantly threatened by oxidizing agents. Oxidized guanines have all been linked to different types of cancers. Theoretical approaches supplement the assorted experimental techniques, and bring new sight and opportunities to investigate the underlying microscopic mechanics. Unfortunately, there is no specific force field to DNA system including oxidized guanines. Taking high level ab initio calculations as benchmark, we developed the ABEEMσπ fluctuating charge force field, which uses multiple fluctuating charges per atom. And it was applied to study the energies, structures and mutations of base pairs containing oxidized guanines. The geometries were obtained in reference to other studies or using B3LYP/6-31+G* level optimization, which is more rational and timesaving among 24 quantum mechanical methods selected and tested by this work. The energies were determined at MP2/aug-cc-pVDZ level with BSSE corrections. Results show that the constructed potential function can accurately simulate the change of H-bond and the buckled angle formed by two base planes induced by oxidized guanine, and it provides reliable information of hydrogen bonding, stacking interaction and the mutation processes. The performance of ABEEMσπ polarizable force field in predicting the bond lengths, bond angles, dipole moments etc. is generally better than those of the common force fields. And the accuracy of ABEEMσπ PFF is close to that of the MP2 method. This shows that ABEEMσπ model is a reliable choice for further research of dynamics behavior of DNA fragment including oxidized guanine.

Development of the ABEEMσπ Polarization Force Field for Base Pairs with Amino Acid Residue Complexes

DNA damage caused by oxidized bases can lead to aging and cancer in living beings. Luckily, a repair enzyme is able to repair the oxidized bases. The key step is to accurately recognize the oxidized bases, which mainly rely on complex hydrogen bond interactions. We have calibrated the charge parameters and torsional parameters of the ABEEMσπ polarization force field (ABEEMσπ PFF) to accurately describe the intermolecular and intramolecular interactions. Taking the experiment and quantum chemical method as the benchmark, a series of properties of base pair-amino acid residue systems, DNA and DNA-protein interaction systems were calculated and compared with those of other force fields. We have done a tremendous amount of tasks in testing, calibrations, and analyses. The ABEEMσπ PFF not only explicitly gives the position and the partial charge of lone-pair electrons but also introduces a function kHB to fit special electrostatic interactions in hydrogen bond interaction regions. Therefore, it can accurately simulate the polarization effect and charge transfer of hydrogen bond interactions, especially for charged systems and sulfur-containing systems, such as the binding energy between amino acid and base pairs (24-28 kcal/mol), which is induced by charge transfer. The RMSD of ABEEMσπ PFF is 1.18 kcal/mol, whereas the RMSD of Amber OL15 is 8.21 kcal/mol. The relative positions of the amino acid residue have significantly changed, and the hydrogen bonds were broken when simulated by fixed charge force fields. In addition, owing to refitting the reasonable torsional parameters, the geometric structures optimized by ABEEMσπ PFF were well consistent with those of the M06-2X/6-311++G** method, but the simulations by fixed force fields have a large rotation of methyl and distortion of the plane of the base pair. After extensive MD simulation with four test DNAs and a DNA-protein system, we conclude that ABEEMσπ PFF shows better agreement when compared to experimental structures, which illustrates the reliability of our model and the transferability of the parameters.

Theoretical reflections on the structural polymorphism of the oxygen-evolving complex in the S2 state and the correlations to substrate water exchange and water oxidation mechanism in photosynthesis

The structural polymorphism of the oxygen-evolving complex is of great significance to photosynthetic water oxidation. Employing density functional theory calculations, we have made further advisement on the interconversion mechanism of O5 transfer in the S2 state, mainly focusing on the potentiality of multi-state reactivity and spin transitions. Then, O5 protonation is proven impossible in S2 for irreversibility of the interconversion, which serves as an auxiliary judgment for that in S1. Besides, the structural polymorphism could also be archived by alternative mechanisms involving Mn3 ligand exchange, one of which with Mn3(III) makes sense to substrate water exchange in S2, although being irresponsible for the derivations of the observed EPR signals. During the water exchange, high-spin states would prevail to facilitate electron transfer between the ferromagnetically coupled Mn centers. In addition, water exchange in S1 could account for the closed-cubane structure as the initial form entering S2 at cryogenic temperatures. With regard to water oxidation, the structural flexibility and variability in both S2 and S3 guarantee smooth W2-O5 coupling in S4, according to the substrate assignments from water exchange kinetics. Within this theoretical framework, the new XFEL findings on S1-S3 can be readily rationalized. Finally, an alternative mechanistic scenario for O-O bond formation with ·OH radical near O4 is presented, followed by water binding to the pivot Mn4(III) from O4 side during S4-S0. This may diversify the substrate sources combined with the Ca channel in water delivery for the forthcoming S-cycle.The structural polymorphism of the oxygen-evolving complex is of great significance to photosynthetic water oxidation. Employing density functional theory calculations, we have made further advisement on the interconversion mechanism of O5 transfer in the S2 state, mainly focusing on the potentiality of multi-state reactivity and spin transitions. Then, O5 protonation is proven impossible in S2 for irreversibility of the interconversion, which serves as an auxiliary judgment for the protonation state of O5 in S1. Besides, the structural polymorphism could also be archived by alternative mechanisms involving Mn3 ligand exchange, one of which with Mn3(III) makes sense to substrate water exchange in S2, although being irresponsible for the derivations of the observed EPR signals. During the water exchange, high-spin states would prevail to facilitate electron transfer between the ferromagnetically coupled Mn centers. In addition, water exchange in S1 could account for the closed-cubane structure as the initial form entering S2 at cryogenic temperatures. With regard to water oxidation, the structural flexibility and variability in both S2 and S3 guarantee smooth W2-O5 coupling in S4, according to the substrate assignments from water exchange kinetics. Within this theoretical framework, the new XFEL findings on S1-S3 can be readily rationalized. Finally, an alternative mechanistic scenario for OO bond formation with ·OH radical near O4 is presented, followed by water binding to the pivot Mn4(III) from O4 side during S4-S0. This may diversify the substrate sources combined with the Ca channel in water delivery for the forthcoming S-cycle.

Morphological transition difference of linear and cyclic block copolymer with polymer blending in a selective solvent by combining dissipative particle dynamics and all-atom molecular dynamics simulations based on the ABEEM polarizable force field

This study describes theoretical simulations of morphological transitions for linear and cyclic block copolymers with polymer blending. Mesoscopic dissipative particle dynamics simulations with reliable interaction parameters from all-atom molecular dynamics simulations based on the ABEEM polarizable force field, can reproduce the self-assembly behavior well and morphological transitions observed by experiments.