Salt Bridges and Conformational Flexibility: Effect on Protein Stability

Proposal of pharmacophore model for HIV reverse transcriptase inhibitors: Combined mutational effect analysis, molecular dynamics, molecular docking and pharmacophore modeling study

OBJECTIVES

Antiretroviral therapy (ART) efficacy is jeopardized by the emergence of drug resistance mutations in HIV, compromising treatment effectiveness. This study aims to propose novel analogs of Effavirenz (EFV) as potential direct inhibitors of HIV reverse transcriptase, employing computer-aided drug design methodologies.

METHODS

Three key approaches were applied: a mutational profile study, molecular dynamics simulations, and pharmacophore development. The impact of mutations on the stability, flexibility, function, and affinity of target proteins, especially those associated with NRTI, was assessed. Molecular dynamics analysis identified G190E as a mutation significantly altering protein properties, potentially leading to therapeutic failure. Comparative analysis revealed that among six first-line antiretroviral drugs, EFV exhibited notably low affinity with viral reverse transcriptase, further reduced by the G190E mutation. Subsequently, a search for EFV-similar inhibitors yielded 12 promising molecules based on their affinity, forming the basis for generating a pharmacophore model.

RESULTS

Mutational analysis pinpointed G190E as a crucial mutation impacting protein properties, potentially undermining therapeutic efficacy. EFV demonstrated diminished affinity with viral reverse transcriptase, exacerbated by the G190E mutation. The search for EFV analogs identified 12 high-affinity molecules, culminating in a pharmacophore model elucidating key structural features crucial for potent inhibition.

CONCLUSION

This study underscores the significance of EFV analogs as potential inhibitors of HIV reverse transcriptase. The findings highlight the impact of mutations on drug efficacy, particularly the detrimental effect of G190E. The generated pharmacophore model serves as a pivotal reference for future drug development efforts targeting HIV, providing essential structural insights for the design of potent inhibitors based on EFV analogs identified in vitro.

patcHwork: a user-friendly pH sensitivity analysis web server for protein sequences and structures

pH regulates protein function and interactions by altering the charge of individual residues causing loss or gain of intramolecular noncovalent bonds, which may lead to structural rearrangements. While tools to analyze residue-specific charge distribution of proteins at a given pH exist, currently no tool is available to investigate noncovalent bond changes at two different pH values. To make protein pH sensitivity analysis more accessible, we developed patcHwork, a web server that combines the identification of amino acids undergoing a charge shift with the determination of affected noncovalent bonds at two user-defined pH values. At the sequence-only level, patcHwork applies the Henderson–Hasselbalch equation to determine pH-sensitive residues. When the 3D protein structure is available, patcHwork can be employed to gain mechanistic understanding of the effect of pH. This is achieved using the PDB2PQR and PROPKA tools and noncovalent bond determination algorithms. A user-friendly interface allows visualizing pH-sensitive residues, affected salt bridges, hydrogen bonds and aromatic (pi–pi and cation–pi) interactions. patcHwork can be used to identify patches, a new concept we propose of pH-sensitive residues in close proximity on the protein, which may have a major impact on function. We demonstrate the attractiveness of patcHwork studying experimentally investigated pH-sensitive proteins (https://patchwork.biologie.uni-freiburg.de/).

Highly Similar Sequence and Structure Yet Different Biophysical Behavior: A Computational Study of Two Triosephosphate Isomerases

Homodimeric triosephosphate isomerases (TIMs) from Trypanosoma cruzi (TcTIM) and Trypanosoma brucei (TbTIM) have markedly similar amino-acid sequences and three-dimensional structures. However, several of their biophysical parameters, such as their susceptibility to sulfhydryl agents and their reactivation speed after being denatured, have significant differences. The causes of these differences were explored with microsecond-scale molecular dynamics (MD) simulations of three different TIM proteins: TcTIM, TbTIM, and a chimeric protein, Mut1. We examined their electrostatic interactions and explored the impact of simulation length on them. The same salt bridge between catalytic residues Lys 14 and Glu 98 was observed in all three proteins, but key differences were found in other interactions that the catalytic amino acids form. In particular, a cation-π interaction between catalytic amino acids Lys 14 and His 96 and both a salt bridge and a hydrogen bond between catalytic Glu 168 and residue Arg 100 were only observed in TcTIM. Furthermore, although TcTIM forms less hydrogen bonds than TbTIM and Mut1, its hydrogen bond network spans almost the entire protein, connecting the residues in both monomers. This work provides new insight into the mechanisms that give rise to the different behavior of these proteins. The results also show the importance of long simulations.

Controlling the Substrate Specificity of an Enzyme through Structural Flexibility by Varying the Salt-Bridge Density

Many enzymes, particularly in one single family, with highly conserved structures and folds exhibit rather distinct substrate specificities. The underlying mechanism remains elusive, the resolution of which is of great importance for biochemistry, biophysics, and bioengineering. Here, we performed a neutron scattering experiment and molecular dynamics (MD) simulations on two structurally similar CYP450 proteins; CYP101 primarily catalyzes one type of ligands, then CYP2C9 can catalyze a large range of substrates. We demonstrated that it is the high density of salt bridges in CYP101 that reduces its structural flexibility, which controls the ligand access channel and the fluctuation of the catalytic pocket, thus restricting its selection on substrates. Moreover, we performed MD simulations on 146 different kinds of CYP450 proteins, spanning distinct biological categories including Fungi, Archaea, Bacteria, Protista, Animalia, and Plantae, and found the above mechanism generally valid. We demonstrated that, by fine changes of chemistry (salt-bridge density), the CYP450 superfamily can vary the structural flexibility of its member proteins among different biological categories, and thus differentiate their substrate specificities to meet the specific biological needs. As this mechanism is well-controllable and easy to be implemented, we expect it to be generally applicable in future enzymatic engineering to develop proteins of desired substrate specificities.

Factors underlying asymmetric pore dynamics of disaggregase and microtubule-severing AAA+ machines

Disaggregation and microtubule-severing nanomachines from the AAA+ (ATPases associated with various cellular activities) superfamily assemble into ring-shaped hexamers that enable protein remodeling by coupling large-scale conformational changes with application of mechanical forces within a central pore by loops protruding within the pore. We probed the asymmetric pore motions and intraring interactions that support them by performing extensive molecular dynamics simulations of single-ring severing proteins and the double-ring disaggregase ClpB. Simulations reveal that dynamic stability of hexameric pores of severing proteins and of the nucleotide-binding domain 1 (NBD1) ring of ClpB, which belong to the same clade, involves a network of salt bridges that connect conserved motifs of central pore loops. Clustering analysis of ClpB highlights correlated motions of domains of neighboring protomers supporting strong interprotomer collaboration. Severing proteins have weaker interprotomer coupling and stronger intraprotomer stabilization through salt bridges involving pore loops. Distinct mechanisms are identified in the NBD2 ring of ClpB involving weaker interprotomer coupling through salt bridges formed by noncanonical loops and stronger intraprotomer coupling. Analysis of collective motions of PL1 loops indicates that the largest amplitude motions in the spiral complex of spastin and ClpB involve axial excursions of the loops, whereas for katanin they involve opening and closing of the central pore. All three motors execute primarily axial excursions in the ring complex. These results suggest distinct substrate processing mechanisms of remodeling and translocation by ClpB and spastin compared to katanin, thus providing dynamic support for the differential action of the two severing proteins. Relaxation dynamics of the distance between the PL1 loops and the center of mass of protomers reveals observation-time-dependent dynamics, leading to predicted relaxation times of tens to hundreds of microseconds on millisecond experimental timescales. For ClpB, the predicted relaxation time is in excellent agreement with the extracted time from smFRET experiments.

A Geometric Definition of Short to Medium Range Hydrogen-Mediated Interactions in Proteins

We present a method to rapidly identify hydrogen-mediated interactions in proteins (e.g., hydrogen bonds, hydrogen bonds, water-mediated hydrogen bonds, salt bridges, and aromatic π-hydrogen interactions) through heavy atom geometry alone, that is, without needing to explicitly determine hydrogen atom positions using either experimental or theoretical methods. By including specific real (or virtual) partner atoms as defined by the atom type of both the donor and acceptor heavy atoms, a set of unique angles can be rapidly calculated. By comparing the distance between the donor and the acceptor and these unique angles to the statistical preferences observed in the Protein Data Bank (PDB), we were able to identify a set of conserved geometries (15 for donor atoms and 7 for acceptor atoms) for hydrogen-mediated interactions in proteins. This set of identified interactions includes every polar atom type present in the Protein Data Bank except OE1 (glutamate/glutamine sidechain) and a clear geometric preference for the methionine sulfur atom (SD) to act as a hydrogen bond acceptor. This method could be readily applied to protein design efforts.

The antioxidant function of Sco proteins depends on a critical surface-exposed residue

BACKGROUND

Besides their role in copper metabolism, Sco proteins from different organisms have been shown to play a defensive role against oxidative stress. In the present study, we set out to identify crucial amino acid residues for the antioxidant activity.

METHODS

Native and mutated Sco proteins from human, Arabidopsis thaliana and the yeast Kluyveromyces lactis were expressed in the model organism Saccharomyces cerevisiae. The oxidative stress resistance of the respective transformants was determined by growth and lipid peroxidation assays.

RESULTS

A functionally important site, located 15 amino acids downstream of the well-conserved copper binding CxxxC motif, was identified. Mutational analysis revealed that a positive charge at this position has a detrimental effect on the antioxidant capacity. Bioinformatic analysis predicts that this site is surface-exposed, and according to Co-IP data it is required for binding of proteins that are connected to known antioxidant pathways.

CONCLUSION

This study shows that the antioxidant capacity of eukaryotic Sco proteins is conserved and depends on the presence of functional site(s) rather than the extent of overall sequence homology.

GENERAL SIGNIFICANCE

These findings provide an insight into the conserved functional sites of eukaryotic Sco proteins that are crucial for combating oxidative stress. This capacity is probably not due to an enzymatic activity but rather is indirectly mediated by interaction with other proteins.

Comparative bioinformatic and structural analyses of pepsin and renin

Pepsin, the archetypal pepsin-like aspartic protease, is irreversibly denatured when exposed to neutral pH conditions whereas renin, a structural homologue of pepsin, is fully stable and optimally active in the same conditions despite sharing highly similar enzyme architecture. To gain insight into the structural determinants of differential aspartic protease pH stability, the present study used comparative bioinformatic and structural analyses. In pepsin, an abundance of polar and aspartic acid residues were identified, a common trait with other acid-stable enzymes. Conversely, renin was shown to have increased levels of basic amino acids. In both pepsin and renin, the solvent exposure of these charged groups was high. Having similar overall acidic residue content, the solvent-exposed basic residues may allow for extensive salt bridge formation in renin, whereas in pepsin, these residues are protonated and serve to form stabilizing hydrogen bonds at low pH. Relative differences in structure and sequence in the turn and joint regions of the β-barrel and ψ-loop in both the N- and C-terminal lobes were identified as regions of interest in defining divergent pH stability. Compared to the structural rigidity of renin, pepsin has more instability associated with the N-terminus, specifically the B/C connector. By contrast, renin exhibits greater C-terminal instability in turn and connector regions. Overall, flexibility differences in connector regions, and amino acid composition, particularly in turn and joint regions of the β-barrel and ψ-loops, likely play defining roles in determining pH stability for renin and pepsin.

The Effects of One Amino Acid Substitutions at the C-Terminal Region of Thermostable L2 Lipase by Computational and Experimental Approach

The substitutions of the amino acid at the predetermined critical point at the C-terminal of L2 lipase may increase its thermostability and enzymatic activity, or even otherwise speed up the unfolding of the protein structure. The C-terminal of most proteins is often flexible and disordered. However, some protein functions are directly related to flexibility and play significant role in enzyme reaction. The critical point for mutation of L2 lipase structure was predicted at the position 385 of the L2 sequence, and the best three mutants were determined based on I-Mutant2.0 software. The best three mutants were S385E, S385I and S385V. The effects of the substitution of the amino acids at the critical point were analysed with molecular dynamics simulation by using Yet Another Scientific Artificial Reality Application software. The predicted mutant L2 lipases were found to have lower root mean square deviation value as compared to L2 lipase. It was indicated that all the three mutants had higher compactness in the structure, consequently enhanced the stability. Root mean square fluctuation analysis showed that the flexibility of L2 lipase was reduced by mutations. Purified S385E lipase had an optimum temperature of 80 °C in Tris-HCl pH 8. The highest enzymatic activity of purified S385E lipase was obtained at 80 °C temperature in Tris-HCl pH 8, while for L2 lipase it was at 70 °C in Glycine-NaOH pH 9. The thermal stability of S385V lipase was enhanced as compared to other protein since that the melting point (T _m) value was at 85.96 °C. S385I lipase was more thermostable compared to recombinant L2 lipase and other mutants at temperature 60 °C within 16 h preincubation.

Potential role of salt-bridges in the hinge-like movement of apicomplexa specific β-hairpin of Plasmodium and Toxoplasma profilins: A molecular dynamics simulation study

Profilin is one of the actin-binding proteins that regulate dynamics of actin polymerization. It plays a key role in cell motility and invasion. It also interacts with several other proteins notably through its poly-L-proline (PLP) binding site. Profilin in apicomplexa is characterized by a unique mini-domain consisting of a large β-hairpin extension and an acidic loop which is relatively longer in Plasmodium species. Profilin is essential for the invasive blood stages of Plasmodium falciparum. In the current study, unbound profilins from Plasmodium falciparum (Pf), Toxoplasma gondii (Tg), and Homo sapiens (Hs) were subjected to molecular dynamics (MD) simulations for a timeframe of 100 ns each to understand the conformational dynamics of these proteins. It was found that the β-hairpin of profilins from Pf and Tg shows a hinge-like movement. This movement in Pf profilin may possibly be driven by the loss of a salt-bridge within profilin. The impact of this conformational change on actin binding was assessed by docking three dimensional (3D) structures of profilin from Pf and Tg with their corresponding actins using ClusPro2.0. The stability of docked Pf profilin-actin complex was assessed through a 50 ns MD simulation. As Hs profilin I does not have the apicomplexa specific mini-domain, MD simulation was performed for this protein and its dynamics was compared to that of profilins from Pf and Tg. Using an immunoinformatics approach, potential epitope regions were predicted for Pf profilin. This has a potential application in the design of vaccines as they mapped to its unique mini-domain.

Relationship between stability and flexibility in the most flexible region of Photinus pyralis luciferase

Firefly luciferase is a protein with a large N-terminal and a small C-terminal domain. B-factor analysis shows that its C-terminal is much more flexible than its N-terminal. Studies on hyperthermophile proteins have been shown that the increased thermal stability of hyperthermophile proteins is due to their enhanced conformational rigidity and the relationship between flexibility, stability and function in most of proteins is on debate. Two mutations (D474K and D476N) in the most flexible region of firefly luciferase were designed. Thermostability analysis shows that D476N mutation doesn't have any significant effect but D474K mutation destabilized protein. On the other hand, flexibility analysis using dynamic quenching and limited proteolysis demonstrates that D474K mutation became much more flexible than wild type although D476N doesn't have any significant difference. Intrinsic and ANS fluorescence studies demonstrate that D476N mutation is brought about by structural changes without significant effect on thermostability and flexibility. Molecular modeling reveals that disruption of a salt bridge between D(474) and K(445) accompanying with some H-bond deletion may be involved in destabilization of D474K mutant.