Crystal structure of the flagellar accessory protein FlaH of Methanocaldococcus jannaschii suggests a regulatory role in archaeal flagellum assembly

Archaeal flagella are unique structures that share functional similarity with bacterial flagella, but are structurally related to bacterial type IV pili. The flagellar accessory protein FlaH is one of the conserved components of the archaeal motility system. However, its function is not clearly understood. Here, we present the 2.2 Å resolution crystal structure of FlaH from the hyperthermophilic archaeon, Methanocaldococcus jannaschii. The protein has a characteristic RecA‐like fold, which has been found previously both in archaea and bacteria. We show that FlaH binds to immobilized ATP—however, it lacks ATPase activity. Surface plasmon resonance analysis demonstrates that ATP affects the interaction between FlaH and the archaeal motor protein FlaI. In the presence of ATP, the FlaH‐FlaI interaction becomes significantly weaker. A database search revealed similarity between FlaH and several DNA‐binding proteins of the RecA superfamily. The closest structural homologs of FlaH are KaiC‐like proteins, which are archaeal homologs of the circadian clock protein KaiC from cyanobacteria. We propose that one of the functions of FlaH may be the regulation of archaeal motor complex assembly.     

Roles of conserved basic amino acid residues and activation mechanism of the hyperthermophilic aspartate racemase at high temperature

X-ray crystallography has revealed two similar alpha/beta domains of the aspartate racemase from the hyperthermophilic archaeon, Pyrococcus horikoshii OT3. The active site is located in the cleft between the two domains where two cysteine residues face each other. This arrangement allows the substrate to enter the cleft and enables the two cysteine residues to act synergistically. However, the distance between their thiolates was estimated to be 9.6 angstroms, which is beyond the distance for cooperative action of them. We examined the molecular mechanism for the racemization reaction of this hyperthermophilic aspartate racemase by mutational analyses and molecular dynamics simulations. The mutational analyses revealed that Arg48 and Lys164 were essential for catalysis in addition to the putative catalytic cysteine residues. The molecular dynamics simulations revealed that the distance between the two active gamma-sulfur atoms of cysteine residues oscillate to periodically become shorter than the predicted cooperative distance at high temperature. In addition, the conformation of Tyr160, which is located at the entrance of the cleft and inhibits the entry of a substrate, changes periodically to open the entrance at 375 K. The opening of the gate is likely to be induced by the motion of the adjacent amino acid, Lys164. The entrance of an aspartate molecule was observed by molecular dynamics (MD) simulations driven by the force of the electrostatic interaction with Arg48, Lys164, and also Asp47. These results provide insights into the roles of amino acid residues at the catalytic site and also the activation mechanism of a hyperthermophilic aspartate racemase at high temperature.     

Inhibition of cysteine proteases by a natural biflavone: behavioral evaluation of fukugetin as papain and cruzain inhibitor

Cruzain is the major cysteine protease of Trypanosoma cruzi, the infectious agent responsible for Chagas disease, and cruzain inhibitors display considerable antitrypanosomal activity. In the present work we elucidated crystallographic data of fukugetin, a biflavone isolated from Garcinia brasiliensis, and investigated the role of this molecule as cysteine protease inhibitor. The kinetic analyses demonstrated that fukugetin inhibited cruzain and papain by a slow reversible type inhibition with K_I of 1.1 and 13.4 µM, respectively. However, cruzain inhibition was about 12 times faster than papain inhibition. Lineweaver–Burk plots demonstrated partial competitive inhibition for cruzain and hyperbolic mixed-type inhibition for papain. Furthermore, the docking results showed that the biflavone binds to ring C′ in the S2 pocket and to ring C in the S3 pocket through hydrophobic interactions and hydrogen bonds. Finally, fukugetin also presented inhibitory activity on proteases of the T. cruzi extract, with IC₅₀ of 7 µM.     

Synthesis of the sulfur amino acids: cysteine and methionine

This review will assess new features reported for the molecular and biochemical aspects of cysteine and methionine biosynthesis in Arabidopsis thaliana with regards to early published data from other taxa including crop plants and bacteria (Escherichia coli as a model). By contrast to bacteria and fungi, plant cells present a complex organization, in which the sulfur network takes place in multiple sites. Particularly, the impact of sulfur amino-acid biosynthesis compartmentalization will be addressed in respect to localization of sulfur reduction. To this end, the review will focus on regulation of sulfate reduction by synthesis of cysteine through the cysteine synthase complex and the synthesis of methionine and its derivatives. Finally, regulatory aspects of sulfur amino-acid biosynthesis will be explored with regards to interlacing processes such as photosynthesis, carbon and nitrogen assimilation.     

Electrostatics of cysteine residues in proteins: Parameterization and validation of a simple model

One of the most popular and simple models for the calculation of pK_as from a protein structure is the semi‐macroscopic electrostatic model MEAD. This model requires empirical parameters for each residue to calculate pK_as. Analysis of current, widely used empirical parameters for cysteine residues showed that they did not reproduce expected cysteine pK_as; thus, we set out to identify parameters consistent with the CHARMM27 force field that capture both the behavior of typical cysteines in proteins and the behavior of cysteines which have perturbed pK_as. The new parameters were validated in three ways: (1) calculation across a large set of typical cysteines in proteins (where the calculations are expected to reproduce expected ensemble behavior); (2) calculation across a set of perturbed cysteines in proteins (where the calculations are expected to reproduce the shifted ensemble behavior); and (3) comparison to experimentally determined pK_a values (where the calculation should reproduce the pK_a within experimental error). Both the general behavior of cysteines in proteins and the perturbed pK_a in some proteins can be predicted reasonably well using the newly determined empirical parameters within the MEAD model for protein electrostatics. This study provides the first general analysis of the electrostatics of cysteines in proteins, with specific attention paid to capturing both the behavior of typical cysteines in a protein and the behavior of cysteines whose pK_a should be shifted, and validation of force field parameters for cysteine residues. Proteins 2012. © 2012 Wiley Periodicals, Inc.     

Molecular dynamics simulation of the substitution of serine for the conserved glycine in the G-loop in the cdc28-srm yeast mutant using the crystal lattice of human CDK2 kinase

Two-nanosecond molecular dynamics simulations of the crystal lattice of an active complex of pT160-CDK2 kinase/cyclin A/ATP-Mg²⁺/substrate were performed. The simulations showed that the structures of the wild-type CDK2 complex and the mutant CDK2 complex involving the substitution G16S-CDK2 corresponding to the yeast substitution G20S-CDC28 differ noticeably and the differences between the structural conformations are most pronounced in the regions that play a key role in the kinase functioning. The results of the computer calculations were used to consider the structural elements that may affect the kinase activity, the regulatory phosphorylation, and the binding of protein kinase with cyclins and substrates.     

Identification of a single and non-essential cysteine residue in dextransucrase of Leuconostoc mesenteroides NRRL B-512F

Amino acid analysis of purified dextransucrase (sucrose: 1,6-α-D-glucan 6-α-D-glucosyltransferase EC 2.4.1.5) from Leuconostoc mesenteroides NRRL B-512F was carried out. The enzyme is virtually devoid of cysteine residue there being only one cysteine residue in the whole enzyme molecule comprising over 1500 amino acid residues. The enzyme is rich in acidic amino acid residues. The number of amino acid residues was calculated based on the molecular weight of 188,000 (Goyal and Katiyar 1994). Amino sugars were not found, implying that the enzyme is not a glycoprotein. It has been shown earlier that the cysteine residue in dextransucrase is not essential for enzyme activity (Goyal and Katiyar 1998). The presence of only one cysteine residue per enzyme molecule illustrates that its tertiary structure is solely dependent on other types of non-covalent interactions such as hydrogen bonding, ionic and nonpolar hydrophobic interactions.     

Simultaneous determination of cysteine sulfinic acid and cysteic acid in rat brain by high-performance liquid chromatography

A sensitive and specific assay method for cysteine sulfinic acid (CSA) and cysteic acid (CA) using high-performance liquid chromatography has been developed. The method includes post-column derivatization of various amino acids with o-phthalaldehyde in the presence of 2-mercaptoethanol. The column packed with cation-exchange resin ($\text{ISC-07}\text{S1504}$, Shimadzu Sci entific instruments, Inc., Kyoto, Japan) was used for obtaining general separation of amino acids except CSA and CA, while the separation of CSA and CA was achieved using a strong-base anion exchange ($\text{ISA-07}\text{S2504}$, Shimadzu Scientific Instruments) column. The fluorescence peak area for CSA was linear between 20 pmol and 5 nmol, whereas that for CA was 10 pmol to 5 nmol. The regional distribution of CSA, CA, and other amino acids in the rat brain was studied using this new assay method.     

The role of a conserved histidine-tyrosine interhelical interaction in the ion channel domain of delta-endotoxins from Bacillus thuringiensis

One of the molecular factors contributing to Leishmania sp. virulence and pathogenesis is the major surface metalloprotease GP63, alternatively called leishmanolysin, MSP, and PSP (EC 3.4.24.36). Here, the molecular dynamics simulation of Leishmania major GP63 in water at pH 7 is reported. Upon solvation, GP63 undergoes a sharp structural relaxation with respect to the crystal structure. The fluctuation pattern occurs essentially in solvent-exposed nonstructured regions. By contrast, the active site turns out to be rigid. Essential dynamics and dynamic-domain analyses, both carried out on the equilibrated portion of GP63, show that the fingerprint fluctuations of GP63 are practically characterized by the motion of a large part of the N-terminal domain. These results appear to be in line with substrate recognition and (pro)enzyme activation played by the N-terminal domain of GP63. A systematic analysis among a series of 10 homologs of GP63 also shows that the residues involved in the interdomain bending result highly conserved. This finding also suggests possible relationship between the maintainance of proteolytic activity and the similarity of the dynamical properties of the related enzymes.     

The concerted action of a positive charge and hydrogen bonds dynamically regulates the pKa of the nucleophilic cysteine in the NrdH‐redoxin family

NrdH‐redoxins shuffle electrons from the NADPH pool in the cell to Class Ib ribonucleotide reductases, which in turn provide the precursors for DNA replication and repair. NrdH‐redoxins have a CVQC active site motif and belong to the thioredoxin‐fold protein family. As for other thioredoxin‐fold proteins, the pK_a of the nucleophilic cysteine of NrdH‐redoxins is of particular interest since it affects the catalytic reaction rate of the enzymes. Recently, the pK_a value of this cysteine in Corynebacterium glutamicum and Mycobacterium tuberculosis NrdH‐redoxins were determined, but structural insights explaining the relatively low pK_a remained elusive. We subjected C. glutamicum NrdH‐redoxin to an extensive molecular dynamics simulation to expose the factors regulating the pK_a of the nucleophilic cysteine. We found that the nucleophilic cysteine receives three hydrogen bonds from residues within the CVQC active site motif. Additionally, a fourth hydrogen bond with a lysine located N‐terminal of the active site further lowers the cysteine pK_a. However, site‐directed mutagenesis data show that the major contribution to the lowering of the cysteine pK_a comes from the positive charge of the lysine and not from the additional Lys‐Cys hydrogen bond. In 12% of the NrdH‐redoxin family, this lysine is replaced by an arginine that also lowers the cysteine pK_a. All together, the four hydrogen bonds and the electrostatic effect of a lysine or an arginine located N‐terminally of the active site dynamically regulate the pK_a of the nucleophilic cysteine in NrdH‐redoxins.     

Spectroscopy and computational studies on the interaction of octyl,dodecyl, and hexadecyl derivatives of anionic and cationic surfactants with adenosine deaminase

Effects of sodium (octyl, dodecyl, hexadecyl) sulfate and their cationic analogous on the structure of adenosine deaminase (ADA) were investigated by fluorescence and circular dichroism spectroscopy as well as molecular dynamics simulation and docking calculation. Root-mean-square derivations, radius of gyration, solvent accessible surface area, and radial distribution function were obtained. The results showed that anionic and cationic surfactants reduce protein stability. Cationic surfactants have more effect on the ADA structure in comparison with anionic surfactants. More concentration and longer surfactants are parallel to higher denaturation. Furthermore, aggregation in the presence of anionic surfactants is more than cationic surfactants. Docking data showed that longer surfactants have more interaction energy and smaller ones bound to the active site.     

Structural basis of pesticide detection by enzymatic biosensing: a molecular docking and MD simulation study

Designing of rapid, facile, selective, and cost-effective biosensor technology is a growing area for the detection of various classes of pesticides. The biosensor with these features can be achieved only through the various bio-components using different transducers. This study, therefore, focuses on the usage of molecular docking, specificity tendencies, and capabilities of proteins for the detection of pesticides. Accordingly, the four transducers, acetylcholinesterase (ACH), cytochromes P450 (CYP), glutathione S-transferase (GST), and protein kinase C (PKC) were selected based on their applications including neurotransmitter, metabolism, detoxification enzyme, and protein phosphorylation. Then after molecular docking of the pesticides, fenobucarb, dichlorodiphenyltrichloroethane (DDT), and parathion onto each enzyme, the conformational behavior of the most stable complexes was further analyzed using 50 ns Molecular Dynamics (MD) simulations carried out under explicit water conditions. In the case of protein kinase C (PKC) and cytochrome P450 3A4 enzyme (CYP), the fenobucarb complex showed the most suitable combination of free energy of binding and inhibition constant ?4.42 kcal/mol (573.73 μM) and ?5.1 kcal/mol (183.49 μM), respectively. Parathion dominated for acetylcholinesterase (ACH) with ?4.57 kcal/mol (448.09 μM) and lastly dichlorodiphenyltrichloroethane for glutathione S-transferase (GST), ?5.43 kcal/mol (103.88 μM). The RMSD variations were critical for understanding the impact of pesticides as they distinctively influence the energetic attributes of the proteins. Overall, the outcomes from the extensive analysis provide an insight into the structural features of the proteins studied, thereby highlighting their potential use as a substrate in biorecognition sensing of pesticide compounds.     

Repurposing FDA-approved drugs as FXR agonists: a structure based in silico pharmacological study

Farnesoid X receptor (FXR) modulates the expression of genes involved in lipid and carbohydrate homeostasis and inflammatory processes. This nuclear receptor is likely a tumor suppressor in several cancers, but its molecular mechanism of suppression is still under study. Several studies reported that FXR agonism increases the survival of colorectal, biliary tract, and liver cancer patients. In addition, FXR expression was shown to be down-regulated in many diseases such as obesity, irritable bowel syndrome, glomerular inflammation, diabetes, proteinuria, and ulcerative colitis. Therefore, development of novel FXR agonists may have significant potential in the prevention and treatment of these diseases. In this scenario, computer-aided drug design procedures can be resourcefully applied for the rapid identification of promising drug candidates. In the present study, we applied the molecular docking method in conjunction with molecular dynamics (MD) simulations to find out potential agonists for FXR based on structural similarity with the drug that is currently used as FXR agonist, obeticholic acid. Our results showed that alvimopan and montelukast could be used as potent FXR activators and outperform the binding affinity of obeticholic acid by forming stable conformation with the protein in silico. However, further investigational studies and validations of the selected drugs are essential to figure out their suitability for preclinical and clinical trials.     

The role of S477N mutation in the molecular behavior of SARS-CoV-2 spike protein: An in-silico perspective

The attachment of SARA-CoV-2 happens between ACE2 and the receptor binding domain (RBD) on the spike protein. Mutations in this domain can affect the binding affinity of the spike protein for ACE2. S477N, one of the most common mutations reported in the recent variants, is located in the RBD. Today's computational approaches in biology, especially during the SARS-CoV-2 pandemic, assist researchers in predicting a protein's behavior in contact with other proteins in more detail. In this study, we investigated the interactions of the S477N-hACE2 in silico to find the impact of this mutation on its binding affinity for ACE2 and immunity responses using dynamics simulation, protein–protein docking, and immunoinformatics methods. Our computational analysis revealed an increased binding affinity of N477 for ACE2. Four new hydrogen and hydrophobic bonds in the mutant RBD-ACE2 were formed (with S19 and Q24 of ACE2), which do not exist in the wild type. Also, the protein spike structure in this mutation was associated with an increase in stabilization and a decrease in its fluctuations at the atomic level. N477 mutation can be considered as the cause of increased escape from the immune system through MHC-II.     

The role of conserved water molecules in the catalytic domain of protein kinases

Protein kinases are essential signaling molecules with a characteristic bilobal shape that has been studied for over 15 years. Despite the number of crystal structures available, little study has been directed away from the prototypical functional elements of the kinase domain. We have performed a structural alignment of 13 active‐conformation kinases and discovered the presence of six water molecules that occur in conserved locations across this group of diverse kinases. Molecular dynamics simulations demonstrated that these waters confer a great deal of stability to their local environment and to a key catalytic residue. Our results highlight the importance of novel elements within the greater kinase family and suggest that conserved water molecules are necessary for efficient kinase function. Proteins 2009. © 2009 Wiley‐Liss, Inc.     

A conserved arginine plays a role in the catalytic cycle of the protein disulphide isomerases

The pK(a) values of the CXXC active-site cysteine residues play a critical role in determining the physiological function of the thioredoxin superfamily. To act as an efficient thiol-disulphide oxidant the thiolate state of the N-terminal cysteine must be stabilised and the thiolate state of the C-terminal cysteine residue destabilised. While increasing the pK(a) value of the C-terminal cysteine residue promotes oxidation of substrates, it has an inhibitory effect on the reoxidation of the enzyme, which is promoted by the formation of a thiolate at this position. Since reoxidation is essential to complete the catalytic cycle, the differential requirement for a high and a low pK(a) value for the C-terminal cysteine residue for different steps in the reaction presents us with a paradox. Here, we report the identification of a conserved arginine residue, located in the loop between beta5 and alpha4 of the catalytic domains of the human protein disulphide isomerase (PDI) family, which is critical for the catalytic function of PDI, ERp57, ERp72 and P5, specifically for reoxidation. An examination of the published NMR structure for the a domain of PDI combined with molecular dynamic studies suggest that the side-chain of this arginine residue moves into and out of the active-site locale and that this has a very marked effect on the pK(a) value of the active-site cysteine residues. This intra-domain motion resolves the apparent dichotomy of the pK(a) requirements for the C-terminal active-site cysteine.     

ApoE: The role of conserved residues in defining function

The amino acid sequences of apolipoprotein E (apoE) from 63 different mammalian species have been downloaded from the protein database. The sequences were compared to human apoE4 to determine conserved and non‐conserved sequences of amino acids. ApoE4 is the major risk factor for the development of late onset Alzheimer's disease while apoE3, which differs from apoE4 by a single amino acid change at position 112, poses little or no risk for the development of this disease. Thus, the two proteins appear to be structurally and functionally different. Seven highly conserved regions, representing approximately 47 amino acids (of 299) have been found. These regions are distributed throughout the protein and reflect ligand binding sites as well as regions proposed to be involved in the propagation of the cysteine–arginine change at position 112 to distant regions of the protein in the N‐ and C‐terminal domains. Highly non‐conserved regions are at the N‐ and C‐terminal ends of the apoE protein.     

Unfolding crystallins: the destabilizing role of a beta-hairpin cysteine in betaB2-crystallin by simulation and experiment

The thermodynamic and kinetic stabilities of the eye lens family of betagamma-crystallins are important factors in the etiology of senile cataract. They control the chance of proteins unfolding, which can lead to aggregation and loss of transparency. betaB2-Crystallin orthologs are of low stability and comprise two typical betagamma-crystallin domains, although, uniquely, the N-terminal domain has a cysteine in one of the conserved folded beta-hairpins. Using high-temperature (500 K) molecular dynamics simulations with explicit solvent on the N-terminal domain of rodent betaB2-crystallin, we have identified in silico local flexibility in this folded beta-hairpin. We have shown in vitro using two-domain human betaB2-crystallin that replacement of this cysteine with a more usual aromatic residue (phenylalanine) results in a gain in conformational stability and a reduction in the rate of unfolding. We have used principal components analysis to visualize and cluster the coordinates from eight separate simulated unfolding trajectories of both the wild-type and the C50F mutant N-terminal domains. These data, representing fluctuations around the native well, show that although the mutant and wild-type appear to behave similarly over the early time period, the wild type appears to explore a different region of conformational space. It is proposed that the advantage of having this low-stability cysteine may be correlated with a subunit-exchange mechanism that allows betaB2-crystallin to interact with a range of other beta-crystallin subunits.     

An integrated approach of network-based systems biology,molecular docking,and molecular dynamics approach to unravel the role of existing antiviral molecules against AIDS-associated cancer

A serious challenge in cancer treatment is to reposition the activity of various already known drug candidates against cancer. There is a need to rewrite and systematically analyze the detailed mechanistic aspect of cellular networks to gain insight into the novel role played by various molecules. Most Human Immunodeficiency Virus infection-associated cancers are caused by oncogenic viruses like Human Papilloma Viruses and Epstein–Bar Virus. As the onset of AIDS-associated cancers marks the severity of AIDS, there might be possible interconnections between the targets and mechanism of both the diseases. We have explored the possibility of certain antiviral compounds to act against major AIDS-associated cancers: Kaposi’s Sarcoma, Non-Hodgkin Lymphoma, and Cervical Cancer with the help of systems pharmacology approach that includes screening for targets and molecules through the construction of a series of drug–target and drug–target–diseases network. Two molecules (Calanolide A and Chaetochromin B) and the target “HRAS” were finally screened with the help of molecular docking and molecular dynamics simulation. The results provide novel antiviral molecules against HRAS target to treat AIDS defining cancers and an insight for understanding the pharmacological, therapeutic aspects of similar unexplored molecules against various cancers.