Molecular dynamics simulations on interaction between bacterial proteins: Implication on pathogenic activities

Decoding the effect of temperatures on conformational stability and order of ligand unbound thermosensing adenine riboswitch using molecular dynamics simulation

The structure-function relationship of the riboswitch is governed mainly by two factors, ligand binding and temperature. Most of the experimental studies shed light on structural dynamics and gene regulation function of Adenine riboswitch from the aspect of ligand instead of temperature. Two unliganded Adenine riboswitch conformations (apoA and apoB) from the thermophile Vibrio vulnificus draw particular attention due to their diverse and polymorphic structures. Ligand-free apoB Adenine riboswitch conformation is not able to interact with the ligand whereas ligand-free apoA Adenine riboswitch conformation adopts ligand-receptive form. The interconversion between apoA and apoB conformation is temperature-dependent and thermodynamically controlled. Therefore Adenine riboswitch is called temperature sensing RNA. The molecular mechanism underlying the thermosensitivity of ligand free Adenine riboswitch is not well known. Hence an attempt is made to examine conformational stability and order of apoA with respect to apoB Adenine riboswitch aptamer computing RMSD, RMSF, RG, principal component analysis, hydrogen bonding interaction and conformational thermodynamics derived from all-atom molecular dynamics trajectories in the temperature range 283K-400K. The temperatures corresponding to the conformational stability and order of apoA adenine riboswitch with respect to apoB adenine riboswitch whole aptamer are shown in descending order 293K∼303K> 313K∼283K>373K>323K. Residue wise and domain wise changes in conformational free energy and entropy of conformational degrees of freedom like pseudo-torsion angle ƞ and θ reflect apoA exhibits pronounced conformational stability compared to apoB at temperatures 293K and 303K, whereas both forms reveal decreased stability at 323K and 400K and may be inactivated, highlighting their role as temperature sensors.

Fluctuation-Dominated Ligand Binding in Molten Globule Protein

A molten globule (MG) state is an intermediate state of protein observed during the unfolding of the native structure. In MG states, milk protein α-lactalbumin (aLA) binds to oleic acid (OLA). This MG-aLA-OLA complex, popularly known as XAMLET, performs cytotoxic activities against cancer cell lines. However, the microscopic understanding of ligand recognition ability in the MG state of the protein has not yet been explored. Motivated by this, we explore the binding of bovine aLA with OLA using all-atom molecular dynamics (MD) simulations. We find the binding mode between MG-aLA and OLA using the conformational thermodynamics method. We also estimate the binding free energy using the umbrella sampling (US) method for both the MG state and the neutral state. We find that the binding free energy obtained from US is comparable with earlier experimental results. We characterize the dihedral fluctuations as the ligand is liberated from the active site of the protein using steered MD. The low energy fluctuations occur near the ligand binding site, which eventually transfer toward the Ca2+-binding site as the ligand is taken away from the protein.

Microscopic insight to specificity of metal ion cofactor in DNA cleavage by restriction endonuclease EcoRV

Restriction endonucleases protect bacterial cells against bacteriophage infection by cleaving the incoming foreign DNA into fragments. In presence of Mg²⁺ ions, EcoRV is able to cleave the DNA but not in presence of Ca²⁺ , although the protein binds to DNA in presence of both metal ions. We make an attempt to understand this difference using conformational thermodynamics. We calculate the changes in conformational free energy and entropy of conformational degrees of freedom, like DNA base pair steps and dihedral angles of protein residues in Mg²⁺ (A)-EcoRV-DNA complex compared to Ca²⁺ (S)-EcoRV-DNA complex using all-atom molecular dynamics (MD) trajectories of the complexes. We find that despite conformational stability and order in both complexes, the individual degrees of freedom behave differently in the presence of two different metal ions. The base pairs in cleavage region are highly disordered in Ca²⁺ (S)-EcoRV-DNA compared to Mg²⁺ (A)-EcoRV-DNA. One of the acidic residues ASP90, coordinating to the metal ion in the vicinity of the cleavage site, is conformationally destabilized and disordered, while basic residue LYS92 gets conformational stability and order in Ca²⁺ (S) bound complex than in Mg²⁺ (A) bound complex. The enhanced fluctuations hinder placement of the metal ion in the vicinity of the scissile phosphate of DNA. Similar loss of conformational stability and order in the cleavage region is observed by the replacement of the metal ion. Considering the placement of the metal ion near scissile phosphate as requirement for cleavage action, our results suggest that the changes in conformational stability and order of the base pair steps and the protein residues lead to cofactor sensitivity of the enzyme. Our method based on fluctuations of microscopic conformational variables can be applied to understand enzyme activities in other protein-DNA systems.

Computer-aided synthesis of dapsone-phytochemical conjugates against dapsone-resistant Mycobacterium leprae

Leprosy continues to be the belligerent public health hazard for the causation of high disability and eventual morbidity cases with stable prevalence rates, even with treatment by the on-going multidrug therapy (MDT). Today, dapsone (DDS) resistance has led to fear of leprosy in more unfortunate people of certain developing countries. Herein, DDS was chemically conjugated with five phytochemicals independently as dapsone-phytochemical conjugates (DPCs) based on azo-coupling reaction. Possible biological activities were verified with computational chemistry and quantum mechanics by molecular dynamics simulation program before chemical synthesis and spectral characterizations viz., proton-HNMR, FTIR, UV and LC-MS. The in vivo antileprosy activity was monitored using the 'mouse-foot-pad propagation method', with WHO recommended concentration 0.01% mg/kg each DPC for 12 weeks, and the host-toxicity testing of the active DPC4 was seen in cultured-human-lymphocytes in vitro. One-log bacilli cells in DDS-resistant infected mice footpads decreased by the DPC4, and no bacilli were found in the DDS-sensitive mice hind pads. Additionally, the in vitro host toxicity study also confirmed that the DCP4 up to 5,000 mg/L level was safety for oral administration, since a minor number of dead cells were found in red color under a fluorescent microscope. Several advanced bioinformatics tools could help locate the potential chemical entity, thereby reducing the time and resources required for in vitro and in vitro tests. DPC4 could be used in place of DDS in MDT, evidenced from in vivo antileprosy activity and in vitro host toxicity study.

Structural and Physiological Exploration of Salmonella Typhi YfdX Uncovers Its Dual Function in Bacterial Antibiotic Stress and Virulence

YfdX is a prokaryotic protein encoded by several pathogenic bacteria including Salmonella enterica serovar Typhi, which causes one of the most fatal infectious diseases, typhoid fever. YfdX is a product of the yfdXWUVE operon and is known to be under the control of EvgA, a regulator protein controlling the expression of several proteins involved in response to environmental stress, in Escherichia coli. Nevertheless, unlike other proteins encoded by the same operon, the structural and physiological aspects of YfdX have been poorly characterized. Here, we identified a previously unknown pH-dependent stoichiometric conversion of S. Typhi YfdX between dimeric and tetrameric states; this conversion was further analyzed via determining its structure by X-ray crystallography at high resolution and by small-angle X-ray scattering in a solution state and via structure-based mutant studies. Biologically, YfdX was proven to be critically involved in Salmonella susceptibility to two β-lactam antibiotics, penicillin G and carbenicillin, as bacterial growth significantly impaired by its deficiency upon treatment with each of the two antibiotics was recovered by chromosomal complementation. Furthermore, by using Galleria mellonella larvae as an in vivo model of Salmonella infection, we demonstrated that Salmonella virulence was remarkably enhanced by YfdX deficiency, which was complemented by a transient expression of the wild-type or dimeric mutant but not by that of the monomeric mutant. The present study work provides direct evidence regarding the participation of YfdX in Salmonella antibiotic susceptibility and in the modulation of bacterial virulence, providing a new insight into this pathogen's strategies for survival and growth.