Relating molecular flexibility to function: a case study of tubulin

Tubulin Vibration Modes Are in the Subterahertz Range, and Their Electromagnetic Absorption Is Affected by Water

Many proteins are thought to coordinate distant sites in their structures through a concerted action of global structural vibrations. However, the direct experimental spectroscopic detection of these vibration modes is rather elusive. We used normal-mode analysis to explore the dominant vibration modes of an all-atom model of the tubulin protein and described their characteristics using a large ensemble of tubulin structures. We quantified the frequency range of the normal vibrational modes to be in the subterahertz band, specifically between ∼40 and ∼160 GHz. Adding water layers to the model increases the frequencies of the low-frequency modes and narrows the frequency variations of the modes among the protein ensemble. We also showed how the electromagnetic absorption of tubulin vibration modes is affected by vibrational damping. These results contribute to our understanding of tubulin's vibrational and electromagnetic properties and provide a foundation for future attempts to control protein behavior via external electromagnetic fields.

Raman Spectroscopy Reveals Photobiomodulation-Induced α-Helix to β-Sheet Transition in Tubulins: Potential Implications for Alzheimer's and Other Neurodegenerative Diseases

In small clinical studies, the application of transcranial photobiomodulation (PBM), which typically delivers low-intensity near-infrared (NIR) to treat the brain, has led to some remarkable results in the treatment of dementia and several neurodegenerative diseases. However, despite the extensive literature detailing the mechanisms of action underlying PBM outcomes, the specific mechanisms affecting neurodegenerative diseases are not entirely clear. While large clinical trials are warranted to validate these findings, evidence of the mechanisms can explain and thus provide credible support for PBM as a potential treatment for these diseases. Tubulin and its polymerized state of microtubules have been known to play important roles in the pathology of Alzheimer's and other neurodegenerative diseases. Thus, we investigated the effects of PBM on these cellular structures in the quest for insights into the underlying therapeutic mechanisms. In this study, we employed a Raman spectroscopic analysis of the amide I band of polymerized samples of tubulin exposed to pulsed low-intensity NIR radiation (810 nm, 10 Hz, 22.5 J/cm2 dose). Peaks in the Raman fingerprint region (300-1900 cm-1)-in particular, in the amide I band (1600-1700 cm-1)-were used to quantify the percentage of protein secondary structures. Under this band, hidden signals of C=O stretching, belonging to different structures, are superimposed, producing a complex signal as a result. An accurate decomposition of the amide I band is therefore required for the reliable analysis of the conformation of proteins, which we achieved through a straightforward method employing a Voigt profile. This approach was validated through secondary structure analyses of unexposed control samples, for which comparisons with other values available in the literature could be conducted. Subsequently, using this validated method, we present novel findings of statistically significant alterations in the secondary structures of polymerized NIR-exposed tubulin, characterized by a notable decrease in α-helix content and a concurrent increase in β-sheets compared to the control samples. This PBM-induced α-helix to β-sheet transition connects to reduced microtubule stability and the introduction of dynamism to allow for the remodeling and, consequently, refreshing of microtubule structures. This newly discovered mechanism could have implications for reducing the risks associated with brain aging, including neurodegenerative diseases like Alzheimer's disease, through the introduction of an intervention following this transition.

Phosphate uptake in PhoX: Molecular mechanisms

PhoX is a high-affinity phosphate binding protein, present in Xanthomonas citri, a phytopathogen responsible for the citrus canker disease. Performing molecular dynamics simulations and different types of computational analyses, we study the molecular mechanisms at play in relation to phosphate binding, revealing the global functioning of the protein: PhoX naturally oscillates along its global normal modes, which allow it to explore both bound and unbound conformations, eventually attracting a nearby negative phosphate ion to the highly positive electrostatic potential on its surface, particularly close to the binding pocket. There, several hydrogen bonds are formed with the two main domains of the structure. Phosphate creates, in this way, a strong bridge that connects the domains, keeping itself between them, in a tight closed conformation, explaining its high binding affinity.

Molecular Dynamics and Machine Learning reveal distinguishing mechanisms of Competitive Ligands to perturb α,β-Tubulin

The mechanisms of action of ligands competing for the Colchicine Binding Site (CBS) of the α,β-Tubulin are non-standard compared to the commonly witnessed ligand-induced inhibition of proteins. This is because their potencies are not solely judged by the binding affinity itself, but also by their capacity to bias the conformational states of the dimer. Regarding the latter requirement, it is observed that ligands competing for the same pocket that binds colchicine exhibit divergence in potential clinical outcomes. Molecular dynamics-based ∼5.2 µs sampling of α,β-Tubulin complexed with four different ligands has revealed that each ligand has its customized way of influencing the complex. Primarily, it is the proportion of twisting and/or bending characteristic of modes of the intrinsic dynamics which is revealed to be 'fundamental' to tune this variation in the mechanism. The milder influence of 'bending' makes a ligand (TUB092), better classifiable under the group of vascular disrupting agents (VDAs), which are phenotypically effective on cytoskeletons; whereas a stronger impact of 'bending' makes the classical ligand Colchicine (COL) a better Anti-Mitotic agent (AMA). Two other ligands BAL27862 (2RR) and Nocodazole (NZO) fall in the intermediate zone as they fail to explicitly induce bending modes. Random Forest Classification method and K-means Clustering is applied to reveal the efficiency of Machine Learning methods in classifying the Tubulin conformations according to their ligand-specific perturbations and to highlight the significance of specific amino acid residues, mostly positioned in the α-β and β-β interfaces involved in the mechanism. These key residues responsible to yield discriminative actions of the ligands are likely to be highly useful in future endeavours to design more precise inhibitors.

Computational analysis of long-range allosteric communications in CFTR

Malfunction of the CFTR protein results in cystic fibrosis, one of the most common hereditary diseases. CFTR functions as an anion channel, the gating of which is controlled by long-range allosteric communications. Allostery also has direct bearings on CF treatment: the most effective CFTR drugs modulate its activity allosterically. Herein, we integrated Gaussian network model, transfer entropy, and anisotropic normal mode-Langevin dynamics and investigated the allosteric communications network of CFTR. The results are in remarkable agreement with experimental observations and mutational analysis and provide extensive novel insight. We identified residues that serve as pivotal allosteric sources and transducers, many of which correspond to disease-causing mutations. We find that in the ATP-free form, dynamic fluctuations of the residues that comprise the ATP-binding sites facilitate the initial binding of the nucleotide. Subsequent binding of ATP then brings to the fore and focuses on dynamic fluctuations that were present in a latent and diffuse form in the absence of ATP. We demonstrate that drugs that potentiate CFTR's conductance do so not by directly acting on the gating residues, but rather by mimicking the allosteric signal sent by the ATP-binding sites. We have also uncovered a previously undiscovered allosteric 'hotspot' located proximal to the docking site of the phosphorylated regulatory (R) domain, thereby establishing a molecular foundation for its phosphorylation-dependent excitatory role. This study unveils the molecular underpinnings of allosteric connectivity within CFTR and highlights a novel allosteric 'hotspot' that could serve as a promising target for the development of novel therapeutic interventions.

Structural basis of microtubule depolymerization by the kinesin-like activity of HIV-1 Rev

HIV-1 Rev is an essential regulatory protein that transports unspliced and partially spliced viral mRNAs from the nucleus to the cytoplasm for the expression of viral structural proteins. During its nucleocytoplasmic shuttling, Rev interacts with several host proteins to use the cellular machinery for the advantage of the virus. Here, we report the 3.5 Å cryo-EM structure of a 4.8 MDa Rev-tubulin ring complex. Our structure shows that Rev's arginine-rich motif (ARM) binds to both the acidic surfaces and the C-terminal tails of α/β-tubulin. The Rev-tubulin interaction is functionally homologous to that of kinesin-13, potently destabilizing microtubules at sub-stoichiometric levels. Expression of Rev in astrocytes and HeLa cells shows that it can modulate the microtubule cytoskeleton within the cellular environment. These results show a previously undefined regulatory role of Rev.

Structural impact of pathogenic SNPs on β-tubulin using molecular dynamics study

Single nucleotide polymorphisms (SNPs) in the TUBB1 (β-tubulin) gene have been implicated as the primary cause of macro thrombocytopenia. Therefore it is essential to identify the potential SNPs which are harmful to cause diseases such as macro thrombocytopenia. The impact caused by these variants on β-tubulin is twofold, both structural and functional. Multiple in-silico tools were used to scrutinise the most deleterious nsSNPs (non-synonymous SNPs) via sequence and structure-based approaches. Further, the β-tubulin protein model incorporating identified mutants was subjected to MD (molecular dynamic) simulations to analyse the impact on protein structure. A total of 2974 SNPs of TUBB1 were retrieved from various sources, and 32 nsSNPs were identified. By screening through sequence-based technique, 13 variants were detected as deleterious and further structure-based filtration was carried out to find thermally destabilising variants. Finally, three variants have been detected as highly destabilising by the mCSM server and chosen for the MD study. All three variants are present in the N-terminal, Intermediate, and C-terminal regions, breaking the spatial arrangement required for microtubule assembly. The spatial arrangement of these variants is in deviation with respect to WT (wild type) β-tubulin. The protein model was subjected to a simulation period of 100 ns. The FEL analysis revealed multiple clusters with minor populations indicating the unstable conformation adapted by the β-tubulin. The normal mode vector analysis exhibited high-intensity flexible motions at the C-terminal end, responsible for binding with MAPs (microtubule-associated proteins), an essential region in microtubule assembly. All these results reveal that the SNP's predicted eventually influence the spatial arrangement of β-tubulin, which would disturb the stacking arrangement of αβ tubulin dimer in microtubule assembly. The present study may set a path to cure the diseases like macro thrombocytopenia.Communicated by Ramaswamy H. Sarma.

Molecular basis of Toxoplasma gondii oryzalin resistance from a novel α-tubulin binding site model

Oryzalin (ORY) is a dinitroaniline derivative that inhibits the microtubule polymerization in plants and parasitic protozoa by selectively binding to the α-tubulin subunit. This herbicidal agent exhibits good antiprotozoal activity against major human parasites, such as Toxoplasma gondii (toxoplasmosis), Leishmania mexicana (leishmaniasis), and Plasmodium falciparum (malaria). Previous chemical mutagenesis assays on T. gondii α-tubulin (TgAT) have identified key mutations that lead to ORY resistance. Herein, we employed alchemical free energy methods and molecular dynamics simulations to determine if the ORY resistance mutations either decrease the TgAT's affinity of the compound or increase the protein stability. Our results here suggest that L136F and V202F mutations significantly decrease the affinity of ORY to TgAT, while T239I and V252L mutations diminish TgAT's flexibility. On the other hand, protein stability predictors determined that R243S mutation reduces TgAT stability due to the loss of its salt bridge interaction with E27. Interestingly, molecular dynamics simulations confirm that the loss of this key interaction leads to ORY binding site closure. Our study provides a better insight into the TgAT-ORY interaction, further supporting our recently proposed ORY-binding site.

Molecular Evolution of Tubulins in Diatoms

Microtubules are formed by α- and β-tubulin heterodimers nucleated with γ-tubulin. Tubulins are conserved eukaryotic proteins. Previously, it was shown that microtubules are involved in diatom silica frustule morphogenesis. Diatom frustules are varied, and their morphology is species-specific. Despite the attractiveness of the problem of elucidating the molecular mechanisms of genetically programmed morphogenesis, the structure and evolution of diatom tubulins have not been studied previously. Based on available genomic and transcriptome data, we analyzed the phylogeny of the predicted amino acid sequences of diatom α-, β- and γ-tubulins and identified five groups for α-tubulins, six for β-tubulins and four for γ-tubulins. We identified characteristic amino acids of each of these groups and also analyzed possible posttranslational modification sites of diatom tubulins. According to our results, we assumed what changes occurred in the diatom tubulin structures during their evolution. We also identified which tubulin groups are inherent in large diatom taxa. The similarity between the evolution of diatom tubulins and the evolution of diatoms suggests that molecular changes in α-, β- and γ-tubulins could be one of the factors in the formation of a high morphological diversity of diatoms.

Mechanisms of influence of the microtubule over-stabilizing ligands on the structure and intrinsic dynamics of α,β-Tubulin

The intervention into the cell cycle progression by administering microtubule over-stabilizing ligands that arrest the mitotic cell division by preventing spindle dissociation, is a promising strategy to fight against cancers. The building blocks of the microtubules and the spindles, i.e. the α,β-tubulin dimer, upon binding of such ligands, stay more comfortably in the microtubular multimeric form; the phenomenon of which is the key to the said over-stabilization. Using two such over-stabilizing ligands, Taxol and Taxotere, the present work reports the collective changes that these ligands induce on the structure and dynamics of the α,β-tubulin dimer which could be reconciled as the molecular basis of the over-stabilization of the microtubules; the trends have been found to be statistically significant across all independent calculations on them. The ligand binding increases the coherence between the residue communities of the two opposite faces of the β-subunit, which in a periodic arrangement in microtubule are knwon to form intermolecular contact with each other. This is likely to create an indirect cooperativity between those structural regions and this is a consequence of the reshuffling of the internal network of interactions upon ligand binding. Such reorganizations are also complemented by the increased contributions of the softer modes of the intrinsic dynamics more, which is likely to increase the plasticity of the system favourable for making structural adjustments in a multimer. Further, the ligands are able to compensate the drawback of lacking one phosphate group in protein-GDP interactions compared to the same for protein-GTP and this is in agreement with the hints form the earlier reports. The findings form a mechanistic basis of the enhanced capacity of the α,β-tubulin dimer to get more favourably accommodated into the microtubule superstructure upon binding either of Taxol and Taxotere.

Atomistic molecular dynamics simulations of tubulin heterodimers explain the motion of a microtubule

Microtubules are essential parts of the cytoskeleton that are built by polymerization of tubulin heterodimers into a hollow tube. Regardless that their structures and functions have been comprehensively investigated in a modern soft matter, it is unclear how properties of tubulin heterodimer influence and promote the self-assembly. A detailed knowledge of such structural mechanisms would be helpful in drug design against neurodegenerative diseases, cancer, diabetes etc. In this work atomistic molecular dynamics simulations were used to investigate the fundamental dynamics of tubulin heterodimers in a sheet and a short microtubule utilizing well-equilibrated structures. The breathing motions of the tubulin heterodimers during assembly show that the movement at the lateral interface between heterodimers (wobbling) dominates in the lattice. The simulations of the protofilament curvature agrees well with recently published experimental data, showing curved protofilaments at polymerization of the microtubule plus end. The tubulin heterodimers exposed at the microtubule minus end were less curved and displayed altered interactions at the site of sheet closure around the outmost heterodimers, which may slow heterodimer binding and polymerization, providing a potential explanation for the limited dynamics observed at the minus end.

The Assembly Switch Mechanism of FtsZ Filament Revealed by All-Atom Molecular Dynamics Simulations and Coarse-Grained Models

Bacterial cytoskeletal protein FtsZ binds and hydrolyzes GTP, and assembles into dynamic filaments that are essential for cell division. Here, we used a multi-scale computational strategy that combined all-atom molecular dynamics (MD) simulations and coarse-grained models to reveal the conformational dynamics of assembled FtsZ. We found that the top end of a filament is highly dynamic and can undergo T-to-R transitions in both GTP- and GDP-bound states. We observed several subcategories of nucleation related dimer species, which leading to a feasible nucleation pathway. In addition, we observed that FtsZ filament exhibits noticeable amounts of twisting, indicating a substantial helicity of the FtsZ filament. These results agree with the previously models and experimental data. Anisotropy network model (ANM) analysis revealed a polymerization enhanced assembly cooperativity, and indicated that the cooperative motions in FtsZ are encoded in the structure. Taken together, our study provides a molecular-level understanding of the diversity of the structural states of FtsZ and the relationships among polymerization, hydrolysis, and cooperative assembly, which should shed new light on the molecular basis of FtsZ’s cooperativity.

Computational Ways to Enhance Protein Inhibitor Design

Two new computational approaches are described to aid in the design of new peptide-based drugs by evaluating ensembles of protein structures from their dynamics and through the assessing of structures using empirical contact potential. These approaches build on the concept that conformational variability can aid in the binding process and, for disordered proteins, can even facilitate the binding of more diverse ligands. This latter consideration indicates that such a design process should be less restrictive so that multiple inhibitors might be effective. The example chosen here focuses on proteins/peptides that bind to hemagglutinin (HA) to block the large-scale conformational change for activation. Variability in the conformations is considered from sets of experimental structures, or as an alternative, from their simple computed dynamics; the set of designe peptides/small proteins from the David Baker lab designed to bind to hemagglutinin, is the large set considered and is assessed with the new empirical contact potentials.

Computational Modeling of Environmental Co-exposure on Oil-Derived Hydrocarbon Overload by Using Substrate-Specific Transport Protein (TodX) with Graphene Nanostructures

BACKGROUND

Bioremediation is a biotechnology field that uses living organisms to remove contaminants from soil and water; therefore, they could be used to treat oil spills from the environment.

METHODS

Herein, we present a new mechanistic approach combining Molecular Docking Simulation and Density Functional Theory to modeling the bioremediation-based nanointeractions of a heterogeneous mixture of oil-derived hydrocarbons by using pristine and oxidized graphene nanostructures and the substrate-specific transport protein (TodX) from Pseudomonas putida.

RESULTS

The theoretical evidences pointing that the binding interactions are mainly based on noncovalent bonds characteristic of physical adsorption mechanism mimicking the "Trojan-horse effect".

CONCLUSION

These results open new horizons to improve bioremediation strategies in over-saturation conditions against oil-spills and expanding the use of nanotechnologies in the context of environmental modeling health and safety.

Mapping the underlying mechanisms of fibrinogen benzothiazole drug interactions using computational and experimental approaches

Three-dimensional conformational crystallographic binding-modes are of paramount importance to understand the docking mechanism of protein-ligand interactions and to identify potential "leading drugs" conformers towards rational drugs-design. Herein, we present an integrated computational-experimental study tackling the problem of multiple binding modes among the ligand 3-(2-Benzothiazolylthio)-propane sulfonic acid (BTS) and the fibrinogen receptor (E-region). Based on molecular docking simulations, we found that the free energy of binding values for nine of different BTS-docking complexes (i.e., BTS-pose_1-9) were very close. We have also identified a docking-mechanism of BTS-interaction mainly based on non-covalent hydrophobic interactions with H-bond contacts stabilizing the fibrinogen-BTS docking complexes. Interestingly, the different BTS-poses_1-9 were found to be able to block the fibrinogen binding site (E-region) by inducing local perturbations in effector and allosteric residues, reducing the degree of collectivity in its flexibility normal modes. As such, we theoretically suggest that the BTS-binding modes can significantly affect the physiological condition of the unoccupied fibrinogen protein structure by bringing global and local perturbations in the frequency domain spectra. The proposed theoretical mechanisms, the interactions involved and the conformational changes suggested, were further corroborated by different experimental techniques such as isothermal titration calorimetry (ITC), zeta potential, UV-vis, fluorescence and small angle X-ray scattering (SAXS). The combined results shall open new avenues towards the application of complex supra-molecular information in rational drugs-design.

When Stiffness Matters: Mechanosensing in Heart Development and Disease

During embryonic morphogenesis, the heart undergoes a complex series of cellular phenotypic maturations (e.g., transition of myocytes from proliferative to quiescent or maturation of the contractile apparatus), and this involves stiffening of the extracellular matrix (ECM) acting in concert with morphogenetic signals. The maladaptive remodeling of the myocardium, one of the processes involved in determination of heart failure, also involves mechanical cues, with a progressive stiffening of the tissue that produces cellular mechanical damage, inflammation, and ultimately myocardial fibrosis. The assessment of the biomechanical dependence of the molecular machinery (in myocardial and non-myocardial cells) is therefore essential to contextualize the maturation of the cardiac tissue at early stages and understand its pathologic evolution in aging. Because systems to perform multiscale modeling of cellular and tissue mechanics have been developed, it appears particularly novel to design integrated mechano-molecular models of heart development and disease to be tested in ex vivo reconstituted cells/tissue-mimicking conditions. In the present contribution, we will discuss the latest implication of mechanosensing in heart development and pathology, describe the most recent models of cell/tissue mechanics, and delineate novel strategies to target the consequences of heart failure with personalized approaches based on tissue engineering and induced pluripotent stem cell (iPSC) technologies.

Mechanical properties of tubulin intra- and inter-dimer interfaces and their implications for microtubule dynamic instability

Thirteen tubulin protofilaments, made of αβ-tubulin heterodimers, interact laterally to produce cytoskeletal microtubules. Microtubules exhibit the striking property of dynamic instability, manifested in their intermittent growth and shrinkage at both ends. This behavior is key to many cellular processes, such as cell division, migration, maintenance of cell shape, etc. Although assembly and disassembly of microtubules is known to be linked to hydrolysis of a guanosine triphosphate molecule in the pocket of β-tubulin, detailed mechanistic understanding of corresponding conformational changes is still lacking. Here we take advantage of the recent generation of in-microtubule structures of tubulin to examine the properties of protofilaments, which serve as important microtubule assembly and disassembly intermediates. We find that initially straight tubulin protofilaments, relax to similar non-radially curved and slightly twisted conformations. Our analysis further suggests that guanosine triphosphate hydrolysis primarily affects the flexibility and conformation of the inter-dimer interface, without a strong impact on the shape or flexibility of αβ-heterodimer. Inter-dimer interfaces are significantly more flexible compared to intra-dimer interfaces. We argue that such a difference in flexibility could be key for distinct stability of the plus and minus microtubule ends. The higher flexibility of the inter-dimer interface may have implications for development of pulling force by curving tubulin protofilaments during microtubule disassembly, a process of major importance for chromosome motions in mitosis.

Structural and energetic evolution of fibrinogen toward to the betablocker interactions

We present a computational analysis coupled with experimental studies, focusing on the binding-interaction between beta-adrenoreceptor blocking agents (acebutolol and propranolol) with fibrinogen protein (E-region). Herein, computational modeling on structural validation and flexibility properties of fibrinogen E-region showed that the E-region interacting residues, which form the funnel-shaped hydrophobic cavity for ligand-binding, can be efficiently modeled. The obtained free energy of binding (FEB) values for the docking complexes, namely acebutolol/fibrinogen E-region and propranolol/fibrinogen E-region, were very close and amounted to - 6.9 kcal/mol and - 6.8 kcal/mol, respectively. They were supported by a high binding-accuracy (R.M.S.D < 2 Å) for the best crystallographic binding-poses in both cases. In this regard, we identify a docking-mechanism of interaction for the propranolol and acebutolol mainly based on non-covalent hydrophobic contacts with the fibrinogen E-region binding-site. Besides, the beta-adrenoreceptor blocking agents are able to induce local perturbations affecting particularly the fibrinogen E-region allosteric residues linked to significant changes in the inter-residue communication and flexibility properties of residue network. In this sense, we show that the key biophysical parameters like frequency and collectivity degree may be compromised in different ways by the interaction with acebutolol and propranolol. Isothermal titration calorimetry, zeta potential and small angle X-ray scattering (SAXS) measurements were performed to complete and corroborate computational analysis. The combined experimental results point out that acebutolol acts to a lesser extent to fibrinogen structure than propranolol.

Tubulin's response to external electric fields by molecular dynamics simulations

Tubulin heterodimers are the building blocks of microtubules and disruption of their dynamics is exploited in the treatment of cancer. Electric fields at certain frequencies and magnitudes are believed to do the same. Here, the tubulin dimer’s response to external electric fields was determined by atomistic simulation. External fields from 50 to 750 kV/cm, applied for 10 ns, caused significant conformational rearrangements that were dependent upon the field’s directionality. Charged and flexible regions, including the α:H1-B2 loop, β:M-loop, and C-termini, were susceptible. Closer inspection of the α:H1-B2 loop in lower strength fields revealed that these effects were consistent and proportional to field strength, and the findings indicate that external electric fields modulate the stability of microtubules through conformational changes to key loops involved in lateral contacts. We also find evidence that tubulin’s curvature and elongation are affected, and external electric fields may bias tubulin towards depolymerization.

Combining multi-scale modelling methods to decipher molecular motions of a branching sucrase from glycoside-hydrolase family 70

Among α-transglucosylases from Glycoside-Hydrolase family 70, the ΔN₁₂₃-GB-CD2 enzyme derived from the bifunctional DSR-E from L. citreum NRRL B-1299 is particularly interesting as it was the first described engineered Branching Sucrase, not able to elongate glucan polymers from sucrose substrate. The previously reported overall structural organization of this multi-domain enzyme is an intricate U-shape fold conserved among GH70 enzymes which showed a certain conformational variability of the so-called domain V, assumed to play a role in the control of product structures, in available X-ray structures. Understanding the role of functional dynamics on enzyme reaction and substrate recognition is of utmost interest although it remains a challenge for biophysical methods. By combining long molecular dynamics simulation (1μs) and multiple analyses (NMA, PCA, Morelet Continuous Wavelet Transform and Cross Correlations Dynamics), we investigated here the dynamics of ΔN₁₂₃-GB-CD2 alone and in interaction with sucrose substrate. Overall, our results provide the detailed picture at atomic level of the hierarchy of motions occurring along different timescales and how they are correlated, in agreement with experimental structural data. In particular, detailed analysis of the different structural domains revealed cooperative dynamic behaviors such as twisting, bending and wobbling through anti- and correlated motions, and also two structural hinge regions, of which one was unreported. Several highly flexible loops surrounding the catalytic pocket were also highlighted, suggesting a potential role in the acceptor promiscuity of ΔN123-GBD-CD2. Normal modes and essential dynamics underlined an interesting two-fold dynamic of the catalytic domain A, pivoting about an axis splitting the catalytic gorge in two parts. The comparison of the conformational free energy landscapes using principal component analysis of the enzyme in absence or in presence of sucrose, also revealed a more harmonic basin when sucrose is bound with a shift population of the bending mode, consistent with the substrate binding event.

Effect of Nucleotide State on the Protofilament Conformation of Tubulin Octamers

At the molecular level, the dynamic instability (random growth and shrinkage) of the microtubule (MT) is driven by the nucleotide state (GTP vs GDP) in the β subunit of the tubulin dimers at the MT cap. Here, we use large-scale molecular dynamics (MD) simulations and normal-mode analysis (NMA) to characterize the effect of a single GTP cap layer on tubulin octamers composed of two neighboring protofilaments (PFs). We utilize recently reported high-resolution structures of dynamic MTs to simulate a GDP octamer both with and without a single GTP cap layer. We perform multiple replicas of long-time atomistic MD simulations (3 replicas, 0.3 μs for each replica, 0.9 μs for each octamer system, and 1.8 μs total) of both octamers. We observe that a single GTP cap layer induces structural differences in neighboring PFs, finding that one PF possesses a gradual curvature, compared to the second PF which possesses a kinked conformation. This results in either curling or splaying between these PFs. We suggest that this is due to asymmetric strengths of longitudinal contacts between the two PFs. Furthermore, using NMA, we calculate mechanical properties of these octamer systems and find that octamer system with a single GTP cap layer possesses a lower flexural rigidity.

The Roles of β-Tubulin Mutations and Isotype Expression in Acquired Drug Resistance

The antitumor drug paclitaxel stabilizes microtubules and reduces their dynamicity, promoting mitotic arrest and eventually apoptosis. Upon assembly of the α/β-tubulin heterodimer, GTP becomes bound to both the α and β-tubulin monomers. During microtubule assembly, the GTP bound to β-tubulin is hydrolyzed to GDP, eventually reaching steady-state equilibrium between free tubulin dimers and those polymerized into microtubules. Tubulin-binding drugs such as paclitaxel interact with β-tubulin, resulting in the disruption of this equilibrium. In spite of several crystal structures of tubulin, there is little biochemical insight into the mechanism by which anti-tubulin drugs target microtubules and alter their normal behavior. The mechanism of drug action is further complicated, as the description of altered β-tubulin isotype expression and/or mutations in tubulin genes may lead to drug resistance as has been described in the literature. Because of the relationship between β-tubulin isotype expression and mutations within β-tubulin, both leading to resistance, we examined the properties of altered residues within the taxane, colchicine and Vinca binding sites. The amount of data now available, allows us to investigate common patterns that lead to microtubule disruption and may provide a guide to the rational design of novel compounds that can inhibit microtubule dynamics for specific tubulin isotypes or, indeed resistant cell lines. Because of the vast amount of data published to date, we will only provide a broad overview of the mutational results and how these correlate with differences between tubulin isotypes. We also note that clinical studies describe a number of predictive factors for the response to anti-tubulin drugs and attempt to develop an understanding of the features within tubulin that may help explain how they may affect both microtubule assembly and stability.

Deformation pattern in vibrating microtubule: Structural mechanics study based on an atomistic approach

The mechanical properties of microtubules are of great importance for understanding their biological function and for applications in artificial devices. Although microtubule mechanics has been extensively studied both theoretically and experimentally, the relation to its molecular structure is understood only partially. Here, we report on the structural analysis of microtubule vibration modes calculated by an atomistic approach. Molecular dynamics was applied to refine the atomic structure of a microtubule and a C_α elastic network model was analyzed for its normal modes. We mapped fluctuations and local deformations up to the level of individual aminoacid residues. The deformation is mode-shape dependent and principally different in α-tubulins and β-tubulins. Parts of the tubulin dimer sequence responding specifically to longitudinal and radial stress are identified. We show that substantial strain within a microtubule is located both in the regions of contact between adjacent dimers and in the body of tubulins. Our results provide supportive evidence for the generally accepted assumption that the mechanics of microtubules, including its anisotropy, is determined by the bonds between tubulins.

Molecular modeling and molecular dynamic simulation of the effects of variants in the TGFBR2 kinase domain as a paradigm for interpretation of variants obtained by next generation sequencing

Variants in the TGFBR2 kinase domain cause several human diseases and can increase propensity for cancer. The widespread application of next generation sequencing within the setting of Individualized Medicine (IM) is increasing the rate at which TGFBR2 kinase domain variants are being identified. However, their clinical relevance is often uncertain. Consequently, we sought to evaluate the use of molecular modeling and molecular dynamics (MD) simulations for assessing the potential impact of variants within this domain. We documented the structural differences revealed by these models across 57 variants using independent MD simulations for each. Our simulations revealed various mechanisms by which variants may lead to functional alteration; some are revealed energetically, while others structurally or dynamically. We found that the ATP binding site and activation loop dynamics may be affected by variants at positions throughout the structure. This prediction cannot be made from the linear sequence alone. We present our structure-based analyses alongside those obtained using several commonly used genomics-based predictive algorithms. We believe the further mechanistic information revealed by molecular modeling will be useful in guiding the examination of clinically observed variants throughout the exome, as well as those likely to be discovered in the near future by clinical tests leveraging next-generation sequencing through IM efforts.

Bending-Twisting Motions and Main Interactions in Nucleoplasmin Nuclear Import

Alpha solenoid proteins play a key role in regulating the classical nuclear import pathway, recognizing a target protein and transporting it into the nucleus. Importin-α (Impα) is the solenoid responsible for cargo protein recognition, and it has been extensively studied by X-ray crystallography to understand the binding specificity. To comprehend the main motions of Impα and to extend the information about the critical interactions during carrier-cargo recognition, we surveyed different conformational states based on molecular dynamics (MD) and normal mode (NM) analyses. Our model of study was a crystallographic structure of Impα complexed with the classical nuclear localization sequence (cNLS) from nucleoplasmin (Npl), which was submitted to multiple 100 ns of MD simulations. Representative conformations were selected for calculating the 87 lowest frequencies NMs of vibration, and a displacement approach was applied along each NM. Based on geometric criteria, using the radius of curvature and inter-repeat angles as the reference metrics, the main motions of Impα were described. Moreover, we determined the salt bridges, hydrogen bonds and hydrophobic interactions in the Impα-NplNLS interface. Our results show the bending and twisting motions participating in the recognition of nuclear proteins, allowing the accommodation and adjustment of a classical bipartite NLS sequence. The essential contacts for the nuclear import were also described and were mostly in agreement with previous studies, suggesting that the residues in the cNLS linker region establish important contacts with Impα adjusting the cNLS backbone. The MD simulations combined with NM analysis can be applied to the Impα-NLS system to help understand interactions between Impα and cNLSs and the analysis of non-classic NLSs.

Relating nucleotide-dependent conformational changes in free tubulin dimers to tubulin assembly

The complex dynamic behavior of microtubules (MTs) is believed to be primarily due to the αβ-tubulin dimer architecture and its intrinsic GTPase activity. Hence, a detailed knowledge of the conformational variations of isolated α-GTP-β-GTP- and α-GTP-β-GDP-tubulin dimers in solution and their implications to interdimer interactions and stability is directly relevant to understand the MT dynamics. An attempt has been made here by combining molecular dynamics (MD) simulations and protein-protein docking studies that unravels key structural features of tubulin dimer in different nucleotide states and correlates their association to tubulin assembly. Results from simulations suggest that tubulin dimers and oligomers attain curved conformations in both GTP and GDP states. Results also indicate that the tubulin C-terminal domain and the nucleotide state are closely linked. Protein-protein docking in combination with MD simulations suggest that the GTP-tubulin dimers engage in relatively stronger interdimer interactions even though the interdimer interfaces are bent in both GTP and GDP tubulin complexes, providing valuable insights on in vitro finding that GTP-tubulin is a better assembly candidate than GDP-tubulin during the MT nucleation and elongation processes.

Computing the Amino Acid Specificity of Fluctuations in Biomolecular Systems

We developed a new amino acid specific method for the computation of spatial fluctuations of proteins around their native structures. We show the consistency with experimental values and the increased performance in comparison to an established model, based on statistical estimates for a set of test proteins. We apply the new method to HIV-1 protease in its wild-type form and to a V82F-I84V mutant that shows resistance to protease inhibitors. We further show how the method can be successfully used to explain the molecular biophysics of drug resistance of the mutant.

Elucidating the mechanism of action of the clinically approved taxanes: a comprehensive comparison of local and allosteric effects

The clinically approved taxanes (paclitaxel, docetaxel and cabazitaxel) target the tubulin protein in microtubules. Despite the clinical success of these agents, the mechanism of action of this class of drugs remains elusive, making rational design of taxanes difficult. Molecular dynamics simulations of these three taxanes with the αβ-tubulin heterodimer examine the similarities and differences in the effects of the drugs on tubulin, probing both local and allosteric effects. Despite their structural similarity, the drugs adopt different conformations in the binding site on β-tubulin. The taxanes similarly increase the helical character of α- and β-tubulins. No correlations are found between microtubule assembly and (i) binding affinity or (ii) the role of the M-loop in enhancing lateral contacts. Instead, changes in intra- and interdimer longitudinal contacts are indicative of the mechanism of action of the taxanes. We find β:H1-S1', and more importantly β:H9 and β:H10, play a role translating the effect of local drug binding in β-tubulin to an allosteric effect in α-tubulin and propose that the displacement of these secondary structures towards α-tubulin may be used as a predictor of the effect of taxanes on the tubulin heterodimers in rational drug design approaches.

Seminal role of deletion of amino acid residues in H1-S2 and S-loop regions in eukaryotic β-tubulin investigated from docking and dynamics perspective

Tubulin is the fundamental unit of microtubules. It is reported to effect different functions like cell division, chromosomal segregation, motility and intracellular transportation. α- and β-tubulin associate laterally and longitudinally to form protofilaments. Both the subunits are structurally identical to each other except for the deletions reported in H1-S2 and S loop regions in eukaryotic β-tubulin. These deletions mimic the ancestral tubulin protein named Latest Common FtsZ-Tubulin Ancestor (LCFTA) with a shorter S-loop region resulting in weak dimerization. However, in eukaryotic beta tubulin, the significance of this shorter region remains elusive till date. The main objective of this study was to model variants of beta tubulin (βmut1, βmut2 and βmut3) with inserts that lengthened the loop, and to compare them with the native α- and β-subunits to understand their biological significance. Further, one more mutant was modeled with the intention of understanding the counter effect of additional deletion of amino acid residues from both H1-S2 and S-loop regions; this mutant was designated as βmut4. Our study confirms that the insertion of amino acid residues considerably increases the protein-protein interactions in βmut1-βmut1, βmut2-βmut2 and βmut3-βmut3 compared to their native β-subunit. Similarly, the binding affinity of GTP also increases in βmut2 and βmut3 as compared to the wild type. However, these deletions result in decreased protein-protein and ligand interactions in wild beta tubulin and βmut4, as compared to βmut1, βmut2,and βmut3. Therefore, we conclude here that residual inserts in the H1-S2 and S loop sub segments bring about conformational changes in regions critically involved in lateral interactions and in the nucleotide binding site, thus altering the binding affinities between the dimers and the ligands. Regarding the biological importance of such deletions in wild beta tubulin, these deletions result in flexible M-loop leading to weak protein-protein interaction. This could be an adaptive feature playing a crucial role in protofilament dissociation during GTP hydrolysis, because of weak dimerization.

Microtubules dual chemo and thermo-responsive depolymerization

The effects of chemotherapeutic agent vinblastine versus low temperature of 277 K were investigated on the structure of αβ-tubulin heterodimer by means of molecular dynamics simulations. Individual experiments have shown that the vinblastine-bound heterodimer, and its apo structure under low temperature of 277 K, both undergo conformational changes toward destabilization of the dimer as compared to the apo tubulin at 300 K. Both factors exhibited weakening of the longitudinal interactions of tubulin heterodimer through displacing dimer interfacial segments, resulting in dominant electrostatic repulsion at the interface of the subunits. The two independent factors of temperature and anti-mitotic agent facilitate alteration of secondary structure in functional segments such as H1-S2 loop, H3, H10 helices, and T7 loop, which are known to be important in either longitudinal or lateral contacts among αβ-heterodimers in MTs protofilaments and their depolymerization mechanism.

Design, Synthesis and Cytotoxicity Evaluationof New 1,2-diaryl-4, 5, 6, 7-Tetrahydro-1H-benzo[d] Imidazolesas Tubulin Inhibitors

A new series of 1,2-diaryl-4,5,6,7-tetrahydro-1H-benzo[d]imidazoles, possessing trimethoxyphenyl pharmacophore, were synthesized to evaluate their biological activities as tubulin inhibitors. Cytotoxic activity of the synthesized compounds 7a-f was assessed against several human cancer cell lines, including MCF-7 (breast cancer cell), HEPG2 (liver hepatocellular cells), A549 (adenocarcinomic human alveolar basal epithelial cells), T47D (Human ductal breast epithelial tumor cell line) and fibroblast. According to our results, HEPG2 seems to be the most sensitive, while MCF7 was the most resistant cell line to the compounds. All the compounds expect 7b, possessed satisfactory activity against HEPG2 with mean IC50 values ranging from 15.60 to 43.81 µM.

Population shift of binding pocket size and dynamic correlation analysis shed new light on the anticooperative mechanism of PII protein

PII protein is one of the largest families of signal transduction proteins in archaea, bacteria, and plants, controlling key processes of nitrogen assimilation. An intriguing characteristic for many PII proteins is that the three ligand binding sites exhibit anticooperative allosteric regulation. In this work, PII protein from Synechococcus elongatus, a model for cyanobacteria and plant PII proteins, is utilized to reveal the anticooperative mechanism upon binding of 2-oxoglutarate (2-OG). To this end, a method is proposed to define the binding pocket size by identifying residues that contribute greatly to the binding of 2-OG. It is found that the anticooperativity is realized through population shift of the binding pocket size in an asymmetric manner. Furthermore, a new algorithm based on the dynamic correlation analysis is developed and utilized to discover residues that mediate the anticooperative process with high probability. It is surprising to find that the T-loop, which is believed to be responsible for mediating the binding of PII with its target proteins, also takes part in the intersubunit signal transduction process. Experimental results of PII variants further confirmed the influence of T-loop on the anticooperative regulation, especially on binding of the third 2-OG. These discoveries extend our understanding of the PII T-loop from being essential in versatile binding of target protein to signal-mediating in the anticooperative allosteric regulation.

Mechanics of severing for large microtubule complexes revealed by coarse-grained simulations

We investigate the mechanical behavior of microtubule (MT) protofilaments under the action of bending forces, ramped up linearly in time, to provide insight into the severing of MTs by microtubule associated proteins (MAPs). We used the self-organized polymer model which employs a coarse-grained description of the protein chain and ran Brownian dynamics simulations accelerated on graphics processing units that allow us to follow the dynamics of a MT system on experimental timescales. Our study focused on the role played in the MT depolymerization dynamics by the inter-tubulin contacts a protofilament experiences when embedded in the MT lattice, and the number of binding sites of MAPs on MTs. We found that proteins inducing breaking of MTs must have at least three attachment points on any tubulin dimer from an isolated protofilament. In contrast, two points of contact would suffice when dimers are located in an intact MT lattice, in accord with experimental findings on MT severing proteins. Our results show that confinement of a protofilament in the MT lattice leads to a drastic reduction in the energy required for the removal of tubulin dimers, due to the drastic reduction in entropy. We further showed that there are differences in the energetic requirements based on the location of the dimer to be removed by severing. Comparing the energy of tubulin dimers removal revealed by our simulations with the amount of energy resulting from one ATP hydrolysis, which is the source of energy for all MAPs, we provided strong evidence for the experimental finding that severing proteins do not bind uniformly along the MT wall.

The critical role of the loops of triosephosphate isomerase for its oligomerization, dynamics, and functionality

Triosephosphate isomerase (TIM) catalyzes the reaction to convert dihydroxyacetone phosphate into glyceraldehyde 3-phosphate, and vice versa. In most organisms, its functional oligomeric state is a homodimer; however, tetramer formation in hyperthermophiles is required for functional activity. The tetrameric TIM structure also provides added stability to the structure, enabling it to function at more extreme temperatures. We apply Principal Component Analysis to find that the TIM structure space is clearly divided into two groups--the open and the closed TIM structures. The distribution of the structures in the open set is much sparser than that in the closed set, showing a greater conformational diversity of the open structures. We also apply the Elastic Network Model to four different TIM structures--an engineered monomeric structure, a dimeric structure from a mesophile--Trypanosoma brucei, and two tetrameric structures from hyperthermophiles Thermotoga maritima and Pyrococcus woesei. We find that dimerization not only stabilizes the structures, it also enhances their functional dynamics. Moreover, tetramerization of the hyperthermophilic structures increases their functional loop dynamics, enabling them to function in the destabilizing environment of extreme temperatures. Computations also show that the functional loop motions, especially loops 6 and 7, are highly coordinated. In summary, our computations reveal the underlying mechanism of the allosteric regulation of the functional loops of the TIM structures, and show that tetramerization of the structure as found in the hyperthermophilic organisms is required to maintain the coordination of the functional loops at a level similar to that in the dimeric mesophilic structure.

The intricate relationship between microtubules and their associated motor proteins during axon growth and maintenance

The hallmarks of neurons are their slender axons which represent the longest cellular processes of animals and which act as the cables that electrically wire the brain, and the brain to the body. Axons extend along reproducible paths during development and regeneration, and they have to be maintained for the lifetime of an organism. Both axon extension and maintenance essentially depend on the microtubule (MT) cytoskeleton. For this, MTs organize into parallel bundles that are established through extension at the leading axon tips within growth cones, and these bundles then form the architectural backbones, as well as the highways for axonal transport essential for supply and intracellular communication. Axon transport over these enormous distances takes days or even weeks and is a substantial logistical challenge. It is performed by kinesins and dynein/dynactin, which are molecular motors that form close functional links to the MTs they walk along. The intricate machinery which regulates MT dynamics, axonal transport and the motors is essential for nervous system development and function, and its investigation has huge potential to bring urgently required progress in understanding the causes of many developmental and degenerative brain disorders. During the last years new explanations for the highly specific properties of axonal MTs and for their close functional links to motor proteins have emerged, and it has become increasingly clear that motors play active roles also in regulating axonal MT networks. Here, I will provide an overview of these new developments.

The interplay of structure and dynamics: insights from a survey of HIV-1 reverse transcriptase crystal structures

HIV-1 reverse transcriptase (RT) is a critical drug target for HIV treatment, and understanding the exact mechanisms of its function and inhibition would significantly accelerate the development of new anti-HIV drugs. It is well known that structure plays a critical role in protein function, but for RT, structural information has proven to be insufficient-despite enormous effort-to explain the mechanism of inhibition and drug resistance of non-nucleoside RT inhibitors. We hypothesize that the missing link is dynamics, information about the motions of the system. However, many of the techniques that give the best information about dynamics, such as solution nuclear magnetic resonance and molecular dynamics simulations, cannot be easily applied to a protein as large as RT. As an alternative, we combine elastic network modeling with simultaneous hierarchical clustering of structural and dynamic data. We present an extensive survey of the dynamics of RT bound to a variety of ligands and with a number of mutations, revealing a novel mechanism for drug resistance to non-nucleoside RT inhibitors. Hydrophobic core mutations restore active-state motion to multiple functionally significant regions of HIV-1 RT. This model arises out of a combination of structural and dynamic information, rather than exclusively from one or the other.

Multiple tubulins: evolutionary aspects and biological implications

Plant tubulin is a dimeric protein that contributes to formation of microtubules, major intracellular structures that are involved in the control of fundamental processes such as cell division, polarity of growth, cell-wall deposition, intracellular trafficking and communications. Because it is a structural protein whose function is confined to the role of microtubule formation, tubulin may be perceived as an uninteresting gene product, but such a perception is incorrect. In fact, tubulin represents a key molecule for studying fundamental biological issues such as (i) microtubule evolution (also with reference to prokaryotic precursors and the formation of cytomotive filaments), (ii) protein structure with reference to the various biochemical features of members of the FstZ/tubulin superfamily, (iii) isoform variations contributed by the existence of multi-gene families and various kinds of post-translational modifications, (iv) anti-mitotic drug interactions and mode of action, (v) plant and cell symmetry, as determined using a series of tubulin mutants, (vi) multiple and sophisticated mechanisms of gene regulation, and (vii) intron molecular evolution. In this review, we present and discuss many of these issues, and offer an updated interpretation of the multi-tubulin hypothesis.

A mass weighted chemical elastic network model elucidates closed form domain motions in proteins

An elastic network model (ENM), usually Cα coarse-grained one, has been widely used to study protein dynamics as an alternative to classical molecular dynamics simulation. This simple approach dramatically saves the computational cost, but sometimes fails to describe a feasible conformational change due to unrealistically excessive spring connections. To overcome this limitation, we propose a mass-weighted chemical elastic network model (MWCENM) in which the total mass of each residue is assumed to be concentrated on the representative alpha carbon atom and various stiffness values are precisely assigned according to the types of chemical interactions. We test MWCENM on several well-known proteins of which both closed and open conformations are available as well as three α-helix rich proteins. Their normal mode analysis reveals that MWCENM not only generates more plausible conformational changes, especially for closed forms of proteins, but also preserves protein secondary structures thus distinguishing MWCENM from traditional ENMs. In addition, MWCENM also reduces computational burden by using a more sparse stiffness matrix.

The effect of complex solvents on the structure and dynamics of protein solutions: The case of Lysozyme in trehalose/water mixtures

We present a Molecular Dynamics simulation study of the effect of trehalose concentration on the structure and dynamics of individual proteins immersed in trehalose/water mixtures. Hen egg-white Lysozyme is used in this study and trehalose concentrations of 0%, 10%, 20%, 30% and 100% by weight are explored. Surprisingly, we have found that changes in trehalose concentration do not change the global structural characteristics of the protein as measured by standard quantities like the mean square deviation, radius of gyration, solvent accessible surface area, inertia tensor and asphericity. Only in the limit of pure trehalose these metrics change significantly. Specifically, we found that the protein is compressed by 2% when immersed in pure trehalose. At the amino acid level there is noticeable rearrangement of the surface residues due to the change in polarity of the surrounding environment with the addition of trehalose. From a dynamic perspective, our computation of the Incoherent Intermediate Scattering Function shows that the protein slows down with increasing trehalose concentration; however, this slowdown is not monotonic. Finally, we also report in-depth results for the hydration layer around the protein including its structure, hydrogen-bonding characteristics and dynamic behavior at different length scales.