Structure fluctuations and conformational changes in protein binding

Efficient delivery of carotenoids to adipocytes with albumin

Carotenoids are very effectively delivered by albumin to adipocytes. The uptake of carotenoids to the cells occurs in the form of self-aggregates that localize in the vicinity of the adipocyte membrane, as shown by high spatial resolution Raman spectroscopy. The binding of carotenoids to albumin and the mechanism of their transport were elucidated with the help of chiroptical spectroscopies, in tandem with molecular docking and molecular dynamics simulations. In particular, apart from the recognized high affinity pocket of albumin that binds a carotenoid monomer in domain I, we have identified a hydrophobic periphery area in domain IIIB that loosely bounds the self-aggregated carotenoid in aqueous media and enables its easy detachment in hydrophobic environments. This explains the effectiveness of albumins as nanocarriers of carotenoids to adipocytes in vitro.

Comparative modeling and enzymatic affinity of novel haloacid dehalogenase from Bacillus megaterium strain BHS1 isolated from alkaline Blue Lake in Turkey

This study presents the initial structural model of L-haloacid dehalogenase (DehLBHS1) from Bacillus megaterium BHS1, an alkalotolerant bacterium known for its ability to degrade halogenated environmental pollutants. The model provides insights into the structural features of DehLBHS1 and expands our understanding of the enzymatic mechanisms involved in the degradation of these hazardous pollutants. Key amino acid residues (Arg40, Phe59, Asn118, Asn176, and Trp178) in DehLBHS1 were identified to play critical roles in catalysis and molecular recognition of haloalkanoic acid, essential for efficient binding and transformation of haloalkanoic acid molecules. DehLBHS1 was modeled using I-TASSER, yielding a best TM-score of 0.986 and an RMSD of 0.53 Å. Validation of the model using PROCHECK revealed that 89.2% of the residues were located in the most favored region, providing confidence in its structural accuracy. Molecular docking simulations showed that the non-simulated DehLBHS1 preferred 2,2DCP over other substrates, forming one hydrogen bond with Arg40 and exhibiting a minimum energy of -2.5 kJ/mol. The simulated DehLBHS1 exhibited a minimum energy of -4.3 kJ/mol and formed four hydrogen bonds with Arg40, Asn176, Asp9, and Tyr11, further confirming the preference for 2,2DCP. Molecular dynamics simulations supported this preference, based on various metrics, including RMSD, RMSF, gyration, hydrogen bonding, and molecular distance. MM-PBSA calculations showed that the DehLBHS1-2,2-DCP complex had a markedly lower binding energy (-21.363 ± 1.26 kcal/mol) than the DehLBHS1-3CP complex (-14.327 ± 1.738 kcal/mol). This finding has important implications for the substrate specificity and catalytic function of DehLBHS1, particularly in the bioremediation of 2,2-DCP in contaminated alkaline environments. These results provide a detailed view of the molecular interactions between the enzyme and its substrate and may aid in the development of more efficient biocatalytic strategies for the degradation of halogenated compounds.Communicated by Ramaswamy H. Sarma.

Characterization of Climbazole-Bovine serum albumin interaction by experimental and in silico approaches

Interactive association of an antifungal drug, climbazole (CBZ) with the carrier protein in bovine circulation, bovine serum albumin (BSA) was explored by fluorescence and absorption spectroscopy along with in silico techniques. The fluorescence and absorption spectral alterations of the protein upon addition of CBZ affirmed the complex foration between CBZ and BSA. The inverse temperature dependence behaviour of the K_SV values as well as the hyperchromic result of the protein's absorption signals characterized CBZ-triggered quenching of BSA fluorescence as the static quenching. A weak binding affinity (K_a = 3.12-1.90-× 10³ M^-1) was reported towards the CBZ-BSA association process. Interpretation of thermodynamic data (entropy change = +14.68 J mol^-1 K^-1 and enthalpy change = -15.07 kJ mol^-1) and in silico analyses anticipated that hydrophobic forces, van der Waals forces and hydrogen bonds were the key intermolecular forces in the complex stabilization. Inclusion of CBZ to BSA produced microenvironmental perturbations around Tyr and Trp residues, and also significantly defended temperature-induced destabilization of BSA. The binding locus of CBZ was detected in the proximity of Sudlow's sites I (subdomain IIA) and II (subdomain IIIA) of BSA, exhibiting greater preference towards site II, as revealed by competitive site-marker displacement investigations and in silico analysis. The stability of the CBZ-BSA complex was further validated by the molecular dynamics simulation assessments.

Study of Allosteric Transitions of Human P-Glycoprotein by Using the Two-State Anisotropic Network Model

Human P-glycoprotein (P-gp) is a kind of ATP-binding cassette (ABC) transporters. Once human P-gp is overexpressed in tumor cells, which can lead to tumor multidrug resistance (MDR). However, the present experimental methods are difficult to obtain the large-scale conformational transition process of human P-gp. In this work, we explored the allosteric pathway of human P-gp from the inward-facing (IF) to the outward-facing (OF) state in the substrate transport process with the two-state anisotropic network model (tANM). These results suggest that the allosteric transitions proceed in a coupled way. The conformational changes of nucleotide-binding domains (NBDs) finally make the transmembrane domains (TMDs) to the OF state via the role of the allosteric propagation of the intracellular helices IH1 and IH2. Additionally, this allosteric pathway is advantageous in energy compared with other methods. This study reveals the conformational transition of P-gp, which contributes to an understanding of the allosteric mechanism of ABC exporters.

Structural modelling of SARS-CoV-2 alpha variant (B.1.1.7) suggests enhanced furin binding and infectivity

The B.1.1.7 SARS-CoV-2 strain that has emerged in the UK in early December presents seven mutations and three deletions on S-protein structure that could lead to a more infective strain. The P681H mutation in the “PRRAR” furin cleavage site might affect the binding affinity to furin enzyme and hence its infectivity.

Therefore, in this study, various structural bioinformatics approaches were used to model the S-protein structure with the B.1.1.7 variant amino acid substitutions and deletions. In addition to modelling the binding of furin to the cleavage site of the wild-type and the B.1.1.7 variant.

Conclusively the B.1.1.7 variant resulted in dynamic stability, conformational changes and variations in binding energies in the S-protein structure, resulting in a more favourable binding of furin enzyme to the SARS-CoV-2 S-protein.

An experimental approach probing the conformational transitions and energy landscape of antibodies: a glimmer of hope for reviving lost therapeutic candidates using ionic liquid

Understanding protein folding in different environmental conditions is fundamentally important for predicting protein structures and developing innovative antibody formulations. While the thermodynamics and kinetics of folding and unfolding have been extensively studied by computational methods, experimental methods for determining antibody conformational transition pathways are lacking. Motivated to fill this gap, we prepared a series of unique formulations containing a high concentration of a chimeric immunoglobin G4 (IgG4) antibody with different excipients in the presence and absence of the ionic liquid (IL) choline dihydrogen phosphate. We determined the effects of different excipients and IL on protein thermal and structural stability by performing variable temperature circular dichroism and bio-layer interferometry analyses. To further rationalise the observations of conformational changes with temperature, we carried out molecular dynamics simulations on a single antibody binding fragment from IgG4 in the different formulations, at low and high temperatures. We developed a methodology to study the conformational transitions and associated thermodynamics of biomolecules, and we showed IL-induced conformational transitions. We showed that the increased propensity for conformational change was driven by preferential binding of the dihydrogen phosphate anion to the antibody fragment. Finally, we found that a formulation containing IL with sugar, amino acids and surfactant is a promising candidate for stabilising proteins against conformational destabilisation and aggregation. We hope that ultimately, we can help in the quest to understand the molecular basis of the stability of antibodies and protein misfolding phenomena and offer new candidate formulations with the potential to revive lost therapeutic candidates.

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.

Glycosylation effects on the structure and dynamics of a full-length Cel7A cellulase

Fungi cellulases are used to degrade cellulose-containing biomass for bioethanol production. Industrial cellulases such as Cel7A from Trichoderma reesei (TrCel7A) are critical in this process. Thus, the understanding of structure and dynamics is crucial for engineering variants with improved cellulolytic activity. This cellulase consists of two domains connected by a flexible and highly glycosylated linker. However, the linker flexibility has hindered the determination of Cel7A complete structure. Herein, based on atomic and sparse data, we applied integrative modelling to build a model of the complete enzyme structure. Next, through simulations, we studied the glycosylation effects on the structure and dynamics of a solubilized TrCel7A. Essential dynamics analysis showed that O-glycosylation in the linker led to the stabilization of protein overall dynamics. O-linked glycans seem to restrict protein dihedral angles distribution in this region, selecting more elongated conformations. Besides the reduced flexibility, functional interdomain motions occurred in a more concerted way in the glycosylated system. In contrast, in the absence of glycosylation, we observed vast conformational plasticity with the functional domains frequently collapsing. We report here evidence that targeting Cel7A linker flexibility by point mutations including modification of glycosylation sites could be a promising design strategy to improve cellulase activity.

B-factor Guided Proline Substitutions in Chromobacterium violaceum Amine Transaminase: Evaluation of the Proline Rule as a Method for Enzyme Stabilization

Biocatalysis is attracting interest in the chemical industry as a sustainable alternative in large-scale chemical transformations. However, low operational stability of naturally evolved enzymes is a challenge and major efforts are required to engineer protein stability, usually by directed evolution. The development of methods for protein stabilization based on rational design is of great interest, as it would minimize the efforts needed to generate stable enzymes. Here we present a rational design strategy based on proline substitutions in flexible areas of the protein identified by analyzing B-factors. Several proline substitutions in the amine transaminase from Chromobacterium violaceum were shown to have a positive impact on stability with increased half-life at 60 °C by a factor of 2.7 (variant K69P/D218P/K304P/R432P) as well as increased melting temperature by 8.3 °C (variant K167P). Finally, the presented method utilizing B-factor analysis in combination with the proline rule was deemed successful at increasing the stability of this enzyme.

Substrate binding versus escape dynamics in a pH-affected fungal beta-glucosidase revealed by molecular dynamics simulations

The cellulolytic ability of fungal species is important to both natural and engineered biocycling of plant matter. One essential step is the conversion of cellobiose into glucose catalyzed by beta-glucosidases. Mutagenesis studies have implicated altering the substrate binding pocket to influence the pH-activity profile of this enzyme. However, structural understanding of the pH-affected substrate binding environment is lacking. Here we conducted molecular dynamics simulations of fully hydrated TrBgl2, a beta-glucosidase of Trichoderma reesei, equilibrated at its optimal pH (pH 6) and two unfavorable pHs (pH 5 and pH 7.5). We identified structural arrangement of specific residues that facilitated substrate escape from the catalytic site at pH 5 but locked the bound substrate in an unfavorable orientation at pH 7.5. For comparative analysis, we also performed simulations of a mutated TrBgl2 with previously demonstrated improved catalysis as a function of pH. We captured the responsible conformational changes in the engineered substrate binding pocket.

Modeling of Protein Structural Flexibility and Large-Scale Dynamics: Coarse-Grained Simulations and Elastic Network Models

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Insulin adsorption onto zinc oxide nanoparticle mediates conformational rearrangement into amyloid-prone structure with enhanced cytotoxic propensity

BACKGROUND

Injection localized amyloidosis is one of the most prevalent disorders in type II diabetes mellitus (TIIDM) patients relying on insulin injections. Previous studies have reported that nanoparticles can play a role in the amyloidogenic process of proteins. Hence, the present study deals with the effect of zinc oxide nanoparticles (ZnONP) on the amyloidogenicity and cytotoxicity of insulin.

METHODS

ZnONP is synthesised and characterized using XRD, Zeta Sizer, UV-Visible spectroscope and TEM. The characterization is followed by ZnONP interaction with insulin, which is studied employing fluorescence spectroscopes, isothermal titration calorimetry and molecular dynamics simulations. The interaction leads insulin conformational rearrangement into amyloid-like fibril, which is studied using thioflavin T dye binding assay, circular dichroism spectroscopy and TEM, followed by cytotoxicity propensity using Alamar Blue dye reduction assay.

RESULTS

Insulin has very weak interaction with ZnONP interface. Insulin at studied concentration forms amorphous aggregates at physiological pH, whereas in presence of ZnONP interface amyloid-like fibrils are formed. While the amyloid-like fibrils are cytotoxic to MIN6 and THP-1 cell lines, insulin and ZnONP individual solutions and their fresh mixtures enhance the cells proliferation.

CONCLUSIONS

The presence of ZnONP interface enhances insulin fibrillation at physiological pH by providing a favourable template for the nucleation and growth of insulin amyloids.

GENERAL SIGNIFICANCE

The studied protein-nanoparticle system from protein conformational dynamics point of view throws caution over nanoparticle use in biological applications, especially in vivo applications, considering the amyloidosis a very slow but non-curable degenerative disease.

Metal-catalyzed oxidation of human serum albumin does not alter the interactive binding to the two principal drug binding sites

It is well known that various physiological factors such as pH, endogenous substances or post-translational modifications can affect the conformational state of human serum albumin (HSA). In a previous study, we reported that both pH- and long chain fatty acid-induced conformational changes can alter the interactive binding of ligands to the two principal binding sites of HSA, namely, site I and site II. In the present study, the effect of metal-catalyzed oxidation (MCO) caused by ascorbate/oxygen/trace metals on HSA structure and the interactive binding between dansyl-L-asparagine (DNSA; a site I ligand) and ibuprofen (a site II ligand) at pH 6.5 was investigated. MCO was accompanied by a time-dependent increase in carbonyl content in HSA, suggesting that the HSA was being oxidized. In addition, The MCO of HSA was accompanied by a change in net charge to a more negative charge and a decrease in thermal stability. SDS-PAGE patterns and α-helical contents of the oxidized HSAs were similar to those of native HSA, indicating that the HSA had not been extensively structurally modified by MCO. MCO also caused a selective decrease in ibuprofen binding. In spite of the changes in the HSA structure and ligand that bind to site II, no change in the interactive binding between DNSA and ibuprofen was observed. These data indicated that amino acid residues in site II are preferentially oxidized by MCO, whereas the spatial relationship between sites I and II (e.g. the distance between sites), the flexibility or space of each binding site are not altered. The present findings provide insights into the structural characteristics of oxidized HSA, and drug binding and drug-drug interactions on oxidized HSA.

Specific stereochemistry of OP-1074 disrupts estrogen receptor alpha helix 12 and confers pure antiestrogenic activity

Complex tissue-specific and cell-specific signaling by the estrogen receptor (ER) frequently leads to the development of resistance to endocrine therapy for breast cancer. Pure ER antagonists, which completely lack tissue-specific agonist activity, hold promise for preventing and treating endocrine resistance, however an absence of structural information hinders the development of novel candidates. Here we synthesize a small panel of benzopyrans with variable side chains to identify pure antiestrogens in a uterotrophic assay. We identify OP-1074 as a pure antiestrogen and a selective ER degrader (PA-SERD) that is efficacious in shrinking tumors in a tamoxifen-resistant xenograft model. Biochemical and crystal structure analyses reveal a structure activity relationship implicating the importance of a stereospecific methyl on the pyrrolidine side chain of OP-1074, particularly on helix 12.

Estrogen receptor alpha (ERα) plays critical roles in the etiology and treatment of breast cancer. Here the authors synthesize benzopyrans with variable side chains to identify antiestrogenic compounds and characterize OP-1074, a compound that exhibits pure antiestrogenic activity by inducing the degradation of ERα and possesses greater potency than tamoxifen or fulvestrant in a xenograft model.

Alanine mutation of the catalytic sites of Pantothenate Synthetase causes distinct conformational changes in the ATP binding region

The enzyme Pantothenate synthetase (PS) represents a potential drug target in Mycobacterium tuberculosis. Its X-ray crystallographic structure has demonstrated the significance and importance of conserved active site residues including His44, His47, Asn69, Gln72, Lys160 and Gln164 in substrate binding and formation of pantoyl adenylate intermediate. In the current study, molecular mechanism of decreased affinity of the enzyme for ATP caused by alanine mutations was investigated using molecular dynamics (MD) simulations and free energy calculations. A total of seven systems including wild-type + ATP, H44A + ATP, H47A + ATP, N69A + ATP, Q72A + ATP, K160A + ATP and Q164A + ATP were subjected to 50 ns MD simulations. Docking score, MM-GBSA and interaction profile analysis showed weak interactions between ATP (substrate) and PS (enzyme) in H47A and H160A mutants as compared to wild-type, leading to reduced protein catalytic activity. However, principal component analysis (PCA) and free energy landscape (FEL) analysis revealed that ATP was strongly bound to the catalytic core of the wild-type, limiting its movement to form a stable complex as compared to mutants. The study will give insight about ATP binding to the PS at the atomic level and will facilitate in designing of non-reactive analogue of pantoyl adenylate which will act as a specific inhibitor for PS.

Searching for potential mTOR inhibitors: Ligand-based drug design, docking and molecular dynamics studies of rapamycin binding site

The PI3K/Akt/mTOR pathway is an important intracellular signaling pathway in cell cycle regulation and its dysregulation is associated with various types of diseases. mTOR (mechanistic or mammalian target of rapamycin) is the main enzyme that performs intermediate control of the signaling pathway through a phosphotransfer process. The classical inhibition of the mTOR pathway is effected by rapamycin and its analogous blocking allosterically the catalytic phosphorylation site, avoiding the deleterious side effects induced by ATP-competitive inhibitors. We employed ligand-based drug design strategies such as pharmacophore searching and analysis, molecular docking, absorption, distribution, metabolism, excretion and toxicity (ADMETox) properties filtering, and molecular dynamics to select potential molecules to become non-ATP competitive inhibitors of the mTOR complex. According to our findings, we propose eight novel potential mTOR inhibitors with similar or better properties than the classic inhibitor complex, rapamycin.

Combining protein sequence, structure, and dynamics: A novel approach for functional evolution analysis of PAS domain superfamily

PAS domains are widespread in archaea, bacteria, and eukaryota, and play important roles in various functions. In this study, we aim to explore functional evolutionary relationship among proteins in the PAS domain superfamily in view of the sequence-structure-dynamics-function relationship. We collected protein sequences and crystal structure data from RCSB Protein Data Bank of the PAS domain superfamily belonging to three biological functions (nucleotide binding, photoreceptor activity, and transferase activity). Protein sequences were aligned and then used to select sequence-conserved residues and build phylogenetic tree. Three-dimensional structure alignment was also applied to obtain structure-conserved residues. The protein dynamics were analyzed using elastic network model (ENM) and validated by molecular dynamics (MD) simulation. The result showed that the proteins with same function could be grouped by sequence similarity, and proteins in different functional groups displayed statistically significant difference in their vibrational patterns. Interestingly, in all three functional groups, conserved amino acid residues identified by sequence and structure conservation analysis generally have a lower fluctuation than other residues. In addition, the fluctuation of conserved residues in each biological function group was strongly correlated with the corresponding biological function. This research suggested a direct connection in which the protein sequences were related to various functions through structural dynamics. This is a new attempt to delineate functional evolution of proteins using the integrated information of sequence, structure, and dynamics.

NOLB: Nonlinear Rigid Block Normal-Mode Analysis Method

We present a new conceptually simple and computationally efficient method for nonlinear normal-mode analysis called NOLB. It relies on the rotations-translations of blocks (RTB) theoretical basis developed by Y.-H. Sanejouand and colleagues [ Durand et al. Biopolymers 1994 , 34 , 759 - 771 . Tama et al. Proteins: Struct., Funct., Bioinf . 2000 , 41 , 1 - 7 ]. We demonstrate how to physically interpret the eigenvalues computed in the RTB basis in terms of angular and linear velocities applied to the rigid blocks and how to construct a nonlinear extrapolation of motion out of these velocities. The key observation of our method is that the angular velocity of a rigid block can be interpreted as the result of an implicit force, such that the motion of the rigid block can be considered as a pure rotation about a certain center. We demonstrate the motions produced with the NOLB method on three different molecular systems and show that some of the lowest frequency normal modes correspond to the biologically relevant motions. For example, NOLB detects the spiral sliding motion of the TALE protein, which is capable of rapid diffusion along its target DNA. Overall, our method produces better structures compared to the standard approach, especially at large deformation amplitudes, as we demonstrate by visual inspection, energy, and topology analyses and also by the MolProbity service validation. Finally, our method is scalable and can be applied to very large molecular systems, such as ribosomes. Standalone executables of the NOLB normal-mode analysis method are available at https://team.inria.fr/nano-d/software/nolb-normal-modes/ . A graphical user interface created for the SAMSON software platform will be made available at https://www.samson-connect.net .

Knowledge-based entropies improve the identification of native protein structures

Evaluating protein structures requires reliable free energies with good estimates of both potential energies and entropies. Although there are many demonstrated successes from using knowledge-based potential energies, computing entropies of proteins has lagged far behind. Here we take an entirely different approach and evaluate knowledge-based conformational entropies of proteins based on the observed frequencies of contact changes between amino acids in a set of 167 diverse proteins, each of which has two alternative structures. The results show that charged and polar interactions break more often than hydrophobic pairs. This pattern correlates strongly with the average solvent exposure of amino acids in globular proteins, as well as with polarity indices and the sizes of the amino acids. Knowledge-based entropies are derived by using the inverse Boltzmann relationship, in a manner analogous to the way that knowledge-based potentials have been extracted. Including these new knowledge-based entropies almost doubles the performance of knowledge-based potentials in selecting the native protein structures from decoy sets. Beyond the overall energy-entropy compensation, a similar compensation is seen for individual pairs of interacting amino acids. The entropies in this report have immediate applications for 3D structure prediction, protein model assessment, and protein engineering and design.

New potential inhibitors of mTOR: a computational investigation integrating molecular docking, virtual screening and molecular dynamics simulation

The mTOR (mammalian or mechanistic Target Of Rapamycin), a complex metabolic pathway that involves multiple steps and regulators, is a major human metabolic pathway responsible for cell growth control in response to multiple factors and that is dysregulated in various types of cancer. The classical inhibition of the mTOR pathway is performed by rapamycin and its analogs (rapalogs). Considering that rapamycin binds to an allosteric site and performs a crucial role in the inhibition of the mTOR complex without causing the deleterious side effects common to ATP-competitive inhibitors, we employ ligand-based drug design strategies, such as virtual screening methodology, computational determination of ADME/Tox properties of selected molecules, and molecular dynamics in order to select molecules with the potential to become non-ATP-competitive inhibitors of the mTOR enzymatic complex. Our findings suggest five novel potential mTOR inhibitors, with similar or better properties than the classic inhibitor complex, rapamycin.

PreFRP: Prediction and visualization of fluctuation residues in proteins

Aim:

The PreFRP web server extracts sequence and basic information of a protein structure and groups amino acid residues in a protein into three important types such as high, moderate, and weak fluctuating residues.

Materials and Methods:

The server takes a protein data bank file or an amino acid sequence as input and prints the probability of amino acid residues to fluctuate. The server also provides a link to Jmol, a molecular visualization program to visualize the high, moderate, and weak fluctuating residues in three different colors.

Results:

Prediction and visualization of fluctuating amino acid residues in proteins may help to understand the complex three-dimensional structure of proteins and may further help in docking and mutation experiments.

Availability:

The web server is freely accessible through the web page of the author's institution http://www.mpi.edu.in/prefrp/link.html.

Local packing modulates diversity of iron pathways and cooperative behavior in eukaryotic and prokaryotic ferritins

Ferritin-like molecules show a remarkable combination of the evolutionary conserved activity of iron uptake and release that engage different pores in the conserved ferritin shell. It was hypothesized that pore selection and iron traffic depend on dynamic allostery with no conformational changes in the backbone. In this study, we detect the allosteric networks in Pseudomonas aeruginosa bacterioferritin (BfrB), bacterial ferritin (FtnA), and bullfrog M and L ferritins (Ftns) by a network-weaving algorithm (NWA) that passes threads of an allosteric network through highly correlated residues using hierarchical clustering. The residue-residue correlations are calculated in the packing-on elastic network model that introduces atom packing into the common packing-off model. Applying NWA revealed that each of the molecules has an extended allosteric network mostly buried inside the ferritin shell. The structure of the networks is consistent with experimental observations of iron transport: The allosteric networks in BfrB and FtnA connect the ferroxidase center with the 4-fold pores and B-pores, leaving the 3-fold pores unengaged. In contrast, the allosteric network directly links the 3-fold pores with the 4-fold pores in M and L Ftns. The majority of the network residues are either on the inner surface or buried inside the subunit fold or at the subunit interfaces. We hypothesize that the ferritin structures evolved in a way to limit the influence of functionally unrelated events in the cytoplasm on the allosteric network to maintain stability of the translocation mechanisms. We showed that the residue-residue correlations and the resultant long-range cooperativity depend on the ferritin shell packing, which, in turn, depends on protein sequence composition. Switching from the packing-on to the packing-off model reduces correlations by 35%-38% so that no allosteric network can be found. The influence of the side-chain packing on the allosteric networks explains the diversity in mechanisms of iron traffic suggested by experimental approaches.

B-factor Analysis and Conformational Rearrangement of Aldose Reductase

The NADPH-dependent reduction of glucose reaction that is catalyzed by Aldose Reductase (AR) follows a sequential ordered kinetic mechanism in which the co-factor NADPH binds to the enzyme prior to the aldehyde substrate. The kinetic/structural experiments have found a conformational change involving a hinge-like movement of a surface loop (residues 213-224) which is anticipated to take place upon the binding of the diphosphate moiety of NADPH. The reorientation of this loop, expected to permit the release of NADP⁺, represents the rate-limiting step of the catalytic mechanism. This study reveals: 1) The Translation/Libration/Screw (TLS) analysis of absolute B-factors of apo AR crystal structures indicates that the 212-224 loop might move as a rigid group. 2) Residues that make the flexible loop slide in the AR binary and ternary complexes. 3) The normalized B-factors separate this segment into three different clusters with fewer residues.