Biomedical subjects

Qiang Cui

Publications and source records attributed to Qiang Cui.

At least 19 recordsLinked to original sources

PubMed2006

"Proton holes" in long-range proton transfer reactions in solution and enzymes: A theoretical analysis.

Proton transfers are fundamental to chemical processes in solution and biological systems. Often, the well-known Grotthuss mechanism is assumed where a series of sequential "proton hops" initiates from the donor and combines to produce the net transfer of a positive charge over a long distance. Although direct experimental evidence for the sequential proton hopping has been obtained recently, alternative mechanisms may be possible in complex molecular systems. To understand these events, all accessible protonation states of the mediating groups should be considered. This is exemplified by transfers through water where the individual water molecules can exist in three protonation states (water, hydronium, and hydroxide); as a result, an alternative to the Grotthuss mechanism for a proton transfer through water is to generate a hydroxide by first protonating the acceptor and then transfer the hydroxide toward the donor through water. The latter mechanism can be most generally described as the transfer of a "proton hole" from the acceptor to the donor where the "hole" characterizes the deprotonated state of any mediating molecule. This pathway is distinct and is rarely considered in the discussion of proton-transfer processes. Using a calibrated quantum mechanical/molecular mechanical (QM/MM) model and an effective sampling technique, we study proton transfers in two solution systems and in Carbonic Anhydrase II. Although the relative weight of the "proton hole" and Grotthuss mechanisms in a specific system is difficult to determine precisely using any computational approach, the current study establishes an energetics motivated framework that hinges on the donor/acceptor pKa values and electrostatics due to the environment to argue that the "proton hole" transfer is likely as important as the classical Grotthuss mechanism for proton transport in many complex molecular systems.

Acetic Acid↗

PubMed2006

The temperature dependence of salt-protein association is sequence specific.

Molecular dynamics (MD) simulations are used to probe the origin of the unexpected temperature dependence of salt accumulation in the C-terminal region of the protein human lymphotactin. As in previous MD simulations, sodium ions accumulate in an enhanced manner near the C-terminal helix at the lower temperature, while the temperature dependence of chloride accumulation is much weaker and slightly positive. In a designed mutant in which all positively charged residues in the C-terminal helix are replaced with neutral polar groups (Ser), the unexpected temperature dependence of the sodium ions is no longer observed. Therefore, these simulations convincingly verified the previous hypothesis that the temperature dependence of ion-protein association is sensitive to the local sequence. This is explained qualitatively in terms of the entropy of association between charged species in solution. These findings have general implications for the interpretation of thermodynamic quantities associated with binding events where ion release is important, such as protein-DNA interactions.

Amino Acid Sequence↗

PubMed2006

Optimization and evaluation of a coarse-grained model of protein motion using x-ray crystal data.

Simple coarse-grained models, such as the Gaussian network model, have been shown to capture some of the features of equilibrium protein dynamics. We extend this model by using atomic contacts to define residue interactions and introducing more than one interaction parameter between residues. We use B-factors from 98 ultra-high resolution (<or=1.0 A) x-ray crystal structures to optimize the interaction parameters. The average correlation between Gaussian network-model fluctuation predictions and the B-factors is 0.64 for the data set, consistent with a previous large-scale study. By separating residue interactions into covalent and noncovalent, we achieve an average correlation of 0.74, and addition of ligands and cofactors further improves the correlation to 0.75. However, further separating the noncovalent interactions into nonpolar, polar, and mixed yields no significant improvement. The addition of simple chemical information results in better prediction quality without increasing the size of the coarse-grained model.

Calmodulin↗

PubMed2006

Effects of temperature and salt concentration on the structural stability of human lymphotactin: insights from molecular simulations.

Extensive molecular dynamics (MD) simulations ( approximately 70 ns total) with explicit solvent molecules and salt ions are carried out to probe the effects of temperature and salt concentration on the structural stability of the human Lymphotactin (hLtn). The distribution of ions near the protein surface and the stability of various structural motifs are observed to exhibit interesting dependence on the local sequence and structure. Whereas chloride association to the protein is overall enhanced as the temperature increases, the sodium distribution in the C-terminal helical region and, to a smaller degree, the chloride distribution in the same region are found higher at the lower temperature. The similar trend is also observed in nonlinear Poisson-Boltzmann calculations with a temperature-dependent water dielectric constant, once conformational averaging over a series of MD snapshots is done. The unexpected temperature dependence in the ion distribution is explained on the basis of the cancellation of association entropy for ion-side chain pairs of opposite-charge and like-charge characters, which have positive and negative contributions, respectively. The C-terminal helix is observed to partially melt whereas a short beta strand forms at the higher temperature with little salt dependence. The N-terminal region, by contrast, develops partial helical structure at a higher salt concentration. These observed behaviors are consistent with solvent and salt screening playing an important role in stabilizing the canonical chemokine fold of hLtn.

Chlorides↗

PubMed2006

Reconciling the "old" and "new" views of protein allostery: a molecular simulation study of chemotaxis Y protein (CheY).

A combination of thirty-two 10-ns-scale molecular dynamics simulations were used to explore the coupling between conformational transition and phosphorylation in the bacteria chemotaxis Y protein (CheY), as a simple but representative example of protein allostery. Results from these simulations support an activation mechanism in which the beta4-alpha4 loop, at least partially, gates the isomerization of Tyr106. The roles of phosphorylation and the conserved Thr87 are deemed indirect in that they stabilize the active configuration of the beta4-alpha4 loop. The indirect role of the activation event (phosphorylation) and/or conserved residues in stabilizing, rather than causing, specific conformational transition is likely a feature in many signaling systems. The current analysis of CheY also helps to make clear that neither the "old" (induced fit) nor the "new" (population shift) views for protein allostery are complete, because they emphasize the kinetic (mechanistic) and thermodynamic aspects of allosteric transitions, respectively. In this regard, an issue that warrants further analysis concerns the interplay of concerted collective motion and sequential local structural changes in modulating cooperativity between distant sites in biomolecules.

Allosteric Regulation↗

PubMed2006

A finite element framework for studying the mechanical response of macromolecules: application to the gating of the mechanosensitive channel MscL.

The gating pathways of mechanosensitive channels of large conductance (MscL) in two bacteria (Mycobacterium tuberculosis and Escherichia coli) are studied using the finite element method. The phenomenological model treats transmembrane helices as elastic rods and the lipid membrane as an elastic sheet of finite thickness; the model is inspired by the crystal structure of MscL. The interactions between various continuum components are derived from molecular-mechanics energy calculations using the CHARMM all-atom force field. Both bacterial MscLs open fully upon in-plane tension in the membrane and the variation of pore diameter with membrane tension is found to be essentially linear. The estimated gating tension is close to the experimental value. The structural variations along the gating pathway are consistent with previous analyses based on structural models with experimental constraints and biased atomistic molecular-dynamics simulations. Upon membrane bending, neither MscL opens substantially, although there is notable and nonmonotonic variation in the pore radius. This emphasizes that the gating behavior of MscL depends critically on the form of the mechanical perturbation and reinforces the idea that the crucial gating parameter is lateral tension in the membrane rather than the curvature of the membrane. Compared to popular all-atom-based techniques such as targeted or steered molecular-dynamics simulations, the finite element method-based continuum-mechanics framework offers a unique alternative to bridge detailed intermolecular interactions and biological processes occurring at large spatial scales and long timescales. It is envisioned that such a hierarchical multiscale framework will find great value in the study of a variety of biological processes involving complex mechanical deformations such as muscle contraction and mechanotransduction.

Cell Membrane↗

PubMed2006

Development of effective quantum mechanical/molecular mechanical (QM/MM) methods for complex biological processes.

Motivated by the long-term goal of understanding vectorial biological processes such as proton transport (PT) in biomolecular ion pumps, a number of developments were made to establish combined quantum mechanical/molecular mechanical (QM/MM) methods suitable for studying chemical reactions involving significant charge separation in the condensed phase. These developments were summarized and discussed with representative problems. Specifically, free energy perturbation and boundary potential methods for treating long-range electrostatics were implemented to test the robustness of QM/MM results for protein systems. It was shown that consistent models with sufficient sampling were able to produce quantitatively satisfactory results, such as pKa for titritable groups in the interior of T4-lysozyme, while an inconsistent treatment of electrostatics or lack of sufficient sampling may produce incorrect results. Modifications were made to an approximate density functional theory (SCC-DFTB) to improve the description of proton affinity and hydrogen-bonding, which are crucial for the treatment of PT in polar systems. Test calculations on water autoionization showed clearly that both improvements are necessary for quantitatively reliable results. Finally, the newly established SCC-DFTB/MM-GSBP protocol was used to explore mechanistic issues in carbonic anhydrase (CA). Preliminary results suggest that PT in CA occurs mainly through short water wires containing two water molecules in a thermally activated fashion. Although longer water wires occur with similar frequencies, PT along those pathways, on average, has substantially higher barriers, a result not expected based on previous studies. The fluctuations of water molecules peripheral to the water wire were found to make a larger impact on the PT energetics compared to polar protein residues in the active site, which are largely pre-organized and therefore have less tendency to reorganize during the reaction.

Computer Simulation↗

PubMed2006

The use of a generalized born model for the analysis of protein conformational transitions: a comparative study with explicit solvent simulations for chemotaxis Y protein (CheY).

To investigate whether implicit solvent models are appropriate for mechanistic studies of conformational transition in proteins, a recently developed generalized Born model (GBSW) was applied to a small signaling protein, chemotaxis protein Y (CheY), with different combinations of the phosphorylation state and conformation of the system; the results were compared to explicit solvent simulations using a stochastic boundary condition. The subtle but distinct conformational transitions involved in CheY activation makes the system ideally suited for comparing implicit and explicit solvent models because these conformational transitions are potentially accessible in both types of simulations. The structural and dynamical properties analyzed include not only those localized to the active site region but also throughout the protein, such as sidechain methyl group order parameters, backbone hydrogen bonding lifetime and occupancy as well as principal components of the trajectories. Overall, many properties were well reproduced by the GBSW simulations when compared with the explicit solvent calculations, although a number of observations consistently point to the suggestion that the current parameterization of the GBSW model tends to overestimate hydrogen-bonding interactions involving both charged groups and (charge-neutral) backbone atoms. This deficiency led to overstabilization of certain secondary structural motifs and more importantly, qualitatively different behaviors for the active site groups (Thr 87, Ala 88, the beta4-alpha4 loop) in response to phosphorylation, when compared with explicit solvent simulations. The current study highlights the value of carrying out both explicit and implicit solvent simulations for complementary mechanistic insights in the analysis of conformational transition in biomolecules.

Bacterial Proteins↗

PubMed2006

Interpreting correlated motions using normal mode analysis.

With the increased popularity of normal mode analyses in structural biology, it is important to carefully consider how to best utilize the results for gaining biological insights without over interpretation. The discussion in this article argues that for the purpose of identifying correlated motions in biomolecules, a case separate from concomitant conformational changes of structural motifs, it is generally important to use a large number of normal modes. This is illustrated through three increasingly complex examples. The simplest case includes two bilinearly coupled harmonic oscillators and serves as a straightforward problem where the important considerations are explicit and transparent. The argument is then generalized to include a system of N-coupled harmonic oscillators and finally to a realistic biomolecule. Although a small number of normal modes are useful for probing structural flexibility, it is clear that a much larger number of modes are required for properly investigating correlated motions in biomolecules.

Algorithms↗

PubMed2005

Reduction of imines by hydroxycyclopentadienyl ruthenium hydride: intramolecular trapping evidence for hydride and proton transfer outside the coordination sphere of the metal.

Reduction of imines by [2,5-Ph2-3,4-Tol2(eta(5)-C4COH)]Ru(CO)2H (2) produces kinetically stable ruthenium amine complexes. Reduction of an imine by 2 in the presence of an external amine trap gives only the complex of the newly generated amine. Reaction of 2 with H2N-p-C6H4N=CHPh (11), which contains an intramolecular amine trap, gave a 1:1 mixture of [2,5-Ph2-3,4-Tol2(eta(4)-C4CO)](CO)2RuNH(CH2Ph)(C6H4-p-NH2) (8), formed by coordination of the newly generated amine to the ruthenium center, and [2,5-Ph2-3,4-Tol2(eta(4)-C4CO)](CO)2RuNH2C6H4-p-NHCH2Ph (9), formed by coordination of the amine already present in the substrate. These results require transfer of hydrogen to the imine outside the coordination sphere of the metal to give a coordinatively unsaturated intermediate that can be trapped inside the initial solvent cage. Amine diffusion from the solvent cage must be much slower than coordination to the metal center. Mechanisms requiring prior coordination of the substrate to ruthenium would have led only to 8 and can be eliminated.

Crystallography, X-Ray↗

PubMed2005

pKa calculations in solution and proteins with QM/MM free energy perturbation simulations: a quantitative test of QM/MM protocols.

The accuracy of biological simulations depends, in large part, on the treatment of electrostatics. Due to the availability of accurate experimental values, calculation of pKa provides stringent evaluation of computational methods. The generalized solvent boundary potential (GSBP) and Ewald summation electrostatic treatments were recently implemented for combined quantum mechanical and molecular mechanics (QM/MM) simulations by our group. These approaches were tested by calculating pKa shifts due to differences in electronic structure and electrostatic environment; the shifts were determined for a series of small molecules in solution, using various electrostatic treatments, and two residues (His 31, Lys 102) in the M102K T4-lysozyme mutant with large pKa shifts, using the GSBP approach. The calculations utilized a free energy perturbation scheme with the QM/MM potential function involving the self-consistent charge density functional tight binding (SCC-DFTB) and CHARMM as the QM and MM methods, respectively. The study of small molecules demonstrated that inconsistent electrostatic models produced results that were difficult to correct in a robust manner; by contrast, extended electrostatics, GSBP, and Ewald simulations produced consistent results once a bulk solvation contribution was carefully chosen. In addition to the electrostatic treatment, the pKa shifts were also sensitive to the level of the QM method and the scheme of treating QM/MM Coulombic interactions; however, simple perturbative corrections based on SCC-DFTB/CHARMM trajectories and higher level single point energy calculations were found to give satisfactory results. Combining all factors gave a root-mean-square difference of 0.7 pKa units for the relative pKa values of the small molecules compared to experiment. For the residues in the lysozyme, an accurate pKa shift was obtained for His 31 with multiple nanosecond simulations. For Lys 102, however, the pKa shift was estimated to be too large, even after more than 10 nanosecond simulations for each lambda window; the difficulty was due to the significant, but slow, reorganization of the protein and water structure when Lys 102 was protonated. The simulations support that Lys 102 is deprotonated in the X-ray structure and the protein is highly destabilized when this residue is protonated.

Algorithms↗

PubMed2005

Comparison of mode analyses at different resolutions applied to nucleic acid systems.

More than two decades of different types of mode analyses has shown that these techniques can be useful in describing large-scale motions in protein systems. A number of mode analyses are available and include quasiharmonics, classical normal mode, block normal mode, and the elastic network model. Each of these methods has been validated for protein systems and this variety allows researchers to choose the technique that gives the best compromise between computational cost and the level of detail in the calculation. These same techniques have not been systematically tested for nucleic acid systems, however. Given the differences in interactions and structural features between nucleic acid and protein systems, the validity of these techniques in the protein regime cannot be directly translated into validity in the nucleic acid realm. In this work, we investigate the usefulness of the above mode analyses as applied to two RNA systems, i.e., the hammerhead ribozyme and a guanine riboswitch. We show that classical normal-mode analysis can match the magnitude and direction of residue fluctuations from the more detailed, anharmonic technique, quasiharmonic analysis of a molecular dynamics trajectory. The block normal-mode approximation is shown to hold in the nucleic acid systems studied. Only the mode analysis at the lowest level of detail, the elastic network model, produced mixed results in our calculations. We present data that suggest that the elastic network model, with the popular parameterization, is not best suited for systems that do not have a close packed structure; this observation also hints at why the elastic network model has been found to be valid for many globular protein systems. The different behaviors of block normal-mode analysis and the elastic network model, which invoke similar degrees of coarse-graining to the dynamics but use different potentials, suggest the importance of applying a heterogeneous potential function in a robust analysis of the dynamics of biomolecules, especially those that are not closely packed. In addition to these comparisons, we briefly discuss insights into the conformational space available to the hammerhead ribozyme.

Algorithms↗

PubMed2005

Benchmark calculations of proton affinities and gas-phase basicities of molecules important in the study of biological phosphoryl transfer.

Benchmark calculations of proton affinities and gas-phase basicities of molecules most relevant to biological phosphoryl transfer reactions are presented and compared with available experimental results. The accuracy of proton affinity and gas-phase basicity results obtained from several multi-level model chemistries (CBS-QB3, G3B3, and G3MP2B3) and density-functional quantum models (PBE0, B1B95, and B3LYP) are assessed and compared. From these data, a set of empirical bond enthalpy, entropy, and free energy corrections are introduced that considerably improve the accuracy and predictive capability of the methods. These corrections are applied to the prediction of proton affinity and gas-phase basicity values of important biological phosphates and phosphoranes for which experimental data does not currently exist. Comparison is made with results from semiempirical quantum models that are commonly employed in hybrid quantum mechanical/molecular mechanical simulations. Data suggest that the design of improved semiempirical quantum models with increased accuracy for relative proton affinity values is necessary to obtain quantitative accuracy for phosphoryl transfer reactions in solution, enzymes, and ribozymes.

Algorithms↗

PubMed2005

Reliable treatment of electrostatics in combined QM/MM simulation of macromolecules.

A robust approach for dealing with electrostatic interactions for spherical boundary conditions has been implemented in the QM/MM framework. The development was based on the generalized solvent boundary potential (GSBP) method proposed by Im et al. [J. Chem. Phys. 114, 2924 (2001)], and the specific implementation was applied to the self-consistent-charge density-functional tight-binding approach as the quantum mechanics (QM) level, although extension to other QM methods is straightforward. Compared to the popular stochastic boundary-condition scheme, the new protocol offers a balanced treatment between quantum mechanics/molecular mechanics (QM/MM) and MM/MM interactions; it also includes the effect of the bulk solvent and macromolecule atoms outside of the microscopic region at the Poisson-Boltzmann level. The new method was illustrated with application to the enzyme human carbonic anhydrase II and compared to stochastic boundary-condition simulations using different electrostatic treatments. The GSBP-based QM/MM simulations were most consistent with available experimental data, while conventional stochastic boundary simulations yielded various artifacts depending on different electrostatic models. The results highlight the importance of carefully treating electrostatics in QM/MM simulations of biomolecules and suggest that the commonly used truncation schemes should be avoided in QM/MM simulations, especially in simulations that involve extensive conformational samplings. The development of the GSBP-based QM/MM protocol has opened up the exciting possibility of studying chemical events in very complex biomolecular systems in a multiscale framework.

Carbonic Anhydrases↗

PubMed2005

Hydrogen elimination from a hydroxycyclopentadienyl ruthenium(II) hydride: study of hydrogen activation in a ligand-metal bifunctional hydrogenation catalyst.

At high temperatures in toluene, [2,5-Ph(2)-3,4-Tol(2)(eta(5)-C(4)COH)]Ru(CO)(2)H (3) undergoes hydrogen elimination in the presence of PPh(3) to produce the ruthenium phosphine complex [2,5-Ph(2)-3,4-Tol(2)-(eta(4)-C(4)CO)]Ru(PPh(3))(CO)(2) (6). In the absence of alcohols, the lack of RuH/OD exchange, a rate law first order in Ru and zero order in phosphine, and kinetic deuterium isotope effects all point to a mechanism involving irreversible formation of a transient dihydrogen ruthenium complex B, loss of H(2) to give unsaturated ruthenium complex A, and trapping by PPh(3) to give 6. DFT calculations showed that a mechanism involving direct transfer of a hydrogen from the CpOH group to form B had too high a barrier to be considered. DFT calculations also indicated that an alcohol or the CpOH group of 3 could provide a low energy pathway for formation of B. PGSE NMR measurements established that 3 is a hydrogen-bonded dimer in toluene, and the first-order kinetics indicate that two molecules of 3 are also involved in the transition state for hydrogen transfer to form B, which is the rate-limiting step. In the presence of ethanol, hydrogen loss from 3 is accelerated and RuD/OH exchange occurs 250 times faster than in its absence. Calculations indicate that the transition state for dihydrogen complex formation involves an ethanol bridge between the acidic CpOH and hydridic RuH of 3; the alcohol facilitates proton transfer and accelerates the reversible formation of dihydrogen complex B. In the presence of EtOH, the rate-limiting step shifts to the loss of hydrogen from B.

Catalysis↗

PubMed2005

Direct determination of reaction paths and stationary points on potential of mean force surfaces.

A simulation approach is introduced for directly determining reaction paths and stationary points on potential of mean force (PMF) surfaces associated with molecular events that occur in complex environments. The nudged elastic band approach was employed to search for steepest descent paths on the PMF surface using the relevant PMF derivatives from a series of local simulations. The steepest descent path on the PMF surface corresponds to the minimum PMF path (i.e. the path with the lowest local PMF barrier), which contains important information about stationary points (e.g. saddle points) on the PMF surface, which in turn can provide useful insights into the thermodynamics and kinetics for the process of interest. By working with the PMF defined in a low-dimensional sub-space rather than a potential energy function of full molecular dimensionality, the main features of the process under study are concisely represented and the orthogonal degrees of freedom are adequately sampled with the appropriate canonical distribution at the desired temperature (e.g. 300 K). Therefore, minimum PMF paths carry statistically meaningful mechanistic information and are complementary to reaction paths of full molecular dimensionality proposed in previous studies. The NEB based path optimization method is direct in the sense that no information regarding the global PMF surface is necessary for the determination of the local reaction path and stationary points along this path. Since only low-dimensional quantities (paths) are searched for, the PMF-path method is expected to scale better in terms of dimension of the PMF sub-space than those aims to fully explore multi-dimensional PMF surfaces. Test applications on simple molecular systems, the alanine di-peptide in vacuum and in solution and a micro-solvated proton-wire, indicate that reliable PMF paths can be determined for both conformational isomerization and chemical reaction processes. However, highly accurate PMF derivatives are required for determining more quantitative observables, such as the free energy profile along the minimum PMF path. Therefore, effective numerical algorithms for calculating local PMF derivatives and systematic protocols for defining the relevant sub-space are the main focus in the near future. Finally, we emphasize that the minimum PMF path defined here includes thermal (e.g. entropic) effects associated with the orthogonal degrees of freedom, but finite kinetic energies associated with the PMF degrees of freedom are not included; this can be improved by adopting a different definition of the reaction path, such as the maximum flux path, on the PMF surface, or thermally sampling all degrees of freedom orthogonal to the one-dimensional path.

Alanine↗

PubMed2004

QM/MM studies of the enzyme-catalyzed dechlorination of 4-chlorobenzoyl-CoA provide insight into reaction energetics.

The conversion of 4-chlorobenzoyl-CoA to 4-hydroxybenzoyl-CoA catalyzed by 4-chlorobenzoyl-CoA dehalogenase is investigated using combined QM/MM approaches. The calculated potential of mean force at the PM3/CHARMM level supports the proposed nucleophilic aromatic substitution mechanism. In particular, a Meisenheimer intermediate was found, stabilized by hydrogen bonds between the benzoyl carbonyl of the ligand and two backbone amide NHs at positions 64 and 114. Mutation of Gly113 to Ala significantly increases the barrier by disrupting the hydrogen bond with the Gly114 backbone. The formation of the Meisenheimer complex is accompanied by significant charge redistribution and structural changes in the substrate benzoyl moiety, consistent with experimental observations. Theoretical results suggest that the reaction rate is limited by the formation of the Meisenheimer complex, rather than by its decomposition. A kinetic model based on the calculated free energy profile is found to be consistent with the experimental time course data.

Acyl Coenzyme A↗

PubMed2004

Normal-mode analysis suggests protein flexibility modulation throughout RNA polymerase's functional cycle.

To explore the domain-scale flexibility of bacterial RNA polymerase (RNAP) throughout its functional cycle, block normal-mode analyses (BNM) were performed on several important functional states, including the holoenzyme, the core complex, a model of RNAP bound to primarily duplex DNA, and a model of the ternary elongation complex. The calculations utilized a molecular mechanics (MM) force field with physical interactions; this is made possible by the use of BNM and the implementation of a sparse-matrix diagonalization routine. The use of homology models necessitated the MM force field rather than the simpler elastic network model (ENM). From the MM/BNM, we have systematically and semiquantitatively calculated the atomic fluctuations in the four functional states without bias due to crystal packing or other artifactual forces. We have observed that both alpha subunits and the omega subunit are rigid, in line with their roles as structural motifs that are not mechanistically involved in RNAP's functional cycle. It has been observed that the beta subunit has two highly mobile domains; these are commonly known as the beta1 and beta2 domains. Our calculations suggest that the flexibility of these domains is modulated throughout the functional cycle and that they move entirely independently of each other unless DNA is bound. From an energetic perspective, we have shown the beta2 domain can flex into and out of the cleft, forming interactions with DNA in the TEC as has been previously proposed. Our calculations also confirm that the beta' subunit's likely flexibility into and out of the DNA binding cleft is energetically allowed. These two observations validate that both of the RNAP crab claw's pincers are mobile, as both beta and beta' have substantial flexibility.

Bacterial Proteins↗