Nucleotide-mediated conformational changes of monomeric actin and Arp3 studied by molecular dynamics simulations

Members of the actin family of proteins exhibit different biochemical properties when ATP, ADP-P_i, ADP, or no nucleotide is bound. We used molecular dynamics simulations to study the effect of nucleotides on the behavior of actin and actin-related protein 3 (Arp3). In all of the actin simulations, the nucleotide cleft stayed closed, as in most crystal structures. ADP was much more mobile within the cleft than ATP, despite the fact that both nucleotides adopt identical conformations in actin crystal structures. The nucleotide cleft of Arp3 opened in most simulations with ATP, ADP, and no bound nucleotide. Deletion of a C-terminal region of Arp3 that extends beyond the conserved actin sequence reduced the tendency of the Arp3 cleft to open. When the Arp3 cleft opened, we observed multiple instances of partial release of the nucleotide. Cleft opening in Arp3 also allowed us to observe correlated movements of the phosphate clamp, cleft mouth, and barbed-end groove, providing a way for changes in the nucleotide state to be relayed to other parts of Arp3. The DNase binding loop of actin was highly flexible regardless of the nucleotide state. The conformation of Ser14/Thr14 in the P1 loop was sensitive to the presence of the γ-phosphate, but other changes observed in crystal structures were not correlated with the nucleotide state on nanosecond timescales. The divalent cation occupied three positions in the nucleotide cleft, one of which was not previously observed in actin or Arp2/3 complex structures. In sum, these simulations show that subtle differences in structures of actin family proteins have profound effects on their nucleotide-driven behavior.     

Effect O6‐guanine alkylation on DNA flexibility studied by comparative molecular dynamics simulations

Alkylation of guanine at the O6 atom is a highly mutagenic DNA lesion because it alters the coding specificity of the base causing G:C to A:T transversion mutations. Specific DNA repair enzymes, e.g. O⁶‐alkylguanin‐DNA‐Transferases (AGT), recognize and repair such damage after looping out the damaged base to transfer it into the enzyme active site. The exact mechanism how the repair enzyme identifies a damaged site within a large surplus of undamaged DNA is not fully understood. The O⁶‐alkylation of guanine may change the deformability of DNA which may facilitate the initial binding of a repair enzyme at the damaged site. In order to characterize the effect of O⁶‐methyl‐guanine (O⁶‐MeG) containing base pairs on the DNA deformability extensive comparative molecular dynamics (MD) simulations on duplex DNA with central G:C, O⁶‐MeG:C or O⁶‐MeG:T base pairs were performed. The simulations indicate significant differences in the helical deformability due to the presence of O⁶‐MeG compared to regular undamaged DNA. This includes enhanced base pair opening, shear and stagger motions and alterations in the backbone fine structure caused in part by transient rupture of the base pairing at the damaged site and transient insertion of water molecules. It is likely that the increased opening motions of O⁶‐MeG:C or O⁶‐MeG:T base pairs play a decisive role for the induced fit recognition or for the looping out of the damaged base by repair enzymes. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 23–32, 2015.     

Aggregation of small peptides studied by molecular dynamics simulations

Peptides and proteins tend to aggregate under appropriate conditions. The amyloid fibrils that are ubiquitously found among these structures are associated with major human diseases like Alzheimer's disease, type II diabetes, and various prion diseases. Lately, it has been observed that even very short peptides like tetra and pentapeptides can form ordered amyloid structures. Here, we present aggregation studies of three such small polypeptide systems, namely, the two amyloidogenic peptides DFNKF and FF, and a control (nonamyloidogenic) one, the AGAIL. The respective aggregation process is studied by all-atom Molecular Dynamics simulations, which allow to shed light on the fine details of the association and aggregation process. Our analysis suggests that naturally aggregating systems exhibit significantly diverse overall cluster shape properties and specific intermolecular interactions. Additional analysis was also performed on the previously studied NFGAIL system.     

Virus capsid dissolution studied by microsecond molecular dynamics simulations

Dissolution of many plant viruses is thought to start with swelling of the capsid caused by calcium removal following infection, but no high-resolution structures of swollen capsids exist. Here we have used microsecond all-atom molecular simulations to describe the dynamics of the capsid of satellite tobacco necrosis virus with and without the 92 structural calcium ions. The capsid expanded 2.5% upon removal of the calcium, in good agreement with experimental estimates. The water permeability of the native capsid was similar to that of a phospholipid membrane, but the permeability increased 10-fold after removing the calcium, predominantly between the 2-fold and 3-fold related subunits. The two calcium binding sites close to the icosahedral 3-fold symmetry axis were pivotal in the expansion and capsid-opening process, while the binding site on the 5-fold axis changed little structurally. These findings suggest that the dissociation of the capsid is initiated at the 3-fold axis.     

Protein flexibility using constraints from molecular dynamics simulations

Proteins are held together in the native state by hydrophobic interactions, hydrogen bonds and interactions with the surrounding water, whose strength as well as spatial and temporal distribution affects protein flexibility and hence function. We study these effects using 10 ns molecular dynamics simulations of pure water and of two proteins, the glutamate receptor ligand binding domain and barnase. We find that most of the noncovalent interactions flicker on and off over typically nanoseconds, and so we can obtain good statistics from the molecular dynamics simulations. Based on this information, a topological network of rigid bonds corresponding to a protein structure with covalent and noncovalent bonds is constructed, with account being taken of the influence of the flickering hydrogen bonds. We define the duty cycle for the noncovalent interactions as the percentage of time a given interaction is present, which we use as an input to investigate flexibility/rigidity patterns, in the algorithm FIRST which constructs and analyses topological networks.     

Structural determinants of MscL gating studied by molecular dynamics simulations

The mechanosensitive channel of large conductance (MscL) in prokaryotes plays a crucial role in exocytosis as well as in the response to osmotic downshock. The channel can be gated by tension in the membrane bilayer. The determination of functionally important residues in MscL, patch-clamp studies of pressure-conductance relationships, and the recently elucidated crystal structure of MscL from Mycobacterium tuberculosis have guided the search for the mechanism of MscL gating. Here, we present a molecular dynamics study of the MscL protein embedded in a fully hydrated POPC bilayer. Simulations totaling 3 ns in length were carried out under conditions of constant temperature and pressure using periodic boundary conditions and full electrostatics. The protein remained in the closed state corresponding to the crystal structure, as evidenced by its impermeability to water. Analysis of equilibrium fluctuations showed that the protein was least mobile in the narrowest part of the channel. The gating process was investigated through simulations of the bare protein under conditions of constant surface tension. Under a range of conditions, the transmembrane helices flattened as the pore widened. Implications for the gating mechanism in light of these and experimental results are discussed.     

Disaccharide conformational flexibility. II. Molecular dynamics simulations of sucrose

Molecular dynamics simulations have been used to study the motions in vacuum of the disaccharide sucrose. Ensembles of trajectories were calculated for each of the five local minimum energy conformations identified in the adiabatic conformational energy mapping of this molecule. The model sucrose molecules were found to exhibit a variety of motions, although the global minimum energy conformation was found to be dynamically stable, and no transitions away from this structure were observed to occur spontaneously. In all but one of these vacuum trajectories, the intramolecular hydrogen bond between residues was maintained, in accord with recent nmr studies of this molecule in aqueous solution. Considerable flexibility of the furanoid ring was found in the trajectories. No "flips" to the opposite puckering for this ring were found in the simulations starting from the global minimum, although such a transition was observed for a trajectory initiated with one of the higher local minimum energy conformations. Overall, the observed structural fluctuations were consistent with the experimental picture of sucrose as a relatively rigid molecule.     

Ribonucleotide reductase inhibition by p-alkoxyphenols studied by molecular docking and molecular dynamics simulations

Ribonucleotide reductase (RNR) is necessary for production of the precursor deoxyribonucleotides for DNA synthesis. Class Ia RNR functions via a stable free radical in one of the two components protein R2. The enzyme mechanism involves long range (proton coupled) electron transfer between protein R1 and the tyrosyl radical in protein R2. Earlier experimental studies showed that p-alkoxyphenols inhibit RNR. Here, molecular docking and molecular dynamics simulations involving protein R2 suggest an inhibition mechanism for p-alkoxyphenols . A low energy binding pocket is identified in protein R2. The preferred configuration provides a structural basis explaining their specific binding to the Escherichia coli and mouse R2 proteins. Trp48 (E. coli numbering), on the electron transfer pathway, is involved in the interactions with the inhibitors. The relative order of the binding energies calculated for the phenol derivatives to protein R2 is correlated with earlier experimental data on inhibition efficiency, in turn related to increasing size of the hydrophobic alkyl substituents. Using the configuration identified by molecular docking as a starting point for molecular dynamics simulations, we find that the p-allyloxyphenol interrupts the catalytic electron transfer pathway of the R2 protein by forming hydrogen bonds with Trp48 and Asp237, thus explaining the inhibitory activity of p-alkoxyphenols.     

Influence of point mutations on the flexibility of cytochrome b5: molecular dynamics simulations of holoproteins

Two membrane-bound isoforms of cytochrome b5 have been identified in mammals, one associated with the outer mitochondrial membrane (OM b5) and the other with the endoplasmic reticulum (microsomal, or Mc b5). The soluble heme binding domains of OM and Mc b5 have highly similar three-dimensional structures but differ significantly in physical properties, with OM b5 exhibiting higher stability due to stronger heme association. In this study, we present results of 8.5-ns length molecular dynamics simulations for rat Mc b5, bovine Mc b5, and rat OM b5, as well as for two rat OM b5 mutants that were anticipated to exhibit properties intermediate between those of rat OM b5 and the two Mc proteins: the A18S/I32L/L47R triple mutant (OM3M) and the A18S/I25L/I32L/L47R/L71S quintuple mutant (OM5M). Analysis of the structure, fluctuations, and interactions showed that the five b5 variants used in this study differed in organization of their molecular surfaces and heme binding cores in a way that could be used to explain certain experimentally observed physical differences. Overall, our simulations provided qualitative microscopic explanations of many of the differences in physical properties between OM and Mc b5 and two mutants in terms of localized changes in structure and flexibility. They also reveal that opening of a surface cleft between hydrophobic cores 1 and 2 in bovine Mc b5, observed in two previously reported simulations (E. M. Storch and V. Daggett, Biochemistry, 1995, Vol. 34, pp. 9682-9693; A. Altuve, Biochemistry, 2001, Vol. 40, pp. 9469-9483), probably resulted from removal of crystal contacts and likely does not occur on the nanosecond time scale. Finally, the MD simulations of OM5M b5 verify that stability and dynamic properties of cytochrome b5 are remarkably resistant to mutations that dramatically alter the stability and structure of the apoprotein.     

Sequence-dependent DNA deformability studied using molecular dynamics simulations

Proteins recognize specific DNA sequences not only through direct contact between amino acids and bases, but also indirectly based on the sequence-dependent conformation and deformability of the DNA (indirect readout). We used molecular dynamics simulations to analyze the sequence-dependent DNA conformations of all 136 possible tetrameric sequences sandwiched between CGCG sequences. The deformability of dimeric steps obtained by the simulations is consistent with that by the crystal structures. The simulation results further showed that the conformation and deformability of the tetramers can highly depend on the flanking base pairs. The conformations of xATx tetramers show the most rigidity and are not affected by the flanking base pairs and the xYRx show by contrast the greatest flexibility and change their conformations depending on the base pairs at both ends, suggesting tetramers with the same central dimer can show different deformabilities. These results suggest that analysis of dimeric steps alone may overlook some conformational features of DNA and provide insight into the mechanism of indirect readout during protein–DNA recognition. Moreover, the sequence dependence of DNA conformation and deformability may be used to estimate the contribution of indirect readout to the specificity of protein–DNA recognition as well as nucleosome positioning and large-scale behavior of nucleic acids.     

Molecular dynamics simulations of Trp side-chain conformational flexibility in the gramicidin A channel

Gramicidin A (gA) is prototypical peptide antibiotic and a model ion channel former. Configured in the solid-state NMR beta(6.5)-helix channel conformation, gA was subjected to 1-ns molecular dynamics (MD) gas phase simulations using the all-atom charmm22 force field to ascertain the conformational stability of the Trp side chains as governed by backbone and neighboring side-chain contacts. Three microcanonical trajectories were computed using different initial atomic velocities for each of twenty different initial structures. For each set, one of the four Trp side chains in each monomer was initially positioned in one of the five non-native conformations (A. E. Dorigo et al., Biophysical Journal, 1999, Vol. 76, 1897-1908), the other Trps being positioned in the native state, o1. In three additional control simulations, all Trps were initiated in the native conformation. After equilibration, constraints were removed and subsequent conformational changes of the initially constrained Trp were measured. The chi(1) was more flexible than chi(2.1). The energetically optimal orientation, o1 (Dorigo et al., 1999), was the most stable in all four Trp positions (9, 11, 13, 15) and remained unchanged for the entire 1 ns simulation in 19 of 24 trials. Changes in chi(1) from each of the 5 suboptimal states occur readily. Two of the non-native conformations reverted readily to o1, whereas the other three converted to an intermediate state, i2. There were frequent interconversions between i2 and o1. We speculate that experimentally observed Trp stability is caused by interactions with the lipid-water interface, and that stabilization of one of the suboptimal conformations in gA, such as i2, by lipid headgroups could produce a secondary, metastable conformational state. This could explain recent experimental studies of differences in the channel conductance dispersity between gA and a Trp-to-Phe gA analog, gramicidin M (gM, J. C. Markham et al., Biochimica et Biophysica Acta, 2001, Vol. 1513, 185-192).     

Potassium and sodium ions in a potassium channel studied by molecular dynamics simulations

We have performed simulations of both a single potassium ion and a single sodium ion within the pore of the bacterial potassium channel KcsA. For both ions there is a dehydration energy barrier at the cytoplasmic mouth suggesting that the crystal structure is a closed conformation of the channel. There is a potential energy barrier for a sodium ion in the selectivity filter that is not seen for potassium. Radial distribution functions for both ions with the carbonyl oxygens of the selectivity filter indicate that sodium may interact more tightly with the filter than does potassium. This suggests that the key to the ion selectivity of KcsA is the greater dehydration energy of Na⁺ ions, and helps to explain the block of KcsA by internal Na⁺ ions.     

Early events in protein aggregation: molecular flexibility and hydrophobicity/charge interaction in amyloid peptides as studied by molecular dynamics simulations

In a previous article (Zbilut et al., Biophys J 2003;85:3544-3557), we demonstrated how an aggregation versus folding choice could be approached considering hydrophobicity distribution and charge. In this work, our aim is highlighting the mutual interaction of charge and hydrophobicity distribution in the aggregation process. Use was made of two different peptides, both derived from a transmembrane protein (amyloid precursor protein; APP), namely, Abeta(1-28) and Abeta(1-40). Abeta(1-28) has a much lower aggregation propensity than Abeta(1-40). The results obtained by means of molecular dynamics simulations show that, when submitted to the most "aggregation-prone" environment, corresponding to the isoelectric point and consequently to zero net charge, both peptides acquire their maximum flexibility, but Abeta(1-40) has a definitely higher conformational mobility than Abeta(1-28). The absence of a hydrophobic "tail," which is the most mobile part of the molecule in Abeta(1-40), is the element lacking in Abeta(1-28) for obtaining a "fully aggregating" phenotype. Our results suggest that conformational flexibility, determined by both hydrophobicity and charge effect, is the main mechanistic determinant of aggregation propensity.     

Influence of the environment in the conformation of alpha-helices studied by protein database search and molecular dynamics simulations

The influence of the solvent on the main-chain conformation (phi and Psi dihedral angles) of alpha-helices has been studied by complementary approaches. A first approach consisted in surveying crystal structures of both soluble and membrane proteins. The residues of analysis were further classified as exposed to either the water (polar solvent) or the lipid (apolar solvent) environment or buried to the core of the protein (intermediate polarity). The statistical results show that the more polar the environment, the lower the value of phi(i) and the higher the value of Psi(i) are. The intrahelical hydrogen bond distance increases in water-exposed residues due to the additional hydrogen bond between the peptide carbonyl oxygen and the aqueous environment. A second approach involved nanosecond molecular dynamics simulations of poly-Ala alpha-helices in environments of different polarity: water to mimic hydrophilic environments that can form hydrogen bonds with the peptide carbonyl oxygen and methane to mimic hydrophobic environments without this hydrogen bond capabilities. These simulations reproduce similar effects in phi and Psi angles and intrahelical hydrogen bond distance and angle as observed in the protein survey analysis. The magnitude of the intrahelical hydrogen bond in the methane environment is stronger than in the water environment, suggesting that alpha-helices in membrane-embedded proteins are less flexible than in soluble proteins. There is a remarkable coincidence between the phi and Psi angles obtained in the analysis of residues exposed to the lipid in membrane proteins and the results from computer simulations in methane, which suggests that this simulation protocol properly mimic the lipidic cell membrane and reproduce several structural characteristics of membrane-embedded proteins. Finally, we have compared the phi and Psi torsional angles of Pro kinks in membrane protein crystal structures and in computer simulations.     

Mechanically induced titin kinase activation studied by force-probe molecular dynamics simulations

The conversion of mechanical stress into a biochemical signal in a muscle cell requires a force sensor. Titin kinase, the catalytic domain of the elastic muscle protein titin, has been suggested as a candidate. Its activation requires major conformational changes resulting in the exposure of its active site. Here, force-probe molecular dynamics simulations were used to obtain insight into the tension-induced activation mechanism. We find evidence for a sequential mechanically induced opening of the catalytic site without complete domain unfolding. Our results suggest the rupture of two terminal beta-sheets as the primary unfolding steps. The low force resistance of the C-terminal relative to the N-terminal beta-sheet is attributed to their different geometry. A subsequent rearrangement of the autoinhibitory tail is seen to lead to the exposure of the active site, as is required for titin kinase activity. These results support the hypothesis of titin kinase as a force sensor.     

Influence of a cis,syn‐cyclobutane pyrimidine dimer damage on DNA conformation studied by molecular dynamics simulations

The photo‐induced formation of cis‐syn‐cyclobutane pyrimidine dimers (CPD) is a highly mutagenic and cancerogenic DNA lesion. In bacteria photolyases can efficiently reverse the dimer formation employing a light‐driven reaction after looping out the CPD damaged bases into the enzyme active site. The exact mechanism how the repair enzyme identifies a damaged site within a large surplus of undamaged DNA is not fully understood. The CPD damage may alter the DNA structure and dynamics already in the absence of the repair enzyme which can facilitate the initial binding of a photolyase repair enzyme. To characterize the effect of a CPD damage, extensive comparative molecular dynamics (MD) simulations on duplex DNA with central regular or CPD damaged nucleotides were performed supplemented with simulations of the DNA‐photolyase complex. Although no spontaneous flipping out transitions of the damaged bases were observed, the simulations showed significant differences in the conformational states of regular and CPD damage DNA. The isolated damaged DNA adopted transient conformations which resembled the global shape of the repair enzyme bound conformation more closely compared to regular B‐DNA. In particular, these conformational changes were observed in most of helical and structural parameters where the protein bound DNA differs drastically from regular B‐DNA. It is likely that the transient overlap of isolated DNA with the enzyme bound DNA conformation plays a decisive role for the specific and rapid initial recognition by the repair enzyme prior to the looping out process of the damaged DNA. © 2014 Wiley Periodicals, Inc. Biopolymers 103: 215–222, 2015.     

Orientation and conformation of a lipase at an interface studied by molecular dynamics simulations

Electron density profiles calculated from molecular dynamics trajectories are used to deduce the orientation and conformation of Thermomyces lanuginosa lipase and a mutant adsorbed at an air-water interface. It is demonstrated that the profiles display distinct fine structures, which uniquely characterize enzyme orientation and conformation. The density profiles are, on the nanosecond timescale, determined by the average enzyme conformation. We outline a computational scheme that from a single molecular dynamics trajectory allows for extraction of electron density profiles referring to different orientations of the lipase relative to an implicit interface. Profiles calculated for the inactive and active conformations of the lipase are compared with experimental electron density profiles measured by x-ray reflectivity for the lipase adsorbed at an air-water interface. The experimental profiles contain less fine structural information than the calculated profiles because the resolution of the experiment is limited by the intrinsic surface roughness of water. Least squares fits of the calculated profiles to the experimental profiles provide areas per adsorbed enzyme and suggest that Thermomyces lanuginosa lipase adsorbs to the air-water interface in a semiopen conformation with the lid oriented away from the interface.     

Relative flexibility of DNA and RNA: a molecular dynamics study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.     

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.     

Characterization of the conformational and orientational dynamics of ganglioside GM1 in a dipalmitoylphosphatidylcholine bilayer by molecular dynamics simulations

The structure and dynamics of a single GM1 (Gal5-β1,3-GalNAc4-β1,4-(NeuAc3-α2,3)-Gal2-β1,4-Glc1-β1,1-Cer) embedded in a DPPC bilayer have been studied by MD simulations. Eleven simulations, each of 10 ns productive run, were performed with different initial conformations of GM1. Simulations of GM1-Os in water and of a DPPC bilayer were also performed to delineate the effects of the bilayer and GM1 on the conformational and orientational dynamics of each other. The conformation of the GM1 headgroup observed in the simulations is in agreement with those reported in literature; but the headgroup is restricted when embedded in the bilayer. NeuAc3 is the outermost saccharide towards the water phase. Glc1 and Gal2 prefer a parallel, and NeuAc3, GalNac4 and Gal5 prefer a perpendicular, orientation with respect to the bilayer normal. The overall characteristics of the bilayer are not affected by the presence of GM1; however, GM1 does influence the DPPC molecules in its immediate vicinity. The implications of these observations on the specific recognition and binding of GM1 embedded in a lipid bilayer by exogenous proteins as well as proteins embedded in lipids have been discussed.