Selectivity principle of the ligand escape process from a two-gate tunnel in myoglobin: molecular dynamics simulation

We proposed a selectivity principle for the ligand escape process from two fluctuating bottlenecks in a cavity with a multigate inside a myoglobin pocket. Our previous analytical theory proposed a fluctuating bottleneck model for a Brownian particle passing through two gates on a cavity surface of an enzyme protein and has determined the escape rate in terms of the time-dependent gate function and the competition effect. It illustrated that with two (or more than two) gates on a cavity surface the gate modulation, which is controlled by protein fluctuation, dominates the ligand escape pathway. We have performed a molecular dynamics simulation to investigate the selectivity principle of the ligand escape process from two-gate tunnel in myoglobin. The simulation results confirm our theoretical conjecture. It indicates that the escape process is actually entropy driven, and the ligand escape pathway is chosen via the gate modulation. This suggests an interesting intrinsic property, that is, the oxymyoglobin tertiary structure is favorable to the departure of the ligand from one direction rather than through a biased random walk.     

Two-state folding observed in individual protein molecules

The folding dynamics of small proteins are often described in terms of a simple two-state kinetic model. Within this notion, the behavior of individual molecules is expected to be stochastic, with a protein molecule residing in either the unfolded or the folded state for extended periods of time, with intermittent rapid jumps across the free energy barrier. However, a direct observation of this bistable behavior has not been made to date. Rather, previous reports of folding trajectories of individual proteins have shown an unexpected degree of complexity. This raises the question whether the simple kinetic properties derived from classical experiments on large ensembles of molecules are reflected in the folding paths taken by individual proteins. Here we report single-molecule folding/unfolding trajectories observed by fluorescence resonance energy transfer for a protein that meets all criteria of a two state-system. The trajectories, measured on molecules immobilized in lipid vesicles, demonstrate the anticipated bistable behavior, with steplike transitions between folded and unfolded conformations. They further allow us to put an upper bound on the barrier crossing time.     

甘油水溶液氢键特性的分子动力学模拟

为了研究低温保护剂溶液的结构和物理化学特性, 以甘油为保护剂, 采用分子动力学方法, 对不同浓度的甘油和水的二元体系进行了模拟. 得到了不同浓度的甘油水溶液在2 ns内的分子动力学运动轨迹, 通过对后1 ns内运动轨迹的分析, 得到了各个原子对的径向分布函数和甘油分子的构型分布. 根据氢键的图形定义, 分析了氢键的结构和动力学特性. 计算了不同浓度下体系中平均每个原子(O和H)和分子(甘油和水)参与氢键个数的百分比分布及其平均值. 同时还计算了所有氢键、水分子之间的氢键以及甘油与水分子之间的氢键的生存周期.     

Modeling heme proteins using atomistic simulations

Heme proteins are found in all living organisms, and perform a wide variety of tasks ranging from electron transport, to the oxidation of organic compounds, to the sensing and transport of small molecules. In this work we review the application of classical and quantum-mechanical atomistic simulation tools to the investigation of several relevant issues in heme proteins chemistry: (i) conformational analysis, ligand migration, and solvation effects studied using classical molecular dynamics simulations; (ii) electronic structure and spin state energetics of the active sites explored using quantum-mechanics (QM) methods; (iii) the interaction of heme proteins with small ligands studied through hybrid quantum mechanics-molecular mechanics (QM-MM) techniques; (iv) and finally chemical reactivity and catalysis tackled by a combination of quantum and classical tools.     

Cooperative motions of protein and hydration water molecules: molecular dynamics study of scytalone dehydratase

Two molecular dynamics (MD) simulations totaling 25 ns of simulation time of monomeric scytalone dehydratase (SD) were performed. The enzyme has a ligand-binding pocket containing a cone-shaped alpha+beta barrel, and the C-terminal region covers the binding pocket. Our simulations clarified the difference in protein dynamics and conformation between the liganded protein and the unliganded protein. The liganded protein held the ligand molecule tightly and the initial structure was maintained during the simulation. The unliganded protein, on the other hand, fluctuated dynamically and its structure changed largely from the initial structure. In the equilibrium state, the binding pocket was fully solvated by opening of the C-terminal region, and the protein dynamics was connected with hydration water molecules entry into and release from the binding pocket. In addition, the cooperative motions of the unliganded protein and the hydration water molecules produced the path through the protein interior for ligand binding.     

Multiple binding modes for dicationic Hoechst 33258 to DNA

The binding of dicationic Hoechst 33258 (ligand) to DNA was characterized by means of the fluorescence spectra, fluorescence intensity titration, time-resolved fluorescence decay, light scattering, circular dichroism, and fluorescence thermal denaturation measurements, and two binding modes were distinguished by the experimental results. Type 1 binding has the stoichiometry of one ligand to more than 12 base pairs, and it is defined as quasi-minor groove binding which has the typical prolonged fluorescence lifetime of about 4.4 ns. In type 1 binding, planar conformation of the ligand is favorable. Type 2 binding with phosphate to ligand ratio (P/L) < 2.5 has the stoichiometry of one ligand to two phosphates. It is defined as a highly dense and orderly stacked binding with DNA backbone as the template. Electrostatic interactions between doubly protonated ligands and negatively charged DNA backbone play a predominant role in the type 2 binding mode. The characteristics of this type of binding result in a twisted conformation of the ligand that has a fluorescence lifetime of less than 1 ns. The results also indicate that the binding is in a cooperative manner primarily by stacking of the aromatic rings of the neighboring ligands. Type 1 binding is only observed for double-stranded DNA (dsDNA) with affinity constant of 1.83 x 10(7) M-1. In the type 2 binding mode, the binding affinity constants are 4.9 x 10(6) and 4.3 x 10(6) M-1 for dsDNA and single-stranded DNA (ssDNA), respectively. The type 2 binding is base pair independent while the type 1 binding is base pair related. The experiments described in this paper revealed that the dication bindings are different from the monocation bindings reported by previous study. The dication binding leads to stronger aggregation at low ligand concentration and results in orderly arrangements of the ligands along DNA chains. Furthermore the dication binding is demonstrated to be beneficial for enhancing the DNA's stability.     

Nonequilibrium dynamics simulations of nitric oxide release: comparative study of nitrophorin and myoglobin

Nitrophorin 4 (NP4) is a heme protein that reversibly binds nitric oxide (NO), with release rates modulated by pH change. High-resolution structures of NP4 revealed that pH changes and NO binding induce a large conformational rearrangement in two loops that serve to protect the heme-bound NO molecule from solvent. We used extended (110 ns) molecular dynamics simulations of NP4 at pH 5 and pH 7, modeled by selective deprotonation of acidic groups. Conformational and dynamic changes were observed, consistent with those found in the crystal. Further, major solvent movement and NO escape were observed at pH 7, while the ligand remained in the heme binding pocket at pH 5. As a control, we also performed molecular dynamics (MD) simulations of sperm whale myoglobin, where NO migration into the interior cavities of the protein was observed, consistent with previous reports. We constructed a kinetic model of ligand escape to quantitatively relate the microscopic rate constants to the observed rates, and tested the predictions against the experimental data. The results suggest that release rates of diatomic molecules from heme proteins can be varied by several orders of magnitude through modest adjustments in geminate rebinding and gating behavior.     

Molecular dynamics simulation of entropy driven ligand escape process in heme pocket

Molecular dynamics simulations were performed to investigate the gate effect of protein motion on the escape of O(2) from the heme pocket. The existing geometric entropy in a spherical cavity pushes the ligand toward the cavity surface, and then the ligand escape along the cavity surface is controlled by the gate size and gate modulation, i.e., protein dynamics regulate the gating behavior, which is an inherent feature of proteins such as myoglobin. Our simulation results confirm that the ligand escape process is basically entropy driven.     

How does a drug molecule find its target binding site?

Although the thermodynamic principles that control the binding of drug molecules to their protein targets are well understood, detailed experimental characterization of the process by which such binding occurs has proven challenging. We conducted relatively long, unguided molecular dynamics simulations in which a ligand (the cancer drug dasatinib or the kinase inhibitor PP1) was initially placed at a random location within a box that also contained a protein (Src kinase) to which that ligand was known to bind. In several of these simulations, the ligand correctly identified its target binding site, forming a complex virtually identical to the crystallographically determined bound structure. The simulated trajectories provide a continuous, atomic-level view of the entire binding process, revealing persistent and noteworthy intermediate conformations and shedding light on the role of water molecules. The technique we employed, which does not assume any prior knowledge of the binding site's location, may prove particularly useful in the development of allosteric inhibitors that target previously undiscovered binding sites.     

Statistical potential for modeling and ranking of protein-ligand interactions

Applications in structural biology and medicinal chemistry require protein-ligand scoring functions for two distinct tasks: (i) ranking different poses of a small molecule in a protein binding site and (ii) ranking different small molecules by their complementarity to a protein site. Using probability theory, we developed two atomic distance-dependent statistical scoring functions: PoseScore was optimized for recognizing native binding geometries of ligands from other poses and RankScore was optimized for distinguishing ligands from nonbinding molecules. Both scores are based on a set of 8,885 crystallographic structures of protein-ligand complexes but differ in the values of three key parameters. Factors influencing the accuracy of scoring were investigated, including the maximal atomic distance and non-native ligand geometries used for scoring, as well as the use of protein models instead of crystallographic structures for training and testing the scoring function. For the test set of 19 targets, RankScore improved the ligand enrichment (logAUC) and early enrichment (EF(1)) scores computed by DOCK 3.6 for 13 and 14 targets, respectively. In addition, RankScore performed better at rescoring than each of seven other scoring functions tested. Accepting both the crystal structure and decoy geometries with all-atom root-mean-square errors of up to 2 ? from the crystal structure as correct binding poses, PoseScore gave the best score to a correct binding pose among 100 decoys for 88% of all cases in a benchmark set containing 100 protein-ligand complexes. PoseScore accuracy is comparable to that of DrugScore(CSD) and ITScore/SE and superior to 12 other tested scoring functions. Therefore, RankScore can facilitate ligand discovery, by ranking complexes of the target with different small molecules; PoseScore can be used for protein-ligand complex structure prediction, by ranking different conformations of a given protein-ligand pair. The statistical potentials are available through the Integrative Modeling Platform (IMP) software package (http://salilab.org/imp) and the LigScore Web server (http://salilab.org/ligscore/).     

The effect of point mutation on the equilibrium structural fluctuations of ferric Myoglobin

The ultrafast equilibrium fluctuations of the Fe(III)-NO complex of a single point mutation of Myoglobin (H64Q) have been studied using Fourier Transform 2D-IR spectroscopy. Comparison with data from wild type Myoglobin (wt-Mb) shows the presence of two conformational substates of the mutant haem pocket where only one exists in the wild type form. One of the substates of the mutant exhibits an almost identical NO stretching frequency and spectral diffusion dynamics to wt-Mb while the other is distinctly different in both respects. The remarkably contrasting dynamics are largely attributable to interactions between the NO ligand and a nearby distal side chain which provides a basis for understanding the roles of these side chains in other ferric haem proteins.     

Lanthanide(III) complexes of novel mixed carboxylic-phosphorus acid derivatives of diethylenetriamine: a step towards more efficient MRI contrast agents

Three novel phosphorus-containing analogues of H(5)DTPA (DTPA = diethylenetriaminepentaacetate) were synthesised (H6L1, H5L2, H5L3). These compounds have a -CH2-P(O)(OH)-R function (R = OH, Ph, CH2NBn2) attached to the central nitrogen atom of the diethylenetriamine backbone. An NMR study reveals that these ligands bind to lanthanide(III) ions in an octadentate fashion through the three nitrogen atoms, a P-O oxygen atom and four carboxylate oxygen atoms. The complexed ligand occurs in several enantiomeric forms due to the chirality of the central nitrogen atom and the phosphorus atom upon coordination. All lanthanide complexes studied have one coordinated water molecule. The residence times (tau(M)298) of the coordinated water molecules in the gadolinium(III) complexes of H6L1 and H5L2 are 88 and 92 ns, respectively, which are close to the optimum. This is particularly important upon covalent and noncovalent attachment of these Gd(3+) chelates to polymers. The relaxivity of the complexes studied is further enhanced by the presence of at least two water molecules in the second coordination sphere of the Gd(3+) ion, which are probably bound to the phosphonate/phosphinate moiety by hydrogen bonds. The complex [Gd(L3)(H2O)](2-) shows strong binding ability to HSA, and the adduct has a relaxivity comparable to MS-325 (40 s(-1) mM(-1) at 40 MHz, 37 degrees C) even though it has a less favourable tau(M) value (685 ns). Transmetallation experiments with Zn(2+) indicate that the complexes have a kinetic stability that is comparable to-or better than-those of [Gd(dtpa)(H2O)](2-) and [Gd(dtpa-bma)(H2O)].     

Effective simulations of gas diffusion through kinetically accessible tunnels in multisubunit proteins: O2 pathways and escape routes in T-state deoxyhemoglobin

The diffusion of small gases to special binding sites within polypeptide matrices pivotally defines the biochemical specificity and reactivity of proteins. We investigate here explicit O(2) diffusion in adult human hemoglobin (HbA) as a case study employing the recently developed temperature-controlled locally enhanced sampling (TLES) method and vary the parameters to greatly increase the simulation efficiency. The method is carefully validated against standard molecular dynamics (MD) simulations and available experimental structural and kinetic data on ligand diffusion in T-state deoxyHbA. The methodology provides a viable alternative approach to traditional MD simulations and/or potential of mean force calculations for: (i) characterizing kinetically accessible diffusion tunnels and escape routes for light ligands in porous proteins; (ii) very large systems when realistic simulations require the inclusion of multiple subunits of a protein; and (iii) proteins that access short-lived conformations relative to the simulation time. In the case of T-state deoxyHbA, we find distinct ligand diffusion tunnels consistent with the experimentally observed disparate Xe cavities in the α- and β-subunits. We identify two distal barriers including the distal histidine (E7) that control access to the heme. The multiple escape routes uncovered by our simulations call for a review of the current popular hypothesis on ligand escape from hemoglobin. Larger deviations from the crystal structure during simulated diffusion in isolated α- and β-subunits highlight the dampening effects of subunit interactions and the importance of including all subunits of multisubunit proteins to map realistic kinetically accessible diffusion tunnels and escape routes.     

Coarse‐grained molecular dynamics simulations of protein–ligand binding

Coarse‐grained molecular dynamics (CGMD) simulations with the MARTINI force field were performed to reproduce the protein–ligand binding processes. We chose two protein–ligand systems, the levansucrase–sugar (glucose or sucrose), and LinB–1,2‐dichloroethane systems, as target systems that differ in terms of the size and shape of the ligand‐binding pocket and the physicochemical properties of the pocket and the ligand. Spatial distributions of the Coarse‐grained (CG) ligand molecules revealed potential ligand‐binding sites on the protein surfaces other than the real ligand‐binding sites. The ligands bound most strongly to the real ligand‐binding sites. The binding and unbinding rate constants obtained from the CGMD simulation of the levansucrase–sucrose system were approximately 10 times greater than the experimental values; this is mainly due to faster diffusion of the CG ligand in the CG water model. We could obtain dissociation constants close to the experimental values for both systems. Analysis of the ligand fluxes demonstrated that the CG ligand molecules entered the ligand‐binding pockets through specific pathways. The ligands tended to move through grooves on the protein surface. Thus, the CGMD simulations produced reasonable results for the two different systems overall and are useful for studying the protein–ligand binding processes. © 2014 Wiley Periodicals, Inc.     

Exploring inhibitor release pathways in histone deacetylases using random acceleration molecular dynamics simulations

Molecular channel exploration perseveres to be the prominent solution for eliciting structure and accessibility of active site and other internal spaces of macromolecules. The volume and silhouette characterization of these channels provides answers for the issues of substrate access and ligand swapping between the obscured active site and the exterior of the protein. Histone deacetylases (HDACs) are metal-dependent enzymes that are involved in the cell growth, cell cycle regulation, and progression, and their deregulations have been linked with different types of cancers. Hence HDACs, especially the class I family, are widely recognized as the important cancer targets, and the characterizations of their structures and functions have been of special interest in cancer drug discovery. The class I HDACs are known to possess two different protein channels, an 11 ? and a 14 ? (named channels A and B1, respectively), of which the former is a ligand or substrate occupying tunnel that leads to the buried active site zinc ion and the latter is speculated to be involved in product release. In this work, we have carried out random acceleration molecular dynamics (RAMD) simulations coupled with the classical molecular dynamics to explore the release of the ligand, N-(2-aminophenyl) benzamide (LLX) from the active sites of the recently solved X-ray crystal structure of HDAC2 and the computationally modeled HDAC1 proteins. The RAMD simulations identified significant structural and dynamic features of the HDAC channels, especially the key 'gate-keeping' amino acid residues that control these channels and the ligand release events. Further, this study identified a novel and unique channel B2, a subchannel from channel B1, in the HDAC1 protein structure. The roles of water molecules in the LLX release from the HDAC1 and HDAC2 enzymes are also discussed. Such structural and dynamic properties of the HDAC protein channels that govern the ligand escape reactions will provide further mechanistic insights into the HDAC enzymes, which, in the long run, have a potential to bring new ideas for developing more promising HDAC inhibitors as well as extend our atomic level understandings on their mechanisms of action.     

Versatile coordination behaviour of an asymmetric half-salen ligand bearing a dansyl fluorophore

The coordinative chemistry of the tridentate half-salen ligand 5-(dimethylamino)-N-(2-((2-hydroxybenzylidene)amino)phenyl)naphthalene-1-sulfonamide (H(2)L, 1) has been studied by means of an electrochemical method. All of the complexes have been characterised using analytical and spectroscopic techniques. Ligand 1 and two nickel (6 and 7), copper (9), zinc (12) and cadmium (14) metal complexes have been studied by crystallography. Complexes 6 and 7 are octahedral and tetrahedral nickel(ii) complexes, respectively, and both contain an [L](2-) molecule that behaves in an [N(2)O] tridentate manner. Nickel(ii) completes its coordination kernel with three water molecules in complex 6, whereas in complex 7 the nickel ion is further bound to a molecule of dansylamine arising from a hydrolysis process. The copper(ii) complex 9 is a monomeric compound that contains a bideprotonated ligand thread and a dimethylsulfoxide molecule coordinated through the sulfur atom. The zinc complex 12 is an unusual pentanuclear cluster compound whose structure consists of four anionic ligand units and two hydroxo anions bound to five zinc(ii) centres. The appearance of the hydroxo anions in this complex provides new evidence for water reduction electrochemically promoted by zinc metal under mild conditions. The cadmium complex 14 is a dimeric compound that comprises two molecules of the anionic ligand and two dimethylsulfoxide molecules. The great structural variety exhibited by all these complexes demonstrates that the introduction of asymmetry in a salen skeleton by incorporating a dansyl pendant increases the versatility of the resulting ligand on coordination. All complexes are luminescent in solution at room temperature in acetonitrile solutions.     

Protein folding pathways extracted by OFLOOD: Outlier FLOODing method

The Outlier FLOODing method (OFLOOD) is proposed as an efficient conformational sampling method to extract biologically rare events such as protein folding. In OFLOOD, sparse distributions (outliers in the conformational space) were regarded as relevant states for the transitions. Then, the transitions were enhanced through conformational resampling from the outliers. This evidence indicates that the conformational resampling of the sparse distributions might increase chances for promoting the transitions from the outliers to other meta‐stable states, which resembles a conformational flooding from the outliers to the neighboring clusters. OFLOOD consists of (i) detections of outliers from conformational distributions and (ii) conformational resampling from the outliers by molecular dynamics (MD) simulations. Cycles of (i) and (ii) are simply repeated. As demonstrations, OFLOOD was applied to folding of Chignolin and HP35. In both cases, OFLOOD automatically extracted folding pathways from unfolded structures with ns‐order computational costs, although µs‐order canonical MD failed to extract them. © 2014 Wiley Periodicals, Inc.     

Transfer of photoenergy in pi-conjugated polymers. Two types of photoluminescence that involve energy transfer along a polymer chain

Two types of energy transfer in pi-conjugated polymers have been investigated using time-resolved photoluminescence (PL) techniques: type i, perpendicular-type energy transfer from the 2,3-di(p-tolyl)quinoxaline unit to the pi-conjugated main chain of poly[2,3-di(p-tolyl)quinoxaline-5,8-diyl], and type ii, parallel-type energy transfer from the oligo(pyridine-2,5-diyl) (O-Py) unit to the oligo(selenophene-2,5-diyl) (O-Se) unit in a block-type copolymer of O-Py and O-Se. Both types of energy transfer were very fast with a time constant shorter than approximately 0.1 ns; in particular, the type ii energy transfer took place with a time constant of approximately 5 ps. Both pi-conjugated polymers were considered to contain segments with various effective pi-conjugation lengths, and the energy transfer to the segment with a larger effective pi-conjugation length and a smaller pi-pi* transition energy required a longer transition time. A polarizing film was obtained by utilizing the perpendicular-type energy transfer.     

Protonation of the oxygen axial ligand in galactose oxidase model compounds as seen with high resolution X-ray emission experiments and FEFF simulations

X-ray Emission Spectroscopy (XES) crossover peaks were shown to be sensitive to the protonation state of solvent molecules in the Zn protein carbonic anhydrase and its model compounds. Here we extend such studies to galactose oxidase models i.e. Cu(ii) open d-shell systems, illustrating that XES combined with FEFF8 simulations reflect changes in the protonation state of the phenolate ligand for the copper center.