Dynamic mechanism of fatty acid transport across cellular membranes through FadL: molecular dynamics simulations

FadL is an important member of the family of fatty acid transport proteins within membranes. In this study, 11 conventional molecular dynamics (CMD) and 25 steered molecular dynamics (SMD) simulations were performed to investigate the dynamic mechanism of transport of long-chain fatty acids (LCFAs) across FadL. The CMD simulations addressed the intrinsically dynamic behavior of FadL. Both the CMD and SMD simulations revealed that a fatty acid molecule can move diffusively to a high-affinity site (HAS) from a low-affinity site (LAS). During this process, the swing motion of the L3 segment and the hydrophobic interaction between the fatty acid and FadL could play important roles. Furthermore, 22 of the SMD simulations revealed that fatty acids can pass through the gap between the hatch domain and the transmembrane domain (TMD) by different pathways. SMD simulations identified nine possible pathways for dodecanoic acid (DA) threading the barrel of FadL. The binding free energy profiles between DA and FadL along the MD trajectories indicate that all of the possible pathways are energetically favorable for the transport of fatty acids; however, one pathway (path VI) might be the most probable pathway for DA transport. The reasonability and reliability of this study were further demonstrated by correlating the MD simulation results with the available mutagenesis results. On the basis of the simulations, a mechanism for the full-length transport process of DA from the extracellular side to the periplasmic space mediated by FadL is proposed.     

氯化钠对朊病毒结构稳定性影响的分子动力学研究

朊病毒疾病是由正常构象的PrPC转化为致病构象的PrPSc引起的一类可传染的蛋白质构象病.采用分子动力学模拟的方法研究了0～500mmol/L的NaCl溶液体系对人朊病毒构象影响并深入探讨了其分子机制.研究发现NaCl可以降低朊病毒的结构稳定性,并引起其α-螺旋含量的急剧降低.进一步的研究表明高浓度NaCl溶液体系能够显著破坏朊病毒螺旋1内部的重要盐桥Asp144-Arg148和Asp147-Arg151,同时明显降低其主要氢键Arg151 N:Asp147 O,Tyr150 N:Glu146 O,Tyr149 N:Tyr145 O和Arg148 N:Asp144 O的稳定性,并诱导朊病毒的疏水核心发生明显扩张,促使朊病毒整体稳定性的下降,这些可能是NaCl促进朊病毒构象转换的重要原因.     

A quantum mechanics/molecular mechanics study of the protein-ligand interaction for inhibitors of HIV-1 integrase

Human immunodeficiency virus type-1 integrase (HIV-1 IN) is an essential enzyme for effective viral replication. Diketo acids such as L-731,988 and S-1360 are potent and selective inhibitors of HIV-1 IN. In this study, we used molecular dynamics simulations, within the hybrid quantum mechanics/molecular mechanics (QM/MM) approach, to determine the protein-ligand interaction energy between HIV-1 IN and L-731,988 and 10 of its derivatives and analogues. This hybrid methodology has the advantage that it includes quantum effects such as ligand polarisation upon binding, which can be very important when highly polarisable groups are embedded in anisotropic environments, as for example in metal-containing active sites. Furthermore, an energy decomposition analysis was performed to determine the contributions of individual residues to the enzyme-inhibitor interactions on averaged structures obtained from rather extensive conformational sampling. Analysis of the results reveals first that there is a correlation between protein-ligand interaction energy and experimental strand transfer into human chromosomes and secondly that the Asn-155, Lys-156 and Lys-159 residues and the Mg(2+) ion are crucial to anti-HIV IN activity. These results may explain the available experimental data.     

Car–Parrinello molecular dynamics on the SN2 reaction Cl+CH3Br in water

The results of the analysis of ab initio molecular dynamics simulations of the transition state of a bimolecular nuclear substitution reaction in water are discussed and a lower limit for the life-time of the hydrogen bonded interactions between the halogens and the surrounding water molecules is estimated. A new function is introduced to compute a lower estimate of the H-bond life-time and compared with the results of the traditional one.     

EGFR和4-苯胺喹唑啉类抑制剂之间相互作用模式的研究

采用分子动力学和MM/PBSA相结合的方法预测了表皮生长因子受体和4-苯胺喹 啉类抑制剂的相互作用模式。在分子动力学采样的基础上,采用MM/PBSA的方法分 别预测了四种可能结合模式下表皮生长因子受体和4-苯胺喹唑啉类抑制剂间的结合 自由能。在MM/PBSA计算中,受体和抑制剂之间的非键相互作用能采用分子力学 (MM)的方法得到;溶剂效应中极性部分对自由能的贡献通过解Possion- Boltzmanne (PB)方程的方法得到;溶液效应中非极性部分对自由能的贡献则通过 分子表面积计算(SA)的方法得到。计算表明,在四种结合模式下,表皮生长因子受 体和4-苯胺喹唑啉类抑制剂之间的结合自由能有较大的差别。在最佳的相互作用模 式中,抑制剂的苯胺部分位于活性口袋的底部,能够与受体残基的非极性侧链产生 很强的范德华和疏水相互作用。抑制剂喹唑啉环上的N(1)原子能够和Met-769上的 NH形成稳定的氢键,而抑制剂上的N(3)原子则和周围的一个水分子形成氢键。同时 ,抑制剂双环上的取代基团也能和活性口袋外部的部分残基形成一定的范德华和疏 水相互作用。最佳结合模式能够很好地解释已有抑制剂结构和活性间的关系。     

Absolute free energies of binding of peptide analogs to the HIV-1 protease from molecular dynamics simulations

The constants of binding of five peptide analogs to the active site of the HIV-1 aspartic-protease are calculated based on a novel sampling scheme that is efficient and does not introduce any approximations in addition to the energy function used to describe the system. The results agree with experiments. The squared correlation coefficient of the calculated vs. the measured values is 0.79. The sampling scheme consists of a series of molecular dynamics integrations with biases. The biases are selected based on an estimate of the probability density function of the system in a way to explore the conformational space and to reduce the statistical error in the calculated binding constants. The molecular dynamics integrations are done with a vacuum potential using a short cutoff scheme. To estimate the probability density of the simulated system, the results of the molecular dynamics integrations are combined using an extension of the weighted histogram analysis method (C. Bartels, Chem. Phys. Letters 331 (2000) 446-454). The probability density of the solvated ligand-protein system is obtained by applying a correction for the use of the short cutoffs in the simulations and by taking into account solvation with an electrostatic term and a hydrophobic term. The electrostatic part of the solvation is determined by finite difference Poisson-Boltzmann calculations; the hydrophobic part of the solvation is set proportional to the solvent accessible surface area. Setting the hydrophobic surface tension parameter equal to 8 mol(-1) K(-1) A(-2), absolute binding constants are in the muM to nM range. This is in agreement with experiments. The standard errors determined from eight repeated binding constant determinations are a factor of 14 to 411. A single determination of a binding constant is done with 499700 steps of molecular dynamics integration and 4500 finite difference Poisson-Boltzmann calculations. The simulations can be analyzed with respect to conformational changes of the active site of the HIV-1 protease or the ligands upon binding and provide information that complements experiments and can be used in the drug development process.     

Computer simulations of trypanosomal nucleoside hydrolase: determination of the protonation state of the bound transition-state analogue

Inosine-uridine nucleoside hydrolase (IU-NH) catalyzes the hydrolysis of nucleosides into base and ribose moieties via a ribooxocarbenium ion transition state, which has been characterized using kinetic isotope effects. Protozoan parasites lack de novo purine and pyrimidine biosynthesis and depend on the purine salvage from the host. Vern Schramm and co-workers characterized p-aminophenyliminoribitol (pAPIR) to be a potent inhibitor of IU-NH from Crithidia fasciculata with K(d) of 30 nM. The cyclic amine function of the iminoribitol ring can be either protonated (pAPIRH(+)) or unprotonated (pAPIR). pAPIRH(+) resembles the charge and geometry of the ribooxocarbenium ion transition state and can be looked upon as a transition-state analogue inhibitor; however, it is known that the pAPIR species is initially bound to the enzyme. We have characterized the pAPIRH(+) species as resident of the active site using ab initio calculations and molecular dynamics simulations. This is a novel use of molecular dynamics to investigate the protonation state of the bound ligand to the active site. Nanosecond molecular dynamics simulations reveal a short hydrogen-bonding network between pAPIRH(+)-O2'-Asp14-His241 triad, which is not seen in the crystal structure. Other features discussed are: hydrogen bonding between pAPIRH(+) and Asn168, unusual geometry of the iminoribitol ring, and hydrophobic interactions.     

纯水溶液中的氢键交换反应路径

利用分子动力学模拟方法对纯水溶液的氢键转化动力学性质进行了深入的微观探讨,溶液中非氢键构型为寿命较短(0.1～0.2 ps)的过渡态构型,我们发现氢键交换通过两种过渡构型完成,氢键角度扭曲激发后与氢键第一壳层水分子沿路径1交换,氢键径向拉伸激发后与氢键第二壳层水分子沿路径2交换,过渡态路径的选择具有温度依赖性.氢键转化需在旧氢键受体氢键过量和新氢键受体氢键不足,同时满足交换反应空间结构要求下才能完成.氢键交换反应对水分子平动和转动行为起着决定作用.     

Drug release by pH‐responsive molecular tweezers: Atomistic details from molecular modeling

pH‐responsive molecular tweezers have been proposed as an approach for targeting drug‐delivery to tumors, which tend to have a lower pH than normal cells. We performed a computational study of a pH‐responsive molecular tweezer using ab initio quantum chemistry in the gas‐phase and molecular dynamics (MD) simulations in solution. The binding free energy in solution was calculated using steered MD. We observe, in atomistic detail, the pH‐induced conformational switch of the tweezer and the resulting release of the drug molecule. Even when the tweezer opens, the drug molecule remains near a hydrophobic arm of the molecular tweezer. Drug release cannot occur, it seems, unless the tweezer is in a hydrophobic environment with low pH. © 2014 Wiley Periodicals, Inc.     

Surface protonation at the rutile (110) interface: explicit incorporation of solvation structure within the refined MUSIC model framework

The detailed solvation structure at the (110) surface of rutile (alpha-TiO2) in contact with bulk liquid water has been obtained primarily from experimentally verified classical molecular dynamics (CMD) simulations of the ab initio-optimized surface in contact with SPC/E water. The results are used to explicitly quantify H-bonding interactions, which are then used within the refined MUSIC model framework to predict surface oxygen protonation constants. Quantum mechanical molecular dynamics (QMD) simulations in the presence of freely dissociable water molecules produced H-bond distributions around deprotonated surface oxygens very similar to those obtained by CMD with nondissociable SPC/E water, thereby confirming that the less computationally intensive CMD simulations provide accurate H-bond information. Utilizing this H-bond information within the refined MUSIC model, along with manually adjusted Ti-O surface bond lengths that are nonetheless within 0.05 A of those obtained from static density functional theory (DFT) calculations and measured in X-ray reflectivity experiments (as well as bulk crystal values), give surface protonation constants that result in a calculated zero net proton charge pH value (pHznpc) at 25 degrees C that agrees quantitatively with the experimentally determined value (5.4+/-0.2) for a specific rutile powder dominated by the (110) crystal face. Moreover, the predicted pHznpc values agree to within 0.1 pH unit with those measured at all temperatures between 10 and 250 degrees C. A slightly smaller manual adjustment of the DFT-derived Ti-O surface bond lengths was sufficient to bring the predicted pHznpcvalue of the rutile (110) surface at 25 degrees C into quantitative agreement with the experimental value (4.8+/-0.3) obtained from a polished and annealed rutile (110) single crystal surface in contact with dilute sodium nitrate solutions using second harmonic generation (SHG) intensity measurements as a function of ionic strength. Additionally, the H-bond interactions between protolyzable surface oxygen groups and water were found to be stronger than those between bulk water molecules at all temperatures investigated in our CMD simulations (25, 150 and 250 degrees C). Comparison with the protonation scheme previously determined for the (110) surface of isostructural cassiterite (alpha-SnO2) reveals that the greater extent of H-bonding on the latter surface, and in particular between water and the terminal hydroxyl group (Sn-OH) results in the predicted protonation constant for that group being lower than for the bridged oxygen (Sn-O-Sn), while the reverse is true for the rutile (110) surface. These results demonstrate the importance of H-bond structure in dictating surface protonation behavior, and that explicit use of this solvation structure within the refined MUSIC model framework results in predicted surface protonation constants that are also consistent with a variety of other experimental and computational data.     

Structural dissimilarity sampling with dynamically self‐guiding selection

Structural dissimilarity sampling (SDS) has been proposed as an enhanced conformational sampling method for reproducing the structural transitions of a given protein. SDS consists of cycles of two steps: (1) Selections of initial structures with structural dissimilarities by referring to a measure. (2) Conformational resampling by restarting short‐time molecular dynamics (MD) simulations from the initial structures. In the present study, an efficient measure is proposed as a dynamically self‐guiding selection to accelerate the structural transitions from a reactant state to a product state as an extension to the original SDS. In the extended SDS, the inner product (IP ) between the reactant and the snapshots generated by short‐time MD simulations are evaluated and ranked according to the IP s at every cycle. Then, the snapshots with low IP s are selected as initial structures for the short‐time MD simulations. This scheme enables one to choose dissimilar and distant initial structures from the reactant, and thus the initial structures dynamically head towards the product, promoting structural transitions from the reactant. To confirm the conformational sampling efficiency, the extended SDS was applied to maltodextrin binding protein (MBP), and we successfully reproduced the structural transition from the open to closed states with submicrosecond‐order simulation times. However, a conventional long‐time MD simulation failed to reproduce the same structural transition. We also compared the performance with that obtained by the ordinary SDS and other sampling techniques that have been developed by us to characterize the possible utility of the extended SDS for actual applications. © 2017 Wiley Periodicals, Inc.     

Study of H-bond characteristics in sub- and supercritical methanol

Classical molecular dynamics simulations of various methanol phase lines near the saturation curve and the critical point have been performed to study the changes in H-bonded clusters structure at transition of methanol to supercritical state. Analysis of H-bonds statistics with combined distance-energy H-bond criterion showed that the correlations between topological characteristics of H-bonds and the mole fraction of H-bonded molecules have unique functional representation despite the phase path applied. In the present study, an attempt has been also made to evaluate the degree of hydrogen bonding by combining the DFT computations on classical MD configurations with the natural bond orbital analysis of the waves functions obtained.     

多肽抑制剂抑制淀粉质多肽42构象转换的分子动力学模拟和结合自由能计算

应用分子动力学模拟和结合自由能计算方法研究了多肽抑制剂KLVFF、VVIA和LPFFD抑制淀粉质多肽42 (Aβ42)构象转换的分子机理. 结果表明, 三种多肽抑制剂均能够有效抑制Aβ42的二级结构由α-螺旋向β-折叠的构象转换. 另外, 多肽抑制剂降低了Aβ42分子内的疏水相互作用, 减少了多肽分子内远距离的接触, 有效抑制了Aβ42的疏水塌缩, 从而起到稳定其初始构象的作用. 这些抑制剂与Aβ42之间的疏水和静电相互作用(包括氢键)均有利于它们抑制Aβ42的构象转换. 此外, 抑制剂中的带电氨基酸残基可以增强其和Aβ42之间的静电相互作用(包括氢键), 并降低抑制剂之间的聚集, 从而大大增强对Aβ42构象转换的抑制能力. 但脯氨酸的引入会破坏多肽的线性结构, 从而大大降低其与Aβ42 之间的作用力. 上述分子模拟的结果揭示了多肽抑制剂KLVFF、VVIA和LPFFD抑制Aβ42构象转换的分子机理, 对于进一步合理设计Aβ的高效短肽抑制剂具有非常重要的理论指导意义.     

Side-Chain Interactions in d/l Peptide Nanotubes: Studies by Crystallography,NMR Spectroscopy and Molecular Dynamics

Our understanding of the factors affecting the stability of cyclic d/l peptide (CP) nanotubes remains underdeveloped. In this work, we investigate the impact of side chain alignment, hydrophobicity and charge on CP nanotube stability through X-ray crystallography, NMR spectroscopy and molecular dynamics (MD) simulations. We characterise the distinct CP-CP alignments that can form and identify stable and unstable dimers by MD simulation. We measure H-bond half-lives of synthesised CPs by ¹H−D exchange experiments and find good correlation with predicted CP-CP stabilities. We find that hydrophobic amino acids improve CP dimer stability but experimentally reduce solubility. Charged amino acids either increase or decrease CP dimer stability depending on the relative orientation and composition of charged groups. X-ray crystal structures are solved for two CPs, revealing non-tubular folded conformations. Ultimately, this work will assist the educated design of stable tubular structures for potential applications in biomedicine.     

Understanding interactions between poly(styrene-co-sodium styrene sulfonate) and single-walled carbon nanotubes

Understanding formation mechanisms of hybrids of carbon nanotubes (CNTs) wrapped by polymers and their interactions is critical in modifying solubility of CNTs in aqueous solution and developing new nanotube-based polymer materials. In the present work, we investigate the structural details of poly(styrene-co-sodium styrene sulfonate) (PSS) wrapping around the CNT and the interactions between the PSS chain and the CNT using molecular dynamics (MD) simulations. The fraction of sulfonated groups significantly influences the wrapping conformations of the PSS chain. Due to limited time scale in the MD simulations, two different initial conformations of the chains are introduced to explore the effect of the initial state on the wrapping behavior. When the chains initially wrap around the CNT in a perfect helix manner, more compact pseudo-helical conformations are obtained. For initial straight line arrangement of the chain monomers, the chains adopt looser wrapping conformations. The free-energy analysis and binding interaction of the PSS chain on the CNT surface take a glance on the relationship between the conformational transition of the chain and the energy evolution.     

Gating Mechanism of the Voltage-Gated Proton Channel Studied by Molecular Dynamics Simulations

The voltage-gated proton channel Hv1 has important roles in proton extrusion, pH homeostasis, sperm motility, and cancer progression. The Hv1 channel has also been found to be highly expressed in cell lines and tissue samples from patients with breast cancer. A high-resolution closed-state structure has been reported for the mouse Hv1 chimera channel (mHv1cc), solved by X-ray crystallography, but the open-state structure of Hv1 has not been solved. Since Hv1 is a promising drug target, various groups have proposed open conformations by molecular modeling and simulation studies. However, the gating mechanism and the open-state conformation under the membrane potential are still debate. Here, we present a molecular dynamics study considering membrane potential and pH conditions. The closed-state structure of mHv1cc was used to run molecular dynamics (MD) simulations with respect to electric field and pH conditions in order to investigate the mechanism of proton transfer. We observed a continuous hydrogen bond chain of water molecules called a water-wire to be formed through the channel pore in the channel opening, triggered by downward displacement of the S2 helix and upward movement of the S4 helix relative to other helices. Due to the movement of the S2 and S4 helices, the internal salt bridge network was rearranged, and the hydrophobic gating layers were destroyed. In line with previous experimental and simulation observations, our simulation results led us to propose a new gating mechanism for the Hv1 proton channel, and may provide valuable information for novel drug discovery.     

First-principles simulation of the photoreaction of a capped azobenzene: the rotational pathway is feasible.

We present first-principles molecular dynamics simulations of azobenzene and a sterically hindered derivative in the first excited state. The restricted open-shell Kohn-Sham (ROKS) approach is employed to describe the motion in the lowest excited state. The rotational pathway is observed in the molecular dynamics simulations for both azobenzene and its azacrown ether capped derivative.