scholarly journals Effects of Cardiolipin on the Conformational Dynamics of Membrane-Anchored Bcl-xL

2021 ◽  
Vol 22 (17) ◽  
pp. 9388
Author(s):  
Vivek Tyagi ◽  
Victor Vasquez-Montes ◽  
J. Alfredo Freites ◽  
Alexander Kyrychenko ◽  
Douglas J. Tobias ◽  
...  

The anti-apoptotic protein Bcl-xL regulates apoptosis by preventing the permeation of the mitochondrial outer membrane by pro-apoptotic pore-forming proteins, which release apoptotic factors into the cytosol that ultimately lead to cell death. Two different membrane-integrated Bcl-xL constructs have been identified: a membrane-anchored and a membrane-inserted conformation. Here, we use molecular dynamics simulations to study the effect of the mitochondrial specific lipid cardiolipin and the protein protonation state on the conformational dynamics of membrane-anchored Bcl-xL. The analysis reveals that the protonation state of the protein and cardiolipin content of the membrane modulate the orientation of the soluble head region (helices α1 through α7) and hence the exposure of its BH3-binding groove, which is required for its interaction with pro-apoptotic proteins.

Author(s):  
Balaji Selvam ◽  
Ya-Chi Yu ◽  
Liqing Chen ◽  
Diwakar Shukla

<p>The SWEET family belongs to a class of transporters in plants that undergoes large conformational changes to facilitate transport of sugar molecules across the cell membrane. However, the structures of their functionally relevant conformational states in the transport cycle have not been reported. In this study, we have characterized the conformational dynamics and complete transport cycle of glucose in OsSWEET2b transporter using extensive molecular dynamics simulations. Using Markov state models, we estimated the free energy barrier associated with different states as well as 1 for the glucose the transport mechanism. SWEETs undergoes structural transition to outward-facing (OF), Occluded (OC) and inward-facing (IF) and strongly support alternate access transport mechanism. The glucose diffuses freely from outside to inside the cell without causing major conformational changes which means that the conformations of glucose unbound and bound snapshots are exactly same for OF, OC and IF states. We identified a network of hydrophobic core residues at the center of the transporter that restricts the glucose entry to the cytoplasmic side and act as an intracellular hydrophobic gate. The mechanistic predictions from molecular dynamics simulations are validated using site-directed mutagenesis experiments. Our simulation also revealed hourglass like intermediate states making the pore radius narrower at the center. This work provides new fundamental insights into how substrate-transporter interactions actively change the free energy landscape of the transport cycle to facilitate enhanced transport activity.</p>


2020 ◽  
Vol 16 (4) ◽  
pp. 451-459 ◽  
Author(s):  
Fortunatus C. Ezebuo ◽  
Ikemefuna C. Uzochukwu

Background: Sulfotransferase family comprises key enzymes involved in drug metabolism. Oxamniquine is a pro-drug converted into its active form by schistosomal sulfotransferase. The conformational dynamics of side-chain amino acid residues at the binding site of schistosomal sulfotransferase towards activation of oxamniquine has not received attention. Objective: The study investigated the conformational dynamics of binding site residues in free and oxamniquine bound schistosomal sulfotransferase systems and their contribution to the mechanism of oxamniquine activation by schistosomal sulfotransferase using molecular dynamics simulations and binding energy calculations. Methods: Schistosomal sulfotransferase was obtained from Protein Data Bank and both the free and oxamniquine bound forms were subjected to molecular dynamics simulations using GROMACS-4.5.5 after modeling it’s missing amino acid residues with SWISS-MODEL. Amino acid residues at its binding site for oxamniquine was determined and used for Principal Component Analysis and calculations of side-chain dihedrals. In addition, binding energy of the oxamniquine bound system was calculated using g_MMPBSA. Results: The results showed that binding site amino acid residues in free and oxamniquine bound sulfotransferase sampled different conformational space involving several rotameric states. Importantly, Phe45, Ile145 and Leu241 generated newly induced conformations, whereas Phe41 exhibited shift in equilibrium of its conformational distribution. In addition, the result showed binding energy of -130.091 ± 8.800 KJ/mol and Phe45 contributed -9.8576 KJ/mol. Conclusion: The results showed that schistosomal sulfotransferase binds oxamniquine by relying on hybrid mechanism of induced fit and conformational selection models. The findings offer new insight into sulfotransferase engineering and design of new drugs that target sulfotransferase.


2021 ◽  
Vol 22 (9) ◽  
pp. 4769
Author(s):  
Pablo Maturana ◽  
María S. Orellana ◽  
Sixto M. Herrera ◽  
Ignacio Martínez ◽  
Maximiliano Figueroa ◽  
...  

Agmatine is the product of the decarboxylation of L-arginine by the enzyme arginine decarboxylase. This amine has been attributed to neurotransmitter functions, anticonvulsant, anti-neurotoxic, and antidepressant in mammals and is a potential therapeutic agent for diseases such as Alzheimer’s, Parkinson’s, and cancer. Agmatinase enzyme hydrolyze agmatine into urea and putrescine, which belong to one of the pathways producing polyamines, essential for cell proliferation. Agmatinase from Escherichia coli (EcAGM) has been widely studied and kinetically characterized, described as highly specific for agmatine. In this study, we analyze the amino acids involved in the high specificity of EcAGM, performing a series of mutations in two loops critical to the active-site entrance. Two structures in different space groups were solved by X-ray crystallography, one at low resolution (3.2 Å), including a guanidine group; and other at high resolution (1.8 Å) which presents urea and agmatine in the active site. These structures made it possible to understand the interface interactions between subunits that allow the hexameric state and postulate a catalytic mechanism according to the Mn2+ and urea/guanidine binding site. Molecular dynamics simulations evaluated the conformational dynamics of EcAGM and residues participating in non-binding interactions. Simulations showed the high dynamics of loops of the active site entrance and evidenced the relevance of Trp68, located in the adjacent subunit, to stabilize the amino group of agmatine by cation-pi interaction. These results allow to have a structural view of the best-kinetic characterized agmatinase in literature up to now.


Glycobiology ◽  
2021 ◽  
Author(s):  
Margrethe Gaardløs ◽  
Sergey A Samsonov ◽  
Marit Sletmoen ◽  
Maya Hjørnevik ◽  
Gerd Inger Sætrom ◽  
...  

Abstract Mannuronan C-5 epimerases catalyse the epimerization of monomer residues in the polysaccharide alginate, changing the physical properties of the biopolymer. The enzymes are utilized to tailor alginate to numerous biological functions by alginate-producing organisms. The underlying molecular mechanisms that control the processive movement of the epimerase along the substrate chain is still elusive. To study this, we have used an interdisciplinary approach combining molecular dynamics simulations with experimental methods from mutant studies of AlgE4, where initial epimerase activity and product formation were addressed with NMR spectroscopy, and characteristics of enzyme-substrate interactions were obtained with isothermal titration calorimetry and optical tweezers. Positive charges lining the substrate-binding groove of AlgE4 appear to control the initial binding of poly-mannuronate, and binding also seems to be mediated by both electrostatic and hydrophobic interactions. After the catalytic reaction, negatively charged enzyme residues might facilitate dissociation of alginate from the positive residues, working like electrostatic switches, allowing the substrate to translocate in the binding groove. Molecular simulations show translocation increments of two monosaccharide units before the next productive binding event resulting in MG-block formation, with the epimerase moving with its N-terminus towards the reducing end of the alginate chain. Our results indicate that the charge pair R343-D345 might be directly involved in conformational changes of a loop that can be important for binding and dissociation. The computational and experimental approaches used in this study complement each other, allowing for a better understanding of individual residues’ roles in binding and movement along the alginate chains.


Processes ◽  
2021 ◽  
Vol 9 (7) ◽  
pp. 1085
Author(s):  
Paula Mihaljević-Jurič ◽  
Sérgio F. Sousa

Statins are important drugs in the regulation of cholesterol levels in the human body that have as a primary target the enzyme β-hydroxy-β-methylglutaryl-CoA reductase (HMGR). This enzyme plays a crucial role in the mevalonate pathway, catalyzing the four-electron reduction of HMG-CoA to mevalonate. A second reduction step of this reaction mechanism has been the subject of much speculation in the literature, with different conflicting theories persisting to the present day. In this study, the different mechanistic hypotheses were evaluated with atomic-level detail through a combination of molecular dynamics simulations (MD) and quantum mechanics/molecular mechanics (QM/MM) calculations. The obtained Gibbs free activation and Gibbs free reaction energy (15 kcal mol−1 and −40 kcal mol−1) show that this hydride step takes place with the involvement of a cationic His405 and Lys639, and a neutral Glu98, while Asp715 remains in an anionic state. The results provide an atomic-level portrait of this step, clearly demonstrating the nature and protonation state of the amino acid residues involved, the energetics associated, and the structure and charge of the key participating atoms in the several intermediate states, finally elucidating this missing step.


2021 ◽  
Vol 12 (1) ◽  
Author(s):  
Carlos G. Acevedo-Rocha ◽  
Aitao Li ◽  
Lorenzo D’Amore ◽  
Sabrina Hoebenreich ◽  
Joaquin Sanchis ◽  
...  

AbstractMultidimensional fitness landscapes provide insights into the molecular basis of laboratory and natural evolution. To date, such efforts usually focus on limited protein families and a single enzyme trait, with little concern about the relationship between protein epistasis and conformational dynamics. Here, we report a multiparametric fitness landscape for a cytochrome P450 monooxygenase that was engineered for the regio- and stereoselective hydroxylation of a steroid. We develop a computational program to automatically quantify non-additive effects among all possible mutational pathways, finding pervasive cooperative signs and magnitude epistasis on multiple catalytic traits. By using quantum mechanics and molecular dynamics simulations, we show that these effects are modulated by long-range interactions in loops, helices and β-strands that gate the substrate access channel allowing for optimal catalysis. Our work highlights the importance of conformational dynamics on epistasis in an enzyme involved in secondary metabolism and offers insights for engineering P450s.


2017 ◽  
Vol 114 (42) ◽  
pp. E8830-E8836 ◽  
Author(s):  
Chang Yun Son ◽  
Arun Yethiraj ◽  
Qiang Cui

Cytochrome c oxidase (CcO) is a transmembrane protein that uses the free energy of O2 reduction to generate the proton concentration gradient across the membrane. The regulation of competitive proton transfer pathways has been established to be essential to the vectorial transport efficiency of CcO, yet the underlying mechanism at the molecular level remains lacking. Recent studies have highlighted the potential importance of hydration-level change in an internal cavity that connects the proton entrance channel, the site of O2 reduction, and the putative proton exit route. In this work, we use atomistic molecular dynamics simulations to investigate the energetics and timescales associated with the volume fluctuation and hydration-level change in this central cavity. Extensive unrestrained molecular dynamics simulations (accumulatively ∼4 μs) and free energy computations for different chemical states of CcO support a model in which the volume and hydration level of the cavity are regulated by the protonation state of a propionate group of heme a3 and, to a lesser degree, the redox state of heme a and protonation state of Glu286. Markov-state model analysis of ∼2-μs trajectories suggests that hydration-level change occurs on the timescale of 100–200 ns before the proton-loading site is protonated. The computed energetic and kinetic features for the cavity wetting transition suggest that reversible hydration-level change of the cavity can indeed be a key factor that regulates the branching of proton transfer events and therefore contributes to the vectorial efficiency of proton transport.


2020 ◽  
Vol 10 (1) ◽  
pp. 4944-4955 ◽  

Alzheimer’s disease (AD) is characterized by the presence of Amyloid-beta (Aβ) peptide, which has the propensity to fold into β-sheets under stress forming aggregated amyloid plaques. Nowadays many studies have focused on the development of novel, specific therapeutic strategies to slow down Aβ aggregation or control preformed aggregates. Albumin, the most abundant protein in the cerebrospinal fluid, was reported to bind Aβ impeding its aggregation. Recently, it has been reported that C-terminal (CTerm) of Human Albumin binds with Aβ1-42, impairs Aβ aggregation and promotes disassembly of Aβ aggregates protecting neurons. In this computational study, we have investigated the effect of CTerm on the conformational dynamics and the aggregation propensity of Aβ1-42 peptide. We have performed molecular dynamics simulations on the Aβ1-42-Aβ1-42 homodimer and Aβ1-42-CTerm of albumin heterodimer using the AMBER force field ff99SBildn. From the Potential of mean force (PMF) study and Binding free energy (BFE) analysis, we observed the association of Aβ1-42 peptide monomer with itself in the form of homodimer to be stronger than its association with the CTerm in the heterodimer complex. The difference in the number of residues in the Aβ1-42 peptide monomer (42 AAs) and CTerm (35 AAs) may be probable reason for the difference in association between the monomeric units in corresponding homodimer and heterodimer complexes. But even then CTerm shows a significant effect on the dimerization of Aβ1-42 peptide. Our findings therefore suggest that CTerm can be used for the disassembly of Aβ1-42 peptide monomer.


eLife ◽  
2015 ◽  
Vol 4 ◽  
Author(s):  
C Keith Cassidy ◽  
Benjamin A Himes ◽  
Frances J Alvarez ◽  
Jun Ma ◽  
Gongpu Zhao ◽  
...  

Chemotactic responses in bacteria require large, highly ordered arrays of sensory proteins to mediate the signal transduction that ultimately controls cell motility. A mechanistic understanding of the molecular events underlying signaling, however, has been hampered by the lack of a high-resolution structural description of the extended array. Here, we report a novel reconstitution of the array, involving the receptor signaling domain, histidine kinase CheA, and adaptor protein CheW, as well as a density map of the core-signaling unit at 11.3 Å resolution, obtained by cryo-electron tomography and sub-tomogram averaging. Extracting key structural constraints from our density map, we computationally construct and refine an atomic model of the core array structure, exposing novel interfaces between the component proteins. Using all-atom molecular dynamics simulations, we further reveal a distinctive conformational change in CheA. Mutagenesis and chemical cross-linking experiments confirm the importance of the conformational dynamics of CheA for chemotactic function.


Sign in / Sign up

Export Citation Format

Share Document