Prediction of protein conformational freedom from distance constraints

A method is presented that generates random protein structures that fulfil a set of upper and lower interatomic distance limits. These limits depend on distances measured in experimental structures and the strength of the interatomic interaction. Structural differences between generated structures are similar to those obtained from experiment and from MD simulation. Although detailed aspects of dynamical mechanisms are not covered and the extent of variations are only estimated in a relative sense, applications to an IgG-binding domain, an SH3 binding domain, HPr, calmodulin, and lysozyme are presented which illustrate the use of the method as a fast and simple way to predict structural variability in proteins. The method may be used to support the design of mutants, when structural fluctuations for a large number of mutants are to be screened. The results suggest that motional freedom in proteins is ruled largely by a set of simple geometric constraints. Proteins 29:240–251, 1997. © 1997 Wiley-Liss, Inc.     

海三棱蔗草种群的密度与生物量动态

本文研究了上海市南汇县东海农场海堤外侧滩涂上海三棱藨草种群的密度动态,高度生长动态、生物量动态以及它们之间及其与环境之间的相互关系。研究结果表明：在环境条件相对稳定的地带A和B内,海三棱藨草种群的高度、高度生长和生物量在生长期内符合Logisfic增长。种群生物量动态与密度动态可分为3个阶段,其中阶段Ⅱ符合Yoda等提出的-3/2自疏定律。地带B为海三棱藨草种群生长的最适地带。地带C内生境条件极不稳定,种群的数量动态变化亦相当剧烈。在不同环境条件下,密度制约因素和非密度制约因素对种群数量动态的相对作用是不同的。在环境条件较稳定的生境中(地带A和B),密度制约因素是决定种群数量动态的主要因素;在环境条件变化剧烈的生境中(地带C),非密度制约因素是决定种群数量动态的主要因素。     

Role of methyl groups in dynamics and evolution of biomolecules

Recent studies have discovered strong differences between the dynamics of nucleic acids (RNA and DNA) and proteins, especially at low hydration and low temperatures. This difference is caused primarily by dynamics of methyl groups that are abundant in proteins, but are absent or very rare in RNA and DNA. In this paper, we present a hypothesis regarding the role of methyl groups as intrinsic plasticizers in proteins and their evolutionary selection to facilitate protein dynamics and activity. We demonstrate the profound effect methyl groups have on protein dynamics relative to nucleic acid dynamics, and note the apparent correlation of methyl group content in protein classes and their need for molecular flexibility. Moreover, we note the fastest methyl groups of some enzymes appear around dynamical centers such as hinges or active sites. Methyl groups are also of tremendous importance from a hydrophobicity/folding/entropy perspective. These significant roles, however, complement our hypothesis rather than preclude the recognition of methyl groups in the dynamics and evolution of biomolecules.     

海三棱藨草种群的密度与生物量动态

本文研究了上海市南汇县东海农场海堤外侧滩涂上海三棱藨草种群的密度动态、高度生长动态、生物量动态以及它们之间及其与环境之间的相互关系。研究结果表明:在环境条件相对稳定的地带A和B内,海三棱藨草种群的高度、高度生长和生物量在生长期内符合Logisfis增长。种群生物量动态与密度动态可分为3个阶段,其中阶段Ⅱ符合Yoda等提出的-3/2自疏定律。地带B为海三棱藨草种群生长的最适地带。地带C内生境条件极不稳定,种群的数量动态变化亦相当剧烈。在不同环境条件下,密度制约因素和非密度制约因素对种群数量动态的相对作用是不同的。在环境条件较稳定的生境中(地带A和B),密度制约因素是决定种群数量动态的主要因素;在环境条件变化剧烈的生境中(地带C),非密度制约因素是决定种群数量动态的主要因素。     

Spatial and temporal dynamics of the cardiac mitochondrial proteome

Mitochondrial proteins alter in their composition and quantity drastically through time and space in correspondence to changing energy demands and cellular signaling events. The integrity and permutations of this dynamism are increasingly recognized to impact the functions of the cardiac proteome in health and disease. This article provides an overview on recent advances in defining the spatial and temporal dynamics of mitochondrial proteins in the heart. Proteomics techniques to characterize dynamics on a proteome scale are reviewed and the physiological consequences of altered mitochondrial protein dynamics are discussed. Lastly, we offer our perspectives on the unmet challenges in translating mitochondrial dynamics markers into the clinic.     

Multiple time step diffusive Langevin dynamics for proteins

We present an algorithm for simulating the long time scale dynamics of proteins and other macromolecules. Our method applies the concept of multiple time step integration to the diffusive Langevin equation, in which short time scale dynamics are replaced by friction and noise. The macromolecular force field is represented at atomic resolution. Slow motions are modeled by constrained Langevin dynamics with very large time steps, while faster degrees of freedom are kept in local thermal equilibrium. In the limit of a sufficiently large molecule, our algorithm is shown to reduce the CPU time required by two orders of magnitude. We test the algorithm on two systems, alanine dipeptide and bovine pancreatic trypsin inhibitor (BPTI), and find that it accurately calculates a variety of equilibrium and dynamical properties. In the case of BPTI, the CPU time required is reduced by nearly a factor of 60 compared to a conventional, unconstrained Langevin simulation using the same force field. Proteins 30:215–227, 1998. © 1998 Wiley-Liss, Inc.     

Cooperative Thermal Denaturation of the Assembly Origin Region of TMV RNA

Abstract

The assembly origin (AO) region of the tobacco mosaic virus RNA melts in an unusually narrow(2.5°C) temperature range. In an 0.01 M phosphate buffer the melting temperature of AO was found to be 41.5°C. This value corresponds to the regions with the most stable secondary/tertiary structure of the whole TMV RNA molecule. It is assumed that the AO region has a specific tertiary structure, which is maintained by the long-range interactions as well as by interactions of the pseudoknot type.     

Nudged elastic band calculation of minimal energy paths for the conformational change of a GG non-canonical pair

The nudged elastic band (NEB) technique has been implemented in AMBER to calculate low-energy paths for conformational changes. A novel simulated annealing protocol that does not require an initial hypothesis for the path is used to sample low-energy paths. This was used to study the conformational change of an RNA cis Watson-Crick/Hoogsteen GG non-canonical pair, with one G syn around the glycosidic bond and the other anti. A previous solution structure, determined by NMR-constrained modeling, demonstrated that the GG pairs change from (syn)G-(anti)G to (anti)G-(syn)G in the context of duplex r(GCAGGCGUGC) on the millisecond timescale. The set of low-energy paths found by NEB show that each G flips independently around the glycosidic bond, with the anti G flipping to syn first. Guanine bases flip without opening adjacent base-pairs by protruding into the major groove, accommodated by a transient change by the ribose to C2'-exo sugar pucker. Hydrogen bonds between bases and the backbone, which lower the energetic barrier to flipping, are observed along the path. The results show the plasticity of RNA base-pairs in helices, which is important for biological processes, including mismatch repair, protein recognition, and translation. The modeling of the GG conformational change also demonstrates that NEB can be used to discover non-trivial paths for macromolecules and therefore NEB can be used as an exploratory method for predicting putative conformational change paths.     

Molecular dynamics studies of alanine racemase: a structural model for drug design

Alanine racemase (AlaR) is a bacterial enzyme that catalyzes the interconversion of L- and D-alanine, which is an essential constituent of the peptidoglycan layer of the bacterial cell wall and requires pyridoxal 5'-phosphate (PLP) as a cofactor. The enzyme is universal to bacteria, including mycobacteria, making it an attractive target for drug design. To investigate the effects of flexibility on the binding modes of the substrate and an inhibitor and to analyze how the active site is affected by the presence of the substrate versus inhibitor, a molecular dynamics simulation on the full AlaR dimer from Bacillus stearothermophilus (pdb code: 1SFT) with a D-alanine molecule in one active site and the noncovalent inhibitor, propionate, in the second site has been carried out. Within the time scale of the simulation, we show that the active site becomes more stabilized in the presence of substrate versus inhibitor. The results of this simulation are in agreement with the proposed mechanism of alanine racemase reaction in which the substrate carboxyl group directly participates in the catalysis by acting cooperatively with Tyr 265' and Lys 39. A structural water molecule in contact with both substrate and inhibitor (i.e., in both active sites) and bridging residues in both active sites was identified. It shows a remarkably low mobility and does not exchange with bulk water. This water molecule can be taken into account for the design of specific AlaR inhibitors by either utilizing it as a bridging group or displacing it with an inhibitor atom. The results presented here provide insights into the dynamics of the alanine racemase in the presence of substrate/inhibitor, which will be used for the rational design of novel inhibitors.     

High pressure reveals structural determinants for globin hexacoordination: Neuroglobin and myoglobin cases

The influence of pressure on the equilibrium between five‐(5c) and six‐coordination (6c) forms in neuroglobin (Ngb) and myoglobin (Mb) has been examined by means of molecular dynamics (MD) simulations at normal and high pressure. The results show that the main effect of high pressure is to reduce the protein mobility without altering the structure in a significant manner. Moreover, our data suggest that the equilibrium between 5c and 6c states in globins is largely controlled by the structure and dynamics of the C‐D region. Finally, in agreement with the available experimental data, the free energy profiles obtained from steered MD for both proteins indicate that high pressure enhances hexacoordination. In Ngb, the shift in equilibrium is mainly related to an increase in the 6c→5c transition barrier, whereas in Mb such a shift is primarily due to a destabilization of the 5c state. Proteins 2009. © 2008 Wiley‐Liss, Inc.     

Investigation of flap flexibility of β-secretase using molecular dynamic simulations

Flap motif and its dynamics were extensively reported in aspartate proteases, e.g. HIV proteases and plasmepsins. Herein, we report the first account of flap dynamics amongst different conformations of β-secretase using molecular dynamics simulation. Various parameters were proposed and a selected few were picked which could appropriately describe the flap motion. Three systems were studied, namely Free (BACE_Free) and two ligand-bound conformations, which belonged to space groups P6₁22 (BACE_Bound1) and C222₁ (BACE_Bound2), respectively and four parameters (distance between the flaps tip residue, Thr72 and Ser325, d₁; dihedral angle, ? (Thr72-Asp32-Asp228-Ser325); TriCα angles, θ₁ (Thr72-Asp32-Ser325), and θ₂ (Thr72-Asp228-Ser325)) were proposed to understand the change in dynamics of flap domain and the extent of flap opening and closing. Analysis of, θ₂, d₁, θ₁ and ? confirmed that the BACE_Free adopted semi-open, open and closed conformations with slight twisting during flap opening. However, BACE_Bound1 (P6122) showed an adaptation to open conformation due to lack of hydrogen bond interaction between the ligand and flap tip residue. A slight flap twisting, ? (lateral twisting) was observed for BACE_Bound1 during flap opening which correlates with the opening of BACE_Free. Contradictory to the BACE_Bound1, the BACE_Bound2 locked the flap in a closed conformation throughout the simulation due to formation of a stable hydrogen bond interaction between the flap tip residue and ligand. Analyses of all three systems highlight that d₁, θ₂ and ? can be precisely used to describe the extent of flap opening and closing concurrently with snapshots along the molecular dynamics trajectory across several conformations of β-secretase.     

Ligand unbinding pathways from the vitamin D receptor studied by molecular dynamics simulations

Molecular dynamics simulation techniques have been used to study the unbinding pathways of 1α,25-dihydroxyvitamin D₃ from the ligand-binding pocket of the vitamin D receptor (VDR). The pathways observed in a large number of relatively short (<200 ps) random acceleration molecular dynamics (RAMD) trajectories were found to be in fair agreement, both in terms of pathway locations and deduced relative preferences, compared to targeted molecular dynamics (TMD) and streered molecular dynamics simulations (SMD). However, the high-velocity ligand expulsions of RAMD tend to favor straight expulsion trajectories and the observed relative frequencies of different pathways were biased towards the probability of entering a particular exit channel. Simulations indicated that for VDR the unbinding pathway between the H1–H2 loop and the β-sheet between H5 and H6 is more favorable than the pathway located between the H1–H2 loop and H3. The latter pathway has been suggested to be the most likely unbinding path for thyroid hormone receptors (TRs) and a likely path for retinoic acid receptor. Ligand entry/exit through these two pathways would not require displacement of H12 from its agonistic position. Differences in the packing of the H1, H2, H3 and β-sheet region explain the changed relative preference of the two unbinding pathways in VDR and TRs. Based on the crystal structures of the ligand binding domains of class 2 nuclear receptors, whose members are VDR and TRs, this receptor class can be divided in two groups according to the packing of the H1, H2, H3 and β-sheet region. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Breaths,Twists, and Turns of Atomistic Nucleosomes

Gene regulation programs establish cellular identity and rely on dynamic changes in the structural packaging of genomic DNA. The DNA is packaged in chromatin, which is formed from arrays of nucleosomes displaying different degree of compaction and different lengths of inter-nucleosomal linker DNA. The nucleosome represents the repetitive unit of chromatin and is formed by wrapping 145–147 basepairs of DNA around an octamer of histone proteins. Each of the four histones is present twice and has a structured core and intrinsically disordered terminal tails. Chromatin dynamics are triggered by inter- and intra-nucleosome motions that are controlled by the DNA sequence, the interactions between the histone core and the DNA, and the conformations, positions, and DNA interactions of the histone tails. Understanding chromatin dynamics requires studying all these features at the highest possible resolution. For this, molecular dynamics simulations can be used as a powerful complement or alternative to experimental approaches, from which it is often very challenging to characterize the structural features and atomic interactions controlling nucleosome motions. Molecular dynamics simulations can be performed at different resolutions, by coarse graining the molecular system with varying levels of details. Here we review the successes and the remaining challenges of the application of atomic resolution simulations to study the structure and dynamics of nucleosomes and their complexes with interacting partners.     

Flap opening dynamics in HIV-1 protease explored with a coarse-grained model

We present a one-bead coarse-grained model that enables dynamical simulations of proteins on the time scale of tens of microseconds. The parameterization of the force field includes accurate conformational terms that allow for fast and reliable exploration of the configurational space. The model is applied to the dynamics of flap opening in HIV-1 protease. The experimental structure of the recently crystallized semi-open conformation of HIV-1 protease is well reproduced in the simulation, which supports the accuracy of our model. Thanks to very long simulations and extensive sampling of opening and closing events, we also investigate the thermodynamics and kinetics of the opening process. We have shown that the effect of the solvent slows down the dynamics to the experimentally observed time scales. The model is found to be reliable for application to substrate docking simulations, which are currently in progress.     

The role of disulfide bonds and N‐terminus in the structural properties of hepcidins: Insights from molecular dynamics simulations

The main purpose of this work is to analyse, by means of molecular dynamics (MD) simulations both structural and mechanical‐dynamical differences between Hepcidin‐20 and Hepcidin‐25 in both oxidized and reduced states in aqueous solution. Results indicate that the presence of disulfide bonds is essential, in both peptides, for maintaining their β‐hairpin motif. As a matter of fact, the lack of this intra‐peptide covalent interactions produces an almost immediate deviation from the oxidized, plausibly active, structure in both the systems. Interestingly, reduced Hepcidin‐25 turns out to be characterized by a highly fluctuating structure which is found to rapidly span a large number of configurations at equilibrium. On the other hand, loss of disulfide bonds in the shorter peptide, results in a more compact and relatively rigid double‐turn structure. Comparison of mechanical–dynamical properties and sidechains–sidechains interactions in oxidized Hepcidin‐20 and Hepcidin‐25 strongly suggest also the key role of N‐terminus in the aggregation tendency of Hepcidin‐25. © 2010 Wiley Periodicals, Inc. Biopolymers 93: 917–926, 2010.     

Emerging computational approaches for the study of protein allostery

Allosteric regulation of protein function is key in controlling cellular processes so its underlying mechanisms are of primary concern to research in areas spanning protein engineering and drug design. However, due to the complex nature of allosteric mechanisms, a clear and predictive understanding of the relationship between protein structure and allosteric function remains elusive. Well established experimental approaches are available to offer a limited degree of characterization of mechanical properties within proteins, but the analytical capabilities of computational methods are evolving rapidly in their ability to accurately define the subtle and concerted structural dynamics that comprise allostery. This review includes a brief overview of allostery in proteins and an exploration of relevant experimental methods. An explanation of the transition from experimental toward computational methods for allostery is discussed, followed by a review of existing and emerging methods.     

Variations on a theme by Debye and Waller: From simple crystals to proteins

Debye and Waller showed how to adjust scattering intensities in diffraction experiments for harmonic motions of atoms about an average, static reference configuration. However, many motions, particularly in biological molecules as compared to simple crystals, are far from harmonic. We show how, using a variety of simple anharmonic, multiconformational models, it is possible to construct a variety of Generalized Debye-Waller Factors, and understand their meaning. A central result for these cases is that, in principle, intensity factors cannot be obtained from true total mean square displacements of the atoms. We make the distinction between the intensity factors for unimodal quasiharmonic motions and those for the anharmonic, multimodal (valley hopping) motions. Only the former affect the conventional B factors. Proteins 29:153–160, 1997. Published 1997 Wiley-Liss, Inc.

1 This article is a US Government work and, as such, is in the public domain in the United States of America.