Multivariate Analysis of Conserved Sequence–Structure Relationships in Kinesins: Coupling of the Active Site and a Tubulin-binding Sub-domain

Wordom update 2: A user-friendly program for the analysis of molecular structures and conformational ensembles

We present the second update of Wordom, a user-friendly and efficient program for manipulation and analysis of conformational ensembles from molecular simulations. The actual update expands some of the existing modules and adds 21 new modules to the update 1 published in 2011. The new adds can be divided into three sets that: 1) analyze atomic fluctuations and structural communication; 2) explore ion-channel conformational dynamics and ionic translocation; and 3) compute geometrical indices of structural deformation. Set 1 serves to compute correlations of motions, find geometrically stable domains, identify a dynamically invariant core, find changes in domain-domain separation and mutual orientation, perform wavelet analysis of large-scale simulations, process the output of principal component analysis of atomic fluctuations, perform functional mode analysis, infer regions of mechanical rigidity, analyze overall fluctuations, and perform the perturbation response scanning. Set 2 includes modules specific for ion channels, which serve to monitor the pore radius as well as water or ion fluxes, and measure functional collective motions like receptor twisting or tilting angles. Finally, set 3 includes tools to monitor structural deformations by computing angles, perimeter, area, volume, β-sheet curvature, radial distribution function, and center of mass. The ring perception module is also included, helpful to monitor supramolecular self-assemblies. This update places Wordom among the most suitable, complete, user-friendly, and efficient software for the analysis of biomolecular simulations. The source code of Wordom and the relative documentation are available under the GNU general public license at http://wordom.sf.net.

A Potential Route of Capsaicin to Its Binding Site in the TRPV1 Ion Channel

Transient receptor potential (TRP) ion channels are important pharmacological targets because of their role in the perception of pain, and so, understanding their chemical regulation is essential for the development of analgesic drugs. Among the currently known TRP channel chemical agonists, capsaicin, the active compound of chili pepper, is probably the most exhaustively studied. The availability of the three-dimensional structure of the vanilloid receptor 1 (TRPV1) has fueled computational studies revealing the molecular details of capsaicin binding modes. Although this is a significant step, a comprehensible binding mechanism or pathway is invaluable for targeting TRP channels in modern pharmacology. In the present work, free-energy and enhanced sampling techniques have been used to explore a possible membrane-mediated pathway for capsaicin to enter the TRPV1 binding pocket where capsaicin accesses the protein starting at the extracellular milieu through the outer leaflet and into its binding site in the protein. The main states visited along this route have been characterized and include (i) a bound state in agreement with the binding mode "head-down, tail-up" and (ii) an alternative state corresponding to a "head-up, tail-down" binding mode. In agreement with previous reports, binding is mediated by both hydrogen bonds and van der Waals interactions, and residue Y511 is crucial for stabilizing the bound state and during the binding process. Together, these results provide a foundation to further understand TRPV channels, and they could be used to guide therapeutic design of selective inhibitors potentially leading to novel avenues for pharmacological applications targeting the TRPV1 channel.

Scaffold-based selective SHP2 inhibitors design using core hopping, molecular docking, biological evaluation and molecular simulation

PTPN11 (coding the gene of SHP2), a classic non-receptor protein tyrosine phosphatase, is implicated in multiple cell signaling pathway. Abnormal activation of SHP2 has been shown to contribute to a variety of human diseases, including Juvenile myelomonocytic leukemia (JMML), Noonan syndrome and tumors. Thus, the SHP2 inhibitors have important therapeutic value. Here, based on the compound PubChem CID 8,478,960 (IC₅₀ = 45.01 μM), a series of thiophene [2,3-d] pyrimidine derivatives (IC₅₀ = 0.4-37.87 μM) were discovered as novel and efficient inhibitors of SHP2 through powerful "core hopping" and CDOCKER technology. Furthermore, the SHP2-PTP phosphatase activity assay indicated that Comp#5 (IC₅₀ = 0.4 μM) was the most active SHP2 inhibitor. Subsequently, the effects of Comp#5 on the structure and function of SHP2 were investigated through molecular dynamics (MD) simulation and post-kinetic analysis. The result indicated that Comp#5 enhanced the interaction of residues THR357, ARG362, LYS366, PRO424, CYS459, SER460, ALA461, ILE463, ARG465, THR507 and GLN510 with the surrounding residues, improving the stability of the catalytic active region and the entrance of catalytic active region. In particular, the Comp#5 conjugated with residue ARG362, elevating the efficient and selectivity of SHP2 protein. The study here may pave the way for discovering the novel SHP2 inhibitors for suffering cancer patients.

The Bio3D packages for structural bioinformatics

Bio3D is a family of R packages for the analysis of biomolecular sequence, structure, and dynamics. Major functionality includes biomolecular database searching and retrieval, sequence and structure conservation analysis, ensemble normal mode analysis, protein structure and correlation network analysis, principal component, and related multivariate analysis methods. Here, we review recent package developments, including a new underlying segregation into separate packages for distinct analysis, and introduce a new method for structure analysis named ensemble difference distance matrix analysis (eDDM). The eDDM approach calculates and compares atomic distance matrices across large sets of homologous atomic structures to help identify the residue wise determinants underlying specific functional processes. An eDDM workflow is detailed along with an example application to a large protein family. As a new member of the Bio3D family, the Bio3D-eddm package supports both experimental and theoretical simulation-generated structures, is integrated with other methods for dissecting sequence-structure-function relationships, and can be used in a highly automated and reproducible manner. Bio3D is distributed as an integrated set of platform independent open source R packages available from: http://thegrantlab.org/bio3d/.

Navigating the conformational landscape of G protein-coupled receptor kinases during allosteric activation

G protein-coupled receptors (GPCRs) are essential for transferring extracellular signals into carefully choreographed intracellular responses controlling diverse aspects of cell physiology. The duration of GPCR-mediated signaling is primarily regulated via GPCR kinase (GRK)-mediated phosphorylation of activated receptors. Although many GRK structures have been reported, the mechanisms underlying GRK activation are not well-understood, in part because it is unknown how these structures map to the conformational landscape available to this enzyme family. Unlike most other AGC kinases, GRKs rely on their interaction with GPCRs for activation and not phosphorylation. Here, we used principal component analysis of available GRK and protein kinase A crystal structures to identify their dominant domain motions and to provide a framework that helps evaluate how close each GRK structure is to being a catalytically competent state. Our results indicated that disruption of an interface formed between the large lobe of the kinase domain and the regulator of G protein signaling homology domain (RHD) is highly correlated with establishment of the active conformation. By introducing point mutations in the GRK5 RHD-kinase domain interface, we show with both in silico and in vitro experiments that perturbation of this interface leads to higher phosphorylation activity. Navigation of the conformational landscape defined by this bioinformatics-based study is likely common to all GPCR-activated GRKs.

Structure-Based Design of Scaffolds Targeting PDE10A by INPHARMA-NMR

Phosphodiesterases (PDE) hydrolyze both cyclic AMP and GMP (cAMP/cGMP) and are responsible for the regulation of their levels in a multitude of cellular functions. PDE10A is expressed in the brain and is a validated target for both schizophrenia and Huntington disease. Here, we address the identification of novel chemical scaffolds that may bind PDE10A via structure-based drug design. For this task, we use INPHARMA, an NMR-based method that measures protein-mediated interligand NOEs between pairs of weakly, competitively binding ligands. INPHARMA is applied to a combination of four chemically diverse PDE10A binding fragments, with the aim of merging their pharmacophoric features into a larger, tighter binding molecule. All four ligands bind the PDE10A cAMP binding domain with affinity in the micromolar range. The application of INPHARMA to identify the correct docking poses of these ligands is challenging due to the nature of the binding pocket and the high content of water-mediated intermolecular contacts. Nevertheless, ensemble docking in the presence of conserved water molecules generates docking poses that are in agreement with all sets of INPHARMA data. These poses are used to build a pharmacophore model with which we search the ZINC database.

Allostery Wiring Map for Kinesin Energy Transduction and Its Evolution

How signals between the kinesin active and cytoskeletal binding sites are transmitted is an open question and an allosteric question. By extracting correlated evolutionary changes within 700+ sequences, we built a model of residues that are energetically coupled and that define molecular routes for signal transmission. Typically, these coupled residues are located at multiple distal sites and thus are predicted to form a complex, non-linear network that wires together different functional sites in the protein. Of note, our model connected the site for ATP hydrolysis with sites that ultimately utilize its free energy, such as the microtubule-binding site, drug-binding loop 5, and necklinker. To confirm the calculated energetic connectivity between non-adjacent residues, double-mutant cycle analysis was conducted with 22 kinesin mutants. There was a direct correlation between thermodynamic coupling in experiment and evolutionarily derived energetic coupling. We conclude that energy transduction is coordinated by multiple distal sites in the protein rather than only being relayed through adjacent residues. Moreover, this allosteric map forecasts how energetic orchestration gives rise to different nanomotor behaviors within the superfamily.

DNA Polymerase Conformational Dynamics and the Role of Fidelity-Conferring Residues: Insights from Computational Simulations

Herein we investigate the molecular bases of DNA polymerase I conformational dynamics that underlie the replication fidelity of the enzyme. Such fidelity is determined by conformational changes that promote the rejection of incorrect nucleotides before the chemical ligation step. We report a comprehensive atomic resolution study of wild type and mutant enzymes in different bound states and starting from different crystal structures, using extensive molecular dynamics (MD) simulations that cover a total timespan of ~5 ms. The resulting trajectories are examined via a combination of novel methods of internal dynamics and energetics analysis, aimed to reveal the principal molecular determinants for the (de)stabilization of a certain conformational state. Our results show that the presence of fidelity-decreasing mutations or the binding of incorrect nucleotides in ternary complexes tend to favor transitions from closed toward open structures, passing through an ensemble of semi-closed intermediates. The latter ensemble includes the experimentally observed ajar conformation which, consistent with previous experimental observations, emerges as a molecular checkpoint for the selection of the correct nucleotide to incorporate. We discuss the implications of our results for the understanding of the relationships between the structure, dynamics, and function of DNA polymerase I at the atomistic level.

Conformational Sampling and Binding Site Assessment of Suppression of Tumorigenicity 2 Ectodomain

Suppression of Tumorigenicity 2 (ST2), a member of the interleukin-1 receptor (IL-1R) family, activates type 2 immune responses to pathogens and tissue damage via binding to IL-33. Dysregulated responses contribute to asthma, graft-versus-host and autoinflammatory diseases and disorders. To study ST2 structure for inhibitor development, we performed the principal component (PC) analysis on the crystal structures of IL1-1R1, IL1-1R2, ST2 and the refined ST2 ectodomain (ST2ECD) models, constructed from previously reported small-angle X-ray scattering data. The analysis facilitates mapping of the ST2ECD conformations to PC subspace for characterizing structural changes. Extensive coverage of ST2ECD conformations was then obtained using the accelerated molecular dynamics simulations started with the IL-33 bound ST2ECD structure as instructed by their projected locations on the PC subspace. Cluster analysis of all conformations further determined representative conformations of ST2ECD ensemble in solution. Alignment of the representative conformations with the ST2/IL-33 structure showed that the D3 domain of ST2ECD (containing D1-D3 domains) in most conformations exhibits no clashes with IL-33 in the crystal structure. Our experimental binding data informed that the D1-D2 domain of ST2ECD contributes predominantly to the interaction between ST2ECD and IL-33 underscoring the importance of the D1-D2 domain in binding. Computational binding site assessment revealed one third of the total detected binding sites in the representative conformations may be suitable for binding to potent small molecules. Locations of these sites include the D1-D2 domain ST2ECD and modulation sites conformed to ST2ECD conformations. Our study provides structural models and analyses of ST2ECD that could be useful for inhibitor discovery.

Dynamic Coupling and Allosteric Networks in the α Subunit of Heterotrimeric G Proteins

G protein α subunits cycle between active and inactive conformations to regulate a multitude of intracellular signaling cascades. Important structural transitions occurring during this cycle have been characterized from extensive crystallographic studies. However, the link between observed conformations and the allosteric regulation of binding events at distal sites critical for signaling through G proteins remain unclear. Here we describe molecular dynamics simulations, bioinformatics analysis, and experimental mutagenesis that identifies residues involved in mediating the allosteric coupling of receptor, nucleotide, and helical domain interfaces of Gαi. Most notably, we predict and characterize novel allosteric decoupling mutants, which display enhanced helical domain opening, increased rates of nucleotide exchange, and constitutive activity in the absence of receptor activation. Collectively, our results provide a framework for explaining how binding events and mutations can alter internal dynamic couplings critical for G protein function.

Enrichment of druggable conformations from apo protein structures using cosolvent-accelerated molecular dynamics

Here we describe the development of an improved workflow for utilizing experimental and simulated protein conformations in the structure-based design of inhibitors for anti-apoptotic Bcl-2 family proteins. Traditional structure-based approaches on similar targets are often constrained by the sparsity of available structures and difficulties in finding lead compounds that dock against flat, flexible protein-protein interaction surfaces. By employing computational docking of known small molecule inhibitors, we have demonstrated that structural ensembles derived from either accelerated MD (aMD) or MD in the presence of an organic cosolvent generally give better scores than those assessed from analogous conventional MD. Furthermore, conformations obtained from combined cosolvent aMD simulations started with the apo-Bcl-xL structure yielded better average and minimum docking scores for known binders than an ensemble of 72 experimental apo- and ligand-bound Bcl-xL structures. A detailed analysis of the simulated conformations indicates that the aMD effectively enhanced conformational sampling of the flexible helices flanking the main Bcl-xL binding groove, permitting the cosolvent acting as small ligands to penetrate more deeply into the binding pocket and shape ligand-bound conformations not evident in conventional simulations. We believe this approach could be useful for identifying inhibitors against other protein-protein interaction systems involving highly flexible binding sites, particularly for targets with less accumulated structural data.

Nucleotide-dependent structural fluctuations and regulation of microtubule-binding affinity of KIF1A

Molecular motors such as kinesin regulate affinity to a rail protein during the ATP hydrolysis cycle. The regulation mechanism, however, is yet to be determined. To understand this mechanism, we investigated the structural fluctuations of the motor head of the single-headed kinesin called KIF1A in different nucleotide states using molecular dynamics simulations of a Gō-like model. We found that the helix α4 at the microtubule (MT) binding site intermittently exhibits a large structural fluctuation when MT is absent. Frequency of this fluctuation changes systematically according to the nucleotide states and correlates strongly with the experimentally observed binding affinity to MT. We also showed that thermal fluctuation enhances the correlation and the interaction with the nucleotide suppresses the fluctuation of the helix α4. These results suggest that KIF1A regulates affinity to MT by changing the flexibility of the helix α4 during the ATP hydrolysis process: the binding site becomes more flexible in the strong binding state than in the weak binding state.

Decrypting the structural, dynamic, and energetic basis of a monomeric kinesin interacting with a tubulin dimer in three ATPase states by all-atom molecular dynamics simulation

We have employed molecular dynamics (MD) simulation to investigate, with atomic details, the structural dynamics and energetics of three major ATPase states (ADP, APO, and ATP state) of a human kinesin-1 monomer in complex with a tubulin dimer. Starting from a recently solved crystal structure of ATP-like kinesin-tubulin complex by the Knossow lab, we have used flexible fitting of cryo-electron-microscopy maps to construct new structural models of the kinesin-tubulin complex in APO and ATP state, and then conducted extensive MD simulations (total 400 ns for each state), followed by flexibility analysis, principal component analysis, hydrogen bond analysis, and binding free energy analysis. Our modeling and simulation have revealed key nucleotide-dependent changes in the structure and flexibility of the nucleotide-binding pocket (featuring a highly flexible and open switch I in APO state) and the tubulin-binding site, and allosterically coupled motions driving the APO to ATP transition. In addition, our binding free energy analysis has identified a set of key residues involved in kinesin-tubulin binding. On the basis of our simulation, we have attempted to address several outstanding issues in kinesin study, including the possible roles of β-sheet twist and neck linker docking in regulating nucleotide release and binding, the structural mechanism of ADP release, and possible extension and shortening of α4 helix during the ATPase cycle. This study has provided a comprehensive structural and dynamic picture of kinesin's major ATPase states, and offered promising targets for future mutational and functional studies to investigate the molecular mechanism of kinesin motors.

Integrating protein structural dynamics and evolutionary analysis with Bio3D

Background

Popular bioinformatics approaches for studying protein functional dynamics include comparisons of crystallographic structures, molecular dynamics simulations and normal mode analysis. However, determining how observed displacements and predicted motions from these traditionally separate analyses relate to each other, as well as to the evolution of sequence, structure and function within large protein families, remains a considerable challenge. This is in part due to the general lack of tools that integrate information of molecular structure, dynamics and evolution.

Results

Here, we describe the integration of new methodologies for evolutionary sequence, structure and simulation analysis into the Bio3D package. This major update includes unique high-throughput normal mode analysis for examining and contrasting the dynamics of related proteins with non-identical sequences and structures, as well as new methods for quantifying dynamical couplings and their residue-wise dissection from correlation network analysis. These new methodologies are integrated with major biomolecular databases as well as established methods for evolutionary sequence and comparative structural analysis. New functionality for directly comparing results derived from normal modes, molecular dynamics and principal component analysis of heterogeneous experimental structure distributions is also included. We demonstrate these integrated capabilities with example applications to dihydrofolate reductase and heterotrimeric G-protein families along with a discussion of the mechanistic insight provided in each case.

Conclusions

The integration of structural dynamics and evolutionary analysis in Bio3D enables researchers to go beyond a prediction of single protein dynamics to investigate dynamical features across large protein families. The Bio3D package is distributed with full source code and extensive documentation as a platform independent R package under a GPL2 license from http://thegrantlab.org/bio3d/.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-014-0399-6) contains supplementary material, which is available to authorized users.

Conserved mechanisms of microtubule-stimulated ADP release, ATP binding, and force generation in transport kinesins

Kinesins are a superfamily of microtubule-based ATP-powered motors, important for multiple, essential cellular functions. How microtubule binding stimulates their ATPase and controls force generation is not understood. To address this fundamental question, we visualized microtubule-bound kinesin-1 and kinesin-3 motor domains at multiple steps in their ATPase cycles—including their nucleotide-free states—at ∼7 Å resolution using cryo-electron microscopy. In both motors, microtubule binding promotes ordered conformations of conserved loops that stimulate ADP release, enhance microtubule affinity and prime the catalytic site for ATP binding. ATP binding causes only small shifts of these nucleotide-coordinating loops but induces large conformational changes elsewhere that allow force generation and neck linker docking towards the microtubule plus end. Family-specific differences across the kinesin–microtubule interface account for the distinctive properties of each motor. Our data thus provide evidence for a conserved ATP-driven mechanism for kinesins and reveal the critical mechanistic contribution of the microtubule interface.

DOI:http://dx.doi.org/10.7554/eLife.03680.001

The interior of a cell is a hive of activity, filled with proteins and other items moving from one location to another. A network of filaments called microtubules forms tracks along which so-called motor proteins carry these items. Kinesins are one group of motor proteins, and a typical kinesin protein has one end (called the ‘motor domain’) that can attach itself to the microtubules. The other end links to the cargo being carried, and a ‘neck’ connects the two. When two of these proteins work together, flexible regions of the neck allow the two motor domains to move past one another, which enable the kinesin to essentially walk along a microtubule in a stepwise manner.

To take these steps along microtubules, each kinesin motor domain in the pair must undergo alternating cycles of tight association and release from their tracks. This cycle is coordinated by binding and breaking down a molecule called ATP, which also provides the energy needed to take the next step. How the cycle of loose and tight microtubule attachment is coordinated with the release of the breakdown products of ATP, and how the energy from the ATP molecule is converted into the force that moves the motor along the microtubule, has been unclear.

Atherton et al. use a technique called cryo-electron microscopy to study—in more detail than previously seen—the structure of the motor domains of two types of kinesin called kinesin-1 and kinesin-3. Images were taken at different stages of the cycle used by the motor domains to extract the energy from ATP molecules. Although the two kinesins have been thought to move along the microtubule tracks in different ways, Atherton et al. find that the core mechanism used by their motor domains is the same.

When a motor domain binds to the microtubule, its shape changes, first stimulating release of the breakdown products of ATP from the previous cycle. This release makes room for a new ATP molecule to bind. The structural changes caused by ATP binding are relatively small but produce larger changes in the flexible neck region that enable individual motor domains within a kinesin pair to co-ordinate their movement and move in a consistent direction. This mechanism involves tight coupling between track binding and fuel usage and makes kinesins highly efficient motors.

The structures uncovered by Atherton et al. reveal a mechanism that links microtubule binding, the energy supplied to the motor domain and the force that moves the kinesin along a microtubule. Future work will clarify whether the key features observed in the motor domains of kinesin-1 and kinesin-3 are also found in other types of kinesin motors.

DOI:http://dx.doi.org/10.7554/eLife.03680.002

Mapping the structural and dynamical features of kinesin motor domains

Kinesin motor proteins drive intracellular transport by coupling ATP hydrolysis to conformational changes that mediate directed movement along microtubules. Characterizing these distinct conformations and their interconversion mechanism is essential to determining an atomic-level model of kinesin action. Here we report a comprehensive principal component analysis of 114 experimental structures along with the results of conventional and accelerated molecular dynamics simulations that together map the structural dynamics of the kinesin motor domain. All experimental structures were found to reside in one of three distinct conformational clusters (ATP-like, ADP-like and Eg5 inhibitor-bound). These groups differ in the orientation of key functional elements, most notably the microtubule binding α4–α5, loop8 subdomain and α2b-β4-β6-β7 motor domain tip. Group membership was found not to correlate with the nature of the bound nucleotide in a given structure. However, groupings were coincident with distinct neck-linker orientations. Accelerated molecular dynamics simulations of ATP, ADP and nucleotide free Eg5 indicate that all three nucleotide states could sample the major crystallographically observed conformations. Differences in the dynamic coupling of distal sites were also evident. In multiple ATP bound simulations, the neck-linker, loop8 and the α4–α5 subdomain display correlated motions that are absent in ADP bound simulations. Further dissection of these couplings provides evidence for a network of dynamic communication between the active site, microtubule-binding interface and neck-linker via loop7 and loop13. Additional simulations indicate that the mutations G325A and G326A in loop13 reduce the flexibility of these regions and disrupt their couplings. Our combined results indicate that the reported ATP and ADP-like conformations of kinesin are intrinsically accessible regardless of nucleotide state and support a model where neck-linker docking leads to a tighter coupling of the microtubule and nucleotide binding regions. Furthermore, simulations highlight sites critical for large-scale conformational changes and the allosteric coupling between distal functional sites.

Kinesin motor proteins transport cargo along microtubule tracks to support essential cellular functions including cell growth, axonal signaling and the separation of chromosomes during cell division. All kinesins contain one or more conserved motor domains that modulate binding and movement along microtubules via cycles of ATP hydrolysis. Important conformational transitions occurring during this cycle have been characterized with extensive crystallographic studies. However, the link between the observed conformations and the mechanisms involved in conformational change and microtubule interaction modulation remain unclear. Here we describe a comprehensive principal component analysis of available motor domain crystallographic structures supplemented with extensive unbiased conventional and accelerated molecular dynamics simulations that together characterize the response of kinesin motor domains to ATP binding and hydrolysis. Our studies reveal atomic details of conformational transitions, as well as novel nucleotide-dependent dynamical couplings, of distal regions and residues potentially important for the allosteric link between nucleotide and microtubule binding sites.

Structure of a kinesin-tubulin complex and implications for kinesin motility

The typical function of kinesins is to transport cargo along microtubules. Binding of ATP to microtubule-attached motile kinesins leads to cargo displacement. To better understand the nature of the conformational changes that lead to the power stroke that moves a kinesin's load along a microtubule, we determined the X-ray structure of human kinesin-1 bound to αβ-tubulin. The structure defines the mechanism of microtubule-stimulated ATP hydrolysis, which releases the kinesin motor domain from microtubules. It also reveals the structural linkages that connect the ATP nucleotide to the kinesin neck linker, a 15-amino acid segment C terminal to the catalytic core of the motor domain, to result in the power stroke. ATP binding to the microtubule-bound kinesin favors neck-linker docking. This biases the attachment of kinesin's second head in the direction of the movement, thus initiating each of the steps taken.

The amino acid substitutions rtP177G and rtF249A in the reverse transcriptase domain of hepatitis B virus polymerase reduce the susceptibility to tenofovir

Long term antiviral therapy with nucleoside/nucleotide analogs have been routinely used to treat chronic hepatitis B virus (HBV) infection but may lead to the emergence of drug-resistant viral mutants. However, the HBV resistance mutations for tenofovir (TDF) remain controversial. It is speculated that the genetic barrier for TDF resistance may be high for HBV. We asked whether selected amino acid substitutions in HBV polymerase may reduce susceptibility to TDF. A series of amino acids in HBV polymerase were selected based on bioinformatics analysis for mutagenesis. The replication competence and susceptibility to TDF of the mutated HBV clones were determined both in vitro and in vivo. nineteen mutations in HBV polymerase were included and impaired the replication competence of HBV genome in different degrees. The mutations at rtL77F (sS69C), rtF88L (sF80Y), and rtP177G (sR169G) also significantly affected HBsAg expression. The HBV mutants with rtP177G and rtF249A were found to have reduced susceptibility to TDF in vitro with a resistance index of 2.53 and 12.16, respectively. The testing in in vivo model based on the hydrodynamic injection revealed the antiviral effect of TDF against wild type and mutated HBV genomes and confirmed the reduced the susceptibility of mutant HBV to TDF.

Principal component and clustering analysis on molecular dynamics data of the ribosomal L11·23S subdomain

With improvements in computer speed and algorithm efficiency, MD simulations are sampling larger amounts of molecular and biomolecular conformations. Being able to qualitatively and quantitatively sift these conformations into meaningful groups is a difficult and important task, especially when considering the structure-activity paradigm. Here we present a study that combines two popular techniques, principal component (PC) analysis and clustering, for revealing major conformational changes that occur in molecular dynamics (MD) simulations. Specifically, we explored how clustering different PC subspaces effects the resulting clusters versus clustering the complete trajectory data. As a case example, we used the trajectory data from an explicitly solvated simulation of a bacteria’s L11·23S ribosomal subdomain, which is a target of thiopeptide antibiotics. Clustering was performed, using K-means and average-linkage algorithms, on data involving the first two to the first five PC subspace dimensions. For the average-linkage algorithm we found that data-point membership, cluster shape, and cluster size depended on the selected PC subspace data. In contrast, K-means provided very consistent results regardless of the selected subspace. Since we present results on a single model system, generalization concerning the clustering of different PC subspaces of other molecular systems is currently premature. However, our hope is that this study illustrates a) the complexities in selecting the appropriate clustering algorithm, b) the complexities in interpreting and validating their results, and c) by combining PC analysis with subsequent clustering valuable dynamic and conformational information can be obtained.

The landscape of the prion protein's structural response to mutation revealed by principal component analysis of multiple NMR ensembles

Prion Proteins (PrP) are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive principle components analysis (PCA) on the tertiary structures of PrP globular proteins to discern PrP subdomains that exhibit conformational change in response to point mutations and clade-specific evolutionary sequence mutation trends. This is to our knowledge the first such large-scale analysis of multiple NMR ensembles of protein structures, and the first study of its kind for PrPs. We conducted PCA on human (n = 11), mouse (n = 14), and wildtype (n = 21) sets of PrP globular structures, from which we identified five conformationally variable subdomains within PrP. PCA shows that different non-local patterns and rankings of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we have observed the conformational clustering of divergent TSE-non-susceptible species pairs; these non-phylogenetic clusterings indicate structural solutions towards TSE resistance that do not necessarily coincide with evolutionary divergence. We discuss the novelty of our approach and the importance of PrP subdomains in structural conversion during disease.

Prion Proteins (PrP) cause a variety of incurable TSE diseases, and are among a small number of proteins for which large numbers of NMR ensembles have been resolved for sequence mutants and diverse species. Here, we perform a comprehensive PCA study to assess conformational variation and discern the landscape of the PrP structural response to sequence mutation. This is to our knowledge the first large-scale analysis of multiple NMR ensembles for a specific protein, and the first study to perform a multivariate PCA on the native globular structures of PrP. We conducted exhaustive PCA on three PrP subsets: human and mouse subsets that include structures of sequence mutants, and the set of wild-type PrP (16 PrP species). PCA shows that different non-local patterns of variable subdomains arise for different pathogenic mutants. These subdomains may thus be key areas for initiating PrP conversion during disease. Furthermore, we observed that some evolutionarily divergent species that are non-susceptible to TSEs have surprising structural similarities in their PrPs. We discuss the novelty of our approach with respect to prions, and the advantage of this analysis as a fast, reliable starting point to identify interesting domains that may warrant further experimental and computational analysis.

Identification of potential small molecule binding pockets on Rho family GTPases

Rho GTPases are conformational switches that control a wide variety of signaling pathways critical for eukaryotic cell development and proliferation. They represent attractive targets for drug design as their aberrant function and deregulated activity is associated with many human diseases including cancer. Extensive high-resolution structures (>100) and recent mutagenesis studies have laid the foundation for the design of new structure-based chemotherapeutic strategies. Although the inhibition of Rho signaling with drug-like compounds is an active area of current research, very little attention has been devoted to directly inhibiting Rho by targeting potential allosteric non-nucleotide binding sites. By avoiding the nucleotide binding site, compounds may minimize the potential for undesirable off-target interactions with other ubiquitous GTP and ATP binding proteins. Here we describe the application of molecular dynamics simulations, principal component analysis, sequence conservation analysis, and ensemble small-molecule fragment mapping to provide an extensive mapping of potential small-molecule binding pockets on Rho family members. Characterized sites include novel pockets in the vicinity of the conformationaly responsive switch regions as well as distal sites that appear to be related to the conformations of the nucleotide binding region. Furthermore the use of accelerated molecular dynamics simulation, an advanced sampling method that extends the accessible time-scale of conventional simulations, is found to enhance the characterization of novel binding sites when conformational changes are important for the protein mechanism.

All-atom structural investigation of kinesin-microtubule complex constrained by high-quality cryo-electron-microscopy maps

In this study, we have performed a comprehensive structural investigation of three major biochemical states of a kinesin complexed with microtubule under the constraint of high-quality cryo-electron-microscopy (EM) maps. In addition to the ADP and ATP state which were captured by X-ray crystallography, we have also modeled the nucleotide-free or APO state for which no crystal structure is available. We have combined flexible fitting of EM maps with regular molecular dynamics simulations, hydrogen-bond analysis, and free energy calculation. Our APO-state models feature a subdomain rotation involving loop L2 and α6 helix of kinesin, and local structural changes in active site similar to a related motor protein, myosin. We have identified a list of hydrogen bonds involving key residues in the active site and the binding interface between kinesin and microtubule. Some of these hydrogen bonds may play an important role in coupling microtubule binding to ATPase activities in kinesin. We have validated our models by calculating the binding free energy between kinesin and microtubule, which quantitatively accounts for the observation of strong binding in the APO and ATP state and weak binding in the ADP state. This study will offer promising targets for future mutational and functional studies to investigate the mechanism of kinesin motors.

Dynamics, flexibility and ligand-induced conformational changes in biological macromolecules: a computational approach

Biomolecules possess important dynamical properties that enable them to adapt and alternate their conformation as a response to environmental stimuli. Recent advancements in computational resources and methodology allow a higher capability to mimic in vitro conditions and open up the possibility of studying large systems over longer timescales. Here, we describe commonly used computational approaches for studying the dynamic properties of proteins. We review a selected set of simulation studies on ligand-induced changes in the chaperonin GroEL-GroES, a molecular folding machine, maltose-binding protein, a prototypical member of the periplasmic binding proteins, and the bacterial ribosomal A-site, focusing on aminoglycoside antibiotic recognition. We also discuss a recent quantitative reconstruction of the binding process of benzamidine and trypsin. These studies contribute to the understanding and further development of the medicinal regulation of large biomolecular systems.

Electrostatically biased binding of kinesin to microtubules

An electrostatic field rotates, slides, and guides the kinesin head to bind the microtubule at a site a short distance ahead, thus determining the direction of movement of the motor.

The minimum motor domain of kinesin-1 is a single head. Recent evidence suggests that such minimal motor domains generate force by a biased binding mechanism, in which they preferentially select binding sites on the microtubule that lie ahead in the progress direction of the motor. A specific molecular mechanism for biased binding has, however, so far been lacking. Here we use atomistic Brownian dynamics simulations combined with experimental mutagenesis to show that incoming kinesin heads undergo electrostatically guided diffusion-to-capture by microtubules, and that this produces directionally biased binding. Kinesin-1 heads are initially rotated by the electrostatic field so that their tubulin-binding sites face inwards, and then steered towards a plus-endwards binding site. In tethered kinesin dimers, this bias is amplified. A 3-residue sequence (RAK) in kinesin helix alpha-6 is predicted to be important for electrostatic guidance. Real-world mutagenesis of this sequence powerfully influences kinesin-driven microtubule sliding, with one mutant producing a 5-fold acceleration over wild type. We conclude that electrostatic interactions play an important role in the kinesin stepping mechanism, by biasing the diffusional association of kinesin with microtubules.

Animal and plant cells contain a molecular-scale “railway” network, in which the tracks, called microtubules, radiate out from the cell centre and locomotive proteins, called kinesins, haul their molecular cargoes along the microtubule tracks. This railway system transports many different cargoes to where they are needed, so it is crucial for the cell's organization and function. Breakdowns in this transport system can cause diseases like Alzheimer's, and drugs that temporarily halt transport make powerful anti-cancer agents. Precisely how kinesin motor proteins move along their microtubule tracks is an important question in biology. We know that some kinesins have twin “heads” that alternately bind to and step along microtubules in a coordinated walking action. But more usually, kinesins have only one head. How single-headed kinesins produce force and movement is poorly understood. In this study, we address this question and show that electrical attraction between single kinesin heads and microtubules is a critical factor deciding the direction of movement: each time the head approaches a microtubule, it slides forwards by the electrical attraction between the engine and the track.

Novel allosteric sites on Ras for lead generation

Aberrant Ras activity is a hallmark of diverse cancers and developmental diseases. Unfortunately, conventional efforts to develop effective small molecule Ras inhibitors have met with limited success. We have developed a novel multi-level computational approach to discover potential inhibitors of previously uncharacterized allosteric sites. Our approach couples bioinformatics analysis, advanced molecular simulations, ensemble docking and initial experimental testing of potential inhibitors. Molecular dynamics simulation highlighted conserved allosteric coupling of the nucleotide-binding switch region with distal regions, including loop 7 and helix 5. Bioinformatics methods identified novel transient small molecule binding pockets close to these regions and in the vicinity of the conformationally responsive switch region. Candidate binders for these pockets were selected through ensemble docking of ZINC and NCI compound libraries. Finally, cell-based assays confirmed our hypothesis that the chosen binders can inhibit the downstream signaling activity of Ras. We thus propose that the predicted allosteric sites are viable targets for the development and optimization of new drugs.

Structural comparison of F1-ATPase: interplay among enzyme structures, catalysis, and rotations

F(1)-ATPase, a rotary motor powered by adenosine triphosphate hydrolysis, has been extensively studied by various methods. Here, we performed a systematic comparison of 29 X-ray crystal structures of F(1)-complexes, finding fine interplay among enzyme structures, catalysis, and rotations. First, analyzing the 87 structures of enzymatic αβ-subunits, we confirmed that the two modes, the hinge motion of β-subunit and the loose/tight motion of the αβ-interface, dominate the variations. The structural ensemble was nearly contiguous bridging three clusters, αβ(TP), αβ(DP), and αβ(E). Second, the catalytic site analysis suggested the correlation between the phosphate binding and the tightening of the αβ-interface. Third, addressing correlations of enzymatic structures with the orientations of the central stalk γ, we found that the γ rotation highly correlates with loosening of αβ(E)-interface and β(DP) hinge motions. Finally, calculating the helix 6 angle of β, we identified the recently observed partially closed conformation being consistent with β(HC).

Conformational selection in G-proteins: lessons from Ras and Rho

The induced fit model has traditionally been invoked to describe the activating conformational change of the monomeric G-proteins, such as Ras and Rho. With this scheme, the presence or absence of the γ-phosphate of GTP leads to an instantaneous switch in conformation. Here we describe atomistic molecular simulations that demonstrate that both Ras and Rho superfamily members harbor an intrinsic susceptibility to sample multiple conformational states in the absence of nucleotide ligand. By comparing the distribution of conformers in the presence and absence of nucleotide, we show that conformational selection is the dominant mechanism by which Ras and Rho undergo nucleotide-dependent conformational changes. Furthermore, the pattern of correlated motions revealed by these simulations predicts a preserved allosteric coupling of the nucleotide-binding site with the membrane interacting C-terminus in both Rho and Ras.

Three clusters of conformational states in p450cam reveal a multistep pathway for closing of the substrate access channel

Conformational changes in the substrate access channel have been observed for several forms of cytochrome P450, but the extent of conformational plasticity exhibited by a given isozyme has not been completely characterized. Here we present crystal structures of P450cam bound to a library of 12 active site probes containing a substrate analogue tethered to a variable linker. The structures provide a unique view of the range of protein conformations accessible during substrate binding. Principal component analysis of a total of 30 structures reveals three discrete clusters of conformations: closed (P450cam-C), intermediate (P450cam-I), and fully open (P450cam-O). Relative to P450cam-C, the P450cam-I state results predominantly from a retraction of helix F, while both helices F and G move in concert to reach the fully open P450cam-O state. Both P450cam-C and P450cam-I are well-defined states, while P450cam-O shows evidence of a somewhat broader distribution of conformations and includes the open form recently seen in the absence of substrate. The observed clustering of protein conformations over a wide range of ligand variants suggests a multistep closure of the enzyme around the substrate that begins by conformational selection from an ensemble of open conformations and proceeds through a well-defined intermediate, P450cam-I, before full closure to the P450cam-C state in the presence of small substrates. This multistep pathway may have significant implications for a full understanding of substrate specificity, kinetics, and coupling of substrate binding to P450 function.

Mechanism of processive movement of monomeric and dimeric kinesin molecules

Kinesin molecules are motor proteins capable of moving along microtubule by hydrolyzing ATP. They generally have several forms of construct. This review focuses on two of the most studied forms: monomers such as KIF1A (kinesin-3 family) and dimers such as conventional kinesin (kinesin-1 family), both of which can move processively towards the microtubule plus end. There now exist numerous models that try to explain how the kinesin molecules convert the chemical energy of ATP hydrolysis into the mechanical energy to "power" their processive movement along microtubule. Here, we attempt to present a comprehensive review of these models. We further propose a new hybrid model for the dimeric kinesin by combining the existing models and provide a framework for future studies in this subject.

Insight into the molecular mechanism of the multitasking kinesin-8 motor

Members of the kinesin-8 motor class have the remarkable ability to both walk towards microtubule plus-ends and depolymerise these ends on arrival, thereby regulating microtubule length. To analyse how kinesin-8 multitasks, we studied the structure and function of the kinesin-8 motor domain. We determined the first crystal structure of a kinesin-8 and used cryo-electron microscopy to calculate the structure of the microtubule-bound motor. Microtubule-bound kinesin-8 reveals a new conformation compared with the crystal structure, including a bent conformation of the α4 relay helix and ordering of functionally important loops. The kinesin-8 motor domain does not depolymerise stabilised microtubules with ATP but does form tubulin rings in the presence of a non-hydrolysable ATP analogue. This shows that, by collaborating, kinesin-8 motor domain molecules can release tubulin from microtubules, and that they have a similar mechanical effect on microtubule ends as kinesin-13, which enables depolymerisation. Our data reveal aspects of the molecular mechanism of kinesin-8 motors that contribute to their unique dual motile and depolymerising functions, which are adapted to control microtubule length.

Kinesin-14: the roots of reversal

Kinesin-14 motor proteins step towards microtubule minus ends, in the opposite direction to other kinesins. Work on the still-enigmatic kinesin-14 mechanism published in BMC Structural Biology shows that the carboxyl terminus of the motor head undergoes a dock-undock cycle, like that of plus-end-directed kinesins.

See research article: http://www.biomedcentral.com/1472-6807/10/19

Structure, evolutionary conservation, and conformational dynamics of Homo sapiens fascin-1, an F-actin crosslinking protein

Eukaryotes have several highly conserved actin-binding proteins that crosslink filamentous actin into compact ordered bundles present in distinct cytoskeletal processes, including microvilli, stereocilia and filopodia. Fascin is an actin-binding protein that is present predominantly in filopodia, which are believed to play a central role in normal and aberrant cell migration. An important outstanding question regards the molecular basis for the unique localization and functional properties of fascin compared with other actin crosslinking proteins. Here, we present the crystal structure of full-length Homo sapiens fascin-1, and examine its packing, conformational flexibility, and evolutionary sequence conservation. The structure reveals a novel arrangement of four tandem beta-trefoil domains that form a bi-lobed structure with approximate pseudo 2-fold symmetry. Each lobe has internal approximate pseudo 2-fold and pseudo 3-fold symmetry axes that are approximately perpendicular, with beta-hairpin triplets located symmetrically on opposite sides of each lobe that mutational data suggest are actin-binding domains. Sequence conservation analysis confirms the importance of hydrophobic core residues that stabilize the beta-trefoil fold, as well as interfacial residues that are likely to stabilize the overall fascin molecule. Sequence conservation also indicates highly conserved surface patches near the putative actin-binding domains of fascin, which conformational dynamics analysis suggests to be coupled via an allosteric mechanism that might have important functional implications for F-actin crosslinking by fascin.

Ras conformational switching: simulating nucleotide-dependent conformational transitions with accelerated molecular dynamics

Ras mediates signaling pathways controlling cell proliferation and development by cycling between GTP- and GDP-bound active and inactive conformational states. Understanding the complete reaction path of this conformational change and its intermediary structures is critical to understanding Ras signaling. We characterize nucleotide-dependent conformational transition using multiple-barrier-crossing accelerated molecular dynamics (aMD) simulations. These transitions, achieved for the first time for wild-type Ras, are impossible to observe with classical molecular dynamics (cMD) simulations due to the large energetic barrier between end states. Mapping the reaction path onto a conformer plot describing the distribution of the crystallographic structures enabled identification of highly populated intermediate structures. These structures have unique switch orientations (residues 25–40 and 57–75) intermediate between GTP and GDP states, or distinct loop3 (46–49), loop7 (105–110), and α5 C-terminus (159–166) conformations distal from the nucleotide-binding site. In addition, these barrier-crossing trajectories predict novel nucleotide-dependent correlated motions, including correlations of α2 (residues 66–74) with α3-loop7 (93–110), loop2 (26–37) with loop10 (145–151), and loop3 (46–49) with α5 (152–167). The interconversion between newly identified Ras conformations revealed by this study advances our mechanistic understanding of Ras function. In addition, the pattern of correlated motions provides new evidence for a dynamic linkage between the nucleotide-binding site and the membrane interacting C-terminus critical for the signaling function of Ras. Furthermore, normal mode analysis indicates that the dominant collective motion that occurs during nucleotide-dependent conformational exchange, and captured in aMD (but absent in cMD) simulations, is a low-frequency motion intrinsic to the structure.

The Ras family of enzymes mediate signaling pathways controlling cell proliferation and development by cycling between active and inactive conformational states. Mutations that affect the ability to switch between states are associated with a variety of cancers. However, details of how the structural changes occur and how mutations affect the fidelity of this process remain to be determined. Here we employ an advanced computational technique, termed accelerated molecular dynamics, to characterize structural transitions and identify novel highly populated transient conformations. Several spatially distant structural regions were found to undergo correlated motions, highlighting a dynamic linkage between the sites of enzymatic reaction and the membrane-interacting C-terminus. In addition, our results indicate that the major motion occurring during the conformational exchange is a low-frequency motion intrinsic to the structure. Hence, features of the characterized transitions likely apply to a large number of structurally similar but functionally diverse nucleotide triphosphatases. These results provide fresh insights into how oncogenic mutations might modulate conformational transitions in Ras.

9-Angström structure of a microtubule-bound mitotic motor

Kinesin-5 (K5) motors are important components of the microtubule (MT)-based cell division machinery and are targets for small-molecule inhibitors currently in cancer clinical trials. However, the nature of the K5-MT interaction and the regulatory mechanisms that control it remain unclear. Using cryo-electron microscopy and image processing, we calculated the structure of a K5 motor bound to MTs at 9 A resolution, providing insight into this important interaction. Our reconstruction reveals the K5 motor domain in an ATP-like conformation in which MT binding induces the conserved nucleotide-sensing switch I and II loops to form a compact subdomain around the bound nucleotide. Our reconstruction also reveals a novel conformation for the K5-specific drug-binding loop 5, suggesting a possible role for it in switching K5s between force generation and diffusional modes of MT binding. Our data thus shed light on regulation of the interaction between spindle components important for chromosome segregation.

Force generation in kinesin hinges on cover-neck bundle formation

In kinesin motors, a fundamental question concerns the mechanism by which ATP binding generates the force required for walking. Analysis of available structures combined with molecular dynamics simulations demonstrates that the conformational change of the neck linker involves the nine-residue-long N-terminal region, the cover strand, as an element that is essential for force generation. Upon ATP binding, it forms a beta sheet with the neck linker, the cover-neck bundle, which induces the forward motion of the neck linker, followed by a latch-type binding to the motor head. The estimated stall force and anisotropic response to external loads calculated from the model agree with force-clamp measurements. The proposed mechanism for force generation by the cover-neck bundle formation appears to apply to several kinesin families. It also elucidates the design principle of kinesin as the smallest known processive motor.

The role of microtubules in processive kinesin movement

Kinesins are microtubule-based motors that are important for various intracellular transport processes. To understand the mechanism of kinesin movement, X-ray crystallography has been used to study the atomic structures of kinesin. However, as crystal structures of kinesin alone accumulate, it is becoming clear that kinesin structures should also be investigated with the microtubule to understand the contribution of the microtubule track to the nucleotide-induced conformational changes of kinesin. Recently, several high-resolution structures of kinesin with microtubules were obtained using cryo-electron microscopy. Comparison with X-ray crystallographic structures revealed the importance of the microtubule in determining the conformation of kinesin. Together with recent biophysical data, we describe different structural models of processive kinesin movement and provide a framework for future experiments.

Molecular motors: not quite like clockwork

Models commonly used to explain the mechanism of myosin motors typically include a power stroke that is attributed to a conformational change in the motor domain and amplified by a long lever arm that connects the motor domain to the cargo. Similar models have proved less enlightening in the case of microtubule motors, for which it may be more helpful to consider models involving thermally driven mechanisms.

Processivity of single-headed kinesin motors

The processive movement of single-headed kinesins is studied by using a ratchet model of non-Markov process, which is built on the experimental evidence that the strong binding of kinesin to microtubule in rigor state induces a large apparent change in the local microtubule conformation. In the model, the microtubule plays a crucial active role in the kinesin movement, in contrast to the previous belief that the microtubule only acts as a passive track for the kinesin motility. The unidirectional movement of single-headed kinesin is resulted from the asymmetric periodic potential between kinesin and microtubule while its processivity is determined by its binding affinity for microtubule in the weak ADP state. Using the model, various experimental results for monomeric kinesin KIF1A, such as the mean step size, the step-size distribution, the long run length and the mean velocity versus load, can be well explained quantitatively. This local conformational change of the microtubule may also play important roles in the processive movement of conventional two-headed kinesins. An experiment to verify the model is suggested.