Structure and Disorder in an Unfolded State under Nondenaturing Conditions from Ensemble Models Consistent with a Large Number of Experimental Restraints

Obtaining detailed structural models of disordered states of proteins under nondenaturing conditions is important for a better understanding of both functional intrinsically disordered proteins and unfolded states of folded proteins. Extensive experimental characterization of the drk N-terminal SH3 domain unfolded state has shown that, although it appears to be highly disordered, it possesses significant nonrandom secondary and tertiary structure. In our previous attempts to generate structural models of the unfolded state using the program ENSEMBLE, we were limited by insufficient experimental restraints and conformational sampling. In this study, we have vastly expanded our experimental restraint set to include ¹H-¹⁵N residual dipolar couplings, small-angle X-ray scattering measurements, nitroxide paramagnetic relaxation enhancements, O₂-induced ¹³C paramagnetic shifts, hydrogen-exchange protection factors, and ¹⁵N R₂ data, in addition to the previously used nuclear Overhauser effects, amino terminal Cu²⁺-Ni²⁺ binding paramagnetic relaxation enhancements, J-couplings, chemical shifts, hydrodynamic radius, and solvent accessibility restraints. We have also implemented a new ensemble calculation methodology that uses iterative conformational sampling and seeks to calculate the simplest possible ensemble models. As a result, we can now generate ensembles that are consistent with much larger experimental data sets than was previously possible. Although highly heterogeneous and having broad molecular size distributions, the calculated drk N-terminal SH3 domain unfolded-state ensembles have very different properties than expected for random or statistical coils and possess significant nonnative α-helical structure and both native-like and nonnative tertiary structure.     

Simulating structural and thermodynamic properties of carcinogen-damaged DNA

A pair of stereoisomeric covalent adducts to guanine in double-stranded DNA, derived from the reaction of mutagenic and tumorigenic metabolites of benzo[a]pyrene, have been well characterized structurally and thermodynamically. Both high-resolution NMR solution structures and an array of thermodynamic data are available for these 10S (+)- and 10R (-)-trans-anti -[BP]-N(2)-dG adducts in double-stranded deoxyoligonucleotides. The availability of experimentally well-characterized duplexes containing these two stereoisomeric guanine adducts provides an opportunity for evaluating the molecular mechanics Poisson-Boltzmann surface area (MM-PBSA) method for computing thermodynamic properties from molecular dynamics ensembles. We have carried out 3-ns molecular dynamics simulations, using NMR solution structures as the starting models for the 10S (+)- and 10R (-)-trans-anti-dG adducts in a DNA duplex 11-mer using AMBER 6.0. We employed the MM-PBSA method to compute the free energies, enthalpies, and entropies of the two adducts. Our complete thermodynamic analysis agrees quite well with the full experimental thermodynamic characterization of these adducts, showing essentially equal stabilities of the two adducts. We also calculated the nuclear Overhauser effect (NOE) distances from the molecular dynamics trajectories, and compared them against the experimental NMR-derived NOE distances. Our results showed that the simulated structures are in good agreement with the NMR experimental NOE data. Furthermore, the molecular dynamics simulations provided new structural and biological insights. Specifically, the puzzling observation that the BP aromatic ring system in the 10S (+)-trans-anti-dG adduct is more exposed to the aqueous solvent than the 10R (-)-trans-anti-dG adduct, is rationalized in terms of the adduct structures. The structural and thermodynamic features of these stereoisomeric adducts are also discussed in relation to their reported low susceptibilities to nucleotide excision repair.     

The structure of the Alzheimer amyloid beta 10-35 peptide probed through replica-exchange molecular dynamics simulations in explicit solvent

The conformational states sampled by the Alzheimer amyloid beta (10-35) (Abeta 10-35) peptide were probed using replica-exchange molecular dynamics (REMD) simulations in explicit solvent. The Abeta 10-35 peptide is a fragment of the full-length Abeta 40/42 peptide that possesses many of the amyloidogenic properties of its full-length counterpart. Under physiological temperature and pressure, our simulations reveal that the Abeta 10-35 peptide does not possess a single unique folded state. Rather, this peptide exists as a mixture of collapsed globular states that remain in rapid dynamic equilibrium with each other. This conformational ensemble is dominated by random coil and bend structures with insignificant presence of an alpha-helical or beta-sheet structure. The 3D structure of Abeta 10-35 is seen to be defined by a salt bridge formed between the side-chains of K28 and D23. This salt bridge is also observed in Abeta fibrils and our simulations suggest that monomeric conformations of Abeta 10-35 contain pre-folded structural motifs that promote rapid aggregation of this peptide.     

Characterization of the conformational and orientational dynamics of ganglioside GM1 in a dipalmitoylphosphatidylcholine bilayer by molecular dynamics simulations

The structure and dynamics of a single GM1 (Gal5-β1,3-GalNAc4-β1,4-(NeuAc3-α2,3)-Gal2-β1,4-Glc1-β1,1-Cer) embedded in a DPPC bilayer have been studied by MD simulations. Eleven simulations, each of 10 ns productive run, were performed with different initial conformations of GM1. Simulations of GM1-Os in water and of a DPPC bilayer were also performed to delineate the effects of the bilayer and GM1 on the conformational and orientational dynamics of each other. The conformation of the GM1 headgroup observed in the simulations is in agreement with those reported in literature; but the headgroup is restricted when embedded in the bilayer. NeuAc3 is the outermost saccharide towards the water phase. Glc1 and Gal2 prefer a parallel, and NeuAc3, GalNac4 and Gal5 prefer a perpendicular, orientation with respect to the bilayer normal. The overall characteristics of the bilayer are not affected by the presence of GM1; however, GM1 does influence the DPPC molecules in its immediate vicinity. The implications of these observations on the specific recognition and binding of GM1 embedded in a lipid bilayer by exogenous proteins as well as proteins embedded in lipids have been discussed.     

BPPred: a Web-based computational tool for predicting biophysical parameters of proteins

We exploit the availability of recent experimental data on a variety of proteins to develop a Web-based prediction algorithm (BPPred) to calculate several biophysical parameters commonly used to describe the folding process. These parameters include the equilibrium m-values, the length of proteins, and the changes upon unfolding in the solvent-accessible surface area, in the heat capacity, and in the radius of gyration. We also show that the knowledge of any one of these quantities allows an estimate of the others to be obtained, and describe the confidence limits with which these estimations can be made. Furthermore, we discuss how the kinetic m-values, or the Beta Tanford values, may provide an estimate of the solvent-accessible surface area and the radius of gyration of the transition state for protein folding. Taken together, these results suggest that BPPred should represent a valuable tool for interpreting experimental measurements, as well as the results of molecular dynamics simulations.     

Computational Design of a PAK1 Binding Protein

We describe a computational protocol, called DDMI, for redesigning scaffold proteins to bind to a specified region on a target protein. The DDMI protocol is implemented within the Rosetta molecular modeling program and uses rigid-body docking, sequence design, and gradient-based minimization of backbone and side-chain torsion angles to design low-energy interfaces between the scaffold and target protein. Iterative rounds of sequence design and conformational optimization were needed to produce models that have calculated binding energies that are similar to binding energies calculated for native complexes. We also show that additional conformation sampling with molecular dynamics can be iterated with sequence design to further lower the computed energy of the designed complexes. To experimentally test the DDMI protocol, we redesigned the human hyperplastic discs protein to bind to the kinase domain of p21-activated kinase 1 (PAK1). Six designs were experimentally characterized. Two of the designs aggregated and were not characterized further. Of the remaining four designs, three bound to the PAK1 with affinities tighter than 350 μM. The tightest binding design, named Spider Roll, bound with an affinity of 100 μM. NMR-based structure prediction of Spider Roll based on backbone and ¹³C^β chemical shifts using the program CS-ROSETTA indicated that the architecture of human hyperplastic discs protein is preserved. Mutagenesis studies confirmed that Spider Roll binds the target patch on PAK1. Additionally, Spider Roll binds to full-length PAK1 in its activated state but does not bind PAK1 when it forms an auto-inhibited conformation that blocks the Spider Roll target site. Subsequent NMR characterization of the binding of Spider Roll to PAK1 revealed a comparably small binding ‘on-rate’ constant (? 10⁵ M^− 1 s^− 1). The ability to rationally design the site of novel protein-protein interactions is an important step towards creating new proteins that are useful as therapeutics or molecular probes.     

Rate of loop formation in peptides: a simulation study

Experimental techniques with high temporal and spatial resolution extend our knowledge of how biological macromolecules self-organise and function. Here, we provide an illustration of the convergence between simulation and experiment made possible by techniques such as triplet-triplet energy transfer and fluorescence quenching with long-lifetime and fast-quenching fluorescent probes. These techniques have recently been used to determine the average time needed for two residues in a peptide or protein segment to form a contact. The timescale of this process is accessible to computer simulation, providing a microscopic interpretation of the data and yielding new insight into the disordered state of proteins. Conversely, such experimental data also provide a test of the validity of alternative choices for the molecular models used in simulations, indicating their possible deficiencies. We carried out simulations of peptides of various composition and length using several models. End-to-end contact formation rates and their dependence on peptide length agree with experimental estimates for some sequences and some force fields but not for others. The deviations are due to artefactual structuring of some peptides, which is not observed when an atomistic model for the solvation water is used. Simulations show that the observed experimental rates are compatible with considerably different distributions of the end-to-end distance; for realistic models, these are never Gaussian but indicative of a rugged energy landscape.     

Binding Hot Spot in the Weak Protein Complex of Physiological Redox Partners Yeast Cytochrome c and Cytochrome c Peroxidase

Transient protein interactions mediate many vital cellular processes such as signal transduction or intermolecular electron transfer. However, due to difficulties associated with their structural characterization, little is known about the principles governing recognition and binding in weak transient protein complexes. In particular, it has not been well established whether binding hot spots, which are frequently found in strong static complexes, also govern transient protein interactions. To address this issue, we have investigated an electron transfer complex of physiological partners from yeast: yeast iso-1-cytochrome c (Cc) and yeast cytochrome c peroxidase (CcP). Using isothermal titration calorimetry and NMR spectroscopy, we show that Cc R13 is a hot-spot residue, as R13A mutation has a strong destabilizing effect on binding. Furthermore, we employ a double-mutant cycle to illustrate that Cc R13 interacts with CcP Y39. The present results, in combination with those of earlier mutational studies, have enabled us to outline the extent of the energetically important Cc-CcP binding region. Based on our analysis, we propose that binding energy hot spots, which are prevalent in static protein complexes, could also govern transient protein interactions.     

Calculation of Proteins’ Total Side-Chain Torsional Entropy and Its Influence on Protein-Ligand Interactions

Despite the high density within a typical protein fold, the ensemble of sterically permissible side-chain repackings is vast. Here, we examine the extent of this variability that survives energetic biases due to van der Waals interactions, hydrogen bonding, salt bridges, and solvation. Monte Carlo simulations of an atomistic model exhibit thermal fluctuations among a diverse set of side-chain arrangements, even with the peptide backbone fixed in its crystallographic conformation. We have quantified the torsional entropy of this native-state ensemble, relative to that of a noninteracting reference system, for 12 small proteins. The reduction in entropy per rotatable bond due to each kind of interaction is remarkably consistent across this set of molecules. To assess the biophysical importance of these fluctuations, we have estimated side-chain entropy contributions to the binding affinity of several peptide ligands with calmodulin. Calculations for our fixed-backbone model correlate very well with experimentally determined binding entropies over a range spanning more than 80 kJ/(mol·308 K).