Instantaneous Normal Modes as an Unforced Reaction Coordinate for Protein Conformational Transitions

We present a novel sampling approach to explore large protein conformational transitions by determining unique substates from instantaneous normal modes calculated from an elastic network model, and applied to a progression of atomistic molecular dynamics snapshots. This unbiased sampling scheme allows us to direct the path sampling between the conformational end states over simulation timescales that are greatly reduced relative to the known experimental timescales. We use adenylate kinase as a test system to show that instantaneous normal modes can be used to identify substates that drive the structural fluctuations of adenylate kinase from its closed to open conformations, in which we observe 16 complete transitions in 4 μs of simulation time, reducing the timescale over conventional simulation timescales by two orders of magnitude. Analysis shows that the unbiased determination of substates is consistent with known pathways determined experimentally.     

Enhanced resolution of molecular recognition to distinguish structurally similar molecules by different conformational responses of a protein upon ligand binding

MutT distinguishes substrate 8-oxo-dGTP from dGTP and also 8-oxo-dGMP from dGMP despite small differences of chemical structures between them. In this paper we show by the method of molecular dynamics simulation that the transition between conformational substates of MutT is a key mechanism for a high-resolution molecular recognition of the differences between the very similar chemical compounds. (1) The native state MutT has two conformational substates with similar free energies, each characterized by either open or closed of two loops surrounding the substrate binding active site. Between the two substates, the open substate is more stable in free MutT and in dGMP-MutT complex, and the closed substate is more stable in 8-oxo-dGMP-MutT complex. (2) Conformational fluctuation of the open substate is much larger than that of the closed substate. An estimate of associated entropy difference was found to be consistent with the experimentally found difference of entropy contribution to the binding free energies of the two molecules. (3) A hydrogen bond between H7 atom of 8-oxo-dGMP and the sidechain of Asn119 plays a crucial role for maintaining the closed substate in 8-oxo-dGMP-MutT complex. When this hydrogen bond is absent in the H7-deficient dGMP-MutT complex, the closed substate is no more maintained and transition to the more entropically-favored open substate is induced. (4) Thus, this mechanism of the hydrogen bond controlling the relative stabilities of the drastically different two conformational substates enhances the resolution to recognize the small difference of the chemical structures between the two molecules, dGMP and 8-oxo-dGMP.     

Mapping the transition state for DNA bending by IHF

How DNA-bending proteins recognize their specific sites on DNA remains elusive, particularly for proteins that use indirect readout, which relies on sequence-dependent variations in DNA flexibility/bendability. The question remains as to whether the protein bends the DNA (protein-induced bending) or, alternatively, "prebent" DNA conformations are thermally accessible, which the protein captures to form the specific complex (conformational capture). To distinguish between these mechanisms requires characterization of reaction intermediates and, in particular, snapshots of the transition state along the recognition pathway. We present such a snapshot, from measurements of DNA bending dynamics in complex with Escherichia coli integration host factor (IHF), an architectural protein that bends specific sites on λ-DNA in a U-turn by creating two sharp kinks in DNA. Fluorescence resonance energy transfer measurements in response to laser temperature-jump perturbation monitor DNA bending. We find that nicks or mismatches that enhance DNA flexibility at the site of the kinks show 3- to 4-fold increase in DNA bending rates that reflect a 4- to 11-fold increase in binding affinities, while sequence modifications away from the kink sites, as well as mutations in IHF designed to destabilize the complex, have negligible effect on DNA bending rates despite >250-fold decrease in binding affinities. These results support the scenario that the bottleneck in the recognition step for IHF is spontaneous kinking of cognate DNA to adopt a partially prebent conformation and point to conformational capture as the underlying mechanism of initial recognition, with additional protein-induced bending occurring after the transition state.     

Role of methyl groups in dynamics and evolution of biomolecules

Recent studies have discovered strong differences between the dynamics of nucleic acids (RNA and DNA) and proteins, especially at low hydration and low temperatures. This difference is caused primarily by dynamics of methyl groups that are abundant in proteins, but are absent or very rare in RNA and DNA. In this paper, we present a hypothesis regarding the role of methyl groups as intrinsic plasticizers in proteins and their evolutionary selection to facilitate protein dynamics and activity. We demonstrate the profound effect methyl groups have on protein dynamics relative to nucleic acid dynamics, and note the apparent correlation of methyl group content in protein classes and their need for molecular flexibility. Moreover, we note the fastest methyl groups of some enzymes appear around dynamical centers such as hinges or active sites. Methyl groups are also of tremendous importance from a hydrophobicity/folding/entropy perspective. These significant roles, however, complement our hypothesis rather than preclude the recognition of methyl groups in the dynamics and evolution of biomolecules.     

Indirect DNA Readout on the Protein Side: Coupling between Histidine Protonation, Global Structural Cooperativity, Dynamics, and DNA Binding of the Human Papillomavirus Type 16 E2C Domain

DNA sequence recognition by the homodimeric C-terminal domain of the human papillomavirus type 16 E2 protein (E2C) is known to involve both direct readout and DNA-dependent indirect readout mechanisms, while protein-dependent indirect readout has been deduced but not directly observed. We have investigated coupling between specific DNA binding and the dynamics of the unusual E2C fold, using pH as an external variable. Nuclear magnetic resonance and isothermal titration calorimetry show that pH titration of His318 in the complex interface and His288 in the core of the domain is coupled to both binding and the dynamics of the β-barrel core of E2C, with a tradeoff between dimer stability and function. Specific DNA binding is, in turn, coupled to the slow dynamics and amide hydrogen exchange in the entire β-barrel, reaching residues far apart from the DNA recognition elements but not affecting the two helices of each monomer. The changes are largest in the dimerization interface, suggesting that the E2C β-barrel acts as a hinge that regulates the relative position of the DNA recognition helices. In conclusion, the cooperative dynamics of the human papillomavirus type 16 E2C β-barrel is coupled to sequence recognition in a protein-dependent indirect readout mechanism. The patterns of residue substitution in genital papillomaviruses support the importance of the protonation states of His288 and His318 and suggest that protein-dependent indirect readout and histidine pH titration may regulate DNA binding in the cell.     

A simple model of backbone flexibility improves modeling of side-chain conformational variability

The considerable flexibility of side-chains in folded proteins is important for protein stability and function, and may have a role in mediating allosteric interactions. While sampling side-chain degrees of freedom has been an integral part of several successful computational protein design methods, the predictions of these approaches have not been directly compared to experimental measurements of side-chain motional amplitudes. In addition, protein design methods frequently keep the backbone fixed, an approximation that may substantially limit the ability to accurately model side-chain flexibility. Here, we describe a Monte Carlo approach to modeling side-chain conformational variability and validate our method against a large dataset of methyl relaxation order parameters derived from nuclear magnetic resonance (NMR) experiments (17 proteins and a total of 530 data points). We also evaluate a model of backbone flexibility based on Backrub motions, a type of conformational change frequently observed in ultra-high-resolution X-ray structures that accounts for correlated side-chain backbone movements. The fixed-backbone model performs reasonably well with an overall rmsd between computed and predicted side-chain order parameters of 0.26. Notably, including backbone flexibility leads to significant improvements in modeling side-chain order parameters for ten of the 17 proteins in the set. Greater accuracy of the flexible backbone model results from both increases and decreases in side-chain flexibility relative to the fixed-backbone model. This simple flexible-backbone model should be useful for a variety of protein design applications, including improved modeling of protein-protein interactions, design of proteins with desired flexibility or rigidity, and prediction of correlated motions within proteins.     

Influence of 8-Oxoguanosine on the Fine Structure of DNA Studied with Biasing-Potential Replica Exchange Simulations

Chemical modification or radiation can cause DNA damage, which plays a crucial role for mutagenesis of DNA, carcinogenesis, and aging. DNA damage can also alter the fine structure of DNA that may serve as a recognition signal for DNA repair enzymes. A new, advanced sampling replica-exchange method has been developed to specifically enhance the sampling of conformational substates in duplex DNA during molecular dynamics (MD) simulations. The approach employs specific biasing potentials acting on pairs of pseudodihedral angles of the nucleic acid backbone that are added in the replica simulations to promote transitions of the most common substates of the DNA backbone. The sampled states can exchange with a reference simulation under the control of the original force field. The application to 7,8-dihydro-8oxo-guanosine, one of the most common oxidative damage in DNA indicated better convergence of sampled states during 10 ns simulations compared to 20 times longer standard MD simulations. It is well suited to study systematically the fine structure and dynamics of large nucleic acids under realistic conditions, including explicit solvent and ions. The biasing potential-replica exchange MD simulations indicated significant differences in the population of nucleic acid backbone substates in the case of 7,8-dihydro-8oxo-guanosine compared to a regular guanosine in the same sequence context. This concerns both the ratio of the B-DNA substates B_I and B_II associated with the backbone dihedral angles ε and ζ but also coupled changes in the backbone dihedral angles α and γ. Such differences may play a crucial role in the initial recognition of damaged DNA by repair enzymes.     

Solvent participation in Serratia marcescens endonuclease complexes

The monomer and dimer of the bacterium Serratia marcescens endonuclease (SMnase) are each catalytically active and the two subunits of the dimer function independently of each other. Specific interfacial waters may play a role in stability, complex formation, and functionality. We performed molecular dynamics simulations of both a subunit of SMnase and its model built complex with DNA and analyzed the relation of the hydration sites to the catalytic mechanism. It was found that the binding of DNA has little influence on the global hydration properties of the protein, including occupancy and water residence time distributions. DNA and protein recognition in our model mainly involves direct contacts by hydrogen bond and hydrophobic interactions. Water-mediated contacts exist, but are less common. Three interior water clusters were identified for SMnase. One cluster around the active site in the monomer-DNA complex shows relatively strong interactions between hydration sites as well as between the sites and the biomolecules. The simulated cluster properties agreed well with experimental data. The magnesium ion shows ligand exchange. Although Mg2+ keeps six ligands during the entire simulation, upon the binding of DNA, Asn119 loses its coordination with Mg2+, while one nonbridging oxygen of the phosphate of a DNA residue and two oxygen atoms of the side chain of Glu127 become the ligands. Waters in a nearby cluster exchange and participate in the resolvation of groups in the presence of DNA. Water thus not only participates in the cleavage of DNA but also can stabilize the transition state and the leaving groups in our model.     

PABMB Lecture. Protein dynamics, folding and misfolding: from basic physical chemistry to human conformational diseases.

Proteins exhibit a variety of motions ranging from amino acid side-chain rotations to the motions of large domains. Recognition of their conformational flexibility has led to the view that protein molecules undergo fast dynamic interconversion between different conformational substates. This proposal has received support from a wide variety of experimental techniques and from computer simulations of protein dynamics. More recently, studies of the subunit dissociation of oligomeric proteins induced by hydrostatic pressure have shown that the characteristic times for subunit exchange between oligomers and for interconversion between different conformations may be rather slow (hours or days). In such cases, proteins cannot be treated as an ensemble of rapidly interconverting conformational substates, but rather as a persistently heterogeneous population of different long-lived conformers. This is reminiscent of the deterministic behavior exhibited by macroscopic bodies, and may have important implications for our understanding of protein folding and biological functions. Here, we propose that the deterministic behavior of proteins may be closely related to the genesis of conformational diseases, a class of pathological conditions that includes transmissible spongiform encephalopathies, Alzheimer's disease and other amyloidosis.     

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.     

The impact of protein flexibility on protein-protein docking

Accounting for protein flexibility in protein-protein docking algorithms is challenging, and most algorithms therefore treat proteins as rigid bodies or permit side-chain motion only. While the consequences are obvious when there are large conformational changes upon binding, the situation is less clear for the modest conformational changes that occur upon formation of most protein-protein complexes. We have therefore studied the impact of local protein flexibility on protein-protein association by means of rigid body and torsion angle dynamics simulation. The binding of barnase and barstar was chosen as a model system for this study, because the complexation of these 2 proteins is well-characterized experimentally, and the conformational changes accompanying binding are modest. On the side-chain level, we show that the orientation of particular residues at the interface (so-called hotspot residues) have a crucial influence on the way contacts are established during docking from short protein separations of approximately 5 A. However, side-chain torsion angle dynamics simulations did not result in satisfactory docking of the proteins when using the unbound protein structures. This can be explained by our observations that, on the backbone level, even small (2 A) local loop deformations affect the dynamics of contact formation upon docking. Complementary shape-based docking calculations confirm this result, which indicates that both side-chain and backbone levels of flexibility influence short-range protein-protein association and should be treated simultaneously for atomic-detail computational docking of proteins.     

Structural and dynamic determinants of protein-peptide recognition

Protein-peptide interactions play important roles in many cellular processes, including signal transduction, trafficking, and immune recognition. Protein conformational changes upon binding, an ill-defined peptide binding surface, and the large number of peptide degrees of freedom make the prediction of protein-peptide interactions particularly challenging. To address these challenges, we perform rapid molecular dynamics simulations in order to examine the energetic and dynamic aspects of protein-peptide binding. We find that, in most cases, we recapitulate the native binding sites and native-like poses of protein-peptide complexes. Inclusion of electrostatic interactions in simulations significantly improves the prediction accuracy. Our results also highlight the importance of protein conformational flexibility, especially side-chain movement, which allows the peptide to optimize its conformation. Our findings not only demonstrate the importance of sufficient sampling of the protein and peptide conformations, but also reveal the possible effects of electrostatics and conformational flexibility on peptide recognition.     

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.     

Side-chain structure and dynamics at the lipid-protein interface: Val1 of the gramicidin A channel.

High resolution dynamics and structural information has been resolved from 2H solid-state NMR spectra of the Val-1 side-chain of the gramicidin channel in a lipid bilayer. Both powder pattern lineshapes and spectra from uniformly aligned samples of gramicidin in lipid bilayers have been analyzed to achieve a fully consistant interpretation of the data. Torsional motions about the C alpha C beta axis (chi 1) are shown to be three-state jumps in which the occupancy of the states is given by the ratio, 75:15:10 for the chi 1 angles of 184 degrees:304 degrees:64 degrees. The dominant conformer is also the most common conformation observed for valines in well defined protein structures. The distribution of conformational substates that represents the chi 1 dynamics appears to be largely independent of the lipid phase transition and the hydration of the sample. However, there is evidence that the residence time between jumps is dependent on the lipid phase transition. Although this time is shown to be approximately 1 microseconds below the phase transition temperature, it is in the fast exchange limit above the transition temperature.     

The B- to A-DNA transition and the reorganization of solvent at the DNA surface

DNA geometry depends on relative humidity. Using the CHARMM22 force field to push B-DNA to A-DNA, a molecular dynamics simulation of a mixed-sequence 24-basepair DNA double-stranded oligomer, starting from B-DNA, was carried out to explore both the mechanism of the transition and the evolution of hydration patterns on the surface of DNA. Over the 11-ns trajectory, the transition recapitulates the slide-first, roll-later mechanism, is opposed by DNA electrostatics, and is favored by an increasing amount of condensed sodium ions. Hydration was characterized by counting the hydrogen bonds between water and DNA, and by the number of water bridges linking two DNA atoms. The number of hydrogen bonds between water and DNA remains constant during the transition, but there is a 40% increase in the number of water bridges, in agreement with the principle of economy of hydration. Water bridges emerge as delicate sensors of both structure and dynamics of DNA. Both local flexibility and the frustration of the water network on the surface of DNA probably account for the low populations and short residence times of the bridges, and for the lubricant role of water in ligand-DNA interactions.     

Protein flexibility directs DNA recognition by the papillomavirus E2 proteins

Although DNA flexibility is known to play an important role in DNA-protein interactions, the importance of protein flexibility is less well understood. Here, we show that protein dynamics are important in DNA recognition using the well-characterized human papillomavirus (HPV) type 6 E2 protein as a model system. We have compared the DNA binding properties of the HPV 6 E2 DNA binding domain (DBD) and a mutant lacking two C-terminal leucine residues that form part of the hydrophobic core of the protein. Deletion of these residues results in increased specific and non-specific DNA binding and an overall decrease in DNA binding specificity. Using (15)N NMR relaxation and hydrogen/deuterium exchange, we demonstrate that the mutation results in increased flexibility within the hydrophobic core and loop regions that orient the DNA binding helices. Stopped-flow kinetic studies indicate that increased flexibility alters DNA binding by increasing initial interactions with DNA but has little or no effect on the structural rearrangements that follow this step. Taken together these data demonstrate that subtle changes in protein dynamics have a major influence on protein-DNA interactions.     

Dynamics and consequences of DNA looping by the FokI restriction endonuclease

Genetic events often require proteins to be activated by interacting with two DNA sites, trapping the intervening DNA in a loop. While much is known about looping equilibria, only a few studies have examined DNA-looping dynamics experimentally. The restriction enzymes that cut DNA after interacting with two recognition sites, such as FokI, can be used to exemplify looping reactions. The reaction pathway for FokI on a supercoiled DNA with two sites was dissected by fast kinetics to reveal, in turn: the initial binding of a protein monomer to each site; the protein–protein association to form the dimer, trapping the loop; the subsequent phosphodiester hydrolysis step. The DNA motion that juxtaposes the sites ought on the basis of Brownian dynamics to take ~2 ms, but loop capture by FokI took 230 ms. Hence, DNA looping by FokI is rate limited by protein association rather than DNA dynamics. The FokI endonuclease also illustrated activation by looping: it cut looped DNA 400 times faster than unlooped DNA.