Molecular simulations of metal-coupled protein folding

Structural studies of catalytic peptides using molecular dynamics simulations

Many self-assembling peptides can form amyloid like structures with different sizes and morphologies. Driven by non-covalent interactions, their aggregation can occur through distinct pathways. Additionally, they can bind metal ions to create enzyme like active sites that allow them to catalyze diverse reactions. Due to the non-crystalline nature of amyloids, it is quite challenging to elucidate their structures using experimental spectroscopic techniques. In this aspect, molecular dynamics (MD) simulations provide a useful tool to derive structures of these macromolecules in solution. They can be further validated by comparing with experimentally measured structural parameters. However, these simulations require a multi-step process starting from the selection of the initial structure to the analysis of MD trajectories. There are multiple force fields, parametrization protocols, equilibration processes, software and analysis tools available for this process. Therefore, it is complicated for non-experts to select the most relevant tools and perform these simulations effectively. In this chapter, a systematic methodology that covers all major aspects of modeling of catalytic peptides is provided in a user-friendly manner. It will be helpful for researchers in this critical area of research.

The nucleotide excision repair proteins through the lens of molecular dynamics simulations

Mutations that affect the proteins responsible for the nucleotide excision repair (NER) pathway can lead to diseases such as xeroderma pigmentosum, trichothiodystrophy, Cockayne syndrome, and Cerebro-oculo-facio-skeletal syndrome. Hence, understanding their molecular behavior is needed to elucidate these diseases' phenotypes and how the NER pathway is organized and coordinated. Molecular dynamics techniques enable the study of different protein conformations, adaptable to any research question, shedding light on the dynamics of biomolecules. However, as important as they are, molecular dynamics studies focused on DNA repair pathways are still becoming more widespread. Currently, there are no review articles compiling the advancements made in molecular dynamics approaches applied to NER and discussing: (i) how this technique is currently employed in the field of DNA repair, focusing on NER proteins; (ii) which technical setups are being employed, their strengths and limitations; (iii) which insights or information are they providing to understand the NER pathway or NER-associated proteins; (iv) which open questions would be suited for this technique to answer; and (v) where can we go from here. These questions become even more crucial considering the numerous 3D structures published regarding the NER pathway's proteins in recent years. In this work, we tackle each one of these questions, revising and critically discussing the results published in the context of the NER pathway.

Single Molecule Force Spectroscopy Studies on Metalloproteins: Opportunities and Challenges

Metalloproteins play important roles in a wide range of biological processes. Elucidating the mechanisms via which metalloproteins fold and constitute their metal centers is critical to the understanding of the functions and dynamics of metalloproteins. Owing to its superior force and length resolution, single-molecule force spectroscopy (SMFS) has evolved into a powerful tool to probe the unfolding and folding mechanisms of metalloproteins at the single level by forcing metalloproteins to unfold and then refold along a reaction coordinate defined by the applied stretching force. The folding of metalloproteins is complex and involves two interwound processes, the folding of the polypeptide chain and the constitution of the metal center. Experimental studies of the folding of metalloproteins are challenging. SMFS studies have allowed researchers to directly probe the folding and unfolding of metalloproteins at the single-molecule level and the effect of metal centers on the folding-unfolding energy landscape of metalloproteins. New mechanistic insights on the folding and unfolding of some metalloproteins have been obtained, demonstrating the power and unique advantages that SMFS techniques may offer. In this Perspective, using calcium-binding proteins and small iron-sulfur proteins as examples, I provide a concise overview of the information and insights that SMFS studies have provided to understand the folding and unfolding of metalloproteins. I also discuss the opportunities and challenges that are present in this fast-progressing area of research.

Effects of External Perturbations on Protein Systems: A Microscopic View

Protein folding can be viewed as the origami engineering of biology resulting from the long process of evolution. Even decades after its recognition, research efforts worldwide focus on demystifying molecular factors that underlie protein structure-function relationships; this is particularly relevant in the era of proteopathic disease. A complex co-occurrence of different physicochemical factors such as temperature, pressure, solvent, cosolvent, macromolecular crowding, confinement, and mutations that represent realistic biological environments are known to modulate the folding process and protein stability in unique ways. In the current review, we have contextually summarized the substantial efforts in unveiling individual effects of these perturbative factors, with major attention toward bottom-up approaches. Moreover, we briefly present some of the biotechnological applications of the insights derived from these studies over various applications including pharmaceuticals, biofuels, cryopreservation, and novel materials. Finally, we conclude by summarizing the challenges in studying the combined effects of multifactorial perturbations in protein folding and refer to complementary advances in experiment and computational techniques that lend insights to the emergent challenges.

Role of Repeated Conformational Transitions in Substrate Binding of Adenylate Kinase

The catalytic cycle of the enzyme adenylate kinase involves large conformational motions between open and closed states. A previous single-molecule experiment showed that substrate binding tends to accelerate both the opening and the closing rates and that a single turnover event often involves multiple rounds of conformational switching. In this work, we showed that the repeated conformational transitions of adenylate kinase are essential for the relaxation of incorrectly bound substrates into the catalytically competent conformation by combining all-atom and coarse-grained molecular simulations. In addition, free energy calculations based on all-atom and coarse-grained models demonstrated that the enzyme with incorrectly bound substrates has much a lower free energy barrier for domain opening compared to that with the correct substrate conformation, which may explain the the acceleration of the domain opening rate by substrate binding. The results of this work provide mechanistic understanding to previous experimental observations and shed light onto the interplay between conformational dynamics and enzyme catalysis.

Effect of Ca2+ binding states of calmodulin on the conformational dynamics and force responses of myosin lever arm

The mechanochemical coupling and biological function of myosin motors are regulated by Ca2+ concentrations. As one of the regulation pathways, Ca2+ binding induces conformational change of the light chain calmodulin and its binding modes with myosin lever arm, which can affect the stiffness of the lever arm and force transmission. However, the underlying molecular mechanism of the Ca2+ regulated stiffness change is not fully understood. Here we study the effect of Ca2+ binding on the conformational dynamics and stiffness of the myosin VIIa lever arm bound with calmodulin by performing molecular dynamics simulations and dynamic correlation network analysis. The results showed that the calmodulin bound lever arm at apo state can sample three different conformations. In addition to the conformation observed in crystal structure, calmodulin bound lever arm at apo condition can also adopt another two conformations featured by different extents of small-angle bending of the lever arm. However, large-angle bending is strongly prohibited. Such results suggest that the calmodulin bound lever arm without Ca2+ binding is plastic for small-angle deformation but shows high stiffness for large-angle deformation. In comparison, after the binding of Ca2+, although the calmodulin bound lever arm is locally more rigid, it can adopt largely deformed or even unfolded conformations, which may render the lever arm incompetent for force transmission. The conformational plasticity of the lever arm for small-angle deformation at apo condition may be utilized as force buffer to prevent the lever arm from unfolding during the power stroke action of the motor domain.

Role of water-bridged interactions in metal ion coupled protein allostery

Allosteric communication between distant parts of proteins controls many cellular functions, in which metal ions are widely utilized as effectors to trigger the allosteric cascade. Due to the involvement of strong coordination interactions, the energy landscape dictating the metal ion binding is intrinsically rugged. How metal ions achieve fast binding by overcoming the landscape ruggedness and thereby efficiently mediate protein allostery is elusive. By performing molecular dynamics simulations for the Ca²⁺ binding mediated allostery of the calmodulin (CaM) domains, each containing two Ca²⁺ binding helix-loop-helix motifs (EF-hands), we revealed the key role of water-bridged interactions in Ca²⁺ binding and protein allostery. The bridging water molecules between Ca²⁺ and binding residue reduces the ruggedness of ligand exchange landscape by acting as a lubricant, facilitating the Ca²⁺ coupled protein allostery. Calcium-induced rotation of the helices in the EF-hands, with the hydrophobic core serving as the pivot, leads to exposure of hydrophobic sites for target binding. Intriguingly, despite being structurally similar, the response of the two symmetrically arranged EF-hands upon Ca²⁺ binding is asymmetric. Breakage of symmetry is needed for efficient allosteric communication between the EF-hands. The key roles that water molecules play in driving allosteric transitions are likely to be general in other metal ion mediated protein allostery.

Natural proteins often utilize allostery in executing a variety of functions. Metal ions are typical cofactors to trigger the allosteric cascade. In this work, using the Ca²⁺ sensor protein calmodulin as the model system, we revealed crucial roles of water-bridged interactions in the metal ion coupled protein allostery. The coordination of the Ca²⁺ to the binding site involves an intermediate in which the water molecule bridges the Ca²⁺ and the liganding residue. The bridging water reduces the free energy barrier height of ligand exchange, therefore facilitating the ligand exchange and allosteric coupling by acting as a lubricant. We also showed that the response of the two symmetrically arranged EF-hand motifs of CaM domains upon Ca²⁺ binding is asymmetric, which is directly attributed to the differing dehydration process of the Ca²⁺ ions and is needed for efficient allosteric communication.

Protein Assembly by Design

Proteins are nature's primary building blocks for the construction of sophisticated molecular machines and dynamic materials, ranging from protein complexes such as photosystem II and nitrogenase that drive biogeochemical cycles to cytoskeletal assemblies and muscle fibers for motion. Such natural systems have inspired extensive efforts in the rational design of artificial protein assemblies in the last two decades. As molecular building blocks, proteins are highly complex, in terms of both their three-dimensional structures and chemical compositions. To enable control over the self-assembly of such complex molecules, scientists have devised many creative strategies by combining tools and principles of experimental and computational biophysics, supramolecular chemistry, inorganic chemistry, materials science, and polymer chemistry, among others. Owing to these innovative strategies, what started as a purely structure-building exercise two decades ago has, in short order, led to artificial protein assemblies with unprecedented structures and functions and protein-based materials with unusual properties. Our goal in this review is to give an overview of this exciting and highly interdisciplinary area of research, first outlining the design strategies and tools that have been devised for controlling protein self-assembly, then describing the diverse structures of artificial protein assemblies, and finally highlighting the emergent properties and functions of these assemblies.

Phytotoxicity and accumulation of Cu in mature and young leaves of submerged macrophyte Hydrilla verticillata (L.f.) Royle

Phytotoxicity and accumulation of Cu in mature and young leaves of submerged macrophyte Hydrilla verticillata (L.f.) Royle were investigated by analyzing the chlorophyll contents, chloroplast ultrastructure and leaf surface structure under different Cu treatments (0, 0.01, 0.05 and 0.1 f mg L^-1). The results showed that 0.05 and 0.1 mg L^-1 Cu treatment decreased the contents of Chl a and Chl b, and caused damage on leaf surface structure and chloroplast ultrastructure compared with control (0 mg L^-1 Cu treatment). Higher concentration of Cu induced Chlorophyll decreases and the damages on the leaf surface structure and chloroplast ultrastructure were more pronounced in mature than in young leaves. It was observed that leaf Cu concentration increased almost linearly with exposure time and majority of the Cu accumulated in the cell walls. Among different cell wall fractions, the majority of Cu accumulated in cell walls was bound to the hemicellulose 1 and cellulose, followed by the pectin, hemicellulose 2. Mature leaves had significantly higher the concentrations of total Cu and bound-Cu in cell walls due to higher uronic acid content in their cell wall fractions (pectin, hemicellulose 1 and cellulose) than young leaves. Distinct cell wall composition might partially contribute to the different Cu toxicity and accumulation between mature and young leaves of submerged macrophyte H. verticillata. Our results show that mature leaves are more efficient in the uptake and accumulation of Cu than young leaves, which might explain why mature leaves sustain more severe damage.

Multimodal Response to Copper Binding in Superoxide Dismutase Dynamics

Copper/zinc superoxide dismutase (SOD) is a homodimeric metalloenzyme that has been extensively studied as a benchmark for structure-function relationships in proteins, in particular because of its implication in the familial form of the neurodegenerative disease amyotrophic lateral sclerosis. Here, we investigate microcrystalline preparations of two differently metalated forms of SOD, namely, the fully mature functional Cu,Zn state and the E,Zn-SOD state in which the Cu site is empty. By using solid-state NMR with fast magic-angle spinning (MAS) at high magnetic fields (¹H Larmor frequency of 800-1000 MHz), we quantify motions spanning a dynamic range from ns to ms. We determine that metal ion uptake does not act as a rigidification element but as a switch redistributing motional processes on different time scales, with coupling of the dynamics of histidine side chains and those of remote key backbone elements of the protein.

Role of the copper ion in pseudoazurin during the mechanical unfolding process

Metalloproteins require the corresponding metal cofactors to exert their proper function. The presence of metal cofactors in the metalloprotein makes it more difficult to investigate its folding and unfolding process. In this study, we employed atomic-force-microscopy-based single-molecule force spectroscopy to reveal the unfolding process of pseudoazurin (PAZ) that belongs to blue copper proteins. Our study shows that holo-PAZ requires a higher rupture force for mechanical unfolding comparing with the apo-PAZ. This result demonstrates that the copper atom not only enables PAZ access to transfer electron, but should also have an influence on its stability. The results also suggest that the electronic configuration of the metal cofactors has a striking effect on the strength of the organometallic bonds. Moreover, the results also reveal that there is an intermediate state during the unfolding process of PAZ. This study provides insight into the characteristics of metalloproteins and leads to a better knowledge of their interaction at the individual molecule level.

Investigations of the underlying mechanisms of HIF-1α and CITED2 binding to TAZ1

The TAZ1 domain of CREB binding protein is crucial for transcriptional regulation and recognizes multiple targets. The interactions between TAZ1 and its specific targets are related to the cellular hypoxic negative feedback regulation. Previous experiments reported that one of the TAZ1 targets, CITED2, is an efficient competitor of another target, HIF-1α. Here, by developing the structure-based models of TAZ1 complexes, we have uncovered the underlying mechanisms of the competitions between the two intrinsic disordered proteins (IDPs) HIF-1α and CITED2 binding to TAZ1. Our results support the experimental hypothesis on the competition mechanisms and the apparent affinity. Furthermore, the simulations locate the dominant position of forming TAZ1-CITED2 complex in both thermodynamics and kinetics. For thermodynamics, TAZ1-CITED2 is the lowest basin located on the free energy surface of binding in the ternary system. For kinetics, the results suggest that CITED2 binds to TAZ1 faster than HIF-1α. In addition, the analysis of contact map and Φ values is important for guiding further experimental studies to understand the biomolecular functions of IDPs.

Comparing Nonbonded Metal Ion Models in the Divalent Cation Binding Protein PsaA

Divalent metal cations are essential for many biological processes; however, accurately modeling divalent metal ions has proved a significant challenge for molecular dynamics force fields. Here we show that the choice of ion model influences the observed dynamics in PsaA, a metal binding protein from Streptococcus pneumoniae. We conduct extensive unbiased simulations and free energy calculations of PsaA bound to its cognate ligand Mn²⁺ and inhibitory ligand Zn²⁺ using three nonbonded ion models: a 12-6 model, a 12-6-4 model, and a multisite model. The observed coordination geometries and metal binding dynamics are sensitive to the choice of ion model, with the most dramatic differences observed in free energy calculations of ion release. We show that the conformational ensemble of Mn-bound PsaA is more similar to the crystallographic metal bound open state. This work extends the current model of PsaA metal binding and provides a framework for the rationalization of experimentally determined metal binding behavior. Our findings support the use of the 12-6-4 ion model for further simulations of divalent cation binding proteins.

X-ray snapshots reveal conformational influence on active site ligation during metalloprotein folding

Cytochrome c (cyt c) has long been utilized as a model system to study metalloprotein folding dynamics and the interplay between active site ligation and tertiary structure. However, recent reports regarding the weakness of the native Fe(ii)-S bond (Fe-Met80) call into question the role of the active site ligation in the protein folding process. In order to investigate the interplay between protein conformation and active site structures, we directly tracked the evolution of both during a photolysis-induced folding reaction using X-ray transient absorption spectroscopy and time-resolved X-ray solution scattering techniques. We observe an intermediate Fe-Met80 species appearing on ∼2 μs timescale, which should not be sustained without stabilization from the folded protein structure. We also observe the appearance of a new active site intermediate: a weakly interacting Fe-H2O state. As both intermediates require stabilization of weak metal-ligand interactions, we surmise the existence of a local structure within the unfolded protein that protects and limits the movement of the ligands, similar to the entatic state found in the native cyt c fold. Furthermore, we observe that in some of the unfolded ensemble, the local stabilizing structure is lost, leading to expansion of the unfolded protein structure and misligation to His26/His33 residues.

Importance of C-Terminal Extension in Thermophilic 1,4-α-Glucan Branching Enzyme from Geobacillus thermoglucosidans STB02

By sequence comparison, the majority of 1,4-α-glucan-branching enzymes (GBEs) consist of an N-terminal carbohydrate-binding domain, a TIM-barrel catalytic domain, and a C-terminal all-beta domain. Among these structures, the GBE from Geobacillus thermoglucosidans STB02 uniquely has a highly charged 26-amino-acid C-terminal extension, whose functional roles are the least understood. In this research, the functional significance of the C-terminal domain in GBE from G. thermoglucosidans STB02 and its extension were assessed using a C-terminal deletion analysis. Mutants lacking of more than 7 residues of the C-terminal all-beta domain could not be detected in lysates of their Escherichia coli expression strains, suggesting that an intact all-beta domain is required for structural stability. In contrast, truncation of the C-terminal extension resulted in greater stability and solubility than the wild type, as well as a lower sensitivity to the presence of added metal ions. Comparison of this mutant with the wild type suggests that the interaction of metal ions with the C-terminal extension influences performance of this enzyme.

Folding of copper proteins: role of the metal?

Copper is a redox-active transition metal ion required for the function of many essential human proteins. For biosynthesis of proteins coordinating copper, the metal may bind before, during or after folding of the polypeptide. If the metal binds to unfolded or partially folded structures of the protein, such coordination may modulate the folding reaction. The molecular understanding of how copper is incorporated into proteins requires descriptions of chemical, thermodynamic, kinetic and structural parameters involved in the formation of protein-metal complexes. Because free copper ions are toxic, living systems have elaborate copper-transport systems that include particular proteins that facilitate efficient and specific delivery of copper ions to target proteins. Therefore, these pathways become an integral part of copper protein folding in vivo. This review summarizes biophysical-molecular in vitro work assessing the role of copper in folding and stability of copper-binding proteins as well as protein-protein copper exchange reactions between human copper transport proteins. We also describe some recent findings about the participation of copper ions and copper proteins in protein misfolding and aggregation reactions in vitro.

Isophthaloyl-Based Selective Fluorescence Receptor for Zn (II) Ion in Semi-Aqueous Medium

A novel Isophthaloyl-based symmetrical (12E,21E)-N1',N3'-bis(2-hydroxybenzylidene) isophthalohydrazide, receptor (1) was synthesized and characterized using various spectroscopic technique. The reorganization ability of receptor (1) was evaluated in semi-aqueous medium and shows significant enhancement in fluorescence intensity for Zn (II) ion over various metal ions in CH₃CN:H₂O (1:1, v/v). The 1:2 binding stoichiometry between receptor (1) and Zn (II) ion was established using Job's plot and the proposed complex structure was calculated by applying Density Functional Theory (DFT) method. The binding constant (K_a) of receptor (1) with Zn (II) ion was established with the Benesi-Hildebrand plot, Scatchard and Connor's plot and the values are 1.00 × 10⁴ M^-1, 1.05× 10⁴ M^-1 and 1.05× 10⁴ M^-1 respectively. The limit of detection (LOD) and limit of quantification (LOQ) of receptor (1) and Zn (II) ion was 0.292 μM and 0.974 μM respectively. The binding mode was due to photo-induced electron transfer (PET) and the coordination of Zn (II) ion with C = N hydroxyl group of receptor (1). Electrochemical analysis of metal free receptor (1) and with Zn (II) ion also confirmed the formation of complex.

Binding of Copper Ions with Octapeptide Region in Prion Protein: Simulations with Charge Transfer Model

Copper ions are important cofactors of many metalloproteins. The binding dynamics of proteins to the copper ion is important for biological functions but is less understood at the microscopic level. What are the key factors determining the recognition and the stabilization of the copper ion during the binding? Our work investigates the binding dynamics of the copper ion with a simple system (the N-terminus of PrP) using simulation methods. To precisely characterize the protein?ion interaction, we build up an effective copper?peptide force field based on quantum chemistry calculations. In our model, the effects of charge transfer, protonation/deprotonation, and induced polarization are considered. With this force field, we successfully characterize the local structures and the complex interactions of the octapeptide around the copper ion. Furthermore, using an enhanced sampling method, the binding/unbinding processes of the copper ion with the octapeptide are simulated. Free-energy landscapes are generated in consequence, and multiple binding pathways are characterized. It is observed that various native ligands contribute differently to the binding processes. Some residues are related to the capture of the ion (behaving like ?arm?s), and some others contribute to the stabilization of the coordination structure (acting like ?core?s). These different interactions induce various pathways. Besides, a nonnative binding ligand is determined, and it has essential contributions and modulations to the binding pathways. With all these results, the picture of copper?octapeptide binding is outlined. These features are believed to happen in many ion?peptide interactions, such as the cooperative stabilization between the coordinations with neighboring backbone nitrogens and an auxiliary intermediate coordination with the neighboring oxygen from the N-terminal direction. We believe that our studies are valuable to understand the complicated ion?peptide binding processes.

Decoding the Time-Dependent Response of Bioluminescent Metal-Detecting Whole-Cell Bacterial Sensors

The signal produced by aqueous dispersions of bioluminescent, metal-responsive whole-cell bacterial sensors is indicative of the concentration of bioavailable metal ions in solution. The conventional calibration-based strategy followed for measuring this concentration is however inadequate to provide any quantitative prediction of the cell response over time as a function of, e.g., their growth features, their defining metal accumulation properties, or the physicochemical medium composition. Such an evaluation is still critically needed for assessing on a mechanistic level the performance of biosensors in terms of metal bioavailability and toxicity monitoring. Herein we report a comprehensive formalism unraveling how the dependence of bioluminescence on time is governed by the dynamics of metal biouptake, by the activation kinetics of lux-based reporter gene, and by the ensuing rate of luciferase production, the kinetics of light emission, and quenching. It is shown that the bioluminescence signal corresponds to the convolution product between two time-dependent functions, one detailing the dynamic interplay of the above micro- and nanoscale processes, and the other pertaining to the change in concentration of photoactive cell sensors over time. Numerical computations illustrate how the shape and magnitude of the bioluminescence peak(s) are intimately connected to the dependence of the photoactive cell concentration on time and to the magnitudes of Deborah numbers that compare the relevant time scales of the biointerfacial and intracellular events controlling light emission. Explicit analytical expressions are further derived for practical situations where bioluminescence is proportional to the concentration of metal ions in solution. The theory is further quantitatively supported by experiments performed on luminescent cadmium-responsive lux-based Escherichia coli biosensors.

Can multiscale simulations unravel the function of metallo-enzymes to improve knowledge-based drug discovery?

Metallo-enzymes are a large class of biomolecules promoting specialized chemical reactions. Quantum-classical quantum mechanics/molecular mechanics molecular dynamics, describing the metal site at quantum mechanics level, while accounting for the rest of system at molecular mechanics level, has an accessible time-scale limited by its computational cost. Hence, it must be integrated with classical molecular dynamics and enhanced sampling simulations to disentangle the functions of metallo-enzymes. In this review, we provide an overview of these computational methods and their capabilities. In particular, we will focus on some systems such as CYP19A1 a Fe-dependent enzyme involved in estrogen biosynthesis, and on Mg2+-dependent DNA/RNA processing enzymes/ribozymes and the spliceosome, a protein-directed ribozyme. This information may guide the discovery of drug-like molecules and genetic manipulation tools.

Real-time monitoring β-lactam/β-lactamase inhibitor (BL/BLI) mixture towards the bacteria porin pathway at single molecule level

Multidrug-resistant bacteria are a great concern and a problem that must be addressed. Extended-spectrum β-lactamases are a common defence mechanism of bacteria to make β-lactam (BL) antibiotics ineffective. β-Lactamase inhibitors (BLIs) are consequently designed and are often clinically prescribed with a BL antibiotic to hinder degradation. Current studies focusing on how BL antibiotics or BLIs interact solely with the bacterial outer membrane nanopores (porins) on reaching the periplasmic side using a nanopore-based sensing technique. In electrochemical studies, the bias voltage allows real-time monitoring of BL antibiotics, BLIs and their mixture through the porin pathway at the single-molecule level. Here we consider the most abundant membrane protein from Escherichia coli (i.e. OmpF), purify and reconstitute the membrane protein in an artificial lipid bilayer and then study its ex vivo electrochemical behaviour. We show the piperacillin/tazobactam mixture interacts with OmpF, whereas the substrate interacts under the maximum bandwidth. The power spectrum analysis of the ionic current trace demonstrates the ampicillin/sulbactram mixture requires more energy than ampicillin alone to pass through the porin pathway. Our results demonstrate that clinically relevant combinations (e.g. piperacillin/tazobactam and ampicillin/sulbactam) interact more strongly with OmpF than either the BL antibiotic or the BLI alone. We suggest a quick and relatively cheap screening method to test the ability of BL antibiotics/BLIs to cross the bacterial cellular membrane.

Calcium-Mediated Allostery of the EGF Fold

The epidermal growth factor (EGF)-like domain is one of the most abundant disulfide-containing domains in nature and is involved in many cellular processes critical to life. Although many EGF-like domains participate in calcium-dependent functions by responding to the local calcium concentration, little is known about how this responsiveness is programmed at the molecular level. Here, we reveal the structural and environmental determinants underpinning the folding of a synthetic analogue of the EGF-A domain (from the low-density lipoprotein receptor). We show that calcium sensitivity is enabled by an allosteric folding pathway, in which calcium binding is connected to the peptide core through local inter-residue interactions. In the absence of calcium, the fold favors disorder because the inherently weak core is insufficient to stabilize the active form, resulting in substantial loss in activity of 2 orders of magnitude. The EGF-A fold, which can freely transition between active and disordered states, is volatile, and we found it to be intolerant of mutations, unlike other disulfide-rich peptides that have been used as stabilizing frameworks. This volatility is beneficial for modularity/plasticity and appears to have evolved for such a purpose, allowing cellular pathways to sense and respond to environmental cues.

Probing Cytochrome c Folding Transitions upon Phototriggered Environmental Perturbations Using Time-Resolved X-ray Scattering

Direct tracking of protein structural dynamics during folding-unfolding processes is important for understanding the roles of hierarchic structural factors in the formation of functional proteins. Using cytochrome c (cyt c) as a platform, we investigated its structural dynamics during folding processes triggered by local environmental changes (i.e., pH or heme iron center oxidation/spin/ligation states) with time-resolved X-ray solution scattering measurements. Starting from partially unfolded cyt c, a sudden pH drop initiated by light excitation of a photoacid caused a structural contraction in microseconds, followed by active site restructuring and unfolding in milliseconds. In contrast, the reduction of iron in the heme via photoinduced electron transfer did not affect conformational stability at short timescales (<1 ms), despite active site coordination geometry changes. These results demonstrate how different environmental perturbations can change the nature of interaction between the active site and protein conformation, even within the same metalloprotein, which will subsequently affect the folding structural dynamics.

Cadmium-induced conformational changes in type 2 metallothionein of medicinal plant Coptis japonica: insights from molecular dynamics studies of apo, partially and fully metalated forms

Plants play an important role in the removal of excess heavy metals from soil and water. Medicinal plants can also have non-traditional use in phytoremediation technologies. Among the heavy metals, Cadmium (Cd) is the most abundant and readily taken up by the crop plants. Plant metallothioneins (MTs) are small proteins having cysteine-rich residues and appear to play key roles in metal homoeostasis. Plant metallothionein 2 (MT 2) from Coptis japonica (Gold-thread; CjMT 2) is a typical member of this subfamily and features two cysteine-rich regions containing eight and six cysteine residues, respectively, separated by 42 amino acids long linker region. In-silico analysis of MT 2 protein sequences of C. japonica was performed. In this study, ab initio methods were utilised for the prediction of three-dimensional structure of CjMT 2. After structure validation, heavy metal-binding sites were predicted for the selected modelled structures of CjMT 2. To obtain Cd_i-CjMT 2 (i = 1-7), metalated complex individual docking experiments were performed. The stability of the metalated docked structures was assessed by molecular dynamics (MD) simulation studies. Our study showed that CjMT 2 binds up to 4 Cd²⁺ ions in two distinct domains: a N-terminal β-domain that binds to 2 Cd²⁺ ions and a C-terminal α-domain that binds with 2 Cd²⁺ ions. Our analysis revealed that Cys residues of alpha and beta domain and some residues of spacer region of CjMT 2 protein might be important for the cadmium interaction. MD simulation studies provided insight into metal-induced conformational changes and mechanism of metalation of CjMT 2, an intrinsically disordered protein. This study provides useful insights into mechanism of cadmium-type 2 metallothionein interaction.

Metal cofactor modulated folding and target recognition of HIV-1 NCp7

The HIV-1 nucleocapsid 7 (NCp7) plays crucial roles in multiple stages of HIV-1 life cycle, and its biological functions rely on the binding of zinc ions. Understanding the molecular mechanism of how the zinc ions modulate the conformational dynamics and functions of the NCp7 is essential for the drug development and HIV-1 treatment. In this work, using a structure-based coarse-grained model, we studied the effects of zinc cofactors on the folding and target RNA(SL3) recognition of the NCp7 by molecular dynamics simulations. After reproducing some key properties of the zinc binding and folding of the NCp7 observed in previous experiments, our simulations revealed several interesting features in the metal ion modulated folding and target recognition. Firstly, we showed that the zinc binding makes the folding transition states of the two zinc fingers less structured, which is in line with the Hammond effect observed typically in mutation, temperature or denaturant induced perturbations to protein structure and stability. Secondly, We showed that there exists mutual interplay between the zinc ion binding and NCp7-target recognition. Binding of zinc ions enhances the affinity between the NCp7 and the target RNA, whereas the formation of the NCp7-RNA complex reshapes the intrinsic energy landscape of the NCp7 and increases the stability and zinc affinity of the two zinc fingers. Thirdly, by characterizing the effects of salt concentrations on the target RNA recognition, we showed that the NCp7 achieves optimal balance between the affinity and binding kinetics near the physiologically relevant salt concentrations. In addition, the effects of zinc binding on the inter-domain conformational flexibility and folding cooperativity of the NCp7 were also discussed.

Unified understanding of folding and binding mechanisms of globular and intrinsically disordered proteins

Extensive experimental and theoretical studies have advanced our understanding of the mechanisms of folding and binding of globular proteins, and coupled folding and binding of intrinsically disordered proteins (IDPs). The forces responsible for conformational changes and binding are common in both proteins; however, these mechanisms have been separately discussed. Here, we attempt to integrate the mechanisms of coupled folding and binding of IDPs, folding of small and multi-subdomain proteins, folding of multimeric proteins, and ligand binding of globular proteins in terms of conformational selection and induced-fit mechanisms as well as the nucleation–condensation mechanism that is intermediate between them. Accumulating evidence has shown that both the rate of conformational change and apparent rate of binding between interacting elements can determine reaction mechanisms. Coupled folding and binding of IDPs occurs mainly by induced-fit because of the slow folding in the free form, while ligand binding of globular proteins occurs mainly by conformational selection because of rapid conformational change. Protein folding can be regarded as the binding of intramolecular segments accompanied by secondary structure formation. Multi-subdomain proteins fold mainly by the induced-fit (hydrophobic collapse) mechanism, as the connection of interacting segments enhances the binding (compaction) rate. Fewer hydrophobic residues in small proteins reduce the intramolecular binding rate, resulting in the nucleation–condensation mechanism. Thus, the folding and binding of globular proteins and IDPs obey the same general principle, suggesting that the coarse-grained, statistical mechanical model of protein folding is promising for a unified theoretical description of all mechanisms.

Metal Ion Modeling Using Classical Mechanics

Metal ions play significant roles in numerous fields including chemistry, geochemistry, biochemistry, and materials science. With computational tools increasingly becoming important in chemical research, methods have emerged to effectively face the challenge of modeling metal ions in the gas, aqueous, and solid phases. Herein, we review both quantum and classical modeling strategies for metal ion-containing systems that have been developed over the past few decades. This Review focuses on classical metal ion modeling based on unpolarized models (including the nonbonded, bonded, cationic dummy atom, and combined models), polarizable models (e.g., the fluctuating charge, Drude oscillator, and the induced dipole models), the angular overlap model, and valence bond-based models. Quantum mechanical studies of metal ion-containing systems at the semiempirical, ab initio, and density functional levels of theory are reviewed as well with a particular focus on how these methods inform classical modeling efforts. Finally, conclusions and future prospects and directions are offered that will further enhance the classical modeling of metal ion-containing systems.

Protein folding in the cell envelope of Escherichia coli

While the entire proteome is synthesized on cytoplasmic ribosomes, almost half associates with, localizes in or crosses the bacterial cell envelope. In Escherichia coli a variety of mechanisms are important for taking these polypeptides into or across the plasma membrane, maintaining them in soluble form, trafficking them to their correct cell envelope locations and then folding them into the right structures. The fidelity of these processes must be maintained under various environmental conditions including during stress; if this fails, proteases are called in to degrade mislocalized or aggregated proteins. Various soluble, diffusible chaperones (acting as holdases, foldases or pilotins) and folding catalysts are also utilized to restore proteostasis. These responses can be general, dealing with multiple polypeptides, with functional overlaps and operating within redundant networks. Other chaperones are specialized factors, dealing only with a few exported proteins. Several complex machineries have evolved to deal with binding to, integration in and crossing of the outer membrane. This complex protein network is responsible for fundamental cellular processes such as cell wall biogenesis; cell division; the export, uptake and degradation of molecules; and resistance against exogenous toxic factors. The underlying processes, contributing to our fundamental understanding of proteostasis, are a treasure trove for the development of novel antibiotics, biopharmaceuticals and vaccines.

Investigation of the early stages of human γD-crystallin aggregation process

Cataract, a major cause of visual impairment worldwide, is a common disease of the eye lens related to protein aggregation. Several factors including the exposure of ultraviolet irradiation and possibly acidic condition may induce the unfolding and subsequent aggregation of the crystallin proteins leading to crystalline lens opacification. Human γD-crystallin (HγDC), a 173 residue monomeric protein, abundant in the nucleus of the human eye lens, has been shown to aggregate and form amyloid fibrils under acidic conditions and that this aggregation route is thought to be a potential initiation pathway for the onset of age-related nuclear cataract. However, the underlying mechanism of fibril formation remains elusive. This report is aimed at examining the structural changes and possible amyloid fibril formation pathway of HγDC using molecular dynamics and molecular docking simulations. Our findings demonstrated that incubation of HγDC under the acidic condition redistributes the protein surface charges and affects the protein interaction with its surrounding solvent environment. This brings about a twist motion in the overall tertiary structure that gives rise to newly formed anti-parallel β-strands in the C-terminal flexible loop regions. The change in protein structural conformation also involves an alteration in specific salt-bridge interactions. Altogether, these findings revealed a plausible mechanism for amyloid fibril formation of HγDC that is important to the early stages of HγDC aggregation involved in cataractogenesis.

Effects of heme binding on myoglobin folding: Coarse grained molecular simulations

Many proteins contain cofactors, such as heme, ATP and metal ions. Binding of cofactors is not only essential for their biological functions, but also can reshape the intrinsic energy landscape of protein molecules and modulate the folding and stability. However, the molecular mechanism of cofactor coupled protein folding is not well understood. In this work, we study the cofactor coupled folding of myoglobin, which is a typical cofactor (heme) containing protein, by performing molecular dynamics simulations with a structure-based protein model developed based on the energy landscape theory. We showed that the heme binding increases the stability of the myoglobin. More importantly, the heme binding tends to increase the protein folding cooperativity, and switch the folding process from a “three-state” mechanism to a “two-state” mechanism. We also showed that the folding pathways of the myoglobin can be modulated by the heme binding. By performing comparative simulations, we revealed that the above effects of heme binding are resulted from the heme induced folding of F-helix, which is otherwise unstructured at apo state, and the heme mediated contacting interactions around the heme binding site. The simulation results are consistent with available experimental data, and provide insights into the molecular mechanism of the effects of cofactor binding on the protein folding and stability.

The Evolution of New Catalytic Mechanisms for Xenobiotic Hydrolysis in Bacterial Metalloenzymes

An increasing number of bacterial metalloenzymes have been shown to catalyse the breakdown of xenobiotics in the environment, while others exhibit a variety of promiscuous xenobiotic-degrading activities. Several different evolutionary processes have allowed these enzymes to gain or enhance xenobiotic-degrading activity. In this review, we have surveyed the range of xenobiotic-degrading metalloenzymes, and discuss the molecular and catalytic basis for the development of new activities. We also highlight how our increased understanding of the natural evolution of xenobiotic-degrading metalloenzymes can be been applied to laboratory enzyme design.

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

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