The dynamic behaviours of protein BMP-2 on hydroxyapatite nanoparticles

The interaction between proteins and nanoscale inorganic particles is one of the most important topics in many fields. In this study, the dynamic behaviours of protein bone morphogenetic protein-2 (BMP-2) (with six different orientations) on hydroxyapatite (HAP) (001) surface were studied using the molecular dynamics and steered molecular dynamics simulation. The results show that the orientation of protein BMP-2 has obvious influence on its adsorption–desorption behaviours. Among the six systems studied in this article, system I exhibits the strongest interaction with the HAP (001) surface, and the number of the adsorbed residues is more than any one of the other five systems correspondingly. These findings suggest that there will be a preferential orientation when a protein is adsorbed onto a nanoscale interface. For protein BMP-2 interacting with the HAP (001) surface, the preferential orientation is the orientation in system I.     

Modeling the interaction of seven bisphosphonates with the hydroxyapatite(100) face

The interaction of seven pamidronate bisphosphonate (Pami-BPs) and its analogs with the hydroxyapatite (HAP) (100) surface was studied using density functional theory (DFT) and molecular dynamic (MD) methods. Partial Mulliken oxygen atomic charges in protonated structures were calculated at the level of B3LYP/6-31G*. The MD simulation was performed using the Discover module of Material Studio by compass force field. The results indicate the abilities of donating electrons of the oxygen atoms of the phosphate groups that are closely associated with the antiresorptive potency. The binding energies, including vdw and electrostatic, are used to discuss the mechanism of antiresorption. The results of calculations show that the strength of interaction of the HAP (100) face with the bisphosphonates is N(4)?>?N(6)?>?N(7)?>?N(5)?>?N(3)?>?N(2)?>?N(1) according to their experimental pIC(50) values.     

Computational analysis on conformational dynamics of bone morphogenetic protein-2 (BMP-2)

BMP-2 is widely used for bone regeneration because of its ability to induce osteoblast differentiation and proliferation. The pharmaceutical application of BMP-2 as bone implant makes the studies on stability and conformational dynamics very relevant as proteins are functional only in their native three-dimensional state. Knowing the factors affecting BMP-2 structure becomes essential for designing bone implants activated by BMP-2. In order to explore the influence of temperature and hydration on protein conformation, we have performed the molecular dynamics (MD) simulations at the time scale of 100 ns with two different force fields. We have examined the dynamic behaviour of BMP-2 monomer and dimer in aqueous medium as well as in vacuum at four different temperatures (300, 350, 400 and 450 K). MD simulation of BMP-2 monomer and dimer in water and vacuum environments shows the major contribution of water in structure stabilization. Temperature of the system affects the secondary structure differently in case of monomer and dimer simulation and the dynamics also depends on the environment viz. vacuum and aqueous. Vacuum simulations show very early loss of the major secondary structure content. On the other hand, BMP-2 monomer and dimer in aqueous environment show the unfolding of α-helix with increasing temperature. This unfolded α-helix is converted into β-sheet at 400 K in monomer of BMP-2. Contrary to this, we did not observe β-sheet formation in dimer BMP-2 even at 450 K indicating that monomers are more aggregation prone entity as compared to dimers of BMP-2.     

Red fluorescent protein DsRed: Parametrization of its chromophore as an amino acid residue for computer modeling in the OPLS-AA force field

Topology of the neutral form of the DsRed fluorescent protein chromophore as a residue of [(4-cis)-2-[(1-cis)-4-amino-4-oxobutanimidoyl]-4-(4-hydroxybenzylidene)-5-oxo-4,5-dihydro-1H-imidazol-1-yl]acetic acid was calculated with OPLS-AA force field. Use of this topology and molecular dynamics simulation allows calculating the parameters of proteins that contain such residue in their polypeptide chains. The chromophore parameters were obtained by ab initio (RHF/6-31G**) quantum chemical calculations applying density functional theory (B3LYP). Using this chromophore, we have calculated the molecular dynamics trajectory of tetrameric fluorescent protein DsRed in solution at 300 K (4 nsec). Correctness of the chromophore parametrization was revealed by comparison of quantitative characteristics of the chromophore structure obtained from the molecular dynamic simulations of DsRed protein with the quantitative characteristics of the chromophore based on the crystallographic X-ray data of fluorescent protein DsRed (PDB ID: 1ZGO, 1G7K, and 1GGX), and also with the quantitative characteristics of the chromophore obtained by quantum chemical calculations. Inclusion of the neutral form of DsRed protein chromophore topology into the OPLS-AA force field yielded the extended force field OPLS-AA/DsRed. This force field can be used for molecular dynamics calculations of proteins containing the DsRed chromophore. The parameter set presented in this study can be applied for similar extension in any other force fields.     

Serine proteases: an ab initio molecular dynamics study

In serine proteases (SPs), the H-bond between His57 and Asp102 and that between Gly193 and the transition state intermediate play a crucial role in enzymatic function. To shed light on the nature of these interactions, we have carried out ab initio molecular dynamics simulations on complexes representing adducts between the reaction intermediate and elastase (one protein belonging to the SP family). Our calculations indicate the presence of a low-barrier H-bond between His57 and Asp102, in complete agreement with NMR experiments on enzyme-transition state analogue complexes. Comparison with an ab initio molecular dynamics simulation on a model of the substrate-enzyme adduct indicates that the Gly193-induced strong stabilization of the intermediate is accomplished by charge/dipole interactions and not by H-bonding as previously suggested. Inclusion of the protein electric field in the calculations does not affect significantly the charge distribution.     

Understanding interaction and dynamics of water molecules in the epoxy via molecular dynamics simulation

The robust structural integrity of the epoxy plays an important role in ensuring the long-term service life of its applications, which is affected by the absorbed moisture. In order to understand the mechanism of the moisture effect, the knowledge of the interaction and dynamics of the water molecules inside the epoxy is of great interest. Molecular dynamics simulation is used in this work to investigate the structure and bonding behaviour of the water molecules in the highly cross-linked epoxy network. When the moisture concentration is low, the water molecules are well dispersed in the cross-linked structure and located in the vicinity of the epoxy functional groups, which predominantly form the hydrogen bond (H-bond) with the epoxy network, resulting in the low water mobility in the epoxy. At the high concentration, the water favourably forms the large cluster due to the predominant water–water H-bond interaction, and the water molecules diffuse primarily inside the cluster, which leads to the high water mobility and the accelerated H-bond dynamics. The variation of the bonding behaviour and dynamics of the water molecules reported here could be exploited to understand the material change and predict the long-term performance of the epoxy-based products during the intended service life.     

Long-term molecular dynamics simulation of copper plastocyanin in water

A long molecular dynamics simulation (1.1 ns) of fully hydrated plastocyanin has been performed and analysed to relate protein dynamics to structural elements and functional properties. The solvated structure is described in detail by the analysis of H-bond network. During all the simulation, the crystal H-bond network is maintained in the beta-sheet regions, while several H-bonds are broken or formed on the external surface of the protein. To evaluate whether such changes could be due to conformational rearrangements or to solvent competition, we have examined the average number of H-bonds between protein atoms and water molecules, and the root mean square deviations from crystal structure as a function of protein residues. Protein mobility and flexibility have been examined by positional and dihedral angle rms fluctuations. Finally, cross-correlation maps have revealed the existence of correlated motions among residues connected by hydrogen bonds.     

Multiple pH regime molecular dynamics simulation for pK calculations

Ionisation equilibria in proteins are influenced by conformational flexibility, which can in principle be accounted for by molecular dynamics simulation. One problem in this method is the bias arising from the fixed protonation state during the simulation. Its effect is mostly exhibited when the ionisation behaviour of the titratable groups is extrapolated to pH regions where the predetermined protonation state of the protein may not be statistically relevant, leading to conformational sampling that is not representative of the true state. In this work we consider a simple approach which can essentially reduce this problem. Three molecular dynamics structure sets are generated, each with a different protonation state of the protein molecule expected to be relevant at three pH regions, and pK calculations from the three sets are combined to predict pK over the entire pH range of interest. This multiple pH molecular dynamics approach was tested on the GCN4 leucine zipper, a protein for which a full data set of experimental data is available. The pK values were predicted with a mean deviation from the experimental data of 0.29 pH units, and with a precision of 0.13 pH units, evaluated on the basis of equivalent sites in the dimeric GCN4 leucine zipper.     

Parameterization of aromatic azido groups: application as photoaffinity probes in molecular dynamics studies

The accuracy of molecular dynamics (MD) simulations is limited by the availability of parameters for the molecular system of interest. In most force fields, parameters of common chemical groups are already present. With the development of novel small organic molecules as probes to study biological systems, more chemical groups require parameterization. An azide group is often used in studies of biological systems but computational studies are still impeded by the lack of parameters. In this paper, we present a set of molecular mechanics (MM) parameters for aromatic and aliphatic azido groups, and their application in MD simulations of a photoaffinity probe currently used in our laboratory for mapping binding modes available in the active site of histone deacetylases. The parameters were developed for the generalized Amber force field (GAFF) using density functional theory (DFT) calculations at B3LYP 6-311G(d) level. The parameters were validated by geometry optimization and MD simulations.     

Exploration of the binding of proton pump inhibitors to human P450 2C9 based on docking and molecular dynamics simulation

Human P450 protein CYP2C9 is one of the major drug-metabolizing isomers, contributing to the oxidation of 16% of the drugs currently in clinical use. To examine the interaction mechanisms between CYP2C9 and proton pump inhibitions (PPIs), we used molecular docking and molecular dynamics (MD) simulation methods to investigate the conformations and interactions around the binding sites of PPIs/CYPP2C9. Results from molecular docking and MD simulations demonstrate that nine PPIs adopt two different conformations (extended and U-bend structures) at the binding sites and position themselves far above the heme of 2C9. The presence of PPIs changes the secondary structures and residue flexibilities of 2C9. Interestingly, at the binding sites of all PPI–CYP2C9 complexes except for Lan/CYP2C9, there are hydrogen-bonding networks made of PPIs, water molecules, and some residues of 2C9. Moreover, there are strong hydrophobic interactions at all binding sites for PPIs/2C9, which indicate that electrostatic interactions and hydrophobic interactions appear to be important for stabilizing the binding sites of most PPIs/2C9. However, in the case of Lan/2C9, the hydrophobic interactions are more important than the electrostatic interactions for stabilizing the binding site. In addition, an interesting conformational conversion from extended to U-bend structures was observed for pantoprazole, which is attributed to an H-bond interaction in the binding pocket, an internal π–π stacking interaction, and an internal electrostatic interaction of pantoprazole.     

Structure determination of a pseudotripeptide zinc complex with the COSMOS-NMR force field and DFT methods

A His-X-His pseudotripeptide zinc complex (X is a N-alkyl glycine derivative) similar to the catalytic center of the carbonic anhydrase was computer designed and experimentally synthesized. Using 2D-NMR techniques, all proton, carbon chemical shifts and nuclear overhauser effect signals were assigned. The three-dimensional structure of the complex was determined with the COSMOS (computer simulation of molecular structures) force field by applying ¹³C bond polarization theory chemical shift pseudo forces and restrictions for NOE distances. From molecular dynamics, simulated annealing simulations and geometry optimizations, the three best force field structures were taken for a final investigation by density functional theory calculations.     

Mechanism of titin unfolding by force: insight from quasi-equilibrium molecular dynamics calculations

We have studied the unfolding by force of one of the immunoglobulin domains of the muscle protein titin using molecular dynamics simulations at 300 K. Previous studies, done at constant pulling rates, showed that under the effect of the force two strands connected to each other by six backbone H-bonds are pulled apart. No details about the mechanism of H-bond breaking were provided. Our simulation protocol "pull and wait" was designed to correspond to very slow pulling, more similar to the rates used in experiments than are the protocols used in previous computational studies. Under these conditions interstrand backbone H-bonds are not "ripped apart" by the application of the force. Instead, small elongations produced by the force weaken specific backbone H-bonds with respect to water-backbone H-bonds. These weakened bonds allow a single water molecule to make H-bonds to the CO and the NH of the same backbone H-bond while they are still bound to each other. The backbone H-bond then breaks (distance > 3.6 A), but its donor and acceptor atoms remain bound to the same water molecule. Further separation of the chains takes place when a second water molecule makes an H-bond with either the protein backbone donor or acceptor atom. Thus, the force does not directly break the main chain H-bonds: it destabilizes them in such a way that they are replaced by H-bonds to water. With this mechanism, the force necessary to break all the H-bonds required to separate the two strands will be strongly dependent on the pulling speed. Further simulations carried out at low forces but long waiting times (> or = 500 ps, < or = 10 ns) show that, given enough time, even a very small pulling force (< 400 pN) is sufficient to destabilize the interstrand H-bonds and allow them to be replaced by H-bonds to two water molecules. As expected, increasing the temperature to 350 K allows the interstrand H-bonds to break at lower forces than those required at 300 K.     

Molecular dynamics and quantum chemical studies of solvent effects on cyclo glycylglycine and glycylalanine dipeptides

Hydrogen bond (H-bond) interactions between the two cyclo dipeptides, cyclo(glycyl-glycine) (CGG) and cyclo(glycyl-alanine) (CGA), and water have been studied using molecular dynamics (MD) and quantum chemical methods. The MD studies have been carried out on CGG and CGA in water using fixed charge force field (AMBER ff03) for over 10 ns with a MD time step of 2 fs. The results of this study show that the solvation pattern influences the conformations of the cyclo dipeptides. Following molecular simulations, post Hartree–Fock and density functional theory methods have been used to explore the molecular properties of the cyclo dipeptides in gaseous and aqueous phase environments. The self-consistent reaction field theory has been used to optimise the cyclopeptides in diethyl ether (? = 4.3) and water (? = 78.5), and the solvent effects have been analysed. A cluster of eight water molecules leads to the formation of first solvation shell of CGG and CGA and the strong H-bonding mainly contributes to the interaction energies. The H-bond interactions have been analysed by the calculation of electron density ρ(r) and its Laplacian ▽²ρ(r) at bond critical points using atoms in molecules theory. The natural bond orbital analysis was carried out to reveal the nature of H-bond interactions. In the solvated complexes, the keto carbons registered the maximum NMR chemical shifts.     

Dynamics of the Hck-SH3 domain: comparison of experiment with multiple molecular dynamics simulations

Molecular dynamics calculations provide a method by which the dynamic properties of molecules can be explored over timescales and at a level of detail that cannot be obtained experimentally from NMR or X-ray analyses. Recent work (Philippopoulos M, Mandel AM, Palmer AG III, Lim C, 1997, Proteins 28:481-493) has indicated that the accuracy of these simulations is high, as measured by the correspondence of parameters extracted from these calculations to those determined through experimental means. Here, we investigate the dynamic behavior of the Src homology 3 (SH3) domain of hematopoietic cell kinase (Hck) via 5N backbone relaxation NMR studies and a set of four independent 4 ns solvated molecular dynamics calculations. We also find that molecular dynamics simulations accurately reproduce fast motion dynamics as estimated from generalized order parameter (S2) analysis for regions of the protein that have experimentally well-defined coordinates (i.e., stable secondary structural elements). However, for regions where the coordinates are not well defined, as indicated by high local root-mean-square deviations among NMR-determined structural family members or high B-factors/low electron density in X-ray crystallography determined structures, the parameters calculated from a short to moderate length (less than 5-10 ns) molecular dynamics trajectory are dependent on the particular coordinates chosen as a starting point for the simulation.     

Specific stimulation of human apurinic/apyrimidinic endonuclease by heat shock protein 70

We previously demonstrated the stimulation of human apurinic/apyrimidinic endonuclease 1 (HAP1) by heat shock protein 70 (HSP70). In this work, we further defined the functional interaction between these proteins. Digestion of HSP70 by trypsin released 48 and 43 kDa amino terminal fragments that retained the ability to stimulate HAP1. In agreement with this result, an HSP70 N-terminal deletion mutant protein containing amino acids 1-385 was comparable to the full-length protein in its ability to enhance HAP1 activity. HSP70 mutants containing carboxy terminal amino acids 386-640 stimulated HAP1 only slightly, as did unrelated proteins. These results implicate the amino terminal portion of HSP70 in stimulating the activity of HAP1.     

Molecular dynamics simulation of proteins under high pressure: Structure,function and thermodynamics

BackgroundMolecular dynamics (MD) simulation is well-recognized as a powerful tool to investigate protein structure, function, and thermodynamics. MD simulation is also used to investigate high pressure effects on proteins. For conducting better MD simulation under high pressure, the main issues to be addressed are: (i) protein force fields and water models were originally developed to reproduce experimental properties obtained at ambient pressure; and (ii) the timescale to observe the pressure effect is often much longer than that of conventional MD simulations.Scope of reviewFirst, we describe recent developments in MD simulation methodologies for studying the high-pressure structure and dynamics of protein molecules. These developments include force fields for proteins and water molecules, and enhanced simulation techniques. Then, we summarize recent studies of MD simulations of proteins in water under high pressure.Major conclusionsRecent MD simulations of proteins in solution under pressure have reproduced various phenomena identified by experiments using high pressure, such as hydration, water penetration, conformational change, helix stabilization, and molecular stiffening.General significanceMD simulations demonstrate differences in the properties of proteins and water molecules between ambient and high-pressure conditions. Comparing the results obtained by MD calculations with those obtained experimentally could reveal the mechanism by which biological molecular machines work well in collaboration with water molecules.     

Enhancing Protein Adsorption Simulations by Using Accelerated Molecular Dynamics

The atomistic modeling of protein adsorption on surfaces is hampered by the different time scales of the simulation ( s) and experiment (up to hours), and the accordingly different ‘final’ adsorption conformations. We provide evidence that the method of accelerated molecular dynamics is an efficient tool to obtain equilibrated adsorption states. As a model system we study the adsorption of the protein BMP-2 on graphite in an explicit salt water environment. We demonstrate that due to the considerably improved sampling of conformational space, accelerated molecular dynamics allows to observe the complete unfolding and spreading of the protein on the hydrophobic graphite surface. This result is in agreement with the general finding of protein denaturation upon contact with hydrophobic surfaces.     

Tumour necrosis factor alpha-stimulated gene-6 inhibits osteoblastic differentiation of human mesenchymal stem cells induced by osteogenic differentiation medium and BMP-2

To better understand the molecular pathogenesis of OPLL (ossification of the posterior longitudinal ligament) of the spine, an ectopic bone formation disease, we performed cDNA microarray analysis on cultured ligament cells from OPLL patients. We found that TSG-6 (tumour necrosis factor alpha-stimulated gene-6) is down-regulated during osteoblastic differentiation. Adenovirus vector-mediated overexpression of TSG-6 inhibited osteoblastic differentiation of human mesenchymal stem cells induced by BMP (bone morphogenetic protein)-2 or OS (osteogenic differentiation medium). TSG-6 suppressed phosphorylation and nuclear accumulation of Smad 1/5 induced by BMP-2, probably by inhibiting binding of the ligand to the receptor, since interaction between TSG-6 and BMP-2 was observed in vitro. TSG-6 has two functional domains, a Link domain (a hyaluronan binding domain) and a CUB domain implicated in protein interaction. The inhibitory effect on osteoblastic differentiation was completely lost with exogenously added Link domain-truncated TSG-6, while partial inhibition was retained by the CUB domain-truncated protein. In addition, the inhibitory action of TSG-6 and the in vitro interaction of TSG-6 with BMP-2 were abolished by the addition of hyaluronan. Thus, TSG-6, identified as a down-regulated gene during osteoblastic differentiation, suppresses osteoblastic differentiation induced by both BMP-2 and OS and is a plausible target for therapeutic intervention in OPLL.     

Improved Ligand-Protein Binding Affinity Predictions Using Multiple Binding Modes

Accurate ligand-protein binding affinity prediction, for a set of similar binders, is a major challenge in the lead optimization stage in drug development. In general, docking and scoring functions perform unsatisfactorily in this application. Docking calculations, followed by molecular dynamics simulations and free energy calculations can be applied to improve the predictions. However, for targets with large, flexible binding sites, with no experimentally determined binding modes for a set of ligands, insufficient sampling can decrease the accuracy of the free energy calculations. Cytochrome P450s, a protein family of major importance for drug metabolism, is an example of a challenging target for binding affinity predictions. As a result, the choice of starting structure from the docking solutions becomes crucial. In this study, an iterative scheme is introduced that includes multiple independent molecular dynamics simulations to obtain weighted ensemble averages to be used in the linear interaction energy method. The proposed scheme makes the initial pose selection less crucial for further simulation, as it automatically calculates the relative weights of the various poses. It also properly takes into account the possibility that multiple binding modes contribute similarly to the overall affinity, or of similar compounds occupying very different poses. The method was applied to a set of 12 compounds binding to cytochrome P450 2C9 and it displayed a root mean-square error of 2.9 kJ/mol.