All-atom empirical potential for molecular modeling and dynamics studies of proteins

Structure, interaction, dynamics and solvent effects on the DNA-EcoRI complex in aqueous solution from molecular dynamics simulation

A 0.7-ns molecular dynamics simulation of the DNA-EcoRI complex in a 7.0-A solvent shell indicated a stable behavior of the system. No significant evaporation or smearing of the solvent's outer boundary occurred. The structure and the intermolecular interactions were found to be well maintained during the simulation. The interaction pattern in the simulation was found to be very similar to that in the crystal structure. Most of the specific interactions between the DNA and the protein were found to be enhanced in the simulation compared to that in the crystal structure as a result of improved interaction geometry. The nonspecific interactions were found to be stronger than the specific ones. The specific interactions between the N7 atoms of Gua(4) or Ade(5) or Ade(6) and the protein were found to be present over almost the entire time of the simulation, whereas hydrogen bonds involving the amino groups of the Ade(5) and Ade(6) with the protein were found to be relatively weaker, with lower probability and shorter lifetime. The time evolution of the root mean square deviations of the DNA and the protein were highly correlated even at the later part of the simulation, showing the tight binding between them. Several long-lived water bridges were found between the DNA backbone atoms and the protein and also between the two protein monomers, which increased the overall stability of the complex. The two protein monomers were found to interact strongly with each other. The energy of the DNA kink deformation was estimated as approximately 31 kcal/mol.

Steered molecular dynamics simulation of the Rieske subunit motion in the cytochrome bc(1) complex

Crystallographic structures of the mitochondrial ubiquinol/cytochrome c oxidoreductase (cytochrome bc(1) complex) suggest that the mechanism of quinol oxidation by the bc(1) complex involves a substantial movement of the soluble head of the Rieske iron-sulfur protein (ISP) between reaction domains in cytochrome b and cytochrome c(1) subunits. In this paper we report the results of steered molecular dynamics simulations inducing, through an applied torque within 1 ns, a 56 degrees rotation of the soluble domain of ISP. For this purpose, a solvated structure of the bc(1) complex in a phospholipid bilayer (a total of 206,720 atoms) was constructed. A subset of 91,061 atoms was actually simulated with 45,131 moving atoms. Point charge distributions for the force field parametrization of heme groups and the Fe(2)S(2) cluster of the Rieske protein included in the simulated complex were determined. The simulations showed that rotation of the soluble domain of ISP is actually feasible. Several metastable conformations of the ISP during its rotation were identified and the interactions stabilizing the initial, final, and intermediate positions of the soluble head of the ISP domain were characterized. A pathway for proton conduction from the Q(o) site to the solvent via a water channel has been identified.

Evidence for the role of PrP(C) helix 1 in the hydrophilic seeding of prion aggregates

Prions are mammalian proteins (PrPs) with a unique pathogenic property: a nonendogenous isoform PrP(Sc) can catalyze conversion of the endogenous PrP(C) isoform into additional PrP(Sc). In this work, we demonstrate that PrP(C) helix 1 has certain properties (hydrophilicity, charge distribution) that make it unique among all naturally occurring alpha-helices, and which are indicative of a highly specific model of prion infectivity. The beta-nucleation model proposes that PrP(Sc) is an aggregate with a hydrophilic core, consisting of a beta-sheet-like arrangement of constituent helix 1 components. It is suggested by using structural arguments, and confirmed by using CHARMM energy calculations, that aggregate formation from two PrP(C) molecules is highly unfavorable, but the addition of chains to an existing aggregate is favorable. The beta-nucleation model is shown to be consistent with the prion species-barrier, as well as with infectivity data. Sequence analysis of all known protein structures indicates that PrP is uniquely suited to beta-nucleation, in contrast to the many proteins that readily form less favorable (often nonspecific) hydrophobic aggregates.

The binding site of sodium in the gramicidin A channel

The available information concerning the structure and location of the main binding site for sodium in the gramicidin A channel is reviewed and discussed. Results from molecular dynamics simulations using an atomic model of the channel embedded in a lipid bilayer are compared with experimental observations. The combined information from experiment and simulation suggests that the main binding sites for sodium are near the channel's mouth, approximately 9.2 A from the centre of the dimer channel, although the motion along the axis could be as large as 1 to 2 A. In the binding site, the sodium ion is lying off axis, making contact with two carbonyl oxygens and two single-file water molecules. The main channel ligand is provided by the carbonyl group of the Leu10-Trp11 peptide linkage, which exhibits the largest deflection from the ion-free channel structure.

Molecular dynamics simulations of the complex between human U1A protein and hairpin II of U1 small nuclear RNA and of free RNA in solution

RNA-protein interactions are essential to a wide range of biological processes. In this paper, a 0.6-ns molecular dynamics simulation of the sequence-specific interaction of human U1A protein with hairpin II of U1 snRNA in solution, together with a 1.2-ns simulation of the free RNA hairpin, is reported. Compared to the findings in the x-ray structure of the complex, most of the interactions remained stable. The nucleotide U8, one of the seven conserved nucleotides AUUGCAC in the loop region, was unusually flexible during the simulation, leading to a loss of direct contacts with the protein, in contrast to the situation in the x-ray structure. Instead the sugar-phosphate backbone of nucleotide C15 was found to form several interactions with the protein. Compared to the NMR structure of U1A protein complexed with the 3'-untranslated region of its own pre-mRNA, the protein core kept the same conformation, and in the two RNA molecules the conserved AUUGCAC of the loop and the closest CG base pair were located in very similar positions and orientations, and underwent very similar interactions with the protein. Therefore, a common sequence-specific interaction mechanism was suggested for the two RNA substrates to bind to the U1A protein. Conformational analysis of the RNA hairpin showed that the conformational changes of the RNA primarily occurred in the loop region, which is just involved in the sites of binding to the protein and in agreement with experimental observation. Both the loop and stem of the RNA became more ordered upon binding to the protein. It was also demonstrated that the molecular dynamics method could be successfully used to simulate the dynamical behavior of a large RNA-protein complex in aqueous solution, thus opening a path for the exploration of the complex biological processes involving RNA at a molecular level.

Probing the role of structural water in a duplex oligodeoxyribonucleotide containing a water-mimicking base analog

Molecular dynamics simulations were performed on models of the dodecamer DNA double-stranded segment, [d(CGCGAATTCGCG)](2), in which each of the adenine residues, individually or jointly, was replaced by the water-mimicking analog 2'-deoxy-7-(hydroxy-methyl)-7-deazaadenosine (hm(7)c(7)dA) [Rockhill, J.K., Wilson,S.R. and Gumport,R.I. (1996) J. Am. Chem. Soc.,118, 10065-10068]. The simulations, when compared with those of the dodecamer itself, show that incorporation of the analog affects neither the overall DNA structure nor its hydrogen-bonding and stacking interactions when it replaces a single individual base. Furthermore, the water molecules near the bases in the singly-substituted oligonucleotides are similarly unaffected. Double substitutions lead to differences in all the aforementioned parameters with respect to the reference sequence. The results suggest that the analog provides a good mimic of specific 'ordered' water molecules observed in contact with DNA itself and at the interface between protein and DNA in specific complexes.

Structure and dynamics of alpha-MSH using DRISM integral equation theory and stochastic dynamics

The structural and dynamical features of the hormone alpha-MSH in solution have been examined over a 100 ns time scale by using free energy molecular mechanics models at room temperature. The free energy surface has been modeled using methods from integral equation theory and the dynamics by the Langevin equation. A modification of the accessible surface area friction drag model was used to calculate the atomic friction coefficients. The molecule shows a stable beta-turn conformation in the message region and a close interaction between the side chains of His6, Phe7, and Trp9. A salt bridge between Glu5 and Arg8 was found not to be a preferred interaction, whereas a Glu5 and Lys11 salt bridge was not sampled, presumably due to relatively high free energy barriers. The message region was more conformationally rigid than the N-terminal region. Several structural features observed here agree well with experimental results. The conformational features suggest a receptor-hormone interaction model where the hydrophobic side chains of Phe7 and Trp9 interact with the transmembrane portion of the MC1 receptor. Also, the positively charged side chain of Arg8 and the imidazole side chain of His6 may interact with the negatively charged portions of the receptor which may even be on the receptor's extracellular loops.

The interactions of substituted pyrido[1,2-e]purines with oligonucleotides depend on the amphiphilic properties of their side chain

Three pyrido[1,2-e]purines of increasing hydrophilicity have been synthesized to evaluate as anticancer agents. These drugs interact quite differently with a synthetic oligodeoxynucleotide d(CGATCG)2. [1] is very hydrophobic due to a phenyl residue in its side chain. It only shows limited interactions with the minihelix without any evidence of intercalation. [2] and [3], on the other hand, have one ([2]) or two ([3]) hydroxyl groups in their acyl chain and present rather amphiphilic properties. The result is a similar intercalation of these derivatives between C and G base pairs as revealed by intermolecular nOe, 1H and 31P chemical shift variations. Models for the intercalation of [2] are proposed using energy minimizations and molecular dynamics (MD) calculations subject to restraints from nOe connectivities. Simulations and experiments indicate weak stability and thus fast exchange of [2] in its intercalation site.

The cavity and pore helices in the KcsA K+ channel: electrostatic stabilization of monovalent cations

The electrostatic influence of the central cavity and pore alpha helices in the potassium ion channel from Streptomyces lividans (KcsA K+ channel) was analyzed by solving the finite difference Poisson equation. The cavity and helices overcome the destabilizing influence of the membrane and stabilize a cation at the membrane center. The electrostatic effect of the pore helices is large compared to that described for water-soluble proteins because of the low dielectric membrane environment. The combined contributions of the ion self-energy and the helix electrostatic field give rise to selectivity for monovalent cations in the water-filled cavity. Thus, the K+ channel uses simple electrostatic principles to solve the fundamental problem of ion destabilization by the cell membrane lipid bilayer.

Consequences of breaking the Asp-His hydrogen bond of the catalytic triad: effects on the structure and dynamics of the serine esterase cutinase

The objective of this study has been to investigate the effects on the structure and dynamics that take place with the breaking of the Asp-His hydrogen bond in the catalytic triad Asp175-His188-Ser120 of the serine esterase cutinase in the ground state. Four molecular dynamics simulations were performed on this enzyme in solution. The starting structures in two simulations had the Asp175-His188 hydrogen bond intact, and in two simulations the Asp175-His188 hydrogen bond was broken. Conformations of the residues comprising the catalytic triad are well behaved during both simulations containing the intact Asp175-His188 hydrogen bond. Short contacts of less than 2.6 A were observed in 1.2% of the sampled distances between the carboxylate oxygens of Asp175 and the NE2 of His188. The simulations showed that the active site residues exhibit a great deal of mobility when the Asp175-His188 hydrogen bond is broken. In the two simulations in which the Asp175-His188 hydrogen bond is not present, the final geometries for the residues in the catalytic triad are not in catalytically productive conformations. In both simulations, Asp175 and His188 are more than 6 A apart in the final structure from dynamics, and the side chains of Ser120 and Asp175 are in closer proximity to the NE2 of His188 than to ND1. Nonlocal effects on the structure of cutinase were observed. A loop formed by residues 26-31, which is on the opposite end of the protein relative to the active site, was greatly affected. Further changes in the dynamics of cutinase were determined from quasiharmonic mode analysis. The frequency of the second lowest mode was greatly reduced when the Asp175-His188 hydrogen bond was broken, and several higher modes showed lower frequencies. All four simulations showed that the oxyanion hole, composed of residues Ser42 and Gln121, is stable. Only one of the hydrogen bonds (Ser42 OG to Gln121 NE2) observed in the crystal structure that stabilize the conformation of Ser42 OG persisted throughout the simulations. This hydrogen bond appears to be enough for the oxyanion hole to retain its structural integrity.

Binding of aldose reductase inhibitors: correlation of crystallographic and mass spectrometric studies

Aldose reductase is a NADP(H)-dependent enzyme, believed to be strongly implicated in the development of degenerative complications of Diabetes Mellitus. The search for specific inhibitors of this enzyme has thus become a major pharmaceutic challenge. In this study, we applied both X-ray crystallography and mass spectrometry to characterize the interactions between aldose reductase and four representative inhibitors: AminoSNM, Imirestat, LCB3071, and IDD384. If crystallography remains obviously the only way to get an extensive description of the contacts between an inhibitor and the enzymatic site, the duration of the crystallographic analysis makes this technique incompatible with high throughput screenings of inhibitors. On the other hand, dissociation experiments monitored by mass spectrometry permitted us to evaluate rapidly the relative gas-phase stabilities of the aldose reductase-inhibitor noncovalent complexes. In our experiments, dissociation in the gas-phase was provoked by increasing the accelerating voltage of the ions (Vc) in the source-analyzer interface region: the Vc value needed to dissociate 50% of the noncovalent complex initially present (Vc50) was taken as a gas-phase stability parameter of the enzyme-inhibitor complex. Interestingly, the Vc50 were found to correlate with the energy of the electrostatic and H-bond interactions involved in the contact aldose reductase/inhibitor (Eel-H), computed from the crystallographic model. This finding may be specially interesting in a context of drug development. Actually, during a drug design optimization phase, the binding of the drug to the target enzyme is often optimized by modifying its interatomic electrostatic and H-bond contacts; because they usually depend on a single atom change on the drug, and are easier to introduce than the hydrophobic interactions. Therefore, the Vc50 may help to monitor the chemical modifications introduced in new inhibitors. X-ray crystallography is clearly needed to get the details of the contacts and to rationalize the design. Nevertheless, once the cycle of chemical modification is engaged, mass spectrometry can be used to select a priori the drug candidates which are worthy of further crystallographic investigation. We thus propose to use the two techniques in a complementary way, to improve the screening of large collections of inhibitors.

Statistical mechanical equilibrium theory of selective ion channels

A rigorous statistical mechanical formulation of the equilibrium properties of selective ion channels is developed, incorporating the influence of the membrane potential, multiple occupancy, and saturation effects. The theory provides a framework for discussing familiar quantities and concepts in the context of detailed microscopic models. Statistical mechanical expressions for the free energy profile along the channel axis, the cross-sectional area of the pore, and probability of occupancy are given and discussed. In particular, the influence of the membrane voltage, the significance of the electric distance, and traditional assumptions concerning the linearity of the membrane electric field along the channel axis are examined. Important findings are: 1) the equilibrium probabilities of occupancy of multiply occupied channels have the familiar algebraic form of saturation properties which is obtained from kinetic models with discrete states of denumerable ion occupancy (although this does not prove the existence of specific binding sites; 2) the total free energy profile of an ion along the channel axis can be separated into an intrinsic ion-pore free energy potential of mean force, independent of the transmembrane potential, and other contributions that arise from the interfacial polarization; 3) the transmembrane potential calculated numerically for a detailed atomic configuration of the gramicidin A channel embedded in a bilayer membrane with explicit lipid molecules is shown to be closely linear over a distance of 25 A along the channel axis. Therefore, the present analysis provides some support for the constant membrane potential field approximation, a concept that has played a central role in the interpretation of flux data based on traditional models of ion permeation. It is hoped that this formulation will provide a sound physical basis for developing nonequilibrium theories of ion transport in selective biological channels.

Structural characterization by computer experiments of the lipid-free LDL-receptor-binding domain of apolipoprotein E

The structure and dynamics of the lipid-free LDL-receptor-binding domain of apolipoprotein E (apoE-RBD) has been investigated by Molecular Dynamics Simulations. ApoE-RBD in its monomeric lipid-free form is a singular four-helix bundle made up of four elongated amphipathic helices. Analysis of one 1.5 ns molecular dynamics trajectory of apoE-RBD performed in water indicates that the lipid-free domain adopts a structure that exhibits characteristics found in native proteins: it has very stable helices and presents a compact structure. Yet its interior exhibits a larger number of transient atomic-size cavities relative to that found in other proteins of similar size and its apolar side chains are more mobile. The latter features distinguish the elongated four-helix bundle as a slightly disordered structure, which shows a structural likeness with some de novo designed four-helix bundle proteins and shares with the latter a leucine-rich residue composition. We anticipate that these unique properties compared with other native helix bundles may be related to the postulated ability of apoE-RBD to undergo an opening of its bundle upon interaction with phospholipids. The distribution of empty cavities computed along the trajectory in the interface regions between the different pairs of helices reveals that the tertiary contacts in one of the interfaces are weaker suggesting that this particular interface could be more easily ruptured upon lipid association.

A method for computational combinatorial peptide design of inhibitors of Ras protein

A computational combinatorial approach is proposed for the design of a peptide inhibitor of Ras protein. The procedure involves three steps. First, a 'Multiple Copy Simultaneous Search' identifies the location of specific functional groups on the Ras surface. This search method allowed us to identify an important binding surface consisting of two beta strands (residues 5-8 and 52-56), in addition to the well known Ras effector loop and switch II region. The two beta strands had not previously been reported to be involved in Ras-Raf interaction. Second, after constructing the peptide inhibitor chain based on the location of N-methylacetamide (NMA) minima, functional groups are selected and connected to the main chain Calpha atom. This step generates a number of possible peptides with different sequences on the Ras surface. Third, potential inhibitors are designed based on a sequence alignment of the peptides generated in the second step. This computational approach reproduces the conserved pattern of hydrophobic, hydrophilic and charged amino acids identified from the Ras effectors. The advantages and limitations of this approach are discussed.

Structural details of urea binding to barnase: a molecular dynamics analysis

BACKGROUND

The molecular mechanism of urea-induced protein unfolding has not been established. It is generally thought that denaturation results from the stabilizing interactions of urea with portions of the protein that are buried in the native state and become exposed upon unfolding of the protein.

RESULTS

We have performed molecular dynamics simulations of barnase (a 110 amino acid RNase from Bacillus amyloliquefaciens) with explicit water and urea molecules at 300 K and 360 K. The native conformation was unaffected in the 300 K simulations at neutral and low pH. Two of the three runs at 360 K and low pH showed some denaturation, with partial unfolding of the hydrophobic core 2. The first solvation shell has a much higher density of urea molecules (water/urea ratio ranging from 2.07 to 2.73) than the bulk (water/urea ratio of 4.56). About one half of the first-shell urea molecules are involved in hydrogen bonds with polar or charged groups on the barnase surface, and between 15% and 18% of the first-shell urea molecules participate in multiple hydrogen bonds with barnase. The more stably bound urea molecules tend to be in crevices or pockets on the barnase surface.

CONCLUSIONS

The simulation results indicate that an aqueous urea solution solvates the surface of a polypeptide chain more favorably than pure water. Urea molecules interact more favorably with nonpolar groups of the protein than water does, and the presence of urea improves the interactions of water molecules with the hydrophilic groups of the protein. The results suggest that urea denaturation involves effects on both nonpolar and polar groups of proteins.

Dramatic Influence of Vβ Gene Polymorphism on an Antigen-Specific CD8+ T Cell Response In Vivo

According to recent crystallographic studies, the TCR-αβ contacts MHC class I-bound antigenic peptides via the polymorphic V gene-encoded complementarity-determining region 1β (CDR1β) and the hypervariable (D)J-encoded CDR3β and CDR3α domains. To evaluate directly the relative importance of CDR1β polymorphism on the fine specificity of T cell responses in vivo, we have taken advantage of congenic Vβa and Vβb mouse strains that differ by a CDR1 polymorphism in the Vβ10 gene segment. The Vβ10-restricted CD8+ T cell response to a defined immunodominant epitope was dramatically reduced in Vβa compared with Vβb mice, as measured either by the expansion of Vβ10+ cells or by the binding of MHC-peptide tetramers. These data indicate that Vβ polymorphism has an important impact on TCR-ligand binding in vivo, presumably by modifying the affinity of CDR1β-peptide interactions.

Comparison of the potentials of mean force for alanine tetrapeptide between integral equation theory and simulation

The dielectrically consistent reference interaction site model (DRISM) integral equation theory is applied to determine the potential of mean force (PMF) for an alanine tetramer. A stochastic dynamics simulation of the alanine tetramer using this PMF is then compared with an explicit water molecular dynamics simulation. In addition, comparison is also done with simulations using other solvent models like the extended reference interaction site model (XRISM) theory, constant dielectric and linear distance-dependent dielectric models. The results show that the DRISM method offers a fairly accurate and computationally inexpensive alternative to explicit water simulations for studies on small peptides.

Identification of substrate binding site of cyclin-dependent kinase 5

Cyclin-dependent kinase 5 (CDK5), unlike other CDKs, is active only in neuronal cells where its neuron-specific activator p35 is present. However, it phosphorylates serines/threonines in S/TPXK/R-type motifs like other CDKs. The tail portion of neurofilament-H contains more than 50 KSP repeats, and CDK5 has been shown to phosphorylate S/T specifically only in KS/TPXK motifs, indicating highly specific interactions in substrate recognition. CDKs have been shown to have a high preference for a basic residue (lysine or arginine) as the n+3 residue, n being the location in the primary sequence of a phosphoacceptor serine or threonine. Because of the lack of a crystal structure of a CDK-substrate complex, the structural basis for this specific interaction is unknown. We have used site-directed mutagenesis ("charged to alanine") and molecular modeling techniques to probe the recognition interactions for substrate peptide (PKTPKKAKKL) derived from histone H1 docked in the active site of CDK5. The experimental data and computer simulations suggest that Asp86 and Asp91 are key residues that interact with the lysines at positions n+2 and/or n+3 of the substrates.

Effect of G40R mutation on the binding of human SRY protein to DNA: a molecular dynamics view

Molecular dynamics simulation was conducted to investigate the reason why the mutant G40R of hSRY protein has a low affinity for DNA. Compared with the previous dynamics results of the wild-type hSRY-HMG-DNA complex, the results of molecular dynamics simulation on the mutant G40R hSRY-HMG-DNA system demonstrated that the whole structure of DNA (especially the second strand) had a major deviation away from the short arm of the HMG box. Consequently, the DNA and the mutant protein could not specifically recognize each other, that is, very different, and low-occupancy, direct, and water-mediated hydrogen bonds were detected at the protein-DNA interface, no conformational changes occurred at the loop region around Met9 during the simulation, and residue IIe13 did not intercalate between the bases of A5 and A6. These results indicated that the mutant G40R did not form a specific complex with the DNA target, hence led to complete gonadal dysgenesis. From the simulation, we realized that the residue Gly40 played a critical structural role in the hSRY-DNA recognition. It might be a structural supporting point of DNA binding because of the absence of a side chain. The reason for the difficulty of the mutant G40R to form a complex with DNA might be that the long and positively charged side chain of Arg40 by its bulk and positive charge hindered the DNA's access to the active sites of the protein.

Towards a better description and understanding of biomolecular solvation

We introduce a flexible framework for the correct description of the solvation of biological macromolecules, the dielectric field equation (DFE). The formalism permits the use of any combination of quantum mechanical (QM), molecular mechanical (MM) and continuum electrostatic (CE) based techniques. For the CE region a method that yields the electric field rather than the potential is outlined. The DFE formalism makes clear the need to consider and to calibrate a dielectric boundary region surrounding the simulation system. The details of how to do this are presented for the case of the Ewald summation method; the effects are demonstrated by calculations of the dielectric properties and the spatially resolved Kirkwood G-factor, G(K)(r), of TIP3P water. Computing the dielectric properties of a multi-component system provides a sensitive method to better understand the solvation of biological macromolecules. Towards this goal a rigorous analysis of the dielectric properties of solvated biomolecules based on a decomposition of the frequency-dependent dielectric constant (or susceptibility) of the full system is presented. The meaning of our approach is investigated, and the results of a first application are reported. Using the method of Voronoi polyhedra, the dielectric properties of the first two solvation shells and bulk water are compared by re-analyzing a 12-ns trajectory of a zinc finger peptide in water [Löffler et al. J. Mol. Biol. 270 (1997) 520]. It is found that the first shell behaves considerably different; in addition, there is a non-negligible contribution to the total susceptibility of the system from coupling between the protein and the bulk water phase, i.e. the water molecules not in the immediate vicinity of the solute.

Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin

In dark-adapted bacteriorhodopsin (bR) the retinal moiety populates two conformers: all-trans and (13,15)cis. Here we examine factors influencing the thermodynamic equilibrium and conformational transition between the two forms, using molecular mechanics and dynamics calculations. Adiabatic potential energy mapping indicates that whereas the twofold intrinsic torsional potentials of the C13==C14 and C15==N16 double bonds favor a sequential torsional pathway, the protein environment favors a concerted, bicycle-pedal mechanism. Which of these two pathways will actually occur in bR depends on the as yet unknown relative weight of the intrinsic and environmental effects. The free energy difference between the conformers was computed for wild-type and modified bR, using molecular dynamics simulation. In the wild-type protein the free energy of the (13,15)cis retinal form is calculated to be 1.1 kcal/mol lower than the all-trans retinal form, a value within approximately kBT of experiment. In contrast, in isolated retinal the free energy of the all-trans state is calculated to be 2.1 kcal/mol lower than (13,15)cis. The free energy differences are similar to the adiabatic potential energy differences in the various systems examined, consistent with an essentially enthalpic origin. The stabilization of the (13,15)cis form in bR relative to the isolated retinal molecule is found to originate from improved protein-protein interactions. Removing internal water molecules near the Schiff base strongly stabilizes the (13,15)cis form, whereas a double mutation that removes negative charges in the retinal pocket (Asp85 to Ala; Asp212 to Ala) has the opposite effect.

Molecular dynamics study of substance P peptides partitioned in a sodium dodecylsulfate micelle

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in an explicit sodium dodecylsulfate (SDS) micelle. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) during the equilibration stage of the simulation. It was shown that when SP-Y8 was initially placed in an insertion (perpendicular) configuration, the peptide equilibrated to a surface-bound (parallel) configuration in approximately 450 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained from the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbones of all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions come from the hydrophobic interaction between the side chains of Lys-3, Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, from Lys-3, Phe-7, Leu-10, and Met-11 in SP-Y8, and the micellar interior. Significant interactions, electrostatic and hydrogen bonding, between the N-terminal residues, Arg-Pro-Lys, and the micellar headgroups were observed. These latter interactions served to affect both the structure and, especially, the flexibility, of the N-terminus. The results from simulation of the same peptides in a water/CCl4 biphasic cell were compared with the results of the present study, and the validity of using the biphasic system as an approximation for peptide-micelle or peptide-bilayer systems is discussed.

Molecular dynamics study of substance P peptides in a biphasic membrane mimic

Two neuropeptides, substance P (SP) and SP-tyrosine-8 (SP-Y8), have been studied by molecular dynamics (MD) simulation in a TIP3P water/CCl4 biphasic solvent system as a mimic for the water-membrane system. Initially, distance restraints derived from NMR nuclear Overhauser enhancements (NOE) were incorporated in the restrained MD (RMD) in the equilibration stage of the simulation. The starting orientation/position of the peptides for the MD simulation was either parallel to the water/CCl4 interface or in a perpendicular/insertion mode. In both cases the peptides equilibrated and adopted a near-parallel orientation within approximately 250 ps. After equilibration, the conformation and orientation of the peptides, the solvation of both the backbone and the side chain of the residues, hydrogen bonding, and the dynamics of the peptides were analyzed from trajectories obtained in the RMD or the subsequent free MD (where the NOE restraints were removed). These analyses showed that the peptide backbone of nearly all residues are either solvated by water or are hydrogen-bonded. This is seen to be an important factor against the insertion mode of interaction. Most of the interactions with the hydrophobic phase come from the hydrophobic interactions of the side chains of Pro-4, Phe-7, Phe-8, Leu-10, and Met-11 for SP, and Phe-7, Leu-10, Met-11 and, to a lesser extent, Tyr-8 in SP-Y8. Concerted conformational transitions took place in the time frame of hundreds of picoseconds. The concertedness of the transition was due to the tendency of the peptide to maintain the necessary secondary structure to position the peptide properly with respect to the water/CCl4 interface.