驱动蛋白结构与运动机制

驱动蛋白是一类能够利用ATP水解释放的化学能驱动其所携带的“货物”分子沿着微管（microtubule,MT）定向运动的分子马达,在细胞器运输、有丝分裂、轴突运输等方面有着重要的生理作用。随着驱动蛋白结合ADP、ATP和未结合核苷酸（APO）三种特征状态的晶体结构的解析,驱动蛋白构象变化的研究得到了进一步发展,而在力产生机制和运动模型方面仍然存在较大争议。本文以kinesin-1家族为例,分析了驱动蛋白三种特征状态结构的特点、状态结构间的构象转变,论述了驱动蛋白的力产生机制和整个迈步过程。并探讨了驱动蛋白的运动模型,同时采用分子动力学模拟比较了驱动蛋白的两种迈步方式,为深入研究驱动蛋白提供了一定的理论计算。最后,基于本课题组对复杂体系的研究,对驱动蛋白体系的控制机制提出了新的假设,并对未来的研究方向进行了展望。     

Toward a detailed description of the pathways of allosteric communication in the GroEL chaperonin through atomistic simulation

GroEL, along with its coprotein GroES, is essential for ensuring the correct folding of unfolded or newly synthesized proteins in bacteria. GroEL is a complex, allosteric molecule, composed of two heptameric rings stacked back to back, that undergoes large structural changes during its reaction cycle. These structural changes are driven by the cooperative binding and subsequent hydrolysis of ATP, by GroEL. Despite numerous previous studies, the precise mechanisms of allosteric communication and the associated structural changes remain elusive. In this paper, we describe a series of all-atom, unbiased, molecular dynamics simulations over relatively long (50-100 ns) time scales of a single, isolated GroEL subunit and also a heptameric GroEL ring, in the presence and absence of ATP. Combined with results from a distance restraint-biased simulation of the single ring, the atomistic details of the earliest stages of ATP-driven structural changes within this complex molecule are illuminated. Our results are in broad agreement with previous modeling studies of isolated subunits and with a coarse-grained, forcing simulation of the single ring. These are the first reported all-atom simulations of the GroEL single-ring complex and provide a unique insight into the role of charged residues K80, K277, R284, R285, and E388 at the subunit interface in transmission of the allosteric signal. These simulations also demonstrate the feasibility of performing all-atom simulations of very large systems on sufficiently long time scales on typical high performance computing facilities to show the origins of the earliest events in biologically relevant processes.     

The rotor in the membrane of the ATP synthase and relatives

In recent years, structural information on the F(1) sector of the ATP synthase has provided an insight into the molecular mechanism of ATP catalysis. The structure strongly supports the proposal that the ATP synthase works as a rotary molecular motor. Insights into the membrane domain have just started to emerge but more detailed structural information is needed if the molecular mechanism of proton translocation coupled to ATP synthesis is to be understood. This review will focus mainly on the ion translocating rotor in the membrane domain of the F-type ATPase, and the related vacuolar and archaeal relatives.     

K and Na Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

K+ and Na+ Conduction in Selective and Nonselective Ion Channels Via Molecular Dynamics Simulations

Generations of scientists have been captivated by ion channels and how they control the workings of the cell by admitting ions from one side of the cell membrane to the other. Elucidating the molecular determinants of ion conduction and selectivity are two of the most fundamental issues in the field of biophysics. Combined with ongoing progress in structural studies, modeling and simulation have been an integral part of the development of the field. As of this writing, the relentless growth in computational power, the development of new algorithms to tackle the so-called rare events, improved force-field parameters, and the concomitant increasing availability of membrane protein structures, allow simulations to contribute even further, providing more-complete models of ion conduction and selectivity in ion channels. In this report, we give an overview of the recent progress made by simulation studies on the understanding of ion permeation in selective and nonselective ion channels.     

Molecular dynamics simulation of ionic conductors: perspectives and limitations

Concepts of modeling ion conductors by means of molecular dynamics simulations are outlined and exemplified by a small series of applications. Computational studies offer atomic level of detail insights into ionic mobility and thus help to rationalize experimental findings in particular for complex ion conduction mechanisms. Standard techniques are discussed within their thermodynamic context, i.e., chemical potentials, conductivity and activation energies. The limitations of molecular dynamics simulations are defined by statistics and the accuracy in evaluating atomic interactions.     

A combined coarse-grained and all-atom simulation of TRPV1 channel gating and heat activation

The transient receptor potential (TRP) channels act as key sensors of various chemical and physical stimuli in eukaryotic cells. Despite years of study, the molecular mechanisms of TRP channel activation remain unclear. To elucidate the structural, dynamic, and energetic basis of gating in TRPV1 (a founding member of the TRPV subfamily), we performed coarse-grained modeling and all-atom molecular dynamics (MD) simulation based on the recently solved high resolution structures of the open and closed form of TRPV1. Our coarse-grained normal mode analysis captures two key modes of collective motions involved in the TRPV1 gating transition, featuring a quaternary twist motion of the transmembrane domains (TMDs) relative to the intracellular domains (ICDs). Our transition pathway modeling predicts a sequence of structural movements that propagate from the ICDs to the TMDs via key interface domains (including the membrane proximal domain and the C-terminal domain), leading to sequential opening of the selectivity filter followed by the lower gate in the channel pore (confirmed by modeling conformational changes induced by the activation of ICDs). The above findings of coarse-grained modeling are robust to perturbation by lipids. Finally, our MD simulation of the ICD identifies key residues that contribute differently to the nonpolar energy of the open and closed state, and these residues are predicted to control the temperature sensitivity of TRPV1 gating. These computational predictions offer new insights to the mechanism for heat activation of TRPV1 gating, and will guide our future electrophysiology and mutagenesis studies.     

Tightening of the ATP-binding sites induces the opening of P2X receptor channels

The opening of ligand-gated ion channels in response to agonist binding is a fundamental process in biology. In ATP-gated P2X receptors, little is known about the molecular events that couple ATP binding to channel opening. In this paper, we identify structural changes of the ATP site accompanying the P2X2 receptor activation by engineering extracellular zinc bridges at putative mobile regions as revealed by normal mode analysis. We provide evidence that tightening of the ATP sites shaped like open 'jaws' induces opening of the P2X ion channel. We show that ATP binding favours jaw tightening, whereas binding of a competitive antagonist prevents gating induced by this movement. Our data reveal the inherent dynamic of the binding jaw, and provide new structural insights into the mechanism of P2X receptor activation.     

Molecular and structural insight into plasmodium falciparum RIO2 kinase

Among approximately 65 kinases of the malarial genome, RIO2 (right open reading frame) kinase belonging to the atypical class of kinase is unique because along with a kinase domain, it has a highly conserved N-terminal winged helix (wHTH) domain. The wHTH domain resembles the wing like domain found in DNA binding proteins and is situated near to the kinase domain. Ligand binding to this domain may reposition the kinase domain leading to inhibition of enzyme function and could be utilized as a novel allosteric site to design inhibitor. In the present study, we have generated a model of RIO2 kinase from Plasmodium falciparum utilizing multiple modeling, simulation approach. A novel putative DNA-binding site is identified for the first time in PfRIO2 kinase to understand the DNA binding events involving wHTH domain and flexible loop. Induced fit DNA docking followed by minimization, molecular dynamics simulation, energetic scoring and binding mode studies are used to reveal the structural basis of PfRIO2-ATP-DNA complex. Ser105 as a potential site of phosphorylation is revealed through the structural studies of ATP binding in PfRIO2. Overall the present study discloses the structural facets of unknown PfRIO2 complex and opens an avenue toward exploration of novel drug target.     

Halting ionic shuttle to disrupt the synthetic machinery—Structural and molecular insights into the inhibitory roles of Bedaquiline towards Mycobacterium tuberculosis ATP synthase in the treatment of tuberculosis

Therapeutic targeting of the adenosine triphosphate (ATP) machinery of Mycobacterium tuberculosis (Mtb) has recently presented a potent and alternative measure to halt the pathogenesis of tuberculosis. This has been potentiated by the development of bedaquiline (BDQ), a novel small molecule inhibitor that selectively inhibits mycobacterial F₁F_o-ATP synthase by targeting its rotor c-ring, resulting in the disruption of ATP synthesis and consequential cell death. Although the structural resolution of the mycobacterial C₉ ring in co`mplex with BDQ provided the first-hand detail of BDQ interaction at the c-ring region of the ATP synthase, there still remains a need to obtain essential and dynamic insights into the mechanistic activity of this drug molecule towards crucial survival machinery of Mtb. As such, for the first time, we report an atomistic model to describe the structural dynamics that explicate the experimentally reported antagonistic features of BDQ in halting ion shuttling by the mycobacterial c-ring, using molecular dynamics simulation and the Molecular Mechanics/Poisson-Boltzmann Surface Area methods. Results showed that BDQ exhibited a considerably high ΔG while it specifically maintained high-affinity interactions with Glu65_B and Asp32_B, blocking their crucial roles in proton binding and shuttling, which is required for ATP synthesis. Moreover, the bulky nature of BDQ induced a rigid and compact conformation of the rotor c-ring, which impedes the essential rotatory motion that drives ion exchange and shuttling. In addition, the binding affinity of a BDQ molecule was considerably increased by the complementary binding of another BDQ molecule, which indicates that an increase in BDQ molecule enhances inhibitory potency against Mtb ATP synthase. Taken together, findings provide atomistic perspectives into the inhibitory mechanisms of BDQ coupled with insights that could enhance the structure-based design of novel ATP synthase inhibitors towards the treatment of tuberculosis.     

Molecular model and ATPase activity of carboxyl-terminal nucleotide binding domain from human P-glycoprotein

ATP binding and hydrolysis are required for P-glycoprotein mediated multidrug resistance. To investigate the molecular mechanism involved in ATP binding and hydrolysis, a three-dimensional model of the carboxyl-terminal nucleotide binding domain (NBD2) was built by homology modeling. Modeling revealed the human P-glycoprotein ATP-binding site and the possible role of conserved Gln1118 residue. Recombinant NBD2 was overexpressed in Escherichia coli and the conserved Gln1118 residue was mutated to an alanine residue. The Vmax for ATP hydrolysis by the mutant NBD2 was approximately 56% of the Vmax of wild-type NBD2. But both proteins displayed similar affinity for ATP, with Km of 479 and 466 microM for mutant and wild-type NBD2, respectively. These results suggest that the possible role of Gln1118 is as an activating residue for ATP hydrolysis. The molecular model also provided structural information about the interactions between NBD2 and the chemosensitizer quercetin. The complex indicated that quercetin was tightly bound to the ATP-binding site and competed for binding. The three-dimensional model of NBD2 can be used to both guide enzymological studies and provide a theoretical basis for the design of potential multidrug resistance reversers.     

Hierarchical approach to predicting permeation in ion channels

A hierarchical computational strategy combining molecular modeling, electrostatics calculations, molecular dynamics, and Brownian dynamics simulations is developed and implemented to compute electrophysiologically measurable properties of the KcsA potassium channel. Models for a series of channels with different pore sizes are developed from the known x-ray structure, using insights into the gating conformational changes as suggested by a variety of published experiments. Information on the pH dependence of the channel gating is incorporated into the calculation of potential profiles for K(+) ions inside the channel, which are then combined with K(+) ion mobilities inside the channel, as computed by molecular dynamics simulations, to provide inputs into Brownian dynamics simulations for computing ion fluxes. The open model structure has a conductance of approximately 110 pS under symmetric 250 mM K(+) conditions, in reasonable agreement with experiments for the largest conducting substate. The dimensions of this channel are consistent with electrophysiologically determined size dependence of quaternary ammonium ion blocking from the intracellular end of this channel as well as with direct structural evidence that tetrabutylammonium ions can enter into the interior cavity of the channel. Realistic values of Ussing flux ratio exponents, distribution of ions within the channel, and shapes of the current-voltage and current-concentration curves are obtained. The Brownian dynamics calculations suggest passage of ions through the selectivity filter proceeds by a "knock-off" mechanism involving three ions, as has been previously inferred from functional and structural studies of barium ion blocking. These results suggest that the present calculations capture the essential nature of K(+) ion permeation in the KcsA channel and provide a proof-of-concept for the integrated microscopic/mesoscopic multitiered approach for predicting ion channel function from structure, which can be applied to other channel structures.     

Inhibition of skeletal muscle ClC-1 chloride channels by low intracellular pH and ATP

Skeletal muscle acidosis during exercise has long been thought to be a cause of fatigue, but recent studies have shown that acidosis maintains muscle excitability and opposes fatigue by decreasing the sarcolemmal chloride conductance. ClC-1 is the primary sarcolemmal chloride channel and has a clear role in controlling muscle excitability, but recombinant ClC-1 has been reported to be activated by acidosis. Following our recent finding that intracellular ATP inhibits ClC-1, we investigated here the interaction between pH and ATP regulation of ClC-1. We found that, in the absence of ATP, intracellular acidosis from pH 7.2 to 6.2 inhibited ClC-1 slightly by shifting the voltage dependence of common gating to more positive potentials, similar to the effect of ATP. Importantly, the effects of ATP and acidosis were cooperative, such that ATP greatly potentiated the effect of acidosis. Adenosine had a similar effect to ATP at pH 7.2, but acidosis did not potentiate this effect, indicating that the phosphates of ATP are important for this cooperativity, possibly due to electrostatic interactions with protonatable residues of ClC-1. A protonatable residue identified by molecular modeling, His-847, was found to be critical for both pH and ATP modulation and may be involved in such electrostatic interactions. These findings are now consistent with, and provide a molecular explanation for, acidosis opposing fatigue by decreasing the chloride conductance of skeletal muscle via inhibition of ClC-1. The modulation of ClC-1 by ATP is a key component of this molecular mechanism.     

The Contribution of the Ankyrin Repeat Domain of TRPV1 as a Thermal Module

The TRPV1 cation nonselective ion channel plays an essential role in thermosensation and perception of other noxious stimuli. TRPV1 can be activated by low extracellular pH, high temperature, or naturally occurring pungent molecules such as allicin, capsaicin, or resiniferatoxin. Its noxious thermal sensitivity makes it an important participant as a thermal sensor in mammals. However, details of the mechanism of channel activation by increases in temperature remain unclear. Here, we used a combination of approaches to try to understand the role of the ankyrin repeat domain (ARD) in channel behavior. First, a computational modeling approach by coarse-grained molecular dynamics simulation of the whole TRPV1 embedded in a phosphatidylcholine and phosphatidylethanolamine membrane provides insight into the dynamics of this channel domain. Global analysis of the structural ensemble shows that the ARD is a region that sustains high fluctuations during dynamics at different temperatures. We then performed biochemical and thermal stability studies of the purified ARD by the means of circular dichroism and tryptophan fluorescence and demonstrate that this region undergoes structural changes at similar temperatures that lead to TRPV1 activation. Our data suggest that the ARD is a dynamic module and that it may participate in controlling the temperature sensitivity of TRPV1.     

Magnesium magnetic isotope effect: A key to the mechanochemistry of phosphorylating enzymes as molecular machines

The discovery of the powerful magnesium isotope effect on enzymatic ATP synthesis provides a new insight into the mechanochemistry of enzymes as molecular machines. The catalytic activities of ATPase, creatine kinase, and glycerophsphate kinase containing a Mg²⁺ ion with magnetic isotope nuclei (²⁵Mg) were found to be two to four times higher than those of the enzymes with spinless, nonmagnetic magnesium cation isotopes (²⁴Mg or ²⁶Mg). This demonstrates unambiguously that ATP synthesis is a spin-selective process involving Mg2+ as the electron-accepting reagent. ATP synthesis proceeds in an ion-radical pair consisting of an ADP oxyradical and Mg²⁺. In this process, the magnesium bivalent cation is the key agent that transforms the mechanics of a protein molecule into chemical processes. This ion is the crucial structural component of enzymes as mechanochemical molecular machines.     

Homology modeling,molecular dynamic simulation and docking studies of cyclin dependent kinase 1

In order to develop promising cyclin dependent kinase 1 inhibitors, homology modeling, docking and molecular dynamic simulation techniques were applied to get insight into the functional and structural properties of cyclin dependent kinase 1 (CDK1). Since there is no reported CDK1 crystal structural data, the three dimensional structure of CDK1 was constructed based on homology modeling. An extensive dynamic simulation was also performed on a Flavopiridol-CDK1 complex for probing the binding pattern of Flavopiridol in the active site of CDK1. The binding modes of other inhibitors to CDK1 were also proposed by molecular docking. The structural requirement for developing more potent CDK1 inhibitors was obtained by the above-mentioned molecular simulations and pharmacophore modeling.     

Molecular structure of low density lipoprotein: current status and future challenges

This review highlights recent advances in structural studies on low density lipoprotein (LDL) with particular emphasis on the apolipoprotein moiety of LDL, apolipoprotein B100 (apoB100). Various molecular aspects of LDL are outlined and obstacles to structure determination are addressed. In this context, the prevailing conceptions of the molecular assembly of LDL and how the synergy of complementary biochemical, biophysical and molecular simulation approaches has lead to the current structural model of LDL are discussed. Evidence is presented that structural heterogeneity and the intrinsic dynamics of LDL are key determinants of the functionality of LDL in both health and disease. Some key research directions, remaining open questions and rapidly emerging new concepts for medical applications of LDL, are furthermore outlined. The article concludes by providing an outlook concerning promising future strategies for the clarification of the molecular details of LDL, in particular of apoB100, combining recent advances in molecular modeling with developments of novel experimental techniques. Although new insights into the molecular organization of LDL are forthcoming, many open questions remain unanswered. The major challenge of the next decade will certainly be the elucidation of the molecular structural and dynamic features of apoB100.     

Structural insights into IKKβ inhibition by natural products staurosporine and quercetin

This Letter describes the results of two combined approaches: homology modeling and molecular docking studies, in order to propose the molecular basis of IKKβ inhibition by staurosporine and quercetin as ATP-competitive inhibitors. The results provides a rationale and structural frameworks for designing potent ATP binding-site inhibitors of IKKβ, which is an attractive drug target for inflammatory diseases and has been found to be responsible for some of the already observed pharmacological effects for marketed drugs.     

Insights into structural mechanisms of gating induced regulation of aquaporins

Aquaporin family comprises of transmembrane channels that are specialized in conducting water and certain small, uncharged molecules across cell membranes. Essential roles of aquaporins in various physiological and pathophysiological conditions have attracted great scientific interest. Pioneering structural studies on aquaporins have almost solved the basic question of mechanism of selective water transport through these channels. Another important structural aspect of aquaporins which seeks attention is that how the flow of water through the channel is regulated by the mechanism of gating. Aquaporins are also regulated at the protein level, i.e. by trafficking which includes changes in their expression levels in the membrane. Availability of high resolution structures along with numerous molecular dynamics simulation studies have helped to gain an understanding of the structural mechanisms by which water flux through aquaporins is controlled. This review will summarize the highlights regarding structural features of aquaporins, mechanisms governing water permeation, proton exclusion and substrate specificity, and describe the structural insights into the mechanisms of aquaporin gating whereby water conduction is regulated by post translational modifications, such as phosphorylation.     

Distinct sensitivities of OmpF and PhoE porins to charged modulators

The inhibition of the anion-selective PhoE porin by ATP and of the cation-selective OmpF porin by polyamines has been previously documented. In the present study, we have extended the comparison of the inhibitor-porin pairs by investigating the effect of anions (ATP and aspartate) and positively charged polyamines (spermine and cadaverine) on both OmpF and PhoE with the patch-clamp technique, and by comparing directly the gating kinetics of the channels modulated by their respective substrates. The novel findings reported here are (1) that the activity of PhoE is completely unaffected by polyamines, and (2) that the kinetic changes induced by ATP on PhoE or polyamines on OmpF suggest different mechanisms of inhibition. ATP induces a high degree of flickering in the PhoE-mediated current and appears to behave as a blocker of ion flow during its presumed transport through PhoE. Polyamines modulate the kinetics of openings and closings of OmpF, in addition to promoting a blocker-like flickering activity. The strong correlation between sensitivity to inhibitors and ion selectivity suggests that some common molecular determinants are involved in these two properties and is in agreement with the hypothesis that polyamines bind inside the pore of cationic porins.