Brownian dynamics simulation of peptides with the University of Houston Brownian Dynamics (UHBD) program

This chapter provides the background theory and a practical protocol for performing Brownian dynamics simulation of peptides. Brownian dynamics simulation represents a complementary approach to Monte Carlo and molecular dynamics methods. Unlike Monte Carlo methods, it could provide dynamical information in a timescale longer than the momentum relaxation time. On the other hand, it is faster than molecular dynamics by approximating the solvent by a continuum and by operating in the over-damped limit. This chapter introduces the use of the University of Houston Brownian Dynamics (UHBD) program [1, 2] to perform Brownian dynamics simulation on peptides.

β-Barrel topology of Alzheimer's β-amyloid ion channels

Emerging evidence supports the ion channel mechanism for Alzheimer's disease pathophysiology wherein small β-amyloid (Aβ) oligomers insert into the cell membrane, forming toxic ion channels and destabilizing the cellular ionic homeostasis. Solid-state NMR-based data of amyloid oligomers in solution indicate that they consist of a double-layered β-sheets where each monomer folds into β-strand-turn-β-strand and the monomers are stacked atop each other. In the membrane, Aβ peptides are proposed to be β-type structures. Experimental structural data available from atomic force microscopy (AFM) imaging of Aβ oligomers in membranes reveal heterogeneous channel morphologies. Previously, we modeled the channels in a non-tilted organization, parallel with the cross-membrane normal. Here, we modeled a β-barrel-like organization. β-Barrels are common in transmembrane toxin pores, typically consisting of a monomeric chain forming a pore, organized in a single-layered β-sheet with antiparallel β-strands and a right-handed twist. Our explicit solvent molecular dynamics simulations of a range of channel sizes and polymorphic turns and comparisons of these with AFM image dimensions support a β-barrel channel organization. Different from the transmembrane β-barrels where the monomers are folded into a circular β-sheet with antiparallel β-strands stabilized by the connecting loops, these Aβ barrels consist of multimeric chains forming double β-sheets with parallel β-strands, where the strands of each monomer are connected by a turn. Although the Aβ barrels adopt the right-handed β-sheet twist, the barrels still break into heterogeneous, loosely attached subunits, in good agreement with AFM images and previous modeling. The subunits appear mobile, allowing unregulated, hence toxic, ion flux.

Computational simulations of the Trichoderma reesei cellobiohydrolase I acting on microcrystalline cellulose Ibeta: the enzyme-substrate complex

Cellobiohydrolases are the dominant components of the commercially relevant Trichoderma reesei cellulase system. Although natural cellulases can totally hydrolyze crystalline cellulose to soluble sugars, the current enzyme loadings and long digestion times required render these enzymes less than cost effective for biomass conversion processes. It is clear that cellobiohydrolases must be improved via protein engineering to reduce processing costs. To better understand cellobiohydrolase function, new simulations have been conducted using charmm of cellobiohydrolase I (CBH I) from T.reesei interacting with a model segment (cellodextrin) of a cellulose microfibril in which one chain from the substrate has been placed into the active site tunnel mimicking the hypothesized configuration prior to final substrate docking (i.e., the +1 and +2 sites are unoccupied), which is also the structure following a catalytic bond scission. No tendency was found for the protein to dissociate from or translate along the substrate surface during this initial simulation, nor to align with the direction of the cellulose chains. However, a tendency for the decrystallized cellodextrin to partially re-anneal into the cellulose surface hints that the arbitrary starting configuration selected was not ideal.

How to lose a kink and gain a helix: pH independent conformational changes of the fusion domains from influenza hemagglutinin in heterogeneous lipid bilayers

We have simulated two conformations of the fusion domain of influenza hemagglutinin (HA) within explicit water, salt, and heterogeneous lipid bilayers composed of POPC:POPG (4:1). Each conformation has seven different starting points in which the initial peptide structure is the same for each conformation, but the location across the membrane normal and lipid arrangement around the peptide are varied, giving a combined total simulation time of 140 ns. For the HA5 conformation (primary structure from recent NMR spectroscopy at pH = 5), the peptide exhibits a stable and less kinked structure in the lipid bilayer compared to that from the NMR studies. The relative fusogenic behavior of the different conformations has been investigated by calculation of the relative free energy of insertion into the hydrophobic region of lipid bilayer as a function of the depth of immersion. For the HA7 conformations (primary structure from recent NMR spectroscopy at pH = 7.4), while the N-terminal helix preserves its initial structure, the flexible C-terminal chain produces a transient helical motif inside the lipid bilayer. This conformational change is pH-independent, and is closely related to the peptide insertion into the lipid bilayer.

Insight into the role of hydration on protein dynamics

The potential energy surface of a protein is rough. This intrinsic energetic roughness affects diffusion, and hence the kinetics. The dynamics of a system undergoing Brownian motion on this surface in an implicit continuum solvent simulation can be tuned via the frictional drag or collision frequency to be comparable to that of experiments or explicit solvent simulations. We show that the kinetic rate constant for a local rotational isomerization in stochastic simulations with continuum solvent and a collision frequency of 2 ps(-1) is about 10(4) times faster than that in explicit water and experiments. A further increase in the collision frequency to 60 ps(-1) slows down the dynamics, but does not fully compensate for the lack of explicit water. We also show that the addition of explicit water does not only slow down the dynamics by increasing the frictional drag, but also increases the local energetic roughness of the energy landscape by as much as 1.0 kcal/mol.

Optimized atomic radii for protein continuum electrostatics solvation forces

Recently, we presented a Green's function approach for the calculation of analytic continuum electrostatic solvation forces based on numerical solutions of the finite-difference Poisson-Botzmann (FDPB) equation [Im et al., Comp. Phys. Comm. 111 (1998) 59]. In this treatment the analytic forces were explicitly defined as the first derivative of the FDPB continuum electrostatic free energy with respect to the coordinates of the solute atoms. A smooth intermediate region for the solute-solvent dielectric boundary needed to be introduced to avoid abrupt discontinuous variations in the solvation free energy and forces as a function of the atomic positions. In the present paper we extend the set of optimized radii, which was previously parametrized from molecular dynamics free energy simulations of the 20 standard amino acids with explicit solvent molecules [Nina et al., J. Phys. Chem. 101 (1997) 5239], to yield accurate solvation free energy by taking the influence of the smoothed dielectric region into account.

On Use of the Amber Potential with the Langevin Dipole Method

Inclusion of solvent effects in biomolecular simulations is most ideally done using explicit methods, as they are able to capture the heterogeneous environment typical of biomolecules and systems involving them (e.g., proteins at solid interfaces). Common explicit methods based on molecular solvent models (e.g., TIP and SPC models) and molecular dynamic or Monte Carlo simulation are computationally expensive and are, therefore, not well-suited to situations where many simulations are required (e.g., in the ab initio structure prediction or design contexts). In such cases, more coarse-grained explicit approaches such as the Langevin dipole (LD) method of Warshel and co-workers are more appropriate. The recent incarnations of the LD method appear to produce good solvation free energy estimates. These incarnations use charges and solute structures obtained from high-level quantum mechanics simulations. As such an approach is clearly not possible for larger solutes or when many structures are to be considered, an alternative must be sought. One possibility is to use structures and charges derived from an existing analytical potential model-we report on such a coupling here with the Amber potential model. The accuracy and computational performance of this hybrid approach, which we term LD-Amber to distinguish it from previous incarnations of the LD method, was assessed by comparing results obtained from the approach with those from experiment and other theoretical methods for the solvation of 18 amino acid analogues and the alanine dipeptide. This comparison shows that the LD-Amber approach can yield results in line with experiment both qualitatively and quantitatively and is as accurate as other explicit methods while being computationally much cheaper.

Conformational Preferences of Non-Prolyl and Prolyl Residues

The conformational study on Ac-Ala-NHMe (the alanine dipeptide) and Ac-Pro-NHMe (the proline dipeptide) is carried out using ab initio HF and density functional methods with the self-consistent reaction field method to explore the differences in the backbone conformational preference and the cis-trans isomerization for the non-prolyl and prolyl residues in the gas phase and in the solutions (chloroform and water). For the alanine and proline dipeptides, with the increase of solvent polarity, the populations of the conformation tC with an intramolecular C(7) hydrogen bond significantly decrease, and those of the polyproline II-like conformation tF and the alpha-helical conformation tA increase, which is in good agreement with the results from circular dichroism and NMR experiments. For both the dipeptides, as the solvent polarity increases, the relative free energy of the cis conformer to the trans conformer decreases and the rotational barrier to the cis-trans isomerization increases. It is found that the cis-trans isomerization proceeds in common through only the clockwise rotation with omega' approximately +120 degrees about the non-prolyl and prolyl peptide bonds in both the gas phase and the solutions. The pertinent distance d(N...H-N(NHMe)) can successfully describe the increase in the rotational barriers for the non-prolyl and prolyl trans-cis isomerization as the solvent polarity increases and the higher barriers for the non-prolyl residue than for the prolyl residue, as seen in experimental and calculated results. By analysis of the contributions to rotational barriers, the cis-trans isomerization for the non-prolyl and prolyl peptide bonds is proven to be entirely enthalpy driven in the gas phase and in the solutions. The calculated cis populations and rotational barriers to the cis-trans isomerization for both the dipeptides in chloroform and/or water accord with the experimental values.

Assessing implicit models for nonpolar mean solvation forces: the importance of dispersion and volume terms

Continuum solvation models provide appealing alternatives to explicit solvent methods because of their ability to reproduce solvation effects while alleviating the need for expensive sampling. Our previous work has demonstrated that Poisson-Boltzmann methods are capable of faithfully reproducing polar explicit solvent forces for dilute protein systems; however, the popular solvent-accessible surface area model was shown to be incapable of accurately describing nonpolar solvation forces at atomic-length scales. Therefore, alternate continuum methods are needed to reproduce nonpolar interactions at the atomic scale. In the present work, we address this issue by supplementing the solvent-accessible surface area model with additional volume and dispersion integral terms suggested by scaled particle models and Weeks-Chandler-Andersen theory, respectively. This more complete nonpolar implicit solvent model shows very good agreement with explicit solvent results and suggests that, although often overlooked, the inclusion of appropriate dispersion and volume terms are essential for an accurate implicit solvent description of atomic-scale nonpolar forces.

Multiple pathways in conformational transitions of the alanine dipeptide: An application of dynamic importance sampling

We present multiple dynamic transition pathways on the two-dimensional dihedral plane between conformational states of the alanine dipeptide. The method used in this study is dynamic importance sampling (DIMS). To perform DIMS, unbiased molecular dynamic simulations are used to generate equilibrium ensembles for the alanine dipeptide within different states. Free energy surfaces on the dihedral plane are calculated from the equilibrium simulations, and four energy minima defined from the surface are used as the starting and ending points for DIMS dynamics. The DIMS method represents an important step towards finding multiple transition pathways within complex biomolecular systems.

A novel method reveals that solvent water favors polyproline II over beta-strand conformation in peptides and unfolded proteins: conditional hydrophobic accessible surface area (CHASA)

In aqueous solution, the ensemble of conformations sampled by peptides and unfolded proteins is largely determined by their interaction with water. It has been a long-standing goal to capture these solute-water energetics accurately and efficiently in calculations. Historically, accessible surface area (ASA) has been used to estimate these energies, but this method breaks down when applied to amphipathic peptides and proteins. Here we introduce a novel method in which hydrophobic ASA is determined after first positioning water oxygens in hydrogen-bonded orientations proximate to all accessible peptide/protein backbone N and O atoms. This conditional hydrophobic accessible surface area is termed CHASA. The CHASA method was validated by predicting the polyproline-II (P(II)) and beta-strand conformational preferences of non-proline residues in the coil library (i.e., non-alpha-helix, non-beta-strand, non-beta-turn library derived from X-ray elucidated structures). Further, the method successfully rationalizes the previously unexplained solvation energies in polyalanyl peptides and compares favorably with published experimentally determined P(II) residue propensities. We dedicate this paper to Frederic M. Richards.

Role of solvent in determining conformational preferences of alanine dipeptide in water

Evidence from a variety of spectroscopic probes indicates that (phi, psi) values corresponding to the left-handed polyproline II helix (P(II)) are preferred for short alanine-based peptides in water. On the basis of results from theoretical studies, it is believed that the observed preference is dictated by favorable peptide-solvent interactions, which are realized through formation of optimal hydrogen-bonding water bridges between peptide donor and acceptor groups. In the present study, we address this issue explicitly by analyzing the hydration structure and thermodynamics of 16 low-energy conformers of the alanine dipeptide (N-acetylalanine-N'-methylamide) in liquid water. Monte Carlo simulations in the canonical ensemble were performed under ambient conditions with all-atom OPLS parameters for the alanine dipeptide and the TIP5P model for water. We find that the number of hydrogen-bonded water molecules connecting the peptide group donor and acceptor atoms has no effect on the solvation thermodynamics. Instead, the latter are determined by the work done to fully hydrate the peptide. This work is minimal for conformations that are characterized by a minimal overlap of the primary hydration shells around the peptide donor and acceptor atoms. As a result, peptide-solvent interactions favor "compact" conformations that do not include P(II)-like geometries. Our main conclusion is that the experimentally observed preference for P(II) does not arise due to favorable direct interactions between the peptide and water molecules. Instead, the latter act to unmask underlying conformational preferences that are a consequence of minimizing intrapeptide steric conflicts.

Brownian dynamics simulation of helix-capping motifs

Helix-capping motifs are believed to play an important role in stabilizing alpha-helices and defining helix start and stop signals. We performed microsecond scale Brownian dynamics simulations to study ten XAAD sequences, with X = (A,E,I,L,N,Q,S,T,V,Y), to examine their propensity to form helix capping motifs and correlate these results with those obtained from analyzing a structural database of proteins. For the widely studied capping box motif S**D, where the asterisk can be any amino acid residue, the simulations suggested that one of the two hydrogen bonds proposed earlier as a stabilizing factor might not be as important. On the other hand, side-chain interactions between the capping residue and the third residue downstream on the polypeptide chain might also play a role in stabilizing this motif. These results are consistent with explicit-solvent molecular dynamics simulations of two capping box motifs found in the proteins BPTI and alpha-dendrotoxin. Principal component analysis of the SAAD trajectory showed that the first three principal components, after those corresponding to translational-rotational motion were removed, accounted for more than half of the conformational fluctuations. The first component separated the conformational space into two parts with the all-helical conformation and the capping box motif lying largely in one part. The second component, on the other hand, could be used to describe conformational transitions between the all-helical form and the capping box motif.

Protocol for MM/PBSA molecular dynamics simulations of proteins

Continuum solvent models have been employed in past years for understanding processes such as protein folding or biomolecular association. In the last decade, several attempts have been made to merge atomic detail molecular dynamics simulations with solvent continuum models. Among continuum models, the Poisson-Boltzmann solvent accessible surface area model is one of the oldest and most fundamental. Notwithstanding its wide usage for simulation of biomolecular electrostatic potential, the Poisson-Boltzmann equation has been very seldom used to obtain solvation forces for molecular dynamics simulation. We propose here a fast and reliable methodology to implement continuum forces in standard molecular mechanics and dynamics algorithms. Results for a totally unrestrained 1 ns molecular dynamics simulation of a small protein are quantitatively similar to results obtained by explicit solvent molecular dynamics simulations.

Comparison of various implicit solvent models in molecular dynamics simulations of immunoglobulin G light chain dimer

The present study tests performance of different solvation models applied to molecular dynamics simulation of a large, dimeric protein molecule. Analytical Continuum Electrostatics (ACE) with two different parameter sets, older V98 and new V01, and Effective Energy Function (EEF) are employed in molecular dynamics simulation of immunoglobulin G (IgG) light chain dimer and variable domain of IgG light chain. Results are compared with explicit solvent and distance dependent dielectric constant (DDE) calculations. The overall analysis shows that the EEF method yields results comparable to explicit solvent simulations; however, the stability of simulations is lower. On the other hand, the ACE_V98 model does not seem to achieve the accuracy or stability expected in nanosecond timescale MD simulation for the studied systems. The ACE_V01 model greatly improves stability of the calculation; nonetheless, changes in radius of gyration and solvent accessible surface of the studied systems may indicate that the parameter set still needs to be improved if the method is supposed to be used for simulations of large, polymeric proteins. Additionally, electrostatic contribution to the solvation free energy calculated in the ACE model is compared with a numerical treatment of the dielectric continuum model. Wall clock time of all simulations is compared. It shows that EEF calculation is six times faster than corresponding ACE and 50 times faster than explicit solvent simulations.

The Poisson-Boltzmann equation for biomolecular electrostatics: a tool for structural biology

Electrostatics plays a fundamental role in virtually all processes involving biomolecules in solution. The Poisson-Boltzmann equation constitutes one of the most fundamental approaches to treat electrostatic effects in solution. The theoretical basis of the Poisson-Boltzmann equation is reviewed and a wide range of applications is presented, including the computation of the electrostatic potential at the solvent-accessible molecular surface, the computation of encounter rates between molecules in solution, the computation of the free energy of association and its salt dependence, the study of pKa shifts and the combination with classical molecular mechanics and dynamics. Theoretical results may be used for rationalizing or predicting experimental results, or for suggesting working hypotheses. An ever-increasing body of successful applications proves that the Poisson-Boltzmann equation is a useful tool for structural biology and complementary to other established experimental and theoretical methodologies.

Protein molecular dynamics with electrostatic force entirely determined by a single Poisson-Boltzmann calculation

The electrostatic force including the intramolecular Coulombic interactions and the electrostatic contribution of solvation effect were entirely calculated by using the finite difference Poisson-Boltzmann method (FDPB), which was incorporated into the GROMOS96 force field to complete a new finite difference stochastic dynamics procedure (FDSD). Simulations were performed on an insulin dimer. Different relative dielectric constants were successively assigned to the protein interior; a value of 17 was selected as optimal for our system. The simulation data were analyzed and compared with those obtained from 500-ps molecular dynamics (MD) simulation with explicit water and a 500-ps conventional stochastic dynamics (SD) simulation without the mean solvent force. The results indicate that the FDSD method with GROMOS96 force field is suitable to study the dynamics and structure of proteins in solution if used with the optimal protein dielectric constant.

New approach to free energy of solvation applying continuum models to molecular dynamics simulation

A new approach to the calculation of the free energy of solvation from trajectories obtained by molecular dynamics simulation is presented. The free energy of solvation is computed as the sum of three contributions originated at the cavitation of the solute by the solvent, the solute-solvent nonpolar (repulsion and dispersion) interactions, and the electrostatic solvation of the solute. The electrostatic term is calculated based on ideas developed for the broadly used continuum models, the cavitational contribution from the excluded volume by the Claverie-Pierotti model, and the Van der Waals term directly from the molecular dynamics simulation. The proposed model is tested for diluted aqueous solutions of simple molecules containing a variety of chemically important functions: methanol, methylamine, water, methanethiol, and dichloromethane. These solutions were treated by molecular dynamics simulations using SPC/E water and the OPLS force field for the organic molecules. Obtained free energies of solvation are in very good agreement with experimental data.

A critical analysis of continuum electrostatics: the screened Coulomb potential--implicit solvent model and the study of the alanine dipeptide and discrimination of misfolded structures of proteins

An analysis of the screened Coulomb potential--implicit solvent model (SCP--ISM) is presented showing that general equations for both the electrostatic and solvation free energy can be derived in a continuum approach, using statistical averaging of the polarization field created by the solvent around the molecule. The derivation clearly shows how the concept of boundary, usually found in macroscopic approaches, is eliminated when the continuum model is obtained from a microscopic treatment using appropriate averaging techniques. The model is used to study the alanine dipeptide in aqueous solution, as well as the discrimination of native protein structures from misfolded conformations. For the alanine dipeptide the free energy surface in the phi--psi space is calculated and compared with recently reported results of a detailed molecular dynamics simulation using an explicit representation of the solvent, and with other available data. The study showed that the results obtained using the SCP--ISM are comparable to those of the explicit water calculation and compares favorably to the FDPB approach. Both transition states and energy minima show a high correlation (r > 0.98) with the results obtained in the explicit water analysis. The study of the misfolded structures of proteins comprised the analysis of three standard decoy sets, namely, the EMBL, Park and Levitt, and Baker's CASP3 sets. In all cases the SCP--ISM discriminated well the native structures of the proteins, and the best-predicted structures were always near-native (cRMSD approximately 2 A).

Evaluation of a fast implicit solvent model for molecular dynamics simulations

A solvation term based on the solvent accessible surface area (SASA) is combined with the CHARMM polar hydrogen force field for the efficient simulation of peptides and small proteins in aqueous solution. Only two atomic solvation parameters are used: one is negative for favoring the direct solvation of polar groups and the other positive for taking into account the hydrophobic effect on apolar groups. To approximate the water screening effects on the intrasolute electrostatic interactions, a distance-dependent dielectric function is used and ionic side chains are neutralized. The use of an analytical approximation of the SASA renders the model extremely efficient (i.e., only about 50% slower than in vacuo simulations). The limitations and range of applicability of the SASA model are assessed by simulations of proteins and structured peptides. For the latter, the present study and results reported elsewhere show that with the SASA model it is possible to sample a significant amount of folding/unfolding transitions, which permit the study of the thermodynamics and kinetics of folding at an atomic level of detail.

Molecular mechanics and dynamics of biomolecules using a solvent continuum model

An easy implementation of molecular mechanics and molecular dynamics simulation using a continuum solvent model is presented that is particularly suitable for biomolecular simulations. The computation of solvation forces is made using the linear Poisson-Boltzmann equation (polar contribution) and the solvent-accessible surface area approach (nonpolar contribution). The feasibility of the methodology is demonstrated on a small protein and a small DNA hairpin. Although the parameters employed in this model must be refined to gain reliability, the performance of the method, with a standard choice of parameters, is comparable with results obtained by explicit water simulations. Copyright 2001 John Wiley & Sons, Inc. J Comput Chem 22: 1830-1842, 2001

Atomistic Brownian dynamics simulation of peptide phosphorylation

We report the implementation of an all-atom Brownian dynamics simulation model of peptides using the constraint algorithm LINCS. The algorithm has been added as a part of UHBD. It uses adaptive time steps to achieve a balance between computational speed and stability. The algorithm was applied to study the effect of phosphorylation on the conformational preference of the peptide Gly-Ser-Ser-Ser. We find that the middle serine residue experiences considerable conformational change from the C(7eq) to the alpha(R) structure upon phosphorylation. NMR (3)J coupling constants were also computed from the Brownian trajectories using the Karplus equation. The calculated (3)J results agree reasonably well with experimental data for phosphorylated peptide but less so for doubly charged phosphorylated one.

Methodology for protein-ligand binding studies: application to a model for drug resistance, the HIV/FIV protease system

A protocol for the rapid energetic analysis of protein-ligand complexes has been developed. This protocol involves the generation of protein-ligand complex ensembles followed by an analysis of the binding free energy components. We apply this methodology toward understanding the origin of binding specificity within the human immunodeficiency virus/feline immunodeficiency virus (HIV/FIV) protease system, a model system for drug resistance studies. A distinct difference in the internal strain of an inhibitor within each protein environment clearly favors the HIV protease complex, as observed experimentally. Our analysis also predicts that residues within the S2-S3 pockets of the FIV protease active site are responsible for this strain. Close examination of the active site residue contributions to interaction energy and desolvation energy identifies specific amino acids that may also play a role in determining the binding preferences of these two enzymes. Proteins 1999;36:318-331.

Unrestrained stochastic dynamics simulations of the UUCG tetraloop using an implicit solvation model

Three unrestrained stochastic dynamics simulations have been carried out on the RNA hairpin GGAC[UUCG] GUCC, using the AMBER94 force field (Cornell et al., 1995. J. Am. Chem. Soc. 117:5179-5197) in MacroModel 5.5 (Mohamadi et al., 1990. J. Comp. Chem. 11:440-467) and either the GB/SA continuum solvation model (Still et al., 1990. J. Am. Chem. Soc. 112:6127-6129) or a linear distance-dependent dielectric (1/R) treatment. The linear distance-dependent treatment results in severe distortion of the nucleic acid structure, restriction of all hydroxyl dihedrals, and collapse of the counterion atmosphere over the course of a 5-ns simulation. An additional vacuum simulation without counterions shows somewhat improved behavior. In contrast, the two GB/SA simulations (1.149 and 3.060 ns in length) give average structures within 1.2 A of the initial NMR structure and in excellent agreement with results of an earlier explicit solvent simulation (Miller and Kollman, 1997. J. Mol. Biol. 270:436-450). In a 3-ns GB/SA simulation starting with the incorrect UUCG tetraloop structure (Cheong et al., 1990. Nature. 346:680-682), this loop conformation converts to the correct loop geometry (Allain and Varani, 1995. J. Mol. Biol. 250:333-353), suggesting enhanced sampling relative to the previous explicit solvent simulation. Thermodynamic effects of 2'-deoxyribose substitutions of loop nucleotides were experimentally determined and are found to correlate with the fraction of time the ribose 2'-OH is hydrogen bonded and the distribution of the hydroxyl dihedral is observed in the GB/SA simulations. The GB/SA simulations thus appear to faithfully represent structural features of the RNA without the computational expense of explicit solvent.

Phosphorylation stabilizes the N-termini of alpha-helices

The role of phosphorylation in stabilizing the N-termini of alpha-helices is examined using computer simulations of model peptides. The models comprise either a phosphorylated or unphosphorylated serine at the helix N-terminus, followed by nine alanines. Monte Carlo/stochastic Dynamics simulations were performed on the model helices. The simulations revealed a distinct stabilization of the helical conformation at the N-terminus after phosphorylation. The stabilization was attributable to favorable electrostatic interactions between the phosphate and the helix backbone. However, direct helix capping by the phosphorylated sidechain was not observed. The results of the calculations are consistent with experimental evidence on the stabilization of helices by phosphates and other anions.