CiteULike: bilab's library [52 articles]

The Solvation Interface is a Determining Factor in Peptide Conformational Preferences

Eric Sorin — 2006-05-10T08:17:56-00:00

Journal of Molecular Biology, Vol. 356, No. 1. (10 February 2006), pp. 248-256.

The 21 residue polyalanine-based Fs peptide was studied using thousands of long, explicit solvent, atomistic molecular dynamics simulations that reached equilibrium at the ensemble level. Peptide conformational preference as a function of hydrophobicity was examined using a spectrum of explicit solvent models, and the peptide length-dependence of the hydrophilic and hydrophobic components of solvent-accessible surface area for several ideal conformational types was considered. Our results demonstrate how the character of the solvation interface induces several conformational preferences, including a decrease in mean helical content with increased hydrophilicity, which occurs predominantly through reduced nucleation tendency and, to a lesser extent, destabilization of helical propagation. Interestingly, an opposing effect occurs through increased propensity for 310-helix conformations, as well as increased polyproline structure. Our observations provide a framework for understanding previous reports of conformational preferences in polyalanine-based peptides including (i) terminal 310-helix prominence, (ii) low [pi]-helix propensity, (iii) increased polyproline conformations in short and unfolded peptides, and (iv) membrane helix stability in the presence and absence of water. These observations provide physical insight into the role of water in peptide conformational equilibria at the atomic level, and expand our view of the complexity of even the most "simple" of biopolymers. Whereas previous studies have focused predominantly on hydrophobic effects with respect to tertiary structure, this work highlights the need for consideration of such effects at the secondary structural level.

Sub-microsecond Protein Folding

Jan Kubelka — 2006-05-02T10:34:41-00:00

Journal of Molecular Biology, Vol. In Press, Corrected Proof

We have investigated the structure, equilibria, and folding kinetics of an engineered 35-residue subdomain of the chicken villin headpiece, an ultrafast-folding protein. Substitution of two buried lysine residues by norleucine residues stabilizes the protein by 1 kcal/mol and increases the folding rate sixfold, as measured by nanosecond laser T-jump. The folding rate at 300 K is (0.7 [mu]s)-1 with little or no temperature dependence, making this protein the first sub-microsecond folder, with a rate only twofold slower than the theoretically predicted speed limit. Using the 70 ns process to obtain the effective diffusion coefficient, the free energy barrier height is estimated from Kramers theory to be less than ~1 kcal/mol. X-ray crystallographic determination at 1 A resolution shows no significant change in structure compared to the single-norleucine-substituted molecule and suggests that the increased stability is electrostatic in origin. The ultrafast folding rate, very accurate X-ray structure, and small size make this engineered villin subdomain an ideal system for simulation by atomistic molecular dynamics with explicit solvent.

Reaction coordinates and rates from transition paths.

RB Best — 2006-02-16T20:16:29-00:00

Proc Natl Acad Sci U S A, Vol. 102, No. 19. (10 May 2005), pp. 6732-6737.

The molecular mechanism of a reaction in solution is reflected in its transition-state ensemble and transition paths. We use a Bayesian formula relating the equilibrium and transition-path ensembles to identify transition states, rank reaction coordinates, and estimate rate coefficients. We also introduce a variational procedure to optimize reaction coordinates. The theory is illustrated with applications to protein folding and the dipole reorientation of an ordered water chain inside a carbon nanotube. To describe the folding of a simple model of a three-helix bundle protein, we variationally optimize the weights of a projection onto the matrix of native and nonnative amino acid contacts. The resulting one-dimensional reaction coordinate captures the folding transition state, with formation and packing of helix 2 and 3 constituting the bottleneck for folding.

Exploring reaction pathways with transition path and umbrella sampling: application to methyl maltoside.

RJ Dimelow — 2006-05-02T07:26:32-00:00

J Chem Phys, Vol. 124, No. 11. (21 March 2006)

The transition path sampling (TPS) method is a powerful approach to study chemical reactions or transitional properties on complex potential energy landscapes. One of the main advantages of the method over potential of mean force methods is that reaction rates can be directly accessed without knowledge of the exact reaction coordinate. We have investigated the complementary nature of these two differing approaches, comparing transition path sampling with the weighted histogram analysis method to study a conformational change in a small model system. In this case study, the transition paths for a transition between two rotational conformers of a model disaccharide molecule, methyl beta-D-maltoside, were compared with a free energy surface constrained by the two commonly used glycosidic (phi,psi) torsional angles. The TPS method revealed a reaction channel that was not apparent from the potential of mean force method, and the suitability of phi and psi as reaction coordinates to describe the isomerization in vacuo was confirmed by examination of the transition path ensemble. Using both transition state theory and transition path sampling methods, the transition rate was estimated. We have estimated a characteristic time between transitions of approximately 160 ns for this rare isomerization event between the two conformations of the carbohydrate. We conclude that transition path sampling can extract subtle information about the dynamics not apparent from the potential of mean force method. However, in calculating the reaction rate, the transition path sampling method required 27.5 times the computational effort than was needed by the potential of mean force method.

In Silico Discovery of beta-Secretase Inhibitors.

D Huang — 2006-05-02T07:23:56-00:00

J Am Chem Soc, Vol. 128, No. 16. (26 April 2006), pp. 5436-5443.

Alzheimer's disease, the most common amyloid-associated disorder, accounts for the majority of the dementia diagnosed after the age of 60. The cleavage of the beta-amyloid precursor protein is initiated by beta-secretase (BACE-1), a membrane-bound aspartic protease, which has emerged as an important but difficult protein target. Here, an in silico screening approach consisting of fragment-based docking, ligand conformational search by a genetic algorithm, and evaluation of free energy of binding was used to identify low-molecular-weight inhibitors of BACE-1. More than 300 000 small molecules were docked and about 15 000 prioritized according to a linear interaction energy model with evaluation of solvation by continuum electrostatics. Eighty-eight compounds were tested in vitro, and 10 of them showed an IC(50) value lower than 100 muM in a BACE-1 enzymatic assay. Interestingly, the 10 active compounds shared a triazine scaffold. Moreover, four of them were active in an assay with mammalian cells (EC(50) < 20 muM), indicating that they are cell-permeable. Therefore, these triazine derivatives are very promising lead candidates for BACE-1 inhibition. The discoveries of this series and two other series of nonpeptidic BACE-1 inhibitors demonstrate the usefulness of our in silico high-throughput screening approach.

Network and graph analyses of folding free energy surfaces

Amedeo Caflisch — 2006-02-14T20:14:38-00:00

Current Opinion in Structural Biology, Vol. 16, No. 1. (February 2006), pp. 71-78.

Protein folding is governed by a complex free energy surface whose entropic contributions are relevant because of the large number of degrees of freedom involved. Such complexity, in particular the conformational heterogeneity of the denatured state, is hidden in projections onto one or two order parameters (e.g. fraction of native contacts and/or radius of gyration), which usually results in relatively smooth surfaces. Recent approaches borrowed from network and graph theory have yielded quantitative unprojected representations of the free energy surfaces of a [beta]-hairpin and a three-stranded [beta]-sheet peptide using equilibrium folding-unfolding molecular dynamics simulations. Interestingly, the network and graph analyses of these structured peptides have revealed a very heterogeneous denatured state ensemble. It includes high-enthalpy, high-entropy conformations with fluctuating non-native secondary structure, as well as low-enthalpy, low-entropy traps.

Coupled folding and binding with alpha-helix-forming molecular recognition elements.

CJ Oldfield — 2006-05-02T07:10:04-00:00

Biochemistry, Vol. 44, No. 37. (20 September 2005), pp. 12454-12470.

Many protein-protein and protein-nucleic acid interactions involve coupled folding and binding of at least one of the partners. Here, we propose a protein structural element or feature that mediates the binding events of initially disordered regions. This element consists of a short region that undergoes coupled binding and folding within a longer region of disorder. We call these features "molecular recognition elements" (MoREs). Examples of MoREs bound to their partners can be found in the alpha-helix, beta-strand, polyproline II helix, or irregular secondary structure conformations, and in various mixtures of the four structural forms. Here we describe an algorithm that identifies regions having propensities to become alpha-helix-forming molecular recognition elements (alpha-MoREs) based on a discriminant function that indicates such regions while giving a low false-positive error rate on a large collection of structured proteins. Application of this algorithm to databases of genomics and functionally annotated proteins indicates that alpha-MoREs are likely to play important roles protein-protein interactions involved in signaling events.

Conformational Diversity of Ligands Bound to Proteins

Gareth Stockwell — 2006-01-29T20:14:56-00:00

Journal of Molecular Biology, Vol. 356, No. 4. (3 March 2006), pp. 928-944.

The phenomenon of molecular recognition, which underpins almost all biological processes, is dynamic, complex and subtle. Establishing an interaction between a pair of molecules involves mutual structural rearrangements guided by a highly convoluted energy landscape, the accurate mapping of which continues to elude us. Increased understanding of the degree to which the conformational space of a ligand is restricted upon binding may have important implications for docking studies, structure refinement and for function prediction methods based on geometrical comparisons of ligands or their binding sites. Here, we present an analysis of the conformational variability exhibited by three of the most ubiquitous biological ligands in nature, ATP, NAD and FAD. First, we demonstrate qualitatively that these ligands bind to proteins in widely varying conformations, including several cases in which parts of the molecule assume energetically unfavourable orientations. Next, by comparing the distribution of bound ligand shapes with the set of all possible molecular conformations, we provide a quantitative assessment of previous observations that ligands tend to unfold when binding to proteins. We show that, while extended forms of ligands are indeed common in ligand-protein structures, instances of ligands in almost maximally compact arrangements can also be found. Thirdly, we compare the conformational variation in two sets of ligand molecules, those bound to homologous proteins, and those bound to unrelated proteins. Although most superfamilies bind ligands in a fairly conserved manner, we find several cases in which significant variation in ligand configuration is observed.

A method for localizing ligand binding pockets in protein structures.

F Glaser — 2006-01-20T20:29:53-00:00

Proteins, Vol. 62, No. 2. (1 February 2006), pp. 479-488.

The accurate identification of ligand binding sites in protein structures can be valuable in determining protein function. Once the binding site is known, it becomes easier to perform in silico and experimental procedures that may allow the ligand type and the protein function to be determined. For example, binding pocket shape analysis relies heavily on the correct localization of the ligand binding site. We have developed SURFNET-ConSurf, a modular, two-stage method for identifying the location and shape of potential ligand binding pockets in protein structures. In the first stage, the SURFNET program identifies clefts in the protein surface that are potential binding sites. In the second stage, these clefts are trimmed in size by cutting away regions distant from highly conserved residues, as defined by the ConSurf-HSSP database. The largest clefts that remain tend to be those where ligands bind. To test the approach, we analyzed a nonredundant set of 244 protein structures from the PDB and found that SURFNET-ConSurf identifies a ligand binding pocket in 75% of them. The trimming procedure reduces the original cleft volumes by 30% on average, while still encompassing an average 87% of the ligand volume. From the analysis of the results we conclude that for those cases in which the ligands are found in large, highly conserved clefts, the combined SURFNET-ConSurf method gives pockets that are a better match to the ligand shape and location. We also show that this approach works better for enzymes than for nonenzyme proteins.

Hot regions in protein--protein interactions: the organization and contribution of structurally conserved hot spot residues.

O Keskin — 2005-03-24T22:35:25-00:00

J Mol Biol, Vol. 345, No. 5. (4 February 2005), pp. 1281-1294.

Structurally conserved residues at protein-protein interfaces correlate with the experimental alanine-scanning hot spots. Here, we investigate the organization of these conserved, computational hot spots and their contribution to the stability of protein associations. We find that computational hot spots are not homogeneously distributed along the protein interfaces; rather they are clustered within locally tightly packed regions. Within the dense clusters, they form a network of interactions and consequently their contributions to the stability of the complex are cooperative; however the contributions of independent clusters are additive. This suggests that the binding free energy is not a simple summation of the single hot spot residue contributions. As expected, around the hot spot residues we observe moderately conserved residues, further highlighting the crucial role of the conserved interactions in the local densely packed environment. The conserved occurrence of these organizations suggests that they are advantageous for protein-protein associations. Interestingly, the total number of hydrogen bonds and salt bridges contributed by hot spots is as expected. Thus, H-bond forming residues may use a "hot spot for water exclusion" mechanism. Since conserved residues are located within highly packed regions, water molecules are easily removed upon binding, strengthening electrostatic contributions of charge-charge interactions. Hence, the picture that emerges is that protein-protein associations are optimized locally, with the clustered, networked, highly packed structurally conserved residues contributing dominantly and cooperatively to the stability of the complex. When addressing the crucial question of "what are the preferred ways of proteins to associate", these findings point toward a critical involvement of hot regions in protein-protein interactions.

Evolutionary trace report maker: a new type of service for comparative analysis of proteins.

I Mihalek — 2006-04-30T14:10:03-00:00

Bioinformatics (27 April 2006)

SUMMARY: Evolutionary trace report_maker offers a new type of service for researchers investigating the function of novel proteins. It pools, from different sources, information about protein sequence, structure and elementary annotation, and to that background superimposes inference about the evolutionary behavior of individual residues, using real-valued evolutionary trace method. As its only input it takes a Protein Data Bank identifier or UniProt accession number, and returns a human-readable document in PDF format, supplemented by the original data needed to reproduce the results quoted in the report. AVAILABILITY: Evolutionary trace reports are freely available for academic users at http://mammoth.bcm.tmc.edu/report_maker.

Ribosome motions modulate electrostatic properties.

J Trylska — 2005-10-12T08:46:40-00:00

Biopolymers, Vol. 74, No. 6. (15 August 2004), pp. 423-431.

The electrostatic properties of the 70S ribosome of Thermus thermophilus were studied qualitatively by solving the Poisson-Boltzmann (PB) equation in aqueous solution and with physiological ionic strength. The electrostatic potential was calculated for conformations of the ribosome derived by recent normal mode analysis (Tama, F., et al. Proc Natl Acad Sci USA 2003 100, 9319-9323) of the ratchet-like reorganization that occurs during translocation (Frank, J.; Agrawal, R. K. Nature 2000 406, 318-322). To solve the PB equation, effective parameters (charges and radii), applicable to a highly charged backbone model of the ribosome, were developed. Regions of positive potential were found at the binding site of the elongation factors G and Tu, as well as where the release factors bind. Large positive potential areas are especially pronounced around the L11 and L6 proteins. The region around the L1 protein is also positively charged, supporting the idea that L1 may interact with the E-site tRNA during its release from the ribosome after translocation. Functional rearrangement of the ribosome leads to electrostatic changes which may help the translocation of the tRNAs during the elongation stage.

Symmetry, Form, and Shape: Guiding Principles for Robustness in Macromolecular Machines.

Florence Tama — 2006-05-01T17:05:55-00:00

Annu Rev Biophys Biomol Struct (13 January 2006)

Computational studies of large macromolecular assemblages have come a long way during the past 10 years. With the explosion of computer power and parallel computing, timescales of molecular dynamics simulations have been extended far beyond the hundreds of picoseconds timescale. However, limitations remain for studies of large-scale conformational changes occurring on timescales beyond nanoseconds especially for large macromolecules. In this review, we describe recent methods based on normal mode analysis that have enabled us to study dynamics on the microsecond timescale for large macromolecules using different levels of coarse graining, from atomically detailed models to those employing only low-resolution structural information. Emerging from such studies is a control principle for robustness in Nature's machines. We discuss this idea in the context of large-scale functionally reorganization of the ribosome, virus particles, and the muscle protein myosin. Expected online publication date for the Annual Review of Biophysics and Biomolecular Structure Volume 35 is May 5, 2006. Please see http://www.annualreviews.org/catalog/pub_dates.asp for revised estimates.

Hidden complexity of free energy surfaces for peptide (protein) folding.

SV Krivov — 2006-02-16T20:20:53-00:00

Proc Natl Acad Sci U S A, Vol. 101, No. 41. (12 October 2004), pp. 14766-14770.

An understanding of the thermodynamics and kinetics of protein folding requires a knowledge of the free energy surface governing the motion of the polypeptide chain. Because of the many degrees of freedom involved, surfaces projected on only one or two progress variables are generally used in descriptions of the folding reaction. Such projections result in relatively smooth surfaces, but they could mask the complexity of the unprojected surface. Here we introduce an approach to determine the actual (unprojected) free energy surface and apply it to the second beta-hairpin of protein G, which has been used as a model system for protein folding. The surface is represented by a disconnectivity graph calculated from a long equilibrium folding-unfolding trajectory. The denatured state is found to have multiple low free energy basins. Nevertheless, the peptide shows exponential kinetics in folding to the native basin. Projected surfaces obtained from the present analysis have a simple form in agreement with other studies of the beta-hairpin. The hidden complexity found for the beta-hairpin surface suggests that the standard funnel picture of protein folding should be revisited.

Folding of Ubiquitin: A Simple Model Describes the Strange Kinetics.

Sergei Chekmarev — 2006-05-01T16:59:03-00:00

J Phys Chem B Condens Matter Mater Surf Interfaces Biophys, Vol. 110, No. 17. (4 May 2006), pp. 8865-8869.

The ubiquitin mutant UbG folding experiments of Sabelko et al. (Proc. Natl. Acad. Sci. U.S.A. 1999, 96, 6031-6036), in which "strange kinetics" were observed, are interpreted in terms of a simple kinetic model. A minimal set of states consisting of a semicompact globule, two off-pathway traps, and the native state are included; the fully unfolded state is not considered because folding to the semicompact globule is fast. Both the low- and the high-temperature experiments of Sabelko et al. are fitted by a system of kinetic equations determining the transitions between these states. It is possible that cold- and heat-denaturated states of UbG are the basis of the off-pathway traps. The fits of the kinetic model to the experimental results provides an estimate of the rate constants for the various reaction channels and show how their contributions vary with temperature. Introduction of an on-pathway intermediate instead of one of the off-pathway traps does not lead to agreement with the experiments.

Molecular dynamics and protein function.

M Karplus — 2006-05-01T16:57:33-00:00

Proc Natl Acad Sci U S A, Vol. 102, No. 19. (10 May 2005), pp. 6679-6685.

A fundamental appreciation for how biological macromolecules work requires knowledge of structure and dynamics. Molecular dynamics simulations provide powerful tools for the exploration of the conformational energy landscape accessible to these molecules, and the rapid increase in computational power coupled with improvements in methodology makes this an exciting time for the application of simulation to structural biology. In this Perspective we survey two areas, protein folding and enzymatic catalysis, in which simulations have contributed to a general understanding of mechanism. We also describe results for the F(1) ATPase molecular motor and the Src family of signaling proteins as examples of applications of simulations to specific biological systems.

De novo protein design: towards fully automated sequence selection.

BI Dahiyat — 2006-05-01T16:52:07-00:00

J Mol Biol, Vol. 273, No. 4. (7 November 1997), pp. 789-796.

Several groups have applied and experimentally tested systematic, quantitative methods to protein design with the goal of developing general design algorithms. We have sought to expand the range of computational protein design by developing quantitative design methods for residues of all parts of a protein: the buried core, the solvent exposed surface, and the boundary between core and surface. Our goal is an objective, quantitative design algorithm that is based on the physical properties that determine protein structure and stability and which is not limited to specific folds or motifs. We chose the betabetaalpha motif typified by the zinc finger DNA binding module to test our design methodology. Using previously published sequence scoring functions developed with a combined experimental and computational approach and the Dead-End Elimination theorem to search for the optimal sequence, we designed 20 out of 28 positions in the test motif. The resulting sequence has less than 40% homology to any known sequence and does not contain any metal binding sites or cysteine residues. The resulting peptide, pda8d, is highly soluble and monomeric and circular dichroism measurements showed it to be folded with a weakly cooperative thermal unfolding transition. The NMR solution structure of pda8d was solved and shows that it is well-defined with a backbone ensemble rms deviation of 0. 55 A. Pda8d folds into the desired betabetaalpha motif with well-defined elements of secondary structure and tertiary organization. Superposition of the pda8d backbone to the design target is excellent, with an atomic rms deviation of 1.04 A.

Engineered antibody Fc variants with enhanced effector function.

GA Lazar — 2006-05-01T16:44:11-00:00

Proc Natl Acad Sci U S A, Vol. 103, No. 11. (14 March 2006), pp. 4005-4010.

Antibody-dependent cell-mediated cytotoxicity, a key effector function for the clinical efficacy of monoclonal antibodies, is mediated primarily through a set of closely related Fcgamma receptors with both activating and inhibitory activities. By using computational design algorithms and high-throughput screening, we have engineered a series of Fc variants with optimized Fcgamma receptor affinity and specificity. The designed variants display >2 orders of magnitude enhancement of in vitro effector function, enable efficacy against cells expressing low levels of target antigen, and result in increased cytotoxicity in an in vivo preclinical model. Our engineered Fc regions offer a means for improving the next generation of therapeutic antibodies and have the potential to broaden the diversity of antigens that can be targeted for antibody-based tumor therapy.

The burial of solvent-accessible surface area is a predictor of polypeptide folding and misfolding as a function of chain elongation.

N Kurt — 2006-05-01T11:39:24-00:00

J Am Chem Soc, Vol. 127, No. 45. (16 November 2005), pp. 15690-15691.

In quest of an empirical potential for protein structure prediction.

Jeffrey Skolnick — 2006-03-12T14:59:59-00:00

Curr Opin Struct Biol (6 March 2006)

Key to successful protein structure prediction is a potential that recognizes the native state from misfolded structures. Recent advances in empirical potentials based on known protein structures include improved reference states for assessing random interactions, sidechain-orientation-dependent pair potentials, potentials for describing secondary or supersecondary structural preferences and, most importantly, optimization protocols that sculpt the energy landscape to enhance the correlation between native-like features and the energy. Improved clustering algorithms that select native-like structures on the basis of cluster density also resulted in greater prediction accuracy. For template-based modeling, these advances allowed improvement in predicted structures relative to their initial template alignments over a wide range of target-template homology. This represents significant progress and suggests applications to proteome-scale structure prediction.

Ensemble versus single-molecule protein unfolding.

R Day — 2006-05-01T11:34:42-00:00

Proc Natl Acad Sci U S A, Vol. 102, No. 38. (20 September 2005), pp. 13445-13450.

Molecular dynamics (MD) simulations are the classic single-molecule "experiments," providing atomic-resolution structural and dynamic information. However, the single-molecule nature of the technique has also been its shortcoming, with frequent criticisms of sampling inadequacies and questions regarding the ensemble behavior of large numbers of molecules. Given the increase in computer power, we now address this issue by performing a large number of simulations and comparing individual and ensemble properties. One hundred independent MD simulations of the protein chymotrypsin inhibitor 2 were carried out for 20 ns each at 498 K in water to more fully describe the potentially diverse routes of protein unfolding and investigate how representative a single trajectory can be. Rapid unfolding was observed in all cases with the trajectories distributed about an average "ensemble" path in which secondary and tertiary structure was lost concomitantly, with tertiary structure loss occurring slightly faster. Individual trajectories did, however, sample conformations far from the average path with very heterogeneous time-dependent properties. Nevertheless, all of the simulations but one followed the average ensemble pathway, such that a small number of simulations (5-10) are sufficient to capture the average properties of these states and the unfolding pathway.

The present view of the mechanism of protein folding.

V Daggett — 2005-09-27T20:49:46-00:00

Nat Rev Mol Cell Biol, Vol. 4, No. 6. (June 2003), pp. 497-502.

We can track the positions and movements of all the atoms in small proteins as they fold and unfold by combining experimental studies with atomic-resolution molecular dynamics simulations. General principles as to how such complex architectures form so rapidly are now emerging from in-depth studies of a few proteins.

Kinetic Definition of Protein Folding Transition State Ensembles and Reaction Coordinates.

Christopher D Snow — 2006-05-01T11:29:32-00:00

Biophys J (14 April 2006)

Using distributed molecular dynamics simulations we located 4 distinct folding transitions for a 39 residue beta-beta-alpha-beta protein fold. To characterize the nature of each room temperature transition, we calculated the probability of transmission for 500 points along each free energy barrier. We introduced a method for determining transition states by employing the transmission probability, Ptrans, and determined which conformations were transition state ensemble members (Ptrans approximately 0.5). The transmission probability may be used to characterize the barrier in several ways. For example, we ran simulations at 82 degrees C, determined the change in Ptrans with temperature for all 2,000 conformations, and quantified Hammond behavior directly using Ptrans correlation. Additionally, we propose that diffusion along Ptrans may provide the configurational diffusion rate at the top of the barrier. Specifically, given a transition state conformation x0 with estimated Ptrans = 0.5, we selected a large set of subsequent conformations from independent trajectories, each exactly a small time delta after x0 (250ps). Calculating Ptrans for the new trial conformations we generated the P(Ptrans|delta=250ps) distribution that reflected diffusion. This approach provides a novel perspective on the diffusive nature of a protein folding transition, and provides a framework for quantitative study of activated relaxation kinetics.

Folding-based molecular simulations reveal mechanisms of the rotary motor F1-ATPase.

N Koga — 2006-05-01T11:26:30-00:00

Proc Natl Acad Sci U S A, Vol. 103, No. 14. (4 April 2006), pp. 5367-5372.

Biomolecular machines fulfill their function through large conformational changes that typically occur on the millisecond time scale or longer. Conventional atomistic simulations can only reach microseconds at the moment. Here, extending the minimalist model developed for protein folding, we propose the "switching Gō model" and use it to simulate the rotary motion of ATP-driven molecular motor F(1)-ATPase. The simulation recovers the unidirectional 120 degrees rotation of the gamma-subunit, the rotor. The rotation was induced solely by steric repulsion from the alpha(3)beta(3) subunits, the stator, which undergoes conformation changes during ATP hydrolysis. In silico alanine mutagenesis further elucidated which residues play specific roles in the rotation. Finally, regarding the mechanochemical coupling scheme, we found that the tri-site model does not lead to successful rotation but that the always-bi-site model produces approximately 30 degrees and approximately 90 degrees substeps, perfectly in accord with experiments. In the always-bi-site model, the number of sites occupied by nucleotides is always two during the hydrolysis cycle. This study opens up an avenue of simulating functional dynamics of huge biomolecules that occur on the millisecond time scales involving large-amplitude conformational change.

Molecular mechanism for stabilizing a short helical peptide studied by generalized-ensemble simulations with explicit solvent.

Y Sugita — 2006-05-01T11:24:04-00:00

Biophys J, Vol. 88, No. 5. (May 2005), pp. 3180-3190.

We study the folding mechanism of an analog of the C-peptide of ribonuclease A in explicit water by a replica-exchange multicanonical molecular dynamics simulation based on all-atom models. The multicanonical weight factor was determined by the combined use of the multicanonical replica-exchange method and the replica-exchange multicanonical algorithm. Using statistically reliable data thus obtained, we have examined the free-energy landscape of the peptide system. The global-minimum free-energy state in the landscape at room temperature has an alpha-helix structure with a distortion near the N-terminus. The state also has a salt bridge between Glu(-)-2 and Arg(+)-10 and an aromatic-aromatic interaction between Phe-8 and His(+)-12, both of which have been observed in x-ray and other experimental measurements. Principal component analysis clearly shows the different roles of these side-chain interactions in the peptide folding. The side-chain interaction between Phe-8 and His(+)-12 greatly enhances the stability of helical structure toward the C-terminal end, whereas the salt bridge between Glu(-)-2 and Arg(+)-10 mainly works as a restraint to prevent the alpha-helix structure from extending to the N-terminus. The free-energy landscape of C-peptide reveals a funnel-like shape where all of these interactions consistently exist only in the global-minimum state. This is the major reason why the native structure of the short helical peptide shows significant stability at low temperatures.

Protonation of the acidic residues in the transmembrane cation-binding sites of the ca(2+) pump.

Y Sugita — 2006-05-01T11:22:33-00:00

J Am Chem Soc, Vol. 127, No. 17. (4 May 2005), pp. 6150-6151.

Atomically detailed description of the unfolding of alpha-lactalbumin by the combined use of experiments and simulations.

T Oroguchi — 2006-05-01T11:20:33-00:00

J Mol Biol, Vol. 354, No. 1. (18 November 2005), pp. 164-172.

The recombinant form of goat alpha-lactalbumin has a significantly faster unfolding rate compared to the authentic form, although the two molecules differ only in an extra methionine at the N terminus of the recombinant. The mechanism of the destabilization caused by this residue was investigated through the combined use of kinetic experiments and molecular dynamics simulations. Unfolding simulations for the authentic and recombinant forms at 398 K (ten trajectories of 5 ns for each form, 100 ns total) precisely reproduced the experimentally observed differences in unfolding behavior. In addition, experiments reproduced the destabilization of a mutant protein, T38A, faithfully as predicted by the simulations. This bidirectional verification between experiments and simulations enabled the atomically detailed description of the role of the extra methionine residue in the unfolding process.

Native protein sequences are designed to destabilize folding intermediates.

Y Isogai — 2006-03-10T18:22:56-00:00

Biochemistry, Vol. 45, No. 8. (28 February 2006), pp. 2488-2492.

Hydrophobic core mutants of sperm whale apomyoglobin were constructed to investigate the amino acid sequence features that determine the folding properties. Replacements of all of the Ile residues with Leu and of all of the Ile and Val residues with Leu decreased the thermodynamic stability of the folded states against the unfolded states but increased the stability of the folding intermediates against the unfolded states, indicating that the amino acid composition of the protein core is important for the protein stability and folding cooperativity. To examine the effect of the arrangement of these hydrophobic residues, mutant proteins were further constructed: 12 sites out of the 18 Leu, 9 Ile, and 8 Val residues of the wild-type myoglobin were randomly replaced with each other so that the amino acid compositions were similar to that of the wild-type protein. Four mutant proteins were obtained without selection of the protein properties. These residue replacements similarly resulted in the stabilization of both the intermediate and folded states against the unfolded states, as compared to the wild-type protein. Thus, the arrangements of the hydrophobic residues in the native amino acid sequence are selected to destabilize the folding intermediate rather than to stabilize the folded state. The present results suggest that the two-state transition of protein folding or the transient formation of the unstable intermediate, which seems to be required for effective production of the functional proteins, has been a major driving force in the molecular evolution of natural globular proteins.

GASH: an improved algorithm for maximizing the number of equivalent residues between two protein structures.

DM Standley — 2006-03-31T18:44:17-00:00

BMC Bioinformatics, Vol. 6 (2005)

BACKGROUND: We introduce GASH, a new, publicly accessible program for structural alignment and superposition. Alignments are scored by the Number of Equivalent Residues (NER), a quantitative measure of structural similarity that can be applied to any structural alignment method. Multiple alignments are optimized by conjugate gradient maximization of the NER score within the genetic algorithm framework. Initial alignments are generated by the program Local ASH, and can be supplemented by alignments from any other program. RESULTS: We compare GASH to DaliLite, CE, and to our earlier program Global ASH on a difficult test set consisting of 3,102 structure pairs, as well as a smaller set derived from the Fischer-Eisenberg set. The extent of alignment crossover, as well as the completeness of the initial set of alignments are examined. The quality of the superpositions is evaluated both by NER and by the number of aligned residues under three different RMSD cutoffs (2,4, and 6A). In addition to the numerical assessment, the alignments for several biologically related structural pairs are discussed in detail. CONCLUSION: Regardless of which criteria is used to judge the superposition accuracy, GASH achieves the best overall performance, followed by DaliLite, Global ASH, and CE. In terms of CPU usage, DaliLite CE and GASH perform similarly for query proteins under 500 residues, but for larger proteins DaliLite is faster than GASH or CE. Both an http interface and a simple object application protocol (SOAP) interface to the GASH program are available at http://www.pdbj.org/GASH/.

The inference of protein-protein interactions by co-evolutionary analysis is improved by excluding the information about the phylogenetic relationships.

T Sato — 2005-11-09T11:01:31-00:00

Bioinformatics, Vol. 21, No. 17. (1 September 2005), pp. 3482-3489.

MOTIVATION: The prediction of protein-protein interactions is currently an important issue in bioinformatics. The mirror tree method uses evolutionary information to predict protein-protein interactions. However, it has been recognized that predictions by the mirror tree method lead to many false positives. The incentive of our study was to solve this problem by improving the method of extracting the co-evolutionary information regarding the protein pairs. RESULTS: We developed a novel method to predict protein-protein interactions from co-evolutionary information in the framework of the mirror tree method. The originality is the use of the projection operator to exclude the information about the phylogenetic relationships among the source organisms from the distance matrix. Each distance matrix was transformed into a vector for the operation. The vector is referred to as a 'phylogenetic vector'. We have proposed three ways to extract the phylogenetic information: (1) using the 16S rRNA from the same source organisms as the proteins under consideration, (2) averaging the phylogenetic vectors and (3) analyzing the principal components of the phylogenetic vectors. We examined the performance of the proposed methods to predict interacting protein pairs from Escherichia coli, using experimentally verified data. Our method was successful, and it drastically reduced the number of false positives in the prediction. AVAILABILITY: The R script for the prediction of protein-protein interactions reported in this manuscript is available at http://timpani.genome.ad.jp/~proj/ CONTACT: sato@kuicr.kyoto-u.ac.jp SUPPLEMENTARY INFORMATION: The information is also available at the same site as the R script.

The 3D profile method for identifying fibril-forming segments of proteins.

MJ Thompson — 2006-05-01T11:13:47-00:00

Proc Natl Acad Sci U S A, Vol. 103, No. 11. (14 March 2006), pp. 4074-4078.

Based on the crystal structure of the cross-beta spine formed by the peptide NNQQNY, we have developed a computational approach for identifying those segments of amyloidogenic proteins that themselves can form amyloid-like fibrils. The approach builds on experiments showing that hexapeptides are sufficient for forming amyloid-like fibrils. Each six-residue peptide of a protein of interest is mapped onto an ensemble of templates, or 3D profile, generated from the crystal structure of the peptide NNQQNY by small displacements of one of the two intermeshed beta-sheets relative to the other. The energy of each mapping of a sequence to the profile is evaluated by using ROSETTADESIGN, and the lowest energy match for a given peptide to the template library is taken as the putative prediction. If the energy of the putative prediction is lower than a threshold value, a prediction of fibril formation is made. This method can reach an accuracy of approximately 80% with a P value of approximately 10(-12) when a conservative energy threshold is used to separate peptides that form fibrils from those that do not. We see enrichment for positive predictions in a set of fibril-forming segments of amyloid proteins, and we illustrate the method with applications to proteins of interest in amyloid research.

Atomic-level description of amyloid beta-dimer formation.

S Gnanakaran — 2006-05-01T11:10:33-00:00

J Am Chem Soc, Vol. 128, No. 7. (22 February 2006), pp. 2158-2159.

Simulating movement of tRNA into the ribosome during decoding.

KY Sanbonmatsu — 2005-12-09T15:24:39-00:00

Proc Natl Acad Sci U S A, Vol. 102, No. 44. (1 November 2005), pp. 15854-15859.

Decoding is the key step during protein synthesis that enables information transfer from RNA to protein, a process critical for the survival of all organisms. We have used large-scale (2.64 x 10(6) atoms) all-atom simulations of the entire ribosome to understand a critical step of decoding. Although the decoding problem has been studied for more than four decades, the rate-limiting step of cognate tRNA selection has only recently been identified. This step, known as accommodation, involves the movement inside the ribosome of the aminoacyl-tRNA from the partially bound "A/T" state to the fully bound "A/A" state. Here, we show that a corridor of 20 universally conserved ribosomal RNA bases interacts with the tRNA during the accommodation movement. Surprisingly, the tRNA is impeded by the A-loop (23S helix 92), instead of enjoying a smooth transition to the A/A state. In particular, universally conserved 23S ribosomal RNA bases U2492, C2556, and C2573 act as a 3D gate, causing the acceptor stem to pause before allowing entrance into the peptidyl transferase center. Our simulations demonstrate that the flexibility of the acceptor stem of the tRNA, in addition to flexibility of the anticodon arm, is essential for tRNA selection. This study serves as a template for simulating conformational changes in large (>10(6) atoms) biological and artificial molecular machines.

Physically realistic homology models built with ROSETTA can be more accurate than their templates.

KM Misura — 2006-05-01T11:08:46-00:00

Proc Natl Acad Sci U S A, Vol. 103, No. 14. (4 April 2006), pp. 5361-5366.

We have developed a method that combines the rosetta de novo protein folding and refinement protocol with distance constraints derived from homologous structures to build homology models that are frequently more accurate than their templates. We test this method by building complete-chain models for a benchmark set of 22 proteins, each with 1 or 2 candidate templates, for a total of 39 test cases. We use structure-based and sequence-based alignments for each of the test cases. All atoms, including hydrogens, are represented explicitly. The resulting models contain approximately the same number of atomic overlaps as experimentally determined crystal structures and maintain good stereochemistry. The most accurate models can be identified by their energies, and in 22 of 39 cases a model that is more accurate than the template over aligned regions is one of the 10 lowest-energy models.

Generation of native-like protein structures from limited NMR data, modern force fields and advanced conformational sampling

Jianhan Chen — 2005-02-09T23:48:42-00:00

Journal of Biomolecular NMR, Vol. 31, No. 1. (January 2005), pp. 59-64.

Refinement of NMR structures using implicit solvent and advanced sampling techniques.

J Chen — 2005-08-04T12:03:21-00:00

J Am Chem Soc, Vol. 126, No. 49. (15 December 2004), pp. 16038-16047.

NMR biomolecular structure calculations exploit simulated annealing methods for conformational sampling and require a relatively high level of redundancy in the experimental restraints to determine quality three-dimensional structures. Recent advances in generalized Born (GB) implicit solvent models should make it possible to combine information from both experimental measurements and accurate empirical force fields to improve the quality of NMR-derived structures. In this paper, we study the influence of implicit solvent on the refinement of protein NMR structures and identify an optimal protocol of utilizing these improved force fields. To do so, we carry out structure refinement experiments for model proteins with published NMR structures using full NMR restraints and subsets of them. We also investigate the application of advanced sampling techniques to NMR structure refinement. Similar to the observations of Xia et al. (J.Biomol. NMR 2002, 22, 317-331), we find that the impact of implicit solvent is rather small when there is a sufficient number of experimental restraints (such as in the final stage of NMR structure determination), whether implicit solvent is used throughout the calculation or only in the final refinement step. The application of advanced sampling techniques also seems to have minimal impact in this case. However, when the experimental data are limited, we demonstrate that refinement with implicit solvent can substantially improve the quality of the structures. In particular, when combined with an advanced sampling technique, the replica exchange (REX) method, near-native structures can be rapidly moved toward the native basin. The REX method provides both enhanced sampling and automatic selection of the most native-like (lowest energy) structures. An optimal protocol based on our studies first generates an ensemble of initial structures that maximally satisfy the available experimental data with conventional NMR software using a simplified force field and then refines these structures with implicit solvent using the REX method. We systematically examine the reliability and efficacy of this protocol using four proteins of various sizes ranging from the 56-residue B1 domain of Streptococcal protein G to the 370-residue Maltose-binding protein. Significant improvement in the structures was observed in all cases when refinement was based on low-redundancy restraint data. The proposed protocol is anticipated to be particularly useful in early stages of NMR structure determination where a reliable estimate of the native fold from limited data can significantly expedite the overall process. This refinement procedure is also expected to be useful when redundant experimental data are not readily available, such as for large multidomain biomolecules and in solid-state NMR structure determination.

GROMACS: fast, flexible, and free.

D Van Der Spoel — 2005-12-15T21:09:13-00:00

J Comput Chem, Vol. 26, No. 16. (December 2005), pp. 1701-1718.

This article describes the software suite GROMACS (Groningen MAchine for Chemical Simulation) that was developed at the University of Groningen, The Netherlands, in the early 1990s. The software, written in ANSI C, originates from a parallel hardware project, and is well suited for parallelization on processor clusters. By careful optimization of neighbor searching and of inner loop performance, GROMACS is a very fast program for molecular dynamics simulation. It does not have a force field of its own, but is compatible with GROMOS, OPLS, AMBER, and ENCAD force fields. In addition, it can handle polarizable shell models and flexible constraints. The program is versatile, as force routines can be added by the user, tabulated functions can be specified, and analyses can be easily customized. Nonequilibrium dynamics and free energy determinations are incorporated. Interfaces with popular quantum-chemical packages (MOPAC, GAMES-UK, GAUSSIAN) are provided to perform mixed MM/QM simulations. The package includes about 100 utility and analysis programs. GROMACS is in the public domain and distributed (with source code and documentation) under the GNU General Public License. It is maintained by a group of developers from the Universities of Groningen, Uppsala, and Stockholm, and the Max Planck Institute for Polymer Research in Mainz. Its Web site is http://www.gromacs.org.

Molecular dynamics studies of the archaeal translocon.

James C Gumbart — 2006-02-27T12:38:01-00:00

Biophys J (13 January 2006)

The translocon is a protein-conducting channel conserved over all domains of life that serves to translocate proteins across or into membranes. While this channel has been well-studied for many years, the recent discovery of a high-resolution crystal structure opens up new avenues of exploration. Taking advantage of this, we performed molecular dynamics simulations of the translocon in a fully-solvated lipid bilayer, examining the translocation abilities of monomeric SecYEbeta by forcing two helices comprised of different amino acid sequences to cross the channel. The simulations revealed that the so-called ;plug' of SecYEbeta swings open during translocation, closing thereafter. Likewise, it was established that the so-called ;pore ring' region of SecYEbeta forms an elastic, yet tight seal around the translocating oligopeptides. The closed state of the channel was found to block permeation of all ions and water molecules; in the open state, ions were blocked. Our results suggest that the SecYEbeta monomer is capable of forming an active channel.

Scalable molecular dynamics with NAMD.

JC Phillips — 2005-12-25T22:19:06-00:00

J Comput Chem, Vol. 26, No. 16. (December 2005), pp. 1781-1802.

NAMD is a parallel molecular dynamics code designed for high-performance simulation of large biomolecular systems. NAMD scales to hundreds of processors on high-end parallel platforms, as well as tens of processors on low-cost commodity clusters, and also runs on individual desktop and laptop computers. NAMD works with AMBER and CHARMM potential functions, parameters, and file formats. This article, directed to novices as well as experts, first introduces concepts and methods used in the NAMD program, describing the classical molecular dynamics force field, equations of motion, and integration methods along with the efficient electrostatics evaluation algorithms employed and temperature and pressure controls used. Features for steering the simulation across barriers and for calculating both alchemical and conformational free energy differences are presented. The motivations for and a roadmap to the internal design of NAMD, implemented in C++ and based on Charm++ parallel objects, are outlined. The factors affecting the serial and parallel performance of a simulation are discussed. Finally, typical NAMD use is illustrated with representative applications to a small, a medium, and a large biomolecular system, highlighting particular features of NAMD, for example, the Tcl scripting language. The article also provides a list of the key features of NAMD and discusses the benefits of combining NAMD with the molecular graphics/sequence analysis software VMD and the grid computing/collaboratory software BioCoRE. NAMD is distributed free of charge with source code at www.ks.uiuc.edu.

Folding Cooperativity in a Three-stranded beta-Sheet Model.

DR Roe — 2005-09-13T13:12:34-00:00

J Mol Biol, Vol. 352, No. 2. (16 September 2005), pp. 370-381.

The thermodynamic behavior of a previously designed three-stranded beta-sheet was studied via several microseconds of standard and replica exchange molecular dynamics simulations. The system is shown to populate at least four thermodynamic minima, including two partially folded states in which only a single hairpin is formed. Simulated melting curves show different profiles for the C and N-terminal hairpins, consistent with differences in secondary structure content in published NMR and CD/FTIR measurements, which probed different regions of the chain. Individual beta-hairpins that comprise the three-stranded beta-sheet are observed to form cooperatively. Partial folding cooperativity between the component hairpins is observed, and good agreement between calculated and experimental values quantifying this cooperativity is obtained when similar analysis techniques are used. However, the structural detail in the ensemble of conformations sampled in the simulations permits a more direct analysis of this cooperativity than has been performed on the basis of experimental data. The results indicate the actual folding cooperativity perpendicular to strand direction is significantly larger than the lower bound obtained previously.

A multistep approach to structure-based drug design: studying ligand binding at the human neutrophil elastase.

T Steinbrecher — 2006-05-01T10:58:46-00:00

J Med Chem, Vol. 49, No. 6. (23 March 2006), pp. 1837-1844.

In this study we show that a combination of different theoretical methods is a viable approach to calculate the binding affinities of new ligands for the human neutrophile elastase. This protease degrades elastin and likely aids neutrophils in fulfilling their immunological functions. Abnormally high human neutrophil elastase (HNE) levels are involved in several diseases; therefore, inhibitors of HNE are of interest as targets for drug design. A recent study has revealed that cinnamic acid and bornyl ester derivatives bind to HNE, but DeltaG0 values from ligand docking results exhibited no correlation with those calculated from the IC50 values. To accurately compute binding affinities, we generated possible protein ligand complex structures by ligand docking calculations. For each of the ligands, the 30 most likely placements were used as starting points of nanosecond length molecular dynamics simulations. The binding free energies for these complex structures were estimated using a continuum solvent (MM-PBSA) approach. These results, along with structural data from the molecular dynamics runs, allowed the identification of a group of similar placements that serve as a model for the natural protein ligand complex structure. This structural model was used to perform thermodynamic integration (TI) calculations to obtain the relative binding free energies of similar ligands to HNE. The TI results were in quantitative agreement with the measured binding affinities. Thus, the presented approach can be used to generate a probable complex structure for known ligands to HNE and to use such a structure to calculate the effects of small ligand modifications on ligand binding, possibly leading to new inhibitors with improved binding affinities.

All-atom structure prediction and folding simulations of a stable protein.

C Simmerling — 2005-09-13T13:24:17-00:00

J Am Chem Soc, Vol. 124, No. 38. (25 September 2002), pp. 11258-11259.

The Amber biomolecular simulation programs.

DA Case — 2005-10-19T23:51:42-00:00

J Comput Chem, Vol. 26, No. 16. (December 2005), pp. 1668-1688.

We describe the development, current features, and some directions for future development of the Amber package of computer programs. This package evolved from a program that was constructed in the late 1970s to do Assisted Model Building with Energy Refinement, and now contains a group of programs embodying a number of powerful tools of modern computational chemistry, focused on molecular dynamics and free energy calculations of proteins, nucleic acids, and carbohydrates. (c) 2005 Wiley Periodicals, Inc. J Comput Chem 26: 1668-1688, 2005.

Effective function annotation through catalytic residue conservation.

Richard A George — 2005-08-05T22:21:55-00:00

Proc Natl Acad Sci U S A (21 July 2005)

This contribution is part of the special series of Inaugural Articles by members of the National Academy of Sciences elected on April 29, 2003. Because of the extreme impact of genome sequencing projects, protein sequences without accompanying experimental data now dominate public databases. Homology searches, by providing an opportunity to transfer functional information between related proteins, have become the de facto way to address this. Although a single, well annotated, close relationship will often facilitate sufficient annotation, this situation is not always the case, particularly if mutations are present in important functional residues. When only distant relationships are available, the transfer of function information is more tenuous, and the likelihood of encountering several well annotated proteins with different functions is increased. The consequence for a researcher is a range of candidate functions with little way of knowing which, if any, are correct. Here, we address the problem directly by introducing a computational approach to accurately identify and segregate related proteins into those with a functional similarity and those where function differs. This approach should find a wide range of applications, including the interpretation of genomics/proteomics data and the prioritization of targets for high-throughput structure determination. The method is generic, but here we concentrate on enzymes and apply high-quality catalytic site data. In addition to providing a series of comprehensive benchmarks to show the overall performance of our approach, we illustrate its utility with specific examples that include the correct identification of haptoglobin as a nonenzymatic relative of trypsin, discrimination of acid-D-amino acid ligases from a much larger ligase pool, and the successful annotation of BioH, a structural genomics target.

An Analysis of Intron Positions in Relation to Nucleotides, Amino Acids, and Protein Secondary Structure.

Gordon S Whamond — 2006-05-01T10:52:51-00:00

J Mol Biol (29 March 2006)

We present an analysis of intron positions in relation to nucleotides, amino acid residues, and protein secondary structure. Previous work has shown that intron sites in proteins are not randomly distributed with respect to secondary structures. Here we show that this preference can be almost totally explained by the nucleotide bias of splice site machinery, and may well not relate to protein stability or conformation at all. Each intron phase is preferentially associated with its own set of residues: phase 0 introns with lysine, glutamine, and glutamic acid before the intron, and valine after; phase 1 introns with glycine, alanine, valine, aspartic acid, and glutamic acid; and phase 2 introns with arginine, serine, lysine, and tryptophan. These preferences can be explained principally on the basis of nucleotide bias at intron locations, which is in accordance with previous literature. Although this work does not prove that introns are inserted into genomes at specific proto-splice sites, it shows that the nucleotide bias surrounding introns, however it originally occurred, explains the observed correlations between introns and protein secondary structure.

Ab initio molecular dynamics with discrete variable representation basis sets: techniques and application to liquid water.

HS Lee — 2006-05-01T10:51:24-00:00

J Phys Chem A Mol Spectrosc Kinet Environ Gen Theory, Vol. 110, No. 16. (27 April 2006), pp. 5549-5560.

Finite temperature ab initio molecular dynamics (AIMD), in which forces are obtained from "on-the-fly" electronic structure calculations, is a widely used technique for studying structural and dynamical properties of chemically active systems. Recently, we introduced an AIMD scheme based on discrete variable representation (DVR) basis sets, which was shown to have improved convergence properties over the conventional plane wave (PW) basis set [Liu,Y.; et al. Phys. Rev. B 2003, 68, 125110]. In the present work, the numerical algorithms for the DVR based AIMD scheme (DVR/AIMD) are provided in detail, and the latest developments of the approach are presented. The accuracy and stability of the current implementation of the DVR/AIMD scheme are tested by performing a simulation of liquid water at ambient conditions. The structural information obtained from the present work is in good agreement with the result of recent AIMD simulations with a PW basis set (PW/AIMD). Advantages of using the DVR/AIMD scheme over the PW/AIMD method are discussed. In particular, it is shown that a DVR/AIMD simulation of liquid water in the complete basis set limit is possible with a relatively small number of grid points.

Potentials of mean force for acetylcholine unbinding from the alpha7 nicotinic acetylcholine receptor ligand-binding domain [j. Am. Chem. Soc. 2006, 128, 3019-3026].

D Zhang — 2006-05-01T10:50:32-00:00

J Am Chem Soc, Vol. 128, No. 13. (5 April 2006)

Coupling hydrophobicity, dispersion, and electrostatics in continuum solvent models.

J Dzubiella — 2006-05-01T10:49:47-00:00

Phys Rev Lett, Vol. 96, No. 8. (3 March 2006)

An implicit solvent model is presented that couples hydrophobic, dispersion, and electrostatic solvation energies by minimizing the system Gibbs free energy with respect to the solvent volume exclusion function. The solvent accessible surface is the output of the theory. The method is illustrated with the solvation of simple solutes on different length scales and captures the sensitivity of hydration to the particular form of the solute-solvent interactions in agreement with recent computer simulations.

Electrostatic properties of cowpea chlorotic mottle virus and cucumber mosaic virus capsids.

R Konecny — 2006-05-01T10:48:54-00:00

Biopolymers, Vol. 82, No. 2. (June 2006), pp. 106-120.

Electrostatic properties of cowpea chlorotic mottle virus (CCMV) and cucumber mosaic virus (CMV) were investigated using numerical solutions to the Poisson-Boltzmann equation. Experimentally, it has been shown that CCMV particles swell in the absence of divalent cations when the pH is raised from 5 to 7. CMV, although structurally homologous, does not undergo this transition. An analysis of the calculated electrostatic potential confirms that a strong electrostatic repulsion at the calcium-binding sites in the CCMV capsid is most likely the driving force for the capsid swelling process during the release of calcium. The binding interaction between the encapsulated genome material (RNA) inside of the capsid and the inner capsid shell is weakened during the swelling transition. This probably aids in the RNA release process, but it is unlikely that the RNA is released through capsid openings due to unfavorable electrostatic interaction between the RNA and capsid inner shell residues at these openings. Calculations of the calcium binding energies show that Ca(2+) can bind both to the native and swollen forms of the CCMV virion. Favorable binding to the swollen form suggests that Ca(2+) ions can induce the capsid contraction and stabilize the native form. (c) 2005 Wiley Periodicals, Inc. Biopolymers 82: 106-120, 2006This article was originally published online as an accepted preprint. The "Published Online" date corresponds to the preprint version. You can request a copy of the preprint by emailing the Biopolymers editorial office at biopolymers@wiley.com.

Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field.

J Chen — 2006-05-01T10:48:09-00:00

J Am Chem Soc, Vol. 128, No. 11. (22 March 2006), pp. 3728-3736.

The efficient and accurate characterization of solvent effects is a key element in the theoretical and computational study of biological problems. Implicit solvent models, particularly generalized Born (GB) continuum electrostatics, have emerged as an attractive tool to study the structure and dynamics of biomolecules in various environments. Despite recent advances in this methodology, there remain limitations in the parametrization of many of these models. In the present work, we demonstrate that it is possible to achieve a balanced implicit solvent force field by further optimizing the input atomic radii in combination with adjusting the protein backbone torsional energetics. This parameter optimization is guided by the potentials of mean force (PMFs) between amino acid polar groups, calculated from explicit solvent free energy simulations, and by conformational equilibria of short peptides, obtained from extensive folding and unfolding replica exchange molecular dynamics (REX-MD) simulations. Through the application of this protocol, the delicate balance between the competing solvation forces and intramolecular forces appears to be better captured, and correct conformational equilibria for a range of both helical and beta-hairpin peptides are obtained. The same optimized force field also successfully folds both beta-hairpin trpzip2 and mini-protein Trp-Cage, indicating that it is quite robust. Such a balanced, physics-based force field will be highly applicable to a range of biological problems including protein folding and protein structural dynamics.