REACH Coarse-Grained Simulation of a Cellulose Fiber

A molecular level understanding of the structure, dynamics and mechanics of cellulose fibers can aid in understanding the recalcitrance of biomass to hydrolysis in cellulosic biofuel production. Here, a residue-scale REACH (Realistic Extension Algorithm via Covariance Hessian) coarse-grained force field was derived from all-atom molecular dynamics (MD) simulations of the crystalline Iβ cellulose fibril. REACH maps the atomistic covariance matrix onto coarse-grained elastic force constants. The REACH force field was found to reproduce the positional fluctuations and low-frequency vibrational spectra from the all-atom model, allowing elastic properties of the cellulose fibril to be characterized using the coarse-grained force field with a speedup of >20 relative to atomistic MD on systems of the same size. The calculated longitudinal/transversal Young's modulus and the velocity of sound are in agreement with experiment. The persistence length of a 36-chain cellulose microcrystal was estimated to be ～380 μm. Finally, the normal-mode analysis with the REACH force field suggests that intrinsic dynamics might facilitate the deconstruction of the cellulose fibril from the hydrophobic surface.     

Dimer Dynamics and Filament Organization of the Bacterial Cell Division Protein FtsA

FtsA is a bacterial actin homolog and one of the core proteins involved in cell division. While previous studies have demonstrated the capability of FtsA to polymerize, little is known about its polymerization state in vivo or if polymerization is necessary for FtsA function. Given that one function of FtsA is to tether FtsZ filaments to the membrane, in vivo polymerization of FtsA imposes geometric constraints and requires a specific polymer curvature direction. Here, we report a series of molecular dynamics simulations probing the structural dynamics of FtsA as a dimer and as a tetrameric single filament. We found that the FtsA polymer exhibits a preferred bending direction that would allow for its placement parallel with FtsZ polymers underneath the cytoplasmic membrane. We also identified key interfacial amino acids that mediate FtsA–FtsA interaction and propose that some amino acids play more critical roles than others. We performed in silico mutagenesis on FtsA and demonstrated that, while a moderate mutation at the polymerization interface does not significantly affect polymer properties such as bending direction and association strength, more drastic mutations change both features and could lead to non-functional FtsA.     

Computational Analysis of C-Reactive Protein for Assessment of Molecular Dynamics and Interaction Properties

Serum C-reactive protein (CRP) is used as a marker of inflammation in several diseases including autoimmune disease and cardiovascular disease. CRP, a member of the pentraxin family, is comprised of five identical subunits. CRP has diverse ligand-binding properties which depend upon different structural states of CRP. However, little is known about the molecular dynamics and interaction properties of CRP. In this study, we used SAPS, SCRATCH protein predictor, PDBsum, ConSurf, ProtScale, Drawhca, ASAView, SCide and SRide server and performed comprehensive analyses of molecular dynamics, protein–protein and residue–residue interactions of CRP. We used 1GNH.pdb file for the crystal structure of human CRP which generated two pentamers (ABCDE and FGHIJ). The number of residues involved in residue–residue interactions between A–B, B–C, C–D, D–E, F–G, G–H, H–I, I–J, A–E and F–J subunits were 12, 11, 10, 11, 12, 11, 10, 11, 10 and 10, respectively. Fifteen antiparallel β sheets were involved in β-sheet topology, and five β hairpins were involved in forming the secondary structure. Analysis of hydrophobic segment distribution revealed deviations in surface hydrophobicity at different cavities present in CRP. Approximately 33 % of all residues were involved in the stabilization centers. We show that the bioinformatics tools can provide a rapid method to predict molecular dynamics and interaction properties of CRP. Our prediction of molecular dynamics and interaction properties of CRP combined with the modeling data based on the known 3D structure of CRP is helpful in designing stable forms of CRP mutants for structure–function studies of CRP and may facilitate in silico drug design for therapeutic targeting of CRP.     

Structure, Dynamics, and RNA Interaction Analysis of the Human SBDS Protein

Shwachman-Bodian-Diamond syndrome is an autosomal recessive genetic syndrome with pleiotropic phenotypes, including pancreatic deficiencies, bone marrow dysfunctions with increased risk of myelodysplasia or leukemia, and skeletal abnormalities. This syndrome has been associated with mutations in the SBDS gene, which encodes a conserved protein showing orthologs in Archaea and eukaryotes. The Shwachman-Bodian-Diamond syndrome pleiotropic phenotypes may be an indication of different cell type requirements for a fully functional SBDS protein. RNA-binding activity has been predicted for archaeal and yeast SBDS orthologs, with the latter also being implicated in ribosome biogenesis. However, full-length SBDS orthologs function in a species-specific manner, indicating that the knowledge obtained from model systems may be of limited use in understanding major unresolved issues regarding SBDS function, namely, the effect of mutations in human SBDS on its biochemical function and the specificity of RNA interaction. We determined the solution structure and backbone dynamics of the human SBDS protein and describe its RNA binding site using NMR spectroscopy. Similarly to the crystal structures of Archaea, the overall structure of human SBDS comprises three well-folded domains. However, significant conformational exchange was observed in NMR dynamics experiments for the flexible linker between the N-terminal domain and the central domain, and these experiments also reflect the relative motions of the domains. RNA titrations monitored by heteronuclear correlation experiments and chemical shift mapping analysis identified a classic RNA binding site at the N-terminal FYSH (fungal, Yhr087wp, Shwachman) domain that concentrates most of the mutations described for the human SBDS.     

Coarse-Grained Description of Protein Internal Dynamics: An Optimal Strategy for Decomposing Proteins in Rigid Subunits

The possibility of accurately describing the internal dynamics of proteins, in terms of movements of a few approximately-rigid subparts, is an appealing biophysical problem with important implications for the analysis and interpretation of data from experiments or numerical simulations. The problem is tackled here by means of a novel variational approach that exploits information about equilibrium fluctuations of interresidues distances, provided, e.g., by atomistic molecular dynamics simulations or coarse-grained models. No contiguity in primary sequence or in space is enforced a priori for amino acids grouped in the same rigid unit. The identification of the rigid protein moduli, or dynamical domains, provides valuable insight into functionally oriented aspects of protein internal dynamics. To illustrate this point, we first discuss the decomposition of adenylate kinase and HIV-1 protease and then extend the investigation to several representatives of the hydrolase enzymatic class. The known catalytic site of these enzymes is found to be preferentially located close to the boundary separating the two primary dynamical subdomains.     

Conformational Dynamics of Protein Side Chains and Enzyme-Substrate Interaction

Abstract

The electric interaction of charged and polar protein side chains with dipole moments of the bonds of substrate molecules bound in enzyme active sites is considered. The conformational motion of the side chains leads to the fact that the electric interaction, besides a constant (electrostatic) part, contains a fluctuating one, which is a random force (noise) exerted on the substrate molecule. On the time scale of enzyme turnovers this noise can be considered as the white one with good approximation. The noise is external, and the explicit expression for its intensity has been obtained. The possible functional role of the noise as an activating factor in enzyme catalysis is discussed.     

Sliding of Proteins Non-specifically Bound to DNA: Brownian Dynamics Studies with Coarse-Grained Protein and DNA Models

DNA binding proteins efficiently search for their cognitive sites on long genomic DNA by combining 3D diffusion and 1D diffusion (sliding) along the DNA. Recent experimental results and theoretical analyses revealed that the proteins show a rotation-coupled sliding along DNA helical pitch. Here, we performed Brownian dynamics simulations using newly developed coarse-grained protein and DNA models for evaluating how hydrodynamic interactions between the protein and DNA molecules, binding affinity of the protein to DNA, and DNA fluctuations affect the one dimensional diffusion of the protein on the DNA. Our results indicate that intermolecular hydrodynamic interactions reduce 1D diffusivity by 30%. On the other hand, structural fluctuations of DNA give rise to steric collisions between the CG-proteins and DNA, resulting in faster 1D sliding of the protein. Proteins with low binding affinities consistent with experimental estimates of non-specific DNA binding show hopping along the CG-DNA. This hopping significantly increases sliding speed. These simulation studies provide additional insights into the mechanism of how DNA binding proteins find their target sites on the genome.     

Protein Surface Dynamics: Interaction with Water and Small Solutes

Previous time resolved measurements had indicated that protons could propagate on the surface of a protein, or a membrane, by a special mechanism that enhances the shuttle of the proton towards a specific site [1]. It was proposed that a proper location of residues on the surface contributes to the proton shuttling function. In the present study, this notion was further investigated using molecular dynamics, with only the mobile charge replaced by Na⁺ and Cl⁻ ions. A molecular dynamics simulation of a small globular protein (the S6 of the bacterial ribosome) was carried out in the presence of explicit water molecules and four pairs of Na⁺ and Cl⁻ ions. A 10 ns simulation indicated that the ions and the protein's surface were in equilibrium, with rapid passage of the ions between the protein's surface and the bulk. Yet it was noted that, close to some domains, the ions extended their duration near the surface, suggesting that the local electrostatic potential prevented them from diffusing to the bulk. During the time frame in which the ions were detained next to the surface, they could rapidly shuttle between various attractor sites located under the electrostatic umbrella. Statistical analysis of molecular dynamics and electrostatic potential/entropy consideration indicated that the detainment state is an energetic compromise between attractive forces and entropy of dilution. The similarity between the motion of free ions next to a protein and the proton transfer on the protein's surface are discussed.     

Coarse-Grained Brownian Dynamics Simulations of the 10-23 DNAzyme

Deoxyribozymes (DNAzymes) are single-stranded DNA that catalyze nucleic acid biochemistry. Although a number of DNAzymes have been discovered by in vitro selection, the relationship between their tertiary structure and function remains unknown. We focus here on the well-studied 10-23 DNAzyme, which cleaves mRNA with a catalytic efficiency approaching that of RNase A. Using coarse-grained Brownian dynamics simulations, we find that the DNAzyme bends its substrate away from the cleavage point, exposing the reactive site and buckling the DNAzyme catalytic core. This hypothesized transition state provides microscopic insights into experimental observations concerning the size of the DNAzyme/substrate complex, the impact of the recognition arm length, and the sensitivity of the enzymatic activity to point mutations of the catalytic core. Upon cleaving the pertinent backbone bond in the substrate, we find that the catalytic core of the DNAzyme unwinds and the overall complex rapidly extends, in agreement with experiments on the related 8-17 DNAzyme. The results presented here provide a starting point for interpreting experimental data on DNAzyme kinetics, as well as developing more detailed simulation models. The results also demonstrate the limitations of using a simple physical model to understand the role of point mutations.     

Conformational Dynamics of Subtilisin-Chymotrypsin Inhibitor 2 Complex by Coarse-Grained Simulations

Abstract

An off-lattice dynamic Monte Carlo (MC) method is used to investigate the conformational dynamics of chymotrypsin inhibitor 2 (CI2) and subtilisin in both free and complex forms over two time windows, referring to short and long time scales. The conformational dynamics of backbone bonds analysed from several independent trajectories reveal that: Both the inhibitor and the enzyme are restricted in their bond rotations, excluding a few bonds, upon binding; the effect being greatest for the loop regions, and for the inhibitor. A cooperativity in the near-neighbor bond rotations are observed on both time scales, whereas the cooperative rotations of the bonds far along the sequence appear only in the long time window, and the latter time window is where most of the interactions between the inhibitor and the enzyme are observed. Upon binding, the cooperatively rotating parts of the inhibitor and the enzyme are readjusted compared to their free forms, and new correlations appear. The binding loop, although it is the closest contact region, is not the only part of the inhibitor involved in the interactions with the enzyme. Loops 3 and 8 and the helices F and G in bound enzyme and the binding loop of the inhibitor contribute at the most to the collective motions of whole structure on the slow time scale and are apparently important for enzyme-inhibitor interactions and function. The results in general provide evidence for the contribution of the loops with cooperative motions to the extensive communication network of the complex.     

A Functional DnaK Dimer Is Essential for the Efficient Interaction with Hsp40 Heat Shock Protein

Highly conserved molecular chaperone Hsp70 heat shock proteins play a key role in maintaining protein homeostasis (proteostasis). DnaK, a major Hsp70 in Escherichia coli, has been widely used as a paradigm for studying Hsp70s. In the absence of ATP, purified DnaK forms low-ordered oligomer, whereas ATP binding shifts the equilibrium toward the monomer. Recently, we solved the crystal structure of DnaK in complex with ATP. There are two molecules of DnaK-ATP in the asymmetric unit. Interestingly, the interfaces between the two molecules of DnaK are large with good surface complementarity, suggesting functional importance of this crystallographic dimer. Biochemical analyses of DnaK protein supported the formation of dimer in solution. Furthermore, our cross-linking experiment based on the DnaK-ATP structure confirmed that DnaK forms specific dimer in an ATP-dependent manner. To understand the physiological function of the dimer, we mutated five residues on the dimer interface. Four mutations, R56A, T301A, N537A, and D540A, resulted in loss of chaperone activity and compromised the formation of dimer, indicating the functional importance of the dimer. Surprisingly, neither the intrinsic biochemical activities, the ATP-induced allosteric coupling, nor GrpE co-chaperone interaction is affected appreciably in all of the mutations except for R56A. Unexpectedly, the interaction with co-chaperone Hsp40 is significantly compromised. In summary, this study suggests that DnaK forms a transient dimer upon ATP binding, and this dimer is essential for the efficient interaction of DnaK with Hsp40.     

Coarse Grained Normal Mode Analysis vs. Refined Gaussian Network Model for Protein Residue-Level Structural Fluctuations

We investigate several approaches to coarse grained normal mode analysis on protein residual-level structural fluctuations by choosing different ways of representing the residues and the forces among them. Single-atom representations using the backbone atoms C _α , C, N, and C _β are considered. Combinations of some of these atoms are also tested. The force constants between the representative atoms are extracted from the Hessian matrix of the energy function and served as the force constants between the corresponding residues. The residue mean-square-fluctuations and their correlations with the experimental B-factors are calculated for a large set of proteins. The results are compared with all-atom normal mode analysis and the residue-level Gaussian Network Model. The coarse-grained methods perform more efficiently than all-atom normal mode analysis, while their B-factor correlations are also higher. Their B-factor correlations are comparable with those estimated by the Gaussian Network Model and in many cases better. The extracted force constants are surveyed for different pairs of residues with different numbers of separation residues in sequence. The statistical averages are used to build a refined Gaussian Network Model, which is able to predict residue-level structural fluctuations significantly better than the conventional Gaussian Network Model in many test cases.     

Simple Topological Features Reflect Dynamics and Modularity in Protein Interaction Networks

The availability of large-scale protein-protein interaction networks for numerous organisms provides an opportunity to comprehensively analyze whether simple properties of proteins are predictive of the roles they play in the functional organization of the cell. We begin by re-examining an influential but controversial characterization of the dynamic modularity of the S. cerevisiae interactome that incorporated gene expression data into network analysis. We analyse the protein-protein interaction networks of five organisms, S. cerevisiae, H. sapiens, D. melanogaster, A. thaliana, and E. coli, and confirm significant and consistent functional and structural differences between hub proteins that are co-expressed with their interacting partners and those that are not, and support the view that the former tend to be intramodular whereas the latter tend to be intermodular. However, we also demonstrate that in each of these organisms, simple topological measures are significantly correlated with the average co-expression of a hub with its partners, independent of any classification, and therefore also reflect protein intra- and inter- modularity. Further, cross-interactomic analysis demonstrates that these simple topological characteristics of hub proteins tend to be conserved across organisms. Overall, we give evidence that purely topological features of static interaction networks reflect aspects of the dynamics and modularity of interactomes as well as previous measures incorporating expression data, and are a powerful means for understanding the dynamic roles of hubs in interactomes.     

The Role of Electrostatic Interaction in Triggering the Unraveling of Stable Helix 1 in Normal Prion Protein. A Molecular Dynamics Simulation Investigation

Abstract

The conversion of normal prion protein (PrP^C) into scrapie isoform (PrP^Sc) is a key event in the pathogenesis of prion diseases. However, the conversion mechanism has given rise to much controversy. For instance, there is much debate on the behavior of helix 1 (H1) in the conversion. A series of experiments demonstrated that H1 in isolated state was very stable under a variety of conditions. But, other experiments indicated that helices 2 and 3 rather than H1 were retained in PrP^Sc. In this paper, molecular dynamics (MD) simulation is employed to investigate the dynamic behavior of H1. It is revealed that although the helix 1 of Human PrP^C (HuPrP^C) is very stable in the isolated state, it becomes unstable when incorporated into native HuPrP^C, which likely results from the long-range electrostatic interaction between Asp147 and Arg208 located in the helices 1 and 3, respectively. This explanation is supported by experimental evaluation and MD simulation on D147N mutant of HuPrP^C that the mutant becomes a little more stable than the wild type HuPrP^C. This finding not only help to reconcile the existing debate on the role of helix 1 in the PrP^C→PrP^Sc transition, but also reveals a possible mechanism for triggering the PrP^C→PrP^Sc conversion.     

Role of Protein Interactions in Defining HIV-1 Viral Capsid Shape and Stability: A Coarse-Grained Analysis

Coarse-grained models of the HIV-1 CA dimer are constructed based on all-atom molecular dynamics simulations. Coarse-grained representations of the capsid shell, which is composed of ∼1500 copies of CA proteins, are constructed and their stability is examined. A key interaction between carboxyl and hexameric amino terminal domains is shown to generate the curvature of the capsid shell. It is demonstrated that variation of the strength of this interaction for different subunits in the lattice can cause formation of asymmetric, conical-shaped closed capsid shells, and it is proposed that variations, in the structure of the additional carboxyl-amino terminal binding interface during self-assembly, are important aspects of capsid cone formation. These results are in agreement with recent structural studies of the capsid hexamer subunit, which suggest that variability in the binding interface is a cause of the differences in subunit environments that exist in a conical structure.     

Anharmonic Normal Mode Analysis of Elastic Network Model Improves the Modeling of Atomic Fluctuations in Protein Crystal Structures

Protein conformational dynamics, despite its significant anharmonicity, has been widely explored by normal mode analysis (NMA) based on atomic or coarse-grained potential functions. To account for the anharmonic aspects of protein dynamics, this study proposes, and has performed, an anharmonic NMA (ANMA) based on the C_α-only elastic network models, which assume elastic interactions between pairs of residues whose C_α atoms or heavy atoms are within a cutoff distance. The key step of ANMA is to sample an anharmonic potential function along the directions of eigenvectors of the lowest normal modes to determine the mean-squared fluctuations along these directions. ANMA was evaluated based on the modeling of anisotropic displacement parameters (ADPs) from a list of 83 high-resolution protein crystal structures. Significant improvement was found in the modeling of ADPs by ANMA compared with standard NMA. Further improvement in the modeling of ADPs is attained if the interactions between a protein and its crystalline environment are taken into account. In addition, this study has determined the optimal cutoff distances for ADP modeling based on elastic network models, and these agree well with the peaks of the statistical distributions of distances between C_α atoms or heavy atoms derived from a large set of protein crystal structures.     

蛋白质相互作用网络进化分析研究进展

近年来,随着高通量实验技术的发展和广泛应用,越来越多可利用的蛋白质相互作用网络数据开始出现．这些数据为进化研究提供了新的视角．从蛋白质、蛋白质相互作用、模体、模块直到整个网络五个层次,综述了近年来蛋白质相互作用网络进化研究领域的主要进展,侧重于探讨蛋白质相互作用、模体、模块直到整个网络对蛋白质进化的约束作用,以及蛋白质相互作用网络不同于随机网络特性的起源和进化等问题．总结了前人工作给学术界的启示,探讨了该领域未来可能的发展方向．