Positive cooperative regulation of double binding sites for human acetylcholinesterase

Does acute exposure to thimerosal, an organic mercury compound, affect the mitochondrial function of an infant model?

BACKGROUND

Thimerosal (TM) is a toxic, organometallic mercury compound (which releases ethyl-mercury-containing compounds in aqueous solutions) used as a preservative in vaccines. Mitochondria are organelle which are highly vulnerable to many chemical compounds, including mercury (Hg) and its derivatives.

METHOD

Wistar rats (at 21 days of age) were used to model a child's TM exposure following childhood vaccination, divided in two groups: TM exposed (20 μg/kg/day) and unexposed controls (saline solution), both for 24 h. Atomic Fluorescence Spectrometry was used to quantify the amounts of mercury in tissues. The electron transport chain (ETC) from isolated mitochondria was evaluated using an oxygen electrode. The mitochondrial membrane potential and H2O2 production were analyzed using selective fluorescence probes. The activity of some enzymes (SOD, CAT, GPx, and AChE) and secondary markers of oxidative stress (GSH, GSSG, total free thiol) were also examined in tissues.

RESULTS

Hg accumulation in the brain and liver was higher in exposed animals when compared to the control. Liver-isolated mitochondria showed that TM improved respiratory control by 23%; however, states 3 and 4 of the ETC presented a decrease of 16% and 37%, respectively. Furthermore, brain-isolated mitochondria presented an improvement of 61% in respiratory control. Brain enzyme activities were significantly impacted in TM-exposed rats compared to unexposed rats as follows: decreases in SOD (32%) and AChE (42%) and increases in GPx (79%) and CAT (100%). GPx enzyme activity in the liver was significantly increased (37%). Among secondary oxidative stress markers, the brain's total reduced thiol (SH) concentration was significantly increased (41%).

CONCLUSION

Acute TM treatment exposure in a Wistar rat model mimicking TM exposure in an infant following childhood vaccination significantly damaged brain bioenergetic pathways. This study supports the ability of TM exposure to preferentially damage the nervous system.

Creatinase: Using Increased Entropy to Improve the Activity and Thermostability

Improving protein thermostability in mutagenesis-based enzyme engineering was often achieved by enhancing interresidue interactions via mutation to increase the enthalpy penalty of unfolding. However, this approach may trade off the functional activity due to the loss of structural flexibility of the biomolecule. Here, by performing X-ray crystallography, enzymatic kinetic experiments, neutron scattering, and thermodynamical measurements, we compared the structures, catalytic behaviors, dynamics, and thermostability between a wild-type creatinase and its four-point mutant. We found that the mutant is an entropy-driven thermostable protein with higher structural flexibility, i.e., higher conformational entropy, in the folded state compared to the wild type. The increased conformational entropy of the mutant in the folded state can reduce the entropy gain during unfolding and thus renders it greater thermostability. Moreover, the increased structural flexibility, particularly around the catalytic site, can broaden the mutant's working temperature range and considerably improve its activity at ambient conditions, which is crucial for its application in diagnosing kidney diseases. Complementary all-atom molecular dynamics simulations indicated that the four mutations replaced several of the strong interresidue interactions (electrostatic interactions and hydrogen bonds) with weak hydrophobic interactions. These substitutions not only release the structural flexibility to promote the thermostability and enzymatic activity of the protein but they also preserve the protein structure from collapsing. Our findings may pave a route for the entropy-driven strategy to design proteins with high thermostability and activity.

Identification of hACE2-interacting sites in SARS-CoV-2 spike receptor binding domain for antiviral drugs screening

The key structure of the interface between the spike protein of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) and human angiotensin-converting enzyme 2 (hACE2) acts as an essential switch for cell entry by the virus and drugs targets. However, this is largely unknown. Here, we tested three peptides of spike receptor binding domain (RBD) and found that peptide 391-465 aa is the major hACE2-interacting sites in SARS-CoV-2 spike RBD. We then identified essential amino acid residues (403R, 449Y, 454R) of peptide 391-465 aa that were critical for the interaction between the RBD and hACE2. Additionally, a pseudotyped virus containing SARS-CoV-2 spike with individual mutation (R454G, Y449F, R403G, N439I, or N440I) was determined to have very low infectivity compared with the pseudotyped virus containing the wildtype (WT) spike from reference strain Wuhan 1, respectively. Furthermore, we showed the key amino acids had the potential to drug screening. For example, molecular docking (Docking) and infection assay showed that Cephalosporin derivatives can bind with the key amino acids to efficiently block infection of the pseudoviruses with wild type spike or new variants. Moreover, Cefixime inhibited live SARS-CoV-2 infection. These results also provide a novel model for drug screening and support further clinical evaluation and development of Cephalosporin derivatives as novel, safe, and cost-effective drugs for prevention/treatment of SARS-CoV-2.

Catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase

Methane is among the most potent of the greenhouse gases, which plays a key role in global climate change. As an excellent carbon and energy source, methane can be utilized by anaerobic methane oxidizing archaea and aerobic methane oxidizing bacteria. The previous work shows that an anaerobic thermophilic enrichment culture composed of dense consortia of archaea and bacteria apparently uses partly similar pathways to oxidize the C4 hydrocarbon butane. However, the catalytic mechanism of butane anaerobic oxidation for alkyl-coenzyme M reductase is still unknown. Therefore, molecular dynamics (MD) simulation was used to investigate the dynamics differences of catalytic mechanism between methane coenzyme M reductase (MCR) and alkyl-coenzyme M reductase (ACR). At first, the binding pocket of ACR is larger than that of MCR. Then, the complex of butane and ACR is more stable than that of methane and ACR. Protein conformation cloud suggests that the position of methane is dynamics and methane escapes from the binding pocket of ACR during most of the simulation time, while butane tightly binds in the pocket of ACR. The hydrophobic interactions between butane and ACR are more and stronger than those between methane and ACR. At the same time, the binding free energy between butane and ACR is significantly lower than that between methane and ACR. The dynamics correlation network indicates that the transformation of information flow for ACR-butane is smoother than that for ACR-methane. The shortest pathway for ACR-butane is from Gln144, Ala141, Hie135, Ile133, Ala160, Arg206, Asp97, Met94, Tyr347 to Phe345 with synergistic effect for two butane molecules. This study can insight into the catalytic mechanism for butane/ACR complex.

Dynamical important residue network (DIRN): network inference via conformational change

MOTIVATION

Protein residue interaction network has emerged as a useful strategy to understand the complex relationship between protein structures and functions and how functions are regulated. In a residue interaction network, every residue is used to define a network node, adding noises in network post-analysis and increasing computational burden. In addition, dynamical information is often necessary in deciphering biological functions.

RESULTS

We developed a robust and efficient protein residue interaction network method, termed dynamical important residue network, by combining both structural and dynamical information. A major departure from previous approaches is our attempt to identify important residues most important for functional regulation before a network is constructed, leading to a much simpler network with the important residues as its nodes. The important residues are identified by monitoring structural data from ensemble molecular dynamics simulations of proteins in different functional states. Our tests show that the new method performs well with overall higher sensitivity than existing approaches in identifying important residues and interactions in tested proteins, so it can be used in studies of protein functions to provide useful hypotheses in identifying key residues and interactions.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Secondary structures transition of tau protein with intrinsically disordered proteins specific force field

Microtubule-associated Tau protein plays a key role in assembling microtubule and modulating the functional organization of the neuron and developing axonal morphology, growth, and polarity. The pathological Tau can aggregate into cross-beta amyloid as one of the hallmarks for Alzheimer's disease (AD). Therefore, one of the top priorities in AD research is to figure out the structural model of Tau aggregation and to screen the inhibitors. The latest generation intrinsically disordered protein specific force field ff14IDPSFF significantly improved the distributions of heterogeneous conformations for intrinsically disordered proteins (IDPs). Here, the molecular dynamics (MD) simulations with three force fields of ff14SB, ff14IDPs, and ff14IDPSFF were employed to investigate the secondary structures transition of Tau (267-312) fragment. The results indicate that ff14IDPSFF can generate more heterogeneous conformers, and the predicted secondary structural distribution is closer to that of the experimental observation. In addition, predicted secondary chemical shifts from ff14IDPSFF are the most approach to those of experiment. Secondary structures transition kinetics for Tau(267-312) with ff14IDPSFF shows that the secondary structures were gradually transformed from α-helix to β-strand and the β-strand located at the regions of the residues 274-280 and residues 305-311. Besides, the driving force for the secondary structures transition of Tau(267-312) is mainly hydrophobic interactions which located at hexa-peptides ²⁷⁵ VQIINK²⁸⁰ and ³⁰⁶ VQIVYK³¹¹ . Secondary structure transition of Tau protein can give insight into the aggregation mechanism for AD.

Order-disorder transition of intrinsically disordered kinase inducible transactivation domain of CREB

Transcription factor cyclic Adenosine monophosphate response-element binding protein plays a critical role in the cyclic AMP response pathway via its intrinsically disordered kinase inducible transactivation domain (KID). KID is one of the most studied intrinsically disordered proteins (IDPs), although most previous studies focus on characterizing its disordered state structures. An interesting question that remains to be answered is how the order-disorder transition occurs at experimental conditions. Thanks to the newly developed IDP-specific force field ff14IDPSFF, the quality of conformer sampling for IDPs has been dramatically improved. In this study, molecular dynamics (MD) simulations were used to study the order-to-disorder transition kinetics of KID based on the good agreement with the experiment on its disordered-state properties. Specifically, we tested four force fields, ff99SBildn, ff99IDPs, ff14IDPSFF, and ff14IDPs in the simulations of KID and found that ff14IDPSFF can generate more diversified disordered conformers and also reproduce more accurate experimental secondary chemical shifts. Kinetics analysis of MD simulations demonstrates that the order-disorder transition of KID obeys the first-order kinetics, and the transition nucleus is I127/L128/L141. The possible transition pathways from the nucleus to the last folded residues were identified as I127-R125-L138-L141-S143-A145 and L128-R125-L138-L141-S143-A145 based on a residue-level dynamical network analysis. These computational studies not only provide testable prediction/hypothesis on the order-disorder transition of KID but also confirm that the ff14IDPSFF force field can be used to explore the correlation between the structure and function of IDPs.

Synergistic regulation mechanism of iperoxo and LY2119620 for muscarinic acetylcholine M2 receptor

Muscarinic acetylcholine receptors are GPCRs that regulate the activity of a diverse array of central and peripheral functions in the human body, including the parasympathetic actions of acetylcholine. The M2 muscarinic receptor subtype plays a key role in modulating cardiac function and many important central processes. The orthosteric agonist and allosteric modulator can bind the pocket of M2. However, the detailed relationship between orthosteric agonist and allosteric modulator of M2 is still unclear. In this study, we intend to elucidate the residue-level regulation mechanism and pathway via a combined approach of dynamical correlation network and molecular dynamics simulation. Specifically computational residue-level fluctuation correlation data was analyzed to reveal detailed dynamics signatures in the regulation process. A hypothesis of "synergistic regulation" is proposed to reveal the cooperation affection between the orthosteric agonist and allosteric modulator, which is subsequently validated by perturbation and mutation analyses. Two possible synergistic regulation pathways of 2CU-I178-Y403-W400-F396-L114-Y440-Nb9 and IXO-V111-F396-L114-Y440-Nb9 were identified by the shortest path algorithm and were confirmed by the mutation of junction node. Furthermore, the efficiency of information transfer of bound M2 is significant higher than any single binding system. Our study shows that targeting the synergistic regulation pathways may better regulate the calcium channel of M2. The knowledge gained in this study may help develop drugs for diseases of the central nervous system and metabolic disorders.