Tunnel-like region observed as a potential allosteric site in Staphylococcus aureus Glyceraldehyde-3-phosphate dehydrogenase.

Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) catalyzing the sixth step of glycolysis has been investigated for allosteric features that might be used as potential target for specific inhibition of Staphylococcus aureus (S.aureus). X-ray structure of bacterial enzyme for which a tunnel-like opening passing through the center previously proposed as an allosteric site has been subjected to six independent 500 ns long Molecular Dynamics simulations. Harmonic bond restraints were employed at key residues to underline the allosteric feature of this region. A noticeable reduction was observed in the mobility of NAD+ binding domains when restrictions were applied. Also, a substantial decrease in cross-correlations between distant Cα fluctuations was detected throughout the structure. Mutual information (MI) analysis revealed a similar decrease in the degree of correspondence in positional fluctuations in all directions everywhere in the receptor. MI between backbone and side-chain torsional variations changed its distribution profile and decreased considerably around the catalytic sites when restraints were employed. Principal component analysis clearly showed that the restrained state sampled a narrower range of conformations than apo state, especially in the first principal mode due to restriction in the conformational flexibility of NAD+ binding domain. Clustering the trajectory based on catalytic site residues displayed a smaller repertoire of conformations for restrained state compared to apo. Representative snapshots subjected to k-shortest pathway analysis revealed the impact of bond restraints on the allosteric communication which displayed distinct optimal and suboptimal pathways for two states, where observed frequencies of critical residues Gln51 and Val283 at the proposed site changed considerably.

Innovative Use of an Injectable, Self-Healing Drug-Loaded Pectin-Based Hydrogel for Micro- and Supermicro-Vascular Anastomoses.

Microvascular surgery plays a crucial role in reconnecting micrometer-scale vessel ends. Suturing remains the gold standard technique for small vessels; however, suturing the collapsed lumen of microvessels is challenging and time-consuming, with the risk of misplaced sutures leading to failure. Although multiple solutions have been reported, the emphasis has predominantly been on resolving challenges related to arteries rather than veins, and none has proven superior. In this study, we introduce an innovative solution to address these challenges through the development of an injectable lidocaine-loaded pectin hydrogel by using computational and experimental methods. To understand the extent of interactions between the drug and the pectin chain, molecular dynamics (MD) simulations and quantum mechanics (QM) calculations were conducted in the first step of the research. Then, a series of experimental studies were designed to prepare lidocaine-loaded injectable pectin-based hydrogels, and their characterization was performed by using Fourier transform infrared spectroscopy (FT-IR), scanning electron microscopy (SEM), and rheological analysis. After all the results were evaluated, the drug-loaded pectin-based hydrogel exhibiting self-healing properties was selected as a potential candidate for in vivo studies to determine its performance during operation. In this context, the hydrogel was injected into the divided vessel ends and perivascular area, allowing for direct suturing through the gel matrix. While our hydrogel effectively prevented vasospasm and facilitated micro- and supermicro-vascular anastomoses, it was noted that it did not cause significant changes in late-stage imaging and histopathological analysis up to 6 months. We strongly believe that pectin-based hydrogel potentially enhanced microlevel arterial, lymphatic, and particularly venous anastomoses.

Repurposing of FDA-approved drugs against active site and potential allosteric drug-binding sites of COVID-19 main protease.

The novel coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) still has serious negative effects on health, social life, and economics. Recently, vaccines from various companies have been urgently approved to control SARS-CoV-2 infections. However, any specific antiviral drug has not been confirmed so far for regular treatment. An important target is the main protease (Mpro ), which plays a major role in replication of the virus. In this study, Gaussian and residue network models are employed to reveal two distinct potential allosteric sites on Mpro that can be evaluated as drug targets besides the active site. Then, Food and Drug Administration (FDA)-approved drugs are docked to three distinct sites with flexible docking using AutoDock Vina to identify potential drug candidates. Fourteen best molecule hits for the active site of Mpro are determined. Six of these also exhibit high docking scores for the potential allosteric regions. Full-atom molecular dynamics simulations with MM-GBSA method indicate that compounds docked to active and potential allosteric sites form stable interactions with high binding free energy (∆Gbind ) values. ∆Gbind values reach -52.06 kcal/mol for the active site, -51.08 kcal/mol for the potential allosteric site 1, and - 42.93 kcal/mol for the potential allosteric site 2. Energy decomposition calculations per residue elucidate key binding residues stabilizing the ligands that can further serve to design pharmacophores. This systematic and efficient computational analysis successfully determines ivermectine, diosmin, and selinexor currently subjected to clinical trials, and further proposes bromocriptine, elbasvir as Mpro inhibitor candidates to be evaluated against SARS-CoV-2 infections.

Local and Global Motions Underlying Antibiotic Binding in Bacterial Ribosome.

The bacterial ribosome is one of the most important targets in the treatment of infectious diseases. As antibiotic resistance in bacteria poses a growing threat, a significant amount of effort is concentrated on exploring new drug-binding sites where testable predictions are of significance. Here, we study the dynamics of a ribosomal complex and 67 small and large subunits of the ribosomal crystal structures (64 antibiotic-bound, 3 antibiotic-free) from Deinococcus radiodurans, Escherichia coli, Haloarcula marismortui, and Thermus thermophilus by the Gaussian network model. Interestingly, a network of nucleotides coupled in high-frequency fluctuations reveals known antibiotic-binding sites. These sites are seen to locate at the interface of dynamic domains that have an intrinsic dynamic capacity to interfere with functional globular motions. The nucleotides and the residues fluctuating in the fast and slow modes of motion thus have promise for plausible antibiotic-binding and allosteric sites that can alter antibiotic binding and resistance. Overall, the present analysis brings a new dynamic perspective to the long-discussed link between small-molecule binding and large conformational changes of the supramolecule.

Noncovalent Pyrene-Polyethylene Glycol Coatings of Carbon Nanotubes Achieve in Vitro Biocompatibility.

Single-walled carbon nanotubes (SWNTs) have become increasingly exploited in biological applications, such as imaging and drug delivery. The application of SWNTs in biological settings requires the surface chemistry to remain through the low solubility in aqueous media. In this research, a facile approach for the preparation of a polyethylene glycol (PEG)-coated SWNT-based nanocarrier was reported. We focused on the effect of PEG chain length and SWNT size on the cytotoxicity of PEG-coated SWNTs as a superior drug delivery nanovector. First, all-atom molecular dynamics (MD) simulations were employed to explore the stability and behavior of SWNT/pyrene-PEG (SWNT/Pyr-PEG) structures at a molecular level that is not attainable with experiments. The MD studies revealed that (i) π-π stacking interactions between the pyrene bearing PEG molecules and SWNTs are maintained in bulky situations, regardless of PEG molecular weight or SWNT size; (ii) pyrene molecules diffuse over the SWNT surface without detaching; and (iii) both short and long dynamic Pyr-PEG chains have the capability of effectively coating the SWNT surface. In light of the simulations, noncovalent (π-π stacking) assemblies of SWNT/Pyr-PEG with different molecular weights of PEG ( Mw = 2000, 5000, and 12000) were successfully fabricated and characterized. For longer PEG chains, more effective coating of SWNTs was obtained, resulting in more biocompatible SWNT/Pyr-PEG nanomaterials. The number of SWNTs coated by Pyr-PEG was highly dependent on the length of pyrene bearing PEG polymers. Moreover, the short SWNTs showed a higher amount of PEG coating with respect to the long SWNTs. Cell viability results demonstrated a dose-dependent cytotoxicity of coated SWNTs. Short SWNTs coated with longer PEG chains have low cytotoxicity to be used in in vivo studies.

Identification of potential allosteric communication pathways between functional sites of the bacterial ribosome by graph and elastic network models.

BACKGROUND: Accumulated evidence indicates that bacterial ribosome employs allostery throughout its structure for protein synthesis. The nature of the allosteric communication between remote functional sites remains unclear, but the contact topology and dynamics of residues may play role in transmission of a perturbation to distant sites. METHODS/RESULTS: We employ two computationally efficient approaches - graph and elastic network modeling to gain insights about the allosteric communication in ribosome. Using graph representation of the structure, we perform k-shortest pathways analysis between peptidyl transferase center-ribosomal tunnel, decoding center-peptidyl transferase center - previously reported functional sites having allosteric communication. Detailed analysis on intact structures points to common and alternative shortest pathways preferred by different states of translation. All shortest pathways capture drug target sites and allosterically important regions. Elastic network model further reveals that residues along all pathways have the ability of quickly establishing pair-wise communication and to help the propagation of a perturbation in long-ranges during functional motions of the complex. CONCLUSIONS: Contact topology and inherent dynamics of ribosome configure potential communication pathways between functional sites in different translation states. Inter-subunit bridges B2a, B3 and P-tRNA come forward for their high potential in assisting allostery during translation. Especially B3 emerges as a potential druggable site. GENERAL SIGNIFICANCE: This study indicates that the ribosome topology forms a basis for allosteric communication, which can be disrupted by novel drugs to kill drug-resistant bacteria. Our computationally efficient approach not only overlaps with experimental evidence on allosteric regulation in ribosome but also proposes new druggable sites.

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation.

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of ß-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Exploring species-specific inhibitors with multiple target sites on S. aureus pyruvate kinase using a computational workflow.

Experimental evidence indicated that bacterial pyruvate kinase of glycolysis can be evaluated as an alternative target to eliminate infections, while antibiotic resistance poses a global threat. Here, we use a computational workflow to reveal and investigate the potential allosteric sites of methicillin-resistant S. aureus PK, which can help in designing species-specific drugs to inhibit activity of this organism. Residue interaction networks point to a known allosteric site at the small C-C interface, a potential allosteric site near the small interface (site #1), and a second potential allosteric site at the large interface (site #2). 2 µs-long molecular dynamics (MD) simulations with AMBER16 generate different conformations of one narrow target site. Known and potential allosteric sites on the selected conformers are investigated using ensemble docking with AutoDock Vina and a library of 2447 FDA-approved drugs. We determine 18 hits, comprising ergot-alkaloids, anti-cancer-agents, antivirals, analgesics, cardiac glycosides, all with a high docking z-score for three sites. 5 selected compounds with high, average and low z-scores are subjected to 50 ns-long MD simulations for MM-GBSA calculations. ΔGbind values up to -49.3 kcal/mol at the C-C interface, up to -32.7 kcal/mol at site #1, and up to -53.3 kcal/mol at site #2 support the docking calculations. We investigate mitapivat and TT-232 as reference compounds under clinical trial, targeting human PK isomers. We suggest 18 FDA-approved hits from the docking calculations and TT-232 as potential inhibitors with multiple target sites on S. aureus PK. This study also proposes pharmacophores models for de novo drug design.Communicated by Ramaswamy H. Sarma.

Potential allosteric sites captured in glycolytic enzymes via residue-based network models: Phosphofructokinase, glyceraldehyde-3-phosphate dehydrogenase and pyruvate kinase.

Likelihood of new allosteric sites for glycolytic enzymes, phosphofructokinase (PFK), glyceraldehyde-3-phosphate dehydrogenase (GADPH) and pyruvate kinase (PK) was evaluated for bacterial, parasitic and human species. Allosteric effect of a ligand binding at a site was revealed on the basis of low-frequency normal modes via Cα-harmonic residue network model. In bacterial PFK, perturbation of the proposed allosteric site outperformed the known allosteric one, producing a high amount of stabilization or reduced dynamics, on all catalytic regions. Another proposed allosteric spot at the dimer interface in parasitic PFK exhibited major stabilization effect on catalytic regions. In parasitic GADPH, the most desired allosteric response was observed upon perturbation of its tunnel region which incorporated key residues for functional regulation. Proposed allosteric site in bacterial PK produced a satisfactory allosteric response on all catalytic regions, whereas in human and parasitic PKs, a partial inhibition was observed. Residue network model based solely on contact topology identified the 'hub residues' with high betweenness tracing plausible allosteric communication pathways between distant functional sites. For both bacterial PFK and PK, proposed sites accommodated hub residues twice as much as the known allosteric site. Tunnel region in parasitic GADPH with the strongest allosteric effect among species, incorporated the highest number of hub residues. These results clearly suggest a one-to-one correspondence between the degree of allosteric effect and the number of hub residues in that perturbation site, which increases the likelihood of its allosteric nature.

Pectin-Zeolite-Based Wound Dressings with Controlled Albumin Release.

Hypoalbuminemia can lead to poor and delayed wound healing, while it is also associated with acute myocardial infarction, heart failure, malignancies, and COVID-19. In elective surgery, patients with low albumin have high risks of postoperative wound complications. Here, we propose a novel cost-effective wound dressing material based on low-methoxy pectin and NaA-zeolite particles with controlled albumin release properties. We focused on both albumin adsorption and release phenomena for wounds with excess exudate. Firstly, we investigated albumin dynamics and calculated electrostatic surfaces at experimental pH values in water by using molecular dynamics methods. Then, we studied in detail pectin-zeolite hydrogels with both adsorption and diffusion into membrane methods using different pH values and albumin concentrations. To understand if uploaded albumin molecules preserved their secondary conformation in different formulations, we monitored the effect of pH and albumin concentration on the conformational changes in albumin after it was released from the hydrogels by using CD-UV spectroscopy analyses. Our results indicate that at pH 6.4, BSA-containing films preserved the protein's folded structure while the protein was being released to the external buffer solutions. In vitro wound healing assay indicated that albumin-loaded hydrogels showed no toxic effects on the fibroblast cells.

Computational assessment of thermostability in miRNA:CNT system using molecular dynamics simulations.

BACKGROUND: Carbon nanotubes (CNTs) show great promise as theranostic agents due to their drug delivery properties, intrinsic near-infrared radiation-responsiveness, and magnetic functionalization. However, temperature elevation caused by these external stimuli during drug delivery should be considered for the evaluation of CNT-based systems loaded with temperature-sensitive biomolecules. METHODS: We examine the thermal stability of a 33 nucleotides long hairpin miRNA encapsulated in (20,20) CNT using all-atom molecular dynamics simulations in explicit water. We systematically increase the temperature as 298, 310, 327, and 343 K, reaching the melting temperature of miRNA. To emphasize the effect of the aromatic confined space, we compare the dynamics of miRNA inside the CNT to its dynamics free in the solution at the same temperatures, reaching a total simulation time of 7.9 µs. RESULTS: miRNA hairpin mostly maintains its double-stranded structure in the confined CNT, even at elevated temperatures. Binding free energies and potential of mean force calculations also underline the strong π-π interactions between the biomolecule and the CNT for 298-343 K. CONCLUSION: The let-7 miRNA mimic, which represents a wide family of RNAi-based therapeutics, can be transported in the CNT under medically applied hyperthermic conditions. GENERAL SIGNIFICANCE: This study shows how the structure and dynamics of miRNA hairpin are affected when encapsulated in an aromatic tube, during a systematic increase of temperature. It also indicates the high potential of CNT-based systems for the delivery of oligonucleotide therapeutics while simultaneous imaging/magnetic field guiding to the target tissue is achieved.

Fmoc-PEG Coated Single-Wall Carbon Nanotube Carriers by Non-covalent Functionalization: An Experimental and Molecular Dynamics Study.

Due to their structural characteristics at the nanoscale level, single-walled carbon nanotubes (SWNTs), hold great promise for applications in biomedicine such as drug delivery systems. Herein, a novel single-walled carbon nanotube (SWNT)-based drug delivery system was developed by conjugation of various Fmoc-amino acid bearing polyethylene glycol (PEG) chains (Mw = 2,000, 5,000, and 12,000). In the first step, full-atom molecular dynamics simulations (MD) were performed to identify the most suitable Fmoc-amino acid for an effective surface coating of SWNT. Fmoc-glycine, Fmoc-tryptophan, and Fmoc-cysteine were selected to attach to the PEG polymer. Here, Fmoc-cysteine and -tryptophan had better average interaction energies with SWNT with a high number of aromatic groups, while Fmoc-glycine provided a non-aromatic control. In the experimental studies, non-covalent modification of SWNTs was achieved by Fmoc-amino acid-bearing PEG chains. The remarkably high amount of Fmoc-glycine-PEG, Fmoc-tryptophan-PEG, and Fmoc-cysteine-PEG complexes adsorbed onto the SWNT surface, as was assessed via thermogravimetric and UV-vis spectroscopy analyses. Furthermore, Fmoc-cysteine-PEG5000 and Fmoc-cysteine-PEG12000 complexes displayed longer suspension time in deionized water, up to 1 and 5 week, respectively, underlying the ability of these surfactants to effectively disperse SWNTs in an aqueous environment. In vitro cell viability assays on human dermal fibroblast cells also showed the low cytotoxicity of these two samples, even at high concentrations. In conclusion, synthesized nanocarriers have a great potential for drug delivery systems, with high loading capacity, and excellent complex stability in water critical for biocompatibility.

Mechanism of cohesin loading onto chromosomes: a conformational dynamics study.

The structure-function relationship of cohesin, an essential chromosome maintenance protein, is investigated by analyzing its collective dynamics and conformational flexibility, enhancing our understanding of the sister chromatid cohesion process. A three-dimensional model of cohesin has been constructed by homology modeling using both crystallographic and electron microscopy image data. The harmonic dynamics of the cohesin structure are calculated with a coarse-grained elastic network model. The model shows that the bending motion of the cohesin ring is able to adopt a head-to-tail conformation, in agreement with experimental data. Low-frequency conformational changes are observed to deform the highly conserved glycine residues at the interface of the cohesin heterodimer. Normal mode analysis further reveals that, near large globular structures such as nucleosome and accessory proteins docked to cohesin, the mobility of the coiled-coil regions is notably affected. Moreover, fully solvated molecular dynamics calculations, performed specifically on the hinge region, indicate that hinge opening starts from one side of the dimerization interface, and is coordinated by highly conserved glycine residues.

Exploring Allosteric Signaling in the Exit Tunnel of the Bacterial Ribosome by Molecular Dynamics Simulations and Residue Network Model.

The bacterial ribosomal tunnel is equipped with numerous sites highly sensitive to the course of the translation process. This study investigates allosteric pathways linking distant functional sites that collaboratively play a role either in translation regulation or recruitment of chaperones. We apply perturbation response scanning (PRS) analysis to 700 ns long and 500 ns long coarse-grained molecular dynamics simulations of E. coli and T. thermophilus large subunits, respectively, to reveal nucleotides/residues with the ability to transmit perturbations by dynamic rationale. We also use the residue network model with the k-shortest pathways method to calculate suboptimal pathways based on the contact topology of the ribosomal tunnel of E. coli crystal structure and 101 ClustENM generated conformers of T. thermophilus large subunit. In the upper part of the tunnel, results suggest that A2062 and A2451 can communicate in both directions for translation stalling, mostly through dynamically coupled C2063, C2064, and A2450. For a similar purpose, U2585 and U2586 are coupled with A2062, while they are also sensitive to uL4 and uL22 at the constriction region through two different pathways at the opposite sides of the tunnel wall. In addition, the constriction region communicates with the chaperone binding site on uL23 at the solvent side but through few nucleotides. Potential allosteric communication pathways between the lower part of the tunnel and chaperone binding site mostly use the flexible loop of uL23, while A1336-G1339 provide a suboptimal pathway. Both species seem to employ similar mechanisms in the long tunnel, where a non-conserved cavity at the bacterial uL23 and 23S rRNA interface is proposed as a novel drug target.

Focused functional dynamics of supramolecules by use of a mixed-resolution elastic network model.

The mixed-resolution elastic network model was introduced previously for computing the motions of a structure, which is described at different levels of detail in different parts, for example, with atomistic and residue-level regions. This method has proved to be an efficient tool to explore the collective dynamics of proteins with some atomistic details, which would be difficult to obtain with either conventional full-atom approaches or fully coarse-grained models. Understanding function often requires atomic detail, but not necessarily for the entire structure. In this study, the calculation of the interaction forces between different resolution regions for the hierarchical levels of coarse-graining is further elaborated on in the new approach by considering explicitly the atomic contacts in the crystal structure. The collective dynamics of the enzyme triosephosphate isomerase and its active site together with loop 6 motions are considered in detail. The supramolecular assemblage ribosome and local atomic motions in its "interesting" functional part-the decoding center-are investigated for the low frequency range of the spectrum with high computational efficiency. This new atom-based mixed coarse-graining approach can be effectively used to generate realistic high-resolution conformations of extremely large protein-DNA or RNA complexes by performing energy minimization on structures deformed along the normal modes of the elastic network model. The new model permits focusing on specific functional parts that move in coordination and response to the remainder of the entire structure.

Collective dynamics of the ribosomal tunnel revealed by elastic network modeling.

The collective dynamics of the nascent polypeptide exit tunnel are investigated with the computationally efficient elastic network model using normal mode analysis. The calculated normal modes are considered individually and in linear combinations with different coefficients mimicking the phase angles between modes, in order to follow the mechanistic motions of tunnel wall residues. The low frequency fluctuations indicate three distinct regions along the tunnel-the entrance, the neck, and the exit-each having distinctly different domain motions. Generally, the lining of the entrance region moves in the exit direction, with the exit region having significantly larger motions, but in a perpendicular direction, whereas the confined neck region has rotational motions. Especially the universally conserved extensions of ribosomal proteins L4 and L22 located at the narrowest and mechanistically strategic region of tunnel undergo generally anti- or non-correlated motions, which may have an important role in nascent polypeptide gating mechanism. These motions appear to be sufficiently robust so as to be unaffected by the presence of a peptide chain in the tunnel.

The ribosome structure controls and directs mRNA entry, translocation and exit dynamics.

The protein-synthesizing ribosome undergoes large motions to effect the translocation of tRNAs and mRNA; here, the domain motions of this system are explored with a coarse-grained elastic network model using normal mode analysis. Crystal structures are used to construct various model systems of the 70S complex with/without tRNA, elongation factor Tu and the ribosomal proteins. Computed motions reveal the well-known ratchet-like rotational motion of the large subunits, as well as the head rotation of the small subunit and the high flexibility of the L1 and L7/L12 stalks, even in the absence of ribosomal proteins. This result indicates that these experimentally observed motions during translocation are inherently controlled by the ribosomal shape and only partially dependent upon GTP hydrolysis. Normal mode analysis further reveals the mobility of A- and P-tRNAs to increase in the absence of the E-tRNA. In addition, the dynamics of the E-tRNA is affected by the absence of the ribosomal protein L1. The mRNA in the entrance tunnel interacts directly with helicase proteins S3 and S4, which constrain the mRNA in a clamp-like fashion, as well as with protein S5, which likely orients the mRNA to ensure correct translation. The ribosomal proteins S7, S11 and S18 may also be involved in assuring translation fidelity by constraining the mRNA at the exit site of the channel. The mRNA also interacts with the 16S 3' end forming the Shine-Dalgarno complex at the initiation step; the 3' end may act as a 'hook' to reel in the mRNA to facilitate its exit.

Exploring allosteric communication in multiple states of the bacterial ribosome using residue network analysis.

Antibiotic resistance is one of the most important problems of our era and hence the discovery of new effective therapeutics is urgent. At this point, studying the allosteric communication pathways in the bacterial ribosome and revealing allosteric sites/residues is critical for designing new inhibitors or repurposing readily approved drugs for this enormous machine. To shed light onto molecular details of the allosteric mechanisms, here we construct residue networks of the bacterial ribosomal complex at four different states of translation by using an effective description of the intermolecular interactions. Centrality analysis of these networks highlights the functional roles of structural components and critical residues on the ribosomal complex. High betweenness scores reveal pathways of residues connecting numerous sites on the structure. Interestingly, these pathways assemble highly conserved residues, drug binding sites, and known allosterically linked regions on the same structure. This study proposes a new residue-level model to test how distant sites on the molecular machine may be linked through hub residues that are critically located on the contact topology to inherently form communication pathways. Findings also indicate intersubunit bridges B1b, B3, B5, B7, and B8 as critical targets to design novel antibiotics.

Conformational dynamics of bacterial trigger factor in apo and ribosome-bound states.

The chaperone trigger factor (TF) binds to the ribosome exit tunnel and helps cotranslational folding of nascent chains (NC) in bacterial cells and chloroplasts. In this study, we aim to investigate the functional dynamics of fully-atomistic apo TF and its complex with 50S. As TF accomodates a high percentage of charged residues on its surface, the effect of ionic strength on TF dynamics is assessed here by performing five independent molecular dynamics (MD) simulations (total of 1.3 micro-second duration) at 29 mM and 150 mM ionic strengths. At both concentrations, TF exhibits high inter- and intra-domain flexibility related to its binding (BD), core (CD), and head (HD) domains. Even though large oscillations in gyration radius exist during each run, we do not detect the so-called 'fully collapsed' state with both HD and BD collapsed upon the core. In fact, the extended conformers with relatively open HD and BD are highly populated at the physiological concentration of 150 mM. More importantly, extended TF snapshots stand out in terms of favorable docking onto the 50S subunit. Elastic network modeling (ENM) indicates significant changes in TF's functional dynamics and domain decomposition depending on its conformation and positioning on the 50S. The most dominant slow motions are the lateral sweeping and vertical opening/closing of HD relative to 50S. Finally, our ENM-based efficient technique -ClustENM- is used to sample atomistic conformers starting with an extended TF-50S complex. Specifically, BD and CD motions are restricted near the tunnel exit, while HD is highly mobile. The atomistic conformers generated without an NC are in agreement with the cryo-EM maps available for TF-ribosome-NC complex.

Collective dynamics of EcoRI-DNA complex by elastic network model and molecular dynamics simulations.

Anisotropic network model (ANM) is used to analyze the collective motions of restriction enzyme EcoRI in free form and in complex with DNA. For comparison, three independent molecular dynamics (MD) simulations, each of 1.5 ns duration, are also performed for the EcoRI-DNA complex in explicit water. Although high mobility (equilibrium fluctuations) of inner and outer loops that surround the DNA is consistent in both methods and experiments, MD runs sample different conformational subspaces from which reliable collective dynamics cannot be extracted. However, ANM employed on different conformations from MD simulations indicates very similar collective motions. The stems of the inner loops are quite immobile even in the free enzyme and form a large, almost fixed, pocket for DNA binding. As a result, the residues that make specific and non-specific interactions with the DNA exhibit very low fluctuations in the free enzyme. The vibrational entropy difference between the EcoRI complex and free protein + unkinked DNA is positive (favorable), which may partially counteract the unfavorable enthalpy difference of DNA kink formation. Dynamic domains in EcoRI complex and cross-correlations between residue fluctuations indicate possible means of communication between the distal active sites.