Structure and atomic fluctuation patterns of potato carboxypeptidase a inhibitor protein

Molecular dynamics (MD) simulation methods were applied to the study of the structural and dynamic fluctuation properties of the potato carboxypeptidase A inhibitor protein (PCI) immersed in a bath of 1259 water molecules. A trajectory of 200 ps was generated at constant temperature and pressure. The crystallographic structure of PCI, as found in its complex with bovine carboxy-peptidase A (CPA), was used to seed the MD simulation. Analyses show that the structure of the PCI core is fairly rigid and stable in itself, and that little deformation is caused by the protein-protein interactions found in the PCI-CPA complex. The N-terminal tail fluctuates to approach the core structure and appears as a relatively disordered region. In contrast, the conformations of the C-terminal tail, which is involved in the inhibitory mechanism, fluctuates in the neighborhood of the X-ray structure in orientations which facilitate CPA binding. Comparison with the structural entries for PCI in water obtained from both 2D-NMR experiments and X-ray data shows that important features of the MD structural results fluctuates between the initial crystal values and those obtained from the NMR solution structure. This fluctuation is not uniform; minor regions move away from the X-ray conformation while they do not approach the NMR conformation. The results reported suggest that the trajectory is long enough to show a behavior that is consistent with the conformational space available to the protein in solution.Abbreviations CPA Carboxypeptidase - DG Distance Geometry - NMR Nuclear Magnetic Resonance - NIS Non Inertial Solvent - MD Molecular Dynamics - PBC Periodic Boundary Conditions - PCI Potato Carboxypeptidase Inhibitor - RMSD Root Mean Square Deviation - a.m.u. Atomic mass units Correspondence to: O. Tapia     

A systematic comparison of three structure determination methods from NMR data: Dependence upon quality and quantity of data

Summary We have systematically examined how the quality of NMR protein structures depends on (1) the number of NOE distance constraints. (2) their assumed precision, (3) the method of structure calculation and (4) the size of the protein. The test sets of distance constraints have been derived from the crystal structures of crambin (5 kDa) and staphylococcal nuclease (17 kDa). Three methods of structure calculation have been compared: Distance Geometry (DGEOM), Restrained Molecular Dynamics (XPLOR) and the Double Iterated Kalman Filter (DIKF). All three methods can reproduce the general features of the starting structure under all conditions tested. In many instances the apparent precision of the calculated structure (as measured by the RMS dispersion from the average) is greater than its accuracy (as measured by the RMS deviation of the average structure from the starting crystal structure). The global RMS deviations from the reference structures decrease exponentially as the number of constraints is increased, and after using about 30% of all potential constraints, the crrors asymptotically approach a limiting value. Increasing the assumed precision of the constraints has the same qualitative effect as increasing the number of constraints. For comparable numbers of constraints/residue, the precision of the calculated structure is less for the larger than for the smaller protein, regardless of the method of calculation. The accuracy of the average structure calculated by Restrained Molecular Dynamics is greater than that of structures obtained by purely geometric methods (DGEOM and DIKF).     

Microcomputers Against Supercomputers? On Geometric Partition of the Computational Box for Vectorized MD Algorithms

Abstract

In the paper three vector algorithms (with original data structures and no nearest neighbour lists) for Molecular Dynamics simulation are compared and their timings presented. The timings have been obtained on ETA 10-P*108 and IBM 3090/150E (with Vector Facility) computers as well as on several microprocessors (e.g. MOTOROLA MC68020/68881, INTEL i80386/387, i486 and INMOS TRANSPUTER T800).     

Grand Canonical Ensemble Monte Carlo Simulation on a Transputer Array

Five methods are described for the distribution over a 3D Transputer array of the calculation of the pair interaction component of particle energy. The most efficient method, expressed in terms of the time to complete a simulation, depends on the size of the simulation and the Transputer array. This dependence is quantified, with emphasis on Grand Canonical Ensemble Monte Carlo simulation, and yields criteria for the optimum strategy for parallel implementation of GCEMC algorithms. The equations derived are generally applicable, and have implications for the programming of Molecular Dynamics simulations.     

Computational simulations on the fish-type-Ⅱ antifreeze protein-ice-solvent system

Based on the computational simulation with the vacuum environment for the fish-type-Ⅱ antifreeze proteinice-solvent (water)system,the multi-complex system of the antifreeze protein-ice-water has been constructed and calculated.We have studied the interaction of such proteinice system with water solvent through the dynamics simulation with 350 ps.By employing the Molecular Dynamics simulation and semi-empirical method calculation,we have further investigated the interface properties of the antifreeze protein and ice crystal combined system.Consequently,a water solvent affects significantly the properties of this combined system.     

Predicting Atomic Details of the Unfolding Pathway for YibK,a Knotted Protein from the SPOUT Superfamily

Abstract

Several protein structures have been reported to contain intricate knots of the polypeptide backbone but the mechanism of the (un)folding process of knotted proteins remains unknown. The members of the SPOUT superfamily of RNA methyltransferases are some of the most intensely studied systems for investigation of the knot formation and function. YibK (whose biochemical function remains unknown) is the representative protein of the SPOUT superfamily. This protein exhibits a deep trefoil knot at the C-terminus.

We conducted an extensive computational analysis of the unfolding process for the monomeric form of YibK. In order to predict the (un)folding pathway of YibK, we have calculated the order of secondary structure disassembly using UNFOLD, and performed thermal unfolding simulations using classical Molecular Dynamics (MD), as well as simulations employing reduced representation of the peptide chain using either MD with the UNRES method or the Monte Carlo (MC) unfolding with the REFINER method.

Results obtained from all methods used in this work are in qualitative agreement. We found that YibK unfolds through four intermediate states. The trefoil knot in YibK disappears at the end of the unfolding process, long after the protein loses its native topology. We observed that the C-terminus leaves the knotting loop folded into a hairpin-like structure, in agreement with the results of coarse-grained simulation reported earlier. We propose that the folding pathway of YibK corresponds to the reversed sequence of events observed in the unfolding pathway elucidated in this study. Thus, we predict that the knot formation is the slowest part of the YibK folding process.     

Towards an Exhaustive Sampling of the Configurational Spaces of the Two Forms of the Peptide Hormone Guanylin

Abstract

The recently introduced Essential Dynamics sampling method is extended such that an exhaustive sampling of the available (backbone) configurational space can be achieved. From an initial Molecular Dynamics simulation an approximated definition of the essential subspace is obtained. This subspace is used to direct subsequent simulations by means of constraint forces. The method is applied to the peptide hormone guanylin, solvated in water, of which the structure was determined recently. The peptide exists in two forms and for both forms, an extensive sampling was produced. The sampling algorithm fills the available space (of the essential coordinates used in the procedure) at a rate that is approximately six to seven times larger than that for traditional Molecular Dynamics. The procedure does not cause any significant perturbation, which is indicated by the fact that free Molecular Dynamics simulations started at several places in the space defined by the Essential Dynamics sample that complete space. Moreover, analyses of the average free Molecular Dynamics step have shown that nowhere except close to the edge of the available space, there are regions where the system shows a drift in a particular direction. This result also shows that in principle, the essential subspace is a constant free energy surface, with well-defined and steep borders, in which the system moves diffusively. In addition, a comparison between two independent essential dynamics sampling runs, of one form of the peptide, shows that the obtained essential subspaces are virtually identical.     

Gibbs-Ensemble Molecular Dynamics: Liquid-Gas Equilibria for Lennard-Jones Spheres and n-Hexane

Abstract

We present a novel method to simulate phase equilibria in atomic and molecular systems. The method is a Molecular Dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles (or molecules) in a thermodynamically consistent fashion. Here we pres the derivation of the generalized equations of motion and discuss the relation of the resulting trajectory averages to the relevant ensemble. We test this Gibbs-Ensemble Molecular Dynamics algorithm by applying it to an atomic and a molecular system, i.e. to the liquid-gas coexistence in a Lennard-Jones fluid and in n-hexane. In both cases our results are in good accord with previous mean field and Gibbs-Ensemble Monte Carlo results as well as with the experimental data in the case of hexane. We also show that our Gibbs-Ensemble Molecular Dynamics algorithm like other Molecular Dynamics techniques can be used to study the dynamics of the system. Self-diffusion coefficients calculated with this method are in agreement with the result of conventional constant temperature Molecular Dynamics.     

Synthesis and conformational analysis of apamin analogues with natural and non-natural cystine/selenocystine connectivities

By replacing two cysteine residues in apamin with selenocysteine, the three possible isomers related to the side-chain connectivities of a bis-cystinyl-peptide were synthesized in regioselective manner exploiting the low redox potential of the diselenide bond. Nuclear magnetic resonance conformational analysis of monoselenocystine analogue apamin with the natural diselenide/disulfide network confirmed the highly isomorphous character of the sulfur replacement with selenium despite its slightly larger atomic radius and increased bond lengths. The comparative conformational analysis of the apamin analogues containing the non-natural side-chain links with wild type apamin clearly revealed retention of the main structural fold and thus the high propensity of these small molecules to adopt the secondary structure elements present in natural apamin. These findings offered interesting hints for a better understanding of the oxidative refolding pathway of the bis-cystinyl peptide that leads exclusively to the correct natural isomer.     

All-atom Monte Carlo simulation of GCAA RNA folding

We report a detailed all-atom simulation of the folding of the GCAA RNA tetraloop. The GCAA tetraloop motif is a very common and thermodynamically stable secondary structure in natural RNAs. We use our simulation methods to study the folding behavior of a 12-base GCAA tetraloop structure with a four-base helix adjacent to the tetraloop proper. We implement an all-atom Monte Carlo (MC) simulation of RNA structural dynamics using a Go potential. Molecular dynamics (MD) simulation of RNA and protein has realistic energetics and sterics, but is extremely expensive in terms of computational time. By coarsely treating non-covalent energetics, but retaining all-atom sterics and entropic effects, all-atom MC techniques are a useful method for the study of protein and now RNA. We observe a sharp folding transition for this structure, and in simulations at room temperature the state histogram shows three distinct minima: an unfolded state (U), a more narrow intermediated state (I), and a narrow folded state (F). The intermediate consists primarily of structures with the GCAA loop and some helix hydrogen bonds formed. Repeated kinetic folding simulations reveal that the number of helix base-pairs forms a simple 1D reaction coordinate for the I-->N transition.     

Practical,Microfabrication-Free Device for Single-Cell Isolation

Microfabricated devices have great potential in cell-level studies, but are not easily accessible for the broad biology community. This paper introduces the Microscale Oil-Covered Cell Array (MOCCA) as a low-cost device for high throughput single-cell analysis that can be easily produced by researchers without microengineering knowledge. Instead of using microfabricated structures to capture cells, MOCCA isolates cells in discrete aqueous droplets that are separated by oil on patterned hydrophilic areas across a relatively more hydrophobic substrate. The number of randomly seeded Escherichia coli bacteria in each discrete droplet approaches single-cell levels. The cell distribution on MOCCA is well-fit with Poisson distribution. In this pioneer study, we created an array of 900-picoliter droplets. The total time needed to seed cells in ∼3000 droplets was less than 10 minutes. Compared to traditional microfabrication techniques, MOCCA dramatically lowers the cost of microscale cell arrays, yet enhances the fabrication and operational efficiency for single-cell analysis.     

Structural heterogeneity and quantitative FRET efficiency distributions of polyprolines through a hybrid atomistic simulation and Monte Carlo approach

Förster Resonance Energy Transfer (FRET) experiments probe molecular distances via distance dependent energy transfer from an excited donor dye to an acceptor dye. Single molecule experiments not only probe average distances, but also distance distributions or even fluctuations, and thus provide a powerful tool to study biomolecular structure and dynamics. However, the measured energy transfer efficiency depends not only on the distance between the dyes, but also on their mutual orientation, which is typically inaccessible to experiments. Thus, assumptions on the orientation distributions and averages are usually made, limiting the accuracy of the distance distributions extracted from FRET experiments. Here, we demonstrate that by combining single molecule FRET experiments with the mutual dye orientation statistics obtained from Molecular Dynamics (MD) simulations, improved estimates of distances and distributions are obtained. From the simulated time-dependent mutual orientations, FRET efficiencies are calculated and the full statistics of individual photon absorption, energy transfer, and photon emission events is obtained from subsequent Monte Carlo (MC) simulations of the FRET kinetics. All recorded emission events are collected to bursts from which efficiency distributions are calculated in close resemblance to the actual FRET experiment, taking shot noise fully into account. Using polyproline chains with attached Alexa 488 and Alexa 594 dyes as a test system, we demonstrate the feasibility of this approach by direct comparison to experimental data. We identified cis-isomers and different static local environments as sources of the experimentally observed heterogeneity. Reconstructions of distance distributions from experimental data at different levels of theory demonstrate how the respective underlying assumptions and approximations affect the obtained accuracy. Our results show that dye fluctuations obtained from MD simulations, combined with MC single photon kinetics, provide a versatile tool to improve the accuracy of distance distributions that can be extracted from measured single molecule FRET efficiencies.     

Conjugated double bonds in lipid bilayers: a molecular dynamics simulation study

Conjugated linoleic acids (CLA) are found naturally in dairy products. Two isomers of CLA, that differ only in the location of cis and trans double bonds, are found to have distinct and different biological effects. The cis 9 trans 11 (C9T11) isomer is attributed to have the anti-carcinogenic effects, while the trans 10 cis 12 (T10C12) isomer is believed to be responsible for the anti-obesity effects. Since dietary CLA are incorporated into membrane phospholipids, we have used Molecular Dynamics (MD) simulations to investigate the comparative effects of the two isomers on lipid bilayer structure. Specifically, simulations of phosphatidylcholine lipid bilayers in which the sn-2 chains contained one of the two isomers of CLA were performed. Force field parameters for the torsional potential of double bonds were obtained from ab initio calculations. From the MD trajectories we calculated and compared structural properties of the two lipid bilayers, including areas per molecule, density profiles, thickness of bilayers, tilt angle of tail chains, order parameters profiles, radial distribution function (RDF) and lateral pressure profiles. The main differences found between bilayers of the two CLA isomers, are (1) the order parameter profile for C9T11 has a dip in the middle of sn-2 chain while the profile for T10C12 has a deeper dip close to terminal of sn-2 chain, and (2) the lateral pressure profiles show differences between the two isomers. Our simulation results reveal localized physical structural differences between bilayers of the two CLA isomers that may contribute to different biological effects through differential interactions with membrane proteins or cholesterol.     

Microsecond scale replica exchange molecular dynamic simulation of villin headpiece: an insight into the folding landscape

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 μs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 ? was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins.     

Gibbs-Ensemble Molecular Dynamics: A New Method for Simulations Involving Particle Exchange

We discuss a novel simulation method suitable for simulating phenomena involving particle exchange. The method is a molecular dynamics version of the Gibbs-Ensemble Monte Carlo technique, which has been developed some years ago for the direct simulation of phase equilibria in fluid systems. The idea is to have two separate simulation boxes, which can exchange particles or molecules in a thermodynamically consistent fashion. We discuss the general idea of the Gibbs-Ensemble Molecular Dynamics technique and present examples for different simple atomic and molecular fluids. Specifically we will discuss Gibbs-Ensemble Molecular Dynamics simulations of gas-liquid and liquid-solid equilibria in Lennard-Jones systems and in hexane as well as an application of the method to adsorption.     

Microsecond Scale Replica Exchange Molecular Dynamic Simulation of Villin Headpiece: An Insight into the Folding Landscape

Abstract

Reaching the experimental time scale of millisecond is a grand challenge for protein folding simulations. The development of advanced Molecular Dynamics techniques like Replica Exchange Molecular Dynamics (REMD) makes it possible to reach these experimental timescales. In this study, an attempt has been made to reach the multi microsecond simulation time scale by carrying out folding simulations on a three helix bundle protein, Villin, by combining REMD and Amber United Atom model. Twenty replicas having different temperatures ranging from 295 K to 390 K were simulated for 1.5 μs each. The lowest Root Mean Square Deviation (RMSD) structure of 2.5 Å was obtained with respect to native structure (PDB code 1VII), with all the helices formed. The folding population landscapes were built using segment-wise RMSD and Principal Components as reaction coordinates. These analyses suggest the two-stage folding for Villin. The combination of REMD and Amber United Atom model may be useful to understand the folding mechanism of various fast folding proteins     

Polypeptide folding with off-lattice Monte Carlo dynamics: the method

The MC dynamics of an off-lattice all-atom protein backbone model with rigid amide planes are studied. The only degrees of freedom are the dihedral angle pairs of the C-atoms. Conformational changes are generated by Monte Carlo (MC) moves. The MC moves considered are single rotations (simple moves, SM's) giving rise to global conformational changes or, alternatively, cooperative rotations in a window of amide planes (window moves, WM's) generating local conformational changes in the window. Outside the window the protein conformation is kept invariant by constraints. These constraints produce a bias in the distribution of dihedral angles. The WM's are corrected for this bias by suitable Jacobians. The energy function used is derived from the CHARMM force field. In a first application to polyalanine it is demonstrated that WM's sample the conformational space more efficiently than SM's.Abbreviations CPU Central Processing Unit - MC Monte Carlo - MCD Monte Carlo Dynamics - MD Molecular Dynamics - RMS Root-Mean-Square - RMSD Root-Mean-Square-Deviation - SM Simple Move - WM Window Move     

Modeling of angiogenin - 3-NMP complex

Angiogenin belongs to the Ribonuclease superfamily and has a weak enzymatic activity that is crucial for its biological function of stimulating blood vessel growth. Structural studies on ligand bound Angiogenin will go a long way in understanding the mechanism of the protein as well as help in designing drugs against it. In this study we present the first available structure of nucleotide ligand bound Angiogenin obtained by computer modeling. The importance of this study in itself notwithstanding, is a precursor to modeling a full dinucleotide substrate onto Angiogenin. Bovine Angiogenin, the structure of which has been solved at a high resolution, was earlier subjected to Molecular Dynamics simulations for a nanosecond. The MD structures offer better starting points for docking as they offer lesser obstruction than the crystal structure to ligand binding. The MD structure with the least serious short contacts was modeled to obtain a steric free Angiogenin - 3' mononucleotide complex structure. The structures were energetically minimized and subjected to a brief spell of Molecular Dynamics. The results of the simulation show that all the ligand-Angiogenin interactions and hydrogen bonds are retained, redeeming the structure and docking procedure. Further, following ligand - protein interactions in the case of the ligands 3'-CMP and 3'-UMP we were able to speculate on how Angiogenin, a predominantly prymidine specific ribonuclease prefers Cytosine to Uracil in the first base position.     

Ewald Summation in the Molecular Dynamics Simulation of Large Ionic Systems

Abstract

The accuracy and efficiency of the direct Ewald summation are discussed in terms of the size of a Molecular Dynamics (MD) ionic system and the ranges of the r-space and q-space summations. The dependence of the convergence parameter α on the size of the system and on the choice of cut-off radius for the short-range potential is given. The possibility of neglecting the q-space term for large ionic systems is discussed in terms of the accuracy and efficiency of the simulation.     

The asparagine-stabilized beta-turn of apamin: contribution to structural stability from dynamics simulation and amide hydrogen exchange analysis

Molecular dynamics simulations of bee venom apamin, and an analogue having an Asn to Ala substitution at residue 2 (apamin-N2A), were analyzed to explore the contribution of hydrogen bonds involving Asn2 to local (beta-turn residues N2, C3, K4, A5) and global stability. The wild-type peptide retained a stable conformation during 2.4 ns of simulation at 67 degrees C, with high beta-turn stability characterized by backbone-side chain hydrogen bonds involving beta-turn residues K4 and A5, with the N2 side chain amide carbonyl. The loss of stabilizing interactions involving the N2 side chain resulted in the loss of the beta-turn conformation in the apamin N2A simulations (27 or 67 degrees C). This loss of beta-turn stability propagates throughout the peptide structure, with destabilization of the C-terminal helix connected to the N-terminal region by two disulfide bonds. Backbone stability in a synthetic peptide analogue (apamin-N2A) was characterized by NMR and amide hydrogen exchange measurements. Consistent with the simulations, loss of hydrogen bonds involving the N2 side chain resulted in destabilization of both the N-terminal beta-turn and the C-terminal helix. Amide exchange protection factors in the C-terminal helix were reduced by 9-11-fold in apamin N2A as compared with apamin, corresponding to free energy (deltaDeltaG(uf)) of around 1.5 kcal M(-1) at 20 degrees C. This is equivalent to the contribution of hydrogen bond interactions involving the N2 side chain to the stability of the beta-turn. Together with additional measures of exchange protection factors, the three main contributions to backbone stability in apamin that account for virtually the full thermodynamic stability of the peptide have been quantitated.