Mechanism of the conformational change of the F1-ATPase β subunit revealed by free energy simulations

F(1)-ATPase is an ATP-driven rotary motor enzyme. The β subunit changes its conformation from an open to a closed form upon ATP binding. The motion in the β subunit is regarded as a major driving force for rotation of the central stalk. In this Article, we explore the conformational change of the β subunit using all-atom free energy simulations with explicit solvent and propose a detailed mechanism for the conformational change. The β subunit conformational change is accomplished roughly in two characteristic steps: changing of the hydrogen-bond network around ATP and the dynamic movement of the C-terminal domain via sliding of the B-helix. The details of the former step agree well with experimental data. In the latter step, sliding of the B-helix enhances the hydrophobic stabilization due to the exclusion of water molecules from the interface and improved packing in the hydrophobic core. This step contributes to a decrease in free energy, leading to the generation of torque in the F(1)-ATPase upon ATP binding.     

Rotation mechanism of F1-ATPase: crucial importance of the water entropy effect

We propose a novel picture of the rotation mechanism of F(1)-ATPase, a rotary-motor protein complex. Entropy, which originates from the translational displacement of water molecules, is treated as the key factor in the proposal. We calculate the water entropy gains upon formation of the α-β, α-γ, and β-γ subunit pairs. The gain is given as the difference between the hydration entropy of a subunit pair and the sum of the hydration entropies of the separate subunits forming the pair. The calculation is made using a hybrid of a statistical-mechanical theory for molecular liquids and morphometric approach. The water entropy gain is considered as a measure of tightness of the packing at each subunit interface. The results are highly correlated with the numbers of stable contacts at the subunit interfaces estimated by a molecular dynamics simulation. We also calculate the hydration entropies of three different subcomplexes comprising the γ subunit, one of the β subunits, and two α subunits adjacent to them. The major finding is that the packing in F(1)-ATPase is highly asymmetrical, and this asymmetry is ascribed to the water entropy effect. We discuss how the rotation of the γ subunit is induced by such chemical processes as ATP binding, ATP hydrolysis, and release of the products. In our picture, the asymmetrical packing plays crucially important roles, and the rotation is driven by the water entropy effect.     

Single-molecule detection of DNA via sequence-specific links between F1-ATPase motors and gold nanorod sensors

We report the construction of a novel biosensing nanodevice to detect single, sequence-specific target DNA molecules. Nanodevice assembly occurs through the association of an immobilized F1-ATPase molecular motor and a functionalized gold nanorod via a single 3',5'-dibiotinylated DNA molecule. Target-dependent 3',5'-dibiotinylated DNA bridges form by combining ligation and exonucleation reactions (LXR), with a specificity capable of selecting against a single nucleotide polymorphism (SNP). Using dark field microscopy to detect gold nanorods, quantitation of assembled nanodevices is sufficient to distinguish the presence of as few as 1800 DNA bridges from nonspecifically bound nanorods. The rotary mechanism of F1-ATPase can drive gold nanorod rotation when the nanorod is attached via the DNA bridge. Therefore, rotation discriminates fully assembled devices from nonspecifically bound nanorods, resulting in a sensitivity limit of one zeptomole (600 molecules).     

Structural characteristics of yeast F1-ATPase before and after 16-degree rotation of the γ subunit: theoretical analysis focused on the water-entropy effect

We have recently proposed a novel picture of the rotation mechanism for F(1)-ATPase [T. Yoshidome, Y. Ito, M. Ikeguchi, and M. Kinoshita, J. Am. Chem. Soc. 133, 4030 (2011)]. In the picture, the asymmetric packing in F(1)-ATPase, originating from the water-entropy effect, plays the key role in the rotation. Here, we analyze the differences between the experimentally determined structures of yeast F(1)-ATPase before and after 16° rotation of the γ subunit with the emphasis on the water-entropy effect. For each of these structures, we calculate the hydration entropies of three sub-complexes comprising the γ subunit, one of the β subunits, and two α subunits adjacent to them. The β(E), β(TP), and β(DP) subunits are involved in sub-complexes I, II, and III, respectively. The calculation is performed using a hybrid of the angle-dependent integral equation theory combined with the molecular model for water and the morphometric approach. The absolute value of the hydration entropy is in the following order: sub-complex I > sub-complex II > sub-complex III. The packing efficiency of the sub-complex follows the opposite order. The rotation gives rise to less efficient packing in sub-complex III and a corresponding water-entropy loss. However, the other two sub-complexes, accompanying water-entropy gains, become more efficiently packed. These results are consistent with our picture of the rotation mechanism, supporting its validity. The water-entropy analysis shows that the interfaces of α(DP)-β(DP) and α(E)-β(E) become more open after the rotation, which is in accord with the experimental observation.     

Photochromic switchable pendant indolyl fulgimide polypyrrole

The photochromic switchable 2-indolyl fulgimide functionalized pyrrole monomer is synthesized and ascertained by spectroscopic techniques. The absorption and emission band of monomer due to intramolecular charge transfer (ICT) state exhibit typical bathochromic shift from nonpolar to polar solvents. The molecular geometry and bandgap energy of monomer in open form as well as closed form are calculated by B3LYP/6-31 G (d) with Gaussian 03 program package. The functionalized pyrrole monomer is electropolymerized by cyclicvoltammetric method on an electrode. The polymer indicated a completely reversible photoswitching between two isomers of open and closed form 2-indolyl fulgimide moiety, which are stable over many cycles of operation by UV/visible light sources. The closed form of polymer exhibit fluorescence quenching attributed to the energy transfer from conductance band of polypyrrole to closed form of 2-indolyl fulgimides unit. This fluorescence quenching system controls the conductance of conducting polymer in molecular electronics and molecular switches.     

Conformational transition pathway of polymerase beta/DNA upon binding correct incoming substrate

The closing conformational transition of wild-type polymerase beta bound to DNA template/primer before the chemical step (nucleotidyl transfer reaction) is simulated using the stochastic difference equation (in length version, "SDEL") algorithm that approximates long-time dynamics. The order of the events and the intermediate states during pol beta's closing pathway are identified and compared to a separate study of pol beta using transition path sampling (TPS) (Radhakrishnan, R.; Schlick, T. Proc. Natl. Acad. Sci. USA 2004, 101, 5970-5975). Results highlight the cooperative and subtle conformational changes in the pol beta active site upon binding the correct substrate that may help explain DNA replication and repair fidelity. These changes involve key residues that differentiate the open from the closed conformation (Asp192, Arg258, Phe272), as well as residues contacting the DNA template/primer strand near the active site (Tyr271, Arg283, Thr292, Tyr296) and residues contacting the beta and gamma phosphates of the incoming nucleotide (Ser180, Arg183, Gly189). This study compliments experimental observations by providing detailed atomistic views of the intermediates along the polymerase closing pathway and by suggesting additional key residues that regulate events prior to or during the chemical reaction. We also show general agreement between two sampling methods (the stochastic difference equation and transition path sampling) and identify methodological challenges involved in the former method relevant to large-scale biomolecular applications. Specifically, SDEL is very quick relative to TPS for obtaining an approximate path of medium resolution and providing qualitative information on the sequence of events; however, associated free energies are likely very costly to obtain because this will require both successful further refinement of the path segments close to the bottlenecks and large computational time.     

On the induced‐fit mechanism of substrate‐enzyme binding structures of nylon‐oligomer hydrolase

We present a detailed computational investigation of the induced‐fit motion in a nylon‐oligomer hydrolase (NylB) upon substrate binding. To this aim, we resort on the recently introduced parallel cascade selection molecular dynamics approach, allowing for an accelerated access to the set of conformational changes from an open‐ to a closed‐state structure to form the enzyme‐substrate complex in a specific induce‐fit mechanism. The structural investigation is quantitatively complemented by free energy analyses within the umbrella sampling algorithm accompanied by weighted histogram analysis. We find that the stabilization free energy is about 1.4 kcal/mol, whereas the highest free energy barrier to be overcome is about 2.3 kcal/mol. Conversely, the energetic contribution for the substrate binding is about 20 kcal/mol, as estimated from Generalized Born/Surface Area. This means that the open‐close induced‐fit motion could occur frequently once the substrate binds to the open state of NylB. © 2014 Wiley Periodicals, Inc.     

Dynamics of closure of zinc bis-porphyrin molecular tweezers with copper(II) ions and electron transfer

Zinc bis-porphyrin molecular tweezers composed of a N(4) spacer bound through pyridyl units to the meso position of porphyrins were synthesized, and the tweezers are closed by the coordination of a copper(II) ion inside the spacer ligand. The effect of the π-π interaction between the porphyrin rings in the closed conformation on the absorption spectra of multi-electron oxidized species and the reduction potentials were clarified by chemical and electrochemical oxidation of the closed form of the zinc bis-porphyrin molecular tweezers in comparison with the open form without copper(II) ion and the corresponding porphyrin monomer. The shifts in redox potentials and absorption spectrum of the porphyrin dication indicate a strong electronic interaction between the two oxidized porphyrins in the closed form, whereas there is little interaction between them in the neutral form. The dynamics of copper(II) ion coordination and subsequent electron transfer was examined by using a stopped-flow UV/Vis spectroscopic technique. It was confirmed that coordination of copper(II) occurs prior to electron-transfer oxidation of the closed form of the zinc bis-porphyrin molecular tweezers.     

Inherent flexibility and protein function: The open/closed conformational transition in the N-terminal domain of calmodulin

The key to understand a protein's function often lies in its conformational dynamics. We develop a coarse-grained variational model to investigate the interplay between structural transitions, conformational flexibility, and function of the N-terminal calmodulin domain (nCaM). In this model, two energy basins corresponding to the "closed" apo conformation and "open" holo conformation of nCaM are coupled by a uniform interpolation parameter. The resulting detailed transition route from our model is largely consistent with the recently proposed EFbeta-scaffold mechanism in EF-hand family proteins. We find that the N-terminal parts of the calcium binding loops shows higher flexibility than the C-terminal parts which form this EFbeta-scaffold structure. The structural transition of binding loops I and II are compared in detail. Our model predicts that binding loop II, with higher flexibility and earlier structural change than binding loop I, dominates the open/closed conformational transition in nCaM.     

Simulation of the opening and closing of Hsp70 chaperones by coarse-grained molecular dynamics

Heat-shock proteins 70 (Hsp70s) are key molecular chaperones which assist in the folding and refolding/disaggregation of proteins. Hsp70s, which consist of a nucleotide-binding domain (NBD, consisting of NBD-I and NBD-II subdomains) and a substrate-binding domain [SBD, further split into the β-sheet (SBD-β) and α-helical (SBD-α) subdomains], occur in two major conformations having (a) a closed SBD, in which the SBD and NBD domains do not interact, (b) an open SBD, in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II. In the SBD-closed conformation, SBD is bound to a substrate protein, with release occurring after transition to the open conformation. While the transition from the closed to the open conformation is triggered efficiently by binding of adenosine triphosphate (ATP) to the NBD, it also occurs, although less frequently, in the absence of ATP. The reverse transition occurs after ATP hydrolysis. Here, we report canonical and multiplexed replica exchange simulations of the conformational dynamics of Hsp70s using a coarse-grained molecular dynamics approach with the UNRES force field. The simulations were run in the following three modes: (i) with the two halves of the NBD unrestrained relative to each other, (ii) with the two halves of the NBD restrained in an "open" geometry as in the SBD-closed form of DnaK (2KHO), and (iii) the two halves of NBD restrained in a "closed" geometry as in known experimental structures of ATP-bound NBD forms of Hsp70. Open conformations, in which the SBD interacted strongly with the NBD, formed spontaneously during all simulations; the number of transitions was largest in simulations carried out with the "closed" NBD domain, and smallest in those carried out with the "open" NBD domain; this observation is in agreement with the experimentally-observed influence of ATP-binding on the transition of Hsp70's from the SBD-closed to the SBD-open form. Two kinds of open conformations were observed: one in which SBD-α interacts with NBD-I and SBD-β interacts with the top parts of NBD-I and NBD-II (as observed in the structures of nucleotide exchange factors), and another one in which this interaction pattern is swapped. A third type of motion, in which SBD-α binds to NBD without dissociating from SBD-β was also observed. It was found that the first stage of interdomain communication (approach of SBD-β, to NBD) is coupled with the rotation of the long axes of NBD-I and NBD-II towards each other. To the best of our knowledge, this is the first successful simulation of the full transition of an Hsp70 from the SBD-closed to the SBD-open conformation.     

Tryptophan synthase: a multienzyme complex with an intramolecular tunnel

Tryptophan synthase is a classic enzyme that channels a metabolic intermediate, indole. The crystal structure of the tryptophan synthase alpha2beta2 complex from Salmonella typhimurium revealed for the first time the architecture of a multienzyme complex and the presence of an intramolecular tunnel. This remarkable hydrophobic tunnel provides a likely passageway for indole from the active site of the alpha subunit, where it is produced, to the active site of the beta subunit, where it reacts with L-serine to form L-tryptophan in a pyridoxal phosphate-dependent reaction. Rapid kinetic studies of the wild type enzyme and of channel-impaired mutant enzymes provide strong evidence for the proposed channeling mechanism. Structures of a series of enzyme-substrate intermediates at the alpha and beta active sites are elucidating enzyme mechanisms and dynamics. These structural results are providing a fascinating picture of loops opening and closing, of domain movements, and of conformational changes in the indole tunnel. Solution studies provide further evidence for ligand-induced conformational changes that send signals between the alpha and beta subunits. The combined results show that the switching of the enzyme between open and closed conformations couples the catalytic reactions at the alpha and beta active sites and prevents the escape of indole.     

Single-molecule Förster resonance energy transfer reveals an innate fidelity checkpoint in DNA polymerase I

Enzymatic reactions typically involve complex dynamics during substrate binding, conformational rearrangement, chemistry, and product release. The noncovalent steps provide kinetic checkpoints that contribute to the overall specificity of enzymatic reactions. DNA polymerases perform DNA replication with outstanding fidelity by actively rejecting noncognate nucleotide substrates early in the reaction pathway. Substrates are delivered to the active site by a flexible fingers subdomain of the enzyme, as it converts from an open to a closed conformation. The conformational dynamics of the fingers subdomain might also play a role in nucleotide selection, although the precise role is currently unknown. Using single-molecule F?rster resonance energy transfer, we observed individual Escherichia coli DNA polymerase I (Klenow fragment) molecules performing substrate selection. We discovered that the fingers subdomain actually samples through three distinct conformations--open, closed, and a previously unrecognized intermediate conformation. We measured the overall dissociation rate of the polymerase-DNA complex and the distribution among the various conformational states in the absence and presence of nucleotide substrates, which were either correct or incorrect. Correct substrates promote rapid progression of the polymerase to the catalytically competent closed conformation, whereas incorrect nucleotides block the enzyme in the intermediate conformation and induce rapid dissociation from DNA. Remarkably, incorrect nucleotide substrates also promote partitioning of DNA to the spatially separated 3'-5' exonuclease domain, providing an additional mechanism to prevent misincorporation at the polymerase active site. These results reveal the existence of an early innate fidelity checkpoint, rejecting incorrect nucleotide substrates before the enzyme encloses the nascent base pair.     

In silico models for the human alpha4beta2 nicotinic acetylcholine receptor

The neuronal alpha4beta2 nicotinic acetylcholine receptor (nAChR) is one of the most widely expressed nAChR subtypes in the brain. Its subunits have high sequence identity (54 and 46% for alpha4 and beta2, respectively) with alpha and beta subunits in Torpedo nAChR. Using the known structure of the Torpedo nAChR as a template, the closed-channel structure of the alpha4beta2 nAChR was constructed through homology modeling. Normal-mode analysis was performed on this closed structure and the resulting lowest frequency mode was applied to it for a "twist-to-open" motion, which increased the minimum pore radius from 2.7 to 3.4 A and generated an open-channel model. Nicotine could bind to the predicted agonist binding sites in the open-channel model but not in the closed one. Both models were subsequently equilibrated in a ternary lipid mixture via extensive molecular dynamics (MD) simulations. Over the course of 11 ns MD simulations, the open channel remained open with filled water, but the closed channel showed a much lower water density at its hydrophobic gate comprised of residues alpha4-V259 and alpha4-L263 and their homologous residues in the beta2 subunits. Brownian dynamics simulations of Na+ permeation through the open channel demonstrated a current-voltage relationship that was consistent with experimental data on the conducting state of alpha4beta2 nAChR. Besides establishment of the well-equilibrated closed- and open-channel alpha4beta2 structural models, the MD simulations on these models provided valuable insights into critical factors that potentially modulate channel gating. Rotation and tilting of TM2 helices led to changes in orientations of pore-lining residue side chains. Without concerted movement, the reorientation of one or two hydrophobic side chains could be enough for channel opening. The closed- and open-channel structures exhibited distinct patterns of electrostatic interactions at the interface of extracellular and transmembrane domains that might regulate the signal propagation of agonist binding to channel opening. A potential prominent role of the beta2 subunit in channel gating was also elucidated in the study.     

Molecular dynamics simulations of ligand-induced backbone conformational changes in the binding site of the periplasmic lysine-, arginine-, ornithine-binding protein

The periplasmic lysine-, arginine-, ornithine-binding protein (LAOBP) traps its ligands by a large hinge bending movement between two globular domains. The overall geometry of the binding site remains largely unchanged between the open (unliganded) and closed (liganded) forms, with only a small number of residues exhibiting limited movement of their side chains. However, in the case of the ornithine-bound structure, the backbone peptide bond between Asp11 and Thr12 undergoes a large rotation. Molecular dynamics simulations have been used to investigate the origin and mechanism of this backbone movement. Simulations allowing flexibility of a limited region and of the whole binding site, with and without bound ligands, suggest that this conformational change is induced by the binding of ornithine, leading to the stabilisation of an energetically favourable alternative conformation.     

Unraveling Dual Aggregation‐Induced Emission Behavior in Steric‐Hindrance Photochromic System for Super Resolution Imaging

Unprecedented dual aggregation‐induced emission (AIE) behavior based on a steric‐hindrance photochromic system is presented, with incorporation one or two bulky aryl groups, resulting in different flexibleness. The dual AIE behavior of open and closed isomers can be explained by restriction of intramolecular rotation (RIR), restriction of intramolecular vibration (RIV), and intermolecular stacking. The large bulky benzothiophene causes restricted rotation, enhancing the emission of open form in solution and weak π–π molecular packing, thereby efficiently enhancing the luminescence performance in the solid state. With incorporation of two large bulky benzothiophene groups, BBTE possesses the most outstanding AIE activity, undergoing highly efficient and reversible off‐to‐on fluorescence in film upon alternating UV and visible light irradiation along with excellent fatigue resistance. The off‐to‐on fluorescent photoswitch is successfully established in super resolution imaging.     

Inhibition of the mitochondrial F1-ATPase by rose bengal mediated photooxidation. Interaction of the Fe2+ chelate of bathophenanthroline with the sensitizer

Rose Bengal mediated photooxidation of mitochondrial F1-ATPase and its beta-subunit resulted in inactivation and loss of about 50 and 60% of their histidine residues, respectively. The beta-subunit was not cleaved upon photooxidation. Photooxidation of histidine probably results in changes in the conformational stability of F1-ATPase leading to its inactivation. The participation of singlet molecular oxygen during the photooxidation process is suggested by the selective loss of histidine residues, while other amino acids, also sensitive to singlet oxygen attack, were not affected. Photochemical damage of F1-ATPase was prevented by various phenanthroline compounds, the order of efficiency being bathophenanthroline-Fe chelate greater than bathophenanthroline greater than orthophenanthroline-Fe chelate greater than bathophenanthroline-sulfonate-Fe chelate. The prevention by bathophenanthroline-Fe chelate of photochemical damage is interpreted on the basis of its interaction with the photosensitizer, Rose Bengal, probably implying a chemical reaction which decreases the actual concentration of the sensitizer and, thereby, the extent of photoinactivation.     

The large scale conformational change of the human DPP III-substrate prefers the "closed" form

Human dipeptidyl peptidase III (DPP III) is a two domain metallo-peptidase from the M49 family. The wide interdomain cleft and broad substrate specificity suggest that this enzyme could experience significant conformational change. Long (>100 ns) molecular dynamics (MD) simulations of DPP III revealed large range conformational changes of the protein, suggesting the pre-existing equilibrium model for a substrate binding. The binding free energy calculations revealed tighter binding of the preferred synthetic substrate Arg-Arg-2-naphtylamide to the "closed" than to the "open" DPP III conformation. Our assumption that Asp372 plays a crucial role in the large scale interdomain closure was proved by the MD simulations of the Asp372Ala variant. During the same simulation time, the variant remained more "open" than the wild type protein. Apparently, Ala was not as efficient as Asp in establishing the interdomain interactions. According to the MM-PBSA calculations, the electrostatic component of the free energy of solvation turned out to be higher for the "closed" protein than for its less compact form. However, the gain in entropy due to water released from the interdomain cleft nicely balanced this negative effect.