CHARMM general force field: A force field for drug-like molecules compatible with the CHARMM all-atom additive biological force fields

Potential Role of the Last Half Repeat in TAL Effectors Revealed by a Molecular Simulation Study

TAL effectors (TALEs) contain a modular DNA-binding domain that is composed of tandem repeats. In all naturally occurring TALEs, the end of tandem repeats is invariantly a truncated half repeat. To investigate the potential role of the last half repeat in TALEs, we performed comparative molecular dynamics simulations for the crystal structure of DNA-bound TALE AvrBs3 lacking the last half repeat and its modeled structure having the last half repeat. The structural stability analysis indicates that the modeled system is more stable than the nonmodeled system. Based on the principle component analysis, it is found that the AvrBs3 increases its structural compactness in the presence of the last half repeat. The comparison of DNA groove parameters of the two systems implies that the last half repeat also causes the change of DNA major groove binding efficiency. The following calculation of hydrogen bond reveals that, by stabilizing the phosphate binding with DNA at the C-terminus, the last half repeat helps to adopt a compact conformation at the protein-DNA interface. It further mediates more contacts between TAL repeats and DNA nucleotide bases. Finally, we suggest that the last half repeat is required for the high-efficient recognition of DNA by TALE.

Predicting Adsorption Affinities of Small Molecules on Carbon Nanotubes Using Molecular Dynamics Simulation

Computational techniques have the potential to accelerate the design and optimization of nanomaterials for applications such as drug delivery and contaminant removal; however, the success of such techniques requires reliable models of nanomaterial surfaces as well as accurate descriptions of their interactions with relevant solutes. In the present work, we evaluate the ability of selected models of naked and hydroxylated carbon nanotubes to predict adsorption equilibrium constants for about 30 small aromatic compounds with a variety of functional groups. The equilibrium constants determined using molecular dynamics coupled with free-energy calculation techniques are directly compared to those derived from experimental measurements. The calculations are highly predictive of the relative adsorption affinities of the compounds, with excellent correlation (r ≥ 0.9) between calculated and measured values of the logarithm of the adsorption equilibrium constant. Moreover, the agreement in absolute terms is also reasonable, with average errors of less than one decade. We also explore possible effects of surface loading, although we demonstrate that they are negligible for the experimental conditions considered. Given the degree of reliability demonstrated, we move on to employing the in silico techniques in the design of nanomaterials, using the optimization of adsorption affinity for the herbacide atrazine as an example. Our simulations suggest that, compared to other modifications of graphenic carbon, polyvinylpyrrolidone conjugation gives the highest affinity for atrazine-substantially greater than that of graphenic carbon alone-and may be useful as a nanomaterial for delivery or sequestration of atrazine.

Aspheric Solute Ions Modulate Gold Nanoparticle Interactions in an Aqueous Solution: An Optimal Way To Reversibly Concentrate Functionalized Nanoparticles

Nanometer-sized gold particles (AuNPs) are of peculiar interest because their behaviors in an aqueous solution are sensitive to changes in environmental factors including the size and shape of the solute ions. In order to determine these important characteristics, we performed all-atom molecular dynamics simulations on the icosahedral Au144 nanoparticles each coated with a homogeneous set of 60 thiolates (4-mercaptobenzoate, pMBA) in eight aqueous solutions having ions of varying sizes and shapes (Na(+), K(+), tetramethylamonium cation TMA(+), tris-ammonium cation TRS(+), Cl(-), and OH(-)). For each solution, we computed the reversible work (potential of mean of force) to bring two nanoparticles together as a function of their separation distance. We found that the behavior of pMBA protected Au144 nanoparticles can be readily modulated by tuning their aqueous environmental factors (pH and solute ion combinations). We examined the atomistic details on how the sizes and shapes of solute ions quantitatively factor in the definitive characteristics of nanoparticle-environment and nanoparticle-nanoparticle interactions. We predict that tuning the concentrations of nonspherical composite ions such as TRS(+) in an aqueous solution of AuNPs be an effective means to modulate the aggregation propensity desired in biomedical and other applications of small charged nanoparticles.

Quantum-Classical Path Integral Simulation of Ferrocene-Ferrocenium Charge Transfer in Liquid Hexane

We employ the quantum-classical path integral methodology to simulate the outer sphere charge-transfer process of the ferrocene-ferrocenium pair in liquid hexane with unprecedented accuracy. Comparison of the simulation results to those obtained by mapping the solvent on an effective harmonic bath demonstrates the accuracy of linear response theory in this system.

Single-Point Mutation with a Rotamer Library Toolkit: Toward Protein Engineering

Protein engineers have long been hard at work to harness biocatalysts as a natural source of regio-, stereo-, and chemoselectivity in order to carry out chemistry (reactions and/or substrates) not previously achieved with these enzymes. The extreme labor demands and exponential number of mutation combinations have induced computational advances in this domain. The first step in our virtual approach is to predict the correct conformations upon mutation of residues (i.e., rebuilding side chains). For this purpose, we opted for a combination of molecular mechanics and statistical data. In this work, we have developed automated computational tools to extract protein structural information and created conformational libraries for each amino acid dependent on a variable number of parameters (e.g., resolution, flexibility, secondary structure). We have also developed the necessary tool to apply the mutation and optimize the conformation accordingly. For side-chain conformation prediction, we obtained overall average root-mean-square deviations (RMSDs) of 0.91 and 1.01 Å for the 18 flexible natural amino acids within two distinct sets of over 3000 and 1500 side-chain residues, respectively. The commonly used dihedral angle differences were also evaluated and performed worse than the state of the art. These two metrics are also compared. Furthermore, we generated a family-specific library for kinases that produced an average 2% lower RMSD upon side-chain reconstruction and a residue-specific library that yielded a 17% improvement. Ultimately, since our protein engineering outlook involves using our docking software, Fitted/Impacts, we applied our mutation protocol to a benchmarked data set for self- and cross-docking. Our side-chain reconstruction does not hinder our docking software, demonstrating differences in pose prediction accuracy of approximately 2% (RMSD cutoff metric) for a set of over 200 protein/ligand structures. Similarly, when docking to a set of over 100 kinases, side-chain reconstruction (using both general and biased conformation libraries) had minimal detriment to the docking accuracy.

Computation of Hydration Free Energies Using the Multiple Environment Single System Quantum Mechanical/Molecular Mechanical Method

A recently developed MESS-E-QM/MM method (multiple-environment single-system quantum mechanical molecular/mechanical calculations with a Roothaan-step extrapolation) is applied to the computation of hydration free energies for the blind SAMPL4 test set and for 12 small molecules. First, free energy simulations are performed with a classical molecular mechanics force field using fixed-geometry solute molecules and explicit TIP3P solvent, and then the non-Boltzmann-Bennett method is employed to compute the QM/MM correction (QM/MM-NBB) to the molecular mechanical hydration free energies. For the SAMPL4 set, MESS-E-QM/MM-NBB corrections to the hydration free energy can be obtained 2 or 3 orders of magnitude faster than fully converged QM/MM-NBB corrections, and, on average, the hydration free energies predicted with MESS-E-QM/MM-NBB fall within 0.10-0.20 kcal/mol of full-converged QM/MM-NBB results. Out of five density functionals (BLYP, B3LYP, PBE0, M06-2X, and ωB97X-D), the BLYP functional is found to be most compatible with the TIP3P solvent model and yields the most accurate hydration free energies against experimental values for solute molecules included in this study.

Molecular Dynamics Simulations of the Initial Adsorption Stages of Fibrinogen on Mica and Graphite Surfaces

Fibrinogen, a blood glycoprotein of vertebrates, plays an essential role in blood clotting by polymerizing into fibrin when activated. Upon adsorption on material surfaces, it also contributes to determine their biocompatibility and has been implicated in the onset of thrombosis and inflammation at medical implants. Here we present the first fully atomistic simulations of the initial stages of the adsorption process of fibrinogen on mica and graphite surfaces. The simulations reveal a weak adsorption on mica that allows frequent desorption and reorientation events. This adsorption is driven by electrostatic interactions between the protein and the silicate surface as well as the counterion layer. Preferred adsorption orientations for the globular regions of the protein are identified. The adsorption on graphite is found to be stronger with fewer reorientation and desorption events and shows the onset of denaturation of the protein.

Fighting Cancer with Transition Metal Complexes: From Naked DNA to Protein and Chromatin Targeting Strategies

Many transition metal complexes have unique physicochemical properties that can be efficiently exploited in medicinal chemistry for cancer treatment. Traditionally, double-stranded DNA has been assumed to be the main binding target; however, recent studies have shown that nucleosomal DNA as well as proteins can act as dominant molecular binding partners. This has raised new questions about the molecular determinants that govern DNA versus protein binding selectivity, and has offered new ways to rationalize their biological activity and possible side effects. To address these questions, molecular simulations at an atomistic level of detail have been used to complement, support, and rationalize experimental data. Herein we review some relevant studies-focused on platinum and ruthenium compounds-to illustrate the power of state-of-the-art molecular simulation techniques and to demonstrate how the interplay between molecular simulations and experiments can make important contributions to elucidating the target preferences of some promising transition metal anticancer agents. This contribution aims at providing relevant information that may help in the rational design of novel drug-discovery strategies.

Molecular Dynamics of Membrane-Spanning DNA Channels: Conductance Mechanism, Electro-Osmotic Transport, and Mechanical Gating

DNA self-assembly has emerged as a new paradigm for design of biomimetic membrane channels. Several experimental groups have already demonstrated assembly and insertion of DNA channels into lipid bilayer membranes; however, the structure of the channels and their conductance mechanism have remained undetermined. Here, we report the results of molecular dynamics simulations that characterized the biophysical properties of the DNA membrane channels with atomic precision. We show that, while overall remaining stable, the local structure of the channels undergoes considerable fluctuations, departing from the idealized design. The transmembrane ionic current flows both through the central pore of the channel as well as along the DNA walls and through the gaps in the DNA structure. Surprisingly, we find that the conductance of DNA channels depend on the membrane tension, making them potentially suitable for force-sensing applications. Finally, we show that electro-osmosis governs the transport of druglike molecules through the DNA channels.

Use of Nonequilibrium Work Methods to Compute Free Energy Differences Between Molecular Mechanical and Quantum Mechanical Representations of Molecular Systems

Carrying out free energy simulations (FES) using quantum mechanical (QM) Hamiltonians remains an attractive, albeit elusive goal. Renewed efforts in this area have focused on using "indirect" thermodynamic cycles to connect "low level" simulation results to "high level" free energies. The main obstacle to computing converged free energy results between molecular mechanical (MM) and QM (ΔA(MM→QM)), as recently demonstrated by us and others, is differences in the so-called "stiff" degrees of freedom (e.g., bond stretching) between the respective energy surfaces. Herein, we demonstrate that this problem can be efficiently circumvented using nonequilibrium work (NEW) techniques, i.e., Jarzynski's and Crooks' equations. Initial applications of computing ΔA(NEW)(MM→QM), for blocked amino acids alanine and serine as well as to generate butane's potentials of mean force via the indirect QM/MM FES method, showed marked improvement over traditional FES approaches.

CHARMM-GUI Input Generator for NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM Simulations Using the CHARMM36 Additive Force Field

Proper treatment of nonbonded interactions is essential for the accuracy of molecular dynamics (MD) simulations, especially in studies of lipid bilayers. The use of the CHARMM36 force field (C36 FF) in different MD simulation programs can result in disagreements with published simulations performed with CHARMM due to differences in the protocols used to treat the long-range and 1-4 nonbonded interactions. In this study, we systematically test the use of the C36 lipid FF in NAMD, GROMACS, AMBER, OpenMM, and CHARMM/OpenMM. A wide range of Lennard-Jones (LJ) cutoff schemes and integrator algorithms were tested to find the optimal simulation protocol to best match bilayer properties of six lipids with varying acyl chain saturation and head groups. MD simulations of a 1,2-dipalmitoyl-sn-phosphatidylcholine (DPPC) bilayer were used to obtain the optimal protocol for each program. MD simulations with all programs were found to reasonably match the DPPC bilayer properties (surface area per lipid, chain order parameters, and area compressibility modulus) obtained using the standard protocol used in CHARMM as well as from experiments. The optimal simulation protocol was then applied to the other five lipid simulations and resulted in excellent agreement between results from most simulation programs as well as with experimental data. AMBER compared least favorably with the expected membrane properties, which appears to be due to its use of the hard-truncation in the LJ potential versus a force-based switching function used to smooth the LJ potential as it approaches the cutoff distance. The optimal simulation protocol for each program has been implemented in CHARMM-GUI. This protocol is expected to be applicable to the remainder of the additive C36 FF including the proteins, nucleic acids, carbohydrates, and small molecules.

Study of the Differential Activity of Thrombin Inhibitors Using Docking, QSAR, Molecular Dynamics, and MM-GBSA

Non-peptidic thrombin inhibitors (TIs; 177 compounds) with diverse groups at motifs P1 (such as oxyguanidine, amidinohydrazone, amidine, amidinopiperidine), P2 (such as cyanofluorophenylacetamide, 2-(2-chloro-6-fluorophenyl)acetamide), and P3 (such as phenylethyl, arylsulfonate groups) were studied using molecular modeling to analyze their interactions with S1, S2, and S3 subsites of the thrombin binding site. Firstly, a protocol combining docking and three dimensional quantitative structure-activity relationship was performed. We described the orientations and preferred active conformations of the studied inhibitors, and derived a predictive CoMSIA model including steric, donor hydrogen bond, and acceptor hydrogen bond fields. Secondly, the dynamic behaviors of some selected TIs (compounds 26, 133, 147, 149, 162, and 177 in this manuscript) that contain different molecular features and different activities were analyzed by creating the solvated models and using molecular dynamics (MD) simulations. We used the conformational structures derived from MD to accomplish binding free energetic calculations using MM-GBSA. With this analysis, we theorized about the effect of van der Waals contacts, electrostatic interactions and solvation in the potency of TIs. In general, the contents reported in this article help to understand the physical and chemical characteristics of thrombin-inhibitor complexes.

Trimodal Therapy: Combining Hyperthermia with Repurposed Bexarotene and Ultrasound for Treating Liver Cancer

Repurposing of existing cancer drugs to overcome their physical limitations, such as insolubility, represents an attractive strategy to achieve enhanced therapeutic efficacy and broaden the range of clinical applications. Such an approach also promises to offer substantial cost savings in drug development efforts. Here we repurposed FDA-approved topical agent bexarotene (Targretin), currently in limited use for cutaneous manifestations of T-cell lymphomas, and re-engineer it for use in solid tumor applications by forming self-assembling nanobubbles. Physico-chemical characterization studies of the novel prodrug nanobubbles demonstrated their stability, enhanced target cell internalization capability, and highly controlled release profile in response to application of focused ultrasound energy. Using an in vitro model of hepatocellular carcinoma and an in vivo large animal model of liver ablation, we demonstrate the effectiveness of bexarotene prodrug nanobubbles when used in conjunction with catheter-based ultrasound, thereby highlighting the therapeutic promise of this trimodal approach.

Supramolecular Enhancement of Protein Analysis via the Recognition of Phenylalanine with Cucurbit[7]uril

Mass spectrometry (MS)-based analysis using enzymatic digestion is widely used for protein sequencing and characterization. The large number of peptides generated from proteolysis, however, suppresses the signal of peptides with low ionization efficiency, thus precluding their observation and analysis. This study describes a technique for improved analysis of peptic peptides by adding the synthetic receptor cucurbit[7]uril (CB[7]), which binds selectively to peptides with N-terminal aromatic residues. Capturing the N-terminal phenylalanine (Phe) of peptides using CB[7] enhances the peptide abundances both in electrospray ionization MS and in matrix-assisted laser desorption ionization MS. Moreover, collision-induced dissociation (CID) of the CB[7]·peptide complex ions generates b- and y-type fragment ions with higher sequence coverage than those generated with uncomplexed peptides. The signal enhancement mediated by CB[7] is attributed to an increase in the peptide proton affinities upon CB[7] complexation. The mechanistic details of the fragmentation process are discussed on the basis of the structures of the complex ions obtained from ion mobility (IM) measurements and molecular modeling. This study demonstrates a novel and powerful approach to the enhancement of protein and peptide analysis using a synthetic receptor, without the need for new instrumentation, chemical modifications, or specialized sample preparation. The simplicity and potential generality of this technique should provide a valuable asset in the toolbox of routine protein and peptide analysis.

Structural Degradation and Swelling of Lipid Bilayer under the Action of Benzene

Benzene and other nonpolar organic solvents can accumulate in the lipid bilayer of cellular membranes. Their effect on the membrane structure and fluidity determines their toxic properties and antibiotic action of the organic solvents on the bacteria. We performed molecular dynamics simulations of the interaction of benzene with the dimyristoylphosphatidylcholine (DMPC) bilayer. An increase in the membrane surface area and fluidity was clearly detected. Changes in the acyl chain ordering, tilt angle, and overall bilayer thickness were, however, much less marked. The dependence of all computed quantities on the benzene content showed two regimes separated by the solubility limit of benzene in water. When the amount of benzene exceeded this point, a layer of almost pure benzene started to grow between the membrane leaflets. This process corresponds to the nucleation of a new phase and provides a molecular mechanism for the mechanical rupture of the bilayer under the action of nonpolar compounds.

Hydroxymethyluracil modifications enhance the flexibility and hydrophilicity of double-stranded DNA

Oxidation of a DNA thymine to 5-hydroxymethyluracil is one of several recently discovered epigenetic modifications. Here, we report the results of nanopore translocation experiments and molecular dynamics simulations that provide insight into the impact of this modification on the structure and dynamics of DNA. When transported through ultrathin solid-state nanopores, short DNA fragments containing thymine modifications were found to exhibit distinct, reproducible features in their transport characteristics that differentiate them from unmodified molecules. Molecular dynamics simulations suggest that 5-hydroxymethyluracil alters the flexibility and hydrophilicity of the DNA molecules, which may account for the differences observed in our nanopore translocation experiments. The altered physico-chemical properties of DNA produced by the thymine modifications may have implications for recognition and processing of such modifications by regulatory DNA-binding proteins.

A Delicate Balance When Substituting a Small Hydrophobe onto Low Molecular Weight Polyethylenimine to Improve Its Nucleic Acid Delivery Efficiency

High molecular weight (HMW) polyethylenimine (PEI) is one of the most versatile nonviral gene vectors that was extensively investigated over the past two decades. The cytotoxic profile of HMW PEI, however, encouraged a search for safer alternatives. Because of lack of cytotoxicity of low molecular weight (LMW) PEI, enhancing its performance via hydrophobic modifications has been pursued to this end. Since the performance of modified PEIs depends on the nature and extent of substituents, we systematically investigated the effect of hydrophobic modification of LMW (1.2 kDa) PEI with a short propionic acid (PrA). Moderate enhancements in PEI hydrophobicity resulted in enhanced cellular uptake of polyplexes and siRNA-induced silencing efficacy, whereas further increase in PrA substitution abolished the uptake as well as the silencing. We performed all-atom molecular dynamics simulations to elucidate the mechanistic details behind these observations. A new assembly mechanism was observed by the presence of hydrophobic PrA moieties, where PrA migrated to core of the polyplex. This phenomenon caused higher surface hydrophobicity and surface charge density at low substitutions, and it caused deleterious effects on surface hydrophobicity and cationic charge at higher substitutions. It is evident that an optimal balance of hydrophobicity/hydrophilicity is needed to achieve the desired polyplex properties for an efficient siRNA delivery, and our mechanistic findings should provide valuable insights for the design of improved substituents on nonviral carriers.

The Role of Linkers in the Excited-State Dynamic Planarization Processes of Macrocyclic Oligothiophene 12-Mers

Linkers adjoining chromophores play an important role in modulating the structure of conjugated systems, which is bound up with their photophysical properties. However, to date, the focus of works dealing with linker effects was limited only to linear π-conjugated materials, and there have been no detailed studies on cyclic counterparts. Herein we report the linker effects on the dynamic planarization processes of π-conjugated macrocyclic oligothiophene 12-mers, where the different ratio between ethynylene and vinylene linkers was chosen to control the backbone rigidity. By analyzing transient fluorescence spectra, we demonstrate that the connecting linkers play a crucial role in the excited-state dynamics of cyclic conjugated systems. Faster dynamic planarization, longer exciton delocalization length, and higher degree of planarity were observed in vinylene inserted cyclic oligothiophenes. Molecular dynamics simulations and density functional theory calculations also stress the importance of the role of linkers in modulating the structure of cyclic oligothiophenes.

The Structural Basis for Activation and Inhibition of ZAP-70 Kinase Domain

ZAP–70 (Zeta-chain-associated protein kinase 70) is a tyrosine kinase that interacts directly with the activated T-cell receptor to transduce downstream signals, and is hence a major player in the regulation of the adaptive immune response. Dysfunction of ZAP–70 causes selective T cell deficiency that in turn results in persistent infections. ZAP–70 is activated by a variety of signals including phosphorylation of the kinase domain (KD), and binding of its regulatory tandem Src homology 2 (SH2) domains to the T cell receptor. The present study investigates molecular mechanisms of activation and inhibition of ZAP–70 via atomically detailed molecular dynamics simulation approaches. We report microsecond timescale simulations of five distinct states of the ZAP–70 KD, comprising apo, inhibited and three phosphorylated variants. Extensive analysis of local flexibility and correlated motions reveal crucial transitions between the states, thus elucidating crucial steps in the activation mechanism of the ZAP–70 KD. Furthermore, we rationalize previously observed staurosporine-bound crystal structures, suggesting that whilst the KD superficially resembles an “active-like” conformation, the inhibitor modulates the underlying protein dynamics and restricts it in a compact, rigid state inaccessible to ligands or cofactors. Finally, our analysis reveals a novel, potentially druggable pocket in close proximity to the activation loop of the kinase, and we subsequently use its structure in fragment-based virtual screening to develop a pharmacophore model. The pocket is distinct from classical type I or type II kinase pockets, and its discovery offers promise in future design of specific kinase inhibitors, whilst mutations in residues associated with this pocket are implicated in immunodeficiency in humans.

ZAP–70 is a key protein kinase in the adaptive immune system. It is essential for development and function of T cells and natural killer cells, and associated mutations can lead to conditions such as severe combined immunodeficiency (SCID). Here, simulations of the ZAP–70 kinase domain are used to study its dynamics in response to different mechanistic signals. We identify crucial motions over microsecond timescales, which help to rationalize in atomic detail previous structural and experimental data regarding its biological regulation. We subsequently propose a scheme for the phosphorylation-dependent activation cascade of ZAP–70, and for its ligand-dependent inhibition. Finally, we characterize a novel cryptic pocket adjacent to the active site and activation loop, which is distinct from classical type I or type II kinase sites. The pocket is in close proximity to several residues whose mutations cause SCID in humans, and its identification offers promise in future drug design efforts.

Monoamine transporters: insights from molecular dynamics simulations

The human monoamine transporters (MATs) facilitate the reuptake of the neurotransmitters serotonin, dopamine, and norepinephrine from the synaptic cleft. Imbalance in monoaminergic neurotransmission is linked to various diseases including major depression, attention deficit hyperactivity disorder, schizophrenia, and Parkinson's disease. Inhibition of the MATs is thus an important strategy for treatment of such diseases. The MATs are sodium-coupled transport proteins belonging to the neurotransmitter/Na(+) symporter (NSS) family, and the publication of the first high-resolution structure of a NSS family member, the bacterial leucine transporter LeuT, in 2005, proved to be a major stepping stone for understanding this family of transporters. Structural data allows for the use of computational methods to study the MATs, which in turn has led to a number of important discoveries. The process of substrate translocation across the membrane is an intrinsically dynamic process. Molecular dynamics simulations, which can provide atomistic details of molecular motion on ns to ms timescales, are therefore well-suited for studying transport processes. In this review, we outline how molecular dynamics simulations have provided insight into the large scale motions associated with transport of the neurotransmitters, as well as the presence of external and internal gates, the coupling between ion and substrate transport, and differences in the conformational changes induced by substrates and inhibitors.

Application of a BOSS-Gaussian interface for QM/MM simulations of Henry and methyl transfer reactions

Hybrid quantum mechanics and molecular mechanics (QM/MM) computer simulations have become an indispensable tool for studying chemical and biological phenomena for systems too large to treat with QM alone. For several decades, semiempirical QM methods have been used in QM/MM simulations. However, with increased computational resources, the introduction of ab initio and density function methods into on-the-fly QM/MM simulations is being increasingly preferred. This adaptation can be accomplished with a program interface that tethers independent QM and MM software packages. This report introduces such an interface for the BOSS and Gaussian programs, featuring modification of BOSS to request QM energies and partial atomic charges from Gaussian. A customizable C-shell linker script facilitates the interprogram communication. The BOSS-Gaussian interface also provides convenient access to Charge Model 5 (CM5) partial atomic charges for multiple purposes including QM/MM studies of reactions. In this report, the BOSS-Gaussian interface is applied to a nitroaldol (Henry) reaction and two methyl transfer reactions in aqueous solution. Improved agreement with experiment is found by determining free-energy surfaces with MP2/CM5 QM/MM simulations than previously reported investigations using semiempirical methods.

The true nature of rotary movements in rotaxanes

Reveal the intricate nature of movements within rotaxanes by means of multidimensional free-energy landscapes.

Disentangling the different movements observed in rotaxanes is critical to characterize their function as molecular and biological motors. How to achieve unidirectional rotation is an important question for successful construction of a highly efficient molecular motor. The motions within a rotaxane composed of a benzylic amide ring threaded on a fumaramide moiety were investigated employing atomistic molecular dynamics simulations. The free-energy profiles describing the rotational process of the ring about the thread were determined from multi-microsecond simulations. Comparing the theoretical free-energy barriers with their experimental counterpart, the syn–anti isomerization of the amide bond within the ring was ruled out. The free-energy barriers arise in fact from the disruption of hydrogen bonds between the ring and the thread. Transition path analysis reveals that complete description of the reaction coordinate requires another collective variable. The free-energy landscape spanned by the two variables characterizing the coupled rotational and shuttling processes of the ring in the rotaxane was mapped. The calculated free-energy barrier, amounting to 9.3 kcal mol^–1, agrees well with experiment. Further analysis shows that shuttling is coupled with the isomerization of the ring, which is not limited to a simplistic chair-to-chair transition. This work provides a cogent example that contrary to chemical intuition, molecular motion can result from complex, entangled movements requiring for their accurate description careful modeling of the underlying reaction coordinate. The methodology described here can be used to evaluate the different components of the multifaceted motion in rotaxanes, and constitutes a robust tool for the rational design of molecular machines.

Effects of different self-assembled monolayers on thin-film morphology: a combined DFT/MD simulation protocol

Organic thin film transistors (OTFTs) are multilayer field-effect transistors that employ an organic conjugated material as semiconductor. Several experimental groups have recently demonstrated that the insertion of an organic self-assembled monolayer (SAM) between the dielectric and the semiconductive layer is responsible for a sensible improvement of the OTFT performances in terms of an increased charge carrier mobility caused by a higher degree of order in the organic semiconductor layer. Here, we describe a combined periodic density functional theory (DFT) and classical molecular dynamics (MD) protocol applied to four different SAMs and a pentacene monolayer deposited onto their surfaces. In particular, we investigate the morphology and the surface of the four SAMs and the translational, orientational, and nematic order of the monolayer through the calculation of several distribution functions and order parameters pointing out the differences among the systems and relating them to known experimental results. Our calculations also suggest that small differences in the SAM molecular design will produce remarkable differences in the SAM surface and monolayer order. In particular, our simulations explain how a SAM with a bulky terminal group results in an irregular and rough surface that determines the deposition of a disordered semiconductive monolayer. On the contrary, SAMs with a small terminal group generate smooth surfaces with uninterrupted periodicity, thus favoring the formation of an ordered pentacene monolayer that increases the mobility of charge carriers and improves the overall performances of the OTFT devices. Our results clearly point out that the in silico procedure presented here might be of help in tuning the design of SAMs in order to improve the quality of OTFT devices.

Interactions of L-3,5,3'-Triiodothyronine [corrected], Allopregnanolone, and Ivermectin with the GABAA Receptor: Evidence for Overlapping Intersubunit Binding Modes

Structural mechanisms of modulation of γ-aminobutyric acid (GABA) type A receptors by neurosteroids and hormones remain unclear. The thyroid hormone L-3,5,3’-triiodothyronine (T3) inhibits GABA_A receptors at micromolar concentrations and has common features with neurosteroids such as allopregnanolone (ALLOP). Here we use functional experiments on α₂β₁γ₂ GABA_A receptors expressed in Xenopus oocytes to detect competitive interactions between T3 and an agonist (ivermectin, IVM) with a crystallographically determined binding site at subunit interfaces in the transmembrane domain of a homologous receptor (glutamate-gated chloride channel, GluCl). T3 and ALLOP also show competitive effects, supporting the presence of both a T3 and ALLOP binding site at one or more subunit interfaces. Molecular dynamics (MD) simulations over 200 ns are used to investigate the dynamics and energetics of T3 in the identified intersubunit sites. In these simulations, T3 molecules occupying all intersubunit sites (with the exception of the α-β interface) display numerous energetically favorable conformations with multiple hydrogen bonding partners, including previously implicated polar/acidic sidechains and a structurally conserved deformation in the M1 backbone.

Understanding the roles of Lys33 and Arg45 in the binding-site stability of LjLTP10, an LTP related to drought stress in Lotus japonicus

In Lotus japonicus, as in most plants, long-chain fatty acids are important components of cuticular wax, one of the principal functions of which is to act as a barrier to water loss in response to drought stress. It is thought that lipid transfer proteins (LTPs) are involved in the process of cuticle formation. We previously described LjLTP10 as an LTP involved in cuticle formation during acclimation response to drought stress in L. japonicus. The structural model of LjLTP10 had two residues (K33 and R45) in the hydrophobic cavity, although the role of these residues was unclear. In the present work, we investigated the molecular mechanism involved in the transport of lipid precursors in L. japonicus and clarified the importance of the residues K33 and R45. First, in silico site-directed mutagenesis studies were carried out on the LjLTP10 structure. Structural analysis showed that LjLTP10 mutants possess similar structures but their hydrophobic cavities are somewhat different. Unfavorable energies for the interactions of the mutant proteins with different ligands were found by molecular docking and molecular dynamics simulations. We also examined the contributions of energetic parameters to the free energy of the protein-ligand complex using the MM-GBSA method. Results showed that the different complexes present similar, favorable van der Waals interactions, whereas electrostatic interactions were not favored in the mutant structures. Our study indicates that the residues K33 and R45 play a crucial role in maintaining the binding pocket structure required for lipid transport.

Diameter-Selective Dispersion of Carbon Nanotubes via Polymers: A Competition between Adsorption and Bundling

The mechanism of the selective dispersion of single-walled carbon nanotubes (CNTs) by polyfluorene polymers is studied in this paper. Using extensive molecular dynamics simulations, it is demonstrated that diameter selectivity is the result of a competition between bundling of CNTs and adsorption of polymers on CNT surfaces. The preference for certain diameters corresponds to local minima of the binding energy difference between these two processes. Such minima in the diameter dependence occur due to abrupt changes in the CNT's coverage with polymers, and their calculated positions are in quantitative agreement with preferred diameters reported experimentally. The presented approach defines a theoretical framework for the further understanding and improvement of dispersion/extraction processes.

Binding enthalpy calculations for a neutral host-guest pair yield widely divergent salt effects across water models

Dissolved salts are a part of the physiological milieu and can significantly influence the kinetics and thermodynamics of various biomolecular processes, such as binding and catalysis; thus, it is important for molecular simulations to reliably describe their effects. The present study uses a simple, nonionized host-guest model system to study the sensitivity of computed binding enthalpies to the choice of water and salt models. Molecular dynamics simulations of a cucurbit[7]uril host with a neutral guest molecule show striking differences in the salt dependency of the binding enthalpy across four water models, TIP3P, SPC/E, TIP4P-Ew, and OPC, with additional sensitivity to the choice of parameters for sodium and chloride. In particular, although all of the models predict that binding will be less exothermic with increasing NaCl concentration, the strength of this effect varies by 7 kcal/mol across models. The differences appear to result primarily from differences in the number of sodium ions displaced from the host upon binding the guest rather than from differences in the enthalpy associated with this displacement, and it is the electrostatic energy that contributes most to the changes in enthalpy with increasing salt concentration. That a high sensitivity of salt affecting the choice of water model, as observed for the present host-guest system despite it being nonionized, raises issues regarding the selection and adjustment of water models for use with biological macromolecules, especially as these typically possess multiple ionized groups that can interact relatively strongly with ions in solution.

Biomimetics: From Bioinformatics to Rational Design of Dendrimers as Gene Carriers

Biomimetics, or the use of principles of Nature for developing new materials, is a paradigm that could help Nanomedicine tremendously. One of the current challenges in Nanomedicine is the rational design of new efficient and safer gene carriers. Poly(amidoamine) (PAMAM) dendrimers are a well-known class of nanoparticles, extensively used as non-viral nucleic acid carriers, due to their positively charged end-groups. Yet, there are still several aspects that can be improved for their successful application in in vitro and in vivo systems, including their affinity for nucleic acids as well as lowering their cytotoxicity. In the search of new functional groups that could be used as new dendrimer-reactive groups, we followed a biomimetic approach to determine the amino acids with highest prevalence in protein-DNA interactions. Then we introduced them individually as terminal groups of dendrimers, generating a new class of nanoparticles. Molecular dynamics studies of two systems: PAMAM-Arg and PAMAM-Lys were also performed in order to describe the formation of complexes with DNA. Results confirmed that the introduction of amino acids as terminal groups in a dendrimer increases their affinity for DNA and the interactions in the complexes were characterized at atomic level. We end up by briefly discussing additional modifications that can be made to PAMAM dendrimers to turned them into promising new gene carriers.

Free energy and hidden barriers of the β-sheet structure of prion protein

On-the-fly free-energy parametrization is a new collective variable biasing approach akin to metadynamics with one important distinction: rather than acquiring an accelerated distribution via a history-dependent bias potential, sampling on this distribution is achieved from the beginning of the simulation using temperature-accelerated molecular dynamics. In the present work, we compare the performance of both approaches to compute the free-energy profile along a scalar collective variable measuring the H-bond registry of the β-sheet structure of the mouse Prion protein. Both methods agree on the location of the free-energy minimum, but free-energy profiles from well-tempered metadynamics are subject to a much higher degree of statistical noise due to hidden barriers. The sensitivity of metadynamics to hidden barriers is shown to be a consequence of the history dependence of the bias potential, and we detail the nature of these barriers for the prion β-sheet. In contrast, on-the-fly parametrization is much less sensitive to these barriers and thus displays improved convergence behavior relative to that of metadynamics. While hidden barriers are a frequent and central issue in free-energy methods, on-the-fly free-energy parametrization appears to be a robust and preferable method to confront this issue.

Realistic sampling of amino acid geometries for a multipolar polarizable force field

The Quantum Chemical Topological Force Field (QCTFF) uses the machine learning method kriging to map atomic multipole moments to the coordinates of all atoms in the molecular system. It is important that kriging operates on relevant and realistic training sets of molecular geometries. Therefore, we sampled single amino acid geometries directly from protein crystal structures stored in the Protein Databank (PDB). This sampling enhances the conformational realism (in terms of dihedral angles) of the training geometries. However, these geometries can be fraught with inaccurate bond lengths and valence angles due to artefacts of the refinement process of the X-ray diffraction patterns, combined with experimentally invisible hydrogen atoms. This is why we developed a hybrid PDB/nonstationary normal modes (NM) sampling approach called PDB/NM. This method is superior over standard NM sampling, which captures only geometries optimized from the stationary points of single amino acids in the gas phase. Indeed, PDB/NM combines the sampling of relevant dihedral angles with chemically correct local geometries. Geometries sampled using PDB/NM were used to build kriging models for alanine and lysine, and their prediction accuracy was compared to models built from geometries sampled from three other sampling approaches. Bond length variation, as opposed to variation in dihedral angles, puts pressure on prediction accuracy, potentially lowering it. Hence, the larger coverage of dihedral angles of the PDB/NM method does not deteriorate the predictive accuracy of kriging models, compared to the NM sampling around local energetic minima used so far in the development of QCTFF.

Modeling Vibrational Spectra of Ester Carbonyl Stretch in Water and DMSO Based on Molecular Dynamics Simulation

On the basis of molecular dynamics simulation, we model the ester carbonyl stretch FTIR signals of methyl acetate in D2O and DMSO. An ab initio map is constructed at the B3LYP/6-311++G** level to relate the carbonyl stretch frequency to the external electric field. Using this map, fluctuating Hamiltonian of the carbonyl stretch is constructed from the MD simulation trajectory. The IR spectra calculated based on this Hamiltonian are found to be in good agreement with the experiment. For methyl acetate in D2O, hydrogen bonding on alkoxy oxygen causes a blue shift of frequency, while that on carbonyl oxygen causes a red shift. Two peaks observed in FTIR signals originate from the balance of these two effects. Furthermore, in both D2O and DMSO solutions, correlations are found between the instantaneous electric field on C═O and the frequencies. Broader line width of the signal in D2O suggests a more inhomogeneous electric field distribution due to the complicated hydrogen-bonding environment.

Stepwise substrate translocation mechanism revealed by free energy calculations of doxorubicin in the multidrug transporter AcrB

AcrB is the inner membrane transporter of the tripartite multidrug efflux pump AcrAB-TolC in E. coli, which poses a major obstacle to the treatment of bacterial infections. X-ray structures have identified two types of substrate-binding pockets in the porter domains of AcrB trimer: the proximal binding pocket (PBP) and the distal binding pocket (DBP), and suggest a functional rotating mechanism in which each protomer cycles consecutively through three distinct conformational states (access, binding and extrusion). However, the details of substrate binding and translocation between the binding pockets remain elusive. In this work, we performed atomic simulations to obtain the free energy profile of the translocation of an antibiotic drug doxorubicin (DOX) inside AcrB. Our simulation indicates that DOX binds at the PBP and DBP with comparable affinities in the binding state protomer, and overcomes a 3 kcal/mol energy barrier to transit between them. Obvious conformational changes including closing of the PC1/PC2 cleft and shrinking of the DBP were observed upon DOX binding in the PBP, resulting in an intermediate state between the access and binding states. Taken together, the simulation results reveal a detailed stepwise substrate binding and translocation process in the framework of functional rotating mechanism.

Solution NMR characterization of chemokine CXCL8/IL-8 monomer and dimer binding to glycosaminoglycans: structural plasticity mediates differential binding interactions

Structural plasticity plays a major role in determining differential binding of CXCL8 monomer and dimer to glycosaminoglycans (GAGs) and that dimer is the high-affinity GAG ligand. We propose that these properties play important roles in orchestrating in vivo chemokine-mediated neutrophil function.

Chemokine CXCL8/interleukin-8 (IL-8) plays a crucial role in directing neutrophils and oligodendrocytes to combat infection/injury and tumour cells in metastasis development. CXCL8 exists as monomers and dimers and interaction of both forms with glycosaminoglycans (GAGs) mediate these diverse cellular processes. However, very little is known regarding the structural basis underlying CXCL8–GAG interactions. There are conflicting reports on the affinities, geometry and whether the monomer or dimer is the high-affinity GAG ligand. To resolve these issues, we characterized the binding of a series of heparin-derived oligosaccharides [heparin disaccharide (dp2), heparin tetrasaccharide (dp4), heparin octasaccharide (dp8) and heparin 14-mer (dp14)] to the wild-type (WT) dimer and a designed monomer using solution NMR spectroscopy. The pattern and extent of binding-induced chemical shift perturbation (CSP) varied between dimer and monomer and between longer and shorter oligosaccharides. NMR-based structural models show that different interaction modes coexist and that the nature of interactions varied between monomer and dimer and oligosaccharide length. MD simulations indicate that the binding interface is structurally plastic and provided residue-specific details of the dynamic nature of the binding interface. Binding studies carried out under conditions at which WT CXCL8 exists as monomers and dimers provide unambiguous evidence that the dimer is the high-affinity GAG ligand. Together, our data indicate that a set of core residues function as the major recognition/binding site, a set of peripheral residues define the various binding geometries and that the structural plasticity of the binding interface allows multiplicity of binding interactions. We conclude that structural plasticity most probably regulates in vivo CXCL8 monomer/dimer–GAG interactions and function.

The Internal Dynamics of Fibrinogen and Its Implications for Coagulation and Adsorption

Fibrinogen is a serum multi-chain protein which, when activated, aggregates to form fibrin, one of the main components of a blood clot. Fibrinolysis controls blood clot dissolution through the action of the enzyme plasmin, which cleaves fibrin at specific locations. Although the main biochemical factors involved in fibrin formation and lysis have been identified, a clear mechanistic picture of how these processes take place is not available yet. This picture would be instrumental, for example, for the design of improved thrombolytic or anti-haemorrhagic strategies, as well as, materials with improved biocompatibility. Here, we present extensive molecular dynamics simulations of fibrinogen which reveal large bending motions centered at a hinge point in the coiled-coil regions of the molecule. This feature, likely conserved across vertebrates according to our analysis, suggests an explanation for the mechanism of exposure to lysis of the plasmin cleavage sites on fibrinogen coiled-coil region. It also explains the conformational variability of fibrinogen observed during its adsorption on inorganic surfaces and it is supposed to play a major role in the determination of the hydrodynamic properties of fibrinogen. In addition the simulations suggest how the dynamics of the D region of fibrinogen may contribute to the allosteric regulation of the blood coagulation cascade through a dynamic coupling between the a- and b-holes, important for fibrin polymerization, and the integrin binding site P1.

Fibrinogen, a protein found in the blood of vertebrates, when activated, aggregates and forms fibrin fibers, the basis of a blood clot. Clots are broken down by the enzyme plasmin, which cuts fibrin fibers at specific places, thus helping the regulation of clot persistence. A mechanistic understanding of fibrin degradation by plasmin is still missing. An important determinant of this process might be the flexibility of fibrinogen. The flexible nature of fibrinogen is reported, for example, by the great variety of conformations observed when fibrinogen adsorbs on material surfaces. However, limits in the spatial resolution of these experiments preclude the identification of the atomistic mechanism behind this flexibility. Here, we perform computer simulations that help identifying with atomistic detail large bending motions occurring at a specific hinge on the molecule. We show how these bending motions can explain the variable conformations observed in experiments and how they help exposing sites where plasmin can cut fibrinogen. Furthermore, our simulations let us identify cooperative effects involving several distant parts of fibrinogen that may play a role in the assembly of fibrin fibers. Both the bending and the cooperative effects, thus, represent potential mechanisms for the regulation of blood clotting.

Mechanism of the Glycosidic Bond Cleavage of Mismatched Thymine in Human Thymine DNA Glycosylase Revealed by Classical Molecular Dynamics and Quantum Mechanical/Molecular Mechanical Calculations

Base excision of mismatched or damaged nucleotides catalyzed by glycosylase enzymes is the first step of the base excision repair system, a machinery preserving the integrity of DNA. Thymine DNA glycosylase recognizes and removes mismatched thymine by cleaving the C1'-N1 bond between the base and the sugar ring. Our quantum mechanical/molecular mechanical calculations of this reaction in human thymine DNA glycosylase reveal a requirement for a positive charge in the active site to facilitate C1'-N1 bond scission: protonation of His151 significantly lowers the free energy barrier for C1'-N1 bond dissociation compared to the situation with neutral His151. Shuttling a proton from His151 to the thymine base further reduces the activation free energy for glycosidic bond cleavage. Classical molecular dynamics simulations of the H151A mutant suggest that the mutation to the smaller, neutral, residue increases the water accessibility of the thymine base, rendering direct proton transfer from the bulk feasible. Quantum mechanical/molecular mechanical calculations of the glycosidic bond cleavage reaction in the H151A mutant show that the activation free energy is slightly lower than in the wild-type enzyme, explaining the experimentally observed higher reaction rates in this mutant.

Control of protein orientation on gold nanoparticles

Gold nanoparticles (Au NPs) have attracted much attention due to their potential applications in nano-medicine. While numerous studies have quantified biomolecular adsorption to Au NPs in terms of equilibrium binding constants, far less is known about biomolecular orientation on nanoparticle surfaces. In this study, the binding of the protein α-synuclein to citrate and (16-mercaptohexadecyl) trimethylammonium bromide (MTAB) coated 12 nm Au NPs is examined by heteronuclear single quantum coherence NMR spectroscopy to provide site-specific measurements of protein-nanoparticle binding. Molecular dynamics simulations support the orientation assignments, which show N-terminus binding to the Au NP for citrate-capped NPs, and C-terminus binding for the MTAB-capped NPs.

A simple and transferable all-atom/coarse-grained hybrid model to study membrane processes

We present an efficient all-atom/coarse-grained hybrid model and apply it to membrane processes. This model is an extension of the all-atom/ELBA model applied previously to processes in water. Here, we improve the efficiency of the model by implementing a multiple-time step integrator that allows the atoms and the coarse-grained beads to be propagated at different timesteps. Furthermore, we fine-tune the interaction between the atoms and the coarse-grained beads by computing the potential of mean force of amino acid side chain analogs along the membrane normal and comparing to atomistic simulations. The model was independently validated on the calculation of small-molecule partition coefficients. Finally, we apply the model to membrane peptides. We studied the tilt angle of the Walp23 and Kalp23 helices in two different model membranes and the stability of the glycophorin A dimer. The model is efficient, accurate, and straightforward to use, as it does not require any extra interaction particles, layers of atomistic solvent molecules or tabulated potentials, thus offering a novel, simple approach to study membrane processes.

Inhibition of Peroxidase Activity of Cytochrome c: De Novo Compound Discovery and Validation

Cytochrome c (cyt c) release from mitochondria is accepted to be the point of no return for eliciting a cascade of interactions that lead to apoptosis. A strategy for containing sustained apoptosis is to reduce the mitochondrial permeability pore opening. Pore opening is enhanced by peroxidase activity of cyt c gained upon its complexation with cardiolipin in the presence of reactive oxygen species. Blocking access to the heme group has been proposed as an effective intervention method for reducing, if not eliminating, the peroxidase activity of cyt c. In the present study, using a combination of druggability simulations, pharmacophore modeling, virtual screening, and in vitro fluorescence measurements to probe peroxidase activity, we identified three repurposable drugs and seven compounds that are validated to effectively inhibit the peroxidase activity of cyt c.

Molecular Simulations of Solved Co-crystallized X-Ray Structures Identify Action Mechanisms of PDEδ Inhibitors

PDEδ is a small protein that binds and controls the trafficking of RAS subfamily proteins. Its inhibition protects initiation of RAS signaling, and it is one of the common targets considered for oncological drug development. In this study, we used solved x-ray structures of inhibitor-bound PDEδ targets to investigate mechanisms of action of six independent all-atom MD simulations. An analysis of atomic simulations combined with the molecular mechanic-Poisson-Boltzmann solvent accessible surface area/generalized Born solvent accessible surface area calculations led to the identification of action mechanisms for a panel of novel PDEδ inhibitors. To the best of our knowledge, this study is one of the first in silico investigations on co-crystallized PDEδ protein. A detailed atomic-scale understanding of the molecular mechanism of PDEδ inhibition may assist in the design of novel PDEδ inhibitors. One of the most common side effects for diverse small molecules/kinase inhibitors is their off-target interactions with cardiac ion channels and human-ether-a-go-go channel specifically. Thus, all of the studied PDEδ inhibitors are also screened in silico at the central cavities of hERG1 potassium channels.

Boramino acid as a marker for amino acid transporters

Amino acid transporters (AATs) are a series of integral channels for uphill cellular uptake of nutrients and neurotransmitters. Abnormal expression of AATs is often associated with cancer, addiction, and multiple mental diseases. Although methods to evaluate in vivo expression of AATs would be highly useful, efforts to develop them have been hampered by a lack of appropriate tracers. We describe a new class of AA mimics-boramino acids (BAAs)-that can serve as general imaging probes for AATs. The structure of a BAA is identical to that of the corresponding natural AA, except for an exotic replacement of the carboxylate with -BF3 (-). Cellular studies demonstrate strong AAT-mediated cell uptake, and animal studies show high tumor-specific accumulation, suggesting that BAAs hold great promise for the development of new imaging probes and smart AAT-targeting drugs.