Comparative simulations of the ground state and the M-intermediate state of the sensory rhodopsin II-transducer complex with a HAMP domain model

The complex of sensory rhodopsin II (SRII) and its cognate transducer HtrII (2:2 SRII-HtrII complex) consists of a photoreceptor and its signal transducer, respectively, associated with negative phototaxis in extreme halophiles. In this study to investigate how photoexcitation in SRII affects the structures of the complex, we conducted two series of molecular dynamics simulations of the complex of SRII and truncated HtrII (residues 1-136) of Natronomonas pharaonis linked with a modeled HAMP domain in the lipid bilayer using the two crystal structures of the ground state and the M-intermediate state as the starting structures. The simulation results showed significant enhancements of the structural differences observed between the two crystal structures. Helix F of SRII showed an outward motion, and the C-terminal end of transmembrane domain 2 (TM2) in HtrII rotated by ～10°. The most significant structural changes were observed in the overall orientations of the two SRII molecules, closed in the ground state and open in the M-state. This change was attributed to substantial differences in the structure of the four-helix bundle of the HtrII dimer causing the apparent rotation of TM2. These simulation results established the structural basis for the various experimental observations explaining the structural differences between the ground state and the M-intermediate state.     

Nano breathers and molecular dynamics simulations in hydrogen-bonded chains

Non-linear localization phenomena in biological lattices have attracted a steadily growing interest and their existence has been predicted in a wide range of physical settings. We investigate the non-linear proton dynamics of a hydrogen-bonded chain in a semi-classical limit using the coherent state method combined with a Holstein–Primakoff bosonic representation. We demonstrate that even a weak inherent discreteness in the hydrogen-bonded (HB) chain may drastically modify the dynamics of the non-linear system, leading to instabilities that have no analog in the continuum limit. We suggest a possible localization mechanism of polarization oscillations of protons in a hydrogen-bonded chain through modulational instability analysis. This mechanism arises due to the neighboring proton–proton interaction and coherent tunneling of protons along hydrogen bonds and/or around heavy atoms. We present a detailed analysis of modulational instability, and highlight the role of the interaction strength of neighboring protons in the process of bioenergy localization. We perform molecular dynamics simulations and demonstrate the existence of nanoscale discrete breather (DB) modes in the hydrogen-bonded chain. These highly localized and long-lived non-linear breather modes may play a functional role in targeted energy transfer in biological systems.     

Exploring the structure and dynamics of nano-confined water molecules using molecular dynamics simulations

Confinement effects can lead to drastic changes in the structural and dynamical properties of water molecules. In this work, we have performed classical molecular dynamics simulations of endohedral fullerenes of type (H₂O)_n@C_m (n = 1, 12, 21, 62, 108 and m = 60, 180, 240, 500 and 720) to explore the effects of spherical confinement on water properties. It is shown that these confined water molecules can form distinct solvation pattern depending upon the available space inside the fullerene cavity. For the systems with smaller diameter, cage-like structure is predominant whereas bulk-like structure is observed for larger fullerenes. The orientational relaxation of these confined water molecules showed slower relaxation as the cavity diameter increases except for the (H₂O)₂₁@C₂₄₀. In this case, stable cage-like structure hinders the overall dynamics of the trapped water molecules. Finally, we have calculated the hydrogen bond lifetimes from the hydrogen bond time correlation functions and compared with that of bulk water.     

Molecular dynamics simulations of migration of ions and molecules through the acetylcholine receptor channel

A dynamic model of the channel of an acetylcholine receptor in a closed state has been proposed. The channel is formed by five a-helices of subunit M2 and stabilized by the cyclic hydrocarbon (CH2)105. The migration of charged and unchanged van der Waals particles with a diameter of 7.72 A equivalent to the diameter of a hydrated sodium ion has been studied. The migration occurred by the action of external force applied to the complex along the channel axis. In the closed state, the inhibition of ions is due to two components: electrostatic interaction and steric constraints. The van der Waals channel gate is formed by residues 13'-A-Val255, B-Val261, C-Val269, D-Val255, and E-Ile264, and the negatively changed residues occurring in the upper part of the channel have a great effect on ion selectivity.     

Calculation of NMR-relaxation parameters for flexible molecules from molecular dynamics simulations

Comparatively small molecules such as peptides can show a high internal mobility with transitions between several conformational minima and sometimes coupling between rotational and internal degrees of freedom. In those cases the interpretation of NMR relaxation data is difficult and the use of standard methods for structure determination is questionable. On the other hand, in the case of those system sizes, the timescale of both rotational and internal motions is accessible by molecular dynamics (MD) simulations using explicit solvent. Thus a comparison of distance averages (r ^–6–1/6 or r ^–31/3) over the MD trajectory with NOE (or ROE) derived distances is no longer necessary, the (back)calculation of the complete spectra becomes possible. In the present study we use two 200 ns trajectories of a heptapeptide of -amino acids in methanol at two different temperatures to obtain theoretical ROESY spectra by calculating the exact spectral densities for the interproton vectors and the full relaxation matrix. Those data are then compared with the experimental ones. This analysis permits to test some of the assumptions and approximations that generally have to be made to interpret NMR spectra, and to make a more reliable prediction of the conformational equilibrium that leads to the experimental spectrum.     

Characterization of the ground state dynamics of proteorhodopsin by NMR and optical spectroscopies

We characterized the dynamics of proteorhodopsin (PR), solubilized in diC7PC, a detergent micelle, by liquid-state NMR spectroscopy at T?=?323?K. Insights into the dynamics of PR at different time scales could be obtained and dynamic hot spots could be identified at distinct, functionally relevant regions of the protein, including the BC loop, the EF loop, the N-terminal part of helix F and the C-terminal part of helix G. We further characterize the dependence of the photocycle on different detergents (n-Dodecyl ??-D-maltoside DDM; 1,2-diheptanoyl-sn-glycero-3-phosphocholine diC7PC) by ultrafast time-resolved UV/VIS spectroscopy. While the photocycle intermediates of PR in diC7PC and DDM exhibit highly similar spectral characteristics, significant changes in the population of these intermediates are observed. In-situ NMR experiments have been applied to characterize structural changes during the photocycle. Light-induced chemical shift changes detected during the photocycle in diC7PC are very small, in line with the changes in the population of intermediates in the photocycle of proteorhodopsin in diC7PC, where the late O-intermediate populated in DDM is missing and the population is shifted towards an equilibrium of intermediates states (M, N, O) without accumulation of a single populated intermediate.     

Molecular dynamics simulations of retinal in rhodopsin: from the dark-adapted state towards lumirhodopsin

The formation of photointermediates and conformational changes observed in the retinal chromophore of bilayer-embedded rhodopsin during the early steps of the protein activation have been studied by molecular dynamics (MD) simulation. In particular, the lysine-bound retinal has been examined, focusing on its conformation in the dark-adapted state (10 ns) and on the early steps after the isomerization of the 11-cis bond to trans (up to 10 ns). The parametrization for the chromophore is based on a recent quantum study [Sugihara, M., Buss, V., Entel, P., Elstner, M., and Frauenheim, T. (2002) Biochemistry 41, 15259-15266] and shows good conformational agreement with recent experimental results. The isomerization, induced by switching the function governing the dihedral angle for the C11=C12 bond, was repeated with several different starting conformations. From the repeated simulations, it is shown that the retinal model exhibits a conserved activation pattern. The conformational changes are sequential and propagate outward from the C11=C12 bond, starting with isomerization of the C11=C12 bond, then a rotation of methyl group C20, and followed by increased fluctuations at the beta-ionone ring. The dynamics of these changes suggest that they are linked with photointermediates observed by spectroscopy. The exact moment when these events occur after the isomerization is modulated by the starting conformation, suggesting that retinal isomerizes through multiple pathways that are slightly different. The amplitudes of the structural fluctuations observed for the protein in the dark-adapted state and after isomerization of the retinal are similar, suggesting a subtle mechanism for the transmission of information from the chromophore to the protein.     

Water molecules in hydroxy/acid networks as a competition between dynamics and bonding. Synthesis of a wet hydrophobic pore

In a model formed by hydroxy acids with a general structure (+/-)-1, we found that solid-state structures depend on steric interactions. Thus, with the exception of molecules 1b and 1e, compounds (+/-)-1a-(+/-)-1m, which possess bulky and conformationally rigid substituents, aggregate by forming tapes and sheets by alternating (+) and (-) subunits held together via carboxylic acid to alcohol hydrogen bonds. Homologue (+/-)-1n with conformationally flexible substituents, which allow conformational deformation gives, by way of the incorporation of water molecules, an efficient hexagonal assembly, which extends to the third-dimension to form tubular H-bonding networks. Each puckered channel can be described as being interconnected by closely packed hexagons in chair-like conformations. The ethyl groups presented in (+/-)-1n provided the volume required to lock the inner hexagonal wall into a rigid structure.     

Interactions of liquid crystal-forming molecules with phospholipid bilayers studied by molecular dynamics simulations

Recent experiments have shown that liquid crystals can be used to image mammalian cell membranes and to amplify structural reorganization in phospholipid-laden liquid crystal-aqueous interfaces. In this work, molecular dynamics simulations were employed to explore the interactions between commonly used liquid crystal-forming molecules and phospholipid bilayers. In particular, umbrella sampling was used to obtain the potential of mean force of 4-cyano-4'-pentylbiphenyl (5CB) and 4'-(3,4-difluor-phenyl)-4-pentyl-bicylohexyl (5CF) molecules partitioning into a dipalmitoylphosphatidylcholine bilayer. In addition, results of simulations are presented for systems consisting of a fully hydrated bilayer with 5CB or 5CF molecules at the lowest (4.5 mol %) and highest (20 mol %) concentrations used in recent laboratory experiments. It is found that mesogens preferentially partition from the aqueous phase into the membrane; the potential of mean force exhibits highly favorable free energy differences for partitioning (-18 k(B)T for 5CB and -26 k(B)T for 5CF). The location and orientation of mesogens associated with the most stable free energies in umbrella sampling simulations of dilute systems were found to be consistent with those observed in liquid-crystal-rich bilayers. It is found that the presence of mesogens in the bilayer enhances the order of lipid acyl tails, and changes the spatial and orientational arrangement of lipid headgroup atoms. These effects are more pronounced at higher liquid-crystal concentrations. In comparing the behavior of 5CB and 5CF, a stronger spatial correlation (i.e., possibly leading to aggregation) is observed between 5CB molecules within a bilayer than between 5CF molecules. Also, the range of molecular orientations and positions along the bilayer normal is larger for 5CB molecules. At the same time, 5CF molecules were found to bind more strongly to lipid headgroups, thereby slowing the lateral motion of lipid molecules.     

Molecular dynamics simulations on the aggregation behavior of indole type organic dye molecules in dye-sensitized solar cells

In Ti0₂ nanostructured dye-sensitized solar cells indole based organic dyes D149, D205 exhibits greater power conversion efficiency. Such organic dye molecules are easily undergone for aggregation. Aggregation in dye molecules leads to reduce electron transfer process in dye-sensitized solar cells. Therefore, anti-aggregating agents such as chenodeoxycholic acid are commonly added to organic dye solution in DSSCs. Studying aggregation of such dye molecules in the absence of semiconductors gives a detailed influence of anti-aggregating agents on dye molecules. Atomistic level of molecular dynamics (MD) simulations were performed on aggregation of indole type dye molecules D149, D205 and D205-F with anti-aggregating agent chenodeoxy cholic acid using AMBER program. The trajectories of the MD simulations were analyzed with order parameters such as radial atom pair distribution functions g(r), diffusion coefficients and root mean square deviations values. MD results suggest that addition of chenodeoxy cholic acid to dyes significantly reduces structural arrangement and increases conformational flexibility and mobility of dye molecules. The influence of semi-perfluorinated alkyl chains in indole dye molecules was analyzed. The parameters such as open-circuit voltage (V_oc) and power conversion efficiency (η) of dye-sensitized solar cells are corroborated with flexibility and diffusion values of dye molecules.     

Strongly hydrogen-bonded water molecules in the Schiff base region of rhodopsins.

In many rhodopsins, a positively charged retinal chromophore is stabilized by a negatively charged carboxylate, and the presence of bound water molecules has been found in the Schiff base region by X-ray crystallography of various rhodopsins. Low-temperature Fourier-transform infrared (FTIR) spectroscopy can directly monitor hydrogen-bonding alterations of internal water molecules of rhodopsins. In particular, we found that a bridged water molecule between the Schiff base and Asp 85 in bacteriorhodopsin (BR), a light-driven proton-pump protein, forms an extremely strong hydrogen bond. It is likely that a hydration switch of the water from Asp 85 to Asp 212 plays an important role in the proton transfer in the Schiff base region of BR. Comprehensive studies of archaeal and visual rhodopsins have revealed that strongly hydrogen-bonded water molecules are only found in the proteins exhibiting proton-pump activities. Strongly hydrogen-bonded water molecules and its transient weakening may be essential for the proton-pump function of rhodopsins.     

Molecular dynamics simulations of trehalose as a 'dynamic reducer' for solvent water molecules in the hydration shell

Systematic computational work for a series of 13 disaccharides was performed to provide an atomic-level insight of unique biochemical role of the alpha,alpha-(1-->1)-linked glucopyranoside dimer over the other glycosidically linked sugars. Superior osmotic and cryoprotective abilities of trehalose were explained on the basis of conformational and hydration characteristics of the trehalose molecule. Analyses of the hydration number and radial distribution function of solvent water molecules showed that there was very little hydration adjacent to the glycosidic oxygen of trehalose and that the dynamic conformation of trehalose was less flexible than any of the other sugars due to this anisotropic hydration. The remarkable conformational rigidity that allowed trehalose to act as a sugar template was required for stable interactions with hydrogen-bonded water molecules. Trehalose made an average of 2.8 long-lived hydrogen bonds per each MD step, which was much larger than the average of 2.1 for the other sugars. The stable hydrogen-bond network is derived from the formation of long-lived water bridges at the expense of decreasing the dynamics of the water molecules. Evidence for this dynamic reduction of water by trehalose was also established based on each of the lowest translational diffusion coefficients and the lowest intermolecular coulombic energy of the water molecules around trehalose. Overall results indicate that trehalose functions as a 'dynamic reducer' for solvent water molecules based on its anisotropic hydration and conformational rigidity, suggesting that macroscopic solvent properties could be modulated by changes in the type of glycosidic linkages in sugar molecules.     

Water transport in aquaporins: osmotic permeability matrix analysis of molecular dynamics simulations

Single-channel osmotic water permeability (p(f)) is a key quantity for investigating the transport capability of the water channel protein, aquaporin. However, the direct connection between the single scalar quantity p(f) and the channel structure remains unclear. In this study, based on molecular dynamics simulations, we propose a p(f)-matrix method, in which p(f) is decomposed into contributions from each local region of the channel. Diagonal elements of the p(f) matrix are equivalent to the local permeability at each region of the channel, and off-diagonal elements represent correlated motions of water molecules in different regions. Averaging both diagonal and off-diagonal elements of the p(f) matrix recovers p(f) for the entire channel; this implies that correlated motions between distantly-separated water molecules, as well as adjacent water molecules, influence the osmotic permeability. The p(f) matrices from molecular dynamics simulations of five aquaporins (AQP0, AQP1, AQP4, AqpZ, and GlpF) indicated that the reduction in the water correlation across the Asn-Pro-Ala region, and the small local permeability around the ar/R region, characterize the transport efficiency of water. These structural determinants in water permeation were confirmed in molecular dynamics simulations of three mutants of AqpZ, which mimic AQP1.     

Snazer: the simulations and networks analyzer

Background  

Networks are widely recognized as key determinants of structure and function in systems that span the biological, physical, and social sciences. They are static pictures of the interactions among the components of complex systems. Often, much effort is required to identify networks as part of particular patterns as well as to visualize and interpret them.     

Temperature dependence of the generalized frequency distribution of water molecules: comparison of experiments and molecular dynamics simulations

The results of inelastic neutron scattering experiments on water in the temperature interval 300–623 K along the coexistence curve are compared with data obtained from molecular dynamics simulations. In general, a good agreement between experiments and calculations is observed and it serves as a satisfactory test of the potential models employed. The temperature dependence of the generalized frequency distribution of water molecules obtained by both experiment and computer simulation demonstrates the accordance with the temperature evolution of the water structure, extracted from neutron and X-ray diffraction measurements.     

Structured pathway across the transition state for peptide folding revealed by molecular dynamics simulations

Small globular proteins and peptides commonly exhibit two-state folding kinetics in which the rate limiting step of folding is the surmounting of a single free energy barrier at the transition state (TS) separating the folded and the unfolded states. An intriguing question is whether the polypeptide chain reaches, and leaves, the TS by completely random fluctuations, or whether there is a directed, stepwise process. Here, the folding TS of a 15-residue β-hairpin peptide, Peptide 1, is characterized using independent 2.5 μs-long unbiased atomistic molecular dynamics (MD) simulations (a total of 15 μs). The trajectories were started from fully unfolded structures. Multiple (spontaneous) folding events to the NMR-derived conformation are observed, allowing both structural and dynamical characterization of the folding TS. A common loop-like topology is observed in all the TS structures with native end-to-end and turn contacts, while the central segments of the strands are not in contact. Non-native sidechain contacts are present in the TS between the only tryptophan (W11) and the turn region (P7-G9). Prior to the TS the turn is found to be already locked by the W11 sidechain, while the ends are apart. Once the ends have also come into contact, the TS is reached. Finally, along the reactive folding paths the cooperative loss of the W11 non-native contacts and the formation of the central inter-strand native contacts lead to the peptide rapidly proceeding from the TS to the native state. The present results indicate a directed stepwise process to folding the peptide.     

Reconstruction of the src-SH3 protein domain transition state ensemble using multiscale molecular dynamics simulations

We use an integrated computational approach to reconstruct accurately the transition state ensemble (TSE) for folding of the src-SH3 protein domain. We first identify putative TSE conformations from free energy surfaces generated by importance sampling molecular dynamics for a fully atomic, solvated model of the src-SH3 protein domain. These putative TSE conformations are then subjected to a folding analysis using a coarse-grained representation of the protein and rapid discrete molecular dynamics simulations. Those conformations that fold to the native conformation with a probability (P(fold)) of approximately 0.5, constitute the true transition state. Approximately 20% of the putative TSE structures were found to have a P(fold) near 0.5, indicating that, although correct TSE conformations are populated at the free energy barrier, there is a critical need to refine this ensemble. Our simulations indicate that the true TSE conformations are compact, with a well-defined central beta sheet, in good agreement with previous experimental and theoretical studies. A structured central beta sheet was found to be present in a number of pre-TSE conformations, however, indicating that this element, although required in the transition state, does not define it uniquely. An additional tight cluster of contacts between highly conserved residues belonging to the diverging turn and second beta-sheet of the protein emerged as being critical elements of the folding nucleus. A number of commonly used order parameters to identify the transition state for folding were investigated, with the number of native Cbeta contacts displaying the most satisfactory correlation with P(fold) values.     

Theoretical electrical conductivity of hydrogen-bonded benzamide-derived molecules and single DNA bases

A benzamide molecule is used as a “reader” molecule to form hydrogen bonds with five single DNA bases, i.e., four normal single DNA bases A,T,C,G and one for 5methylC. The whole molecule is then attached to the gold surface so that a meta-molecule junction is formed. We calculate the transmission function and conductance for the five metal–molecule systems, with the implementation of density functional theory-based non-equilibrium Green function method. Our results show that each DNA base exhibits a unique conductance and most of them are on the pS level. The distinguishable conductance of each DNA base provides a way for the fast sequencing of DNA. We also investigate the dependence of conductivity of such a metal–molecule system on the hydrogen bond length between the “reader” molecule and DNA base, which shows that conductance follows an exponential decay as the hydrogen bond length increases, i.e., the conductivity is highly sensitive to the change in hydrogen bond length.