Conformational stability and order of Hoogsteen base pair induced by protein binding

Several experimental studies have shown that Hoogsteen (HG) base pair (bp) stabilizes in the presence of proteins. The molecular mechanism underlying this stabilization is not well known. This leads us to examine the stability of the HG bp in duplex DNA using all-atom molecular dynamics simulation in both the absence and presence of proteins. We use conformational thermodynamics to investigate the stability of a HG bp in duplex DNA at the molecular level. We compute the changes in the conformational free energy and entropy of DNA when DNA adopts a HG bp in its bp sequence rather than a Watson-Crick (WC) bp in both naked DNA and protein-bound DNA complex. We observe that the presence of proteins stabilizes and organizes the HG bp and the entire DNA duplex. Sugar-phosphate, sugar-base, and sugar-pucker torsion angles play key roles in stabilizing and ordering the HG bp in the protein-bound DNA complex.

Computational Studies of Cardiac and Skeletal Troponin

Troponin is a key regulatory protein in muscle contraction, consisting of three subunits troponin C (TnC), troponin I (TnI), and troponin T (TnT). Calcium association to TnC initiates contraction by causing a series of dynamic and conformational changes that allow the switch peptide of TnI to bind and subsequently cross bridges to form between the thin and thick filament of the sarcomere. Owing to its pivotal role in contraction regulation, troponin has been the focus of numerous computational studies over the last decade. These studies elegantly supplemented a large volume of experimental work and focused on the structure, dynamics and function of the whole troponin complex, individual subunits, and even on segments of the thin filament. Molecular dynamics, Brownian dynamics, and free energy simulations have been used to elucidate the conformational dynamics and underlying free energy landscape of troponin, calcium, and switch peptide binding, as well as the effect of disease mutations, small molecules and post-translational modifications such as phosphorylation. Frequently, simulations have been used to confirm or explain experimental observations. Computer-aided drug discovery tools have been employed to identify novel potential calcium sensitizing agents binding to the TnC-TnI interface. Finally, Markov modeling has contributed to simulating contraction within the sarcomere on the mesoscale. Here we are reviewing and classifying the existing computational work on troponin and its subunits, outline current gaps in simulations elucidating troponin's role in contraction and suggest potential future developments in the field.

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.

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

We perform molecular dynamics simulation studies on interaction between bacterial proteins: an outer-membrane protein STY3179 and a yfdX protein STY3178 of Salmonella Typhi. STY3179 has been found to be involved in bacterial adhesion and invasion. STY3178 is recently biophysically characterized. It is a soluble protein having antibiotic binding and chaperon activity capabilities. These two proteins co-occur and are from neighboring gene in Salmonella Typhi-occurrence of homologs of both STY3178 and STY3179 are identified in many Gram-negative bacteria. We show using homology modeling, docking followed by molecular dynamics simulation that they can form a stable complex. STY3178 belongs to aqueous phase, while the beta barrel portion of STY3179 remains buried in DPPC bilayer with extra-cellular loops exposed to water. To understand the molecular basis of interaction between STY3178 and STY3179, we compute the conformational thermodynamics which indicate that these two proteins interact through polar and acidic residues belonging to their interfacial region. Conformational thermodynamics results further reveal instability of certain residues in extra-cellular loops of STY3179 upon complexation with STY3178 which is an indication for binding with host cell protein laminin.

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.

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.