Simulation of High Density Lipoprotein Behavior on a Few Layer Graphene Undergoing Non-Uniform Mechanical Load.

Effect of a nonuniform external mechanical load on high density lipoprotein (HDL) in aqueous medium was investigated using course-grained molecular dynamics simulations. The nonuniform load was achieved by a few layer graphene on one side and closed single-walled carbon nanotube (SWNT) (7, 7) on the opposite side of lipoprotein. The tube had a diameter of 1 nm and was oriented perpendicularly to the graphene. HDL was located between them. The tube was approaching to HDL on graphene deforming it. We considered two cases of the tube movement with velocities of 20 and 5 m/s. Coarse-grained (CG) molecular dynamics with application of the MARTINI force field for HDL and coarse-grained model with an all-atom (AA)/CG mapping ratio of 1.5 for carbon nanotube (CNT) (each CG bead was modeled by the 4-site CG benzene) were used. Coarse-grained model of HDL was received by method of self-assembly. HDL was static but not fixed that gave the possibility to compensate its external influence in some way. It was established that in water medium HDL interacted with graphene substrate. It was established that in water HDL interacts with graphene substrate, slightly flattening but retaining its shape of the whole. It was also observed that during the calculations HDL partially dodged nanotube. Lipoprotein belts unfolded on the graphene substrate in the way of the best compensation for the impact of nanotubes. Finally, we observed that the approaching tube has passed through the less dense medium of dipalmitoylphosphatidylcholine (DPPC) and its pressure on the macromolecule decreased. Inhomogeneity of the external exposure deformed HDL at approximately 10-50%. The character of deformation demonstrated that lipoprotein has viscoelastic properties similar to a fluid. The discovered ability of lipoprotein may help to establish mechanism of interaction of lipoproteins with arterial walls and dynamic behavior of lipoproteins in arterial intima.

Donor/acceptor coupling shortcuts in electron transfer within ruthenium-modified derivatives of cytochrome b(562).

Quantitative theoretical studies of long-range electron transfer are still rare, and reliable computational methods to analyze these reactions are still being developed. We re-examined electron transfer reactions in ruthenium-modified cytochrome b562 derivatives focusing on accurate calculation of statistical average of electron transfer rates that are dominated by a small fraction of accessible protein conformations. We performed a series of ab initio calculations of donor/acceptor interactions over protein fragments sampled from long molecular dynamic trajectories and compared computed electron transfer rates to available experimental data. Our approach takes into account cofactor electronic structure and effects of solvation on the donor-acceptor interactions. It allows predicting absolute values of electron transfer rates in contrast to other computational methodologies that give only qualitative results. Our calculations reproduced with a good accuracy experimental electron transfer rates. We also found that electron transfer in some of the cytochrome b562 derivatives is dominated by "shortcut" conformations, where donor/acceptor interactions are mediated by nonbonded interactions of Ru ligands with protein surface groups. Several derivatives adopt long-lived conformations with the Ru complex interacting with negatively charged protein residues that are characterized by shorter Ru-Fe distances and higher ET rates. We argue that quantitative theoretical analysis is essential for detailed understanding of protein electron transfer and mechanisms of biological redox reactions.

Control of cytochrome c redox reactivity through off-pathway modifications in the protein hydrogen-bonding network.

Measurements of photoinduced Fe(2+)-to-Ru(3+) electron transfer (ET), supported by theoretical analysis, demonstrate that mutations off the dominant ET pathways can strongly influence the redox reactivity of cytochrome c. The effects arise from the change in the protein dynamics mediated by the intraprotein hydrogen-bonding network.

Modeling DNA-bending in the nucleosome: role of AA periodicity.

This paper uses atomistic molecular mechanics within the framework of the JUMNA model to study the bending properties of DNA segments, with emphasis on understanding the role of the 10 bp periodicity associated with AA repeats that has been found to dominate in nucleosomal DNA. The calculations impose a bending potential on 18 bp segments that is consistent with nucleosome structures (i.e., radius of curvature of 4.1 nm), and then determine the energies of the minimum energy structures for different values of the rotational register (a measure of the direction of bending of the DNA) subject to forces derived from the Amber force field (parm99bsc0). The results show that sequences that contain the 10 bp repeats but are otherwise random have a narrow distribution of rotational register values that minimize the energy such that it is possible to combine several minimized structures to give the 147 bp nearly planar loop structure of the nucleosome. The rotational register values that lead to minimum bending energy with 10 bp AA repeats have a narrower minor groove, which points toward the histone interior at the positions of the AA repeats, which is a result that matches the experiments. The calculations also show that these sequences have a relatively flat potential energy landscape for bending to a 4.1 nm radius of curvature. Random sequences that do not have the 10 bp AA repeats have less stable bent structures, and a flat rotational register distribution, such that low energy nearly planar loops are less likely.

Cooperative melting in caged dimers with only two DNA duplexes.

Small molecule-DNA hybrids with only two parallel DNA duplexes (rSMDH2) displayed sharper melting profiles compared to unmodified DNA duplexes, consistent with predictions from neighboring-duplex theory. Using adjusted thermodynamic parameters obtained from a coarse-grain dynamic simulation, the experimental data fit well to an analytical model.

Using DNA to Link Gold Nanoparticles, Polymers and Molecules: a Theoretical Perspective.

This Perspective describes theoretical studies aimed at understanding the structural and thermal properties of materials in which DNA is used to link gold nanoparticles, or polymers or organic molecules. Particularly in the case of gold nanoparticles, the materials derived from this structural motif have proven to be important for biological sensing and other applications, however additional applications may arise as a result of recent advances in the preparation of crystalline materials based on DNA-linked particles. From a theory perspective these are challenging materials to describe due to the large number of atoms, and the polyelectrolyte character of DNA, however there has been important progress recently using all-atom and coarse-grained molecular dynamics, and with analytical theory. Among topics that we discuss are the structure and density of DNA when attached to gold particles, and the size and melting properties of DNA-linked nanoparticles in different environments.

DNA melting in small-molecule-DNA-hybrid dimer structures: experimental characterization and coarse-grained molecular dynamics simulations.

When DNA hybridization is used to link together nanoparticles or molecules, the melting transition of the resulting DNA-linked material often is very sharp. In this paper, we study a particularly simple version of this class of material based on a small-molecule-DNA-hybrid (SMDH) structure that has three DNA strands per 1,3,5-tris(phenylethynyl)benzene core. By varying the concentration of the SMDHs, it is possible to produce either SMDH dimers or bulk aggregates, with the former having highly packed duplex DNA while the latter has an extended network. Melting measurements that we present show that the dimers exhibit sharp melting while the extended aggregates show broad melting. To interpret these results, we have performed coarse-grained molecular dynamics (CGMD) studies of the dimer melting and also of isolated duplex melting using CGMD potentials that have either implicit or explicit ions. Details of the melting simulation technology demonstrate that the simulations properly describe equilibrium transitions in isolated duplexes. The results show that the SMDH dimer has much sharper melting than the isolated duplex. Both implicit and explicit ion calculations show this effect, but the explicit ion results are sharper. An analytical model of the melting thermodynamics is developed which shows that the sharp melting is entropically driven and can be understood primarily in terms of the differences between the effective concentrations of the DNA strands for intracomplex hybridization events compared to intermolecular hybridization.

Coupling coherence distinguishes structure sensitivity in protein electron transfer.

Quantum mechanical analysis of electron tunneling in nine thermally fluctuating cytochrome b562 derivatives reveals two distinct protein-mediated coupling limits. A structure-insensitive regime arises for redox partners coupled through dynamically averaged multiple-coupling pathways (in seven of the nine derivatives) where heme-edge coupling leads to the multiple-pathway regime. A structure-dependent limit governs redox partners coupled through a dominant pathway (in two of the nine derivatives) where axial-ligand coupling generates the single-pathway limit and slower rates. This two-regime paradigm provides a unified description of electron transfer rates in 26 ruthenium-modified heme and blue-copper proteins, as well as in numerous photosynthetic proteins.

Photoselected electron transfer pathways in DNA photolyase.

Cyclobutane dimer photolyases are proteins that bind to UV-damaged DNA containing cyclobutane pyrimidine dimer lesions. They repair these lesions by photo-induced electron transfer. The electron donor cofactor of a photolyase is a two-electron-reduced flavin adenine dinucleotide (FADH(-)). When FADH(-) is photo-excited, it transfers an electron from an excited pi --> pi* singlet state to the pyrimidine dimer lesion of DNA. We compute the lowest excited singlet states of FADH(-) using ab initio (time-dependent density functional theory and time-dependent Hartree-Fock), and semiempirical (INDO/S configuration interaction) methods. The calculations show that the two lowest pi --> pi* singlet states of FADH(-) are localized on the side of the flavin ring that is proximal to the dimer lesion of DNA. For the lowest-energy donor excited state of FADH(-), we compute the conformationally averaged electronic coupling to acceptor states of the thymine dimer. The coupling calculations are performed at the INDO/S level, on donor-acceptor cofactor conformations obtained from molecular dynamics simulations of the solvated protein with a thymine dimer docked in its active site. These calculations demonstrate that the localization of the (1)FADH(-)* donor state on the flavin ring enhances the electronic coupling between the flavin and the dimer by permitting shorter electron-transfer pathways to the dimer that have single through-space jumps. Therefore, in photolyase, the photo-excitation itself enhances the electron transfer rate by moving the electron towards the dimer.

Ab initio based calculations of electron-transfer rates in metalloproteins.

A long-standing challenge in electron-transfer theory is to compute accurate rates of long-distance reactions in proteins. We describe an ab initio Hartree-Fock approach to compute electronic-coupling interactions and electron-transfer rates in proteins that allows the favorable comparison with experiment. The method includes the following key features; each is essential for reliable rate computations: (1) summing contributions over multiple tunneling pathways, (2) averaging couplings over thermally accessible protein conformations, (3) describing donor and acceptor electronic structure explicitly, including solvation effects, and averaging coupling over multiple energy-level crossings of the nearly degenerate donor-acceptor ligand-field states, and (4) eliminating basis set artifacts associated with diffuse basis functions. The strong dependence of coupling on donor-acceptor distance and on pathway interferences causes large variations of the computed electron-coupling values with protein geometry, and the strongest coupled conformers dominate the electron-transfer rate. As such, averaging over thermally accessible conformers of the protein and of the redox cofactors is essential. This approach was tested on six ruthenium-modified azurin derivatives using the high temperature nonadiabatic rate expression and compared with simpler pathways, average barrier, and semiempirical INDO models. Results of ab initio Hartree-Fock calculations with a split-valence basis set are in good agreement with the experimental rates. Predicted rates in the longer-distance derivatives are underestimated by 3-8-fold. This analysis indicates that the key ingredients needed for quantitatively reliable protein electron-transfer rate calculations are accessible.