Exciton/Charge-transfer Electronic Couplings in Organic Semiconductors

Charge transfer (CT) states and excitons are important in energy conversion processes that occur in organic light emitting devices (OLEDS) and organic solar cells. An ab initio density functional theory (DFT) method for obtaining CT−exciton electronic couplings between CT states and excitons is presented. This method is applied to two organic heterodimers to obtain their CT−exciton coupling and adiabatic energy surfaces near their CT−exciton diabatic surface crossings. The results show that the new method provides a new window into the role of CT states in exciton−exciton transitions within organic semiconductors.United States. Dept. of Energy (DEFG02- 07ER46474)David & Lucile Packard Foundation (Fellowship

Functional Rotation of the Transporter AcrB: Insights into Drug Extrusion from Simulations

The tripartite complex AcrAB-TolC is the major efflux system in Escherichia coli. It extrudes a wide spectrum of noxious compounds out of the bacterium, including many antibiotics. Its active part, the homotrimeric transporter AcrB, is responsible for the selective binding of substrates and energy transduction. Based on available crystal structures and biochemical data, the transport of substrates by AcrB has been proposed to take place via a functional rotation, in which each monomer assumes a particular conformation. However, there is no molecular-level description of the conformational changes associated with the rotation and their connection to drug extrusion. To obtain insights thereon, we have performed extensive targeted molecular dynamics simulations mimicking the functional rotation of AcrB containing doxorubicin, one of the two substrates that were co-crystallized so far. The simulations, including almost half a million atoms, have been used to test several hypotheses concerning the structure-dynamics-function relationship of this transporter. Our results indicate that, upon induction of conformational changes, the substrate detaches from the binding pocket and approaches the gate to the central funnel. Furthermore, we provide strong evidence for the proposed peristaltic transport involving a zipper-like closure of the binding pocket, responsible for the displacement of the drug. A concerted opening of the channel between the binding pocket and the gate further favors the displacement of the drug. This microscopically well-funded information allows one to identify the role of specific amino acids during the transitions and to shed light on the functioning of AcrB

A Structural Model of the Pore-Forming Region of the Skeletal Muscle Ryanodine Receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca2+ over K+ in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca2+ in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K+ but not Ca2+. Our simulations with G4898R likewise show a decrease in the preference of Ca2+ over K+ in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level

On the Solvation of the Zn 2+

The solvation properties of the Zn2+ ion in methanol solution have been investigated using a combined approach based on molecular dynamics (MD) simulations and extended X-ray absorption fine structure (EXAFS) experimental results. The quantum mechanical potential energy surface for the interaction of the Zn2+ ion with a methanol molecule has been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective Zn-methanol interactions have been fitted by suitable analytical potentials, and have been utilized in the MD simulation to obtain the structural properties of the solution. The reliability of the whole procedure has been assessed by comparing the theoretical structural results with the EXAFS experimental data. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to the EXAFS experiments

The Structural Role of Mg2+ Ions in a Class I RNA Polymerase Ribozyme. A Molecular Simulation Study

According to the RNA world hypothesis, self-replicating ribozymes, storing the genetic information and being able to perform catalysis, were the constituents of the first living organisms. In particular, RNA polymerase ribozymes, similar to current proteinaceous enzymatic polymerases, may have been able to promote the synthesis of RNA strands in a primitive world. Polymerase catalysis is usually assisted by Mg2+ ions, but it is not always trivial to find out experimentally the number of Mg2+ ions placed in the active site as well as the identity and the number of their coordination ligands. Here, we addressed this issue in an artificial class I ligase ribozyme. On the basis of a recently solved crystal structure, we constructed computational models of reactant and product states of this ribozyme, considering monometallic and bimetallic species. Our models were relaxed by force field based molecular dynamics (MD) simulations and mixed quantum-classical (QM/MM) MD. The structural and dynamical properties of these models were consistent with experimental data and were validated by a comparison with the catalytic sites of proteinaceous DNA and RNA polymerases. Consistently with enzymatic polymerases, our results suggest that class I RNA ligases most probably contain two magnesium ions in the active site and they may, therefore, catalyze the junction of two RNA strands via "a two Mg2+ ions" mechanism

Mechanism of falcipain-2 inhibition by α,β-unsaturated benzo[1,4]diazepin-2-one methyl ester

Falcipain-2 (FP-2) is a papain-family cysteine protease of Plasmodium falciparum whose primary function is to degrade the host red cell hemoglobin, within the food vacuole, in order to provide free amino acids for parasite protein synthesis. Additionally it promotes host cell rupture by cleaving the skeletal proteins of the erythrocyte membrane. Therefore, the inhibition of FP-2 represents a promising target in the search of novel anti-malarial drugs. A potent FP-2 inhibitor, characterized by the presence in its structure of the 1,4-benzodiazepine scaffold and an a,b-unsaturated methyl ester moiety capable to react with the Cys42 thiol group located in the active site of FP-2, has been recently reported in literature. In order to study in depth the inhibition mechanism triggered by this interesting compound, we carried out, through ONIOM hybrid calculations, a computational investigation of the processes occurring when the inhibitor targets the enzyme and eventually leads to an irreversible covalent Michael adduct. Each step of the reaction mechanism has been accurately characterized and a detailed description of each possible intermediate and transition state along the pathway has been reported

Free Energies of Quantum Particles: The Coupled-Perturbed Quantum Umbrella Sampling Method

We introduce a new simulation method called Coupled-Perturbed Quantum Umbrella Sampling that extends the classical umbrella sampling approach to reaction coordinates involving quantum mechanical degrees of freedom. The central idea in our method is to solve coupled-perturbed equations to find the response of the quantum system's wave function along a reaction coordinate of interest. This allows for propagation of the system's dynamics under the influence of a quantum biasing umbrella potential and provides a method to rigorously undo the effects of the bias to compute equilibrium ensemble averages. In this way, one can drag electrons into regions of high free energy where they would otherwise not go, thus enabling chemistry by fiat. We demonstrate the applicability of our method for two condensed-phase systems of interest. First, we consider the interaction of a hydrated electron with an aqueous sodium cation, and we calculate a potential of mean force that shows that an e(-):Na(+) contact pair is the thermodynamically favored product starting from either a neutral sodium atom or the separate cation and electron species. Second, we present the first determination of a hydrated electron's free-energy profile relative to an air/water interface. For the particular model parameters used, we find that the hydrated electron is more thermodynamically stable in the bulk rather than at the interface. Our analysis suggests that the primary driving force keeping the electron away from the interface is the long-range electron-solvent polarization interaction rather than the short-range details of the chosen pseudopotential