Computational study of the transmembrane domain of the acetylcholine receptor

The nicotinic acetylcholine receptor (nAChR) is a ligand-gated ion channel protein whose transmembrane domain (TM-domain) is believed to be responsible for channel gating via a hydrophobic effect. In this work, we perform molecular dynamics and Brownian dynamics simulations to investigate the effect of transmembrane potential on the conformation and water occupancy of TM-domain, and the resulting ion permeation events. The results show that the behavior of the hydrophobic gate is voltage-dependent. Large hyperpolarized membrane potential can change the conformation of TM-domain and water occupancy in this region, which may enable ion conduction. An electrostatic gating mechanism is also proposed from our simulations, which seems to play a role in addition to the well-known hydrophobic effect.     

Determination of the conformation of d(GGAAATTTCC)2 in solution by use of 1H NMR and restrained molecular dynamics

The conformation of the putative bent DNA d(GGAAATTTCC)2 in solution was studied by use of 1H NMR and restrained molecular dynamics. Most of the resonances were assigned sequentially. A total of 182 interproton distance restraints were determined from two-dimensional nuclear Overhauser effect spectra with short mixing times. Torsion angle restraints for each sugar moiety were determined by qualitative analysis of a two-dimensional correlated spectrum. Restrained molecular dynamics was carried out with the interproton distances and torsion angles incorporated into the total energy function of the system in the form of effective potential terms. As initial conformations for restrained molecular dynamics, classical A-DNA and B-DNA were adopted. The root mean square deviation (rmsd) between these two conformations is 5.5 A. The conformations obtained by use of restrained molecular dynamics are very similar to each other, the rmsd being 0.8 A. On the other hand, the conformations obtained by use of molecular dynamics without experimental restraints or restrained energy minimization depended heavily on the initial conformations, and convergence to a similar conformation was not attained. The conformation obtained by use of restrained molecular dynamics exhibits a few remarkable features. The second G residue takes on the BII conformation [Fratini, A. V., Kopka, M. L., Drew, H. R., & Dickerson, R. E. (1982) J. Biol. Chem. 257, 14686-14707] rather than the standard BI conformation. There is discontinuity of the sugar puckering between the eighth T and ninth C. The minor groove of the oligo(dA) tract is rather compressed. As a result, d(GGAAATTTCC)2 is bent.     

Activity and conformational changes of alpha-chymotrypsin in reverse micelles studied by spin labeling.

alpha-Chymotrypsin (CT), spin-labeled at the active site by using an acylating label which constitutes a substrate for this protein, has been investigated in reverse micelles formed by AOT in isooctane. The electron spin resonance spectra provided information on conformation, dynamics and deacylation activity. The dynamics of the label bound to CT appears to be more hindered in reverse micelles than in aqueous solution, probably owing to the effect of the micellar environment on protein conformation. The deacylation rate in reverse micelles does not show the characteristic bell-shaped dependence on water content which is generally found for CT enzymatic activity.     

Theoretical studies of the structure and molecular dynamics of a peptide crystal

Energy minimizations and molecular dynamics simulations have been performed on the cyclic peptide cyclo-(Ala-Pro-D-Phe)2 in both the isolated and crystal states. The results of these calculations have been analyzed, both to investigate our ability to reproduce experimental data (structure and vibrational and NMR spectra) and to investigate the effects of environment on the energy, structure, and dynamics of peptides. Comparison of the minimized and time-averaged crystal systems with the experimental peptide structure shows that the calculations have closely reproduced the experimental structure. Molecular dynamics of the isolated molecule has led to a new conformation, which is approximately equal to 8.5 kcal/mol more stable than the conformation that exists in the crystal, the latter conformation being stabilized by intermolecular (packing) forces. This illustrates the considerable effect that environment can have on the conformation of peptides. The crystal environment has also been shown to significantly reduce the dynamic conformational fluctuations seen for the isolated molecule. The behavior of the peptide during the isolated simulation also supports the experimental NMR observation of a symmetric structure that differs from the asymmetric, instantaneous structures which characterize the molecule during the dynamics. Calculations of vibrational frequencies of the peptide in the crystal and isolated states show the expected shifts in bond-stretching frequencies due to intermolecular interactions. Finally, we have calculated NMR coupling constants from the dynamics simulation of the isolated peptide and have compared these with the experimental values. This has led to a possible reinterpretation of the experimental data.     

Theoretical investigation of the dynamics of the active site lid in Rhizomucor miehei lipase.

Interfacial activation of Rhizomucor miehei lipase is accompanied by a hinge-type motion of a single helix (residues 83-94) that acts as a lid over the active site. Activation of the enzyme involves the displacement of the lid to expose the active site, suggesting that the dynamics of the lid could be of mechanistic and kinetic importance. To investigate possible activation pathways and to elucidate the effect of a hydrophobic environment (as would be provided by a lipid membrane) on the lid opening, we have applied molecular dynamics and Brownian dynamics techniques. Our results indicate that the lipase activation is enhanced in a hydrophobic environment. In nonpolar low-dielectric surroundings, the lid opens in approximately 100 ns in the BD simulations. In polar high-dielectric (aqueous) surroundings, the lid does not always open up in simulations of up to 900 ns duration, but it does exhibit some gating motion, suggesting that the enzyme molecule may exist in a partially active form before the catalytic reaction. The activation is controlled by the charged residues ARG86 and ASP91. In the inactive conformation, ASP91 experiences repulsive forces and pushes the lid toward the open conformation. Upon activation ARG86 approaches ASP61, and in the active conformation, these residues form a salt bridge that stabilizes the open conformation.     

Conformational properties of oxytocin in dimethyl sulfoxide solution: NMR and restrained molecular dynamics studies.

The conformation of oxytocin, the neurohypophyseal nonapeptide hormone, in solution in deuterated dimethyl sulfoxide has been determined by 1H-nmr. The structural determination is based on the experimental data set of nuclear Overhauser effect restraints. Obtained after the restrained molecular dynamics simulation on an initial structure of extended conformation, five resultant structures satisfy the experimental restraints well. These structures resemble that of the crystal structure of deamino-oxytocin, an analogue of oxytocin, in terms of a close correlation observed both at two beta-turn regions of the 20-membered tocin ring and at the tripeptide tail end. Based on this comparison and analysis of restrained molecular dynamics trajectories, we found that, although the turns are stabilized by the formation of hydrogen bonds, the oxytocin molecule possesses a slight twist in DMSO solution relative to the orientation of deamino-oxytocin in the crystalline state. Analyses of oxytocin conformation indicate that the tripeptide tail is more flexible than the tocin ring.     

All-atom molecular dynamics study of EAK16 peptide: the effect of pH on single-chain conformation,dimerization and self-assembly behavior

Single-chain equilibrium conformation and dimerization of the three types of ionic EAK16 peptide are studied under three pH conditions using all-atom molecular dynamics simulations. It is found that both the single-chain conformation and the dimerization process of EAK16-IV are considerably different from those of the two other types, EAK16-I and EAK16-II. The value of pH is found to have a stronger effect on the single-chain conformation and dimerization of EAK16-IV. It is shown that in addition to the charge pattern on the peptide chains, the size of the side chains of the charged amino acids plays role in the conformation of the peptide chains and their dimerization. The results shed light on the pH-dependent self-assembly behavior of EAK16 peptide in the bulk solution, which has been reported in the literature.     

去七肽胰岛素的分子动力学研究

本文用分子动力学的方法对去七肽胰岛素(DHPI)分子的构象进行了研究,首先用分子动力学方法对晶体胰岛素分子的构象能进行了优化,然后除去B链C端的最后七个残基(B24—B30),做分子动力学模拟,得到了DHPI的平衡构象和均方差波动。胰岛素分子的X射线晶体衍射结构和能量优化构象之间的均方根偏差为0.1;所得DHPI构象和胰岛素能量优化构象间C原子间的均方根偏差为1.8。变化最大的区域是A8—A10,A18—A21,B1—B41和B18—B23。     

Solution dynamics of the oligosaccharide moiety of ganglioside GM1: Comparison of solution conformations with the bound state conformation in association with cholera toxin B-pentamer

The solution dynamics of the oligosaccharide moiety of ganglioside G_M1 have been determined by use of a combination of ¹H rotating frame Overhauser effect measurements and restrained molecular dynamics simulations, It is found that the Galβ1-3 and NeuNAc moieties which are primarily recognized by cholera toxin both exhibit considerable torsional flexibility about their respective glycosidic linkages. A comparison with the bound state conformation of the ganglioside in association with cholera toxin B-pentamer, shows that a low energy conformation of the oligosaccharide, which closely approximates the globel minimum, is selected upon binding.     

Ligand-induced conformational and structural dynamics changes in Escherichia coli cyclic AMP receptor protein

Cyclic AMP receptor protein (CRP) regulates the expression of a large number of genes in E. coli. It is activated by cAMP binding, which leads to some yet undefined conformational changes. These changes do not involve significant redistribution of secondary structures. A potential mechanism of activation is a ligand-induced change in structural dynamics. Hence, the cAMP-mediated conformational and structural dynamics changes in the wild-type CRP were investigated using hydrogen-deuterium exchange and Fourier transform infrared spectroscopy. Upon cAMP binding, the two functional domains within the wild-type CRP undergo conformational and structural dynamics changes in two opposite directions. While the smaller DNA-binding domain becomes more flexible, the larger cAMP-binding domain shifts to a less dynamic conformation, evidenced by a faster and a slower amide H-D exchange, respectively. To a lesser extent, binding of cGMP, a nonfunctional analogue of cAMP, also stabilizes the cAMP-binding domain, but it fails to mimic the relaxation effect of cAMP on the DNA-binding domain. Despite changes in the conformation and structural dynamics, cAMP binding does not alter significantly the secondary structural composition of the wild-type CRP. The apparent difference between functional and nonfunctional analogues of cAMP is the ability of cAMP to effect an increase in the dynamic motions of the DNA binding domain.     

The effect of changes in gramicidin conformation on bilayer lipid properties.

The effects of two different gramicidin conformations on lipid phase behaviour and dynamics are compared. Samples of chain-perdeuterated dimyristoylphosphatidylcholine containing gramicidin were first prepared with gramicidin in a state having a circular dichroism spectrum generally identified as corresponding to the non-channel conformation. The effects, on bilayer lipid properties, of gramicidin in this conformation were then determined using deuterium nuclear magnetic resonance measurements of acyl chain orientational order and transverse relaxation times as a function of temperature. These samples were then incubated at 65 degrees C to convert the gramicidin to a state with a circular dichroism spectrum of the type generally identified with the channel conformation. The nuclear magnetic resonance measurements were then repeated. In the gel phase, it was found that transverse relaxation time and chain orientational order of the lipid were insensitive to gramicidin conformation. In the liquid crystalline phase, gramicidin in the channel conformation was found to have a slightly larger effect on transverse relaxation and orientational order than gramicidin in the non-channel conformation. The perturbation of the phase behavior by gramicidin was found to be relatively insensitive to gramicidin conformation.     

Quaternary structure dependence of kinetic hole burning and conformational substates interconversion in hemoglobin

Using a sol-gel encapsulation technique, we have prepared samples of CO saturated human adult hemoglobin locked in the R or T quaternary conformation. We report time-resolved spectra of these samples in the Soret region following flash photolysis, in the time interval ranging from 250 ns to 200 ms and in the temperature interval of 100-170 K. A suitable analysis of the measured difference spectra enables us to obtain the spectral contribution of deoxyHb and HbCO molecules as a function of time and/or of the fraction N(t) of deoxyHb molecules. In our experimental time window geminate CO rebinding to hemoglobin in the T quaternary conformation is about 2 orders of magnitude slower than to hemoglobin in the R conformation: this suggests that the barrier distribution for the CO rebinding, g(H), depends strongly on the protein quaternary structure. In our temperature interval, spectral shifts due to kinetic hole burning (KHB) are present: for HbCO the KHB effect is large in the R conformation and small in the T conformation. For deoxyHb the opposite is true. We attribute the observed behavior to the effect of interconversion between the relevant substates. This effect is stronger for HbCO molecules in the T conformation and for deoxyHb molecules in the R conformation; it confirms the quaternary structure dependence of the hemoglobin energy landscape and suggests enhanced dynamics of ligation intermediate species such as T-state HbCO or R-state deoxyHb.     

A computational study of the open and closed forms of the N-lobe human serum transferrin apoprotein

Human serum transferrin tightly binds ferric ions in the blood stream but is able to release them in cells by a process involving receptor-mediated endocytosis and decrease in pH. Iron binding and release are accompanied by a large conformation change. In this study, we investigate theoretically the open and closed forms of the N-lobe human serum transferrin apoprotein by performing pKa calculations and molecular dynamics and free-energy simulations. In agreement with the hypothesis based on the x-ray crystal structures, our calculations show that there is a shift in the pKa values of the lysines forming the dilysine trigger when the conformation changes. We argue, however, that simple electrostatic repulsion between the lysines is not sufficient to trigger domain opening and, instead, propose an alternative explanation for the dilysine-trigger effect. Analysis of the molecular dynamics and free-energy results indicate that the open form is more mobile than the closed form and is much more stable at pH 5.3, in large part due to entropic effects. Despite a lower free energy, the dynamics simulation of the open form shows that it is flexible enough to sample conformations that are consistent with iron binding.     

Fluorescence lifetime distributions in human superoxide dismutase. Effect of temperature and denaturation.

The internal dynamics of human superoxide dismutase has been studied using time-resolved fluorescence. The fluorescence decay has been analyzed using continuous distribution of lifetime values. The effect of temperature and conformational state on the lifetime distribution has been investigated. The emission of the single tryptophan residue depends on the nature and dynamics of the protein matrix. Conformational changes have been induced by increased concentration of guanidinium hydrochloride. We found that both temperature and conformation strongly effect the width of the lifetime distribution.     

Cofactor assisted gating mechanism in the active site of NADH oxidase from Thermus thermophilus

NADH oxidase (NOX) from Thermus thermophilus is a member of a structurally homologous flavoprotein family of nitroreductases and flavin reductases. The importance of local conformational dynamics in the active site of NOX has been recently demonstrated. The enzyme activity was increased by 250% in the presence of 1 M urea with no apparent perturbation of the native structure of the protein. The present in silico results correlate with the in vitro data and suggest the possible explanation about the effect of urea on NOX activity at the molecular level. Both, X-ray structure and molecular dynamics (MD) simulations, show open conformation of the active site represented by approximately 0.9 nm distance between the indole ring of Trp47 and the isoalloxazine ring of FMN412. In this conformation, the substrate molecule can bind in the active site without sterical restraints. MD simulations also indicate more stable conformation of the active site called "closed" conformation. In this conformation, Trp47 and the isoalloxazine ring of FMN412 are so close to each other (approximately 0.5 nm) that the substrate molecule is unable to bind between them without perturbing this conformation. The open/close transition of the active site between Trp47 and the flavin ring is accompanied by release of the "tightly" bound water molecule from the active site--cofactor assisted gating mechanism. The presence of urea in aqueous solutions of NOX prohibits closing of the active site and even unlocks the closed active site because of the concomitant binding of a urea molecule in the active site cavity. The binding of urea in the active site is stabilized by formation of one/two persistent hydrogen bonds involving the carbonyl group of the urea molecule. Our report represents the first MD study of an enzyme from the novel flavoprotein family of nitroreductases and flavin reductases. The common occurrence of aromatic residues covering the active sites in homologous enzymes suggests the possibility of a general gating mechanism and the importance of local dynamics within this flavoprotein family.     

One gene, two diseases and three conformations: molecular dynamics simulations of mutants of human prion protein at room temperature and elevated temperatures

Fatal familial insomnia (FFI) and Creutzfeldt-Jakob disease (CJD) are associated to the same mutation at codon 178 but differentiate into clinicopathologically distinct diseases determined by this mutation and a naturally occurring methionine-valine polymorphism at codon 129 of the prion protein gene. It has been suggested that the clinical and pathological difference between FFI and CJD is caused by different conformations of the prion protein. Using molecular dynamics (MD), we investigated the effect of the mutation at codon 178 and the polymorphism at codon 129 on prion protein dynamics and conformation at normal and elevated temperatures. Four model structures were examined with a focus on their dynamics and conformational changes. The results showed differences in stability and dynamics between polymorphic variants. Methionine variants demonstrated a higher stability than valine variants. Elongation of existing beta-sheets and formation of new beta-sheets was found to occur more readily in valine polymorphic variants. We also discovered the inhibitory effect of proline residue on existing beta-sheet elongation.     

Characterizing conformation changes in proteins through the torsional elastic response

The relationship between functional conformation changes and thermal dynamics of proteins is investigated with the help of the torsional network model (TNM), an elastic network model in torsion angle space that we recently introduced. We propose and test a null-model of “random” conformation changes that assumes that the contributions of normal modes to conformation changes are proportional to their contributions to thermal fluctuations. Deviations from this null model are generally small. When they are large and significant, they consist in conformation changes that are represented by very few low frequency normal modes and overcome small energy barriers. We interpret these features as the result of natural selection favoring the intrinsic protein dynamics consistent with functional conformation changes. These “selected” conformation changes are more frequently associated to ligand binding, and in particular phosphorylation, than to pairs of conformations with the same ligands. This deep relationship between the thermal dynamics of a protein, represented by its normal modes, and its functional dynamics can reconcile in a unique framework the two models of conformation changes, conformational selection and induced fit. The program TNM that computes torsional normal modes and analyzes conformation changes is available upon request. This article is part of a Special Issue entitled: The emerging dynamic view of proteins: Protein plasticity in allostery, evolution and self-assembly.     

Computer Simulation of Macromolecules

Computer simulation techniques are now an essential part of modern structural molecular biology. They are used in many different ways in order to study the conformation, dynamics and interactions of proteins and nucleic acids. In this paper, I shall review several of these applications and then focus on three specific areas, namely the conformation and dynamics of proteins including the use of free energy perturbation methods to study mutant proteins, the conformation and dynamics of DNA and DNA-drug complexes, and the use of computers with parallel architectures. Although simulation of molecules as large and complex as proteins and nucleic acids may be considered a grand challenge in itself, there are even greater challenges for the future.     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.