CO migration pathways in cytochrome P450cam studied by molecular dynamics simulations

Previous laser flash photolysis investigations between 100 and 300 K have shown that the kinetics of CO rebinding with cytochrome P450(cam)(camphor) consist of up to four different processes revealing a complex internal dynamics after ligand dissociation. In the present work, molecular dynamics simulations were undertaken on the ternary complex P450(cam)(cam)(CO) to explore the CO migration pathways, monitor the internal cavities of the protein, and localize the CO docking sites. One trajectory of 1 nsec with the protein in a water box and 36 trajectories of 1 nsec in the vacuum were calculated. In each trajectory, the protein contained only one CO ligand on which no constraints were applied. The simulations were performed at 200, 300, and 320 K. The results indicate the presence of seven CO docking sites, mainly hydrophobic, located in the same moiety of the protein. Two of them coincide with xenon binding sites identified by crystallography. The protein matrix exhibits eight persistent internal cavities, four of which corresponding to the ligand docking sites. In addition, it was observed that water molecules entering the protein were mainly attracted into the polar pockets, far away from the CO docking sites. Finally, the identified CO migration pathways provide a consistent interpretation of the experimental rebinding kinetics.     

Dissociation mechanism of GDP from Cdc42 via DOCK9 revealed by molecular dynamics simulations

Cell division control protein 42 homolog (Cdc42) influences a variety of cellular responses such as cell migration and polarity. Deregulation of Cdc42 has been associated with several human diseases and developmental disorders. Over-activation of Cdc42 through guanine nucleotide exchange factor (GEF) is a critical event for Cdc42 involved cancer metastasis. Members of DOCK family of GEF are important activators of Cdc42. However, this activation mechanism is still unknown. Molecular dynamics (MD) simulations and molecular mechanics-Poisson-Boltzmann surface area (MM-PBSA) calculations were employed to investigate the central step of the activation of Cdc42: the dissociation mechanism of GDP from Cdc42 via DOCK9. Simulation results show that Mg²⁺ ion has a remarkable influence on the conformational change of switch I of Cdc42 through residue Pro34 which functions as a “clasp” to control the flexibility of switch I. In the GDP dissociation process, the Mg²⁺ ion leave first to result in a suitable conformation of Cdc42 for following DOCK9 binding to. When DOCK9 binds to Cdc42, it changes the orientations of residues Lys16, Thr17, Cys18 and Phe28 of Cdc42 to weaken the interactions between Cdc42 and GDP to release GDP. This study first elucidates the dissociation mechanism of GDP from Cdc42 via DOCK9 and identifies the essential residues of Cdc42 in this process. These simulation results are consistent with the recent findings of biochemical and amino acid mutational studies, and the observations are beneficial to understand the activation mechanism of Cdc42 and to provide insights for designing compounds targeting on Cdc42 related cancer metastasis.     

Mechanisms for allosteric activation of protease DegS by ligand binding and oligomerization as revealed from molecular dynamics simulations

A local perturbation of a protein may lead to functional changes at some distal site, a phenomenon denoted as allostery. Here, we study the allosteric control of a protease using molecular dynamics simulations. The system considered is the bacterial protein DegS which includes a protease domain activated on ligand binding to an adjacent PDZ domain. Starting from crystallographic structures of DegS homo‐trimers, we perform simulations of the ligand‐free and ‐bound state of DegS at equilibrium. Considering a single protomer only, the trimeric state was mimicked by applying restraints on the residues in contact with other protomers in the DegS trimer. In addition, the bound state was also simulated without any restraints to mimic the monomer. Our results suggest that not only ligand release but also disassembly of a DegS trimer inhibits proteolytic activity. Considering various observables for structural changes, we infer allosteric pathways from the interface with other protomers to the active site. Moreover, we study how ligand release leads to (i) catalytically relevant changes involving residues 199–201 and (ii) a transition from a stretched to a bent conformation for residues 217–219 (which prohibits proper substrate binding). Finally, based on ligand‐induced C_α shifts we identify residues in contact with other protomers in the DegS trimer that likely transduce the perturbation from ligand release from a given protomer to adjacent protomers. These residues likely play a key role in the experimentally known effect of ligand release from a protomer on the proteolytic activity of the other protomers. Proteins 2016; 84:1690–1705. © 2016 Wiley Periodicals, Inc.     

Probing the "annealing" mechanism of GroEL minichaperone using molecular dynamics simulations

Although the intact chaperonin machinery is needed to rescue natural substrate proteins (SPs) under non-permissive conditions the "minichaperone" alone, containing only the isolated apical domain of GroEL, can assist folding of a certain class of proteins. To understand the annealing function of the minichaperone, we have carried out molecular dynamics simulations in the NPT ensemble totaling 300ns for four systems; namely, the isolated strongly binding peptide (SBP), the minichaperone, and the SBP and a weakly binding peptide (WBP) in complex with the minichaperone. The SBP, which is structureless in isolation, adopts a beta-hairpin conformation in complex with the minichaperone suggesting that favorable non-specific interactions of the SPs confined to helices H and I of the apical domains can induce local secondary structures. Comparison of the dynamical fluctuations of the apo and the liganded forms of the minichaperone shows that the stability (needed for SP capture) involves favorable hydrophobic interactions and hydrogen bond network formation between the SBP and WBP, and helices H and I. The release of the SP, which is required for the annealing action, involves water-mediated interactions of the charged residues at the ends of H and I helices. The simulation results are consistent with a transient binding release (TBR) model for the annealing action of the minichaperone. According to the TBR model, SP annealing occurs in two stages. In the first stage the SP is captured by the apical domain. This is followed by SP release (by thermal fluctuations) that places it in a different region of the energy landscape from which it can partition rapidly to the native state with probability Phi or be trapped in another misfolded state. The process of binding and release can result in enhancement of the native state yield. The TBR model suggests "that any cofactor that can repeatedly bind and release SPs can be effective in assisting protein folding." By comparing the structures of the non-chaperone alpha-casein (which has no sequence similarity with the apical domain) and the minichaperone and the hydrophobicity profiles we show that alpha-casein has a pair of helices that have similar sequence and structural profiles as H and I. Based on this comparison we identify residues that stabilize (destabilize) alpha-casein-protein complexes. This suggests that alpha-casein assists folding by the TBR mechanism.     

Homology modeling and molecular dynamics simulations of the glycine receptor ligand binding domain

We present a homology based model of the ligand binding domain (LBD) of the homopentameric alpha1 glycine receptor (GlyR). The model is based on multiple sequence alignment with other members of the nicotinicoid ligand gated ion channel superfamily and two homologous acetylcholine binding proteins (AChBP) from the freshwater (Lymnaea stagnalis) and saltwater (Aplysia californica) snails with known high resolution structure. Using two template proteins with known structure to model three dimensional structure of a target protein is especially advantageous for sequences with low homology as in the case presented in this paper. The final model was cross-validated by critical evaluation of experimental and published mutagenesis, functional and other biochemical studies. In addition, a complex structure with strychnine antagonist in the putative binding site is proposed based on docking simulation using Autodock program. Molecular dynamics (MD) simulations with simulated annealing protocol are reported on the proposed LBD of GlyR, which is stable in 5 ns simulation in water, as well as for a deformed LBD structure modeled on the corresponding domain determined in low-resolution cryomicroscopy structure of the alpha subunit of the full-length acetylcholine receptor (AChR). Our simulations demonstrate that the beta-sandwich central core of the protein monomer is fairly rigid in the simulations and resistant to deformations in water.     

Insights on pH-dependent conformational changes of mosquito odorant binding proteins by molecular dynamics simulations

Chemical recognition plays an important role for the survival and reproduction of many insect species. Odorant binding proteins (OBPs) are the primary components of the insect olfactory mechanism and have been documented to play an important role in the host-seeking mechanism of mosquitoes. They are “transport proteins” believed to transport odorant molecules from the external environment to their respective membrane targets, the olfactory receptors. The mechanism by which this transport occurs in mosquitoes remains a conundrum in this field. Nevertheless, OBPs have proved to be amenable to conformational changes mediated by a pH change in other insect species. In this paper, the effect of pH on the conformational flexibility of mosquito OBPs is assessed computationally using molecular dynamics simulations of a mosquito OBP “CquiOBP1” bound to its pheromone 3OG (PDB ID: 3OGN). Conformational twist of a loop, driven by a set of well-characterized changes in intramolecular interactions of the loop, is demonstrated. The concomitant (i) closure of what is believed to be the entrance of the binding pocket, (ii) expansion of what could be an exit site, and (iii) migration of the ligand towards this putative exit site provide preliminary insights into the mechanism of ligand binding and release of these proteins in mosquitoes. The correlation of our results with previous experimental observations based on NMR studies help us provide a cardinal illustration on one of the probable dynamics and mechanism by which certain mosquito OBPs could deliver their ligand to their membrane-bound receptors at specific pH conditions.     

Dynamics of Ca2+-saturated calmodulin D129N mutant studied by multiple molecular dynamics simulations

Fifteen independent 1-nsec MD simulations of fully solvated Ca(2+) saturated calmodulin (CaM) mutant D129N were performed from different initial conditions to provide a sufficient statistical basis to gauge the significance of observed dynamical properties. In all MD simulations the four Ca(2+) ions remained in their binding sites, and retained a single water ligand as observed in the crystal structure. The coordination of Ca(2+) ions in EF-hands I, II, and III was sevenfold. In EF-hand IV, which was perturbed by the mutation of a highly conserved Asp129, an anomalous eightfold Ca(2+) coordination was observed. The Ca(2+) binding loop in EF-hand II was observed to dynamically sample conformations related to the Ca(2+)-free form. Repeated MD simulations implicate two well-defined conformations of Ca(2+) binding loop II, whereas similar effect was not observed for loops I, III, and IV. In 8 out of 15 MD simulations Ca(2+) binding loop II adopted an alternative conformation in which the Thr62 >C=O group was displaced from the Ca(2+) coordination by a water molecule, resulting in the Ca(2+) ion ligated by two water molecules. The alternative conformation of the Ca(2+) binding loop II appears related to the "closed" state involved in conformational exchange previously detected by NMR in the N-terminal domain fragment of CaM and the C-terminal domain fragment of the mutant E140Q. MD simulations suggest that conformations involved in microsecond exchange exist partially preformed on the nanosecond time scale.     

Mutated form (G52E) of inactive diphtheria toxin CRM197: molecular simulations clearly display effect of the mutation to NAD binding

Mutated form (G52E) of diphtheria toxin (DT) CRM197 is an inactive and nontoxic enzyme. Here, we provided a molecular insight using comparative molecular dynamics (MD) simulations to clarify the influence of a single point mutation on overall protein and active-site loop. Post-processing MD analysis (i.e. stability, principal component analysis, hydrogen-bond occupancy, etc.) is carried out on both wild and mutated targets to investigate and to better understand the mechanistic differences of structural and dynamical properties on an atomic scale especially at nicotinamide adenine dinucleotide (NAD) binding site when a single mutation (G52E) happens at the DT. In addition, a docking simulation is performed for wild and mutated forms. The docking scoring analysis and docking poses results revealed that mutant form is not able to properly accommodate the NAD molecule.     

Effects of ligand binding on the dynamics of rice nonspecific lipid transfer protein 1: a model from molecular simulations

Plant nonspecific lipid transfer proteins (nsLTPs) are small, basic proteins constituted mainly of alpha-helices and stabilized by four conserved disulfide bridges. They are characterized by the presence of a tunnel-like hydrophobic cavity, capable of transferring various lipid molecules between lipid bilayers in vitro. In this study, molecular dynamics (MD) simulations were performed at room temperature to investigate the effects of lipid binding on the dynamic properties of rice nsLTP1. Rice nsLTP1, either in the free form or complexed with one or two lipids was subjected to MD simulations. The C-terminal loop was very flexible both before and after lipid binding, as revealed by calculating the root-mean-square fluctuation. After lipid binding, the flexibility of some residues that were not in direct contact with lipid molecules increased significantly, indicating an increase of entropy in the region distal from the binding site. Essential dynamics analysis revealed clear differences in motion between unliganded and liganded rice nsLTP1s. In the free form of rice nsLTP1, loop1 exhibited the largest directional motion. This specific essential motion mode diminished after binding one or two lipid molecules. To verify the origin of the essential motion observed in the free form of rice nsLTP1, we performed multiple sequence alignments to probe the intrinsic motion encoded in the primary sequence. We found that the amino acid sequence of loop1 is highly conserved among plant nsLTP1s, thus revealing its functional importance during evolution. Furthermore, the sequence of loop1 is composed mainly of amino acids with short side chains. In this study, we show that MD simulations, together with essential dynamics analysis, can be used to determine structural and dynamic differences of rice nsLTP1 upon lipid binding.     

Structural dynamics of precursor and product of the RNA enzyme from the hepatitis delta virus as revealed by molecular dynamics simulations

The hepatitis delta virus (HDV) ribozyme is a self-cleaving RNA enzyme involved in the replication of a human pathogen, the hepatitis delta virus. Recent crystal structures of the precursor and product of self-cleavage, together with detailed kinetic analyses, have led to hypotheses on the catalytic strategies employed by the HDV ribozyme. We report molecular dynamics (MD) simulations (approximately 120 ns total simulation time) to test the plausibility that specific conformational rearrangements are involved in catalysis. Site-specific self-cleavage requires cytidine in position 75 (C75). A precursor simulation with unprotonated C75 reveals a rather weak dynamic binding of C75 in the catalytic pocket with spontaneous, transient formation of a H-bond between U-1(O2') and C75(N3). This H-bond would be required for C75 to act as the general base. Upon protonation in the precursor, C75H+ has a tendency to move towards its product location and establish a firm H-bonding network within the catalytic pocket. However, a C75H+(N3)-G1(O5') H-bond, which would be expected if C75 acted as a general acid catalyst, is not observed on the present simulation timescale. The adjacent loop L3 is relatively dynamic and may serve as a flexible structural element, possibly gated by the closing U20.G25 base-pair, to facilitate a conformational switch induced by a protonated C75H+. L3 also controls the electrostatic environment of the catalytic core, which in turn may modulate C75 base strength and metal ion binding. We find that a distant RNA tertiary interaction involving a protonated cytidine (C41) becomes unstable when left unprotonated, leading to disruptive conformational rearrangements adjacent to the catalytic core. A Na ion temporarily compensates for the loss of the protonated hydrogen bond, which is strikingly consistent with the experimentally observed synergy between low pH and high Na+ concentrations in mediating residual self-cleavage of the HDV ribozyme in the absence of divalents.     

Retardation of the unfolding process by single N-glycosylation of ribonuclease A based on molecular dynamics simulations

The conformational characteristics of glycosylated- and unglycosylated bovine pancreatic ribonuclease A (RNaseA) were traced with unfolding molecular dynamics simulations using CHARMM program at 470 K. The glycosylated RNase (Glc_RNase) possesses nearly identical protein structure with RNaseA, differing only by presence of a single acetylglucosamine residue N-linked to Asn34 in the RNaseA. Attaching of monomeric N-acetylglucosamine residue to the Asn34 in RNaseA resulted in a change of denaturing process of Glc_RNase. Simulations showed that the unfolding of RNaseA involved significant weakening of nonlocal interactions whereas the glycosylation led Glc_RNase to preserve the nonlocal interactions even in its denatured form. Even in simulations over 8 ns at 470 K, Glc_RNase remained relatively stable as a less denatured conformation. However, conformation of RNaseA was changed to a fully unfolded state before 3 ns of the simulations at 470 K. This difference was due to fact that formation of hydrogen bond bridges and nonlocal contacts induced by the attached N-acetylglucosamine of Glc_RNase showing in the unfolding simulations. These high-temperature unfolding MD simulations provided a theoretical basis for the previous experimental work in which Glc_RNase showed slower unfolding kinetics compared with unglycosylated RNaseA, suggesting that single N-glycosylation induced retardation of unfolding process of the ribonuclease protein.     

Characterization of the denaturation of human alpha-lactalbumin in urea by molecular dynamics simulations

Molecular dynamics (MD) simulations were used to characterize the non-cooperative denaturation of the molten globule A-state of human alpha-lactalbumin by urea. A solvent of explicit urea and water molecules was used, corresponding to a urea concentration of approximately 6M. Three simulations were performed at temperatures of 293K, 360K and 400K, with lengths of 2 ns, 8 ns and 8 ns respectively. The results of the simulations were compared with experimental data from NMR studies of human alpha-lactalbumin and related peptides. During the simulations, hydrogen bonds were formed from the protein to both urea and water molecules as intra-protein hydrogen bonds were lost. Urea was shown to compete efficiently with water as both a hydrogen bond donor and acceptor. Radial distribution functions of water and urea around hydrophobic side chain atoms showed a significant increase in urea molecules in the solvation shell as the side chains became exposed during denaturation. A considerable portion of the native-like secondary structure persisted throughout the simulations. However, in the simulations at 360K and 400K, there were substantial changes in the packing of aromatic and other hydrophobic side chains in the protein, and many native contacts were lost. The results suggest that during the non-cooperative denaturation of the molten globule, secondary structure elements are stabilized by non-specific, non-native interactions.     

Mutant bovine odorant-binding protein: Temperature affects the protein stability and dynamics as revealed by infrared spectroscopy and molecular dynamics simulations

Bovine odorant-binding protein (bOBP), a member of the lipocalin family, presents the so-called 3D "domain-swapped" protein structure. In fact, in solution, it appears as a dimer in which each monomer is composed by the classical lipocalin fold, with a central beta-barrel followed by a stretch of residues and the alpha-helix domain protruding out of the barrel and crossing the dimer interface. Recently, a deswapped mutant form of bOBP was obtained, in which a Gly residue was inserted after position 121 and the two residues in position 64 and 156 were replaced by Cys residues for restoring the disulfide bridge common to the lipocalin family. In this work, we used Fourier transform infrared spectroscopy and molecular dynamics simulations to investigate the effect of temperature on the structural stability and conformational dynamics of the mutant bOBP. The spectroscopic and molecular simulation data pointed out that the hydrophobic regions of the protein matrix appear to be an important factor for the protein stability and integrity. In addition, it was also found that the mutant bOBP is significantly stabilized by the binding of the ligand, which may have an impact on the biological function of bOBP. The obtained results will allow for a better use of this protein as probe for the design of advanced protein-based biosensors for the detection of compounds used in the fabrication of explosive powders.     

A computational insight into binding modes of inhibitors XD29, XD35, and XD28 to bromodomain-containing protein 4 based on molecular dynamics simulations

In the current work, conformational changes of bromodomain-containing protein 4 (1) (BRD4-1) induced by bindings of inhibitors XD29 (57G), XD35 (57F), and XD28 (L28) were investigated using molecular dynamics (MD) simulations and principal component analysis. The results demonstrate that inhibitor bindings produce significant effect on the motion of ZA loop in BRD4-1. Moreover, to further study binding modes of three inhibitors to BRD4-1, binding free energies of inhibitors to BRD4-1 were also calculated using molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method. The results indicate that van der Waals interactions are main factors to modulate inhibitor bindings. Energy decomposition and hydrogen bond analysis demonstrate that residues Pro82, Leu92, Asn140, and Ile146 play important roles in binding processes of inhibitors to BRD4-1. This study is not only helpful for better understanding function and internal dynamics of BRD4-1, but also can provide a theoretical basis for rational designs of effective anticancer drugs targeting BRD4-1.     

Conformational changes of the p53-binding cleft of MDM2 revealed by molecular dynamics simulations

Two 35-ns molecular dynamics simulations of both ligated [mouse double minute protein 2 (MDM2(p53))] and unligated (MDM2(apo)) structures of human MDM2 bound to the N-terminal domain of the tumor suppressor p53 have been performed. Analysis of the dynamics revealed that the most flexible region of MDM2 was the p53-binding cleft. When MDM2 was bound to p53, a wider and more stable topology of the cleft was obtained, while unligated MDM2 showed a narrower and highly flexible cleft. It was also found that the dynamics involved in the opening/closing motions were due to the movement of different domains of the protein, which is in agreement with recent experimental data. Considering our results, a mechanism in which p53 might be recognized and attached to MDM2 is proposed, and some implications on future directions for in silico anticancer drug design efforts are discussed. In summary, the observations made here would be very useful not only for better understanding of the biological implications of the MDM2 dynamics, but also for future efforts in anticancer drug design and discovery.     

On the structure of molecularly imprinted polymers by modifying charge on functional groups through molecular dynamics simulations

Molecularly imprinted polymers (MIPs) have been widely applied in many fields owing to their advantages. The recognition mechanism between target molecules and MIPs and the influence of dominant factor on the recognition process are still poorly understood. In this paper, a cubic methacrylate-based MIP model was constructed, and the charge on carboxyl group atoms was changed artificially to investigate the recognition process. It is found that the diffusion coefficients of the target molecules (cholesterol) are not affected by polymer network structure. The recognition process is mainly determined by the mesh sizes of the polymer network. In addition, the structure of modified MIP systems was also discussed from the viewpoints of radial distribution function and hydrogen bonds system. These results suggest that the polymer matrix structure would be enhanced with an increase in charge. Thus, it influences the structure of the water molecules in the system a little.     

Different effects of lectins on the ligand binding of the NMDA receptors and sigma sites in rat brain hippocampus synaptic membranes

The effects of the lectins concanavalin A, WGA, ricin, abrin, and the mistletoe lectins from Viscum album MLI, MLII, and MLIII on the binding of ligands of the NMDA and sigma receptors in rat hippocampus synaptic plasma membranes were investigated. Binding of [³H]MK-801, [³H]glutamate, [³H]5,7-DCKA, and [³H]glycine to the membranes was decreased by 40-60% after addition of galactose-specific lectins (mistletoe lectins MLI, MLII, ricin, abrin) at concentrations of 0.01 mg/ml, but was not affected by the glucose- and mannose-specific lectin Con A, an acetylglucosamine-specific lectin WGA, or an acetylgalactosamine-specific lectin MLIII. The binding of [³H]SKF 10047 was decreased only in the presence of MLIII and did not change after addition of the other lectins. It is suggested that lectin-sensitive ligand binding sites of sigma- and NMDA receptors are located separately, and that the carbohydrate side chains of the sigma receptor do not participate in the modulation of the NMDA-receptor.     

Influence of dendrimer surface chemistry and pH on the binding and release pattern of chalcone studied by molecular dynamics simulations

Poly(amidoamine) (PAMAM) dendrimers are promising nanocarriers that can enhance the solubility of hydrophobic drugs. The surface chemistry of dendrimers is of great relevance as end groups of these nanocarriers can be easily modified to improve the bioavailability and sustained release of the cargo. Therefore, a molecular‐level understanding of the host‐guest interactions that can give both qualitative and quantitative information is particularly desirable. In this work, fully atomistic molecular dynamics simulations were used to study the association of a bioactive natural product, ie, chalcone, with amine‐, acetyl‐, and carboxyl‐terminated PAMAM dendrimers at physiological and acidic pH environments. Amine‐ and carboxyl‐terminated PAMAM dendrimers have an open microstructure at low pH that is not able to hold the ligand tightly, resulting in an unfavorable encapsulation of the chalcone molecule. In the case of acetyl‐terminated dendrimer, chalcone molecule diffuses out of the dendritic cavities a few times during the simulation time and prefers to locate close to the surface of dendrimer. Average center of mass distance values at neutral pH showed that the chalcone molecule bounds firmly in the internal pockets of amine‐, acetyl‐, and carboxyl‐terminated dendrimers and forms stable complexes with these nanovectors. The potential of mean force calculations showed that the release of the ligand from the dendrimers occurs at a controlled rate in the body.