Alamethicin channels in a membrane: molecular dynamics simulations

Alamethicin (Alm) is a 20 residue peptide which forms a kinked alpha-helix in membrane and membrane-mimetic environments. Ion channels formed by intramembraneous aggregates of Alm are thought to be formed by bundles of approximately parallel Alm helices surrounding a central bilayer pore. Different channel conductance levels correspond to different numbers of helices per bundle, ranging from N = 5 to N > 8. Calculation of the predicted pKA values of the ring of Glu18 sidechains at the C-terminal mouth of the pore suggests that at neutral pH most or all of these sidechains will remain protonated. Nanosecond molecular dynamics (MD) simulations of N = 5, 6, 7 and 8 bundles of Alm helices in a POPC bilayer have been run, corresponding to a total simulation time of 4 ns. These simulations explore the stability and conformational dynamics of these helix bundle channels when embedded in a full phospholipid bilayer in an aqueous environment. The structural and dynamic properties of water in these model channels are examined. As in earlier in vacuo simulations (J. Breed, R. Sankararamakrishnan, I. D. Kerr and M. S. P. Sansom, Biophys. J., 1996, 70, 1643) the dipole moments of water molecules within the pores are aligned antiparallel to the helix dipoles. This helps to contribute to the stability of the helix bundles.     

固有无序超嗜热菌FlgM蛋白在两种不同温度下的构象变化特征

基于分子动力学模拟方法比较了超嗜热菌FlgM 蛋白在常温(293 K)和生理温度(358 K)下的结构特征.基于GROMACS软件包, 采用OPLS-AA分子力场和TIP3P水模型, 对超嗜热菌FlgM 蛋白在293 和358 K进行了2组独立的长时间分子动力学模拟, 每组体系模拟时间为1500 ns. 主要分析了两种温度下超嗜热菌FlgM蛋白的二级结构特征、整体构象变化及半无序化区域和结构化区域的构象特征. 研究结果表明: 在常温下, N端具有一定程度的螺旋成分, 但在生理温度下, 超嗜热菌FlgM 蛋白的结构变得松散, 螺旋结构减少, 构象稳定性减弱, H1 螺旋散开, FlgM 蛋白构象灵活性增强, 不稳定程度增加. 这些不同温度的结构变化说明: 半无序化区域(N端)在非结合状态下有一定的螺旋结构, 但该段螺旋的稳定性随温度升高而降低. 超嗜热菌FlgM蛋白会通过增加结构的无序程度使结构更加灵活, 以适应高温, 从而使该类固有无序蛋白更好地行使其功能, 如提高同其他成分的结合速率等.     

Surface tension effect on transmembrane channel stability in a model membrane

The effect of surface tension on the lipid bilayer membrane is a question that has drawn considerable research effort. This interest has been driven both by the desire to determine the surface tension effects on the lipid bilayer and from the suggestion that adding finite surface tension to a small membrane system may provide more realistic lipid properties in molecular dynamics simulations. Here, the effect of surface tension on a palmitololelylphosphatidylcholine (POPC) bilayer membrane containing a four-helix transmembrane alamethicin peptide bundle is investigated. Simulations of 10 ns were undertaken for two different ensembles, NPT and NP(z)gammaT with a surface tension, gamma, of 20 mN m(-1) per interface, which is near the pore-forming region. The significance of differences between the tension-free and surface tension simulations was determined using nonparametric statistical analysis on replicate simulations with different initial conditions. The results suggest that, when the membrane is under surface tension, the peptide helical structure is perturbed from that in the tension-free state but that the bundle conformation is more stable than that in the tension-free state, with hydrogen bonding playing an important stabilizing role. Surface tension counteracts the influence of the transmembrane helix bundle on nearby lipid order, making the lipid order more uniform throughout the membrane in the tension state. Conversely, the lipid mobility was less uniform in the tension state, with lipids far from the bundle being significantly more mobile than those near the bundle. One general implication of the results is that surface tension can affect the membrane nonuniformly, in that the properties of lipids near the peptide are different from those further away.     

A generalized Born formalism for heterogeneous dielectric environments: application to the implicit modeling of biological membranes

Reliable computer simulations of complex biological environments such as integral membrane proteins with explicit water and lipid molecules remain a challenging task. We propose a modification of the standard generalized Born theory of homogeneous solvent for modeling the heterogeneous dielectric environments such as lipid/water interfaces. Our model allows the representation of biological membranes in the form of multiple layered dielectric regions with dielectric constants that are different from the solute cavity. The proposed new formalism is shown to predict the electrostatic component of solvation free energy with a relative error of 0.17% compared to exact finite-difference solutions of the Poisson equation for a transmembrane helix test system. Molecular dynamics simulations of melittin and bacteriorhodopsin are carried out and performed over 10 ns and 7 ns of simulation time, respectively. The center of melittin along the membrane normal in these stable simulations is in excellent agreement with the relevant experimental data. Simulations of bacteriorhodopsin started from the experimental structure remained stable and in close agreement with experiment. We also examined the free energy profiles of water and amino acid side chain analogs upon membrane insertion. The results with our implicit membrane model agree well with the experimental transfer free energy data from cyclohexane to water as well as explicit solvent simulations of water and selected side chain analogs.     

What is the time scale for α-helix nucleation?

Helix formation is an elementary process in protein folding, influencing both the rate and mechanism of the global folding reaction. Yet, because helix formation is less cooperative than protein folding, the kinetics are often multiexponential, and the observed relaxation times are not straightforwardly related to the microscopic rates for helix nucleation and elongation. Recent ultrafast spectroscopic measurements on the peptide Ac-WAAAH(+)-NH(2) were best fit by two relaxation modes on the ～0.1-1 ns time scale, (1) apparently much faster than had previously been experimentally inferred for helix nucleation. Here, we use replica-exchange molecular dynamics simulations with an optimized all-atom protein force field (Amber ff03w) and an accurate water model (TIP4P/2005) to study the kinetics of helix formation in this peptide. We calculate temperature-dependent microscopic rate coefficients from the simulations by treating the dynamics between helical states as a Markov process using a recently developed formalism. The fluorescence relaxation curves obtained from simulated temperature jumps are in excellent agreement with the experimentally determined results. We find that the kinetics are multiphasic but can be approximated well by a double-exponential function. The major processes contributing to the relaxation are the shrinking of helical states at the C-terminal end and a faster re-equilibration among coil states. Despite the fast observed relaxation, the helix nucleation time is estimated from our model to be 20-70 ns at 300 K, with a dependence on temperature well described by Arrhenius kinetics.     

Molecular dynamics of nicotinic acetylcholine receptor correlating biological functions

The nicotinic acetylcholine receptor (nAChR) that mediates fast intercellular communication in response to neurotransmitters is a paradigm of ligand-gated ion channels. Molecular dynamics (MD) simulations are valuable in understanding membrane protein function at atomic level, providing useful clues for further experimental/theoretical studies. In this brief review, recent progress in MD simulations of the nAChR has been illustrated, mainly focusing on the latest simulation of the whole transmembrane domain of the receptor. On the basis of MD simulations, asymmetrical and asynchronous motions of five subunits were observed both in the ligand binding and transmembrane domains; a closed-to-open conformational shift of the gate was captured in different simulation systems; the contributions from the lipid molecules and other transmembrane segments rather than M2 to the gate switch as well as the conformational change of the whole channel were assessed; the dynamic behavior and related physical/chemical properties of the water molecules and cations within the ion channel were examined; and an experimentally comparable single-channel conductance and ion selectivity were obtained.     

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.     

Insights into the dynamics of HIV-1 protease: a kinetic network model constructed from atomistic simulations

The conformational dynamics in the flaps of HIV-1 protease plays a crucial role in the mechanism of substrate binding. We develop a kinetic network model, constructed from detailed atomistic simulations, to determine the kinetic mechanisms of the conformational transitions in HIV-1 PR. To overcome the time scale limitation of conventional molecular dynamics (MD) simulations, our method combines replica exchange MD with transition path theory (TPT) to study the diversity and temperature dependence of the pathways connecting functionally important states of the protease. At low temperatures the large-scale flap opening is dominated by a small number of paths; at elevated temperatures the transition occurs through many structurally heterogeneous routes. The expanded conformation in the crystal structure 1TW7 is found to closely mimic a key intermediate in the flap-opening pathways at low temperature. We investigated the different transition mechanisms between the semi-open and closed forms. The calculated relaxation times reveal fast semi-open ? closed transitions, and infrequently the flaps fully open. The ligand binding rate predicted from this kinetic model increases by 38-fold from 285 to 309 K, which is in general agreement with experiments. To our knowledge, this is the first application of a network model constructed from atomistic simulations together with TPT to analyze conformational changes between different functional states of a natively folded protein.     

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.     

Surfaces of alcohol-water mixtures studied by sum-frequency generation vibrational spectroscopy

Sum-frequency generation vibrational spectroscopy was used to investigate the surface molecular structure of binary mixtures of water and alcohol (methanol, ethanol, and propanol) at the air/liquid interface. In this study, it is shown that the sum-frequency signal from the alcohol molecules in the CH-stretch vibration region is always larger for mixtures than that from pure alcohol. For example, the sum-frequency signal from a propanol mixture surface at a 0.1 bulk mole fraction was approximately 3 times larger than that from a pure propanol surface. However, the ratio between the sum-frequency signals taken at different polarization combinations was found to be constant within experimental errors as the bulk alcohol concentration was changed. This suggested that the orientation of surface alcohol molecules does not vary appreciably with the change of concentration and that the origin of the signal enhancement is mainly due to the increase in the surface number density of alcohol molecules contributing to the sum-frequency signal for the alcohol/water mixture as compared to the pure alcohol surface.     

The Aminotriazole Antagonist Cmpd‐1 Stabilises a Distinct Inactive State of the Adenosine 2A Receptor

The widely expressed G‐protein coupled receptors (GPCRs) are versatile signal transducer proteins that are attractive drug targets but structurally challenging to study. GPCRs undergo a number of conformational rearrangements when transitioning from the inactive to the active state but have so far been believed to adopt a fairly conserved inactive conformation. Using ¹⁹F NMR spectroscopy and advanced molecular dynamics simulations we describe a novel inactive state of the adenosine 2A receptor which is stabilised by the aminotriazole antagonist Cmpd‐1. We demonstrate that the ligand stabilises a unique conformation of helix V and present data on the putative binding mode of the compound involving contacts to the transmembrane bundle as well as the extracellular loop 2.     

Mpro-C蛋白三维结构域交换的机理: 来自分子模拟的线索

SARS冠状病毒主蛋白酶(M^pro)在病毒的蛋白酶切过程中发挥着重要作用. M^pro的晶体结构显示它存在两种形式的二聚体: 一种是发生三维结构域交换的形式, 另一种是非交换的形式. M^pro的C端结构域(M^pro-C)单独表达时也能形成与全长M^pro类似的三维结构域交换二聚体. 三维结构域交换通常发生在蛋白质的表面, 但M^pro-C 的结构域交换却发生在疏水核心. 在本文中, 我们利用分子动力学模拟及三维结构域交换预测算法研究了M^pro-C 中被高度埋藏的核心螺旋片段发生交换的机理. 我们发现基于结构与基于序列的已有算法都不能正确预言出M^pro-C和M^pro中发生结构域交换的铰链区位置. 分子模拟结果表明M^pro-C中的交换片段在天然态下埋藏得很好, 但在变性单体中则会被释放并暴露在外面. 因此, 在完全或部分解折叠状态下交换片段的打开有助于促进单体间的相互作用及结构域交换二聚体的形成.     

Characterization of the dynamics of an essential helix in the U1A protein by time-resolved fluorescence measurements

The RNA recognition motif (RRM), one of the most common RNA-binding domains, recognizes single-stranded RNA. A C-terminal helix that undergoes conformational changes upon binding is often an important contributor to RNA recognition. The N-terminal RRM of the U1A protein contains a C-terminal helix (helix C) that interacts with the RNA-binding surface of a beta-sheet in the free protein (closed conformation), but is directed away from this beta-sheet in the complex with RNA (open conformation). The dynamics of helix C in the free protein have been proposed to contribute to binding affinity and specificity. We report here a direct investigation of the dynamics of helix C in the free U1A protein on the nanosecond time scale using time-resolved fluorescence anisotropy. The results indicate that helix C is dynamic on a 2-3 ns time scale within a 20 degrees range of motion. Steady-state fluorescence experiments and molecular dynamics simulations suggest that the dynamical motion of helix C occurs within the closed conformation. Mutation of a residue on the beta-sheet that contacts helix C in the closed conformation dramatically destabilizes the complex (Phe56Ala) and alters the steady-state fluorescence, but not the time-resolved fluorescence anisotropy, of a Trp in helix C. Mutation of Asp90 in the hinge region between helix C and the remainder of the protein to Ala or Gly subtly alters the dynamics of the U1A protein and destabilizes the complex. Together these results show that helix C maintains a dynamic closed conformation that is stable to these targeted protein modifications and does not equilibrate with the open conformation on the nanosecond time scale.     

Molecular dynamics simulations and thermodynamics analysis of DNA-drug complexes. Minor groove binding between 4',6-diamidino-2-phenylindole and DNA duplexes in solution

An extended set of nanosecond-scale molecular dynamics simulations of DNA duplex sequences in explicit solvent interacting with the minor groove binding drug 4',6-diamidino-2-phenylindole (DAPI) are investigated for four different and sequence specific binding modes. Force fields for DAPI have been parametrized to properly reflect its internal nonplanarity. Sequences investigated include the binding modes observed experimentally, that is, AATT in d(CGCGAATTCGCG)(2) and ATTG in d(GGCCAATTGG)(2) and alternative shifted binding modes ATTC and AATT, respectively. In each case, stable MD simulations are obtained, well reproducing specific hydration patterns seen in the experiments. In contrast to the 2.4 A d(CGCGAATTCGCG)(2) crystal structure, the DAPI is nonplanar, consistent with its gas-phase geometry and the higher resolution crystal structure. The simulations also suggest that the DAPI molecule is able to adopt different conformational substates accompanied by specific hydration patterns that include long-residing waters. The MM_PBSA technology for estimating relative free energies was utilized. The most consistent free energy results were obtained with an approach that uses a single trajectory of the DNA-DAPI complex to estimate all free energy terms. It is demonstrated that explicit inclusion of a subset of bound water molecules shifts the calculated relative binding free energies in favor of both crystallographically observed binding modes, underlining the importance of structured hydration.     

W2466.48 Opens a Gate for a Continuous Intrinsic Water Pathway during Activation of the Adenosine A2A Receptor

The question how G‐protein‐coupled receptors transduce an extracellular signal by a sequence of transmembrane conformational transitions into an intracellular response remains to be solved at molecular detail. Herein, we use molecular dynamics simulations to reveal distinct conformational transitions of the adenosine A_2A receptor, and we found that the conserved W246^6.48 residue in transmembrane helix TM6 performs a key rotamer toggle switch. Agonist binding induces the sidechain of W246^6.48 to fluctuate between two distinct conformations enabling the diffusion of water molecules from the bulk into the center of the receptor. After passing the W246^6.48 gate, the internal water molecules induce another conserved residue, Y288^7.53, to switch to a distinct rotamer conformation establishing a continuous transmembrane water pathway. Further, structural changes of TM6 and TM7 induce local structural changes of the adjacent lipid bilayer.     

应用ABEEMσπ／MM模型探讨小α-螺旋在显性水中折叠／展开的可逆过程

应用ABEEMσπ／MM模型进行分子动力学模拟,研究了显性水溶液中小α-螺旋（短肽Ala5）折叠／展开的可逆过程．动力学分析显示,300K下α-螺旋可以保存2ns的时间,该结果支持Margulis等人的结论．每个结构与α-螺旋结构骨架重原子的均方根偏差的时间轨迹指出,“300K下螺旋成核现象在0．1ns内快速发生”的结论是不恰当的．通过对300、400和500K温度下的研究,首次定量地给出各温度下螺旋保存的时间分别为2ns、1～1．5ns和0．8ns,并且增加温度并不改变折叠／展开的方式,只是改变折叠／展开的速率．本文对“转化态集合”结构的分析表明,从螺旋到卷曲的转换,主要通过螺旋端的氢键断裂发生（92％）、尤其是C端的氢键断裂发生（50％）．氢键的破坏和形成在0．1ns的时间内完成．     

Molecular dynamics simulation of the A-DNA to B-DNA transition in aqueous RbCl solution

Unrestrained molecular dynamics (MD) simulations have been carried out to characterize the stability of DNA conformations and the dynamics of A-DNA→B-DNA conformational transitions in aqueous RbCl solutions. The PARM99 force field in the AMBER8 package was used to investigate the effect of RbCl concentration on the dynamics of the A→B conformational transition in the DNA duplex d(CGCGAATTCGCG)2 . Canonical Aand B-form DNA were assumed for the initial conformation and the final conformation had a length per complete turn that matched the canonical B-DNA. The DNA structure was monitored for 3.0 ns and the distances between the C5′ atoms were obtained from the simulations. It was found that all of the double stranded DNA strands of A-DNA converged to the structure of B-form DNA within 1.0 ns during the unrestrained MD simulations. In addition, increasing the RbCl concentration in aqueous solution hindered the A→B conformational transition and the transition in aqueous RbCl solution was faster than that in aqueous NaCl solution for the same electrolyte strength. The effects of the types and concentrations of counterions on the dynamics of the A→B conformational transition can be understood in terms of the variation in water activity and the number of accumulated counterions in the major grooves of A-DNA. The rubidium ion distributions around both fixed A-DNA and B-DNA were obtained using the restrained MD simulations to help explain the effect of RbCl concentration on the dynamics of the A→B conformational transition.     

Evidence of complete hydrophobic coating of bombesin by trifluoroethanol in aqueous solution: an NMR spectroscopic and molecular dynamics study

Bombesin is a tetradecapeptide that possesses a random coil structure in pure water. In the presence of 30 % (v/v) 2,2,2-trifluoroethanol (TFE), it adopts a partial helical conformation involving the C-terminal amino acids 6-14. This conformational change, known as the TFE effect, is studied here in terms of the solvation state of the peptide at different TFE concentrations by means of intermolecular homo- and heteronuclear NOE measurements. When an aqueous solution of bombesin is titrated with TFE, a continual decrease in the water/peptide interactions and a concomitant increase in the TFE/peptide interactions is observed, and at 30 % (v/v) TFE no homonuclear NOEs between water and the peptide can be detected. The conformational transition of the bombesin molecule is thus accompanied by a complete surface covering with TFE. A parallel molecular dynamics (MD) study of the peptide in aqueous solution with the single-point charge (SPC) water model and in a 30 % (v/v) TFE/water mixture with a recently developed TFE model has also been performed. The 10 ns simulations were in agreement with the experimental data. The calculations indicate stabilisation of the alpha-helix in the H(2)O/TFE mixture, in contrast to the situation in pure water, and clustering of the TFE molecules around the peptide.     

Folding processes of the B domain of protein A to the native state observed in all-atom ab initio folding simulations

Reaching the native states of small proteins, a necessary step towards a comprehensive understanding of the folding mechanisms, has remained a tremendous challenge to ab initio protein folding simulations despite the extensive effort. In this work, the folding process of the B domain of protein A (BdpA) has been simulated by both conventional and replica exchange molecular dynamics using AMBER FF03 all-atom force field. Started from an extended chain, a total of 40 conventional (each to 1.0 micros) and two sets of replica exchange (each to 200.0 ns per replica) molecular dynamics simulations were performed with different generalized-Born solvation models and temperature control schemes. The improvements in both the force field and solvent model allowed successful simulations of the folding process to the native state as demonstrated by the 0.80 A C(alpha) root mean square deviation (RMSD) of the best folded structure. The most populated conformation was the native folded structure with a high population. This was a significant improvement over the 2.8 A C(alpha) RMSD of the best nativelike structures from previous ab initio folding studies on BdpA. To the best of our knowledge, our results demonstrate, for the first time, that ab initio simulations can reach the native state of BdpA. Consistent with experimental observations, including Phi-value analyses, formation of helix II/III hairpin was a crucial step that provides a template upon which helix I could form and the folding process could complete. Early formation of helix III was observed which is consistent with the experimental results of higher residual helical content of isolated helix III among the three helices. The calculated temperature-dependent profile and the melting temperature were in close agreement with the experimental results. The simulations further revealed that phenylalanine 31 may play critical to achieve the correct packing of the three helices which is consistent with the experimental observation. In addition to the mechanistic studies, an ab initio structure prediction was also conducted based on both the physical energy and a statistical potential. Based on the lowest physical energy, the predicted structure was 2.0 A C(alpha) RMSD away from the experimentally determined structure.