Exploring the helix-coil transition via all-atom equilibrium ensemble simulations

Modelling proteins: conformational sampling and reconstruction of folding kinetics

In the last decades biomolecular simulation has made tremendous inroads to help elucidate biomolecular processes in-silico. Despite enormous advances in molecular dynamics techniques and the available computational power, many problems involve long time scales and large-scale molecular rearrangements that are still difficult to sample adequately. In this review we therefore summarise recent efforts to fundamentally improve this situation by decoupling the sampling of the energy landscape from the description of the kinetics of the process. Recent years have seen the emergence of many advanced sampling techniques, which permit efficient characterisation of the relevant family of molecular conformations by dispensing with the details of the short-term kinetics of the process. Because these methods generate thermodynamic information at best, they must be complemented by techniques to reconstruct the kinetics of the process using the ensemble of relevant conformations. Here we review recent advances for both types of methods and discuss their perspectives to permit efficient and accurate modelling of large-scale conformational changes in biomolecules. This article is part of a Special Issue entitled: Protein Dynamics: Experimental and Computational Approaches.

SH3 domain-peptide binding energy calculations based on structural ensemble and multiple peptide templates

SH3 domains mediate signal transduction by recognizing short peptides. Understanding of the driving forces in peptide recognitions will help us to predict the binding specificity of the domain-peptide recognition and to understand the molecular interaction networks of cells. However, accurate calculation of the binding energy is a tough challenge. In this study, we propose three ideas for improving our ability to predict the binding energy between SH3 domains and peptides: (1) utilizing the structural ensembles sampled from a molecular dynamics simulation trajectory, (2) utilizing multiple peptide templates, and (3) optimizing the sequence-structure mapping. We tested these three ideas on ten previously studied SH3 domains for which SPOT analysis data were available. The results indicate that calculating binding energy using the structural ensemble was most effective, clearly increasing the prediction accuracy, while the second and third ideas tended to give better binding energy predictions. We applied our method to the five SH3 targets in DREAM4 Challenge and selected the best performing method.

Comparative biochemical and computational study of the role of naturally occurring mutations at Ambler positions 104 and 170 in GES β-lactamases

In GES-type β-lactamases, positions 104 and 170 are occupied by Glu or Lys and by Gly, Asn, or Ser, respectively. Previous studies have indicated an important role of these amino acids in the interaction with β-lactams, although their precise role, especially that of residue 104, remains uncertain. In this study, we constructed GES-1 (Glu104, Gly170), GES-2 (Glu104, Asn170), GES-5 (Glu104, Ser170), GES-6 (Lys104, Ser170), GES-7 (Lys104, Gly170), and GES-13 (Lys104, Asn170) by site-specific mutagenesis and compared their hydrolytic properties. Isogenic comparisons of β-lactam resistance levels conferred by these GES variants were also performed. Data indicated the following patterns: (i) Lys104-containing enzymes exhibited enhanced hydrolysis of oxyimino-cephalosporins and reduced efficiency against imipenem in relation to enzymes possessing Glu104, (ii) Asn170-containing enzymes showed reduced hydrolysis rates of penicillins and older cephalosporins, (iii) Ser170 enabled GES to hydrolyze cefoxitin efficiently, and (iv) Asn170 and Ser170 increased the carbapenemase character of GES enzymes but reduced their activity against ceftazidime. Molecular dynamic simulations of GES apoenzyme models, as well as construction of GES structures complexed with cefoxitin and an achiral ceftazidime-like boronic acid, provided insights into the catalytic behavior of the studied mutants. There were indications that an increased stability of the hydrogen bonding network of Glu166-Lys73-Ser70 and an altered positioning of Trp105 correlated with the substrate spectra, especially with acylation of GES by imipenem. Furthermore, likely effects of Ser170 on GES interactions with cefoxitin and of Lys104 on interactions with oxyimino-cephalosporins were revealed. Overall, the data unveiled the importance of residues 104 and 170 in the function of GES enzymes.

Globular structure of a human immunodeficiency virus-1 protease (1DIFA dimer) in an effective solvent medium by a Monte Carlo simulation

A coarse-grained model is used to study the structure and dynamics of a human immunodeficiency virus-1 protease (1DIFA dimer) consisting of 198 residues in an effective solvent medium on a cubic lattice by Monte Carlo simulations for a range of interaction strengths. Energy and mobility profiles of residues are found to depend on the interaction strength and exhibit remarkable segmental symmetries in two monomers. Lowest energy residues such as Arg(41) and Arg(140) (most electrostatic and polar) are not the least mobile; despite the higher energy, the hydrophobic residues (Ile, Leu, and Val) are least mobile and form the core by pinning down the local segments for the globular structure. Variations in the gyration radius (R(g)) and energy (E(c)) of the protein show nonmonotonic dependence on the interaction strength with the smallest R(g) around the largest value of E(c). Pinning of the conformations by the hydrophobic residues at high interaction strength seems to provide seed for the protein chain to collapse.

Nontemplated approach to tuning the spectral properties of cyanine-based fluorescent nanoGUMBOS

Template-free controlled aggregation and spectral properties in fluorescent organic nanoparticles (FONs) is highly desirable for various applications. Herein, we report a nontemplated method for controlling the aggregation in near-infrared (NIR) cyanine-based nanoparticles derived from a group of uniform materials based on organic salts (GUMBOS). Cationic heptamethine cyanine dye 1,1',3,3,3',3'-hexamethylindotricarbocyanine (HMT) was coupled with five different anions, viz., [NTf(2)(-)], [BETI(-)], [TFPB(-)], [AOT(-)], and [TFP4B(-)], by an ion-exchange method to obtain the respective GUMBOS. The nanoGUMBOS obtained via a reprecipitation method were primarily amorphous and spherical (30-100 nm) as suggested by selected area electron diffraction (SAED) and transmission electron microscopy (TEM). The formation of tunable self-assemblies within the nanoGUMBOS was characterized using absorption and fluorescence spectroscopy in conjunction with molecular dynamics simulations. Counterion-controlled spectral properties observed in the nanoGUMBOS were attributed to variations in J/H ratios with different anions. Association with the [AOT(-)] anion afforded predominant J aggregation enabling the highest fluorescence intensity, whereas [TFP4B(-)] disabled the fluorescence due to predominant H aggregation in the nanoparticles. Analyses of the stacking angle of the cations based on molecular dynamic simulation results in [HMT][NTf(2)], [HMT][BETI], and [HMT][AOT] dispersed in water and a visual analysis of the representative simulation snapshots also imply that the type of aggregation was controlled through the counterion associated with the dye cation.

Force field-dependent structural divergence revealed during long time simulations of Calbindin d9k

The structural and the dynamic features of the Calbindin (CaB) protein in its holo and apo states are compared using molecular dynamics simulations under nine different force fields (FFs) (G43a1, G53a6, Opls-AA, Amber94, Amber99, Amber99p, AmberGS, AmberGSs, and Amber99sb). The results show that most FFs reproduce reasonably well the majority of the experimentally derived features of the CaB protein. However, in several cases, there are significant differences in secondary structure properties, root mean square deviations (RMSDs), root mean square fluctuations (RMSFs), and S(2) order parameters among the various FFs. What is more, in certain cases, these parameters differed from the experimentally derived values. Some of these deviations became noticeable only after 50 ns. A comparison with experimental data indicates that, for CaB, the Amber94 shows overall best agreement with the measured values, whereas several others seem to deviate from both crystal and nuclear magnetic resonance data.

Metadynamics as a tool for mapping the conformational and free-energy space of peptides--the alanine dipeptide case study

There is a need for a fast, accurate, and reliable method of sampling the conformational space of peptides and proteins in order to obtain a balanced free-energy profile which can lead to our understanding of protein structure. We have utilized metadynamics for the conformational study of the solvated alanine dipeptide molecule, and our results show that the method has proven to be competent as a fast, robust, and reliable method for the conformation free-energy calculations of peptides in an explicit solvent, surpassing traditional methods such as umbrella sampling. We have also addressed the issue of the influence of different water models on the resulting free-energy profile in order to consistently decompose the setting of our simulation. All of the explicit water models for the simulation of biomolecules TIP3P, TIP4P, TIP4P/Ew, TIP5P, and SPCE have exhibited similar effects on the conformational preferences of alanine dipeptide with no significant differences. On the other hand, by comparing the potential energy surface in the gas phase and the free-energy surface in a water environment, we have shown that the interaction with water molecules is one of the most important structure-driving elements, with a great influence on the free-energy surface (FES) of the solvated peptide and the conformational preferences of the peptide backbone. All of the tested force fields (ff03, ff99SB, opls-aa, and charmm27) appreciably differ in the population of the individual conformers and the barriers between them. Significant divergence was found on both the potential energy surface (PES) as well as free-energy surface (FES) calculated by charmm27. We have therefore concluded that the differences originate dominantly from the parametrization of the peptide backbone in the given force field rather than from a noncovalent interaction with water molecules.

Theoretical study of molecular mechanism of binding TRAP220 coactivator to Retinoid X Receptor alpha, activated by 9-cis retinoic acid

Study on molecular mechanism of conformational reorientation of RXR-alpha ligand binding domain is presented. We employed CABS--a reduced model of protein dynamics to model folding pathways of binding 9-cis retinoic acid to apo-RXR molecule and TRAP220 peptide fragment to the holo form. Based on obtained results we also propose a sequential model of RXR activation by 9-cis retinoic acid and TRAP220 coactivator. Methodology presented here may be used for investigation of binding pathways of other NR/hormone/cofactor sets.

Enhanced sampling and applications in protein folding in explicit solvent

We report a single-copy tempering method for simulating large complex systems. In a generalized ensemble, the method uses runtime estimate of the thermal average energy computed from a novel integral identity to guide a continuous temperature-space random walk. We first validated the method in a two-dimensional Ising model and a Lennard-Jones liquid system. It was then applied to folding of three small proteins, trpzip2, trp-cage, and villin headpiece in explicit solvent. Within 0.5-1 microsecond, all three systems were reversibly folded into atomic accuracy: the alpha carbon root mean square deviations of the best folded conformations from the native states were 0.2, 0.4, and 0.4 A, for trpzip2, trp-cage, and villin headpiece, respectively.

Structures of apo and GTP-bound molybdenum cofactor biosynthesis protein MoaC from Thermus thermophilus HB8

The first step in the molybdenum cofactor (Moco) biosynthesis pathway involves the conversion of guanosine triphosphate (GTP) to precursor Z by two proteins (MoaA and MoaC). MoaA belongs to the S-adenosylmethionine-dependent radical enzyme superfamily and is believed to generate protein and/or substrate radicals by reductive cleavage of S-adenosylmethionine using an Fe-S cluster. MoaC has been suggested to catalyze the release of pyrophosphate and the formation of the cyclic phosphate of precursor Z. However, structural evidence showing the binding of a substrate-like molecule to MoaC is not available. Here, apo and GTP-bound crystal structures of MoaC from Thermus thermophilus HB8 are reported. Furthermore, isothermal titration calorimetry experiments have been carried out in order to obtain thermodynamic parameters for the protein-ligand interactions. In addition, molecular-dynamics (MD) simulations have been carried out on the protein-ligand complex of known structure and on models of relevant complexes for which X-ray structures are not available. The biophysical, structural and MD results reveal the residues that are involved in substrate binding and help in speculating upon a possible mechanism.

Sodium perchlorate effects on the helical stability of a mainly alanine peptide

Sodium perchlorate salt (NaClO(4)) is commonly used as an internal intensity standard in ultraviolet resonance Raman (UVRR) spectroscopy experiments. It is well known that NaClO(4) can have profound effects on peptide stability. The impact of NaClO(4) on protein stability in UVRR experiments has not yet been fully investigated. It is well known from experiment that protein stability is strongly affected by the solution composition (water, salts, osmolytes, etc.). Therefore, it is of the utmost importance to understand the physical basis on which the presence of salts and osmolytes in the solution impact protein structure and stability. The aim of this study is to investigate the effects of NaClO(4), on the helical stability of an alanine peptide in water. Based upon replica-exchange molecular dynamics data, it was found that NaClO(4) solution strongly stabilizes the helical state and that the number of pure helical conformations found at room temperature is greater than in pure water. A thorough investigation of the anion effects on the first and second solvation shells of the peptide, along with the Kirkwood-Buff theory for solutions, allows us to explain the physical mechanisms involved in the observed specific ion effects. A direct mechanism was found in which ClO(4)(-) ions are strongly attracted to the folded backbone.

Structural and dynamical effects induced by the anticancer drug topotecan on the human topoisomerase I - DNA complex

BACKGROUND

Human topoisomerase I catalyzes the relaxation of DNA supercoils in fundamental cell processes like transcription, replication and chromosomal segregation. It is the only target of the camptothecin family of anticancer drugs. Among these, topotecan has been used to treat lung and ovarian carcinoma for several years. Camptothecins reversibly binds to the covalent intermediate DNA-enzyme, stabilizing the cleavable complex and reducing the religation rate. The stalled complex then collides with the progression of the replication fork, producing lethal double strand DNA breaks and eventually cell death.

METHODOLOGY/PRINCIPAL FINDINGS

Long lasting molecular dynamics simulations of the DNA-topoisomerase I binary complex and of the DNA-topoisomerase-topotecan ternary complex have been performed and compared. The conformational space sampled by the binary complex is reduced by the presence of the drug, as observed by principal component and cluster analyses. This conformational restraint is mainly due to the reduced flexibility of residues 633-643 (the region connecting the linker to the core domain) that causes an overall mobility loss in the ternary complex linker domain. During the simulation, DNA/drug stacking interactions are fully maintained, and hydrogen bonds are maintained with the enzyme. Topotecan keeps the catalytic residue Lys532 far from the DNA, making it unable to participate to the religation reaction. Arg364 is observed to interact with both the B and E rings of topotecan with two stable direct hydrogen bonds. An interesting constrain exerted by the protein on the geometrical arrangement of topotecan is also observed.

CONCLUSIONS/SIGNIFICANCE

Atomistic-scale understanding of topotecan interactions with the DNA-enzyme complex is fundamental to the explaining of its poisonous effect and of the drug resistance observed in several single residue topoisomerase mutants. We observed significant alterations due to topotecan in both short-range interactions and long-range protein domain communications.

Crucial roles of the subnanosecond local dynamics of the flap tips in the global conformational changes of HIV-1 protease

To understand the underlying mechanisms of the open and closed conformational change of HIV-1 protease (HIV-1 PR) at multiple time scales, we performed serial fully unrestrained, extremely long time molecular dynamics simulations with an explicit solvent model. Spontaneous semiopen to closed conformational transition and inhibitor-collision-induced opening of the flaps were simulated in a real time scale. We found that the rapid, local subnanosecond fluctuations of the flap tips might be the mechanisms triggering the global open and close conformational transitions at the 100-ns time scale. The subnanosecond fluctuation is induced by the Phi-Psi rotations of the residues at the flap tips, mainly Psi of Gly49 and Phi of Ile50, disturbing the interactions between the two tips and then their stability. We further showed that the water molecule W301 is helpful for the stability of the PR-inhibitor complex by acting as a collision buffer for the dynamic interaction between flap tips and the inhibitor. These results might help gain a better insight into the dynamics of HIV-1 PR, especially the local dynamics of the flap tips, which may provide important guidelines for design of novel potent inhibitors.

An isoform-specific PDZ-binding motif targets type I PIP5 kinase beta to the uropod and controls polarization of neutrophil-like HL60 cells

Type I phosphatidylinositol 4-phosphate 5-kinase (PIP5KI)-beta participates in establishing polarity during leukocyte chemotaxis. Its final 83 amino acids localize PIP5KIbeta to the uropod of chemotaxing neutrophils and T cells, and interact with ezrin-radixin-moesin (ERM) proteins and EBP50 (4.1-ERM-binding phosphoprotein 50), a scaffold protein with 2 PDZ (PSD-95, disc large, ZO-1) domains. The structural motifs at the PIP5KIbeta C terminus that confer signaling specificity are, nonetheless, unknown. We show that the last 4 residues of PIP5KIbeta constitute an atypical PDZ-binding motif, which steers PIP5KIbeta to the uropod by binding to both EBP50 PDZ domains. Molecular modeling and mutagenesis indicated that PDZ-binding motif is necessary for PIP5KIbeta localization and for chemoattractant-induced neutrophil polarization. Polarity in cells that express PIP5KIbeta mutants lacking the PDZ-binding motif was restored by overexpression of PIP5KIbeta, but not of PIP5KIgamma_i2, another isoform that localizes to the neutrophil uropod. Our results identify an isoform-specific PDZ-binding motif in PIP5KIbeta, which confers specificity for PIP5KIbeta signaling at the uropod during leukocyte chemotaxis.

A crystallographic and modelling study of a human telomeric RNA (TERRA) quadruplex

DNA telomeric repeats in mammalian cells are transcribed to guanine-rich RNA sequences, which adopt parallel-stranded G-quadruplexes with a propeller-like fold. The successful crystallization and structure analysis of a bimolecular human telomeric RNA G-quadruplex, folded into the same crystalline environment as an equivalent DNA oligonucleotide sequence, is reported here. The structural basis of the increased stability of RNA telomeric quadruplexes over DNA ones and their preference for parallel topologies is described here. Our findings suggest that the 2′-OH hydroxyl groups in the RNA quadruplex play a significant role in redefining hydration structure in the grooves and the hydrogen bonding networks. The preference for specific nucleotides to populate the C3′-endo sugar pucker domain is accommodated by alterations in the phosphate backbone, which leads to greater stability through enhanced hydrogen bonding networks. Molecular dynamics simulations on the DNA and RNA quadruplexes are consistent with these findings. The computations, based on the native crystal structure, provide an explanation for RNA G-quadruplex ligand binding selectivity for a group of naphthalene diimide ligands as compared to the DNA G-quadruplex.

Evaluating molecular mechanical potentials for helical peptides and proteins

Multiple variants of the AMBER all-atom force field were quantitatively evaluated with respect to their ability to accurately characterize helix-coil equilibria in explicit solvent simulations. Using a global distributed computing network, absolute conformational convergence was achieved for large ensembles of the capped A₂₁ and F_s helical peptides. Further assessment of these AMBER variants was conducted via simulations of a flexible 164-residue five-helix-bundle protein, apolipophorin-III, on the 100 ns timescale. Of the contemporary potentials that had not been assessed previously, the AMBER-99SB force field showed significant helix-destabilizing tendencies, with beta bridge formation occurring in helical peptides, and unfolding of apolipophorin-III occurring on the tens of nanoseconds timescale. The AMBER-03 force field, while showing adequate helical propensities for both peptides and stabilizing apolipophorin-III, (i) predicts an unexpected decrease in helicity with ALA→ARG⁺ substitution, (ii) lacks experimentally observed 3₁₀ helical content, and (iii) deviates strongly from average apolipophorin-III NMR structural properties. As is observed for AMBER-99SB, AMBER-03 significantly overweighs the contribution of extended and polyproline backbone configurations to the conformational equilibrium. In contrast, the AMBER-99φ force field, which was previously shown to best reproduce experimental measurements of the helix-coil transition in model helical peptides, adequately stabilizes apolipophorin-III and yields both an average gyration radius and polar solvent exposed surface area that are in excellent agreement with the NMR ensemble.

Water film thickness-dependent conformation and diffusion of single-strand DNA on poly(ethylene glycol)-silane surface

In this paper, we investigate, using molecular dynamics simulations, the conformation and diffusion of longer and shorter single-strand DNA (ssDNA) as a function of water film thickness. While the conformation of the shorter ssDNA is significantly affected and the diffusion is suppressed with reduction in water film thickness, the conformation and diffusion of the longer DNA is not influenced. We explain our observations by considering the competition between stacking interaction of bases and solvation tendency of ssDNA. This paper suggests an approach to control the surface motion of ssDNA in nanoscale water films using film thickness.

Equilibrium conformational dynamics in an RNA tetraloop from massively parallel molecular dynamics

Conformational equilibrium within the ubiquitous GNRA tetraloop motif was simulated at the ensemble level, including 10 000 independent all-atom molecular dynamics trajectories totaling over 110 µs of simulation time. This robust sampling reveals a highly dynamic structure comprised of 15 conformational microstates. We assemble a Markov model that includes transitions ranging from the nanosecond to microsecond timescales and is dominated by six key loop conformations that contribute to fluctuations around the native state. Mining of the Protein Data Bank provides an abundance of structures in which GNRA tetraloops participate in tertiary contact formation. Most predominantly observed in the experimental data are interactions of the native loop structure within the minor groove of adjacent helical regions. Additionally, a second trend is observed in which the tetraloop assumes non-native conformations while participating in multiple tertiary contacts, in some cases involving multiple possible loop conformations. This tetraloop flexibility can act to counterbalance the energetic penalty associated with assuming non-native loop structures in forming tertiary contacts. The GNRA motif has thus evolved not only to readily participate in simple tertiary interactions involving native loop structure, but also to easily adapt tetraloop secondary conformation in order to participate in larger, more complex tertiary interactions.

How the projection domains of NF-L and alpha-internexin determine the conformations of NF-M and NF-H in neurofilaments

Making use of a numerical self-consistent field method and polymer brush concepts, we model the solvated corona of neurofilaments (NF) composed of projection domains (unstructured tails) of constituent proteins. Projections are modeled with amino acid resolution. We focus on the importance of the two shortest ones (alpha-internexin and NF-L) in regulating the conformations of the two longer ones (NF-M and NF-H) in an isolated NF. We take the wild-type NF with no alpha-internexin as the reference, for which the phosphorylation-induced translocation of M- and H-tails has been examined previously. We demonstrate that a subbrush of L-tails creates an electrostatic potential profile with an approximately parabolic shape. An experimentally relevant (2:1) ratio of L- to alpha-projections reduces the charge density of the L subbrush and shifts the translocation transition of the H-tails to slightly higher degrees of phosphorylation. Replacing all L-tails by alpha-projections destroys the substructure of the NF corona and this alters the NF response to the phosphorylation of long tails.

Structure of the preamyloid dimer of beta-2-microglobulin from covalent labeling and mass spectrometry

Beta-2-microglobulin (beta2m) self-associates into fibrillar amyloid deposits in the musculoskeletal system of patients undergoing hemodialysis treatment. Previous studies have shown that stoichiometric amounts of Cu(II) at near physiological conditions can cause beta2m to organize into native-like dimers prior to forming amyloid fibrils. Here, we report the results from selective covalent labeling reactions combined with mass spectrometry that provide insight into the amino acid residues that mediate dimer formation in the wild-type protein. Using three complementary covalent labeling reagents, we find that the dimer interface is formed by the antiparallel stacking of ABED beta-sheets from two beta2m monomers. In addition, our data clearly indicate that a dimer interface involving the interactions of D-D strands from separate protein units as seen in the recent crystal structures of two mutant beta2m oligomers is unlikely.

Evidence for a partially structured state of the amylin monomer

Islet amyloid polypeptide (amylin) is the main component in amyloid deposits formed in type II diabetes. We used triplet quenching to probe the dynamics of contact formation between the N-terminal disulfide loop and a C-terminal tryptophan in monomeric amylins from human and rat. Quenching rates measured in the absence of denaturant are four times larger than those in 6 M guanidinium chloride, indicating a decrease in the average end-to-end distance (collapse) at low denaturant concentrations. We were surprised to find an even greater (sevenfold) increase in quenching rates on removal of denaturant for a hydrophilic control peptide containing the disulfide loop compared to the same peptide without the loop (twofold change). These results suggest that collapse is driven by backbone-backbone and backbone-side chain interactions involving the disulfide loop portion of the chain rather than by the formation of side-chain hydrophobic contacts. Molecular dynamics simulations of the control peptide show that the collapse results from hydrogen-bonding interactions between the central residues of the chain and the disulfide loop. The quenching experiments also indicate that the monomer of the human, amyloidogenic form of amylin is more compact than the rat form, which does not form amyloid. We discuss these newly observed differences between human and rat amylin in solution and their possible relation to aggregation and to the physiological function of amylin binding to the calcitonin receptor.

Solution structures of rat amylin peptide: simulation, theory, and experiment

Amyloid deposits of amylin in the pancreas are an important characteristic feature found in patients with Type-2 diabetes. The aggregate has been considered important in the disease pathology and has been studied extensively. However, the secondary structures of the individual peptide have not been clearly identified. In this work, we present detailed solution structures of rat amylin using a combination of Monte Carlo and molecular dynamics simulations. A new Monte Carlo method is presented to determine the free energy of distinct biomolecular conformations. Both folded and random-coil conformations of rat amylin are observed in water and their relative stability is examined in detail. The former contains an alpha-helical segment comprised of residues 7-17. We find that at room temperature the folded structure is more stable, whereas at higher temperatures the random-coil structure predominates. From the configurations and weights we calculate the alpha-carbon NMR chemical shifts, with results that are in reasonable agreement with experiments of others. We also calculate the infrared spectrum in the amide I stretch regime, and the results are in fair agreement with the experimental line shape presented herein.

Free and ATP-bound structures of Ap4A hydrolase from Aquifex aeolicus V5

Asymmetric diadenosine tetraphosphate (Ap(4)A) hydrolases degrade the metabolite Ap(4)A back into ATP and AMP. The three-dimensional crystal structure of Ap(4)A hydrolase (16 kDa) from Aquifex aeolicus has been determined in free and ATP-bound forms at 1.8 and 1.95 A resolution, respectively. The overall three-dimensional crystal structure of the enzyme shows an alphabetaalpha-sandwich architecture with a characteristic loop adjacent to the catalytic site of the protein molecule. The ATP molecule is bound in the primary active site and the adenine moiety of the nucleotide binds in a ring-stacking arrangement equivalent to that observed in the X-ray structure of Ap(4)A hydrolase from Caenorhabditis elegans. Binding of ATP in the active site induces local conformational changes which may have important implications in the mechanism of substrate recognition in this class of enzymes. Furthermore, two invariant water molecules have been identified and their possible structural and/or functional roles are discussed. In addition, modelling of the substrate molecule at the primary active site of the enzyme suggests a possible path for entry and/or exit of the substrate and/or product molecule.

A long DNA segment in a linear nanoscale Paul trap

We study the dynamics of a linearly distributed line charge such as single stranded DNA (ssDNA) in a nanoscale, linear 2D Paul trap in vacuum. Using molecular dynamics simulations we show that a line charge can be trapped effectively in the trap for a well defined range of stability parameters. We investigated (i) a flexible bonded string of charged beads and (ii) a ssDNA polymer of variable length, for various trap parameters. A line charge undergoes oscillations or rotations as it moves, depending on its initial angle, the position of the center of mass and the velocity. The stability region for a strongly bonded line of charged beads is similar to that of a single ion with the same charge to mass ratio. Single stranded DNA as long as 40 nm does not fold or curl in the Paul trap, but could undergo rotations about the center of mass. However, we show that a stretching field in the axial direction can effectively prevent the rotations and increase the confinement stability.

Bias-dependent amino-acid-induced conductance changes in short semi-metallic carbon nanotubes

We study the interaction between short semi-metallic carbon nanotubes and different amino acids using molecular dynamics and ab initio (density functional theory/non-equilibrium Green's function) simulations. We identify two different mechanisms of nanotube conductance change upon adsorption of amino acids: one due to the change of the coordinates of the nanotube arising from van der Waals forces of interaction with the adsorbed amino acid; and one due to electrostatic interactions, which appear only in the case of charged amino acids. We also find that the transport mechanism and the changes in the conductance of the tube upon amino acid adsorption are bias dependent.

Conformational transitions upon ligand binding: holo-structure prediction from apo conformations

Biological function of proteins is frequently associated with the formation of complexes with small-molecule ligands. Experimental structure determination of such complexes at atomic resolution, however, can be time-consuming and costly. Computational methods for structure prediction of protein/ligand complexes, particularly docking, are as yet restricted by their limited consideration of receptor flexibility, rendering them not applicable for predicting protein/ligand complexes if large conformational changes of the receptor upon ligand binding are involved. Accurate receptor models in the ligand-bound state (holo structures), however, are a prerequisite for successful structure-based drug design. Hence, if only an unbound (apo) structure is available distinct from the ligand-bound conformation, structure-based drug design is severely limited. We present a method to predict the structure of protein/ligand complexes based solely on the apo structure, the ligand and the radius of gyration of the holo structure. The method is applied to ten cases in which proteins undergo structural rearrangements of up to 7.1 A backbone RMSD upon ligand binding. In all cases, receptor models within 1.6 A backbone RMSD to the target were predicted and close-to-native ligand binding poses were obtained for 8 of 10 cases in the top-ranked complex models. A protocol is presented that is expected to enable structure modeling of protein/ligand complexes and structure-based drug design for cases where crystal structures of ligand-bound conformations are not available.

Electronic continuum model for molecular dynamics simulations of biological molecules

Electronic polarizability is an important factor in molecular interactions. In the conventional force fields such as AMBER or CHARMM, however, there is inconsistency in how the effect of electronic dielectric screening of Coulombic interactions, inherent for the condensed phase media, is treated. Namely, the screening appears to be accounted for via effective charges only for neutral moieties, whereas the charged residues are treated as if they were in vacuum. As a result, the electrostatic interactions between ionized groups are exaggerated in molecular simulations by the factor of about 2. The discussed here MDEC (Molecular Dynamics in Electronic Continuum) model provides a theoretical framework for modification of the standard non-polarizable force fields to make them consistent with the idea of uniform electronic screening of partial atomic charges. The present theory states that the charges of ionized groups and ions should be scaled; i.e. reduced by a factor of about 0.7. In several examples, including the interaction between Na (+) ions, which is of interest for ion-channel simulations, and the dynamics of an important salt-bride in Cytochrome c Oxidase, we compared the standard non-polarizable MD simulations with MDEC simulations, and demonstrated that MDEC charge scaling procedure results in more accurate interactions. The inclusion of electronic screening for charged moieties is shown to result in significant changes in protein dynamics and can give rise to new qualitative results compared with the traditional non-polarizable force fields simulations.

Simulations of RNA interactions with monovalent ions

RNA folding and binding reactions are mediated by interactions with ions that make up the surrounding aqueous electrolytic milieu. Although Mg(2+) ions are often implicated as being crucial for RNA folding, it is known that folding is feasible in high concentrations of monovalent alkali-halide salts. Experiments have yielded important information regarding the salt dependence of RNA stability. Recent work has shown that molecular simulations based on explicit representations of solvent molecules and monovalent ions can also provide useful insights regarding the ionic atmospheres around model RNA systems. These insights can help rationalize intriguing observations regarding the dependence of RNA stability on cation type providing one pays attention to important considerations that go into the proper design of molecular simulations. These issues are discussed in detail and the methods are applied to an A-form RNA and B-form DNA sequence. The results of these simulations are compared to previous work on a kissing-loop system with analogous sequence. In particular, ionic atmospheres obtained from molecular simulations are compared to those obtained using the nonlinear Poisson Boltzmann model for continuum electrostatics for these three different nucleic acid systems. The comparisons indicate reasonable agreement in terms of coarse-grained observables such as the numbers of counterions accumulated around the solutes. However, details of the ionic atmospheres, captured in terms of spatial free energy density profiles, are quite different between the two approaches. These comparisons suggest the need for improvements in continuum models to capture sequence-specific effects, ion-ion correlation, and the effects of partial dehydration of ions.

Salt dependence of an alpha-helical peptide folding energy landscapes

We used CD, UV resonance Raman spectroscopy, and molecular dynamics simulation to examine the impact of salts on the conformational equilibria and the Ramachandran Psi angle (un)folding Gibbs free energy landscape coordinate of a mainly polyalanine alpha-helical peptide, AP of sequence AAAAA(AAARA)(3)A. NaClO(4) stabilizes alpha-helical-like conformations more than does NaCl, which stabilizes more than Na(2)SO(4) at identical ionic strengths. This alpha-helix stabilization ordering is the reverse of the Hofmeister series of anions in their ability to disorder water hydrogen bonding. Much of the NaClO(4) alpha-helix stabilization results from ClO(4)(-) association with the AP terminal -NH(3)(+) groups and Arg side chains. ClO(4)(-) stabilizes 3(10)-helix conformations but destabilizes turn conformations. The decreased Cl(-) and SO(4)(2-) AP alpha-helix stabilization probably results from a decreased association with the Arg and terminal -NH(3)(+) groups. Cl(-) is expected to have a smaller binding affinity and thus stabilizes alpha-helical conformations intermediately between NaClO(4) and Na(2)SO(4). Electrostatic screening stabilizes pi-bulge conformations.

Secondary structure propensities in peptide folding simulations: a systematic comparison of molecular mechanics interaction schemes

We present a systematic study directed toward the secondary structure propensity and sampling behavior in peptide folding simulations with eight different molecular dynamics force-field variants in explicit solvent. We report on the combinational result of force field, water model, and electrostatic interaction schemes and compare to available experimental characterization of five studied model peptides in terms of reproduced structure and dynamics. The total simulation time exceeded 18 mus and included simulations that started from both folded and extended conformations. Despite remaining sampling issues, a number of distinct trends in the folding behavior of the peptides emerged. Pronounced differences in the propensity of finding prominent secondary structure motifs in the different applied force fields suggest that problems point in particular to the balance of the relative stabilities of helical and extended conformations.

Non-sequence-specific interactions can account for the compaction of proteins unfolded under "native" conditions

Proteins unfolded by high concentrations of chemical denaturants adopt expanded, largely structure-free ensembles of conformations that are well approximated as random coils. In contrast, globular proteins unfolded under less denaturing conditions (via mutations, or transiently unfolded after a rapid jump to native conditions) and molten globules (arising due to mutations or cosolvents) are often compact. Here we explore the origins of this compaction using a truncated equilibrium-unfolded variant of the 57-residue FynSH3 domain. As monitored by far-UV circular dichroism, NMR spectroscopy, and hydrogen-exchange kinetics, CDelta4 (a 4-residue carboxy-terminal deletion variant of FynSH3) appears to be largely unfolded even in the absence of denaturant. Nevertheless, CDelta4 is quite compact under these conditions, with a hydrodynamic radius only slightly larger than that of the native protein. In order to understand the origins of this molten-globule-like compaction, we have characterized a random sequence polypeptide of identical amino acid composition to CDelta4. Notably, we find that the hydrodynamic radius of this random sequence polypeptide also approaches that of the native protein. Thus, while native-like interactions may contribute to the formation of compact "unfolded" states, it appears that non-sequence-specific monomer-monomer interactions can also account for the dramatic compaction observed for molten globules and the "physiological" unfolded state.

Quantifying Correlations Between Allosteric Sites in Thermodynamic Ensembles

Allostery describes altered protein function at one site due to a perturbation at another site. One mechanism of allostery involves correlated motions, which can occur even in the absence of substantial conformational change. We present a novel method, "MutInf", to identify statistically significant correlated motions from equilibrium molecular dynamics simulations. Our approach analyzes both backbone and sidechain motions using internal coordinates to account for the gear-like twists that can take place even in the absence of the large conformational changes typical of traditional allosteric proteins. We quantify correlated motions using a mutual information metric, which we extend to incorporate data from multiple short simulations and to filter out correlations that are not statistically significant. Applying our approach to uncover mechanisms of cooperative small molecule binding in human interleukin-2, we identify clusters of correlated residues from 50 ns of molecular dynamics simulations. Interestingly, two of the clusters with the strongest correlations highlight known cooperative small-molecule binding sites and show substantial correlations between these sites. These cooperative binding sites on interleukin-2 are correlated not only through the hydrophobic core of the protein but also through a dynamic polar network of hydrogen bonding and electrostatic interactions. Since this approach identifies correlated conformations in an unbiased, statistically robust manner, it should be a useful tool for finding novel or "orphan" allosteric sites in proteins of biological and therapeutic importance.

Energy landscape of the trpzip2 peptide

Replica exchange simulations are used to study the energy landscape of trpzip2, a model beta-hairpin system, using the AMBER99sb force field and explicit solvent. The total simulation time is 300 ns per replica (approximately 10 mus total). The trp side chains are observed to adopt multiple packing arrangements with a freezing temperature below 273 K in the simulated system. The secondary structure and native hydrogen bonds melt out cooperatively around 273 K. The residual beta-strand structure and antiparallel bonding persist at high temperature. These results provide a model for the three apparent melting transitions observed experimentally in this system. The dominant folding mechanism of trpzip2 in this model appears to be zipping, which is consistent with recent measurements on similar hairpins. Structures for which the turn is native-like and the termini are non-native-like collectively form a metastable intermediate. Most of the stabilizing enthalpy is gained after the formation of the turn. Equilibrium thermodynamic quantities are compared against experiment. Although the AMBER99sb force field reproduces the native structure with good fidelity, the stability of the native state is significantly underpredicted with a melting temperature near 273 K, and the relative heat capacity is only about one tenth of its experimental value.

Molecular determinants of snurportin 1 ligand affinity and structural response upon binding

The transport of large biomolecules such as proteins and RNA across nuclear pore complexes is a field of strong interest and research. Although the basic mechanisms are fairly well understood, the details of the underlying intermolecular interaction within these transport complexes are still unclear. The recognition dynamics and energetics of cargo binding to the transport receptor are not yet resolved. Here, the binding of dimethylated RNA-caps to snurportin 1 is studied by molecular-dynamics simulations. The simulations reveal a strong structural response of the protein upon RNA-cap release. In particular, major rearrangements occur in regions already intrinsically flexible in the holo structure. Additionally, the difference in free energy of binding to snurportin 1 between the two methylation states of the RNA-cap, responsible for the directionality of the transport is quantified. In particular, desolvation of the ligand is revealed as the key-step in binding to snurportin 1. These findings suggest that the binding of m(3)G-capped RNA is mainly driven by the enhanced water entropy gain of the solvation shell.

Kinetics of helix unfolding: molecular dynamics simulations with milestoning

The unfolding process of a helical heteropeptide is studied by computer simulation in explicit solvent. A combination of a functional optimization to determine the reaction coordinate and short time trajectories between "milestones" is used to study the kinetics of the unfolding. One hundred unfolding trajectories along three different unfolding pathways are computed between all nearby milestones, providing adequate statistics to compute the overall first passage time. The radius of gyration is smaller for intermediate configurations compared to the initial and final states, suggesting that the kinetics (but not the thermodynamics) is sensitive to pressure. The transitions are dominated by local bond rotations (the psi dihedral angle) that are assisted by significant nonmonotonic fluctuations of nearby torsions. The most effective unfolding pathway is via the N-terminal.