Evaluating Molecular Mechanical Potentials for Helical Peptides and Proteins

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

Non-Bulk-Like Solvent Behavior in the Ribosome Exit Tunnel

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities

Amyloid-Associated Nucleic Acid Hybridisation

Nucleic acids promote amyloid formation in diseases including Alzheimer's and Creutzfeldt-Jakob disease. However, it remains unclear whether the close interactions between amyloid and nucleic acid allow nucleic acid secondary structure to play a role in modulating amyloid structure and function. Here we have used a simplified system of short basic peptides with alternating hydrophobic and hydrophilic amino acid residues to study nucleic acid - amyloid interactions. Employing biophysical techniques including X-ray fibre diffraction, circular dichroism spectroscopy and electron microscopy we show that the polymerized charges of nucleic acids concentrate and enhance the formation of amyloid from short basic peptides, many of which would not otherwise form fibres. In turn, the amyloid component binds nucleic acids and promotes their hybridisation at concentrations below their solution Kd, as shown by time-resolved FRET studies. The self-reinforcing interactions between peptides and nucleic acids lead to the formation of amyloid nucleic acid (ANA) fibres whose properties are distinct from their component polymers. In addition to their importance in disease and potential in engineering, ANA fibres formed from prebiotically-produced peptides and nucleic acids may have played a role in early evolution, constituting the first entities subject to Darwinian evolution

SH3 Domain-Peptide Binding Energy Calculations Based on Structural Ensemble and Multiple Peptide Templates

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

Nucleotide Binding Switches the Information Flow in Ras GTPases

The Ras superfamily comprises many guanine nucleotide-binding proteins (G proteins) that are essential to intracellular signal transduction. The guanine nucleotide-dependent intrinsic flexibility patterns of five G proteins were investigated in atomic detail through Molecular Dynamics simulations of the GDP- and GTP-bound states (SGDP and SGTP, respectively). For all the considered systems, the intrinsic flexibility of SGDP was higher than that of SGTP, suggesting that Guanine Exchange Factor (GEF) recognition and nucleotide switch require higher amplitude motions than effector recognition or GTP hydrolysis. Functional mode, dynamic domain, and interaction energy correlation analyses highlighted significant differences in the dynamics of small G proteins and Gα proteins, especially in the inactive state. Indeed, SGDP of Gαt, is characterized by a more extensive energy coupling between nucleotide binding site and distal regions involved in GEF recognition compared to small G proteins, which attenuates in the active state. Moreover, mechanically distinct domains implicated in nucleotide switch could be detected in the presence of GDP but not in the presence of GTP. Finally, in small G proteins, functional modes are more detectable in the inactive state than in the active one and involve changes in solvent exposure of two highly conserved amino acids in switches I and II involved in GEF recognition. The average solvent exposure of these amino acids correlates in turn with the rate of GDP release, suggesting for them either direct or indirect roles in the process of nucleotide switch. Collectively, nucleotide binding changes the information flow through the conserved Ras-like domain, where GDP enhances the flexibility of mechanically distinct portions involved in nucleotide switch, and favors long distance allosteric communication (in Gα proteins), compared to GTP

Studying the Underlying Event in Drell-Yan and High Transverse Momentum Jet Production at the Tevatron

We study the underlying event in proton-antiproton collisions by examining the behavior of charged particles (transverse momentum pT > 0.5 GeV/c, pseudorapidity |\eta| < 1) produced in association with large transverse momentum jets (~2.2 fb-1) or with Drell-Yan lepton-pairs (~2.7 fb-1) in the Z-boson mass region (70 < M(pair) < 110 GeV/c2) as measured by CDF at 1.96 TeV center-of-mass energy. We use the direction of the lepton-pair (in Drell-Yan production) or the leading jet (in high-pT jet production) in each event to define three regions of \eta-\phi space; toward, away, and transverse, where \phi is the azimuthal scattering angle. For Drell-Yan production (excluding the leptons) both the toward and transverse regions are very sensitive to the underlying event. In high-pT jet production the transverse region is very sensitive to the underlying event and is separated into a MAX and MIN transverse region, which helps separate the hard component (initial and final-state radiation) from the beam-beam remnant and multiple parton interaction components of the scattering. The data are corrected to the particle level to remove detector effects and are then compared with several QCD Monte-Carlo models. The goal of this analysis is to provide data that can be used to test and improve the QCD Monte-Carlo models of the underlying event that are used to simulate hadron-hadron collisions.Comment: Submitted to Phys.Rev.

Precision measurement of the top quark mass from dilepton events at CDF II

We report a measurement of the top quark mass, M_t, in the dilepton decay channel of $t\bar{t}\to b\ell'^{+}\nu_{\ell'}\bar{b}\ell^{-}\bar{\nu}_{\ell}$ using an integrated luminosity of 1.0 fb^{-1} of p\bar{p} collisions collected with the CDF II detector. We apply a method that convolutes a leading-order matrix element with detector resolution functions to form event-by-event likelihoods; we have enhanced the leading-order description to describe the effects of initial-state radiation. The joint likelihood is the product of the likelihoods from 78 candidate events in this sample, which yields a measurement of M_{t} = 164.5 \pm 3.9(\textrm{stat.}) \pm 3.9(\textrm{syst.}) \mathrm{GeV}/c^2, the most precise measurement of M_t in the dilepton channel.Comment: 7 pages, 2 figures, version includes changes made prior to publication by journa

Measurement of the Ratios of Branching Fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) and B(Bs -> Ds pi) / B(Bd -> Dd pi)

Using 355 pb^-1 of data collected by the CDF II detector in \ppbar collisions at sqrt{s} = 1.96 TeV at the Fermilab Tevatron, we study the fully reconstructed hadronic decays B -> D pi and B -> D pi pi pi. We present the first measurement of the ratio of branching fractions B(Bs -> Ds pi pi pi) / B(Bd -> Dd pi pi pi) = 1.05 pm 0.10 (stat) pm 0.22 (syst). We also update our measurement of B(Bs -> Ds pi) / B(Bd -> Dd pi) to 1.13 pm 0.08 (stat) pm 0.23 (syst) improving the statistical uncertainty by more than a factor of two. We find B(Bs -> Ds pi) = [3.8 pm 0.3 (stat) pm 1.3 (syst)] \times 10^{-3} and B(Bs -> Ds pi pi pi) = [8.4 pm 0.8 (stat) pm 3.2 (syst)] \times 10^{-3}.Comment: 7 pages, 2 figure

Cross Section Measurements of High-pTp_T Dilepton Final-State Processes Using a Global Fitting Method

We present a new method for studying high-$p_T$ dilepton events ($e^{\pm}e^{\mp}$, $\mu^{\pm}\mu^{\mp}$, $e^{\pm}\mu^{\mp}$) and simultaneously extracting the production cross sections of $p\bar{p} \to t\bar{t}$, $p\bar{p} \to W^+W^-$, and p\bar{p} \to \ztt at a center-of-mass energy of $\sqrt{s} = 1.96$ TeV. We perform a likelihood fit to the dilepton data in a parameter space defined by the missing transverse energy and the number of jets in the event. Our results, which use $360 {\rm pb^{-1}}$ of data recorded with the CDF II detector at the Fermilab Tevatron Collider, are $\sigma(t\bar{t}) = 8.5_{-2.2}^{+2.7}$ pb, $\sigma(W^+W^-) = 16.3^{+5.2}_{-4.4}$ pb, and \sigma(\ztt) =291^{+50}_{-46} pb.Comment: 20 pages, 2 figures, to be submitted to PRD-R

Measurement of the Dipion Mass Spectrum in X(3872) -> J/Psi Pi+ Pi- Decays

We measure the dipion mass spectrum in X(3872)--> J/Psi Pi+ Pi- decays using 360 pb-1 of pbar-p collisions at 1.96 TeV collected with the CDF II detector. The spectrum is fit with predictions for odd C-parity (3S1, 1P1, and 3DJ) charmonia decaying to J/Psi Pi+ Pi-, as well as even C-parity states in which the pions are from Rho0 decay. The latter case also encompasses exotic interpretations, such as a D0-D*0Bar molecule. Only the 3S1 and J/Psi Rho hypotheses are compatible with our data. Since 3S1 is untenable on other grounds, decay via J/Psi Rho is favored, which implies C=+1 for the X(3872). Models for different J/Psi-Rho angular momenta L are considered. Flexibility in the models, especially the introduction of Rho-Omega interference, enable good descriptions of our data for both L=0 and 1.Comment: 7 pages, 4 figures -- Submitted to Phys. Rev. Let