Conformational thermodynamics of metal-ion binding to a protein

In-silicon studies on hydration in EcoRI-cognate DNA complex

Restriction endonucleases (REs) cleave DNA at specific site in presence of Mg2+ ion. Experiments further emphasize the role of hydration in metal ion specificity and sequence specificity of DNA cleavage. However, the relation between hydration and specificity has not been understood till date. This leads us to study via all-atom molecular dynamics (MD) simulations how the hydration around the scissile phosphate group changes in presence of Mg2+ and Ca2+ and depend on the DNA sequence. We observe the least number of hydrogen bonds around the scissile phosphate group in presence of Mg2+ ion. We further find that the hydrogen bonds decrease at the scissile phosphate on mutating one base pair in the cleavage region of the DNA in Mg2+ loaded EcoRI-DNA complex. We also perform steered MD simulations and observe that the rate of decrease of fraction of hydrogen bonds is slower in the mutated complex than the unmutated complex.

New insight into scale inhibition during tea brewing: Ca2+/Mg2+ complexing and alkalinity consumption

Scale not only affects the taste and color of water, but also increases the risks of osteoporosis and cardiovascular diseases associated with drinking it. As a popular beverage, tea is rich many substances that have considerable potential for scale inhibition, including protein, tea polyphenols and organic acids. In this study, the effect of tea brewing on scale formation was explored. It was found that the proteins, catechins and organic acids in tea leaves could be released when the green tea was brewed in water with sufficient hardness and alkalinity. The tea-released protein was able to provide carboxyl groups to chelate with calcium ions (Ca²⁺), preventing the Ca²⁺ from reacting with the carbonate ions (CO₃^2-). The B rings of catechins were another important structure in the complexation of Ca²⁺ and magnesium ions (Mg²⁺). The carboxyl and hydroxyl groups on the organic acids was able to form five-membered chelating rings with Ca²⁺ and Mg²⁺, resulting in a significant decrease in Ca²⁺ from 100.0 to 60.0 mg/L. Additionally, the hydrogen ions (H⁺) provided by the organic acids consumed and decreased the alkalinity of the water from 250.0 to 131.4 mg/L, leading to a remarkable reduction in pH from 8.93 to 7.73. It further prevented the bicarbonate (HCO₃^-) from producing CO₃^2- when the water was heated. The reaction of the tea constituents with the hardness and alkalinity inhibited the formation of scale, leading to a significant decrease in turbidity from 10.6 to 1.4 NTU. Overall, this study provides information to help build towards an understanding of the scale inhibition properties of tea and the prospects of tea for anti-scaling in industrial applications.

Effect on the conformations of the spike protein of SARS-CoV-2 due to mutation

The spike protein of SARS-CoV-2 mediates receptor binding and cell entry and is the key immunogenic target for virus neutralization and the present attention of many vaccine layouts. It exhibits significant conformational flexibility. We study the structural fluctuations of spike protein among the most common mutations that appeared in the variant of concerns (VOC). We report the thermodynamics of conformational changes in mutant spike protein with respect to the wild-type from the distributions of the dihedral angles obtained from the equilibrium configurations generated via all-atom molecular dynamics simulations. We find that the mutation causes the increase in distance between the N-terminal domain and receptor binding domain, leading to an obtuse angle cosine θ distribution in the trimeric structure in spike protein. Thus, an increase in open state is conferred to the more infectious variants of SARS-CoV-2. The thermodynamically destabilized and disordered residues of receptor binding motif among the mutant variants of spike protein are proposed to serve as better binding sites for the host factor. We identify a short stretch of region connecting the N-terminal domain and receptor binding domain forming a linker loop where many residues undergo stabilization in the open state compared to the closed one.

Remediation and Mechanisms of Cadmium Biosorption by a Cadmium-Binding Protein from Lentinula edodes

Cadmium bioremediation with metal-binding proteins is primarily conducted using metallothioneins (MTs). However, in the present study, we investigated a non-MT cadmium-binding protein from Lentinula edodes (LECBP) as a remediation tool for cadmium biosorption in Escherichia coli. The results indicated that the expression of LECBP significantly enhanced the cadmium biosorption capacity of transgenic E. coli. The secondary structure and conformation of LECBP were changed after binding with cadmium as evidenced by circular dichroism and fluorescence spectroscopy. The results of Fourier transform infrared spectroscopy indicated that carboxyl oxygen and amino nitrogen atoms were involved in the interaction between LECBP and cadmium. The results further demonstrated that glutamic acid and histidine residues are the potential binding sites. Our results have thus provided new insights into cadmium bioremediation in an aquatic environment.

Structural and dynamic responses of calcium ion binding loop residues in metallo-proteins

Conformational changes in bio-molecular systems are fundamental to several biological processes. It is important to study changes in responses of underlying microscopic variables, like dihedral angles as conformational change takes place. We perform all-atom simulations and modelling via Langevin equation to illustrate the changes in structural and dynamic responses of dihedral angles of calcium ion binding residues of different proteins in metal ion free (apo) and bound (holo) states. The equilibrium distributions of dihedral angles in apo- and holo-states represent structural response. Our studies show the presence of dihedrals with multiple peaks (isomeric states) separated by barrier heights is more frequent in apo- than in holo-state. The relaxation time-scale of dihedral fluctuations is found to increase linearly with decreasing barrier height due to more frequent barrier re-crossing events. The slow kinetic response of the dihedrals also contributes to slowing down of macro-scale fluctuations, which may be useful to understand kinetics of various bio-molecular processes.

What factors determine the stability of a weak protein-protein interaction in a charged aqueous droplet?

Maintaining the interface of a weak transient protein complex transferred from bulk solution to the gaseous state via evaporating droplets is a critical question in the detection of the complex association (dissociation) constant by using electrospray ionization mass spectrometry (ESI-MS). Here we explore the factors that may affect the stability of a protein-protein interaction (PPI) using atomistic molecular dynamics (MD) modelling of a complex of ubiquitin (Ub) and the ubiquitin-associated domain (UbA) (RCSB PDB code ) and a non-covalent complex of diubiquitin (RCSB PDB code ) in aqueous droplets. A general method is presented to determine the protonation states of the complexes we investigate in particular, and that of a protein in general, under various pH conditions that an evaporating droplet acquires due to its change in size. We find that the combination of high temperature and high charge states of the protein complexes may destabilize the interface by creating new interfaces instead of a direct rupture of the initial stable interface. We provide evidence that highly charged protein complexes are found in droplets that form conical extrusions of the solvent on the surface due to charge-induced instability. This distinct droplet morphology leads to a higher solvent evaporation rate that assists in transferring the complex in the gaseous state without dissociation. The conical solvent protrusions expose on the droplet surface certain amino acids that otherwise would be solvated in a droplet with the protein complex of low charge states. The new vapor-protein interface does not have a direct effect on the stability of the PPI. A common way in experiments to stabilize the protein complexes in droplets is to reduce the protonation state of the proteins. Here we find that weakly bound protein complexes even at high protonation states can be stabilized by the presence of a small number of counterions, without affecting the protonation state of the protein. Our findings may provide guiding principles in ESI-MS experiments to stabilize weak transient PPIs.

Allostery in Orai1 binding to calmodulin revealed from conformational thermodynamics

Here, we study microscopic mechanism of complex formation between Ca²⁺-bound calmodulin (holoCaM) and Orai1 that regulates Ca²⁺-dependent inactivation process in eukaryotic cells. We compute conformational thermodynamic changes in holoCaM with respect to complex of Orai1 bound to C-terminal domain of holoCaM using histograms of dihedral angles of the proteins over trajectories from molecular dynamics simulations. Our analysis shows that the N-terminal domain residues L4, T5, Q41, N42, T44 and E67 of holoCaM get destabilized and disordered due to Orai1 binding to C-terminal domain of calmodulin affect the N-terminal domain residues. Among these residues, polar T44, having maximum destabilization and disorder via backbone fluctuations, shows the largest change in solvent exposure. This suggests that N-terminal domain is allosterically regulated via T44 by the binding of Orai1 to the C-terminal domain.

Molecular dynamics studies on conformational thermodynamics of Orai1-calmodulin complex

Molecular understanding of bio-macromolecular binding is a challenging task due to large sizes of the molecules and presence of variety of interactions. Here, we study the molecular mechanism of calmodulin (CaM) binding to Orai1 that regulates Ca²⁺-dependent inactivation process in eukaryotic cells. Although experimental observations indicate that Orai1 binds to the C-terminal of Ca²⁺-loaded CaM, it is not decisive if N-domain of CaM interacts with Orai1. We address the issue of interaction of different domains of CaM with Orai1 using conformational thermodynamic changes, computed from histograms of dihedral angles over simulated trajectories of CaM, CaM-binding domain of Orai1 and complexes of CaM with Orai1. The changes for all residues of both C and N terminal domains of CaM upon Orai1 binding are compared. Our analysis shows that Orai1binds to both C-terminal and N-terminal domains of CaM, indicating 1:2 stoichiometry. The Orai1 binding to N-terminal domain of CaM is less stable than that to the C-terminal domain. The binding residues are primarily hydrophobic. These observations are in qualitative agreement to the experiments. The conformational thermodynamic changes thus provide a useful computational tool to provide atomic details of interactions in bio-macromolecular binding.

Polyanionic Antimicrobial Membranes: An Experimental and Theoretical Study

Polycationic polymers have been widely used as antimicrobial materials because of their broad spectrum activity and potential use as new antibiotics. Herein, we report the synthesis of polyanionic antimicrobial membranes by in situ photo-cross-linking of a sulfate based anionic monomer, followed by cation-exchange with organic (quaternary ammonium or imidazolium) or metal (Ag⁺, Cu²⁺, Fe³⁺, Zn²⁺, Na⁺, K⁺) cations. The resultant polyanionic membranes show high and broad spectrum antibacterial activities against both bacteria (Escherichia coli, Staphylococcus aureus) and fungi (Candida albicans ). In addition, the polyanionic antimicrobial membranes efficiently inhibited the formation of biofilms by SC5314 and its crk1 gene deleted (Δcrk1) C. albicans strains. Furthermore, the synthesized polyanionic membranes exhibit good blood compatibility, low cytotoxicity and long-term antibacterial stability, demonstrating safe antimicrobial materials in the application of healthcare.

SDS induced dissociation of STY3178 oligomer: experimental and molecular dynamics studies

STY3178 the yfdX Salmonella Typhi protein dissociates reversibly in presence of sodium dodecyl sulphate from trimer to monomer.

Conformational thermodynamics guided structural reconstruction of biomolecular fragments

Conformational thermodynamics compares the modeling protocols to identify the conformation of the missing region leading to a suitable model for metal ion free (apo) skeletal muscle Troponin C.

Quantum chemical studies on nucleophilic sites in calcium ion bound zwitterionic calmodulin loops

Quantum chemical calculation on Ca²⁺ bound zwitterionic calmodulin-loops shows terminal capping contributions localized near HOMO and LUMO, which decay exponentially and presence of nucleophilic site at the phenyl-ring O_h of Y99.

A microscopic insight from conformational thermodynamics to functional ligand binding in proteins

We show that the thermodynamics of metal ion-induced conformational changes aid to understand the functions of protein complexes. This is illustrated in the case of a metalloprotein, alpha-lactalbumin (aLA), a divalent metal ion binding protein. We use the histograms of dihedral angles of the protein, generated from all-atom molecular dynamics simulations, to calculate conformational thermodynamics. The thermodynamically destabilized and disordered residues in different conformational states of a protein are proposed to serve as binding sites for ligands. This is tested for β-1,4-galactosyltransferase (β4GalT) binding to the Ca(2+)-aLA complex, in which the binding residues are known. Among the binding residues, the C-terminal residues like aspartate (D) 116, glutamine (Q) 117, tryptophan (W) 118 and leucine (L) 119 are destabilized and disordered and can dock β4GalT onto Ca(2+)-aLA. No such thermodynamically favourable binding residues can be identified in the case of the Mg(2+)-aLA complex. We apply similar analysis to oleic acid binding and predict that the Ca(2+)-aLA complex can bind to oleic acid through the basic histidine (H) 32 of the A2 helix and the hydrophobic residues, namely, isoleucine (I) 59, W60 and I95, of the interfacial cleft. However, the number of destabilized and disordered residues in Mg(2+)-aLA are few, and hence, the oleic acid binding to Mg(2+)-bound aLA is less stable than that to the Ca(2+)-aLA complex. Our analysis can be generalized to understand the functionality of other ligand bound proteins.