Absolute Protein-Protein Association Rate Constants from Flexible, Coarse-Grained Brownian Dynamics Simulations: The Role of Intermolecular Hydrodynamic Interactions in Barnase-Barstar Association

Theory and computation have long been used to rationalize the experimental association rate constants of protein-protein complexes, and Brownian dynamics (BD) simulations, in particular, have been successful in reproducing the relative rate constants of wild-type and mutant protein pairs. Missing from previous BD studies of association kinetics, however, has been the description of hydrodynamic interactions (HIs) between, and within, the diffusing proteins. Here we address this issue by rigorously including HIs in BD simulations of the barnase-barstar association reaction. We first show that even very simplified representations of the proteins—involving approximately one pseudoatom for every three residues in the protein—can provide excellent reproduction of the absolute association rate constants of wild-type and mutant protein pairs. We then show that simulations that include intermolecular HIs also produce excellent estimates of association rate constants, but, for a given reaction criterion, yield values that are decreased by ∼35–80% relative to those obtained in the absence of intermolecular HIs. The neglect of intermolecular HIs in previous BD simulation studies, therefore, is likely to have contributed to the somewhat overestimated absolute rate constants previously obtained. Consequently, intermolecular HIs could be an important component to include in accurate modeling of the kinetics of macromolecular association events

Physicochemical Properties of Cells and Their Effects on Intrinsically Disordered Proteins (IDPs)

It has long been axiomatic that a protein’s structure determines its function. Intrinsically disordered proteins (IDPs) and disordered protein regions (IDRs) defy this structure–function paradigm. They do not exhibit stable secondary and/or tertiary structures and exist as dynamic ensembles of interconverting conformers with preferred, nonrandom orientations.(1-4) The concept of IDPs and IDRs as functional biological units was initially met with skepticism. For a long time, disorder, intuitively implying chaos, had no place in our perception of orchestrated molecular events controlling cell biology. Over the past years, however, this notion has changed. Aided by findings that structural disorder constitutes an ubiquitous and abundant biological phenomenon in organisms of all phyla,(5-7) and that it is often synonymous with function,(8-11) disorder has become an integral part of modern protein biochemistry. Disorder thrives in eukaryotic signaling pathways(12) and functions as a prominent player in many regulatory processes.(13-15) Disordered proteins and protein regions determine the underlying causes of many neurodegenerative disorders and constitute the main components of amyloid fibrils.(16) They further contribute to many forms of cancer, diabetes and to cardiovascular and metabolic diseases.(17, 18) Research into disordered proteins produced significant findings and established important new concepts. On the structural side, novel experimental and computational approaches identified and described disordered protein ensembles(3, 19, 20) and led to terms such as secondary structure propensities, residual structural features, and transient long-range contacts.(1, 21) The discovery of coupled folding-and-binding reactions defined the paradigm of disorder-to-order transitions(22) and high-resolution insights into the architectures of amyloid fibrils were obtained.(23, 24) On the biological side, we learned about the unexpected intracellular stability of disordered proteins, their roles in integrating post-translational protein modifications in cell signaling and about their functions in regulatory processes ranging from transcription to cell fate decisions.(15, 25, 26) One open question remaining to be addressed is how these in vitro structural insights relate to biological in vivo effects. How do complex intracellular environments modulate the in vivo properties of disordered proteins and what are the implications for their biological functions (Figure 1)?(27-29

Parametrization of Backbone Flexibility in a Coarse-Grained Force Field for Proteins (COFFDROP) Derived from All-Atom Explicit-Solvent Molecular Dynamics Simulations of All Possible Two-Residue Peptides

Recently, we reported the parametrization of a set of coarse-grained (CG) nonbonded potential functions, derived from all-atom explicit-solvent molecular dynamics (MD) simulations of amino acid pairs and designed for use in (implicit-solvent) Brownian dynamics (BD) simulations of proteins; this force field was named COFFDROP (COarse-grained Force Field for Dynamic Representations Of Proteins). Here, we describe the extension of COFFDROP to include bonded backbone terms derived from fitting to results of explicit-solvent MD simulations of all possible two-residue peptides containing the 20 standard amino acids, with histidine modeled in both its protonated and neutral forms. The iterative Boltzmann inversion (IBI) method was used to optimize new CG potential functions for backbone-related terms by attempting to reproduce angle, dihedral, and distance probability distributions generated by the MD simulations. In a simple test of the transferability of the extended force field, the angle, dihedral, and distance probability distributions obtained from BD simulations of 56 three-residue peptides were compared to results from corresponding explicit-solvent MD simulations. In a more challenging test of the COFFDROP force field, it was used to simulate eight intrinsically disordered proteins and was shown to quite accurately reproduce the experimental hydrodynamic radii (<i>R</i>_hydro), provided that the favorable nonbonded interactions of the force field were uniformly scaled downward in magnitude. Overall, the results indicate that the COFFDROP force field is likely to find use in modeling the conformational behavior of intrinsically disordered proteins and multidomain proteins connected by flexible linkers

Molecular Dynamics Simulations of 441 Two-Residue Peptides in Aqueous Solution: Conformational Preferences and Neighboring Residue Effects with the Amber ff99SB-ildn-NMR Force Field

Understanding the intrinsic conformational preferences of amino acids and the extent to which they are modulated by neighboring residues is a key issue for developing predictive models of protein folding and stability. Here we present the results of 441 independent explicit-solvent MD simulations of all possible two-residue peptides that contain the 20 standard amino acids with histidine modeled in both its neutral and protonated states. ³<i>J</i>_HNHα coupling constants and δ_Hα chemical shifts calculated from the MD simulations correlate quite well with recently published experimental measurements for a corresponding set of two-residue peptides. Neighboring residue effects (NREs) on the average ³<i>J</i>_HNHα and δ_Hα values of adjacent residues are also reasonably well reproduced, with the large NREs exerted experimentally by aromatic residues, in particular, being accurately captured. NREs on the secondary structure preferences of adjacent amino acids have been computed and compared with corresponding effects observed in a coil library and the average β-turn preferences of all amino acid types have been determined. Finally, the intrinsic conformational preferences of histidine, and its NREs on the conformational preferences of adjacent residues, are both shown to be strongly affected by the protonation state of the imidazole ring