Many Local Motions Cooperate to Produce the Adenylate Kinase Conformational Transition

Conformational transitions are functionally important in many proteins. In the enzyme adenylate kinase (AK), two small domains (LID and NMP) close over the larger CORE domain; the reverse (opening) motion limits the rate of catalytic turnover. Here, using double-well Gō simulations of Escherichia coli AK, we elaborate on previous investigations of the AK transition mechanism by characterizing the contributions of rigid-body (Cartesian), backbone dihedral, and contact motions to transition-state (TS) properties. In addition, we compare an apo simulation to a pseudo-ligand-bound simulation to reveal insights into allostery. In Cartesian space, LID closure precedes NMP closure in the bound simulation, consistent with prior coarse-grained models of the AK transition. However, NMP-first closure is preferred in the apo simulation. In backbone dihedral space, we find that, as expected, backbone fluctuations are reduced in the O/C transition in parts of all three domains. Among these “quenching” residues, most in the CORE, especially residues 11–13, are rigidified in the TS of the bound simulation, while residues 42–44 in the NMP are flexible in the TS. In contact space, in both apo and bound simulations, one nucleus of closed-state contacts includes parts of the NMP and CORE; CORE–LID contacts are absent in the TS of the apo simulation but formed in the TS of the bound simulation. From these results, we predict mutations that will perturb the opening and/or closing transition rates by changing the entropy of dihedrals and/or the enthalpy of contacts. Furthermore, regarding allostery, the fully closed structure is populated in the apo simulation, but our contact results imply that ligand binding shifts the preferred O/C transition pathway, thus precluding a simple conformational selection mechanism. Finally, the analytical approach and the insights derived from this work may inform the rational design of flexibility and allostery in proteins.     

Crystal structure of human adenylate kinase 4 (L171P) suggests the role of hinge region in protein domain motion

It is well known that motion of LID and NMP-binding (NMP_bind) domains in adenylate kinase (AK) is important in ligand binding and catalysis. However, the nature of such domain motions is poorly characterized. One of the critical hinge regions is hinge IV, which connects the CORE and LID domains. In addition, the hinge IV contains a strictly conserved residue, L171, in the AK family. To investigate the role of hinge IV, crystal structure of human adenylate kinase 4 (AK4) L171P mutant was determined. This mutation dramatically changes the orientation of the LID domain, which could be described as a novel twisted-and-closed conformation in contrast to the open and closed conformations in other AKs. This mutant provides a new example of domain motions in AK family.     

Crystal structure of the zinc-, cobalt-, and iron-containing adenylate kinase from <Emphasis Type="Italic">Desulfovibrio gigas</Emphasis>: a novel metal-containing adenylate kinase from Gram-negative bacteria

Adenylate kinases (AK) from Gram-negative bacteria are generally devoid of metal ions in their LID domain. However, three metal ions, zinc, cobalt, and iron, have been found in AK from Gram-negative bacteria. Crystal structures of substrate-free AK from Desulfovibrio gigas with three different metal ions (Zn²⁺, Zn-AK; Co²⁺, Co-AK; and Fe²⁺, Fe-AK) bound in its LID domain have been determined by X-ray crystallography to resolutions 1.8, 2.0, and 3.0 Å, respectively. The zinc and iron forms of the enzyme were crystallized in space group I222, whereas the cobalt-form crystals were C2. The presence of the metals was confirmed by calculation of anomalous difference maps and by X-ray fluorescence scans. The work presented here is the first report of a structure of a metal-containing AK from a Gram-negative bacterium. The native enzyme was crystallized, and only zinc was detected in the LID domain. Co-AK and Fe-AK were obtained by overexpressing the protein in Escherichia coli. Zn-AK and Fe-AK crystallized as monomers in the asymmetric unit, whereas Co-AK crystallized as a dimer. Nevertheless, all three crystal structures are very similar to each other, with the same LID domain topology, the only change being the presence of the different metal atoms. In the absence of any substrate, the LID domain of all holoforms of AK was present in a fully open conformational state. Normal mode analysis was performed to predict fluctuations of the LID domain along the catalytic pathway.     

Conformational changes in hemoglobin S (betaE6V) imposed by mutation of the beta Glu7-beta Lys132 salt bridge and detected by UV resonance Raman spectroscopy

The impact upon molecular structure of an additional point mutation adjacent to the existing E6V mutation in sickle cell hemoglobin was probed spectroscopically. The UV resonance Raman results show that the conformational consequences of mutating the salt bridge pair, betaGlu(7)-betaLys(132), are dependent on which residue of the pair is modified. The betaK132A mutants exhibit the spectroscopic signatures of the R --> T state transition in both the "hinge" and "switch" regions of the alpha(1)beta(2) interface. Both singly and doubly mutated hemoglobin (Hb) betaepsilon7Alpha exhibit the switch region signature for the R --> T quaternary state transition but not the hinge signature. The absence of this hinge region-associated quaternary change is the likely origin of the observed increased oxygen binding affinity for the Hb betaepsilon7Alpha mutants. The observed large decrease in the W3 alpha14beta15 band intensity for doubly mutated Hb betaepsilon7Alpha is attributed to an enhanced separation in the A helix-E helix tertiary contact of the beta subunits. The results for the Hb A betaGlu(7)-betaLys(132) salt bridge mutants demonstrate that attaining the T state conformation at the hinge region of the alpha(1)beta(2) dimer interface can be achieved through different intraglobin pathways; these pathways are subject to subtle mutagenic manipulation at sites well removed from the dimer interface.     

On the roles of substrate binding and hinge unfolding in conformational changes of adenylate kinase

We characterized the conformational change of adenylate kinase (AK) between open and closed forms by conducting five all-atom molecular-dynamics simulations, each of 100 ns duration. Different initial structures and substrate binding configurations were used to probe the pathways of AK conformational change in explicit solvent, and no bias potential was applied. A complete closed-to-open and a partial open-to-closed transition were observed, demonstrating the direct impact of substrate-mediated interactions on shifting protein conformation. The sampled configurations suggest two possible pathways for connecting the open and closed structures of AK, affirming the prediction made based on available x-ray structures and earlier works of coarse-grained modeling. The trajectories of the all-atom molecular-dynamics simulations revealed the complexity of protein dynamics and the coupling between different domains during conformational change. Calculations of solvent density and density fluctuations surrounding AK did not show prominent variation during the transition between closed and open forms. Finally, we characterized the effects of local unfolding of an important hinge near Pro¹⁷⁷ on the closed-to-open transition of AK and identified a novel mechanism by which hinge unfolding modulates protein conformational change. The local unfolding of Pro¹⁷⁷ hinge induces alternative tertiary contacts that stabilize the closed structure and prevent the opening transition.     

The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase. A highly stable salt bridge D118-K136 at the LID-CORE interface, contributing substantially to the total nonbonded LID-CORE interactions, was identified as a major factor that stabilizes the open conformation. Alternative transition pathways, such as AMPbd opening following LID opening, seem unlikely, as full transition events were not observed along this pathway. The simulation data indicate a high enthalpic penalty, possibly obstructing transitions along this route.     

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.     

Induced fit in guanidino kinases--comparison of substrate-free and transition state analog structures of arginine kinase

Arginine kinase (AK) is a member of the guanidino kinase family that plays an important role in buffering ATP concentration in cells with high and fluctuating energy demands. The AK specifically catalyzes the reversible phosphoryl transfer between ATP and arginine. We have determined the crystal structure of AK from the horseshoe crab (Limulus polyphemus) in its open (substrate-free) form. The final model has been refined at 2.35 A with a final R of 22.3% (R(free) = 23.7%). The structure of the open form is compared to the previously determined structure of the transition state analog complex in the closed form. Classically, the protein would be considered two domain, but dynamic domain (DynDom) analysis shows that most of the differences between the two structures can be considered as the motion between four rigid groups of nonsequential residues. ATP binds near a cluster of positively charged residues of a fixed dynamic domain. The other three dynamic domains close the active site with separate hinge rotations relative to the fixed domain. Several residues of key importance for the induced motion are conserved within the phosphagen kinase family, including creatine kinase. Substantial conformational changes are induced in different parts of the enzyme as intimate interactions are formed with both substrates. Thus, although induced fit occurs in a number of phosphoryl transfer enzymes, the conformational changes in phosphagen kinases appear to be more complicated than in prior examples.     

Explaining an unusually fast parasitic enzyme: folate tail-binding residues dictate substrate positioning and catalysis in Cryptosporidium hominis thymidylate synthase

The essential enzyme TS-DHFR from Cryptosporidium hominis undergoes an unusually rapid rate of catalysis at the conserved TS domain, facilitated by two nonconserved residues, Ala287 and Ser290, in the folate tail-binding region. Mutation of these two residues to their conserved counterparts drastically affects multiple steps of the TS catalytic cycle. We have determined the crystal structures of all three mutants (A287F, S290G, and A287F/S290G) in complex with active site ligands dUMP and CB3717. The structural data show two effects of the mutations: an increased distance between the ligands in the active site and increased flexibility of the folate ligand in the partially open enzyme state that precedes conformational change to the active catalytic state. The latter effect is able to be rescued by the mutants containing the A287F mutation. In addition, the conserved water network of TS is altered in each of the mutants. The structural results point to a role of the folate tail-binding residues in closely positioning ChTS ligands and restricting ligand flexibility in the partially open state to allow for a rapid transition to the active closed state and enhanced rate of catalysis. These results provide an explanation on how folate tail-binding residues at one end of the active site affect long-range interactions throughout the TS active site and validate these residues as targets for species-specific drug design.     

Crystal structure combined with genetic analysis of the Thermus thermophilus ribosome recycling factor shows that a flexible hinge may act as a functional switch

Ribosome recycling factor (RRF), in concert with elongation factor EF-G, is required for disassembly of the posttermination complex of the ribosome after release of polypeptides. The crystal structure of Thermus thermophilus RRF was determined at 2.6 A resolution. It is a tRNA-like L-shaped molecule consisting of two domains: a long three-helix bundle (domain 1) and a three-layer beta/alpha/beta sandwich (domain 2). Although the individual domain structures are similar to those of Thermotoga maritima RRF (Selmer et al., Science, 1999, 286:2349-2352), the interdomain angle differs by 33 degrees in two molecules, suggesting that the hinge between two domains is potentially flexible and responsive to different conditions of crystal packing. The hinge connects hydrophobic junctions of domains 1 and 2. The structure-based genetic analysis revealed the strong correlation between the hinge flexibility and the in vivo function of RRF. First, altering the hinge flexibility by making alanine or serine substitutions for large-size residues conserved at the hinge loop and nearby in domain 1 frequently gave rise to gain of function except a Pro residue conserved at the hinge loop. Second, the hinge defect resulting from a too relaxed hinge structure can be compensated for by secondary alterations in domain 1 that seem to increase the hydrophobic contact between domain 1 and the hinge loop. These results show that the hinge flexibility is vital for the function of RRF and that the steric interaction between the hinge loop and domains 1 and 2 restricts the interdomain angle and/or the hinge flexibility. These results indicate that RRF possesses an architectural difference from tRNA regardless of a resemblance to tRNA shape: RRF has a "gooseneck" elbow, whereas the tRNA elbow is rigid, and the direction of flex of RRF and tRNA is at a nearly right angle to each other. Moreover, surface electrostatic potentials of the two RRF proteins are dissimilar and do not mimic the surface potential of tRNA or EF-G. These properties will add a new insight into RRF, suggesting that RRF is more than a simple tRNA mimic.     

Essential dynamics sampling study of adenylate kinase: comparison to citrate synthase and implication for the hinge and shear mechanisms of domain motions

Essential dynamics sampling simulations of the domain conformations of unliganded Escherichia coli adenylate kinase have been performed to determine whether the ligand-induced closed-domain conformation is accessible to the open unliganded enzyme. Adenylate kinase is a three- domain protein with a central CORE domain and twoflanking domains, the LID and the NMPbind domains. The sampling simulations were applied to the CORE and NMPbind domain pair and the CORE and LID domain pair separately. One aim is to compare the results to those of a similar study on the enzyme citrate synthase to determine whether a similar domain-locking mechanism operates in adenylate kinase. Although for adenylate kinase the simulations suggest that the closed-domain conformation of the unliganded enzyme is at a slightly higher free energy than the open for both domain pairs, the results are radically different to those found for citrate synthase. In adenylate kinase the targeted domain conformations could always be achieved, whereas this was not the case in citrate synthase due to an apparent free-energy barrier between the open and closed conformations. Adenylate kinase has been classified as a protein that undergoes closure through a hinge mechanism, whereas citrate synthase has been assigned to the shear mechanism. This was quantified here in terms of the change in the number of interdomain contacting atoms upon closure which showed a considerable increase in adenylate kinase. For citrate synthase this number remained largely the same, suggesting that the domain faces slide over each other during closure. This suggests that shear and hinge mechanisms of domain closure may relate to the existence or absence of an appreciable barrier to closure for the unliganded protein, as the latter can hinge comparatively freely, whereas the former must follow a more constrained path. In general though it appears a bias toward keeping the unliganded enzyme in the open-domain conformation may be a common feature of domain enzymes.     

Mechanisms of domain closure in proteins

Certain enzymes respond to the binding of substrates and coenzymes by the closure of an active site that lies in a cleft between two domains. We have examined the mechanism of the domain closure in citrate synthase, for which atomic co-ordinates are available for "open" and "closed" forms. We show that the mechanism of domain closure involves small shifts and rotations of packed helices within the two domains and at their interface. Large motions of distant segments of the structure are the cumulative effect of the small relative shifts in intervening pairs of packed segments. These shifts are accommodated not by changes in packing but rather by small conformational changes in side-chains. We call this the helix interface shear mechanism of domain closure. The relative movements of packed helices follow the principles suggested by our recent study of insulin. This mechanism of domain closure is quite different from the hinge mechanisms that allow the rigid body movements of domains in immunoglobulins. The large interface between the domains of citrate synthase precludes a simple hinge mechanism for its conformational change. The helix interface shear mechanism of conformational change occurs in other enzymes that contain extensive domain-domain interfaces.     

Crystal structure of unligated guanylate kinase from yeast reveals GMP-induced conformational changes

The crystal structure of guanylate kinase (GK) from yeast (Saccharomyces cerevisiae) with a non-acetylated N terminus has been determined in its unligated form (apo-GK) as well as in complex with GMP (GK.GMP). The structure of apo-GK was solved with multiwavelength anomalous diffraction data and refined to an R-factor of 0.164 (R(free)=0.199) at 2.3 A resolution. The structure of GK.GMP was determined using the crystal structure of GK with an acetylated N terminus as the search model and refined to an R-factor of 0.156 (R(free)=0.245) at 1.9 A. GK belongs to the family of nucleoside monophosphate (NMP) kinases and catalyzes the reversible phosphoryl transfer from ATP to GMP. Like other NMP kinases, GK consists of three dynamic domains: the CORE, LID, and NMP-binding domains. Dramatic movements of the GMP-binding domain and smaller but significant movements of the LID domain have been revealed by comparing the structures of apo-GK and GK.GMP. apo-GK has a much more open conformation than the GK.GMP complex. Systematic analysis of the domain movements using the program DynDom shows that the large movements of the GMP-binding domain involve a rotation around an effective hinge axis approximately parallel with helix 3, which connects the GMP-binding and CORE domains. The C-terminal portion of helix 3, which connects to the CORE domain, has strikingly higher temperature factors in GK.GMP than in apo-GK, indicating that these residues become more mobile upon GMP binding. The results suggest that helix 3 plays an important role in domain movement. Unlike the GMP-binding domain, which moves toward the active center of the enzyme upon GMP binding, the LID domain moves away from the active center and makes the presumed ATP-binding site more open. Therefore, the LID domain movement may facilitate the binding of MgATP. The structure of the recombinant GK.GMP complex superimposes very well with that of the native GK.GMP complex, indicating that N-terminal acetylation does not have significant impact on the three-dimensional structure of GK.     

Backbone flexibility, conformational change, and catalysis in a phosphohexomutase from Pseudomonas aeruginosa

The enzyme phosphomannomutase/phosphoglucomutase (PMM/PGM) from the bacterium Pseudomonas aeruginosa is involved in the biosynthesis of several complex carbohydrates, including alginate, lipopolysaccharide, and rhamnolipid. Previous structural studies of this protein have shown that binding of substrates produces a rotation of the C-terminal domain, changing the active site from an open cleft in the apoenzyme into a deep, solvent inaccessible pocket where phosphoryl transfer takes place. We report herein site-directed mutagenesis, kinetic, and structural studies in examining the role of residues in the hinge between domains 3 and 4, as well as residues that participate in enzyme-substrate contacts and help form the multidomain "lid" of the active site. We find that the backbone flexibility of residues in the hinge region (e.g., mutation of proline to glycine/alanine) affects the efficiency of the reaction, decreasing k cat by approximately 10-fold and increasing K m by approximately 2-fold. Moreover, thermodynamic analyses show that these changes are due primarily to entropic effects, consistent with an increase in the flexibility of the polypeptide backbone leading to a decreased probability of forming a catalytically productive active site. These results for the hinge residues contrast with those for mutants in the active site of the enzyme, which have profound effects on enzyme kinetics (10 (2)-10 (3)-fold decrease in k cat/ K m) and also show substantial differences in their thermodynamic parameters relative to those of the wild-type (WT) enzyme. These studies support the concept that polypeptide flexibility in protein hinges may evolve to optimize and tune reaction rates.     

Mechanochemical coupling in the myosin motor domain. I. Insights from equilibrium active-site simulations

Although the major structural transitions in molecular motors are often argued to couple to the binding of Adenosine triphosphate (ATP), the recovery stroke in the conventional myosin has been shown to be dependent on the hydrolysis of ATP. To obtain a clearer mechanistic picture for such "mechanochemical coupling" in myosin, equilibrium active-site simulations with explicit solvent have been carried out to probe the behavior of the motor domain as functions of the nucleotide chemical state and conformation of the converter/relay helix. In conjunction with previous studies of ATP hydrolysis with different active-site conformations and normal mode analysis of structural flexibility, the results help establish an energetics-based framework for understanding the mechanochemical coupling. It is proposed that the activation of hydrolysis does not require the rotation of the lever arm per se, but the two processes are tightly coordinated because both strongly couple to the open/close transition of the active site. The underlying picture involves shifts in the dominant population of different structural motifs as a consequence of changes elsewhere in the motor domain. The contribution of this work and the accompanying paper [] is to propose the actual mechanism behind these "population shifts" and residues that play important roles in the process. It is suggested that structural flexibilities at both the small and large scales inherent to the motor domain make it possible to implement tight couplings between different structural motifs while maintaining small free-energy drops for processes that occur in the detached states, which is likely a feature shared among many molecular motors. The significantly different flexibility of the active site in different X-ray structures with variable level arm orientations supports the notation that external force sensed by the lever arm may transmit into the active site and influence the chemical steps (nucleotide hydrolysis and/or binding).     

Investigation of the mechanism of domain closure in citrate synthase by molecular dynamics simulation

Six, 2 ns molecular dynamics simulations have been performed on the homodimeric enzyme citrate synthase. In three, both monomers were started from the open, unliganded X-ray conformation. In the remaining three, both monomers started from a closed, liganded X-ray conformation, with the ligands removed. Projecting the motion from the simulations onto the experimental domain motion revealed that the free-energy profile is rather flat around the open conformation, with steep sides. The most closed conformations correspond to hinge-bending angles of 12-14 compared to the 20 degrees that occurs upon the binding of oxaloacetate. It is also found that the open, unliganded X-ray conformation is situated at the edge of the steep rise in free energy, although conformations that are about 5 degrees more open were sampled. A rigid-body essential dynamics analysis of the combined open trajectories has shown that domain motions in the direction of the closed X-ray conformation are compatible with the natural domain motion of the unliganded protein, which has just two main degrees of freedom. The simulations starting from the closed conformation suggest a free-energy profile with a small barrier in going from the closed to open conformation. A combined essential dynamics and hinge-bending analysis of a trajectory that spontaneously converts from the closed to open state shows an almost exact correspondence to the experimental transition that occurs upon ligand binding. The simulations support the conclusion from an earlier analysis of the experimental transition that the beta-hairpin acts as a mechanical hinge by attaching the small domain to the large domain through a conserved main-chain hydrogen bond and salt-bridges, and allowing rotation to occur via its two flexible termini. The results point to a mechanism of domain closure in citrate synthase that has analogy to the process of closing a door.     

Cloning and characterization of adenylate kinase from Chlamydia pneumoniae

Chlamydiae proliferate only within the infected host cells and are thought to be "energy parasites," because they take up ATP from the host cell as an energy source. In the present study, we isolated from Chlamydia pneumoniae the gene encoding adenylate kinase (AK). Using the enzyme produced in Escherichia coli, its properties were characterized. K(m) values for AMP and for ADP of the purified C. pneumoniae AK (AKcpn) were each 330 microm, which is significantly higher than the reported values of other AKs, whereas K(m) for ATP was 24 microm, which was rather lower than others. AKcpn contains 1 g atom of zinc/mol of 24,000-dalton protein. Mass spectrometric analysis of AKcpn and analysis of properties of mutated AKcpn strongly suggested that zinc is associated with four cysteine residues in the LID domain of the enzyme. The apo-AKcpn that lost zinc retained AK activity, although K(m) for AMP of apo-AKcpn increased about 2-fold and V(max) decreased about one-half from that of holo-AKcpn. The apo-AKcpn was more thermolabile and sensitive to trypsin digestion than the holo-AKcpn. Moreover, the recovery in vitro of the AK activity during the renaturation process of the denatured apo-AKcpn was dependent on zinc. A mutated protein in which cysteine residues in the LID domain were substituted by other amino acids lost both zinc and enzyme activity. The mutated protein was more sensitive to protease than the apo-AKcpn. These results indicate that zinc in AKcpn, although not essential for the catalysis, stabilizes the enzyme and probably plays a crucial role in proper folding of the protein. Furthermore, the catalytic properties of AKcpn suggest a distinctive regulatory mechanism in the metabolism compared with AKs in other organisms.     

Molecular dynamics simulations of DNA curvature and flexibility: helix phasing and premelting

Recent studies of DNA axis curvature and flexibility based on molecular dynamics (MD) simulations on DNA are reviewed. The MD simulations are on DNA sequences up to 25 base pairs in length, including explicit consideration of counterions and waters in the computational model. MD studies are described for ApA steps, A-tracts, for sequences of A-tracts with helix phasing. In MD modeling, ApA steps and A-tracts in aqueous solution are essentially straight, relatively rigid, and exhibit the characteristic features associated with the B'-form of DNA. The results of MD modeling of A-tract oligonucleotides are validated by close accord with corresponding crystal structure results and nuclear magnetic resonance (NMR) nuclear Overhauser effect (NOE) and residual dipolar coupling (RDC) structures of d(CGCGAATTCGCG) and d(GGCAAAAAACGG). MD simulation successfully accounts for enhanced axis curvature in a set of three sequences with phased A-tracts studied to date. The primary origin of the axis curvature in the MD model is found at those pyrimidine/purine YpR "flexible hinge points" in a high roll, open hinge conformational substate. In the MD model of axis curvature in a DNA sequence with both phased A-tracts and YpR steps, the A-tracts appear to act as positioning elements that make the helix phasing more precise, and key YpR steps in the open hinge state serve as curvature elements. Our simulations on a phased A-tract sequence as a function of temperature show that the MD simulations exhibit a premelting transition in close accord with experiment, and predict that the mechanism involves a B'-to-B transition within A-tracts coupled with the prediction of a transition in key YpR steps from the high roll, open hinge, to a low roll, closed hinge substate. Diverse experimental observations on DNA curvature phenomena are examined in light of the MD model with no serious discrepancies. The collected MD results provide independent support for the "non-A-tract model" of DNA curvature. The "junction model" is indicated to be a special case of the non-A-tract model when there is a Y base at the 5' end of an A-tract. In accord with crystallography, the "ApA wedge model" is not supported by MD.     

Modeling of DNA binding to the condensin hinge domain using molecular dynamics simulations guided by atomic force microscopy

The condensin protein complex compacts chromatin during mitosis using its DNA-loop extrusion activity. Previous studies proposed scrunching and loop-capture models as molecular mechanisms for the loop extrusion process, both of which assume the binding of double-strand (ds) DNA to the hinge domain formed at the interface of the condensin subunits Smc2 and Smc4. However, how the hinge domain contacts dsDNA has remained unknown. Here, we conducted atomic force microscopy imaging of the budding yeast condensin holo-complex and used this data as basis for coarse-grained molecular dynamics simulations to model the hinge structure in a transient open conformation. We then simulated the dsDNA binding to open and closed hinge conformations, predicting that dsDNA binds to the outside surface when closed and to the outside and inside surfaces when open. Our simulations also suggested that the hinge can close around dsDNA bound to the inside surface. Based on these simulation results, we speculate that the conformational change of the hinge domain might be essential for the dsDNA binding regulation and play roles in condensin-mediated DNA-loop extrusion.     

Alternative S2 hinge regions of the myosin rod differentially affect muscle function, myofibril dimensions and myosin tail length

Muscle myosin heavy chain (MHC) rod domains intertwine to form alpha-helical coiled-coil dimers; these subsequently multimerize into thick filaments via electrostatic interactions. The subfragment 2/light meromyosin "hinge" region of the MHC rod, located in the C-terminal third of heavy meromyosin, may form a less stable coiled-coil than flanking regions. Partial "melting" of this region has been proposed to result in a helix to random-coil transition. A portion of the Drosophila melanogaster MHC hinge is encoded by mutually exclusive alternative exons 15a and 15b, the use of which correlates with fast (hinge A) or slow (hinge B) muscle physiological properties. To test the functional significance of alternative hinge regions, we constructed transgenic fly lines in which fast muscle isovariant hinge A was switched for slow muscle hinge B in the MHC isoforms of indirect flight and jump muscles. Substitution of the slow muscle hinge B impaired flight ability, increased sarcomere lengths by approximately 13% and resulted in minor disruption to indirect flight muscle sarcomeric structure compared with a transgenic control. With age, residual flight ability decreased rapidly and myofibrils developed peripheral defects. Computational analysis indicates that hinge B has a greater coiled-coil propensity and thus reduced flexibility compared to hinge A. Intriguingly, the MHC rod with hinge B was approximately 5 nm longer than myosin with hinge A, consistent with the more rigid coiled-coil conformation predicted for hinge B. Our study demonstrates that hinge B cannot functionally substitute for hinge A in fast muscle types, likely as a result of differences in the molecular structure of the rod, subtle changes in myofibril structure and decreased ability to maintain sarcomere structure in indirect flight muscle myofibrils. Thus, alternative hinges are important in dictating the distinct functional properties of myosin isoforms and the muscles in which they are expressed.