Structural stability of myoglobin and glycomyoglobin: a comparative molecular dynamics simulation study

Glycoproteins are formed as the result of enzymatic glycosylation or chemical glycation in the body, and produced in vitro in industrial processes. The covalently attached carbohydrate molecule(s) confer new properties to the protein, including modified stability. In the present study, the structural stability of a glycoprotein form of myoglobin, bearing a glucose unit in the N-terminus, has been compared with its native form by the use of molecular dynamics simulation. Both structures were subjected to temperatures of 300 and 500 K in an aqueous environment for 10 ns. Changes in secondary structures and RMSD were then assessed. An overall higher stability was detected for glycomyoglobin, for which the most stable segments/residues were highlighted and compared with the native form. The simple addition of a covalently bound glucose is suggested to exert its stabilizing effect via increased contacts with surrounding water molecules, as well as a different pattern of interactions with neighbor residues.

Electronic supplementary material

The online version of this article (doi:10.1007/s10867-015-9383-2) contains supplementary material, which is available to authorized users.     

NO rebinding to myoglobin: a reactive molecular dynamics study

The rebinding of NO to myoglobin after photolysis is studied using the 'reactive molecular dynamics' method. In this approach the energy of the system is evaluated on two potential energy surfaces that include the heme-ligand interactions which change between liganded and unliganded myoglobin. This makes it possible to take into account in a simple way, the high dimensionality of the transition seam connecting the reactant and product states. The dynamics of the dissociated NO molecules are examined, and the geometrical and energetic properties of the transition seam are studied. Analysis of the frequency of recrossing shows that the height of the effective rebinding barrier is dependent on the time after photodissociation. This effect is due mainly to protein relaxation and may contribute to the experimentally observed non-exponential rebinding rate of NO, as has been suggested previously.     

Cholesterol-sphingomyelin interactions: a molecular dynamics simulation study

Stearoylsphingomyelin (SSM) bilayers containing 0, 22, and 50 mol % cholesterol (Chol) and a pentadecanoyl-stearoylphosphatidylcholine (15SPC) bilayer containing 22 mol % Chol were molecular dynamics simulated at two temperatures (37 degrees C and 60 degrees C). 15SPC is the best PC equivalent of SSM. The Chol effect on the SSM bilayer differs significantly from that on the 15SPC bilayer. At the same temperature and Chol content, H-bonding of Chol with SSM is more extensive than with 15SPC. SSM-Chol H-bonding anchors the OH group of Chol in the lower regions of the SSM-Chol bilayer interface. Such a location strengthens the influence of Chol on the SSM chains. In effect, the phase of the SSM-Chol bilayer containing 22 mol % Chol at 37 degrees C is shifted from the gel to the liquid-ordered phase, and the bilayer displays similar properties below and above the main phase-transition temperature for a pure SSM bilayer of approximately 45 degrees C. In contrast, due to a higher location, Chol is not able to change the phase of the 15SPC-Chol bilayer, which at 37 degrees C remains in the gel phase. Chol affects both the core and interface of the SSM bilayer. With increasing Chol content, the order of SSM chains and hydration of SSM headgroups increase, whereas polar interactions between lipids decrease.     

Conformational dynamics of RNA-peptide binding: a molecular dynamics simulation study

Molecular dynamics simulations of the binding of the heterochiral tripeptide KkN to the transactivation responsive (TAR) RNA of HIV-1 is presented, using an all-atom force field with explicit water. To obtain starting structures for the TAR-KkN complex, semirigid docking calculations were performed that employ an NMR structure of free TAR RNA. The molecular dynamics simulations show that the starting structures in which KkN binds to the major groove of TAR (as it is the case for the Tat-TAR complex of HIV-1) are unstable. On the other hand, the minor-groove starting structures are found to lead to several binding modes, which are stabilized by a complex interplay of stacking, hydrogen bonding, and electrostatic interactions. Although the ligand does not occupy the binding position of Tat protein, it is shown to hinder the interhelical motion of free TAR RNA. The latter is presumably necessary to achieve the conformational change of TAR RNA to bind Tat protein. Considering the time evolution of the trajectories, the binding process is found to be ligand-induced and cooperative. That is, the conformational rearrangement only occurs in the presence of the ligand and the concerted motion of the ligand and a large part of the RNA binding site is necessary to achieve the final low-energy binding state.     

Copper-containing nitrite reductase: a DFT study of nitrite and nitric oxide adducts

Copper-containing nitrite reductases (Cu-NIRs) reduce nitrite to NO. Reported here are DFT (density functional theory) results on models of the Cu-NIR active site bound to nitrite and nitric oxide. The Cu-NIR active site appears to have been designed to exclude N-nitrite binding even though N-O bond cleavage would be equally facile in the N- and O-isomers. The active site also appears to force a side-on coordination of the end-product, nitric oxide. The latter feature has to rely on the sterics of the active site to destabilize, thermodynamically speaking, the Cu-NO adduct; under these conditions, the absence of N-nitrite coordination is proposed to be merely a side-effect. For the Cu(II)-NO adduct, sterical crowding appears to also favour the Cu-NO electromer over Cu(I)-NO+, helping to avoid the potentially damaging chemistry associated with an NO+ moiety. These conclusions are in reasonable agreement with previous conclusions drawn from experiment [Science 304 (2004) 867].     

Prevalent mutations of human prion protein: a molecular modeling and molecular dynamics study

Point mutations in the human prion protein gene, leading to amino acid substitutions in the human prion protein contribute to conversion of PrPC to PrPSc and amyloid formation, resulting in prion diseases such as familial Creutzfeldt-Jakob disease (CJD), Gerstmann-Straussler-Scheinker disease (GSS), and fatal familial insomnia. We have investigated impressions of prevalent mutations including Q217R, D202N, F198S, on the human prion protein and compared the mutant models with wild types. Structural analyses of models were performed with molecular modeling and molecular dynamics simulation methods. According to our results, frequently occurred mutations are observed in conserved and fully conserved sequences of human prion protein and the most fluctuation values occur in the Helix 1 around residues 144-152 and C-terminal end of the Helix 2. Our analysis of results obtained from MD simulation clearly shows that this long-range effect plays an important role in the conformational fluctuations in mutant structures of human prion protein. Results obtained from molecular modeling such as creation or elimination of some hydrogen bonds, increase or decrease of the accessible surface area and molecular surface, loss or accumulation of negative or positive charges on specific positions, and altering the polarity and pKa values, show that amino acid point mutations, though not urgently change the stability of PrP, might have some local impacts on the protein interactions which are required for oligomerization into fibrillar species.     

Pressure denaturation of apomyoglobin: a molecular dynamics simulation study

The effect of pressure on the structure and mobility of Sperm Wale Apomyoglobin was studied by Molecular Dynamics computer simulation at 1 bar and 3 kbar (1 atm=1.01325 bar=101.325 kPa). The results are in good agreement with the available experimental data, allowing further analysis of other features of the effect of pressure on the protein solution. From the analysis of Secondary Structures (SS) along the trajectories it is observed that alpha-helixes are favoured under pressure at the expense of bends, turns and 3-helixes. The studies of mobility show that although the general mobility is restricted under pressure this is not true for some particular residues. The studies of tertiary structure show important conformational changes. The evolution of the Solvent Accessed Surface (SAS) with pressure shows a notorious increase due almost completely to a biased raise in the hydrophobic area exposed, which consequently shows that the hydrophobic interaction is considerably weaker under high hydrostatic pressure conditions.     

Combination of anti-hypertensive drugs: a molecular dynamics simulation study

The anti-hypertensive drugs amlodipine, atenolol and lisinopril, in ordinary and PEGylated forms, with different combined-ratios, were studied by molecular dynamics simulations using GROMACS software. Twenty simulation systems were designed to evaluate the interactions of drug mixtures with a dimyristoylphosphatidylcholine (DMPC) lipid bilayer membrane, in the presence of water molecules. In the course of simulations, various properties of the systems were investigated, including drug location, diffusion and mass distribution in the membrane; drug orientation; the lipid chain disorder as a result of drug penetration into the DMPC membrane; the number of hydrogen bonds; and drug surface area. According to the results obtained, combined drugs penetrate deeper into the DMPC lipid bilayer membrane, and the lipid chains remain ordered. Also, the combined PEGylated drugs, at a combination ratio of 1:1:1, enhance drug penetration into the DMPC membrane, reduce drug agglomeration, orient the drug in a proper angle for easy penetration into the membrane, and decrease undesirable lipotoxicity due to distorted membrane self-assembly and thickness.

Open image in new window

Graphical abstract ?

     

Photoinduced conformational dynamics of a photoswitchable peptide: a nonequilibrium molecular dynamics simulation study

Employing nonequilibrium molecular dynamics simulations, a comprehensive computational study of the photoinduced conformational dynamics of a photoswitchable bicyclic azobenzene octapeptide is presented. The calculation of time-dependent probability distributions along various global and local reaction coordinates reveals that the conformational rearrangement of the peptide is rather complex and occurs on at least four timescales: 1) After photoexcitation, the azobenzene unit of the molecule undergoes nonadiabatic photoisomerization within 0.2 ps. 2) On the picosecond timescale, the cooling (13 ps) and the stretching (14 ps) of the photoexcited peptide is observed. 3) Most reaction coordinates exhibit a 50-100 ps component reflecting a fast conformational rearrangement. 4) The 500-1000 ps component observed in the simulation accounts for the slow diffusion-controlled conformational equilibration of the system. The simulation of the photoinduced molecular processes is in remarkable agreement with time-resolved optical and infrared experiments, although the calculated cooling as well as the initial conformational rearrangements of the peptide appear to be somewhat too slow. Based on an ab initio parameterized vibrational Hamiltonian, the time-dependent amide I frequency shift is calculated. Both intramolecular and solvent-induced contributions to the frequency shift were found to change by < or = 2 cm(-1), in reasonable agreement with experiment. The potential of transient infrared spectra to characterize the conformational dynamics of peptides is discussed in some detail.     

Retraction: Molecular modeling and molecular dynamics simulation studies of Delta-Notch complex

Majumder R, Roy, S, Thakur AR, J Biomol Struct Dyn .2011 Oct;29(2):297-310. According to a letter that the Journal received as an attachment to an email dated September 27, 2011 from the senior author Prof. Ashoke Ranjan Thakur, they have indicated to have their paper retracted, because an original reviewer raised questions (when he saw the published copy) regarding the correctness assignment of a residue number in their calculation. The authors have already started the recalculation. They intend to submit the whole paper back again for consideration for publication in a future issue of the Journal.     

Interfacial water on crystalline silica: a comparative molecular dynamics simulation study

Understanding the properties of interfacial water at solid–liquid interfaces is important in a wide range of applications. Molecular dynamics is becoming a widespread tool for this purpose. Unfortunately, however, the results of such studies are known to strongly depend on the selection of force fields. It is, therefore, of interest to assess the extent by which the implemented force fields can affect the predicted properties of interfacial water. Two silica surfaces, with low and high surface hydroxyl density, respectively, were simulated implementing four force fields. These force fields yield different orientation and flexibility of surface hydrogen atoms, and also different interaction potentials with water molecules. The properties for interfacial water were quantified by calculating contact angles, atomic density profiles, surface density distributions, hydrogen bond density profiles and residence times for water near the solid substrates. We found that at low surface density of hydroxyl groups, the force field strongly affects the predicted contact angle, while at high density of hydroxyl groups, water wets all surfaces considered. From a molecular-level point of view, our results show that the position and intensity of peaks observed from oxygen and hydrogen atomic density profiles are quite different when different force fields are implemented, even when the simulated contact angles are similar. Particularly, the surfaces simulated by the CLAYFF force field appear to attract water more strongly than those simulated by the Bródka and Zerda force field. It was found that the surface density distributions for water strongly depend on the orientation of surface hydrogen atoms. In all cases, we found an elevated number of hydrogen bonds formed between interfacial water molecules. The hydrogen bond density profile does not depend strongly on the force field implemented to simulate the substrate, suggesting that interfacial water assumes the necessary orientation to maximise the number of water–water hydrogen bonds irrespectively of surface properties. Conversely, the residence time for water molecules near the interface strongly depends on the force field and on the flexibility of surface hydroxyl groups. Specifically, water molecules reside for longer times at contact with rigid substrates with high density of hydroxyl groups. These results should be considered when comparisons between simulated and experimental data are attempted.     

Cholesterol effects on a mixed-chain phosphatidylcholine bilayer: a molecular dynamics simulation study

A molecular dynamics simulation of a mono-cis-unsaturated 1-palmitoyl-2-oleoyl-phosphatidylcholine bilayer containing approximately 22 mol% of cholesterol (POPC-Chol) was carried out for 15 ns. An 8-ns trajectory was analysed to determine the effects of Chol on the membrane properties and compare it with that on the fully saturated 1,2-dimyristoyl-phosphatidylcholine bilayer containing approximately 22 mol% of Chol (DMPC-Chol). The study suggests that the experimentally observed weaker effect of Chol on the POPC than DMPC bilayer might result from a different vertical localisation of the Chol hydroxyl group (OH-Chol) in both bilayers: in the POPC-Chol bilayer, OH-Chol is placed approximately 3 A higher in the bilayer interface than in the DMPC-Chol bilayer. Because of the rigid cis double bond in the beta-chain of POPC, Chol fits worse to the POPC-Chol membrane environment and is pushed up, in effect all Chol ring atoms are, on average, located above the double bond. Both in mono-cis-unsaturated and fully saturated PC bilayers, Chol induces stronger van der Waals interactions among the chains, whereas its interactions with the chains are weak. In contrast to DMPC, the smooth alpha-face of the Chol ring lowers the order of POPC chains, whereas the rough beta-face increases the order.     

Homology modeling and molecular dynamics study on N-acetylneuraminate lyase

With homology modeling techniques, molecular mechanics and molecular dynamics methods, a 3D structure model of N-acetylneuraminate lyase from human (hNAL, EC 4.1.3.3) was created and refined. This model was further assessed by Profile-3D and PROCHECK, which confirms that the refined model is reliable. Furthermore, the docking results of the substrates (sialic acid and KDO) into the active site of hNAL indicate that hNAL can cleave the sialic acid and KDO. Thr51 and Tyr143 may be the key amino acids residues as they have strong hydrogen bonding interactions with the substrates, which is in good agreement with the experimental results by Izard et al. (Structure 2:361–369. doi:10.1016/S0969-2126(00)00038-1 (1994)). From the docking studies, we also suggest that Asp176 and Ser218 only form hydrogen bonds with sialic acid, therefore, they may help sialic acid interact with hNAL steadly.     

Molecular modeling and molecular dynamics simulation studies of Delta-Notch complex

Notch is a single-pass transmembrane receptor protein which is composed of a short extracellular region, a single-pass transmembrane domain and a small intracellular region. Notch ligand like Delta, member of the DSL protein family, is also single-pass transmembrane protein. It has been demonstrated that of the 36 EGF repeats of Notch, 11th and 12th are sufficient to mediate interactions with Delta. Crystal structure of mammalian Notch extracellular ligand binding domain contains 11 and 12 EGF-like repeats. Here a portion of the Delta protein of Drosophila, known to interact with Notch extracellular domain (ECD) has been modeled using homology modeling. The structure of the Delta-Notch complex was subsequently modeled by protein docking method using GRAMM. MD simulations of the modeled structures were performed. The structure for Delta-Notch complex has been proposed based on interaction energy parameter and planarity studies.     

Pharmacophore modeling, molecular docking, and molecular dynamics simulation approaches for identifying new lead compounds for inhibiting aldose reductase 2

Aldose reductase 2 (ALR2), which catalyzes the reduction of glucose to sorbitol using NADP as a cofactor, has been implicated in the etiology of secondary complications of diabetes. A pharmacophore model, Hypo1, was built based on 26 compounds with known ALR2-inhibiting activity values. Hypo1 contains important chemical features required for an ALR2 inhibitor, and demonstrates good predictive ability by having a high correlation coefficient (0.95) as well as the highest cost difference (128.44) and the lowest RMS deviation (1.02) among the ten pharmacophore models examined. Hypo1 was further validated by Fisher's randomization method (95%), test set (r = 0.91), and the decoy set shows the goodness of fit (0.70). Furthermore, during virtual screening, Hypo1 was used as a 3D query to screen the NCI database, and the hit leads were sorted by applying Lipinski's rule of five and ADME properties. The best-fitting leads were subjected to docking to identify a suitable orientation at the ALR2 active site. The molecule that showed the strongest interactions with the critical amino acids was used in molecular dynamics simulations to calculate its binding affinity to the candidate molecules. Thus, Hypo1 describes the key structure-activity relationship along with the estimated activities of ALR2 inhibitors. The hit molecules were searched against PubChem to find similar molecules with new scaffolds. Finally, four molecules were found to satisfy all of the chemical features and the geometric constraints of Hypo1, as well as to show good dock scores, PLPs and PMFs. Thus, we believe that Hypo1 facilitates the selection of novel scaffolds for ALR2, allowing new classes of ALR2 inhibitors to be designed.     

Structural modeling and molecular dynamics simulation of the actin filament

Actin is a major structural protein of the eukaryotic cytoskeleton and enables cell motility. Here, we present a model of the actin filament (F-actin) that not only incorporates the global structure of the recently published model by Oda et al. but also conserves internal stereochemistry. A comparison is made using molecular dynamics simulation of the model with other recent F-actin models. A number of structural determents such as the protomer propeller angle, the number of hydrogen bonds, and the structural variation among the protomers are analyzed. The MD comparison is found to reflect the evolution in quality of actin models over the last 6 years. In addition, simulations of the model are carried out in states with both ADP or ATP bound and local hydrogen-bonding differences characterized.     

Principles of carbopeptoid folding: a molecular dynamics simulation study.

The conformational spaces of five oligomers of tetrahydrofuran-based carbopeptoids in chloroform and dimethyl sulfoxide were investigated through nine molecular dynamics simulations. Prompted by nuclear magnetic resonance experiments that indicated various stable folds for some but not all of these carbopeptoids, their folding behaviour was investigated as a function of stereochemistry, chain length and solvent. The conformational distributions of these molecules were analysed in terms of occurrence of hydrogen bonds, backbone torsional-angle distributions, conformational clustering and solute configurational entropy. While a cis-linkage across the tetrahydrofuran ring favours right-handed helical structures, a trans-linkage results in a larger conformational variability. Intra-solute hydrogen bonding is reduced with increasing chain length and with increasing solvent polarity. Solute configurational entropies confirm the picture obtained: they are smaller for cis- than for trans-linked peptides, for chloroform than for dimethyl sulfoxide as solvent and for shorter peptide chains. The simulations provide an atomic picture of molecular conformational variability that is consistent with the available experimental data.     

High-Temperature unfolding of a trp-Cage mini-protein: a molecular dynamics simulation study

Background  

Trp cage is a recently-constructed fast-folding miniprotein. It consists of a short helix, a 3,10 helix and a C-terminal poly-proline that packs against a Trp in the alpha helix. It is known to fold within 4 ns.     

Extended molecular dynamics simulation of the carbon monoxide migration in sperm whale myoglobin

We report the results of an extended molecular dynamics simulation on the migration of photodissociated carbon monoxide in wild-type sperm whale myoglobin. Our results allow following one possible ligand migration dynamics from the distal pocket to the Xe1 cavity via a path involving the other xenon binding cavities and momentarily two additional packing defects along the pathway. Comparison with recent time resolved structural data obtained by Laue crystallography with subnanosecond to millisecond resolution shows a more than satisfactory agreement. In fact, according to time resolved crystallography, CO, after photolysis, can occupy the Xe1 and Xe4 cavities. However, no information on the trajectory of the ligand from the distal pocket to the Xe1 is available. Our results clearly show one possible path within the protein. In addition, although our data refer to a single trajectory, the local dynamics of the ligand in each cavity is sufficiently equilibrated to obtain local structural and thermodynamic information not accessible to crystallography. In particular, we show that the CO motion and the protein fluctuations are strictly correlated: free energy calculations of the migration between adjacent cavities show that the migration is not a simple diffusion but is kinetically or thermodynamically driven by the collective motions of the protein; conversely, the protein fluctuations are influenced by the ligand in such a way that the opening/closure of the passage between adjacent cavities is strictly correlated to the presence of CO in its proximity. The compatibility between time resolved crystallographic experiments and molecular dynamics simulations paves the way to a deeper understanding of the role of internal dynamics and packing defects in the control of ligand binding in heme proteins.     

Non-polar interactions between cholesterol and phospholipids: a molecular dynamics simulation study

A 15-ns molecular dynamics simulation of the fully hydrated dimyristoylphosphatidylcholine-cholesterol (DMPC-Chol) bilayer containing approximately 22 mol% Chol was carried out. An 8-ns trajectory was analysed to investigate the effect of Chol on the chain packing in the bilayer core. While the packing of DMPC chains on the smooth alpha-face side of the Chol ring is similar to that in the pure DMPC bilayer, the packing on the rough beta-face side is less regular and less tight. Two methyl groups located on the Chol beta-face disturb the packing; in effect, van der Waals (vdW) interactions between Chol rings and DMPC chains are weaker than the ones between sole DMPC chains. VdW interactions between an alkyl chain of DMPC and an isooctyl tail of Chol are similarly strong as those between two DMPC chains.