Diphenylhexatriene membrane probes DPH and TMA-DPH: A comparative molecular dynamics simulation study

Fluorescence spectroscopy andmicroscopy have been utilized as tools in membrane biophysics for decades now. Because phospholipids are non-fluorescent, the use of extrinsic membrane probes in this context is commonplace. Among the latter, 1,6-diphenylhexatriene (DPH) and its trimethylammonium derivative (TMA-DPH) have been extensively used. It is widely believed that, owing to its additional charged group, TMA-DPH is anchored at the lipid/water interface and reports on a bilayer region that is distinct from that of the hydrophobic DPH. In this study, we employ atomistic MD simulations to characterize the behavior of DPH and TMA-DPH in 1-palmitoyl-2-oleoyl-sn-glycero-3-phosphocholine (POPC) and POPC/cholesterol (4:1) bilayers. We show that although the dynamics of TMA-DPH in thesemembranes is noticeably more hindered than that of DPH, the location of the average fluorophore of TMA-DPH is only ~3–4 Å more shallow than that of DPH. The hindrance observed in the translational and rotational motions of TMA-DPH compared to DPH is mainly not due to significant differences in depth, but to the favorable electrostatic interactions of the former with electronegative lipid atoms instead. By revealing detailed insights on the behavior of these two probes, our results are useful both in the interpretation of past work and in the planning of future experiments using themasmembrane reporters.

Computing the structural and vibrational properties of polymorphic organic molecular crystals through van der Waals corrected density functional theory and the electronic properties of organic thin films through microelectrostatic calculations.

The present Thesis reports on the various research projects to which I have contributed during my PhD period, working with several research groups, and whose results have been communicated in a number of scientific publications. The main focus of my research activity was to learn, test, exploit and extend the recently developed vdW-DFT (van der Waals corrected Density Functional Theory) methods for computing the structural, vibrational and electronic properties of ordered molecular crystals from first principles. A secondary, and more recent, research activity has been the analysis with microelectrostatic methods of Molecular Dynamics (MD) simulations of disordered molecular systems. While only very unreliable methods based on empirical models were practically usable until a few years ago, accurate calculations of the crystal energy are now possible, thanks to very fast modern computers and to the excellent performance of the best vdW-DFT methods. Accurate energies are particularly important for describing organic molecular solids, since they often exhibit several alternative crystal structures (polymorphs), with very different packing arrangements but very small energy differences. Standard DFT methods do not describe the long-range electron correlations which give rise to the vdW interactions. Although weak, these interactions are extremely sensitive to the packing arrangement, and neglecting them used to be a problem. The calculations of reliable crystal structures and vibrational frequencies has been made possible only recently, thanks to development of some good representations of the vdW contribution to the energy (known as “vdW corrections”).

Multiscale modeling of light-harvesting systems based on rhodamine B aggregates

Rhodamine B (RB) has been successfully exploited in the synthesis of light harvesting systems, but since RB is prone to form dimers acting as quenchers for the fluorescence, high energy transfer efficiencies can be reached only when using bulky and hydrophobic counterions acting as spacers between RBs. In this PhD thesis, a multiscale theoretical study aimed at providing insights into the structural, photophysical and optical properties of RB and its aggregates is presented. At the macroscopic level (no atomistic details) a phenomenological model describing the fluorescence decay of RB networks in presence of both quenching from dimers and exciton-exciton annihiliation is presented and analysed, showing that the quenching from dimers affects the decay only at long times, a feature that can be exploited in global fitting analysis to determine relevant chemical and photophysical information. At the mesoscopic level (atomistic details but no electronic structure) the RB aggregation in water in presence of different counterions is studied with molecular dynamics (MD) simulations. A new force field has been parametrized for describing the RB flexibility and the RB-RB interaction driving the dimerization. Simulations correctly predict the RB/counterion aggregation only in presence of bulky and hydrophobic counterion and its ability to prevent the dimerization. Finally, at the microscopic level, DFT calculations are performed to demonstrate the spacing action of bulky counterions, but standard TDDFT calculations are showed to fail in correctly describing the excited states of RB and its dimers. Moreover, also standard procedures proposed in literature for obtaining ad hoc functionals are showed to not work properly. A detailed analysis on the effect of the exact exchange shows that its short-range contribution is the crucial quantity for ameliorating results, and a new functional containing a proper amount of such an exchange is proposed and successfully tested.

Discrete (MD) and continuum (DD) simulations of void growth in single crystals.

Void growth in ductile materials is an important problem from the fundamental and technological viewpoint. Most of the models developed to quantify and understand the void growth process did not take into account two important factors: the anisotropic nature of plastic flow in single crystals and the size effects that appear when plastic flow is confined into very small regions.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Comparison between molecular dynamics abd Monte Carlo simulations of an ordering process in a binary alloy

Ordering in a binary alloy is studied by means of a molecular-dynamics (MD) algorithm which allows to reach the domain growth regime. Results are compared with Monte Carlo simulations using a realistic vacancy-atom (MC-VA) mechanism. At low temperatures fast growth with a dynamical exponent x>1/2 is found for MD and MC-VA. The study of a nonequilibrium ordering process with the two methods shows the importance of the nonhomogeneity of the excitations in the system for determining its macroscopic kinetics.

Modeling of human P450 oxidoreductase structure by in silico mutagenesis and MD simulation

P450 oxidoreductase (POR) is the obligate electron donor for microsomal cytochrome P450s and mutations in POR cause several metabolic disorders. We have modeled the structure of human P450 oxidoreductase by in silico amino acid replacements in the rat POR crystal structure. The rat POR has 94% homology with human POR and 38 amino acids were replaced to make its sequence identical to human POR. Several rounds of molecular dynamic simulations refined the model and removed structural clashes from side chain alterations of replaced amino acids. This approach has the advantage of keeping the cofactor contacts and structural features of the core enzyme intact which could not be achieved by homology based approaches. The final model from our approach was of high quality and compared well with experimentally determined structures of other PORs. This model will be used for analyzing the structural implications of mutations and polymorphisms in human POR.

Li & Pb MD Sinulations

Pb17Li is today a reference breeder material in diverse fusion R&D programs worldwide. Extracting dynamic and structural properties of liquid LiPb mixtures via molecular dynamics simulations, represent a crucial step for multiscale modeling efforts in order to understand the suitability of this compound for future Nuclear Fusion technologies. At present a Li-Pb cross potential is not available in the literature. Here we present our first results on the validation of two semi-empirical potentials for Li and Pb in liquid phase. Our results represent the establishment of a solid base as a previous but crucial step to implement a LiPb cross potential. Structural and thermodynamical analyses confirm that the implemented potentials for Li and Pb are realistic to simulate both elements in the liquid phase.

The deficit and the economy : an update of long-term simulations : report to the Chairman, Committee on the Budget, U.S. Senate, and the Chairman, Committee on the Budget, House of Representatives /

"B-259450"--P. 1.

Why are MD simulated protein folding times wrong?

The question of significant deviations of protein folding times simulated using molecular dynamics from experimental values is investigated. It is shown that in the framework of Markov State Model (MSM) describing the conformational dynamics of peptides and proteins, the folding time is very sensitive to the simulation model parameters, such as forcefield and temperature. Using two peptides as examples, we show that the deviations in the folding times can reach an order of magnitude for modest variations of the molecular model. We, therefore, conclude that the folding rate values obtained in molecular dynamics simulations have to be treated with care.

Data of the molecular dynamics simulations of mutations in the human connexin46 docking interface

The structure of hCx26 derived from the X-ray analysis was used to generate a homology model for hCx46. Interacting connexin molecules were used as starting model for the molecular dynamics (MD) simulation using NAMD and allowed us to predict the dynamic behavior of hCx46wt and the cataract related mutant hCx46N188T as well as two artificial mutants hCx46N188Q and hCx46N188D. Within the 50 ns simulation time the docked complex composed of the mutants dissociate while hCx46wt remains stable. The data indicates that one hCx46 molecule forms 5-7 hydrogen bonds (HBs) with the counterpart connexin of the opposing connexon. These HBs appear essential for a stable docking of the connexons as shown by the simulation of an entire gap junction channel and were lost for all the tested mutants. The data described here are related to the research article entitled "The cataract related mutation N188T in human connexin46 (hCx46) revealed a critical role for residue N188 in the docking process of gap junction channels" (Schadzek et al., 2015) [1].

INVESTIGATIONS INTO MECHANISMS UNDERLYING EXTREME WAVE FORMATIONS AND COMPUTATIONALLY INTENSIVE SIMULATIONS

Various mechanisms have been proposed to explain extreme waves or rogue waves in an oceanic environment including directional focusing, dispersive focusing, wave-current interaction, and nonlinear modulational instability. The Benjamin-Feir instability (nonlinear modulational instability), however, is considered to be one of the primary mechanisms for rogue-wave occurrence. The nonlinear Schrodinger equation is a well-established approximate model based on the same assumptions as required for the derivation of the Benjamin-Feir theory. Solutions of the nonlinear Schrodinger equation, including new rogue-wave type solutions are presented in the author's dissertation work. The solutions are obtained by using a predictive eigenvalue map based predictor-corrector procedure developed by the author. Features of the predictive map are explored and the influences of certain parameter variations are investigated. The solutions are rescaled to match the length scales of waves generated in a wave tank. Based on the information provided by the map and the details of physical scaling, a framework is developed that can serve as a basis for experimental investigations into a variety of extreme waves as well localizations in wave fields. To derive further fundamental insights into the complexity of extreme wave conditions, Smoothed Particle Hydrodynamics (SPH) simulations are carried out on an advanced Graphic Processing Unit (GPU) based parallel computational platform. Free surface gravity wave simulations have successfully characterized water-wave dispersion in the SPH model while demonstrating extreme energy focusing and wave growth in both linear and nonlinear regimes. A virtual wave tank is simulated wherein wave motions can be excited from either side. Focusing of several wave trains and isolated waves has been simulated. With properly chosen parameters, dispersion effects are observed causing a chirped wave train to focus and exhibit growth. By using the insights derived from the study of the nonlinear Schrodinger equation, modulational instability or self-focusing has been induced in a numerical wave tank and studied through several numerical simulations. Due to the inherent dissipative nature of SPH models, simulating persistent progressive waves can be problematic. This issue has been addressed and an observation-based solution has been provided. The efficacy of SPH in modeling wave focusing can be critical to further our understanding and predicting extreme wave phenomena through simulations. A deeper understanding of the mechanisms underlying extreme energy localization phenomena can help facilitate energy harnessing and serve as a basis to predict and mitigate the impact of energy focusing.

On The Consistency Of Monte Carlo Track Structure Dna Damage Simulations.

Monte Carlo track structures (MCTS) simulations have been recognized as useful tools for radiobiological modeling. However, the authors noticed several issues regarding the consistency of reported data. Therefore, in this work, they analyze the impact of various user defined parameters on simulated direct DNA damage yields. In addition, they draw attention to discrepancies in published literature in DNA strand break (SB) yields and selected methodologies. The MCTS code Geant4-DNA was used to compare radial dose profiles in a nanometer-scale region of interest (ROI) for photon sources of varying sizes and energies. Then, electron tracks of 0.28 keV-220 keV were superimposed on a geometric DNA model composed of 2.7 × 10(6) nucleosomes, and SBs were simulated according to four definitions based on energy deposits or energy transfers in DNA strand targets compared to a threshold energy ETH. The SB frequencies and complexities in nucleosomes as a function of incident electron energies were obtained. SBs were classified into higher order clusters such as single and double strand breaks (SSBs and DSBs) based on inter-SB distances and on the number of affected strands. Comparisons of different nonuniform dose distributions lacking charged particle equilibrium may lead to erroneous conclusions regarding the effect of energy on relative biological effectiveness. The energy transfer-based SB definitions give similar SB yields as the one based on energy deposit when ETH ≈ 10.79 eV, but deviate significantly for higher ETH values. Between 30 and 40 nucleosomes/Gy show at least one SB in the ROI. The number of nucleosomes that present a complex damage pattern of more than 2 SBs and the degree of complexity of the damage in these nucleosomes diminish as the incident electron energy increases. DNA damage classification into SSB and DSB is highly dependent on the definitions of these higher order structures and their implementations. The authors' show that, for the four studied models, different yields are expected by up to 54% for SSBs and by up to 32% for DSBs, as a function of the incident electrons energy and of the models being compared. MCTS simulations allow to compare direct DNA damage types and complexities induced by ionizing radiation. However, simulation results depend to a large degree on user-defined parameters, definitions, and algorithms such as: DNA model, dose distribution, SB definition, and the DNA damage clustering algorithm. These interdependencies should be well controlled during the simulations and explicitly reported when comparing results to experiments or calculations.