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

CSM-lig: a web server for assessing and comparing protein-small molecule affinities

Determining the affinity of a ligand for a given protein is a crucial component of drug development and understanding their biological effects. Predicting binding affinities is a challenging and difficult task, and despite being regarded as poorly predictive, scoring functions play an important role in the analysis of molecular docking results. Here, we present CSM-Lig (http://structure.bioc.cam.ac.uk/csm_lig), a web server tailored to predict the binding affinity of a protein-small molecule complex, encompassing both protein and small-molecule complementarity in terms of shape and chemistry via graph-based structural signatures. CSM-Lig was trained and evaluated on different releases of the PDBbind databases, achieving a correlation of up to 0.86 on 10-fold cross validation and 0.80 in blind tests, performing as well as or better than other widely used methods. The web server allows users to rapidly and automatically predict binding affinities of collections of structures and assess the interactions made. We believe CSM-lig would be an invaluable tool for helping assess docking poses, the effects of multiple mutations, including insertions, deletions and alternative splicing events, in protein-small molecule affinity, unraveling important aspects that drive protein–compound recognition.

Solute-Solvent Energetics Based on Proximal Distribution Functions

We consider the hydration structure and thermodynamic energetics of solutes in aqueous solution. On the basis of the dominant local correlation between the solvent and the chemical nature of the solute atoms, proximal distribution functions (pDF) can be used to quantitatively estimate the hydration pattern of the macromolecules. We extended this technique to study the solute-solvent energetics including the van der Waals terms representing excluded volume and tested the method with butane and propanol. Our results indicate that the pDF-reconstruction algorithm can reproduce van der Waals solute-solvent interaction energies to useful kilocalorie per mole accuracy. We subsequently computed polyalanine-water interaction energies for a variety of conformers, which also showed agreement with the simulated values.

Development of Potent and Selective Inhibitors for Group VIA Calcium-Independent Phospholipase A2 Guided by Molecular Dynamics and Structure-Activity Relationships

The development of inhibitors for phospholipase A2 (PLA2) is important in elucidating the enzymes implication in various biological pathways. PLA2 enzymes are an important pharmacological target implicated in various inflammatory diseases. Computational chemistry, organic synthesis, and in vitro assays were employed to develop potent and selective inhibitors for group VIA calcium-independent PLA2. A set of fluoroketone inhibitors was studied for their binding mode with two human cytosolic PLA2 enzymes: group IVA cPLA2 and group VIA iPLA2. New compounds were synthesized and assayed toward three major PLA2s. This study led to the development of four potent and selective thioether fluoroketone inhibitors as well as a thioether keto-1,2,4-oxadiazole inhibitor for GVIA iPLA2, which will serve as lead compounds for future development and studies. The keto-1,2,4-oxadiazole functionality with a thioether is a novel structure, and it will be used as a lead to develop inhibitors with higher potency and selectivity toward GVIA iPLA2.

Atomistic Simulation of Solubilization of Polycyclic Aromatic Hydrocarbons in a Sodium Dodecyl Sulfate Micelle

Solubilization of two polycyclic aromatic hydrocarbons (PAHs), naphthalene (NAP, 2-benzene-ring PAH) and pyrene (PYR, 4-benzene-ring PAH), into a sodium dodecyl sulfate (SDS) micelle was studied through all-atom molecular dynamics (MD) simulations. We find that NAP as well as PYR could move between the micelle shell and core regions, contributing to their distribution in both regions of the micelle at any PAH concentration. Moreover, both NAP and PYR prefer to stay in the micelle shell region, which may arise from the greater volume of the micelle shell, the formation of hydrogen bonds between NAP and water, and the larger molecular volume of PYR. The PAHs are able to form occasional clusters (from dimer to octamer) inside the micelle during the simulation time depending on the PAH concentration in the solubilization systems. Furthermore, the micelle properties (i.e., size, shape, micelle internal structure, alkyl chain conformation and orientation, and micelle internal dynamics) are found to be nearly unaffected by the solubilized PAHs, which is irrespective of the properties and concentrations of PAHs.

Simulation-Based Approaches for Determining Membrane Permeability of Small Compounds

Predicting the rate of nonfacilitated permeation of solutes across lipid bilayers is important to drug design, toxicology, and signaling. These rates can be estimated using molecular dynamics simulations combined with the inhomogeneous solubility-diffusion model, which requires calculation of the potential of mean force and position-dependent diffusivity of the solute along the transmembrane axis. In this paper, we assess the efficiency and accuracy of several methods for the calculation of the permeability of a model DMPC bilayer to urea, benzoic acid, and codeine. We compare umbrella sampling, replica exchange umbrella sampling, adaptive biasing force, and multiple-walker adaptive biasing force for the calculation of the transmembrane PMF. No definitive advantage for any of these methods in their ability to predict the membrane permeability coefficient Pm was found, provided that a sufficiently long equilibration is performed. For diffusivities, a Bayesian inference method was compared to a generalized Langevin method, both being sensitive to chosen parameters and the slow relaxation of membrane defects. Agreement within 1.5 log units of the computed Pm with experiment is found for all permeants and methods. Remaining discrepancies can likely be attributed to limitations of the force field as well as slowly relaxing collective movements within the lipid environment. Numerical calculations based on model profiles show that Pm can be reliably estimated from only a few data points, leading to recommendations for calculating Pm from simulations.

Augmenting drug-carrier compatibility improves tumour nanotherapy efficacy

A major goal of cancer nanotherapy is to use nanoparticles as carriers for targeted delivery of anti-tumour agents. The drug-carrier association after intravenous administration is essential for efficient drug delivery to the tumour. However, a large number of currently available nanocarriers are self-assembled nanoparticles whose drug-loading stability is critically affected by the in vivo environment. Here we used in vivo FRET imaging to systematically investigate how drug-carrier compatibility affects drug release in a tumour mouse model. We found the drug's hydrophobicity and miscibility with the nanoparticles are two independent key parameters that determine its accumulation in the tumour. Next, we applied these findings to improve chemotherapeutic delivery by augmenting the parent drug's compatibility; as a result, we achieved better antitumour efficacy. Our results help elucidate nanomedicines' in vivo fate and provide guidelines for efficient drug delivery.

Atomistic Molecular Dynamics Simulations of Crude Oil/Brine Displacement in Calcite Mesopores

Unconventional reservoirs such as hydrocarbon-bearing shale formations and ultratight carbonates generate a large fraction of oil and gas production in North America. The characteristic feature of these reservoirs is their nanoscale porosity that provides significant surface areas between the pore walls and the occupying fluids. To better assess hydrocarbon recovery from these formations, it is crucial to develop an improved insight into the effects of wall-fluid interactions on the interfacial phenomena in these nanoscale confinements. One of the important properties that controls the displacement of fluids inside the pores is the threshold capillary pressure. In this study, we present the results of an integrated series of large-scale molecular dynamics (MD) simulations performed to investigate the effects of wall-fluid interactions on the threshold capillary pressures of oil-water/brine displacements in a calcite nanopore with a square cross section. Fully atomistic models are utilized to represent crude oil, brine, and calcite in order to accommodate electrostatic interactions and H-bonding between the polar molecules and the calcite surface. To this end, we create mixtures of various polar and nonpolar organic molecules to better represent the crude oil. The interfacial tension between oil and water/brine and their contact angle on calcite surface are simulated. We study the effects of oil composition, water salinity, and temperature and pressure conditions on these properties. The threshold capillary pressure values are also obtained from the MD simulations for the calcite nanopore. We then compare the MD results against those generated using the Mayer-Stowe-Princen (MSP) method and explain the differences.

Kainate receptor pore-forming and auxiliary subunits regulate channel block by a novel mechanism

KEY POINTS

Kainate receptor heteromerization and auxiliary subunits, Neto1 and Neto2, attenuate polyamine ion-channel block by facilitating blocker permeation. Relief of polyamine block in GluK2/GluK5 heteromers results from a key proline residue that produces architectural changes in the channel pore α-helical region. Auxiliary subunits exert an additive effect to heteromerization, and thus relief of polyamine block is due to a different mechanism. Our findings have broad implications for work on polyamine block of other cation-selective ion channels.

ABSTRACT

Channel block and permeation by cytoplasmic polyamines is a common feature of many cation-selective ion channels. Although the channel block mechanism has been studied extensively, polyamine permeation has been considered less significant as it occurs at extreme positive membrane potentials. Here, we show that kainate receptor (KAR) heteromerization and association with auxiliary proteins, Neto1 and Neto2, attenuate polyamine block by enhancing blocker permeation. Consequently, polyamine permeation and unblock occur at more negative and physiologically relevant membrane potentials. In GluK2/GluK5 heteromers, enhanced permeation is due to a single proline residue in GluK5 that alters the dynamics of the α-helical region of the selectivity filter. The effect of auxiliary proteins is additive, and therefore the structural basis of polyamine permeation and unblock is through a different mechanism. As native receptors are thought to assemble as heteromers in complex with auxiliary proteins, our data identify an unappreciated impact of polyamine permeation in shaping the signalling properties of neuronal KARs and point to a structural mechanism that may be shared amongst other cation-selective ion channels.

Unraveling the dynamics and structure of functionalized self-assembled monolayers on gold using 2D IR spectroscopy and MD simulations

Functionalized self-assembled monolayers (SAMs) are the focus of ongoing investigations because they can be chemically tuned to control their structure and dynamics for a wide variety of applications, including electrochemistry, catalysis, and as models of biological interfaces. Here we combine reflection 2D infrared vibrational echo spectroscopy (R-2D IR) and molecular dynamics simulations to determine the relationship between the structures of functionalized alkanethiol SAMs on gold surfaces and their underlying molecular motions on timescales of tens to hundreds of picoseconds. We find that at higher head group density, the monolayers have more disorder in the alkyl chain packing and faster dynamics. The dynamics of alkanethiol SAMs on gold are much slower than the dynamics of alkylsiloxane SAMs on silica. Using the simulations, we assess how the different molecular motions of the alkyl chain monolayers give rise to the dynamics observed in the experiments.

Novel Analogues of (R)-5-(Methylamino)-5,6-dihydro-4H-imidazo[4,5,1-ij]quinolin-2(1H)-one (Sumanirole) Provide Clues to Dopamine D2/D3 Receptor Agonist Selectivity

Novel 1-, 5-, and 8-substituted analogues of sumanirole (1), a dopamine D2/D3 receptor (D2R/D3R) agonist, were synthesized. Binding affinities at both D2R and D3R were higher when determined in competition with the agonist radioligand [(3)H]7-hydroxy-N,N-dipropyl-2-aminotetralin (7-OH-DPAT) than with the antagonist radioligand [(3)H]N-methylspiperone. Although 1 was confirmed as a D2R-preferential agonist, its selectivity in binding and functional studies was lower than previously reported. All analogues were determined to be D2R/D3R agonists in both GoBRET and mitogenesis functional assays. Loss of efficacy was detected for the N-1-substituted analogues at D3R. In contrast, the N-5-alkyl-substituted analogues, and notably the n-butyl-arylamides (22b and 22c), all showed improved affinity at D2R over 1 with neither a loss of efficacy nor an increase in selectivity. Computational modeling provided a structural basis for the D2R selectivity of 1, illustrating how subtle differences in the highly homologous orthosteric binding site (OBS) differentially affect D2R/D3R affinity and functional efficacy.

Flexible CDOCKER: Development and application of a pseudo-explicit structure-based docking method within CHARMM

Protein-ligand docking is a commonly used method for lead identification and refinement. While traditional structure-based docking methods represent the receptor as a rigid body, recent developments have been moving toward the inclusion of protein flexibility. Proteins exist in an interconverting ensemble of conformational states, but effectively and efficiently searching the conformational space available to both the receptor and ligand remains a well-appreciated computational challenge. To this end, we have developed the Flexible CDOCKER method as an extension of the family of complete docking solutions available within CHARMM. This method integrates atomically detailed side chain flexibility with grid-based docking methods, maintaining efficiency while allowing the protein and ligand configurations to explore their conformational space simultaneously. This is in contrast to existing approaches that use induced-fit like sampling, such as Glide or Autodock, where the protein or the ligand space is sampled independently in an iterative fashion. Presented here are developments to the CHARMM docking methodology to incorporate receptor flexibility and improvements to the sampling protocol as demonstrated with re-docking trials on a subset of the CCDC/Astex set. These developments within CDOCKER achieve docking accuracy competitive with or exceeding the performance of other widely utilized docking programs.

How Far Does a Receptor Influence Vibrational Properties of an Odorant?

The biophysical mechanism of the sense of smell, or olfaction, is still highly debated. The mainstream explanation argues for a shape-based recognition of odorant molecules by olfactory receptors, while recent investigations suggest the primary olfactory event to be triggered by a vibrationally-assisted electron transfer reaction. We consider this controversy by studying the influence of a receptor on the vibrational properties of an odorant in atomistic details as the coupling between electronic degrees of freedom of the receptor and the vibrations of the odorant is the key parameter of the vibrationally-assisted electron transfer. Through molecular dynamics simulations we elucidate the binding specificity of a receptor towards acetophenone odorant. The vibrational properties of acetophenone inside the receptor are then studied by the polarizable embedding density functional theory approach, allowing to quantify protein-odorant interactions. Finally, we judge whether the effects of the protein provide any indications towards the existing theories of olfaction.

Computational Investigation of the Effect of Lipid Membranes on Ion Permeation in Gramicidin A

Membrane proteins are embedded in a lipid bilayer and interact with the lipid molecules in subtle ways. This can be studied experimentally by examining the effect of different lipid bilayers on the function of membrane proteins. Understanding the causes of the functional effects of lipids is difficult to dissect experimentally but more amenable to a computational approach. Here we perform molecular dynamics simulations and free energy calculations to study the effect of two lipid types (POPC and NODS) on the conductance of the gramicidin A (gA) channel. A larger energy barrier is found for the K⁺ potential of mean force in gA embedded in POPC compared to that in NODS, which is consistent with the enhanced experimental conductance of cations in gA embedded in NODS. Further analysis of the contributions to the potential energy of K⁺ reveals that gA and water molecules in gA make similar contributions in both bilayers but there are significant differences between the two bilayers when the lipid molecules and interfacial waters are considered. It is shown that the stronger dipole moments of the POPC head groups create a thicker layer of interfacial waters with better orientation, which ultimately is responsible for the larger energy barrier in the K⁺ PMF in POPC.

High-Throughput Fragment Docking into the BAZ2B Bromodomain: Efficient in Silico Screening for X-Ray Crystallography

Bromodomains are protein modules that bind to acetylated lysine side chains in histones and other proteins. The bromodomain adjacent to zinc finger domain protein 2B (BAZ2B) has been reported to be poorly druggable. Here, we screened an in-house library of 350 fragments by automatic docking to the BAZ2B bromodomain. The top 12 fragments according to the predicted binding energy were selected for experiments of soaking into apo crystals of BAZ2B which yielded the structure of the complex for four of them, which is a hit rate of 33%. Additional crystal structures were solved for BAZ2B and two scaffolds identified by analogy. For three topologically similar fragments, the crystal structures reveal binding modes with different penetration, i.e., with zero, one, and two water molecules, respectively, located between the fragment and the side chain of a conserved tyrosine (Tyr1901) in the bottom of the acetyl lysine pocket of BAZ2B. Furthermore, a remarkable stereoselectivity of the acetyl lysine pocket emerges from the crystal structures of the bromodomains of BAZ2B and SMARCA4 in complex with the chiral diol MPD (2-methyl-2,4-pentanediol).

Computational Study of Binding of μ-Conotoxin GIIIA to Bacterial Sodium Channels NaVAb and NaVRh

Structures of several voltage-gated sodium (NaV) channels from bacteria have been determined recently, but the same feat might not be achieved for the mammalian counterparts in the near future. Thus, at present, computational studies of the mammalian NaV channels have to be performed using homology models based on the bacterial crystal structures. A successful homology model for the mammalian NaV1.4 channel was recently constructed using the extensive mutation data for binding of μ-conotoxin GIIIA to NaV1.4, which was further validated through studies of binding of other μ-conotoxins and ion permeation. Understanding the similarities and differences between the bacterial and mammalian NaV channels is an important issue, and the NaV1.4-GIIIA system provides a good opportunity for such a comparison. To this end, we study the binding of GIIIA to the bacterial channels NaVAb and NaVRh. The complex structures are obtained using docking and molecular dynamics simulations, and the dissociation of GIIIA is studied through umbrella sampling simulations. The results are compared to those obtained from the NaV1.4-GIIIA system, and the differences in the binding modes arising from the changes in the selectivity filters are highlighted.

Constant pH Molecular Dynamics Reveals pH-Modulated Binding of Two Small-Molecule BACE1 Inhibitors

Targeting β-secretase (BACE1) with small-molecule inhibitors offers a promising route for treatment of Alzheimer's disease. However, the intricate pH dependence of BACE1 function and inhibitor efficacy has posed major challenges for structure-based drug design. Here we investigate two structurally similar BACE1 inhibitors that have dramatically different inhibitory activity using continuous constant pH molecular dynamics (CpHMD). At high pH, both inhibitors are stably bound to BACE1; however, within the enzyme active pH range, only the iminopyrimidinone-based inhibitor remains bound, while the aminothiazine-based inhibitor becomes partially dissociated following the loss of hydrogen bonding with the active site and change of the 10s loop conformation. The drastically lower activity of the second inhibitor is due to the protonation of a catalytic aspartate and the lack of a propyne tail. This work demonstrates that CpHMD can be used for screening pH-dependent binding profiles of small-molecule inhibitors, providing a new tool for structure-based drug design and optimization.

De novo reconstruction of DNA origami structures through atomistic molecular dynamics simulation

The DNA origami method has brought nanometer-precision fabrication to molecular biology labs, offering myriads of potential applications in the fields of synthetic biology, medicine, molecular computation, etc. Advancing the method further requires controlling self-assembly down to the atomic scale. Here we demonstrate a computational method that allows the equilibrium structure of a large, complex DNA origami object to be determined to atomic resolution. Through direct comparison with the results of cryo-electron microscopy, we demonstrate de novo reconstruction of a 4.7 megadalton pointer structure by means of fully atomistic molecular dynamics simulations. Furthermore, we show that elastic network-guided simulations performed without solvent can yield similar accuracy at a fraction of the computational cost, making this method an attractive approach for prototyping and validation of self-assembled DNA nanostructures.

The pathway of ligand entry from the membrane bilayer to a lipid G protein-coupled receptor

The binding process through the membrane bilayer of lipid-like ligands to a protein target is an important but poorly explored recognition process at the atomic level. In this work we succeeded in resolving the binding of the lipid inhibitor ML056 to the sphingosine-1-phosphate receptor 1 (S1P1R) using unbiased molecular dynamics simulations with an aggregate sampling of over 800 μs. The binding pathway is a multi-stage process consisting of the ligand diffusing in the bilayer leaflet to contact a "membrane vestibule" at the top of TM 7, subsequently moving from this lipid-facing vestibule to the orthosteric binding cavity through a channel formed by TMs 1 and 7 and the N-terminal of the receptor. Unfolding of the N-terminal alpha-helix increases the volume of the channel upon ligand entry, helping to reach the crystallographic pose that also corresponds to the predicted favorable pose. The relaxation timescales of the binding process show that the binding of the ligand to the "membrane vestibule" is the rate-limiting step in the multi microseconds timescale. We comment on the significance and parallels of the binding process in the context of other binding studies.

Effects of cytosine modifications on DNA flexibility and nucleosome mechanical stability

Cytosine can undergo modifications, forming 5-methylcytosine (5-mC) and its oxidized products 5-hydroxymethylcytosine (5-hmC), 5-formylcytosine (5-fC) and 5-carboxylcytosine (5-caC). Despite their importance as epigenetic markers and as central players in cellular processes, it is not well understood how these modifications influence physical properties of DNA and chromatin. Here we report a comprehensive survey of the effect of cytosine modifications on DNA flexibility. We find that even a single copy of 5-fC increases DNA flexibility markedly. 5-mC reduces and 5-hmC enhances flexibility, and 5-caC does not have a measurable effect. Molecular dynamics simulations show that these modifications promote or dampen structural fluctuations, likely through competing effects of base polarity and steric hindrance, without changing the average structure. The increase in DNA flexibility increases the mechanical stability of the nucleosome and vice versa, suggesting a gene regulation mechanism where cytosine modifications change the accessibility of nucleosomal DNA through their effects on DNA flexibility.

Chiral Graphene Quantum Dots

Chiral nanostructures from metals and semiconductors attract wide interest as components for polarization-enabled optoelectronic devices. Similarly to other fields of nanotechnology, graphene-based materials can greatly enrich physical and chemical phenomena associated with optical and electronic properties of chiral nanostructures and facilitate their applications in biology as well as other areas. Here, we report that covalent attachment of l/d-cysteine moieties to the edges of graphene quantum dots (GQDs) leads to their helical buckling due to chiral interactions at the "crowded" edges. Circular dichroism (CD) spectra of the GQDs revealed bands at ca. 210-220 and 250-265 nm that changed their signs for different chirality of the cysteine edge ligands. The high-energy chiroptical peaks at 210-220 nm correspond to the hybridized molecular orbitals involving the chiral center of amino acids and atoms of graphene edges. Diverse experimental and modeling data, including density functional theory calculations of CD spectra with probabilistic distribution of GQD isomers, indicate that the band at 250-265 nm originates from the three-dimensional twisting of the graphene sheet and can be attributed to the chiral excitonic transitions. The positive and negative low-energy CD bands correspond to the left and right helicity of GQDs, respectively. Exposure of liver HepG2 cells to L/D-GQDs reveals their general biocompatibility and a noticeable difference in the toxicity of the stereoisomers. Molecular dynamics simulations demonstrated that d-GQDs have a stronger tendency to accumulate within the cellular membrane than L-GQDs. Emergence of nanoscale chirality in GQDs decorated with biomolecules is expected to be a general stereochemical phenomenon for flexible sheets of nanomaterials.

Molecular Insights into the Adsorption Mechanism of Human β-Defensin-3 on Bacterial Membranes

Human β-defensin-3 (hBD3) is an endogenous antimicrobial peptide that exhibits broad-spectrum antibacterial activity without eukaryotic cytotoxicity. In this work, we carried out molecular dynamics (MD) simulations to explore its adsorption mechanism on, and the structural and thermodynamic contributions of individual residues to its antibacterial activity with both Gram-negative (GN) and Gram-positive (GP) bacterial membrane. Due to the strong electrostatic interaction of hBD3 with POPG lipids, which are more prevalent on the GP membrane, its adhesion to the GP membrane is stronger than to the GN membrane and stabilized more rapidly. On the surface of both bacterial membranes, the orientation of hBD3 is dominated by an electric dipole. We next analyzed the binding free energy decompositions of the hBD3-membrane complex using the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The results of both the GN and the GP membrane simulations show that Arg17, Arg36, and Arg38 form both polar and nonpolar interactions and are potentially the key residues for hBD3 antibacterial activity. On the other hand, there was a significant difference in the energy contribution of Arg12 between the GP and GN membrane simulations, suggesting that Arg12 is a key factor in the toxicity of hBD3 to specifically GP bacteria. Our findings shed light on the antibacterial activity of hBD3 on bacterial membranes and yield insights useful for the design of potent antimicrobial peptides targeting multidrug resistant bacteria.

Comparing pharmacophore models derived from crystal structures and from molecular dynamics simulations

ABSTRACT

Pharmacophore modeling is a widely used technique in computer-aided drug discovery. Structure-based pharmacophore models of a ligand in complex with a protein have proven to be useful for supporting in silico hit discovery, hit to lead expansion, and lead optimization. As a structure-based approach it depends on the correct interpretation of ligand-protein interactions. There are legitimate concerns about the fidelity of the bound ligand and about non-physiological contacts with parts of the crystal and the solvent effects that influence the protein structure. A possible way to refine the structure of a protein-ligand system is to use the final structure of a given MD simulation. In this study we compare pharmacophore models built using the initial protein-ligand structure obtained from the protein data bank (PDB) with pharmacophore models built with the final structure of a molecular dynamics simulation. We show that the pharmacophore models differ in feature number and feature type and that the pharmacophore models built from the last structure of a MD simulation shows in some cases better ability to distinguish between active and decoy ligand structures.

GRAPHICAL ABSTRACT

C60 fullerene localization and membrane interactions in RAW 264.7 immortalized mouse macrophages

There continues to be a significant increase in the number and complexity of hydrophobic nanomaterials that are engineered for a variety of commercial purposes making human exposure a significant health concern. This study uses a combination of biophysical, biochemical and computational methods to probe potential mechanisms for uptake of C60 nanoparticles into various compartments of living immune cells. Cultures of RAW 264.7 immortalized murine macrophage were used as a canonical model of immune-competent cells that are likely to provide the first line of defense following inhalation. Modes of entry studied were endocytosis/pinocytosis and passive permeation of cellular membranes. The evidence suggests marginal uptake of C60 clusters is achieved through endocytosis/pinocytosis, and that passive diffusion into membranes provides a significant source of biologically-available nanomaterial. Computational modeling of both a single molecule and a small cluster of fullerenes predicts that low concentrations of fullerenes enter the membrane individually and produce limited perturbation; however, at higher concentrations the clusters in the membrane causes deformation of the membrane. These findings are bolstered by nuclear magnetic resonance (NMR) of model membranes that reveal deformation of the cell membrane upon exposure to high concentrations of fullerenes. The atomistic and NMR models fail to explain escape of the particle out of biological membranes, but are limited to idealized systems that do not completely recapitulate the complexity of cell membranes. The surprising contribution of passive modes of cellular entry provides new avenues for toxicological research that go beyond the pharmacological inhibition of bulk transport systems such as pinocytosis.

All the O2 Consumed by Thermus thermophilus Cytochrome ba3 Is Delivered to the Active Site through a Long, Open Hydrophobic Tunnel with Entrances within the Lipid Bilayer

Cytochrome ba3 is a proton-pumping heme-copper oxygen reductase from the extreme thermophile Thermus thermophilus. Despite the fact that the enzyme's active site is buried deep within the protein, the apparent second order rate constant for the initial binding of O2 to the active-site heme has been experimentally found to be 10(9) M(-1) s(-1) at 298 K, at or near the diffusion limit, and 2 orders of magnitude faster than for O2 binding to myoglobin. To provide quantitative and microscopic descriptions of the O2 delivery pathway and mechanism in cytochrome ba3, extensive molecular dynamics simulations of the enzyme in its membrane-embedded form have been performed, including different protocols of explicit ligand sampling (flooding) simulations with O2, implicit ligand sampling analysis, and in silico mutagenesis. The results show that O2 diffuses to the active site exclusively via a Y-shaped hydrophobic tunnel with two 25-Å long membrane-accessible branches that coincide with the pathway previously suggested by the crystallographically identified xenon binding sites. The two entrances of the bifurcated tunnel of cytochrome ba3 are located within the lipid bilayer, where O2 is preferentially partitioned from the aqueous phase. The largest barrier to O2 migration within the tunnel is estimated to be only 1.5 kcal/mol, allowing O2 to reach the enzyme active site virtually impeded by one-dimensional diffusion once it reaches a tunnel entrance at the protein surface. Unlike other O2-utilizing proteins, the tunnel is "open" with no transient barriers observed due to protein dynamics. This unique low-barrier passage through the protein ensures that O2 transit through the protein is never rate-limiting.

Design and evaluation of an apta-nano-sensor to detect Acetamiprid in vitro and in silico

Pesticide detection is a main concern of food safety experts. Therefore, it is urgent to design an accurate, rapid, and cheap test. Biosensors that detect pesticide residues could replace current methods, such as HPLC or GC-MC. This research designs a biosensor based on aptamer (Oligonucleotide ss-DNA) in the receptor role, silver nanoparticles (AgNPs) as optical sensors and salt (NaCl) as the aggregative inducer of AgNPs to detect the presence of Acetamiprid. After optimization, .6 μM aptamer and 100 mM salt were employed. The selectivity and sensitivity of the complex were examined by different pesticides and different Acetamiprid concentrations. To simulate in vitro experimental conditions, bioinformatics software was used as in silico analysis. The results showed the detection of Acetamiprid at the .02 ppm (89.8 nM) level in addition to selectivity. Docking outputs introduced two loops as active sites in aptamer and confirmed aptamer-Acetamiprid bonding. Circular dichroism spectroscopy (CD) confirmed upon Acetamiprid binding, aptamer was folded due to stem-loop formation. Stability of the Apt-Acetamiprid complex in a simulated aqueous media was examined by molecular dynamic studies.

Violation of DNA neighbor exclusion principle in RNA recognition

DNA intercalation has been very useful for engineering DNA-based functional materials.

DNA intercalation has been very useful for engineering DNA-based functional materials. It is generally expected that the intercalation phenomenon in RNA would be similar to that in DNA. Here we note that the neighbor-exclusion principle is violated in RNA by naphthalene-based cationic probes, in contrast to the fact that it is usually valid in DNA. All the intercalation structures are responsible for the fluorescence, where small naphthalene moieties are intercalated in between bases via π–π interactions. The structure is aided by hydrogen bonds between the cationic moieties and the ribose-phosphate backbone, which results in specific selectivity for RNA over DNA. This experimentally observed mechanism is supported by computationally reproducing the fluorescence and CD data. MD simulations confirm the unfolding of RNA due to the intercalation of probes. Elucidation of the mechanism of selective sensing for RNA over DNA would be highly beneficial for dynamical observation of RNA which is essential for studying its biological roles.

A computational assessment of the permeability and salt rejection of carbon nanotube membranes and their application to water desalination

Membranes made from nanomaterials such as nanotubes and graphene have been suggested to have a range of applications in water filtration and desalination, but determining their suitability for these purposes requires an accurate assessment of the properties of these novel materials. In this study, we use molecular dynamics simulations to determine the permeability and salt rejection capabilities for membranes incorporating carbon nanotubes (CNTs) at a range of pore sizes, pressures and concentrations. We include the influence of osmotic gradients and concentration build up and simulate at realistic pressures to improve the reliability of estimated membrane transport properties. We find that salt rejection is highly dependent on the applied hydrostatic pressure, meaning high rejection can be achieved with wider tubes than previously thought; while membrane permeability depends on salt concentration. The ideal size of the CNTs for desalination applications yielding high permeability and high salt rejection is found to be around 1.1 nm diameter. While there are limited energy gains to be achieved in using ultra-permeable CNT membranes in desalination by reverse osmosis, such membranes may allow for smaller plants to be built as is required when size or weight must be minimized. There are diminishing returns in further increasing membrane permeability, so efforts should focus on the fabrication of membranes containing narrow or functionalized CNTs that yield the desired rejection or selection properties rather than trying to optimize pore densities.

Computational studies of the binding profile of phosphoinositide PtdIns (3,4,5) P₃ with the pleckstrin homology domain of an oomycete cellulose synthase

Saprolegnia monoica is a model organism to investigate Saprolegnia parasitica, an important oomycete which causes considerable loss in aquaculture every year. S. monoica contains cellulose synthases vital for oomycete growth. However, the molecular mechanism of the cellulose biosynthesis process in the oomycete growth is still poorly understood. Some cellulose synthases of S. monoica, such as SmCesA2, are found to contain a plecsktrin homology (PH) domain, which is a protein module widely found in nature and known to bind to phosphoinositides, a class of signaling compounds involved in many biological processes. Understanding the molecular interactions between the PH domain and phosphoinositides would help to unravel the cellulose biosynthesis process of oomycetes. In this work, the binding profile of PtdIns (3,4,5) P₃, a typical phosphoinositide, with SmCesA2-PH was studied by molecular docking, molecular dynamics and metadynamics simulations. PtdIns (3,4,5) P₃ is found to bind at a specific site located at β1, β2 and β1-β2 loop of SmCesA2-PH. The high affinity of PtdIns (3,4,5) P₃ to SmCesA2-PH is contributed by the free phosphate groups, which have electrostatic and hydrogen-bond interactions with Lys88, Lys100 and Arg102 in the binding site.

DIRECT-ID: An automated method to identify and quantify conformational variations--application to β2 -adrenergic GPCR

The conformational dynamics of a macromolecule can be modulated by a number of factors, including changes in environment, ligand binding, and interactions with other macromolecules, among others. We present a method that quantifies the differences in macromolecular conformational dynamics and automatically extracts the structural features responsible for these changes. Given a set of molecular dynamics (MD) simulations of a macromolecule, the norms of the differences in covariance matrices are calculated for each pair of trajectories. A matrix of these norms thus quantifies the differences in conformational dynamics across the set of simulations. For each pair of trajectories, covariance difference matrices are parsed to extract structural elements that undergo changes in conformational properties. As a demonstration of its applicability to biomacromolecular systems, the method, referred to as DIRECT-ID, was used to identify relevant ligand-modulated structural variations in the β2 -adrenergic (β2 AR) G-protein coupled receptor. Micro-second MD simulations of the β2 AR in an explicit lipid bilayer were run in the apo state and complexed with the ligands: BI-167107 (agonist), epinephrine (agonist), salbutamol (long-acting partial agonist), or carazolol (inverse agonist). Each ligand modulated the conformational dynamics of β2 AR differently and DIRECT-ID analysis of the inverse-agonist vs. agonist-modulated β2 AR identified residues known through previous studies to selectively propagate deactivation/activation information, along with some previously unidentified ligand-specific microswitches across the GPCR. This study demonstrates the utility of DIRECT-ID to rapidly extract functionally relevant conformational dynamics information from extended MD simulations of large and complex macromolecular systems.

Structure-based design of N-substituted 1-hydroxy-4-sulfamoyl-2-naphthoates as selective inhibitors of the Mcl-1 oncoprotein

Structure-based drug design was utilized to develop novel, 1-hydroxy-2-naphthoate-based small-molecule inhibitors of Mcl-1. Ligand design was driven by exploiting a salt bridge with R263 and interactions with the p2 pocket of the protein. Significantly, target molecules were accessed in just two synthetic steps, suggesting further optimization will require minimal synthetic effort. Molecular modeling using the Site-Identification by Ligand Competitive Saturation (SILCS) approach was used to qualitatively direct ligand design as well as develop quantitative models for inhibitor binding affinity to Mcl-1 and the Bcl-2 relative Bcl-xL as well as for the specificity of binding to the two proteins. Results indicated hydrophobic interactions in the p2 pocket dominated affinity of the most favourable binding ligand (3bl: Ki = 31 nM). Compounds were up to 19-fold selective for Mcl-1 over Bcl-xL. Selectivity of the inhibitors was driven by interactions with the deeper p2 pocket in Mcl-1 versus Bcl-xL. The SILCS-based SAR of the present compounds represents the foundation for the development of Mcl-1 specific inhibitors with the potential to treat a wide range of solid tumours and hematological cancers, including acute myeloid leukemia.

The Location of the Protonated and Unprotonated Forms of Arbidol in the Membrane: A Molecular Dynamics Study

Abstract

Arbidol is a potent broad-spectrum antiviral molecule for the treatment and prophylaxis of many viral infections. Viruses that can be inhibited by arbidol include enveloped and non-enveloped viruses, RNA and DNA viruses, as well as pH-independent and pH-dependent ones. These differences in viral types highlight the broad spectrum of Arb antiviral activity and, therefore, it must affect a common viral critical step. Arbidol incorporates rapidly into biological membranes, and some of its antiviral effects might be related to its capacity to interact with and locate into the membrane. However, no information is available of the molecular basis of its antiviral mechanism/s. We have aimed to locate the protonated (Arp) and unprotonated (Arb) forms of arbidol in a model membrane system. Both Arb and Arp locate in between the hydrocarbon acyl chains of the phospholipids but its specific location and molecular interactions differ from each other. Whereas both Arb and Arp average location in the membrane palisade is a similar one, Arb tends to be perpendicular to the membrane surface, whereas Arp tends to be parallel to it. Furthermore, Arp, in contrast to Arb, seems to interact stronger with POPG than with POPC, implying the existence of a specific interaction between Arp, the protonated from, with negatively charged phospholipids. This data would suggest that the active molecule of arbidol in the membrane is the protonated one, i.e., the positively charged molecule. The broad antiviral activity of arbidol would be defined by the perturbation it exerts on membrane structure and therefore membrane functioning.

Graphical Abstract

Electronic supplementary material

The online version of this article (doi:10.1007/s00232-016-9876-3) contains supplementary material, which is available to authorized users.

Additive CHARMM force field for naturally occurring modified ribonucleotides

More than 100 naturally occurring modified nucleotides have been found in RNA molecules, in particular in tRNAs. We have determined molecular mechanics force field parameters compatible with the CHARMM36 all‐atom additive force field for all these modifications using the CHARMM force field parametrization strategy. Emphasis was placed on fine tuning of the partial atomic charges and torsion angle parameters. Quantum mechanics calculations on model compounds provided the initial set of target data, and extensive molecular dynamics simulations of nucleotides and oligonucleotides in aqueous solutions were used for further refinement against experimental data. The presented parameters will allow for computational studies of a wide range of RNAs containing modified nucleotides, including the ribosome and transfer RNAs. © 2016 The Authors. Journal of Computational Chemistry Published by Wiley Periodicals, Inc.

Pore dilatation increases the bicarbonate permeability of CFTR, ANO1 and glycine receptor anion channels

KEY POINTS

Cellular stimuli can modulate the ion selectivity of some anion channels, such as CFTR, ANO1 and the glycine receptor (GlyR), by changing pore size. Ion selectivity of CFTR, ANO1 and GlyR is critically affected by the electric permittivity and diameter of the channel pore. Pore size change affects the energy barriers of ion dehydration as well as that of size-exclusion of anion permeation. Pore dilatation increases the bicarbonate permeability (P HC O3/ Cl ) of CFTR, ANO1 and GlyR. Dynamic change in P HC O3/ Cl may mediate many physiological and pathological processes.

ABSTRACT

Chloride (Cl(-) ) and bicarbonate (HCO3 (-) ) are two major anions and their permeation through anion channels plays essential roles in our body. However, the mechanism of ion selection by the anion channels is largely unknown. Here, we provide evidence that pore dilatation increases the bicarbonate permeability (P HC O3/ Cl ) of anion channels by reducing energy barriers of size-exclusion and ion dehydration of HCO3 (-) permeation. Molecular, physiological and computational analyses of major anion channels, such as cystic fibrosis transmembrane conductance regulator (CFTR), anoctamin-1(ANO1/TMEM16A) and the glycine receptor (GlyR), revealed that the ion selectivity of anion channels is basically determined by the electric permittivity and diameter of the pore. Importantly, cellular stimuli dynamically modulate the anion selectivity of CFTR and ANO1 by changing the pore size. In addition, pore dilatation by a mutation in the pore-lining region alters the anion selectivity of GlyR. Changes in pore size affected not only the energy barriers of size exclusion but that of ion dehydration by altering the electric permittivity of water-filled cavity in the pore. The dynamic increase in P HC O3/ Cl by pore dilatation may have many physiological and pathophysiological implications ranging from epithelial HCO3 (-) secretion to neuronal excitation.

Ligand Binding Pathways of Clozapine and Haloperidol in the Dopamine D2 and D3 Receptors

The binding of a small molecule ligand to its protein target is most often characterized by binding affinity and is typically viewed as an on/off switch. The more complex reality is that binding involves the ligand passing through a series of intermediate states between the solution phase and the fully bound pose. We have performed a set of 29 unbiased molecular dynamics simulations to model the binding pathways of the dopamine receptor antagonists clozapine and haloperidol binding to the D2 and D3 dopamine receptors. Through these simulations we have captured the binding pathways of clozapine and haloperidol from the extracellular vestibule to the orthosteric binding site and thereby, we also predict the bound pose of each ligand. These are the first long time scale simulations of haloperidol or clozapine binding to dopamine receptors. From these simulations, we have identified several important stages in the binding pathway, including the involvement of Tyr7.35 in a "handover" mechanism that transfers the ligand between the extracellular vestibule and Asp3.32. We have also performed interaction and cluster analyses to determine differences in binding pathways between the D2 and D3 receptors and identified metastable states that may be of use in drug design.

Photo-switchable tweezers illuminate pore-opening motions of an ATP-gated P2X ion channel

P2X receptors function by opening a transmembrane pore in response to extracellular ATP. Recent crystal structures solved in apo and ATP-bound states revealed molecular motions of the extracellular domain following agonist binding. However, the mechanism of pore opening still remains controversial. Here we use photo-switchable cross-linkers as ‘molecular tweezers’ to monitor a series of inter-residue distances in the transmembrane domain of the P2X2 receptor during activation. These experimentally based structural constraints combined with computational studies provide high-resolution models of the channel in the open and closed states. We show that the extent of the outer pore expansion is significantly reduced compared to the ATP-bound structure. Our data further reveal that the inner and outer ends of adjacent pore-lining helices come closer during opening, likely through a hinge-bending motion. These results provide new insight into the gating mechanism of P2X receptors and establish a versatile strategy applicable to other membrane proteins.

DOI:http://dx.doi.org/10.7554/eLife.11050.001

Protein receptors in the cell membrane play an important role transmitting signals from outside to inside the cell. Members of the P2X family of receptors are ion channels that form pores through the membrane. When a molecule of ATP binds to the external region of the receptor, it activates it and causes the receptor to change from a closed to an open shape. Once opened, ions flow through the channel’s pore and trigger a response inside the cell. P2X receptors are found on most animal cells (including nerve cells) and are involved in both normal cellular activity and processes linked to disease, including inflammation and chronic pain.

The P2X receptor has three parts or subunits, and each contributes to the channel’s pore. Recent research using a technique called X-ray crystallography has revealed how ATP binding causes shape changes in the external region of the receptor. But these three-dimensional structures did not reveal details of how the subunits move to open or close the channel’s pore.

Habermacher et al. have now added light-sensitive linkers onto the P2X receptor in a way that meant that different colors of light could be used to force parts of the receptor to come closer together or move apart. This allowed the pore to be opened and closed in response to changes in light. Habermacher et al. then studied the behavior of these modified receptors within a natural membrane and found that the light stimulated movements were similar to those seen with ATP. When the behavior of the receptor and light-sensitive linkers was studied using computer simulations, it led to new models of the P2X pore in the open and closed state. In these models, the open channel was more tightly packed than in the previous structure and an unexpected hinge-bending movement was seen to accompany the opening of the channel. It is hoped that this new approach will also be useful for probing how other membrane proteins change their shape when activated.

DOI:http://dx.doi.org/10.7554/eLife.11050.002

1,3-Propanediol binds inside the water-conducting pore of aquaporin 4: Does this efficacious inhibitor have sufficient potency?

Among the thirteen types of water channel proteins, aquaporins (AQPs), which play various essential roles in human physiology, AQP4 is richly expressed in cells of the central nervous system and implicated in pathological conditions such as brain edema. Therefore, researchers have been looking for ways to inhibit AQP4's water-conducting function. Many small molecules have been investigated for their interactions with the residues that form the AQP4 channel entry vestibule on the extracellular side and their interruption of waters entering into the conducting pore. Conducting all-atom simulations on the basis of CHARMM 36 force field, we study one such inhibitor, 5-acetamido-1,3,4-thiadiazole-2-sulfonamide (AZM), to achieve quantitative agreement between the computed and the experimentally measured values of AZM-AQP4 binding affinity. Using the same method, we examine the possibility of plugging up the AQP4 channel around the Asn-Pro-Ala motifs located near the channel center because a small molecule bound there would totally occlude water conduction through AQP4. We compute the binding affinities of 1,2-ethanediol (EDO) and 1,3-propanediol (PDO) inside the AQP4 conducting pore and identify the specificities of the interactions. The EDO-AQP4 interaction is weak with a dissociation constant of 80 mM. The PDO-AQP4 interaction is rather strong with a dissociation constant of 328 μM, which indicates that PDO is an efficacious AQP4 inhibitor with sufficiently high potency. Considering the fact that PDO is classified by the US Food and Drug Administration as generally safe, we predict that 1,3-propanediol could be an effective drug for brain edema and other AQP4-correlated neurological conditions.

Small Molecule Inhibitors of Ca(2+)-S100B Reveal Two Protein Conformations

The drug pentamidine inhibits calcium-dependent complex formation with p53 ((Ca)S100B·p53) in malignant melanoma (MM) and restores p53 tumor suppressor activity in vivo. However, off-target effects associated with this drug were problematic in MM patients. Structure-activity relationship (SAR) studies were therefore completed here with 23 pentamidine analogues, and X-ray structures of (Ca)S100B·inhibitor complexes revealed that the C-terminus of S100B adopts two different conformations, with location of Phe87 and Phe88 being the distinguishing feature and termed the "FF-gate". For symmetric pentamidine analogues ((Ca)S100B·5a, (Ca)S100B·6b) a channel between sites 1 and 2 on S100B was occluded by residue Phe88, but for an asymmetric pentamidine analogue ((Ca)S100B·17), this same channel was open. The (Ca)S100B·17 structure illustrates, for the first time, a pentamidine analog capable of binding the "open" form of the "FF-gate" and provides a means to block all three "hot spots" on (Ca)S100B, which will impact next generation (Ca)S100B·p53 inhibitor design.

Iron depletion strategy for targeted cancer therapy: utilizing the dual roles of neutrophil gelatinase-associated lipocalin protein

Decreasing iron uptake and increasing iron efflux may result in cell death by oxidative inactivation of vital enzymes. Applying the dual function of neutrophil gelatinase-associated lipocalin (NGAL) could achieve the goal of iron depletion in the cancer cells. Tyr106, Lys125 or Lys134 was the key binding site for NGAL protein to sequester iron-chelating siderophores. In this study, we employed all bioactive peptides in peptide databank to dock with the siderophore-binding sites of NGAL protein by virtual screening. In addition, we performed molecular dynamics (MD) simulation to observe the molecular character and structural variation of ligand-protein interaction. Glu-Glu-Lys-Glu (EEKE), Glu-Glu-Asp-Cys-Lys (EEDCK), and Gly-Glu-Glu-Cys-Asp (GEECD) were selected preliminarily by rigorous scoring functions for further investigation. GEECD was excluded due to higher binding total energy than the others. Moreover, we also excluded EEKE due to larger influence to the stability of binding residues by the information of root mean square fluctuation (RMSF) and principal component analysis (PCA). Thus, we suggested that EEDCK was the potential bioactive peptide which had been proved to inhibit malignant cells for targeted cancer therapy. Graphical Abstract Perspective drug design of occupying the siderophore-binding sites of NGAL outside the cell temporarily by a potential short peptide until NGAL enters into the cell, and releasing the siderophore-binding sites inside the cell.

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.

Molecular switches of the κ opioid receptor triggered by 6'-GNTI and 5'-GNTI

The κ opioid receptor (κOR) is a member of G-protein-coupled receptors, and is considered as a promising drug target for treating neurological diseases. κOR selective 6'-GNTI was proved to be a G-protein biased agonist, whereas 5'-GNTI acts as an antagonist. To investigate the molecular mechanism of how these two ligands induce different behaviors of the receptor, we built two systems containing the 5'-GNTI-κOR complex and the 6'-GNTI-κOR complex, respectively, and performed molecular dynamics simulations of the two systems. We observe that transmembrane (TM) helix 6 of the κOR rotates about 4.6(o) on average in the κOR-6'-GNTI complex. Detailed analyses of the simulation results indicate that E297(6.58) and I294(6.55) play crucial roles in the rotation of TM6. In the simulation of the κOR-5'-GNTI system, it is revealed that 5'-GNTI can stabilize TM6 in the inactive state form. In addition, the kink of TM7 is stabilized by a hydrogen bond between S324(7.47) and the residue V69(1.42) on TM1.

Mechanism of Hydrotropic Action of Hydrotrope Sodium Cumene Sulfonate on the Solubility of Di-t-Butyl-Methane: A Molecular Dynamics Simulation Study

Hydrotropes are special class of amphiphilic molecules that have an ability to solubilize the insoluble or sparingly soluble molecules in water. To find out the mechanism of hydrotropic action of hydrotropes on hydrophobic molecules, we have carried out classical molecular dynamics simulation of hydrophobic solute di-t-butyl-methane (DTBM) and hydrotrope sodium cumene sulfonate (SCS) in water with a regime of SCS concentrations. Our study demonstrates that, above the minimum hydrotrope concentration (MHC), the self-aggregation of SCS starts, and it creates a micellar-like environment in which the hydrophobic tail part of SCS points inward while its hydrophilic sulfonate group points outward to make favorable contact with water molecules. The formation of the hydrophobic core of SCS cluster creates a hydrophobic environment where the hydrophobic DTBM molecules are encapsulated. Interestingly, the determination of average water-SCS hydrogen bonds further suggests that the aggregate formation of SCS molecules has a negligible influence on it. Moreover, the calculations of Flory-Huggins interaction parameters also reveal favorable interactions between hydrotrope SCS and solute DTBM molecules. The implications of these findings on the mechanism of hydrotrope assisted enhanced solubility of hydrophobic molecules are discussed.

Allosteric inhibitors have distinct effects, but also common modes of action, in the HCV polymerase

The RNA-dependent RNA polymerase from the Hepatitis C Virus (gene product NS5B) is a validated drug target because of its critical role in genome replication. There are at least four distinct allosteric sites on the polymerase to which several small molecule inhibitors bind. In addition, numerous crystal structures have been solved with different allosteric inhibitors bound to the polymerase. However, the molecular mechanisms by which these small molecules inhibit the enzyme have not been fully elucidated. There is evidence that allosteric inhibitors alter the intrinsic motions and distribution of conformations sampled by the enzyme. In this study we use molecular dynamics simulations to understand the structural and dynamic changes that result when inhibitors are bound at three different allosteric binding sites on the enzyme. We observe that ligand binding at each site alters the structure and dynamics of NS5B in a distinct manner. Nonetheless, our studies also highlight commonalities in the mechanisms of action of the different inhibitors. Each inhibitor alters the conformational states sampled by the enzyme, either by rigidifying the enzyme and preventing transitions between functional conformational states or by destabilizing the enzyme and preventing functionally relevant conformations from being adequately sampled. By illuminating the molecular mechanisms of allosteric inhibition, these studies delineate the intrinsic functional properties of the enzyme and pave the way for designing novel and more effective polymerase inhibitors. This information may also be important to understand how allosteric regulation occurs in related viral polymerases and other enzymes.