Modeling Hsp70/Hsp40 interaction by multi-scale molecular simulations and co-evolutionary sequence analysis

The interaction between the Heat Shock Proteins 70 and 40 is at the core of the ATPase regulation of the chaperone machinery that maintains protein homeostasis. However, the structural details of this fundamental interaction are still elusive and contrasting models have been proposed for the transient Hsp70/Hsp40 complexes. Here we combine molecular simulations based on both coarsegrained and atomistic models with co-evolutionary sequence analysis to shed light on this problem by focusing on the bacterial DnaK/DnaJ system. The integration of these complementary approaches resulted into a novel structural model that rationalizes previous experimental observations. We identify an evolutionary-conserved interaction surface formed by helix II of the DnaJ J-domain and a groove on lobe IIA of the DnaK nucleotide binding domain, involving the inter-domain linker.

Related papers:

title: The heat shock protein 70 family: Highly homologous proteins with overlapping and distinct functions
abstract: The human heat shock protein 70 (Hsp70) family contains at least eight homologous chaperone proteins. Endoplasmatic reticulum and mitochondria have their specific Hsp70 proteins, whereas the remaining six family members reside mainly in the cytosol and nucleus. The requirement for multiple highly homologous although different Hsp70 proteins is still far from clear, but their individual and tissue‐specific expression suggests that they are assigned distinct biological tasks. This concept is supported by the fact that mice knockout for different Hsp70 genes display remarkably discrete phenotypes. Moreover, emerging data suggest that individual Hsp70 proteins can bring about non‐overlapping and chaperone‐independent functions essential for growth and survival of cancer cells. This review summarizes our present knowledge of the individual members of human Hsp70 family and elaborate on the functional differences between the cytosolic/nuclear representatives.

title: Molecular chaperones in protein folding and proteostasis
abstract: Most proteins must fold into defined three-dimensional structures to gain functional activity. But in the cellular environment, newly synthesized proteins are at great risk of aberrant folding and aggregation, potentially forming toxic species. To avoid these dangers, cells invest in a complex network of molecular chaperones, which use ingenious mechanisms to prevent aggregation and promote efficient folding. Because protein molecules are highly dynamic, constant chaperone surveillance is required to ensure protein homeostasis (proteostasis). Recent advances suggest that an age-related decline in proteostasis capacity allows the manifestation of various protein-aggregation diseases, including Alzheimer's disease and Parkinson's disease. Interventions in these and numerous other pathological states may spring from a detailed understanding of the pathways underlying proteome maintenance.

title: Hsp70 chaperone dynamics and molecular mechanism.
abstract: The chaperone functions of heat shock protein (Hsp)70 involve an allosteric control mechanism between the nucleotide-binding domain (NBD) and polypeptide substrate-binding domain (SBD): ATP binding and hydrolysis regulates the affinity for polypeptides, and polypeptide binding accelerates ATP hydrolysis. These data suggest that Hsp70s exist in at least two conformational states. Although structural information on the conformation with high affinity for polypeptides has been available for several years, the conformation with an open polypeptide binding cleft was elucidated only recently. In addition, other biophysical studies have revealed a more dynamic picture of Hsp70s, shedding light on the molecular mechanism by which Hsp70s assist protein folding. In this review recent insights into the structure and mechanism of Hsp70s are discussed.

title: The HSP70 chaperone machinery: J proteins as drivers of functional specificity
abstract: Heat shock 70 kDa proteins (HSP70s) are ubiquitous molecular chaperones that function in a myriad of biological processes, modulating polypeptide folding, degradation and translocation across membranes, and protein–protein interactions. This multitude of roles is not easily reconciled with the universality of the activity of HSP70s in ATP-dependent client protein-binding and release cycles. Much of the functional diversity of the HSP70s is driven by a diverse class of cofactors: J proteins. Often, multiple J proteins function with a single HSP70. Some target HSP70 activity to clients at precise locations in cells and others bind client proteins directly, thereby delivering specific clients to HSP70 and directly determining their fate.

title: Hsp70 chaperones accelerate protein translocation and the unfolding of stable protein aggregates by entropic pulling.
abstract: Hsp70s are highly conserved ATPase molecular chaperones mediating the correct folding of de novo synthesized proteins, the translocation of proteins across membranes, the disassembly of some native protein oligomers, and the active unfolding and disassembly of stress-induced protein aggregates. Here, we bring thermodynamic arguments and biochemical evidences for a unifying mechanism named entropic pulling, based on entropy loss due to excluded-volume effects, by which Hsp70 molecules can convert the energy of ATP hydrolysis into a force capable of accelerating the local unfolding of various protein substrates and, thus, perform disparate cellular functions. By means of entropic pulling, individual Hsp70 molecules can accelerate unfolding and pulling of translocating polypeptides into mitochondria in the absence of a molecular fulcrum, thus settling former contradictions between the power-stroke and the Brownian ratchet models for Hsp70-mediated protein translocation across membranes. Moreover, in a very different context devoid of membrane and components of the import pore, the same physical principles apply to the forceful unfolding, solubilization, and assisted native refolding of stable protein aggregates by individual Hsp70 molecules, thus providing a mechanism for Hsp70-mediated protein disaggregation.

title: The kinetic parameters and energy cost of the Hsp70 chaperone as a polypeptide unfoldase.
abstract: Hsp70-Hsp40-NEF and possibly Hsp100 are the only known molecular chaperones that can use the energy of ATP to convert stably pre-aggregated polypeptides into natively refolded proteins. However, the kinetic parameters and ATP costs have remained elusive because refolding reactions have only been successful with a molar excess of chaperones over their polypeptide substrates. Here we describe a stable, misfolded luciferase species that can be efficiently renatured by substoichiometric amounts of bacterial Hsp70-Hsp40-NEF. The reactivation rates increased with substrate concentration and followed saturation kinetics, thus allowing the determination of apparent V(max)' and K(m)' values for a chaperone-mediated renaturation reaction for the first time. Under the in vitro conditions used, one Hsp70 molecule consumed five ATPs to effectively unfold a single misfolded protein into an intermediate that, upon chaperone dissociation, spontaneously refolded to the native state, a process with an ATP cost a thousand times lower than expected for protein degradation and resynthesis.

title: Allostery in the Hsp70 chaperone proteins.
abstract: Heat shock 70-kDa (Hsp70) chaperones are essential to in vivo protein folding, protein transport, and protein re-folding. They carry out these activities using repeated cycles of binding and release of client proteins. This process is under allosteric control of nucleotide binding and hydrolysis. X-ray crystallography, NMR spectroscopy, and other biophysical techniques have contributed much to the understanding of the allosteric mechanism linking these activities and the effect of co-chaperones on this mechanism. In this chapter these findings are critically reviewed. Studies on the allosteric mechanisms of Hsp70 have gained enhanced urgency, as recent studies have implicated this chaperone as a potential drug target in diseases such as Alzheimer's and cancer. Recent approaches to combat these diseases through interference with the Hsp70 allosteric mechanism are discussed.

title: Solution conformation of wild-type E. coli Hsp70 (DnaK) chaperone complexed with ADP and substrate
abstract: DnaK is the canonical Hsp70 molecular chaperone protein from Escherichia coli. Like other Hsp70s, DnaK comprises two main domains: a 44-kDa N-terminal nucleotide-binding domain (NBD) that contains ATPase activity, and a 25-kDa substrate-binding domain (SBD) that harbors the substrate-binding site. Here, we report an experimental structure for wild-type, full-length DnaK, complexed with the peptide NRLLLTG and with ADP. It was obtained in aqueous solution by using NMR residual dipolar coupling and spin labeling methods and is based on available crystal structures for the isolated NBD and SBD. By using dynamics methods, we determine that the NBD and SBD are loosely linked and can move in cones of ±35° with respect to each other. The linker region between the domains is a dynamic random coil. Nevertheless, an average structure can be defined. This structure places the SBD in close proximity of subdomain IA of the NBD and suggests that the SBD collides with the NBD at this area to establish allosteric communication.

title: Structure and dynamics of the ATP-bound open conformation of Hsp70 chaperones.
abstract: Central to the chaperone function of Hsp70s is the transition between open and closed conformations of their polypeptide substrate binding domain (SBD), which is regulated through an allosteric mechanism via ATP binding and hydrolysis in their nucleotide binding domain (NBD). Although the structure of the closed conformation of Hsp70s is well studied, the open conformation has remained elusive. Here, we report on the 2.4 Å crystal structure of the ATP-bound open conformation of the Escherichia coli Hsp70 homolog DnaK. In the open DnaK structure, the β sheet and α-helical lid subdomains of the SBD are detached from one another and docked to different faces of the NBD. The contacts between the β sheet subdomain and the NBD reveal the mechanism of allosteric regulation. In addition, we demonstrate that docking of the β sheet and α-helical lid subdomains to the NBD is a sequential process influenced by peptide and protein substrates.

title: Allosteric opening of the polypeptide-binding site when an Hsp70 binds ATP
abstract: The 70-kilodalton (kDa) heat-shock proteins (Hsp70s) are ubiquitous molecular chaperones essential for cellular protein folding and proteostasis. Each Hsp70 has two functional domains: a nucleotide-binding domain (NBD), which binds and hydrolyzes ATP, and a substrate-binding domain (SBD), which binds extended polypeptides. NBD and SBD interact little when in the presence of ADP; however, ATP binding allosterically couples the polypeptide- and ATP-binding sites. ATP binding promotes polypeptide release; polypeptide rebinding stimulates ATP hydrolysis. This allosteric coupling is poorly understood. Here we present the crystal structure of an intact ATP-bound Hsp70 from Escherichia coli at 1.96-Å resolution. The ATP-bound NBD adopts a unique conformation, forming extensive interfaces with an SBD that has changed radically, having its α-helical lid displaced and the polypeptide-binding channel of its β-subdomain restructured. These conformational changes, together with our biochemical assays, provide a structural explanation for allosteric coupling in Hsp70 activity.

title: NMR structure of the J-domain and the Gly/Phe-rich region of the Escherichia coli DnaJ chaperone.
abstract: The recombinant N-terminal 107-amino acid polypeptide fragment 2-108 of the DnaJ molecular chaperone of Escherichia coli, which contains the J-domain (residues 2 to 76) and the Gly/Phe-rich region (residues 77 to 108), was uniformly labeled with nitrogen-15 and carbon-13. The complete NMR solution structure of the J-domain was determined with the program DIANA on the basis of 682 nuclear Overhauser enhancement (NOE) upper distance limits and 180 dihedral angle constraints. It contains three well-defined helices comprising residues 6 to 10, 18 to 32 and 41 to 57, and a fourth helix, consisting of residues 61 to 68, which is well defined as a regular secondary structure but for which the location relative to the remainder of the molecule is not precisely determined. The helices II and III form an antiparallel helical coiled-coil. Helix I is approximately parallel to the plane defined by the helices II and III and runs from the carboxy-terminal end of the helix III to the center of helix II. Helix IV is positioned near the carboxy-terminal end of helix III and is on the same side of the coiled coil as helix I, but it is oriented approximately perpendicular to the plane of the helices II and III. This novel alpha-protein topology leads to formation of a hydrophobic core involving side-chains of all four helices. A strong correlation is seen between the extent of sequence-conservation of hydrophobic residues in the family of J-domain homologues, and the structural organization of the hydrophobic core in these proteins. The residues which have key roles for the specificity of the interaction of DnaJ-like proteins with their corresponding Hsp70 counterparts are located on the outer surfaces of the helices II and III, and in the loop connecting these two helices. Measurements of backbone amide proton exchange rates, 15N spin relaxation times and heteronuclear 15N {1H} NOEs provided additional insights into local conformational equilibria and internal rate processes in DnaJ(2-108). In the Gly/Phe-rich region, which is poorly ordered in the NMR solution structure and does not form a globular core, the polypeptide segment 90 to 103 differs from the segments 77 to 89 and 104 to 108 by reduced local flexibility. Considering that this same segment shows sequence conservation with corresponding segments in the Gly/Phe-rich regions of other DnaJ-like proteins, its reduced flexibility may be directly linked to the formation of the ternary DnaJ-DnaK-polypeptide complex.

title: Role of the J-domain in the cooperation of Hsp40 with Hsp70.
abstract: The Escherichia coli Hsp40 DnaJ and Hsp70 DnaK cooperate in the binding of proteins at intermediate stages of folding, assembly, and translocation across membranes. Binding of protein substrates to the DnaK C-terminal domain is controlled by ATP binding and hydrolysis in the N-terminal ATPase domain. The interaction of DnaJ with DnaK is mediated at least in part by the highly conserved N-terminal J-domain of DnaJ that includes residues 2-75. Heteronuclear NMR experiments with uniformly 15N-enriched DnaJ2-75 indicate that the chemical environment of residues located in helix II and the flanking loops is perturbed on interaction with DnaK or a truncated DnaK molecule, DnaK2-388. NMR signals corresponding to these residues broaden and exhibit changes in chemical shifts in the presence of DnaK(MgADP). Addition of MgATP largely reversed the broadening, indicating that NMR signals of DnaJ2-75 respond to ATP-dependent changes in DnaK. The J-domain interaction is localized to the ATPase domain of DnaK and is likely to be dominated by electrostatic interactions. The results suggest that the J-domain tethers DnaK to DnaJ-bound substrates, which DnaK then binds with its C-terminal peptide-binding domain.

title: Investigation of the interaction between DnaK and DnaJ by surface plasmon resonance spectroscopy.
abstract: Hsp70 chaperones assist protein folding through ATP-regulated transient association with substrates. Substrate binding by Hsp70 is controlled by DnaJ co-chaperones which stimulate Hsp70 to hydrolyze ATP and, consequently, to close its substrate binding cavity allowing trapping of substrates. We analyzed the interaction of the Escherichia coli Hsp70 homologue, DnaK, with DnaJ using surface plasmon resonance (SPR) spectroscopy. Resonance signals of complex kinetic characteristics were detected when DnaK was passed over a sensor chip with coupled DnaJ. This interaction was specific as it was not detected with a functionally defective DnaJ mutant protein, DnaJ259, that carries a mutation in the HPD signature motif of the conserved J-domain. Detectable DnaK-DnaJ interaction required ATP hydrolysis by DnaK and was competitively inhibited by chaperone substrates of DnaK. For DnaK mutant proteins with amino acid substitutions in the substrate binding cavity that affect substrate binding, the strength of detected interaction with DnaJ decreased proportionally with increased strength of the substrate binding defects. These findings indicate that the detected response signals resulted from DnaJ and ATP hydrolysis-dependent association of DnaJ as substrate for DnaK. Although not considered as physiologically relevant, this association allowed us to experimentally unravel the mechanism of DnaJ action. Accordingly, DnaJ stimulates ATP hydrolysis only after association of a substrate with the substrate binding cavity of DnaK. Further analysis revealed that this coupling mechanism required the J-domain of DnaJ and was also functional for natural DnaK substrates, and thus is central to the mechanism of action of the DnaK chaperone system.

title: Interaction of the Hsp70 molecular chaperone, DnaK, with its cochaperone DnaJ.
abstract: Chaperones of the Hsp70 family bind to unfolded or partially folded polypeptides to facilitate many cellular processes. ATP hydrolysis and substrate binding, the two key molecular activities of this chaperone, are modulated by the cochaperone DnaJ. By using both genetic and biochemical approaches, we provide evidence that DnaJ binds to at least two sites on the Escherichia coli Hsp70 family member DnaK: under the ATPase domain in a cleft between its two subdomains and at or near the pocket of substrate binding. The lower cleft of the ATPase domain is defined as a binding pocket for the J-domain because (i) a DnaK mutation located in this cleft (R167H) is an allele-specific suppressor of the binding defect of the DnaJ mutation, D35N and (ii) alanine substitution of two residues close to R167 in the crystal structure, N170A and T173A, significantly decrease DnaJ binding. A second binding determinant is likely to be in the substrate-binding domain because some DnaK mutations in the vicinity of the substrate-binding pocket are defective in either the affinity (G400D, G539D) or rate (D526N) of both peptide and DnaJ binding to DnaK. Binding of DnaJ may propagate conformational changes to the nearby ATPase catalytic center and substrate-binding sites as well as facilitate communication between these two domains to alter the molecular properties of Hsp70.

title: A Conserved HPD Sequence of the J-domain Is Necessary for YDJ1 Stimulation of Hsp70 ATPase Activity at a Site Distinct from Substrate Binding (*)
abstract: The 46-kDa protein YDJ1 is one of several known yeast homologues of the Escherichia coli DnaJ protein. Like all J homologues, it shares homology with the highly conserved NH-terminal “J-domain” of DnaJ. A component of the DnaK (Hsp70) chaperone machinery that mediates protein folding, DnaJ is necessary for survival at elevated temperatures. It stimulates ATP hydrolysis by DnaK and effects the release of DnaK-bound polypeptides. Previous genetic and biochemical studies indicate that the J-domain is necessary for these functions. Using peptides corresponding to J-domain sequence, we show that a peptide containing the highly conserved His-Pro-Asp sequence at positions 34-36 in the J-domain competes off YDJ1 stimulation of Hsp70 ATPase activity. Inhibitory concentrations of peptide do not prevent binding of folding substrates, therefore YDJ1 must interact with Hsp70 at a site distinct from that for substrate binding. This interaction is critical for Hsp70 activity, since a mutant YDJ1 protein harboring a H34Q change (ydj1Q34) stimulates neither Hsp70 ATPase nor substrate release. The importance of the proper function of this region of the protein is supported by the poor growth and temperature-sensitive phenotype of yeast expressing ydj1Q34.

title: Hsp70 chaperones are non-equilibrium machines that achieve ultra-affinity by energy consumption
abstract: 70-kDa Heat shock proteins are ATP-driven molecular chaperones that perform a myriad of essential cellular tasks. Although structural and biochemical studies have shed some light on their functional mechanism, the fundamental issue of the role of energy consumption, due to ATP-hydrolysis, has remained unaddressed. Here we establish a clear connection between the non-equilibrium nature of Hsp70, due to ATP hydrolysis, and the determining feature of its function, namely its high affinity for its substrates. Energy consumption can indeed decrease the dissociation constant of the chaperone-substrate complex by several orders of magnitude with respect to an equilibrium scenario. We find that the biochemical requirements for observing such ultra-affinity coincide with the physiological conditions in the cell. Our results rationalize several experimental observations and pave the way for further analysis of non-equilibrium effects underlying chaperone functions.

title: The conformational dynamics of the mitochondrial Hsp70 chaperone.
abstract: Heat shock proteins 70 (Hsp70) represent a ubiquitous and conserved family of molecular chaperones involved in a plethora of cellular processes. The dynamics of their ATP hydrolysis-driven and cochaperone-regulated conformational cycle are poorly understood. We used fluorescence spectroscopy to analyze, in real time and at single-molecule resolution, the effects of nucleotides and cochaperones on the conformation of Ssc1, a mitochondrial member of the family. We report that the conformation of its ADP state is unexpectedly heterogeneous, in contrast to a uniform ATP state. Substrates are actively involved in determining the conformation of Ssc1. The J protein Mdj1 does not interact transiently with the chaperone, as generally believed, but rather is released slowly upon ATP hydrolysis. Analysis of the major bacterial Hsp70 revealed important differences between highly homologous members of the family, possibly explaining tuning of Hsp70 chaperones to meet specific functions in different organisms and cellular compartments.

title: Heat shock protein 70 kDa chaperone/DnaJ cochaperone complex employs an unusual dynamic interface
abstract: The heat shock protein 70 kDa (Hsp70)/DnaJ/nucleotide exchange factor system assists in intracellular protein (re)folding. Using solution NMR, we obtained a three-dimensional structure for a 75-kDa Hsp70–DnaJ complex in the ADP state, loaded with substrate peptide. We establish that the J domain (residues 1–70) binds with its positively charged helix II to a negatively charged loop in the Hsp70 nucleotide-binding domain. The complex shows an unusual “tethered” binding mode which is stoichiometric and saturable, but which has a dynamic interface. The complex represents part of a triple complex of Hsp70 and DnaJ both bound to substrate protein. Mutagenesis data indicate that the interface is also of relevance for the interaction of Hsp70 and DnaJ in the ATP state. The solution complex is completely different from a crystal structure of a disulfide-linked complex of homologous proteins [Jiang, et al. (2007) .

title: Scanning mutagenesis identifies amino acid residues essential for the in vivo activity of the Escherichia coli DnaJ (Hsp40) J-domain.
abstract: The DnaJ (Hsp40) cochaperone regulates the DnaK (Hsp70) chaperone by accelerating ATP hydrolysis in a cycle closely linked to substrate binding and release. The J-domain, the signature motif of the Hsp40 family, orchestrates interaction with the DnaK ATPase domain. We studied the J-domain by creating 42 mutant E. coli DnaJ variants and examining their phenotypes in various separate in vivo assays, namely, bacterial growth at low and high temperatures, motility, and propagation of bacteriophage lambda. Most mutants studied behaved like wild type in all assays. In addition to the (33)HisProAsp(35) (HPD) tripeptide found in all known functional J-domains, our study uncovered three new single substitution mutations (Y25A, K26A, and F47A) that totally abolish J-domain function. Furthermore, two glycine substitution mutants in an exposed flexible loop (R36G, N37G) showed partial loss of J-domain function alone and complete loss of function as a triple (RNQ-GGG) mutant coupled with the phenotypically silent Q38G. Interestingly, all the essential residues map to a small region on the same solvent-exposed face of the J-domain. Engineered mutations in the corresponding residues of the human Hdj1 J-domain grafted in E. coli DnaJ also resulted in loss of function, suggesting an evolutionarily conserved interaction surface. We propose that these clustered residues impart critical sequence determinants necessary for J-domain catalytic activity and reversible contact interface with the DnaK ATPase domain.

title: Structural Features Required for the Interaction of the Hsp70 Molecular Chaperone DnaK with Its Cochaperone DnaJ*
abstract: Hsp70 family members together with their Hsp40 cochaperones function as molecular chaperones, using an ATP-controlled cycle of polypeptide binding and release to mediate protein folding. Hsp40 plays a key role in the chaperone reaction by stimulating the ATPase activity and activating the substrate binding of Hsp70. We have explored the interaction between the Escherichia coli Hsp70 family member, DnaK, and its cochaperone partner DnaJ. Our data show that the binding of ATP, subsequent conformational changes in DnaK, and DnaJ-stimulated ATP hydrolysis are all required for the formation of a DnaK-DnaJ complex as monitored by Biacore analysis. In addition, our data imply that the interaction of the J-domain with DnaK depends on the substrate binding state of DnaK.

title: Structural basis of J cochaperone binding and regulation of Hsp70.
abstract: The many protein processing reactions of the ATP-hydrolyzing Hsp70s are regulated by J cochaperones, which contain J domains that stimulate Hsp70 ATPase activity and accessory domains that present protein substrates to Hsp70s. We report the structure of a J domain complexed with a J responsive portion of a mammalian Hsp70. The J domain activates ATPase activity by directing the linker that connects the Hsp70 nucleotide binding domain (NBD) and substrate binding domain (SBD) toward a hydrophobic patch on the NBD surface. Binding of the J domain to Hsp70 displaces the SBD from the NBD, which may allow the SBD flexibility to capture diverse substrates. Unlike prokaryotic Hsp70, the SBD and NBD of the mammalian chaperone interact in the ADP state. Thus, although both nucleotides and J cochaperones modulate Hsp70 NBD:linker and NBD:SBD interactions, the intrinsic persistence of those interactions differs in different Hsp70s and this may optimize their activities for different cellular roles.

title: Evaluation of competing J domain:Hsp70 complex models in light of existing mutational and NMR data
abstract: The work by Ahmad et al. (1) presented an NMR-based model for a bacterial DnaJ J domain:DnaK(Hsp70):ADP complex that differs from our crystal structure of a disulfide-linked bovine Hsc70:auxilin J domain complex (2). The work by Ahmad et al. (1) claimed that their model can better account for published mutational data, that their model is in better agreement with a previous NMR study of a bacterial DnaJ J domain:DnaK complex (3), and that the cross-link in our structure renders it irrelevant for understanding J domain:Hsp70 interactions or mechanism. However, bacterial DnaK and mammalian Hsc70 bind J proteins very differently. Bacterial DnaK interacts primarily with helix II of DnaJ (3), whereas bovine Hsc70 interacts primarily with helix III and the histidine-proline-aspartate (HPD) loop (4). In our structure, helix III and the HPD loop of the bovine auxilin J domain form the primary interaction surface with bovine Hsc70, and our structure is in excellent agreement with an NMR study of a bovine Hsc70:polyoma virus T-ag J domain complex (4), which is a more appropriate comparison than the bacterial complex. Our structure can better account for the mutational data: 8 of 9 DnaJ J domain residues with mutations that affect J function, 12 of 14 residues in viral T-ag J domains with mutations that affect J function, and 10 of 12 residues/regions in DnaK/Hsc70 with mutations that affect interaction with J correspond to residues that are close (<8 A) to the binding partner in our complex. In contrast, in the model by Ahmad et al. (1), only 1 of 9 functionally critical DnaJ residues, 1 of 12 DnaK residues/regions, and 0 of 14 viral T-ag J domain residues are within 15 A of residues in the binding partner that were spin-labeled to provide NMR distance constraints. The latter includes a spin-labeled DnaK residue (V210) that is in the center of the loop proposed in the model by Ahmad et al. (1) to form the binding surface for DnaJ. Cross-linking is frequently used to capture transient or dynamic complexes for crystallization, with several validated examples in the Protein Data Bank (i.e., 1RTD, 1KLN, and 1S6V). Evidence for the relevance of our structure is that (i) the J domain has ATPase stimulating effects on the cross-linked Hsc70 that exactly recapitulate the ATPase stimulation seen for the two uncross-linked species (2), (ii) the structure predicts the location of Hsc70 mutations that disrupt its ability to be stimulated by the J protein (2), (iii) the structure indicates how the J domain stimulates Hsc70 ATPase activity through conformational effects on the linker between the Hsc70 nucleotide and protein binding domains (2), a mechanism that is gaining increasing acceptance, and (iv) the structure is consistent with a cryo-EM structure of a clathrin cage:auxilin:Hsc70:ATP complex. Space and reference limitations preclude additional elaboration of these issues or listing of all relevant references here. However, interested readers can access a detailed discussion along with tables summarizing the NMR and mutational data and references in ref. 5.

title: Reply to Sousa et al.: Evaluation of competing J domain:Hsp70 complex models in light of methods used
abstract: We respond to Sousa et al. (1) who compare a crystal structure (2) of a disulfide-linked adduct of human Hsp70 nucleotide binding domain (NBD) and auxilin J domain, with the solution complex (3) of Escherichia coli Hsp70 NBD, substrate binding domain (SBD), substrate and DnaJ J domain. The molecular models are very different despite >60% sequence homology (the coordinates in ref. 3 are available on request).

title: Coarse-grained models for simulations of multiprotein complexes: application to ubiquitin binding.
abstract: We develop coarse-grained models and effective energy functions for simulating thermodynamic and structural properties of multiprotein complexes with relatively low binding affinity (K(d) >1 microM) and apply them to binding of Vps27 to membrane-tethered ubiquitin. Folded protein domains are represented as rigid bodies. The interactions between the domains are treated at the residue level with amino-acid-dependent pair potentials and Debye-Hückel-type electrostatic interactions. Flexible linker peptides connecting rigid protein domains are represented as amino acid beads on a polymer with appropriate stretching, bending, and torsion-angle potentials. In simulations of membrane-attached protein complexes, interactions between amino acids and the membrane are described by residue-dependent short-range potentials and long-range electrostatics. We parameterize the energy functions by fitting the osmotic second virial coefficient of lysozyme and the binding affinity of the ubiquitin-CUE complex. For validation, extensive replica-exchange Monte Carlo simulations are performed of various protein complexes. Binding affinities for these complexes are in good agreement with the experimental data. The simulated structures are clustered on the basis of distance matrices between two proteins and ranked according to cluster population. In approximately 70% of the complexes, the distance root-mean-square is less than 5 A from the experimental structures. In approximately 90% of the complexes, the binding interfaces on both proteins are predicted correctly, and in all other cases at least one interface is correct. Transient and nonspecifically bound structures are also observed. With the validated model, we simulate the interaction between the Vps27 multiprotein complex and a membrane-tethered ubiquitin. Ubiquitin is found to bind preferentially to the two UIM domains of Vps27, but transient interactions between ubiquitin and the VHS and FYVE domains are observed as well. These specific and nonspecific interactions are found to be positively cooperative, resulting in a substantial enhancement of the overall binding affinity beyond the approximately 300 microM of the specific domains. We also find that the interactions between ubiquitin and Vps27 are highly dynamic, with conformational rearrangements enabling binding of Vps27 to diverse targets as part of the multivesicular-body protein-sorting pathway.

title: Replica exchange simulations of transient encounter complexes in protein–protein association
abstract: Recent paramagnetic relaxation enhancement (PRE) studies on several weakly interacting protein complexes have unequivocally demonstrated the existence of transient encounter complexes. Here, we present a computational method to study protein–protein binding by creating equilibrium ensembles that include both specific and nonspecific protein complexes. In a joint analysis of simulation and experiment we explore the physical nature and underlying physicochemical characteristics of encounter complexes involving three protein–protein interactions of the bacterial phosphotransferase system. Replica exchange Monte Carlo simulations using a coarse-grained energy function recover the structures of the specific complexes and produce binding affinities in good agreement with experiment. Together with the specific complex, a relatively small number of distinct nonspecific complexes largely accounts for the measured PRE data. The combined relative population of the latter is less than ∼10%. The binding interfaces of the specific and nonspecific complexes differ primarily in size but exhibit similar amino acid compositions. We find that the overall funnel-shaped energy landscape of complex formation is dominated by the specific complex, a small number of structured nonspecific complexes, and a diffuse cloud of loosely bound complexes connecting the specific and nonspecific binding sites with each other and the unbound state. Nonspecific complexes may not only accelerate the binding kinetics by enhancing the rate of success of random diffusional encounters but also play a role in protein function as alternative binding modes.

title: SAXS ensemble refinement of ESCRT-III CHMP3 conformational transitions.
abstract: We developed and implemented an ensemble-refinement method to study dynamic biomolecular assemblies with intrinsically disordered segments. Data from small angle X-ray scattering (SAXS) experiments and from coarse-grained molecular simulations were combined by using a maximum-entropy approach. The method was applied to CHMP3 of ESCRT-III, a protein with multiple helical domains separated by flexible linkers. Based on recent SAXS data by Lata et al. ( 2008), we constructed ensembles of CHMP3 at low- and high-salt concentration to characterize its closed autoinhibited state and open active state. At low salt, helix α(5) is bound to the tip of helices α(1) and α(2), in excellent agreement with a recent crystal structure. Helix α(6) remains free in solution and does not appear to be part of the autoinhibitory complex. The simulation-based ensemble refinement is general and effectively increases the resolution of SAXS beyond shape information to atomically detailed structures.

title: The J-domain of Hsp40 couples ATP hydrolysis to substrate capture in Hsp70.
abstract: The Escherichia coli Hsp40 DnaJ uses its J-domain to target substrate polypeptides for binding to the Hsp70 DnaK, but the mechanism of J-domain function has been obscured by a substrate-like interaction between DnaJ and DnaK. ATP hydrolysis in DnaK is associated with a conformational change that captures the substrate, and both DnaJ and substrate can stimulate ATP hydrolysis. However, substrates cannot trigger capture by DnaK in the presence of ATP, and substrates stimulate a DnaK conformational change that is uncoupled from ATP hydrolysis. The role of the J-domain was examined using the fluorescent derivative of a fusion protein composed of the J-domain and a DnaK-binding peptide. In the absence of ATP, DnaK-binding affinity of the fusion protein is similar to that of the unfused peptide. However, in the presence of ATP, the affinity of the fusion protein is dramatically increased, which is opposite to the decrease in DnaK affinity typically exhibited by peptides. Binding of a fusion protein that contains a defective J-domain is insensitive to ATP. According to results from isothermal titration calorimetry, the J-domain binds to the DnaK ATPase domain with weak affinity (K(D) = 23 microM at 20 degrees C). The interaction is characterized by a positive enthalpy, small heat capacity change (DeltaC(p)= -33 kcal mol(-1)), and increasing binding affinity for increasing temperatures in the physiological range. In conditions that support binding of the J-domain to the ATPase domain, the J-domain accelerates ATP hydrolysis and a simultaneous conformational change in DnaK that is associated with peptide capture. The defective J-domain is inactive, despite the fact that it binds to the DnaK ATPase domain with higher than wild-type affinity. The results are most consistent with an allosteric mechanism of J-domain action in which the J-domain couples ATP hydrolysis to peptide capture by accelerating ATP hydrolysis and delaying DnaK closure until ATP is hydrolyzed.

title: Emerging methods in protein co-evolution
abstract: Co-evolution is a fundamental component of the theory of evolution and is essential for understanding the relationships between species in complex ecological networks. A wide range of co-evolution-inspired computational methods has been designed to predict molecular interactions, but it is only recently that important advances have been made. Breakthroughs in the handling of phylogenetic information and in disentangling indirect relationships have resulted in an improved capacity to predict interactions between proteins and contacts between different protein residues. Here, we review the main co-evolution-based computational approaches, their theoretical basis, potential applications and foreseeable developments.

title: Protein structure prediction from sequence variation
abstract: Genomic sequences contain rich evolutionary information about functional constraints on macromolecules such as proteins. This information can be efficiently mined to detect evolutionary couplings between residues in proteins and address the long-standing challenge to compute protein three-dimensional structures from amino acid sequences. Substantial progress has recently been made on this problem owing to the explosive growth in available sequences and the application of global statistical methods. In addition to three-dimensional structure, the improved understanding of covariation may help identify functional residues involved in ligand binding, protein-complex formation and conformational changes. We expect computation of covariation patterns to complement experimental structural biology in elucidating the full spectrum of protein structures, their functional interactions and evolutionary dynamics.

title: Direct-coupling analysis of residue coevolution captures native contacts across many protein families
abstract: The similarity in the three-dimensional structures of homologous proteins imposes strong constraints on their sequence variability. It has long been suggested that the resulting correlations among amino acid compositions at different sequence positions can be exploited to infer spatial contacts within the tertiary protein structure. Crucial to this inference is the ability to disentangle direct and indirect correlations, as accomplished by the recently introduced direct-coupling analysis (DCA). Here we develop a computationally efficient implementation of DCA, which allows us to evaluate the accuracy of contact prediction by DCA for a large number of protein domains, based purely on sequence information. DCA is shown to yield a large number of correctly predicted contacts, recapitulating the global structure of the contact map for the majority of the protein domains examined. Furthermore, our analysis captures clear signals beyond intradomain residue contacts, arising, e.g., from alternative protein conformations, ligand-mediated residue couplings, and interdomain interactions in protein oligomers. Our findings suggest that contacts predicted by DCA can be used as a reliable guide to facilitate computational predictions of alternative protein conformations, protein complex formation, and even the de novo prediction of protein domain structures, contingent on the existence of a large number of homologous sequences which are being rapidly made available due to advances in genome sequencing.

title: Identification of direct residue contacts in protein–protein interaction by message passing
abstract: Understanding the molecular determinants of specificity in protein–protein interaction is an outstanding challenge of postgenome biology. The availability of large protein databases generated from sequences of hundreds of bacterial genomes enables various statistical approaches to this problem. In this context covariance-based methods have been used to identify correlation between amino acid positions in interacting proteins. However, these methods have an important shortcoming, in that they cannot distinguish between directly and indirectly correlated residues. We developed a method that combines covariance analysis with global inference analysis, adopted from use in statistical physics. Applied to a set of >2,500 representatives of the bacterial two-component signal transduction system, the combination of covariance with global inference successfully and robustly identified residue pairs that are proximal in space without resorting to ad hoc tuning parameters, both for heterointeractions between sensor kinase (SK) and response regulator (RR) proteins and for homointeractions between RR proteins. The spectacular success of this approach illustrates the effectiveness of the global inference approach in identifying direct interaction based on sequence information alone. We expect this method to be applicable soon to interaction surfaces between proteins present in only 1 copy per genome as the number of sequenced genomes continues to expand. Use of this method could significantly increase the potential targets for therapeutic intervention, shed light on the mechanism of protein–protein interaction, and establish the foundation for the accurate prediction of interacting protein partners.

title: Sequence co-evolution gives 3D contacts and structures of protein complexes
abstract: Protein–protein interactions are fundamental to many biological processes. Experimental screens have identified tens of thousands of interactions, and structural biology has provided detailed functional insight for select 3D protein complexes. An alternative rich source of information about protein interactions is the evolutionary sequence record. Building on earlier work, we show that analysis of correlated evolutionary sequence changes across proteins identifies residues that are close in space with sufficient accuracy to determine the three-dimensional structure of the protein complexes. We evaluate prediction performance in blinded tests on 76 complexes of known 3D structure, predict protein–protein contacts in 32 complexes of unknown structure, and demonstrate how evolutionary couplings can be used to distinguish between interacting and non-interacting protein pairs in a large complex. With the current growth of sequences, we expect that the method can be generalized to genome-wide elucidation of protein–protein interaction networks and used for interaction predictions at residue resolution.

title: Robust and accurate prediction of residue–residue interactions across protein interfaces using evolutionary information
abstract: Do the amino acid sequence identities of residues that make contact across protein interfaces covary during evolution? If so, such covariance could be used to predict contacts across interfaces and assemble models of biological complexes. We find that residue pairs identified using a pseudo-likelihood-based method to covary across protein–protein interfaces in the 50S ribosomal unit and 28 additional bacterial protein complexes with known structure are almost always in contact in the complex, provided that the number of aligned sequences is greater than the average length of the two proteins. We use this method to make subunit contact predictions for an additional 36 protein complexes with unknown structures, and present models based on these predictions for the tripartite ATP-independent periplasmic (TRAP) transporter, the tripartite efflux system, the pyruvate formate lyase-activating enzyme complex, and the methionine ABC transporter.

title: Allosteric Regulation of Hsp70 Chaperones Involves a Conserved Interdomain Linker*
abstract: The 70-kDa heat shock proteins (Hsp70) are essential members of the cellular chaperone machinery that assists protein-folding processes. To perform their functions Hsp70 chaperones toggle between two nucleotide-controlled conformational states. ATP binding to the ATPase domain triggers the transition to the low affinity state of the substrate-binding domain, while substrate binding to the substrate-binding domain in synergism with the action of a J-domain-containing cochaperone stimulates ATP hydrolysis and thereby transition to the high affinity state. Thus, ATPase and substrate-binding domains mutually affect each other through an allosteric control mechanism, the basis of which is largely unknown. In this study we identified two positively charged, surface-exposed residues in the ATPase domain and a negatively charged residue in the linker connecting both domains that are important for interdomain communication. Furthermore, we demonstrate that the linker alone is sufficient to stimulate the ATPase activity, an ability that is lost upon amino acid replacement. The linker therefore is most likely the lever that is wielded by the substrate-binding domain and the cochaperone onto the ATPase domain to induce a conformation favorable for ATP hydrolysis. Based on our results we propose a mechanism of interdomain communication.

title: An Interdomain Energetic Tug-of-War Creates the Allosterically Active State in Hsp70 Molecular Chaperones
abstract: The allosteric mechanism of Hsp70 molecular chaperones enables ATP binding to the N-terminal nucleotide-binding domain (NBD) to alter substrate affinity to the C-terminal substrate-binding domain (SBD) and substrate binding to enhance ATP hydrolysis. Cycling between ATP-bound and ADP/substrate-bound states requires Hsp70s to visit a state with high ATPase activity and fast on/off kinetics of substrate binding. We have trapped this "allosterically active" state for the E. coli Hsp70, DnaK, and identified how interactions among the NBD, the β subdomain of the SBD, the SBD α-helical lid, and the conserved hydrophobic interdomain linker enable allosteric signal transmission between ligand-binding sites. Allostery in Hsp70s results from an energetic tug-of-war between domain conformations and formation of two orthogonal interfaces: between the NBD and SBD, and between the helical lid and the β subdomain of the SBD. The resulting energetic tension underlies Hsp70 functional properties and enables them to be modulated by ligands and cochaperones and "tuned" through evolution.

title: Structure and energetics of an allele-specific genetic interaction between dnaJ and dnaK: correlation of nuclear magnetic resonance chemical shift perturbations in the J-domain of Hsp40/DnaJ with binding affinity for the ATPase domain of Hsp70/DnaK.
abstract: The molecular chaperone machine composed of Escherichia coli Hsp70/DnaK and Hsp40/DnaJ binds and releases client proteins in cycles of ATP-dependent protein folding, membrane translocation, disassembly, and degradation. The J-domain of DnaJ simultaneously stimulates ATP hydrolysis in the ATPase domain and capture of the client protein in the peptide-binding domain of DnaK. ATP-dependent binding of DnaJ to DnaK mimics DnaJ-dependent capture of a client protein. The dnaJ mutation that replaces aspartate-35 with asparagine (D35N) in the J-domain causes a defect in binding of DnaJ to DnaK. The dnaK mutation that replaces arginine-167 with alanine (R167A) in the ATPase domain of DnaK(R167A) restores binding of DnaJ(D35N). This genetic interaction was said to be allele-specific because wild-type DnaJ does not bind to DnaK(R167A). The J-domain of DnaJ binds to the ATPase domain of DnaK in its capacity as modulator of DnaK ATPase activity and conformational behavior. Surprisingly, the mutations affect the domainwise interaction in an almost opposite manner. D35N increases the affinity of the J-domain for the ATPase domain. R167A has no affect on the affinity of the ATPase domain for the D35N mutant J-domain, but it reduces the affinity for the wild-type J-domain. Previous amide ((1)H, (15)N) NMR chemical shift perturbation mapping in the J-domain suggested that the ATPase domain binds to J-domain helix II and the flanking loops. In the D35N mutant J-domain, chemical shift perturbations include additional effects at amides in the flexible loop II-III and helix III, which have been proposed to undergo an induced fit conformational change upon binding to DnaK. The integrated magnitudes of chemical shift perturbations for the various J-domain and ATPase domain pairs correlate with the free energies of binding. Thus, the J-domain structure can be described as a dynamic ensemble of conformations that is constrained by binding to the ATPase domain. J-domain helix II bends upon binding to the ATPase domain. D35N increases helix II bending, but less so in combination with R167A in the ATPase domain. Taken together, the results suggest that D35N overstabilizes an induced fit conformational change in loop II-III and helix III that is necessary for the J-domain to couple ATP hydrolysis with a conformational change in DnaK, and R167A destabilizes the induced conformation. Conclusions from this work have implications for understanding mechanisms of protein-protein interaction that are involved in allosteric regulation and genetic suppression.

title: Nucleotide-Dependent Interactions within a Specialized Hsp70/Hsp40 Complex Involved in Fe–S Cluster Biogenesis
abstract: The structural mechanism by which Hsp70-type chaperones interact with Hsp40-type co-chaperones has been of great interest, yet still remains a matter of debate. Here, we used solution NMR spectroscopy to investigate the ATP-/ADP-dependent interactions between Escherichia coli HscA and HscB, the specialized Hsp70/Hsp40 molecular chaperones that mediate iron–sulfur cluster transfer. We observed that NMR signals assigned to amino acid residues in the J-domain and its “HPD” motif of HscB broadened severely upon the addition of ATP-bound HscA, but these signals were not similarly broadened by ADP-bound HscA or the isolated nucleotide binding domain of HscA complexed with either ATP or ADP. An HscB variant with an altered HPD motif, HscB(H32A,P33A,D34A), failed to manifest WT-like NMR signal perturbations and also abolished WT-like stimulation of ATP hydrolysis by HscA. In addition, residues 153–171 in the C-terminal region of HscB exhibited NMR signal perturbations upon interaction with HscA, alone or complexed with ADP or ATP. These results demonstrate that the HPD motif in the J-domain of HscB directly interacts with ATP-bound HscA and suggest that a second, less nucleotide-dependent binding site for HscA resides in the C-terminal region of HscB.

title: Hsc70 contacts helix III of the J domain from polyomavirus T antigens: addressing a dilemma in the chaperone hypothesis of how they release E2F from pRb.
abstract: Hsc70's expected binding site on helix II of the J domain of T antigens appears to be blocked in its structure bound to tumor suppressor pRb. We used NMR to map where mammalian Hsc70 binds the J domain of murine polyomavirus T antigens (PyJ). The ATPase domain of Hsc70 unexpectedly has its biggest effects on the NMR peak positions of the C-terminal end of helix III of PyJ. The Hsc70 ATPase domain protects the C-terminal end of helix III of PyJ from an uncharged paramagnetic probe of chelated Gd(III), clearly suggesting the interface. Effects on the conserved HPD loop and helix II of PyJ are smaller. The NMR results are supported by a novel assay of Hsc70's ATP hydrolysis showing that mutations of surface residues in PyJ helix III impair PyJ-dependent stimulation of Hsc70 activity. Evolutionary trace analysis of J domains suggests that helix III usually may join helix II in contributing specificities for cognate hsp70s. Our novel evidence implicating helix III differs from evidence that Escherichia coli DnaK primarily affects helix II and the HPD loop of DnaJ. We find the pRb-binding fragment of E2F1 to be intrinsically unfolded and a good substrate for Hsc70 in vitro. This suggests that E2F1 could be a substrate for Hsc70 recruited by T antigen to an Rb family member. Importantly, our results strengthen the chaperone hypothesis for E2F release from an Rb family member by Hsc70 recruited by large T antigen. That is, it now appears that Hsc70 can freely access helix III and the HPD motif of large T antigen bound to an Rb family member.

title: Peptide Folding: When Simulation Meets Experiment
abstract: Mol. dynamics simulation studies on the folding of beta-peptides H-beta3-HVal-beta3-HAla-beta3-HLeu-(S,S)-beta3-HAla(alphaMe)-beta3-HVal-beta3-HAla-beta3-HLeu-OH and H-beta2-HVal-beta3-HAla-beta2-HLeu-beta3-HVal-beta2-HAla-beta3-HLeu-OH were carried out. Despite the small differences in sequence between the two peptides studied, the simulations correctly predict a left-handed 31-helical fold for the beta-heptapeptide and a right-handed helical fold for the beta-hexapeptide.

title: HMMER web server: interactive sequence similarity searching
abstract: HMMER is a software suite for protein sequence similarity searches using probabilistic methods. Previously, HMMER has mainly been available only as a computationally intensive UNIX command-line tool, restricting its use. Recent advances in the software, HMMER3, have resulted in a 100-fold speed gain relative to previous versions. It is now feasible to make efficient profile hidden Markov model (profile HMM) searches via the web. A HMMER web server (http://hmmer.janelia.org) has been designed and implemented such that most protein database searches return within a few seconds. Methods are available for searching either a single protein sequence, multiple protein sequence alignment or profile HMM against a target sequence database, and for searching a protein sequence against Pfam. The web server is designed to cater to a range of different user expertise and accepts batch uploading of multiple queries at once. All search methods are also available as RESTful web services, thereby allowing them to be readily integrated as remotely executed tasks in locally scripted workflows. We have focused on minimizing search times and the ability to rapidly display tabular results, regardless of the number of matches found, developing graphical summaries of the search results to provide quick, intuitive appraisement of them.

title: Fast pseudolikelihood maximization for direct-coupling analysis of protein structure from many homologous amino-acid sequences
abstract: Direct-coupling analysis is a group of methods to harvest information about coevolving residues in a protein family by learning a generative model in an exponential family from data. In protein families of realistic size, this learning can only be done approximately, and there is a trade-off between inference precision and computational speed. We here show that an earlier introduced l 2 -regularized pseudolikelihood maximization method called plmDCA can be modified as to be easily parallelizable, as well as inherently faster on a single processor, at negligible difference in accuracy. We test the new incarnation of the method on 143 protein family/structure-pairs from the Protein Families database (PFAM), one of the larger tests of this class of algorithms to date.

title: Benchmarking and Analysis of Protein Docking Performance in Rosetta v3.2
abstract: RosettaDock has been increasingly used in protein docking and design strategies in order to predict the structure of protein-protein interfaces. Here we test capabilities of RosettaDock 3.2, part of the newly developed Rosetta v3.2 modeling suite, against Docking Benchmark 3.0, and compare it with RosettaDock v2.3, the latest version of the previous Rosetta software package. The benchmark contains a diverse set of 116 docking targets including 22 antibody-antigen complexes, 33 enzyme-inhibitor complexes, and 60 ‘other’ complexes. These targets were further classified by expected docking difficulty into 84 rigid-body targets, 17 medium targets, and 14 difficult targets. We carried out local docking perturbations for each target, using the unbound structures when available, in both RosettaDock v2.3 and v3.2. Overall the performances of RosettaDock v2.3 and v3.2 were similar. RosettaDock v3.2 achieved 56 docking funnels, compared to 49 in v2.3. A breakdown of docking performance by protein complex type shows that RosettaDock v3.2 achieved docking funnels for 63% of antibody-antigen targets, 62% of enzyme-inhibitor targets, and 35% of ‘other’ targets. In terms of docking difficulty, RosettaDock v3.2 achieved funnels for 58% of rigid-body targets, 30% of medium targets, and 14% of difficult targets. For targets that failed, we carry out additional analyses to identify the cause of failure, which showed that binding-induced backbone conformation changes account for a majority of failures. We also present a bootstrap statistical analysis that quantifies the reliability of the stochastic docking results. Finally, we demonstrate the additional functionality available in RosettaDock v3.2 by incorporating small-molecules and non-protein co-factors in docking of a smaller target set. This study marks the most extensive benchmarking of the RosettaDock module to date and establishes a baseline for future research in protein interface modeling and structure prediction.

title: GROMACS: High performance molecular simulations through multi-level parallelism from laptops to supercomputers
abstract: GROMACS is one of the most widely used open-source and free software codes in chemistry, used primarily for dynamical simulations of biomolecules. It provides a rich set of calculation types, prepa ...

title: Comparison of simple potential functions for simulating liquid water
abstract: Classical Monte Carlo simulations have been carried out for liquid water in the NPT ensemble at 25 °C and 1 atm using six of the simpler intermolecular potential functions for the water dimer: Bernal–Fowler (BF), SPC, ST2, TIPS2, TIP3P, and TIP4P. Comparisons are made with experimental thermodynamic and structural data including the recent neutron diffraction results of Thiessen and Narten. The computed densities and potential energies are in reasonable accord with experiment except for the original BF model, which yields an 18% overestimate of the density and poor structural results.  TIPS2 and TIP4P potentials yield oxygen–oxygen partial structure functions in good agreement with the neutron diffraction results. The accord with the experimental OH and HH partial structure functions is poorer; however, the computed results for these functions are similar for all the potential functions. Consequently, the discrepancy may be due to the correction terms needed in processing the neutron data or to an effect uniformly neglected in the computations. Comparisons are also made for self‐diffusion coefficients obtained from molecular dynamics simulations. Overall, the SPC, ST2, TIPS2, and TIP4P models give reasonable structural and thermodynamic descriptions of liquid water and they should be useful in simulations of aqueous solutions. The simplicity of the SPC, TIPS2, and TIP4P functions is also attractive from a computational standpoint.

title: Canonical sampling through velocity rescaling.
abstract: The authors present a new molecular dynamics algorithm for sampling the canonical distribution. In this approach the velocities of all the particles are rescaled by a properly chosen random factor. The algorithm is formally justified and it is shown that, in spite of its stochastic nature, a quantity can still be defined that remains constant during the evolution. In numerical applications this quantity can be used to measure the accuracy of the sampling. The authors illustrate the properties of this new method on Lennard-Jones and TIP4P water models in the solid and liquid phases. Its performance is excellent and largely independent of the thermostat parameter also with regard to the dynamic properties.

title: LINCS: A linear constraint solver for molecular simulations
abstract: In this article, we present a new LINear Constraint Solver (LINCS) for molecular simulations with bond constraints. The algorithm is inherently stable, as the constraints themselves are reset instead of derivatives of the constraints, thereby eliminating drift. Although the derivation of the algorithm is presented in terms of matrices, no matrix matrix multiplications are needed and only the nonzero matrix elements have to be stored, making the method useful for very large molecules. At the same accuracy, the LINCS algorithm is three to four times faster than the SHAKE algorithm. Parallelization of the algorithm is straightforward.