Coarse-grained molecular simulations of large biomolecules

Recently, we saw a dramatic increase in the number of researches that rely on coarse-grained (CG) simulations for large biomolecules. Here, first, we briefly describe recently developed and used CG models for proteins and nucleic acids. Balance between structure-based and physico-chemical terms is a key issue. We also discuss the multiscale algorithms used to derive CG parameters. Next, we comment on the dynamics used in CG simulations with an emphasis on the importance of hydrodynamic interactions. We then discuss the pros and cons of CG simulations. Finally, we overview recent exciting applications of CG simulations. Publicly available tools and software for CG simulations are also summarized.     

Molecular dynamics simulation of a phosphatidylglycerol membrane

Although molecular dynamics simulations are an important tool for studying membrane systems, relatively few simulations have used anionic lipids. This paper reports the first simulation of a pure phosphatidylglycerol (PG) bilayer. The properties of this equilibrated palmitoyloleoylphosphatidylglycerol membrane agree with experimental observations of PG membranes and with previous simulations of monolayers and mixed bilayers containing PG lipids. These simulations also provide interesting insights into hydrogen bonding interactions in PG membranes. This equilibrated membrane will be a useful starting point for simulations of membrane proteins interacting with PG lipids.     

Ion flux onto a nonconducting surface in a beam-plasma discharge: Simulation and experiment

Computer simulations of a plasma sheath formed near a nonconducting surface in a beamplasma discharge (BPD) at a low magnetic field in the reactor chamber are performed. The results of computer simulations are compared with experimental data. The optimal mode of GaAs etching in a BPD-based plasma processing reactor is chosen using results of computer simulations. Simulations for carbon (DLC) film sputtering modes in the framework of a similar model are also performed.     

Detailed molecular dynamics simulations of model biological membranes containing cholesterol

Detailed molecular dynamics simulations performed to study the nature of lipid raft domains that appear in model membranes are reviewed in this paper. The described simulations were performed on hydrated bilayers containing binary mixtures of cholesterol with phospholipids and also on ternary mixtures containing cholesterol, a phospholipid with a high main transition temperature T_m, and a phospholipid with a low transition temperature T_m. These simulations provide qualitative and semi-quantitative information about cholesterol-lipid interactions and also a testing ground for major assumptions made to explain the nature of lipid rafts in model membranes.     

Characterizing rare-event property distributions via replicate molecular dynamics simulations of proteins

As computational resources increase, molecular dynamics simulations of biomolecules are becoming an increasingly informative complement to experimental studies. In particular, it has now become feasible to use multiple initial molecular configurations to generate an ensemble of replicate production-run simulations that allows for more complete characterization of rare events such as ligand-receptor unbinding. However, there are currently no explicit guidelines for selecting an ensemble of initial configurations for replicate simulations. Here, we use clustering analysis and steered molecular dynamics simulations to demonstrate that the configurational changes accessible in molecular dynamics simulations of biomolecules do not necessarily correlate with observed rare-event properties. This informs selection of a representative set of initial configurations. We also employ statistical analysis to identify the minimum number of replicate simulations required to sufficiently sample a given biomolecular property distribution. Together, these results suggest a general procedure for generating an ensemble of replicate simulations that will maximize accurate characterization of rare-event property distributions in biomolecules.     

Molecular dynamics simulations of arachidonic acid complexes with COX-1 and COX-2: insights into equilibrium behavior

The cyclooxygenase (COX) enzymes are responsible for the committed step in prostaglandin biosynthesis, the generation of prostaglandin H(2). As a result, these enzymes are pharmacologically important targets for nonsteroidal antiinflammatory drugs, such as aspirin and newer COX-2 selective inhibitors. The cyclooxygenases are functional homodimers, and each subunit contains both a cyclooxygenase and a peroxidase active site. These enzymes are quite interesting mechanistically, as the conversion of arachidonic acid to prostaglandin H(2) requires two oxygenation and two cyclization reactions, resulting in the formation of five new chiral centers with nearly absolute regio- and stereochemical fidelity. We have used molecular dynamics (MD) simulations to investigate the equilibrium behavior of both COX-1 and COX-2 enzyme isoforms with bound arachidonate. These simulations were compared with reference simulations of arachidonate in solution to explore the effect of enzyme on substrate conformation and positioning in the active site. The simulations suggest that the substrate has greater conformational freedom in the COX-2 active site, consistent with the larger COX-2 active site volume observed in X-ray crystal structures. The simulations reveal different conformational behavior for arachidonate in each subunit over the course of extended equilibrium MD simulations. The simulations also provide detailed information for several protein channels that might be important for oxygen and water transport to or from active sites or for intermediate trafficking between the cyclooxygenase and peroxidase active sites. The detailed comparisons for COX-1 versus COX-2 active site structural fluctuations may also provide useful information for design of new isozyme-selective inhibitors.     

Molecular dynamics studies of a DNA-binding protein: 1. A comparison of the trp repressor and trp aporepressor aqueous simulations.

The results of two 30-ps molecular dynamics simulations of the trp repressor and trp aporepressor proteins are presented in this paper. The simulations were obtained using the AMBER molecular mechanical force field and in both simulations a 6-A shell of TIP3P waters surrounded the proteins. The trp repressor protein is a DNA-binding regulatory protein and it utilizes a helix-turn-helix (D helix-turn-E helix) motif to interact with DNA. The trp aporepressor, lacking two molecules of the L-tryptophan corepressor, cannot bind specifically to DNA. Our simulations show that the N- and C-termini and the residues in and near the helix-turn-helix motifs are the most mobile regions of the proteins, in agreement with the X-ray crystallographic studies. Our simulations also find increased mobility of the residues in the turn-D helix-turn regions of the proteins. We find the average distance separating the DNA-binding motifs to be larger in the repressor as compared to the aporepressor. In addition to examining the protein residue fluctuations and deviations with respect to X-ray structures, we have also focused on backbone dihedral angles and corepressor hydrogen-bonding patterns in this paper.     

Computer simulations of membrane proteins

Computer simulations are rapidly becoming a standard tool to study the structure and dynamics of lipids and membrane proteins. Increasing computer capacity allows unbiased simulations of lipid and membrane-active peptides. With the increasing number of high-resolution structures of membrane proteins, which also enables homology modelling of more structures, a wide range of membrane proteins can now be simulated over time spans that capture essential biological processes. Longer time scales are accessible by special computational methods. We review recent progress in simulations of membrane proteins.     

Quaternary association in β-prism I fold plant lectins: Insights from X-ray crystallography, modelling and molecular dynamics

Dimeric banana lectin and calsepa, tetrameric artocarpin and octameric heltuba are mannose-specific beta-prism I fold lectins of nearly the same tertiary structure. MD simulations on individual subunits and the oligomers provide insights into the changes in the structure brought about in the protomers on oligomerization, including swapping of the N-terminal stretch in one instance. The regions that undergo changes also tend to exhibit dynamic flexibility during MD simulations. The internal symmetries of individual oligomers are substantially retained during the calculations. Energy minimization and simulations were also carried out on models using all possible oligomers by employing the four different protomers. The unique dimerization pattern observed in calsepa could be traced to unique substitutions in a peptide stretch involved in dimerization. The impossibility of a specific mode of oligomerization involving a particular protomer is often expressed in terms of unacceptable steric contacts or dissociation of the oligomer during simulations. The calculations also led to a rationale for the observation of a heltuba tetramer in solution although the lectin exists as an octamer in the crystal, in addition to providing insights into relations among evolution, oligomerization and ligand binding.     

Spatial and Temporal Simulation of Human Evolution. Methods,Frameworks and Applications

Analyses of human evolution are fundamental to understand the current gradients of human diversity. In this concern, genetic samples collected from current populations together with archaeological data are the most important resources to study human evolution. However, they are often insufficient to properly evaluate a variety of evolutionary scenarios, leading to continuous debates and discussions. A commonly applied strategy consists of the use of computer simulations based on, as realistic as possible, evolutionary models, to evaluate alternative evolutionary scenarios through statistical correlations with the real data. Computer simulations can also be applied to estimate evolutionary parameters or to study the role of each parameter on the evolutionary process. Here we review the mainly used methods and evolutionary frameworks to perform realistic spatially explicit computer simulations of human evolution. Although we focus on human evolution, most of the methods and software we describe can also be used to study other species. We also describe the importance of considering spatially explicit models to better mimic human evolutionary scenarios based on a variety of phenomena such as range expansions, range shifts, range contractions, sex-biased dispersal, long-distance dispersal or admixtures of populations. We finally discuss future implementations to improve current spatially explicit simulations and their derived applications in human evolution.     

Molecular dynamics simulations of a double unit cell in a protein crystal: volume relaxation at constant pressure and correlation of motions between the two unit cells

Eight molecular dynamics simulations of a double crystal unit cell of ubiquitin were performed to investigate the effects of simulating at constant pressure and of simulating two unit cells compared to a single unit cell. To examine the influence of different simulation conditions, the constant-pressure and constant-volume simulations were each performed with and without counterions and using two different treatments of the long-range electrostatic interactions (lattice-sum and reaction-field methods). The constant-pressure simulations were analyzed in terms of unit cell deformation and accompanying protein deformations. Energetic and structural properties of the proteins in the simulations of the double unit cell were compared to the results of previously reported one-unit-cell simulations. Correlation between the two unit cells was also investigated based on relative translational and rotational movements of the proteins and on dipole fluctuations. The box in the constant-pressure simulations is found to deform slowly to reach convergence only after 5-10 ns. This deformation does not result from a distortion in the structure of the proteins but rather from changes in protein packing within the unit cell. The results of the double-unit-cell simulations are closely similar to the results of the single-unit-cell simulations, and little motional correlation is found between the two unit cells.     

Event detection and sub‐state discovery from biomolecular simulations using higher‐order statistics: Application to enzyme adenylate kinase

Biomolecular simulations at millisecond and longer time‐scales can provide vital insights into functional mechanisms. Because post‐simulation analyses of such large trajectory datasets can be a limiting factor in obtaining biological insights, there is an emerging need to identify key dynamical events and relating these events to the biological function online, that is, as simulations are progressing. Recently, we have introduced a novel computational technique, quasi‐anharmonic analysis (QAA) (Ramanathan et al., PLoS One 2011;6:e15827), for partitioning the conformational landscape into a hierarchy of functionally relevant sub‐states. The unique capabilities of QAA are enabled by exploiting anharmonicity in the form of fourth‐order statistics for characterizing atomic fluctuations. In this article, we extend QAA for analyzing long time‐scale simulations online. In particular, we present HOST4MD—a higher‐order statistical toolbox for molecular dynamics simulations, which (1) identifies key dynamical events as simulations are in progress, (2) explores potential sub‐states, and (3) identifies conformational transitions that enable the protein to access those sub‐states. We demonstrate HOST4MD on microsecond timescale simulations of the enzyme adenylate kinase in its apo state. HOST4MD identifies several conformational events in these simulations, revealing how the intrinsic coupling between the three subdomains (LID, CORE, and NMP) changes during the simulations. Further, it also identifies an inherent asymmetry in the opening/closing of the two binding sites. We anticipate that HOST4MD will provide a powerful and extensible framework for detecting biophysically relevant conformational coordinates from long time‐scale simulations. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

Effect of lipid composition on buforin II structure and membrane entry

Buforin II is a 21-amino acid polycationic antimicrobial peptide derived from a peptide originally isolated from the stomach tissue of the Asian toad Bufo bufo gargarizans. It is hypothesized to target a wide range of bacteria by translocating into cells without membrane permeabilization and binding to nucleic acids. Previous research found that the structure and membrane interactions of buforin II are related to lipid composition. In this study, we used molecular dynamics (MD) simulations along with lipid vesicle experiments to gain insight into how buforin II interacts differently with phosphatidylcholine (PC), phosphatidylglycerol (PG), and phosphatidylethanolamine (PE) lipids. Fluorescent spectroscopic measurements agreed with the previous assertion that buforin II does not interact with pure PC vesicles. Nonetheless, the reduced entry of the peptide into anionic PG membranes versus neutral PC membranes during simulations correlates with the experimentally observed reduction in BF2 translocation through pure PG membranes. Simulations showing membrane entry into PC also provide insight into how buforin II may initially penetrate cell membranes. Our MD simulations also allowed us to consider how neutral PE lipids affect the peptide differently than PC. In particular, the peptide had a more helical secondary structure in simulations with PE lipids. A change in structure was also apparent in circular dichroism measurements. PE also reduced membrane entry in simulations, which correlates with decreased translocation in the presence of PE observed in previous studies. Together, these results provide molecular-level insight into how lipid composition can affect buforin II structure and function and will be useful in efforts to design peptides with desired antimicrobial and cell-penetrating properties.     

Effects of wall distensibility in hemodynamic simulations of an arteriovenous fistula

Arteriovenous fistulae are created surgically to provide adequate access for dialysis patients suffering from end-stage renal disease. It has long been hypothesized that the rapid blood vessel remodeling occurring after fistula creation is in part a process to restore the mechanical stresses to some preferred level, i.e., mechanical homeostasis. The current study presents fluid–structure interaction (FSI) simulations of a patient-specific model of a mature arteriovenous fistula reconstructed from 3D ultrasound scans. The FSI results are compared with previously published data of the same model but with rigid walls. Ultrasound-derived wall motion measurements are also used to validate the FSI simulations of the wall motion. Very large time-averaged shear stresses, 10–15 Pa, are calculated at the fistula anastomosis in the FSI simulations, values which are much larger than what is typically thought to be the normal homeostatic shear stress in the peripheral vasculature. Although this result is systematically lower by as much as 50 % compared to the analogous rigid-walled simulations, the inclusion of distensible vessel walls in hemodynamic simulations does not reduce the high anastomotic shear stresses to “normal” values. Therefore, rigid-walled analyses may be acceptable for identifying high shear regions of arteriovenous fistulae.     

Modeling simultaneous glucose and xylose uptake in Saccharomyces cerevisiae from kinetics and gene expression of sugar transporters

A kinetic model for glucose and xylose co-substrate uptake in Saccharomyces cerevisiae is presented. The model couples the enzyme kinetics with the glucose-dependent genetic expression of the individual transport proteins. This novel approach implies several options for optimizing the co-substrate utilization. Interestingly, the simulations predict a maximum xylose uptake rate at a glucose concentration >0 g/L, which suggests that the genetic expressions of the considered transport proteins are of importance when optimizing the xylose uptake. This was also evident in fed-batch simulations, where a distinct optimal glucose addition rate >0 g/L x h was found. Strategies for improving the co-substrate utilization by genetic engineering of the transport systems are furthermore suggested based on simulations.     

Comparison of the observed with the simulated distributions of the parental genome contribution in two marker-assisted backcross programs in rice

Computer simulations are useful tools to optimize marker-assisted breeding programs. The objective of our study was to investigate the closeness of computer simulations of the recurrent parent genome recovery with experimental data obtained in two marker-assisted backcrossing programs in rice (Orzya sativa L.). We simulated the breeding programs as they were practically carried out. In the simulations we estimated the frequency distributions of the recurrent parent genome proportion in the backcross populations. The simulated distributions were in good agreement with those obtained practically. The simulation results were also observed to be robust with respect to the choice of the mapping function and the accuracy of the linkage map. We conclude that computer simulations are a useful tool for pre-experiment estimation of selection response in marker-assisted backcrossing. Vanessa Prigge and Hans Peter Maurer contributed equally to this work.     

Insights into prion protein function from atomistic simulations

Computer simulations are a powerful tool for studies of biological systems. They have often been used to study prion protein (PrP), a protein responsible for neurodegenerative diseases, which include “mad cow disease” in cattle and Creutzfeldt-Jacob disease in humans. An important aspect of the prion protein is its interaction with copper ion, which is thought to be relevant for PrP’s yet undetermined function and also potentially play a role in prion diseases. For studies of copper attachment to the prion protein, computer simulations have often been used to complement experimental data and to obtain binding structures of Cu-PrP complexes. This paper summarizes the results of recent ab initio calculations of copper-prion protein interactions focusing on the recently discovered concentration-dependent binding modes in the octarepeat region of this protein. In addition to determining the binding structures, computer simulations were also used to make predictions about PrP’s function and the role of copper in prion diseases. The results demonstrate the predictive power and applicability of ab initio simulations for studies of metal-biomolecular complexes.     

Computer simulations of de novo designed helical proteins.

In the context of reduced protein models, Monte Carlo simulations of three de novo designed helical proteins (four-member helical bundle) were performed. At low temperatures, for all proteins under consideration, protein-like folds having different topologies were obtained from random starting conformations. These simulations are consistent with experimental evidence indicating that these de novo designed proteins have the features of a molten globule state. The results of Monte Carlo simulations suggest that these molecules adopt four-helix bundle topologies. They also give insight into the possible mechanism of folding and association, which occurs in these simulations by on-site assembly of the helices. The low-temperature conformations of all three sequences have the features of a molten globule state.     

Molecular dynamics of water in the neighborhood of aquaporins

Present knowledge obtained by molecular dynamics (MD) simulation studies regarding the dynamics of water, both in the vicinity of biological membranes and within the proteinaceous water channels, also known as aquaporins (AQPs), is reviewed. A brief general summary of the water models most extensively employed in MD simulations (SPC, SPC/E, TIP3P, TIP4P), indicating their most relevant pros and cons, is likewise provided. Structural considerations of water are also discussed, based on different order parameters, which can be extracted from MD simulations as well as from experiments. Secondly, the behaviour of water in the neighbourhood of membranes by means of molecular dynamics simulations is addressed. Consequently, the comparison with previous experimental evidence is pointed out. In living cells, water is transported across the plasma membrane through the lipid bilayer and the aforementioned AQPs, which motivates this review to focus mostly on MD simulation studies of water within AQPs. Relevant contributions explaining peculiar properties of these channels are discussed, such as selectivity and gating. Water models used in these studies are also summarised. Finally, based on the information presented here, further MD studies are encouraged.