Guangju Chen

A comparison between artificial and natural water oxidation

Two artificial water oxidation catalysts, the blue dimer and the Llobet catalyst, have been studied using hybrid DFT methods. The results are compared to those for water oxidation in the natural photosystem II enzyme. Studies on the latter system have now reached a high level of understanding, at present much higher than the one for the artificial systems. A recent high resolution X-ray structural investigation of PSII has confirmed the main features of the structure of the oxygen evolving complex (OEC) suggested by previous DFT cluster studies. The O-O bond formation mechanism suggested is of direct coupling (DC) type between an oxygen radical and a bridging oxo ligand. A similar DC mechanism is found for the Llobet catalyst, while an acid-base (AB) mechanism is preferred for the blue dimer. All of them require at least one oxygen radical. Full energy diagrams, including both redox and chemical steps, have been constructed illustrating similarities and differences to the natural system. Unlike previous DFT studies, the results of the present study suggest that the blue dimer is rate-limited by the initial redox steps, and the Llobet catalyst by O(2) release. The results could be useful for further improvement of the artificial systems.

Structures and Stabilities of Pbn (n ≤ 20) Clusters

We have performed global structural optimizations for neutral lead clusters Pb(n) (n = 2-20) by using a genetic algorithm (GA) coupled with a tight-binding (TB) potential. The low-energy structures identified from a GA/TB search were further optimized at the DFT-PBE level. The calculated results show that the Pb(n) (14 < n </= 20) clusters favor compact spherical structures with hexagon and pentagon rings. These structures are different from those of Si(n), Ge(n), and Sn(n) clusters which favor prolates in the same size range. The binding energies, second differences in energy, and fragmentation behaviors of the Pb(n) clusters were also discussed. Pb(n) (n = 4, 7, 10, 13, 15, and 17) clusters are found to be special stable clusters, which is in good agreement with the experimental results.

Computational Determination of Binding Structures and Free Energies of Phosphodiesterase-2 with Benzo[1,4]diazepin-2-one Derivatives

Phosphodiesterase-2 (PDE2) is a key enzyme catalyzing hydrolysis of both cyclic adenosine monophosphate (cAMP) and cyclic guanosine monophosphate (cGMP) that serve as intracellular second messengers. PDE2 has been recognized as an attractive drug target, and selective inhibitors of PDE2 are expected to be promising candidates for the memory enhancer, antidepressant, and anxiolytic agent. In the present study, we examined the detailed binding structures and free energies for PDE2 interacting with a promising series of inhibitors, i.e., benzo[1,4]diazepin-2-one derivatives, by carrying out molecular docking, molecular dynamics (MD) simulations, binding free energy calculations, and binding energy decompositions. The computational results provide valuable insights into the detailed enzyme-inhibitor binding modes including important intermolecular interactions, e.g., the π-π stacking interactions with the common benzo[1,4]diazepin-2-one scaffold of the inhibitors, hydrogen bonding and hydrophobic interactions with the substituents on the benzo[1,4]diazepin-2-one scaffold. Future rational design of new, more potent inhibitors of PDE2 should carefully account for all of these favorable intermolecular interactions. By use of the MD-simulated binding structures, the calculated binding free energies are in good agreement with the experimental activity data for all of the examined benzo[1,4]diazepin-2-one derivatives. The enzyme-inhibitor binding modes determined and the agreement between the calculated and experimental results are expected to be valuable for future rational design of more potent inhibitors of PDE2.

Analysis of Ice-Binding Sites in Fish Type II Antifreeze Protein by Quantum Mechanics

Many organisms living in cold environments can survive subzero temperatures by producing antifreeze proteins (AFPs) or antifreeze glycoproteins. In this paper we investigate the ice-binding surface of type II AFP by quantum mechanical methods, which, to the best of our knowledge, represents the first time that molecular orbital computational approaches have been applied to AFPs. Molecular mechanical approaches, including molecular docking, energy minimization, and molecular dynamics simulation, were used to obtain optimal systems for subsequent quantum mechanical analysis. We selected 17 surface patches covering the entire surface of the type II AFP and evaluated the interaction energy between each of these patches and two different ice planes using semi-empirical quantum mechanical methods. We have demonstrated the weak orbital overlay phenomenon and the change of bond orders in ice. These results consistently indicate that a surface patch containing 19 residues (K37, L38, Y20, E22, Y21, I19, L57, T56, F53, M127, T128, F129, R17, C7, N6, P5, G10, Q1, and W11) is the most favorable ice-binding site for both a regular ice plane and an ice plane where water O atoms are randomly positioned. Furthermore, for the first time the computation results provide new insights into the weakening of the ice lattice upon AFP binding, which may well be a primary factor leading to AFP-induced ice growth inhibition.

Interaction identification of Zif268 and TATAZF proteins with GC-/AT-rich DNA sequence: A theoretical study

Molecular dynamics (MD) simulations for Zif268 (a zinc‐finger‐protein binding specifically to the GC‐rich DNA)‐d(A1G2C3G4T5G6G7G8C9A10C11)2 and TATAZF (a zinc‐finger‐protein recognizing the AT‐rich DNA)‐d(A1C2G3C4T5A6T7A8A9A10A11G12G13)2 complexes have been performed for investigating the DNA binding affinities and specific recognitions of zinc fingers to GC‐rich and AT‐rich DNA sequences. The binding free energies for the two systems have been further analyzed by using the molecular mechanics Poisson‐Boltzmann surface area (MM‐PBSA) method. The calculations of the binding free energies reveal that the affinity energy of Zif268‐DNA complex is larger than that of TATAZF‐DNA one. The affinity between the zinc‐finger‐protein and DNA is mainly driven by more favorable van‐der‐Waals and nonpolar/solvation interactions in both complexes. However, the affinity energy difference of the two binding systems is mainly caused by the difference of van‐der‐Waals interactions and entropy components. The decomposition analysis of MM‐PBSA free energies on each residue of the proteins predicts that the interactions between the residues with the positive charges and DNA favor the binding process; while the interactions between the residues with the negative charges and DNA behave in the opposite way. The interhydrogen‐bonds at the protein‐DNA interface and the induced intrafinger hydrogen bonds between the residues of protein for the Zif268‐DNA complex have been identified at some key contact sites. However, only the interhydrogen‐bonds between the residues of protein and DNA for TATAZF‐DNA complex have been found. The interactions of hydrogen‐bonds, electrostatistics and van‐der‐Waals type at some new contact sites have been identified. Moreover, the recognition characteristics of the two studied zinc‐finger‐proteins have also been discussed.

Theoretical investigation of the self-assembly of cyclo[(-β3-HGly)4-]

Abstract Cyclo[(-β3-HGly)4-] and its oligomers were studied using the density function B3LYP method. Their energetics and structural characteristics were calculated and analyzed. Significantly, we observed that the average interaction energy between two adjacent monomers in the β3-cyclopeptide oligomers showed marked increase upon addition of more monomers. However, the enhancement will gradually become weaker and eventually reach a plateau, giving rise to a stable nanotube. The same phenomenon was observed for their geometric structures and dipole moments of the assembly. Based on these new insights, we suggest that the synergetic effect resulted from addition of monomers will facilitate the enhancement of favorable interaction of the nanotube, which is the primary driving force of self-assembly.

Computational Study on the Function of Water within a β-Helix Antifreeze Protein Dimer and in the Process of Ice-Protein Binding

Antifreeze proteins (AFPs) help many organisms protect themselves from freezing in subzero temperatures. The most active AFPs found to date are those from insects, which possess exceptionally regular beta-helical structures. On the ice-binding surface of these proteins, regularly arrayed water molecules are observed within the repeating Thr-Xxx-Thr motif, but the exact role of these water molecules remains unknown. In this work, we have employed a number of computational methods to examine the role of these water molecules in an AFP from Tenebrio molitor (TmAFP). Our investigation involved a combination of molecular and quantum mechanical approaches. Properties such as stability, interaction energy, orbital overlap, and conformational analysis of various systems, including TmAFP-water, TmAFP-water-ice, and TmAFP-ice, were systematically evaluated and compared. The regularly arrayed water molecules were found to remain associated with TmAFP before ice binding, demonstrating that they are an intrinsic part of the protein. These water molecules may assist TmAFP in the process of ice recognition and binding. However, after facilitating the initial stages of ice recognition and binding, these water molecules are excluded in the final formation of the AFP-ice complex. The departure of these water molecules enables a better two-dimensional match between TmAFP and ice. These results agree with experimental observations showing that although these water molecules are aligned with the ice-binding hydroxyl groups of Thr residues in one dimension, they are in fact positioned slightly off in the second dimension, making a good two-dimensional match impossible.

Molecular dynamics simulation studies on the positive cooperativity of the Kemptide substrate with protein kinase A induced by the ATP ligand.

The positive cooperativity of the Kemptide substrate or the ATP molecule with the PKA catalytic subunit has been studied by dynamics simulations and free energy calculations on a series of binary and ternary models. The results revealed that the first ATP binding to the PKA catalytic subunit is energetically favorable for the successive Kemptide binding, confirming the positive cooperativity. The key residues Thr51, Glu170, and Phe187 in PKA contributing to the positive cooperativity have been found. The binding of ATP to PKA induces the positive cooperativity through one direct allosteric communication network in PKA from the ATP binding sites in the catalytic loop of the large lobe to the Kemptide binding sites in the activation segment of the large lobe, two indirect ones from those in the glycine-rich loop and the β3 strand of the small lobe, and from those in the catalytic loop to those in the activation segment via the αF helix media. The Tyr204Ala mutation in the activation segment of PKA causes both the decoupling of the cooperativity and the disruption of the corresponding allosteric network through the αF helix media.

Differences in conformational dynamics of [Pt3(HPTAB)]6+–DNA adducts with various cross-linking modes

We present here molecular dynamics simulations and DNA conformational dynamics for a series of trinuclear platinum [Pt3(HPTAB)]6+-DNA adducts [HPTAB = N,N,N′,N′,N′′,N′′-hexakis (2-pyridyl-methyl)-1,3,5-tris(aminomethyl) benzene], including three types of bifunctional crosslinks and four types of trifunctional crosslinks. Our simulation results reveal that binding of the trinuclear platinum compound to a DNA duplex induces the duplex unwinding in the vicinity of the platination sites, and causes the DNA to bend toward the major groove. As a consequence, this produces a DNA molecule whose minor groove is more widened and shallow compared to that of an undamaged bare-DNA molecule. Notably, for trifunctional crosslinks, we have observed extensive DNA conformational distortions, which is rarely seen for normal platinum–DNA adducts. Our findings, in this study, thus provide further support for the idea that platinum compounds with trifunctional intra-strand or long-range-inter-strand cross-linking modes can generate larger DNA conformational distortions than other types of cross-linking modes.

Theoretical study on stabilities of multiple hydrogen bonded dimers

A series of self-constituted multiple hydrogen bonded (MHB) complexes has been investigated systematically by density functional theory (PBE1PBE /6-31G**), the Morokuma energy decomposition method (HF/6-31G**) and MP2 (6-31G** and 6-311++G**) calculation. We have discovered that (i) for doubly hydrogen bonded (DHB) complexes, both the interaction energy and stability increase with the charge transfer energy; (ii) for quadruple hydrogen bonded (QHB) complexes, cooperativity is the most important factor determining stability of the complex: stronger cooperative energy correlates well with larger interaction energy and thus more stable complex and vice versa; (iii) correlation energy plays an important role in intermolecular interactions. The correlation energy, mainly consisting of dispersive energy, also exhibits cooperativity in MHB dimers: positive for M-aadd and generally negative for other complexes.