Hydrophobic hydration of amphipathic peptides. (65/16154)

Biomolecular surfaces and interfaces are commonly found with apolar character. The hydrophobic effect thus plays a crucial role in processes involving association with biomolecular surfaces in the cellular environment. By computer simulation, we compared the hydrogen bonding structures and energetics of the proximal hydration shells of the monomer and dimer from a recent study of an extrinsic membrane peptide, melittin. The two peptides were studied in their amphipathic alpha-helical forms, which possess extended hydrophobic surfaces characterized by different topography. The topography of the peptide-water interface was found to be critical in determining the enthalpic nature of hydrophobic hydration. This topographical dependence has far-reaching implications in the regulation of bioactivities in the presence of amphipathicity. This result also engenders reconsideration of the validity of using free energy parameters that depend solely on the chemical nature of constituent moieties in characterizing hydrophobic hydration of proteins and biomolecules in general.  (+info)

An alamethicin channel in a lipid bilayer: molecular dynamics simulations. (66/16154)

We present the results of 2-ns molecular dynamics (MD) simulations of a hexameric bundle of Alm helices in a 1-palmitoyl-2-oleoylphosphatidylcholine bilayer. These simulations explore the dynamic properties of a model of a helix bundle channel in a complete phospholipid bilayer in an aqueous environment. We explore the stability and conformational dynamics of the bundle in a phospholipid bilayer. We also investigate the effect on bundle stability of the ionization state of the ring of Glu18 side chains. If all of the Glu18 side chains are ionised, the bundle is unstable; if none of the Glu18 side chains are ionized, the bundle is stable. pKA calculations suggest that either zero or one ionized Glu18 is present at neutral pH, correlating with the stable form of the helix bundle. The structural and dynamic properties of water in this model channel were examined. As in earlier in vacuo simulations (Breed et al., 1996 .Biophys. J. 70:1643-1661), the dipole moments of water molecules within the pore were aligned antiparallel to the helix dipoles. This contributes to the stability of the helix bundle.  (+info)

Molecular dynamics of human methemoglobin: the transmission of conformational information between subunits in an alpha beta dimer. (67/16154)

Spectroscopic studies indicate an interaction of the distal histidine with the heme iron as well as the transmission of distal heme perturbations across the alpha1beta1 interface. Molecular dynamics simulations have been used to explain the molecular basis for these processes. Using a human methemoglobin alpha beta dimer, it has been shown that at 235 K after 61 ps, a rearrangement occurs in the alpha-chain corresponding to the formation of a bond with the distal histidine. This transition does not take place in the beta-chain during a 100-ps simulation and is reversed at 300 K. The absence of the distal histidine transition in the isolated chains and with the interface frozen indicate the involvement of the alphabeta interface. A detailed analysis of the simulation has been performed in terms of RMS fluctuations, domain cross-correlation maps, the disruption of helix hydrogen bonds, as well changes in electrostatic interactions and dihedral angles. This analysis shows that the rearrangements in the alpha-chain necessary to bring the histidine closer to the iron involve alterations primarily in the CD loop and at the interface. Communication to the beta-chain distal pocket is propagated by increased interactions of the alpha-chain B helix with the beta-chain G-GH-H segment and the flexibility in the EF loop. The G helices shown to be involved in propagation of perturbation across the alpha1beta1 interface extend into the alpha1beta2 interfaces, providing a mechansim whereby distal interactions can modulate the T<==>R transition in hemoglobin.  (+info)

Cardiac sodium channel Markov model with temperature dependence and recovery from inactivation. (68/16154)

A Markov model of the cardiac sodium channel is presented. The model is similar to the CA1 hippocampal neuron sodium channel model developed by Kuo and Bean (1994. Neuron. 12:819-829) with the following modifications: 1) an additional open state is added; 2) open-inactivated transitions are made voltage-dependent; and 3) channel rate constants are exponential functions of enthalpy, entropy, and voltage and have explicit temperature dependence. Model parameters are determined using a simulated annealing algorithm to minimize the error between model responses and various experimental data sets. The model reproduces a wide range of experimental data including ionic currents, gating currents, tail currents, steady-state inactivation, recovery from inactivation, and open time distributions over a temperature range of 10 degrees C to 25 degrees C. The model also predicts measures of single channel activity such as first latency, probability of a null sweep, and probability of reopening.  (+info)

Defining the transmembrane helix of M2 protein from influenza A by molecular dynamics simulations in a lipid bilayer. (69/16154)

Integral membrane proteins containing at least one transmembrane (TM) alpha-helix are believed to account for between 20% and 30% of most genomes. There are several algorithms that accurately predict the number and position of TM helices within a membrane protein sequence. However, these methods tend to disagree over the beginning and end residues of TM helices, posing problems for subsequent modeling and simulation studies. Molecular dynamics (MD) simulations in an explicit lipid and water environment are used to help define the TM helix of the M2 protein from influenza A virus. Based on a comparison of the results of five different secondary structure prediction algorithms, three different helix lengths (an 18mer, a 26mer, and a 34mer) were simulated. Each simulation system contained 127 POPC molecules plus approximately 3500-4700 waters, giving a total of approximately 18,000-21,000 atoms. Two simulations, each of 2 ns duration, were run for the 18mer and 26mer, and five separate simulations were run for the 34mer, using different starting models generated by restrained in vacuo MD simulations. The total simulation time amounted to 11 ns. Analysis of the time-dependent secondary structure of the TM segments was used to define the regions that adopted a stable alpha-helical conformation throughout the simulation. This analysis indicates a core TM region of approximately 20 residues (from residue 22 to residue 43) that remained in an alpha-helical conformation. Analysis of atomic density profiles suggested that the 18mer helix revealed a local perturbation of the lipid bilayer. Polar side chains on either side of this region form relatively long-lived H-bonds to lipid headgroups and water molecules.  (+info)

Noncontact dipole effects on channel permeation. II. Trp conformations and dipole potentials in gramicidin A. (70/16154)

The four Trp dipoles in the gramicidin A (gA) channel modulate channel conductance, and their side chain conformations should therefore be important, but the energies of different conformations are unknown. A conformational search for the right-handed helix based on molecular mechanics in vacuo yielded 46 conformations within 20 kcal/mol of the lowest energy conformation. The two lowest energy conformations correspond to the solid-state and solution-state NMR conformations, suggesting that interactions within the peptide determine the conformation. For representative conformations, the electrostatic potential of the Trp side chains on the channel axis was computed. A novel application of the image-series method of. Biophys. J. 9:1160-1170) was introduced to simulate the polarization of bulk water by the Trp side chains. For the experimentally observed structures, the CHARm toph19 potential energy (PE) of a cation in the channel center is -1.65 kcal/mol without images. With images, the PE is -1.9 kcal/mol, demonstrating that the images further enhance the direct dipole effect. Nonstandard conformations yielded less favorable PEs by 0.4-1.1 kcal/mol.  (+info)

Simulation analysis of the retinal conformational equilibrium in dark-adapted bacteriorhodopsin. (71/16154)

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.  (+info)

A possible resolution of the gating paradox. (72/16154)

We introduce a Markov model for the gating of membrane channels. The model features a possible solution to the so-called gating current paradox, namely that the bell-shaped curve that describes the voltage dependence of the kinetics is broader than expected from, and shifted relative to, the sigmoidal curve that describes the voltage dependence of the activation. The model also predicts some temperature dependence of this shift, but presence of the latter has not been tested experimentally so far.  (+info)