We combine molecular dynamics simulations and new high-field NMR experiments to describe the solution structure of the Aβ [21][22][23][24][25][26][27][28][29][30] peptide fragment that may be relevant for understanding structural mechanisms related to Alzheimer's disease. By using two different empirical force-field combinations, we provide predictions of the three-bond scalar coupling constants ( 3 J H N H α), chemical-shift values, 13 C relaxation parameters, and rotating-frame nuclear Overhauser effect spectroscopy (ROESY) crosspeaks that can then be compared directly to the same observables measured in the corresponding NMR experiment of Aβ [21][22][23][24][25][26][27][28][29][30] . We find robust prediction of the 13 C relaxation parameters and medium-range ROESY crosspeaks by using new generation TIP4P-Ew water and Amber ff99SB protein force fields, in which the NMR validates that the simulation yields both a structurally and dynamically correct ensemble over the entire Aβ 21-30 peptide. Analysis of the simulated ensemble shows that all medium-range ROE restraints are not satisfied simultaneously and demonstrates the structural diversity of the Aβ 21-30 conformations more completely than when determined from the experimental medium-range ROE restraints alone. We find that the structural ensemble of the Aβ 21-30 peptide involves a majority population (~60%) of unstructured conformers, lacking any secondary structure or persistent hydrogen-bonding networks. However, the remaining minority population contains a substantial percentage of conformers with a β-turn centered at Val24 and Gly25, as well as evidence of the Asp23 to Lys28 salt bridge important to the fibril structure. This study sets the stage for robust theoretical work on Aβ 1-40 and Aβ 1-42 , for which collection of detailed NMR data on the monomer will be more challenging because of aggregation and fibril formation on experimental timescales at physiological conditions. In addition, we believe that the interplay of modern molecular simulation and high-quality NMR experiments has reached a fruitful stage for characterizing structural ensembles of disordered peptides and proteins in general.

Myoglobin is a globular protein involved in oxygen storage and transport. No consensus yet exists on the atomic level mechanism by which oxygen and other small nonpolar ligands move between the myoglobin's buried heme, which is the ligand binding site, and surrounding solvent. This study uses room temperature molecular dynamics simulations to provide a complete atomic level picture of ligand migration in myoglobin. Multiple trajectories-providing a cumulative total of 7 s of simulation-are analyzed. Our simulation results are consistent with and tie together previous experimental findings. Specifically, we characterize: (i) Explicit full trajectories in which the CO ligand shuttles between the internal binding site and the solvent and (ii) pattern and structural origins of transient voids available for ligand migration. The computations are performed both in sperm whale myoglobin wild-type and in sperm whale V68F myoglobin mutant, which is experimentally known to slow ligand-binding kinetics. On the basis of these independent, but mutually consistent ligand migration and transient void computations, we find that there are two discrete dynamical pathways for ligand migration in myoglobin. Trajectory hops between these pathways are limited to two bottleneck regions. Ligand enters and exits the protein matrix in common identifiable portals on the protein surface. The pathways are located in the ''softer'' regions of the protein matrix and go between its helices and in its loop regions. Localized structural fluctuations are the primary physical origin of the simulated CO migration pathways inside the protein.proteins ͉ molecular dynamics ͉ kinetics ͉ structure-function relationships M yoglobin is a small globular protein used by muscle cells for oxygen storage and transport. Historically, myoglobin has been an important starting point for structural and biophysical characterization of larger proteins. It has therefore been dubbed the ''Hydrogen Atom'' of molecular biology (1). However, despite nearly half a century of experimental and theoretical studies, the precise manner in which the oxygen ligand moves between heme iron and solvent in myoglobin remains unresolved. None of the over 250 static crystallographic structures of myoglobin in the Protein Data Bank show an obvious static path between the external solvent and the heme iron. It was recognized long ago that thermal fluctuations create a web of transient structural voids and, therefore, a dynamical pathway through which small nonpolar ligands such as O 2 , CO or NO can migrate inside the molecule. Numerous experimental and theoretical works have elucidated many aspects of the ligand migration process (2-17). Some of these studies have provided evidence for a single major pathway for ligand migration that would require movement of the distal histidine ("His gate") (16,(18)(19)(20)(21)(22)(23). However, other investigators have suggested the existence of multiple pathways for ligand diffusion in myoglobin (4,5,(24)(25)(26).Reconciling these two seemingly contr...

This work investigates statistical prevalence and overall physical origins of changes in charge states of receptor proteins upon ligand binding. These changes are explored as a function of the ligand type (small molecule, protein, and nucleic acid), and distance from the binding region. Standard continuum solvent methodology is used to compute, on an equal footing, pK changes upon ligand binding for a total of 5899 ionizable residues in 20 protein-protein, 20 protein-small molecule, and 20 protein-nucleic acid high-resolution complexes. The size of the data set combined with an extensive error and sensitivity analysis allows us to make statistically justified and conservative conclusions: in 60% of all protein-small molecule, 90% of all protein-protein, and 85% of all protein-nucleic acid complexes there exists at least one ionizable residue that changes its charge state upon ligand binding at physiological conditions (pH = 6.5). Considering the most biologically relevant pH range of 4-8, the number of ionizable residues that experience substantial pK changes (DeltapK > 1.0) due to ligand binding is appreciable: on average, 6% of all ionizable residues in protein-small molecule complexes, 9% in protein-protein, and 12% in protein-nucleic acid complexes experience a substantial pK change upon ligand binding. These changes are safely above the statistical false-positive noise level. Most of the changes occur in the immediate binding interface region, where approximately one out of five ionizable residues experiences substantial pK change regardless of the ligand type. However, the physical origins of the change differ between the types: in protein-nucleic acid complexes, the pK values of interface residues are predominantly affected by electrostatic effects, whereas in protein-protein and protein-small molecule complexes, structural changes due to the induced-fit effect play an equally important role. In protein-protein and protein-nucleic acid complexes, there is a statistically significant number of substantial pK perturbations, mostly due to the induced-fit structural changes, in regions far from the binding interface.

Molecular dynamics simulations of the nucleosome core particle and its isolated DNA free in solution are reported. The simulations are based on the implicit solvent methodology and provide insights into the nature of large-scale structural fluctuations and flexibility of the nucleosomal DNA. In addition to the kinked regions previously identified in the x-ray structure of the nucleosome, the simulations support the existence of a biochemically identified distorted region of the DNA. Comparison of computed relative free energies shows that formation of the kinks is associated with little, if any, energy cost relative to a smooth, ideal conformation of the DNA superhelix. Isolated nucleosomal DNA is found to be considerably more flexible than expected for a 147 bp stretch of DNA based on its canonical persistence length of 500 A. Notably, the significant bending of the DNA observed in our simulations occurs without breaking of Watson-Crick bonds. The computed relative stability of bent conformations is sensitive to the ionic strength of the solution in the physiological range; the sensitivity suggests possible experiments that might provide further insights into the structural origins of the unusual flexibility of the DNA.