Molecular dynamics simulation and essential dynamics study of mutated plastocyanin: structural, dynamical and functional effects of a disulfide bridge insertion at the protein surface

A molecular dynamics simulation (1.1 ns) at 300 K, of fully hydrated Ile21Cys, Glu25Cys plastocyanin mutant has been performed to investigate the structural, dynamical and functional effects of a disulfide bridge insertion at the surface of the protein. A detailed analysis of the root mean square fluctuations, H-bonding pattern and dynamical cross-correlation map has been performed. An essential dynamics method has also been applied as complementary analysis to identify concerted motions (essential modes), that could be relevant to the electron transfer function. The results have been compared with those previously obtained for wild-type plastocyanin and have revealed that the mutant shows a different pattern of H-bonds, with several interactions lost and a higher flexibility, especially around the electron transfer copper site. The analysis of dynamical cross-correlation map and of essential modes, has shown that the mutant performs different functional concerted motions, which might be related to the binding recognition with its electron transfer partners in comparison with the wild-type protein.     

Concerted motions in copper plastocyanin and azurin: an essential dynamics study

Essential dynamics analysis of molecular dynamics simulation trajectories (1.1 ns) of two copper containing electron transfer proteins, plastocyanin and azurin, has been performed. The protein essential modes have been analysed in order to identify large concerted motions which could be relevant for the electron transfer function exerted by these proteins. The analysis, conducted for temporal windows of different lengths along the protein trajectories, shows a rapid convergence and indicates that for both the proteins the predominant internal motions occur in a subspace of only a few degrees of freedom. Moreover, it is found that for both the proteins the likely binding sites (i.e. the hydrophobic and negative patches) with the reaction partners move in a concerted fashion with a few structural regions far from the active site. Such results are discussed in connection with the possible involvement of large concerted motions in the recognition and binding interaction with physiological electron transfer partners.     

Evidence of reduced flexibility in disulfide bridge-depleted azurin: a molecular dynamics simulation study.

Two molecular dynamics simulations have been performed for 2 ns, at room temperature, on fully hydrated wild type and Cys3Ala/Cys26Ala double-mutant azurin, to investigate the role of the unique disulfide bridge on the structure and dynamics of the protein. The results show that the removal of the [bond]SS[bond] bond does not affect the structural features of the protein, whereas alterations of the dynamical properties are observed. The root mean square fluctuations of the atomic positions are, on average, considerably reduced in the azurin mutant with respect to the wild type form. The number of intramolecular hydrogen bonds between protein backbone atoms that are lost during the simulation, with respect to the starting configuration, are reduced in the absence of the disulfide bond. The analysis of the dynamical cross-correlation map, characterising the protein co-ordinated internal motions, demonstrates in the mutated azurin a significant decrease in anti-correlated displacements between protein residues, with the only exception occurring in the region of the mutation sites. The overall findings show a relevant reduction in flexibility as a consequence of the disulfide bridge depletion in azurin, suggesting that the [bond]SS[bond] bond is a structural element which significantly contributes to the dynamic properties of the native protein.     

Structural, dynamical and functional aspects of the inner motions in the blue copper protein azurin

Molecular dynamics was applied to dissect out the internal motions of azurin, a copper protein performing electron transfer. Simulations of 16.5 ns were analyzed in search of coordinated displacements of amino acid residues that are important for the protein function. A region with high conformational instability was found in the 'southern' end of the molecule, far away from the copper site and the binding sites for the redox partners of azurin. By excluding the 'southern' region from the subsequent analysis, correlated motions were identified in the hydrophobic patch that surrounds the protein active site. The simulation results are in excellent agreement with recent NMR data on azurin in solution [A. V. Zhuravleva, D. M. Korzhnev, E. Kupce, A. S. Arseniev, M. Billeter, V. Y. Orekhov, Gated electron transfers and electron pathways in azurin: a NMR dynamic study at multiple fields and temperatures, J. Mol. Biol. 342 (2004) 1599-1611] and suggest a rationale for cooperative displacements of protein residues that are thought to be critical for the electron transfer process. A number of other structural and dynamic features of azurin are discussed in the context of the blue copper protein family and an explanation is proposed to account for the variability/conservation of some regions in the cupredoxins.     

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.     

A molecular dynamics simulation study of the solvent isotope effect on copper plastocyanin

The effect of heavy water on the structure and dynamics of copper plastocyanin as well as on some aspects of the solvent dynamics at the protein-solvent interfacial region have been investigated by molecular dynamics simulation. The simulated system has been analyzed in terms of the atomic root mean square deviation and fluctuations, intraprotein H-bond pattern, dynamical cross-correlation map and the results have been compared with those previously obtained for plastocyanin in H2O (Ciocchetti et al. Biophys. Chem. 69 (1997), 185-198). The simulated plastocyanin structure in the two solvents, averaging 1 ns, is very similar along the beta-structure regions, while the most significant differences are registered, analogous to the turns and the regions likely involved in the electron transfer pathway. Moreover, plastocyanin in D2O shows an increase in the number of both the intraprotein H-bonds and the residues involved in correlated motions. An analysis of the protein-solvent coupling evidenced that D2O makes the H-bond formation more difficult with the solvent molecules for positively charged and polar residues, while an opposite trend is observed for negatively charged residues. On the other hand, the frequency of exchange of the solvent molecules involved in the protein-solvent H-bond formation is significantly depressed in D2O. The results are discussed also in connection with protein functionality and briefly with some experimental results connected with the thermostability of proteins in D2O.     

An analysis of the influence of protein intrinsic dynamical properties on its thermal unfolding behavior

The influence of the protein topology-encoded dynamical properties on its thermal unfolding motions was studied in the present work. The intrinsic dynamics of protein topology was obtained by the anisotropic network model (ANM). The ANM has been largely used to investigate protein collective functional motions, but it is not well elucidated if this model can also reveal the preferred large-scale motions during protein unfolding. A small protein barnase is used as a typical case study to explore the relationship between protein topology-encoded dynamics and its unfolding motions. Three thermal unfolding simulations at 500 K were performed for barnase and the entire unfolding trajectories were sampled and partitioned into several windows. For each window, the preferred unfolding motions were investigated by essential dynamics analysis, and then associated with the intrinsic dynamical properties of the starting conformation in this window, which is detected by ANM. The results show that only a few slow normal modes imposed by protein structure are sufficient to give a significant overlap with the preferred unfolding motions. Especially, the large amplitude unfolding movements, which imply that the protein jumps out of a local energy basin, can be well described by a single or several ANM slow modes. Besides the global motions, it is also found that the local residual fluctuations encoded in protein structure are highly correlated with those in the protein unfolding process. Furthermore, we also investigated the relationship between protein intrinsic flexibility and its unfolding events. The results show that the intrinsic flexible regions tend to unfold early. Several early unfolding events can be predicted by analysis of protein structural flexibility. These results imply that protein structure-encoded dynamical properties have significant influences on protein unfolding motions.     

Molecular dynamics simulations of hemoglobin A in different states and bound to DPG: effector-linked perturbation of tertiary conformations and HbA concerted dynamics

Recent functional studies reported on human adult hemoglobin (HbA) show that heterotropic effector-linked tertiary structural changes are primarily responsible for modulating the oxygen affinity of hemoglobin. We present the results of 6-ns molecular dynamics simulations performed to gain insights into the dynamical and structural details of these effector-linked tertiary changes. All-atom simulations were carried out on a series of models generated for T- and R-state HbA, and for 2,3-diphosphoglycerate-bound models. Cross-correlation analyses identify both intra- and intersubunit correlated motions that are perturbed by the presence of the effector. Principal components analysis was used to decompose the covariance matrix extracted from the simulations and reconstruct the trajectories along the principal coordinates representative of functionally important collective motions. It is found that HbA in both quaternary states exists as ensembles of tertiary conformations that introduce dynamic heterogeneity in the protein. 2,3-Diphosphoglycerate induces significant perturbations in the fluctuations of both HbA states that translate into the protein visiting different tertiary conformations within each quaternary state. The analysis reveals that the presence of the effector affects the most important components of HbA motions and that heterotropic effectors modify the overall dynamics of the quaternary equilibrium via tertiary changes occurring in regions where conserved functionally significant residues are located, namely in the loop regions between helices C and E, E and F, and F and G, and in concerted helix motions. The changes are not apparent when comparing the available x-ray crystal structures in the presence and absence of effector, but are striking when comparing the respective dynamic tertiary conformations of the R and T tetramers.     

Conformational dynamics of cytochrome c: correlation to hydrogen exchange.

We study the dynamical fluctuations of horse heart cytochrome c by molecular dynamics (MD) simulations in aqueous solution, at four temperatures: 300 K, 360 K, 430 K, and 550 K. Each simulation covers a production time of at least 1.5 nanoseconds (ns). The conformational dynamics of the system is analyzed in terms of collective motions that involve the whole protein, and local motions that involve the formation and breaking of intramolecular hydrogen bonds. The character of the MD trajectories can be described within the framework of rugged energy landscape dynamics. The MD trajectories sample multiple conformational minima, with basins in protein conformational space being sampled for a few hundred picoseconds. The trajectories of the system in configurational space can be described in terms of diffusion of a particle in real space with a waiting time distribution due to partial trapping in shallow minima. As a consequence of the hierarchical nature of the dynamics, the mean square displacement autocorrelation function, <|x(t) - x(0)|2>, exhibits a power law dependence on time, with an exponent of around 0.5 for times shorter than 100 ps, and an exponent of 1.75 for longer times. This power law behavior indicates that the system exhibits suppressed diffusion (sub-diffusion) in sampling of configurational space at time scales shorter than 100 ps, and enhanced (super-diffusion) at longer time scales. The multi-basin feature of the trajectories is present at all temperatures simulated. Structural changes associated with inter-basin displacements correspond to collective motions of the Omega loops and coiled regions and relative motions of the alpha-helices as rigid bodies. Similar motions may be involved in experimentally observed amide hydrogen exchange. However, some groups showing large correlated motions do not expose the amino hydrogens to the solvent. We show that large fluctuations are not necessarily correlated to hydrogen exchange. For example, regions of the proteins forming alpha helices and turns show significant fluctuations, but as rigid bodies, and the hydrogen bonds involved in the formation of these structures do not break in proportion to these fluctuations. Proteins 1999;36:175-191. Published 1999 Wiley-Liss, Inc.     

Insight derived from molecular dynamics simulation into dynamics and molecular motions of cuticle-degrading serine protease Ver112

Cuticle-degrading serine protease Ver112, which derived from a nematophagous fungus Lecanicillium psalliotae, has been exhibited to have high cuticle-degrading and nematicidal activities. We have performed molecular dynamics (MD) simulation based on the crystal structure of Ver112 to investigate its dynamic properties and large-scale concerted motions. The results indicate that the structural core of Ver112 shows a small fluctuation amplitude, whereas the substrate binding sites, and the regions close to and opposite the substrate binding sites experience significant conformational fluctuations. The large concerted motions obtained from essential dynamics (ED) analysis of MD trajectory can lead to open or close of the substrate binding sites, which are proposed to be linked to the functional properties of Ver112, such as substrate binding, orientation, catalytic, and release. The significant motion in the loop regions that is located opposite the binding sites are considered to play an important role in modulating the dynamics of the substrate binding sites. Furthermore, the bottom of free energy landscape (FEL) of Ver112 are rugged, which is mainly caused by the fluctuations of substrate binding regions and loops located opposite the binding site. In addition, the mechanism underlying the high flexibility and catalytic activity of Ver112 was also discussed. Our simulation study complements the biochemical and structural studies, and provides insight into the dynamics-function relationship of cuticle-degrading serine protease Ver112.     

Insights derived from molecular dynamics simulation into the molecular motions of serine protease proteinase K

Serine protease proteinase K, a member of the subtilisin family of enzymes, is of significant industrial, agricultural and biotechnological importance. Despite the wealth of structural information about proteinase K provided by static X-ray structures, a full understanding of the enzymatic mechanism requires further insight into the dynamic properties of this enzyme. Molecular dynamics simulations and essential dynamics (ED) analysis were performed to investigate the molecular motions in proteinase K. The results indicate that the internal core of proteinase K is relatively rigid, whereas the surface-exposed loops, most notably the substrate-binding regions, exhibit considerable conformational fluctuations. Further ED analysis reveals that the large concerted motions in the substrate-binding regions cause opening/closing of the substrate-binding pockets, thus supporting the proposed induced-fit mechanism of substrate binding. The distinct electrostatic/hydrogen-bonding interactions between Asp39 and His69 and between His69 and Ser224 within the catalytic triad lead to different thermal motions and orientations of these three catalytic residues, which can be related to their different functional roles in the catalytic process. Statistical analyses of the geometrical/functional properties as well as evolutionary conservation of the glycines in proteinase K-like proteins reveal that glycines may play an important role in determining the folding architecture and structural flexibility of this class of enzymes. Our simulation study complements the biochemical and structural studies and provides new insights into the dynamic structural basis of the functional properties of this class of enzymes.     

A molecular dynamics analysis of protein structural elements

The relation between protein secondary structure and internal motions was examined by using molecular dynamics to calculate positional fluctuations of individual helix, beta-sheet, and loop structural elements in free and substrate-bound hen egg-white lysozyme. The time development of the fluctuations revealed a general correspondence between structure and dynamics; the fluctuations of the helices and beta-sheets converged within the 101 psec period of the simulation and were lower than average in magnitude, while the fluctuations of the loop regions were not converged and were mostly larger than average in magnitude. Notable exceptions to this pattern occurred in the substrate-bound simulation. A loop region (residues 101-107) of the active site cleft had significantly reduced motion due to interactions with the substrate. Moreover, part of a loop and a 3(10) helix (residues of 67-88) not in contact with the substrate showed a marked increase in fluctuations. That these differences in dynamics of free and substrate-bound lysozyme did not result simply from sampling errors was established by an analysis of the variations in the fluctuations of the two halves of the 101 psec simulation of free lysozyme. Concerted transitions of four to five mainchain phi and psi angles between dihedral wells were shown to be responsible for large coordinate shifts in the loops. These transitions displaced six or fewer residues and took place either abruptly, in 1 psec or less, or with a diffusive character over 5-10 psec. Displacements of rigid secondary structures involved longer timescale motions in bound lysozyme; a 0.5 A rms change in the position of a helix occurred over the 55 psec simulation period. This helix reorientation within the protein appears to be a response to substrate binding. There was little correlation between the solvent accessible surface area and the dynamics of the different structural elements.     

Effects of cavity-forming mutations on the internal dynamics of azurin

The effects of two single-point cavity-forming mutations, F110S and I7S, on the internal dynamics of azurin from Pseudomonas aeruginosa were probed by the phosphorescence emission of Trp-48, deeply buried in the compact hydrophobic core of the macromolecule. Changes in flexibility of the protein matrix around the chromophore were monitored by the intrinsic phosphorescence lifetime (tau(0)) whereas more general effects on structural fluctuations were deduced from the phosphorescence acrylamide quenching rate constant (k(q)), which measures the diffusion of the solute through the protein fold. The results show a spectacular, 4-5 orders of magnitude, increase of k(q) emphasizing that large amplitude structural fluctuations permitting acrylamide migration to the protein core have been drastically enhanced in each azurin mutant. The large, 12-15 kcal/mol, decrease in the activation enthalpy associated to k(q) suggests that the rate enhancement is caused, rather than through a generalized increase of protein flexibility, by the elimination of an inner barrier to the diffusion process. According to tau(0) the chromophore environment is more fluid with I7S but strikingly more rigid with F110S, demonstrating that when internal cavities are formed local effects on the mobility at the mutation site are unpredictable. Both tau(0) and k(q) reveal a structure tightening role of bound Cd(2+) that correlates with the increase in stability from apo- to holo-azurin. While these alterations in internal dynamics of azurin do not seem to play a role on electron transfer through the central region, the enhanced migration of acrylamide emphasizes that cavities may be critical for the rapid diffusion of substrates to buried, solvent inaccessible sites of enzymes.     

Relating molecular flexibility to function: a case study of tubulin

Microtubules (MT), along with a variety of associated motor proteins, are involved in a range of cellular functions including vesicle movement, chromosome segregation, and cell motility. MTs are assemblies of heterodimeric proteins, alpha beta-tubulins, the structure of which has been determined by electron crystallography of zinc-induced, pacilitaxel-stabilized tubulin sheets. These data provide a basis for examining relationships between structural features and protein function. Here, we study the fluctuation dynamics of the tubulin dimer with the aim of elucidating its functional motions relevant to substrate binding, polymerization/depolymerization and MT assembly. A coarse-grained model, harmonically constrained according to the crystal structure, is used to explore the global dynamics of the dimer. Our results identify six regions of collective motion, comprised of structurally close but discontinuous sequence fragments, observed only in the dimeric form, dimerization being a prerequisite for domain identification. Boundaries between regions of collective motions appear to act as linkages, found primarily within secondary-structure elements that lack sequence conservation, but are located at minima in the fluctuation curve, at positions of hydrophobic residues. Residue fluctuations within these domains identify the most mobile regions as loops involved in recognition of the adjacent regions. The least mobile regions are associated with nucleotide binding sites where lethal mutations occur. The functional coupling of motions between and within regions identifies three global motions: torsional and wobbling movements, en bloc, between the alpha- and beta-tubulin monomers, and stretching longitudinally. Further analysis finds the antitumor drug pacilitaxel (TaxotereR) to reduce flexibility in the M loop of the beta-tubulin monomer; an effect that may contribute to tightening lateral interactions between protofilaments assembled into MTs. Our analysis provides insights into relationships between intramolecular tubulin movements of MT organization and function.     

Restrained and unrestrained molecular dynamics simulations in the NVT ensemble of alamethicin

Molecular dynamics simulations on the transmembrane antibiotic peptide alamethicin have been performed in the NVT ensemble (i.e., the number of particles N, the volume V, and the temperature T of the system are kept constant). Results on the structure and conformational flexibility of this molecule are discussed and compared with previous experimental CD, x-ray, nmr data and theoretical computations on fragments analogues. An extensive study of structural and dynamic properties from H-bonding pattern analysis is presented. Evidences for a largely alpha-helix structure with some extent of freedom in the C-terminal region are found. Further, a partition of the molecule into three regions on the base of structural features and dynamic behavior has been proposed, and the correlation among the motions of the three regions is described.     

MD simulation of a plastocyanin mutant adsorbed onto a gold surface

MD simulation of plastocyanin, an electron transfer protein, adsorbed onto a gold surface, has been performed for 10 ns. Starting from the crystallographic structure of a poplar plastocyanin mutant engineered with the insertion of a disulfide bridge, the protein has been anchored to a gold substrate modeled by a cluster of three layers in the Au<111> configuration. A number of significant structural and dynamical properties of the protein molecule, covalently bound through either one or two sulfur atoms to the gold surface, has been extracted and compared with those of the free protein. Attention has been paid to investigate the dynamical aspects putatively related to the electron transfer process. In particular, the cross-correlation function between specific active site vibrations and all the other protein atom motions and the principal component analysis have been calculated in order to put into evidence dynamical correlation of some functional relevance. The results are discussed also in connection with related experiments.     

Topological and dynamical properties of Azurin anchored to a gold substrate as investigated by molecular dynamics simulation

A classical molecular dynamics study of the electron transfer protein azurin, covalently bound to a gold substrate through its native disulphide group, is carried out at full hydration. With the aim of investigating the effects on the protein structure and dynamics as induced by the presence of an electric field, simulations are performed on neutral, positively and negatively charged substrates. A number of parameters, such as the average structure, the root mean square deviations and fluctuations, the intraprotein hydrogen bonds and solvent accessible surface of the protein, are monitored during 10 ns of run. The orientation, the height and the lateral size of the protein, with respect to the substrate are evaluated and compared with the experimental data obtained by scanning probe nanoscopies. The electron transfer properties between the copper redox center and the disulphide bridge bound to the substrate are investigated and briefly discussed.     

Differences in internal dynamics of actin under different structural states detected by neutron scattering

F-actin, a helical polymer formed by polymerization of the monomers (G-actin), plays crucial roles in various aspects of cell motility. Flexibility of F-actin has been suggested to be important for such a variety of functions. Understanding the flexibility of F-actin requires characterization of a hierarchy of dynamical properties, from internal dynamics of the actin monomers through domain motions within the monomers and relative motions between the monomers within F-actin to large-scale motions of F-actin as a whole. As a first step toward this ultimate purpose, we carried out elastic incoherent neutron scattering experiments on powders of F-actin and G-actin hydrated with D₂O and characterized the internal dynamics of F-actin and G-actin. Well established techniques and analysis enabled the extraction of mean-square displacements and their temperature dependence in F-actin and in G-actin. An effective force constant analysis with a model consisting of three energy states showed that two dynamical transitions occur at ∼150 K and ∼245 K, the former of which corresponds to the onset of anharmonic motions and the latter of which couples with the transition of hydration water. It is shown that behavior of the mean-square displacements is different between G-actin and F-actin, such that G-actin is “softer” than F-actin. The differences in the internal dynamics are detected for the first time between the different structural states (the monomeric state and the polymerized state). The different behavior observed is ascribed to the differences in dynamical heterogeneity between F-actin and G-actin. Based on structural data, the assignment of the differences observed in the two samples to dynamics of specific loop regions involved in the polymerization of G-actin into F-actin is proposed.     

Collective motions in proteins: a covariance analysis of atomic fluctuations in molecular dynamics and normal mode simulations.

A method is described for identifying collective motions in proteins from molecular dynamics trajectories or normal mode simulations. The method makes use of the covariances of atomic positional fluctuations. It is illustrated by an analysis of the bovine pancreatic trypsin inhibitor. Comparison of the covariance and cross-correlation matrices shows that the relative motions have many similar features in the different simulations. Many regions of the protein, especially regions of secondary structure, move in a correlated manner. Anharmonic effects, which are included in the molecular dynamics simulations but not in the normal analysis, are of some importance in determining the larger scale collective motions, but not the more local fluctuations. Comparisons of molecular dynamics simulations in the present and absence of solvent indicate that the environment is of significance for the long-range motions.