Functionally important correlated motions in the single-stranded DNA-binding protein encoded by filamentous phage Pf3

The Dynamic Behavior of the P2X4 Ion Channel in the Closed Conformation

We present the results of a detailed molecular dynamics study of the closed form of the P2X₄ receptor. The fluctuations observed in the simulations were compared with the changes that occur in the transition from the closed to the open structure. To get further insight on the opening mechanism, the actual displacements were decomposed into interchain motions and intrachain deformations. This analysis revealed that the iris-like expansion of the transmembrane helices mainly results from interchain motions that already take place in the closed conformation. However, these movements cannot reach the amplitude required for the opening of the channel because they are impeded by interactions occurring around the ATP binding pocket. This suggests that the union of ATP produces distortions in the chains that eliminate the restrictions on the interchain displacements, leading to the opening of the pore.

Structure and assembly of filamentous bacteriophages

Filamentous bacteriophages are interesting paradigms in structural molecular biology, in part because of the unusual mechanism of filamentous phage assembly. During assembly, several thousand copies of an intracellular DNA-binding protein bind to each copy of the replicating phage DNA, and are then displaced by membrane-spanning phage coat proteins as the nascent phage is extruded through the bacterial plasma membrane. This complicated process takes place without killing the host bacterium. The bacteriophage is a semi-flexible worm-like nucleoprotein filament. The virion comprises a tube of several thousand identical major coat protein subunits around a core of single-stranded circular DNA. Each protein subunit is a polymer of about 50 amino-acid residues, largely arranged in an α-helix. The subunits assemble into a helical sheath, with each subunit oriented at a small angle to the virion axis and interdigitated with neighbouring subunits. A few copies of "minor" phage proteins necessary for infection and/or extrusion of the virion are located at each end of the completed virion. Here we review both the structure of the virion and aspects of its function, such as the way the virion enters the host, multiplies, and exits to prey on further hosts. In particular we focus on our understanding of the way the components of the virion come together during assembly at the membrane. We try to follow a basic rule of empirical science, that one should chose the simplest theoretical explanation for experiments, but be prepared to modify or even abandon this explanation as new experiments add more detail.

New insights in the activation of human cholesterol esterase to design potent anti-cholesterol drugs

Primary hypercholesterolemia is the root cause for major health issues like coronary heart disease and atherosclerosis. Regulating plasma cholesterol level, which is the product of biosynthesis as well as dietary intake, has become one of the major therapeutic strategies to effectively control these diseases. Human cholesterol esterase (hCEase) is an interesting target involved in the regulation of plasma cholesterol level and thus inhibition of this enzyme is highly effective in the treatment of hypercholesterolemia. This study was designed to understand the activation mechanism that enables the enzyme to accommodate long chain fatty acids and to identify the structural elements for the successful catalysis. Primarily the activation efficiencies of three different bile salts were studied and compared using molecular dynamics simulations. Based on the conformations of major surface loops, hydrogen bond interactions, and distance analyses, taurocholate was concluded as the preferred activator of the enzyme. Furthermore, the importance of two bile salt binding sites (proximal and remote) and the crucial role of 7α-OH group of the bile salts in the activation of hCEase was examined and evidenced. The results of our study explain the structural insights of the activation mechanism and show the key features of the bile salts responsible for the enzyme activation which are very useful in hypolipidemic drug designing strategies.

Exploring residue component contributions to dynamical network models of allostery

Allosteric regulation in biological systems is of considerable interest given the vast number of proteins that exhibit such behavior. Network models obtained from molecular dynamics simulations have been shown to be powerful tools for the analysis of allostery. In this work, different coarse-grain residue representations (nodes) are used together with a dynamical network model to investigate models of allosteric regulation. This model assumes that allosteric signals are dependent on positional correlations of protein substituents, as determined through molecular dynamics simulations, and uses correlated motion to generate a signaling weight between two given nodes. We examine four types of network models using different node representations in Cartesian coordinates: the (i) residue alpha-carbons, (ii) sidechain center of mass, (iii) backbone center of mass, and the entire (iv) residue center of mass. All correlations are filtered by a dynamic contact map that defines the allowable interactions between nodes based on physical proximity. We apply the four models to imidazole glycerol phosphate synthase (IGPS), which provides a well-studied experimental framework in which allosteric communication is known to persist across disparate protein domains (e.g. a protein dimer interface). IGPS is modeled as a network of nodes and weighted edges. Optimal allosteric pathways are traced using the Floyd Warshall algorithm for weighted networks, and community analysis (a form of hierarchical clustering) is performed using the Girvan-Newman algorithm. Our results show that dynamical information encoded in the residue center of mass must be included in order to detect residues that are experimentally known to play a role in allosteric communication for IGPS. More broadly, this new method may be useful for predicting pathways of allosteric communication for any biomolecular system in atomic detail.

The molecular origin of the MMR-dependent apoptosis pathway from dynamics analysis of MutSα-DNA complexes

The cellular response to DNA damage signaling by mismatch-repair (MMR) proteins is incompletely understood. It is generally accepted that MMR-dependent apoptosis pathway in response to DNA damage detection is independent of MMR's DNA repair function. In this study, we investigate correlated motions in response to the binding of mismatched and platinum cross-linked DNA fragments by MutSα, as derived from 50 ns molecular dynamics simulations. The protein dynamics in response to the mismatched and damaged DNA recognition suggests that MutSα signals their recognition through independent pathways providing evidence for the molecular origin of the MMR-dependent apoptosis. MSH2 subunit is indicated to play a key role in signaling both mismatched and damaged DNA recognition; localized and collective motions within the protein allow identifying sites on the MSH2 surface possible involved in recruiting proteins responsible for downstream events. Unlike in the mismatch complex, predicted key communication sites specific for the damage recognition are on the list of known cancer-causing mutations or deletions. This confirms MSH2's role in signaling DNA damage-induced apoptosis and suggests that defects in MMR alone is sufficient to trigger tumorigenesis, supporting the experimental evidence that MMR-damage response function could protect from the early occurrence of tumors. Identifying these particular communication sites may have implications for the treatment of cancers that are not defective for MMR, but are unable to function optimally for MMR-dependent responses following DNA damage such as the case of resistance to cisplatin.

Quantifying Correlations Between Allosteric Sites in Thermodynamic Ensembles

Allostery describes altered protein function at one site due to a perturbation at another site. One mechanism of allostery involves correlated motions, which can occur even in the absence of substantial conformational change. We present a novel method, "MutInf", to identify statistically significant correlated motions from equilibrium molecular dynamics simulations. Our approach analyzes both backbone and sidechain motions using internal coordinates to account for the gear-like twists that can take place even in the absence of the large conformational changes typical of traditional allosteric proteins. We quantify correlated motions using a mutual information metric, which we extend to incorporate data from multiple short simulations and to filter out correlations that are not statistically significant. Applying our approach to uncover mechanisms of cooperative small molecule binding in human interleukin-2, we identify clusters of correlated residues from 50 ns of molecular dynamics simulations. Interestingly, two of the clusters with the strongest correlations highlight known cooperative small-molecule binding sites and show substantial correlations between these sites. These cooperative binding sites on interleukin-2 are correlated not only through the hydrophobic core of the protein but also through a dynamic polar network of hydrogen bonding and electrostatic interactions. Since this approach identifies correlated conformations in an unbiased, statistically robust manner, it should be a useful tool for finding novel or "orphan" allosteric sites in proteins of biological and therapeutic importance.

Solution structure of the DNA-binding domain of RPA from Saccharomyces cerevisiae and its interaction with single-stranded DNA and SV40 T antigen

Replication protein A (RPA) is a three-subunit complex with multiple roles in DNA metabolism. DNA-binding domain A in the large subunit of human RPA (hRPA70A) binds to single-stranded DNA (ssDNA) and is responsible for the species-specific RPA–T antigen (T-ag) interaction required for Simian virus 40 replication. Although Saccharomyces cerevisiae RPA70A (scRPA70A) shares high sequence homology with hRPA70A, the two are not functionally equivalent. To elucidate the similarities and differences between these two homologous proteins, we determined the solution structure of scRPA70A, which closely resembled the structure of hRPA70A. The structure of ssDNA-bound scRPA70A, as simulated by residual dipolar coupling-based homology modeling, suggested that the positioning of the ssDNA is the same for scRPA70A and hRPA70A, although the conformational changes that occur in the two proteins upon ssDNA binding are not identical. NMR titrations of hRPA70A with T-ag showed that the T-ag binding surface is separate from the ssDNA-binding region and is more neutral than the corresponding part of scRPA70A. These differences might account for the species-specific nature of the hRPA70A–T-ag interaction. Our results provide insight into how these two homologous RPA proteins can exhibit functional differences, but still both retain their ability to bind ssDNA.

Preorganization and protein dynamics in enzyme catalysis

Recently, an alternative has been offered to the concept of transition state (TS) stabilization as an explanation for rate enhancements in enzyme-catalyzed reactions. Instead, most of the rate increase has been ascribed to preorganization of the enzyme active site to bind substrates in a geometry close to that of the TS, which then transit the activation barrier impelled by motions along the reaction coordinate. The question as to how an enzyme achieves such preorganization and concomitant TS stabilization as well as potential coupled motions along the reaction coordinate leads directly to the role of protein dynamic motion. Dihydrofolate reductase (DHFR) is a paradigm in which the role of dynamics in catalysis continues to be unraveled by a wealth of kinetic, structural, and computational studies. DHFR has flexible loop regions adjacent to the active site whose motions modulate passage through the kinetically preferred pathway. The participation of residues distant from the DHFR active site in enhancing the rate of hydride transfer, however, is unanticipated and may signify the importance of long range protein motions. The general significance of protein dynamics in understanding other biological processes is briefly discussed.

Characterization of binding-induced changes in dynamics suggests a model for sequence-nonspecific binding of ssDNA by replication protein A

Single-stranded-DNA-binding proteins (SSBs) are required for numerous genetic processes ranging from DNA synthesis to the repair of DNA damage, each of which requires binding with high affinity to ssDNA of variable base composition. To gain insight into the mechanism of sequence-nonspecific binding of ssDNA, NMR chemical shift and (15)N relaxation experiments were performed on an isolated ssDNA-binding domain (RPA70A) from the human SSB replication protein A. The backbone (13)C, (15)N, and (1)H resonances of RPA70A were assigned for the free protein and the d-CTTCA complex. The binding-induced changes in backbone chemical shifts were used to map out the ssDNA-binding site. Comparison to results obtained for the complex with d-C(5) showed that the basic mode of binding is independent of the ssDNA sequence, but that there are differences in the binding surfaces. Amide nitrogen relaxation rates (R(1) and R(2)) and (1)H-(15)N NOE values were measured for RPA70A in the absence and presence of d-CTTCA. Analysis of the data using the Model-Free formalism and spectral density mapping approaches showed that the structural changes in the binding site are accompanied by some significant changes in flexibility of the primary DNA-binding loops on multiple timescales. On the basis of these results and comparisons to related proteins, we propose that the mechanism of sequence-nonspecific binding of ssDNA involves dynamic remodeling of the binding surface.

Independent tyrosyl contributions to the CD of Ff gene 5 protein and the distinctive effects of Y41H and Y41F mutants on protein-protein cooperative interactions

The gene 5 protein (g5p) of the Ff virus contains five Tyr, individual mutants of which have now all been characterized by CD spectroscopy. The protein has a dominant tyrosyl 229-nm L(a) CD band that is shown to be approximately the sum of the five individual Tyr contributions. Tyr41 is particularly important in contributing to the high cooperativity with which the g5p binds to ssDNA, and Y41F and Y41H mutants are known to differ in dimer-dimer packing interactions in crystal structures. We compared the solution structures and binding properties of the Y41F and Y41H mutants using CD spectroscopy. Secondary structures of the mutants were similar by CD analyses and close to those derived from the crystal structures. However, there were significant differences in the binding properties of the two mutant proteins. The Y41H protein had an especially low binding affinity and perturbed the spectrum of poly[d(A)] in 2 mM Na(+) much less than did Y41F and the wild-type gene 5 proteins. Moreover, a change in the Tyr 229 nm band, assigned to the perturbation of Tyr34 at the dimer-dimer interface, was absent in titrations with the Y41H mutant under low salt conditions. In contrast, titrations with the Y41H mutant in 50 mM Na(+) exhibited typical CD changes of both the nucleic acid and the Tyr 229-nm band. Thus, protein-protein and g5p-ssDNA interactions appeared to be mutually influenced by ionic strength, indicative of correlated changes in the ssDNA binding and cooperativity loops of the protein or of indirect structural constraints.

Reorientational eigenmode dynamics: a combined MD/NMR relaxation analysis method for flexible parts in globular proteins

An approach is presented for the interpretation of heteronuclear NMR spin relaxation data in mobile protein parts in terms of reorientational eigenmode dynamics. The method is based on the covariance matrix of the spatial functions of the nuclear spin interactions that cause relaxation expressed as spherical harmonics of rank 2. The approach was applied to characterize the dynamics of a loop region of ubiquitin. The covariance matrix was determined from a conformational ensemble generated by a 5 ns molecular dynamics simulation. It was found that the time correlation functions of the dominant eigenmodes decay in good approximation with a single correlation time. From the reorientational eigenmodes, their eigenvalues, and correlation times, NMR relaxation data were calculated in accordance with Bloch-Wangsness-Redfield relaxation theory and directly compared with experimental (15)N relaxation parameters. Using a fitting procedure, agreement between calculated and experimental data was improved significantly by adjusting eigenvalues and correlation times of the dominant modes. The presented procedure provides detailed information on correlated reorientational dynamics of flexible parts in globular proteins. The covariance matrix was linked to the covariance matrix of backbone dihedral angle fluctuations, allowing one to study the motional behavior of these degrees of freedom on nano- and subnanosecond time scales.

Concerted motions in copper plastocyanin and azurin: an essential dynamics study

Essential dynamics analysis of molecular dynamics simulation trajectories (1.1 ns) of two copper containing electron transfer proteins, plastocyanin and azurin, has been performed. The protein essential modes have been analysed in order to identify large concerted motions which could be relevant for the electron transfer function exerted by these proteins. The analysis, conducted for temporal windows of different lengths along the protein trajectories, shows a rapid convergence and indicates that for both the proteins the predominant internal motions occur in a subspace of only a few degrees of freedom. Moreover, it is found that for both the proteins the likely binding sites (i.e. the hydrophobic and negative patches) with the reaction partners move in a concerted fashion with a few structural regions far from the active site. Such results are discussed in connection with the possible involvement of large concerted motions in the recognition and binding interaction with physiological electron transfer partners.

A space-time structure determination of human CD2 reveals the CD58-binding mode

We describe a procedure for a space-time description of protein structures. The method is capable of determining populations of conformational substates, and amplitudes and directions of internal protein motions. This is achieved by fitting static and dynamic NMR data. The approach is based on the jumping-among-minima concept. First, a wide conformational space compatible with structural NMR data is sampled to find a large set of substates. Subsequently, intrasubstate motions are sampled by using molecular dynamics calculations with force field energy terms. Next, the populations of substates are fitted to NMR relaxation data. By diagonalizing a second moment matrix, directions and amplitudes of motions are identified. The method was applied to the adhesion domain of human CD2. We found that very few substates can account for most of the experimental data. Furthermore, only two types of collective motions have high amplitudes. They represent transitions between a concave (closed) and flat (open) binding face and resemble the change upon counter-receptor (CD58) binding.