Multiple molecular dynamics simulations of TEM beta-lactamase: dynamics and water binding of the omega-loop

Structure and dynamics of whole-sequence homology model of ORF3a protein of SARS-CoV-2: An insight from microsecond molecular dynamics simulations

The ORF3a is a large accessory protein in severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which plays an important role in virulence and viral replication; especially in inflammasome activation and apoptosis. However,, the existing cryo-EM structure of SARS-CoV-2 ORF3a is incomplete, . making it challenging to understand its structural and functional features. The aim of this study is to investigate the dynamic behaviors of the full-sequence homology model of ORF3a and compare it with the cryo-EM structure using microsecond molecular dynamics simulations. The previous studies indicated that the unresolved residues of the cryo-EM structure are not only involved in the pathogenesis of the SARS-CoV-2 but also exhibit a significant antigenicity. The dynamics scenario of homology model revealed higher RMSD, Rg, and SASA values with stable pattern when compared to the cryo-EM structure. Moreover, the RMSF analysis demonstrated higher fluctuations at specific positions (1-43, 97-110, 172-180, 219-243) in the model structure, whereas the cryo-EM structure displayed lower overall drift (except 1-43) in comparison to the model structure.Secondary structural features indicated that a significant unfolding in the transmembrane domains and β-strand at positions 166 to 172, affecting the stability and compactness of the cryo-EM structure , whereas the model exhibited noticeable unfolding in transmembrane domains and small-coiled regions in the N-terminal. , The results from molecular docking and steered molecular dynamics investigations showed the model structure had a greater number of non-bonding interactions, leading to enhanced stability when compared to the cryo-EM structure. Consequently, higher forces were necessary for unbinding of the baricitinib and ruxolitinib inhibitors from the model structure.. Our findings can help better understanding of the significance of unresolved residues at the molecular level. Additionally, this information can guide researchers for experimental endeavors aimed at completing the full-sequence structure of the ORF3a.Communicated by Ramaswamy H. Sarma.

Molecular cloning, characterization and 3D modelling of spotted snakehead fbn1 C-terminal region encoding asprosin and expression analysis of fbn1

The FBN1 gene encodes profibrillin protein that is cleaved by the enzyme furin to release fibrillin-1 and a glucogenic hormone, asprosin. Asprosin is implicated in diverse metabolic functions as well as pathological conditions in mammals. However, till date, there are no studies on asprosin in any non-mammalian vertebrate. In this study, we have retrieved the spotted snakehead Channa punctata fbn1 gene (ss fbn1) from the testicular transcriptome data and validated it. The transcript is predicted to encode 2817 amino acid long putative profibrillin protein. Amino acid sequence alignment of deduced ss profibrillin with human profibrillin revealed that the furin cleavage site in profibrillin is well conserved in C. punctata. Further, differential expression of ss fbn1 was observed in various tissues with the highest expression in gonads. Prominent expression of furin was also observed in the gonads suggesting the possibility of proteolytic cleavage of profibrillin protein and secretion of asprosin in C. punctata. In addition, the C-terminal of the fbn1 gene of C. punctata that codes for asprosin protein has been cloned. Using in silico approach, physicochemical properties of the putative ss asprosin were characterized and post-translational changes were predicted. The putative ss asprosin protein sequence is predicted to consist of 142 amino acid residues, with conserved glycosylation sites. Further, the 3D model of ss asprosin was predicted followed by MD (molecular dynamics) simulation for energy minimization. Thus, the current study, for the first time in non-mammalian vertebrates, predicts and characterizes the novel protein asprosin using in silico approach.

Detailed investigation of catalytically important residues of class A β-lactamase

An increasing global health challenge is antimicrobial resistance. Bacterial infections are often treated by using β-lactam antibiotics. But several resistance mechanisms have evolved in clinically mutated bacteria, which results in resistance against such antibiotics. Among which production of novel β-lactamase is the major one. This results in bacterial resistance against penicillin, cephalosporin, and carbapenems, which are considered to be the last resort of antibacterial treatment. Hence, β-lactamase enzymes produced by such bacteria are called extended-spectrum β-lactamase and carbapenemase enzymes. Further, these bacteria have developed resistance against many β-lactamase inhibitors as well. So, investigation of important residues that play an important role in altering and expanding the spectrum activity of these β-lactamase enzymes becomes necessary. This review aims to gather knowledge about the role of residues and their mutations in class A β-lactamase, which could be responsible for β-lactamase mediated resistance. Class A β-lactamase enzymes contain most of the clinically significant and expanded spectrum of β-lactamase enzymes. Ser70, Lys73, Ser130, Glu166, and Asn170 residues are mostly conserved and have a role in the enzyme's catalytic activity. In-depth investigation of 69, 130, 131, 132, 164, 165, 166, 170, 171, 173, 176, 178, 179, 182, 237, 244, 275 and 276 residues were done along with its kinetic analysis for knowing its significance. Further, detailed information from many previous studies was gathered to know the effect of mutations on the kinetic activity of class A β-lactamase enzymes with β-lactam antibiotics.Communicated by Ramaswamy H. Sarma.

Unveiling the mechanistic roles of chlorine substituted phthalazinone-based compounds containing chlorophenyl moiety towards the differential inhibition of poly (ADP-ribose) polymerase-1 in the treatment of lung cancer

PARP-1 has become an attractive target in cancer treatment owing to its significant role in breast and ovarian cancers. The design of highly selective and effective poly (ADP ribose) polymerase-1 inhibitors has significant therapeutic advantages and has remained the core of several PARP-1-based drug discovery research. The pharmacophoric relevance of a chlorine substituent in a recent study led to the design of compounds 11c (meta-chlorophenyl) and 11d (para-chlorophenyl). In this study, we resolved the mechanistic effects of the changes in chlorine positional orientation, which underlie the inhibitory potencies and selectivity exhibited disparately by 11c and 11d. Compared to 11d, among other multiple higher-affinity complementary interactions with key site residues, the meta-Cl positioning in 11c facilitated its optimal motion and orientation towards conserved residues Arg878 and Asp766 with consistent pi-cation and pi-anion interactions, respectively, thereby favoring the stability of the ligand towards PARP-1. These could account for the higher inhibitory potency exhibited by 11c relative to 11d against PARP-1. The thermodynamics calculation revealed that 11c had a relatively higher total binding energy (ΔGbind) than 11d. We also observed that 11d displayed high deviations, compared to 11c, indicative of its unstable binding orientation. Furthermore, we reported in this study that the high involvement of electrostatic and van der Waal effects potentiated the binding affinity and strength of 11c (ΔE_vdW = -50.58 and ΔE_ele = -27.20) relative to 11d (ΔE_vdW = -49.46 and ΔE_ele = -19.96) at PARP-1 binding pocket. We believe the findings in this current study would provide valuable insights into the design of selective PARP-1 inhibitors containing chlorine substituent for cancer treatment, including lung cancer.Communicated by Ramaswamy H. Sarma.

Talazoparib Dual-targeting on Poly (ADP-ribose) Polymerase-1 and -16 Enzymes Offers a Promising Therapeutic Strategy in Small Cell Lung Cancer Therapy: Insight from Biophysical Computations

In recent times, inhibition of poly (ADP-ribose) polymerase (PARP) enzymes by pharmacological drugs has attracted much attention as an anticancer therapy. As reported, PARP-16 has been discovered as a novel anticancer target for small cell lung cancer, and that the inhibition of both PARP-16 and PARP-1 by talazoparib can increase the overall effectiveness of talazoparib in the SCLC treatment. In this study, we employed computational approaches to investigate the differential inhibitory potency of Talazoparib on PARP-1 and PARP-16. Talazoparib has excellent PARP-1 and PARP-16 binding activities, as revealed by the ΔG_bind (total binding energy). Pp16-tpb had binding energy of -34.85 kcal/mol, while pp1-tpb had a binding energy of -26.36 kcal/mol. The binding activity of Talazoparib on both PARP-1 and PARP-16 was significantly influenced by van der Waal and electrostatic interactions. Correspondingly, according to the findings of this study, binding residues with total binding energy greater than 1.00 kcal/mol contributed considerably to the Talazoparib's binding activities on PARP-1 and PARP-16. We believe the findings of this study will pave the way for developing dual targeting of PARP enzymes as a strategy for small-cell lung cancer treatment.

Prospecting the therapeutic edge of a novel compound (B12) over berberine in the selective targeting of Retinoid X Receptor in colon cancer

The Retinoid X Receptor (RXR) is an attractive target in the treatment of colon cancer. Different therapeutic binders with high potency have been used to specifically target RXR. Among these compounds is a novel analogue of berberine, B12. We provided structural and molecular insights into the therapeutic activity properties of B12 relative to its parent compound, berberine, using force field estimations and thermodynamic calculations. Upon binding of B12 to RXR, the high instability elicited by RXR was markedly reduced; similar observation was seen in the berberine-bound RXR. However, our analysis revealed that B12 could have a more stabilizing effect on RXR when compared to berberine. Interestingly, the mechanistic behaviour of B12 in the active site of RXR opposed its impact on RXR protein. This disparity could be due to the bond formation and breaking elicited between B12/berberine and the active site residues. We observed that B12 and berberine could induce a disparate conformational change in regions Gly250-Asp258 located on the His-RXRα/LBD domain. Comparatively, the high agonistic and activation potential reported for B12 compared to berberine might be due to its superior binding affinity as evidenced in the thermodynamic estimations. The total affinity for B12 (-25.76 kcal/mol) was contributed by electrostatic interactions from Glu243 and Glu239. Also, Arg371, which plays a crucial role in the activity of RXR, formed a strong hydrogen bond with B12; however, a weak interaction was elicited between Arg371 and berberine. Taken together, our study has shown the RXRα activating potential of B12, and findings from this study could provide a framework in the future design of RXRα binders specifically tailored in the selective treatment of colon cancer.

Insight into the Therapeutic Potential of a Bicyclic Hydroxypyridone Compound 2-[(2,4-Dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one as COMT Inhibitor in the Treatment of Parkinson's Disease: A Molecular Dynamic Simulation Approach

Parkinson's disease (PD) is one of the most targeted neurodegenerative diseases in clinical research. Awareness of research is due to its increasing number of affected people worldwide. The pathology of PD has been linked to several key proteins upregulation such as the catechol O-Methyltransferase (COMT). Hence, the synthesis of compounds possessing inhibitory capacity has been the frontline of research in recent years. Several compounds have been synthesized among which is the nitrocatechol. However, major limitations associated with the nitrocatechol scaffold include the inability to possess adequate CNS penetration properties and hepatic toxicity associated with the compounds. However, a series of bicyclic hydroxypyridones compounds were synthesized to evaluate their inhibitory potentials on COMT protein with compound 38 (c38) 2-[(2,4-dichlorophenyl)methyl]-7-hydroxy-1,2,3,4-tetrahydro-8H-pyrido[1,2-a]pyrazin-8-one shown to have a 40 fold increase level coverage in its IC₅₀ over brain exposure when compared to the other synthesized compound. The molecular dynamics method was employed to understand the nature of interaction exhibited by c38. Molecular mechanics of c38 revealed a disruptive effect on the secondary structure of COMT protein. Per residue decomposition analysis revealed similar crucial residues involved in the favorable binding of c38 and tolcapone implicated its increased inhibitory capacity on COMT in preventing PD. Free binding energy (ΔG_bind ) of c38 further revealed the inhibitory capacity towards COMT protein in comparison to the FDA approved tolcapone. Ligand mobility analysis of both compounds showed a timewise different mobility pattern across the simulation time frame at the active site pocket of the protein connoting the different inhibitory potency exhibited by c38 and tolcapone. Findings from this study revealed optimization of c38 could facilitate the discovery of new compounds with enhanced inhibitory properties towards COMT in treating PD.

Piecing the fragments together: Dynamical insights into the enhancement of BRD4-BD1 (BET protein) druggability in cancer chemotherapy using novel 8-methyl-pyrrolo[1,2-a]pyrazin-1(2H)-one derivatives

BACKGROUND

Fragment-based drug discovery in recent times has been explored in the design of highly potent therapeutics.

METHODS

In this study, we explored the inhibitory dynamics of Compound 38 (Cpd38), a newly synthesized Bromodomain-containing protein 4 bromodomain 1 (BRD4-BD1) protein inhibitor derived from the synthetic coupling of Fragment 47 (Fgt47) into ABBV-075 scaffold. Using dynamic simulation methods, we unraveled the augmentative effects of chemical fragmentation on improved BRD4-BD1 inhibition.

RESULTS

Findings from this study revealed that although Fgt47 exhibited a considerable ΔGbind, its incorporation into the difluoro-phenoxy pyridine scaffold (Cpd38) notably enhanced the binding affinity. Time-based analyses of interaction dynamics further revealed that the bulkiness of Cpd38 favored its interaction at the BRD4-BD1 active site relative to the fragment. Strikingly, when compared to Fgt47, Cpd38 demonstrated high mobility, which could have enabled it to bind optimally and complementarily with key residues of the active site such as Ile146, Asn140, Cys136, Tyr98, Leu94, Val87, Phe83 and Trp81.

DISCUSSION

On the contrary, majority of these interactions were gradually lost in Fgt47 which could further indicate the essence of coupling it with the difluoro-phenoxy pyridine scaffold. Furthermore, Cpd38 had a more altering effect on BRD4-BDI relative to Fgt47 which could also be as a result of its higher inhibitory activity.

CONCLUSION

Conclusively, the design of highly potent therapeutics could be facilitated by the incorporation of pharmacologically active small molecule fragments into the scaffold of existing drugs.

East to West not North-West: Structure-Based Mechanistic Resolution of 8-Hydroxyl Replacement and Resulting Effects on the Activities of Imidazole-Based Heme Oxygenase-1 Inhibitors

Upregulation of Heme Oxygenase-1 (HO-1) has been widely implicated in cancer growth and chemoresistance. This explains the numerous drug discovery efforts aimed at mitigating its pro-carcinogenic roles till date. In a recent study, two selective azole-based HO-1 inhibitors (Cpd1 and Cpd2) were synthesized, which exhibited differential inhibitory potencies of ~200μm. Interestingly, variations in the affinities of these compounds were determined by their positioning across specific regions of the HO-1 binding domain, pin-pointing a pharmacological interrelationship that remains unresolved. Therefore, in this study, using molecular dynamics simulations and binding free energy calculations, we investigate how dynamical orientations of these compounds influence their binding affinities at the active HO-1 domain. Findings revealed favorable binding for the bromobenzene and imidazole substituents of Cpd1 at the western and eastern regions of the HO-1 active domain. The constituent hydroxyl group was coordinated by residues Asp140 and Arg136 over the simulation period. On the contrary, stable binding of the bromobenzene and imidazole substituents were negated by the optimal orientations of the benzyl substituent, which extended into the northeastern region. These were supported by the displacement of Asp140 and Arg136, crucial for hydrogen bond formation in Cpd1. Also, we observed that Cpd2 exhibited high deviations indicative of an unstable binding relative to Cpd1. This further supports the presumption that Cpd2 was systematically oriented away from the active HO-1 region, a phenomenon that was due to the optimal motions of the benzyl group at the northeastern regions. The highlight of our findings is that the benzyl substituent in Cpd2 elicited negative effects on HO-1, vis a vis, instability, displacement of crucial residues, and low binding energy when compared to Cpd1. Findings pave the way for future drug discovery efforts related to HO-1 inhibition in cancer therapy.

Investigating the Mechanistic Inhibitory Discrepancies of Novel Halogen and Alkyl Di-Substituted Oxadiazole-Based Dibenzo-Azepine-Dione Derivatives on Poly (ADP-Ribose) Polymerase-1

Numerous studies have established the involvement of Poly (ADP-ribose) Polymerase-1 (PARP-1) in cancer presenting it as an important therapeutic target over recent years. Although homology among the PARP protein family makes selective targeting difficult, two compounds [d11 (0.939 μM) and d21 (0.047 μM)] with disparate inhibitory potencies against PARP-1 were recently identified. In this study, free energy calculations and molecular simulations were used to decipher underlying mechanisms of differential PARP-1 inhibition exhibited by the two compounds. The thermodynamics calculation revealed that compound d21 had a relatively higher ΔG_bind than d11. High involvement of van der Waal and electrostatic effects potentiated the affinity of d21 at PARP-1 active site. More so, incorporated methyl moiety in d11 accounted for steric hindrance which, in turn, prevented complementary interactions of key site residues such as TYR889, MET890, TYR896, TYR907. Conformational studies also revealed that d21 is more stabilized for interactions in the active site compared to d11. We believe that findings from this study would provide an important avenue for the development of selective PARP-1 inhibitors.

Therapeutic Path to Double Knockout: Investigating the Selective Dual-Inhibitory Mechanisms of Adenosine Receptors A1 and A2 by a Novel Methoxy-Substituted Benzofuran Derivative in the Treatment of Parkinson's Disease

The dual inhibition of adenosine receptors A1 (A₁ AR) and A2 (A_2A AR) has been considered as an efficient strategy in the treatment of Parkinson's disease (PD). This led to the recent development of a series of methoxy-substituted benzofuran derivatives among which compound 3j exhibited dual-inhibitory potencies in the micromolar range. Therefore, in this study, we seek to resolve the mechanisms by which this novel compound elicits its selective dual targeting against A₁ AR and A_2A AR. Unique to the binding of 3j in both proteins, from our findings, is the ring-ring interaction elicited by _A1Phe275 (→ _A2Phe170) with the benzofuran ring of the compound. As observed, this π-stacking interaction contributes notably to the stability of 3j at the active sites of A₁ and A_2A AR. Besides, conserved active site residues in the proteins such as _A1Ala170 (→ _A2Ala65), _A1Ile173 (→ _A2Ile68), _A1Val191 (→ _A2Val86), _A1Leu192 (→ _A2Leu87), _A1Ala195 (→ _A2Ala90), _A1Met284 (→ _A2Met179), _A1Tyr375 (→ _A2Tyr369), _A1Ile378 (→ _A2Ile372), and _A1His382 (→ _A2His376) were commonly involved with other ring substituents which further complement the dual binding and stability of 3j. This reflects a similar interaction mechanism that involved aromatic (π) interactions. Consequentially, vdW energies contributed immensely to the dual binding of the compound, which culminated in high ΔG_binds that were homogenous in both proteins. Furthermore, 3j commonly disrupted the stable and compact conformation of A₁ and A_2A AR, coupled with their active sites where Cα deviations were relatively high. Ligand mobility analysis also revealed that both compounds exhibited a similar motion pattern at the active site of the proteins relative to their optimal dual binding. We believe that findings from this study with significantly aid the structure-based design of highly selective dual-inhibitors of A₁ and A_2A AR.

Can We Exploit β-Lactamases Intrinsic Dynamics for Designing More Effective Inhibitors?

β-lactamases (BLs) represent the most frequent cause of antimicrobial resistance in Gram-negative bacteria. Despite the continuous efforts in the development of BL inhibitors (BLIs), new BLs able to hydrolyze the last developed antibiotics rapidly emerge. Moreover, the insurgence rate of effective mutations is far higher than the release of BLIs able to counteract them. This results in a shortage of antibiotics that is menacing the effective treating of infectious diseases. The situation is made even worse by the co-expression in bacteria of BLs with different mechanisms and hydrolysis spectra, and by the lack of inhibitors able to hit them all. Differently from other targets, BL flexibility has not been deeply exploited for drug design, possibly because of the small protein size, for their apparent rigidity and their high fold conservation. In this mini-review, we discuss the evidence for BL binding site dynamics being crucial for catalytic efficiency, mutation effect, and for the design of new inhibitors. Then, we report on identified allosteric sites in BLs and on possible allosteric inhibitors, as a strategy to overcome the frequent occurrence of mutations in BLs and the difficulty of competing efficaciously with substrates. Nevertheless, allosteric inhibitors could work synergistically with traditional inhibitors, increasing the chances of restoring bacterial susceptibility towards available antibiotics.

Ligand Binding Mechanism and Its Relationship with Conformational Changes in Adenine Riboswitch

Riboswitches are naturally occurring RNA aptamers that control the expression of essential bacterial genes by binding to specific small molecules. The binding with both high affinity and specificity induces conformational changes. Thus, riboswitches were proposed as a possible molecular target for developing antibiotics and chemical tools. The adenine riboswitch can bind not only to purine analogues but also to pyrimidine analogues. Here, long molecular dynamics (MD) simulations and molecular mechanics Poisson–Boltzmann surface area (MM-PBSA) computational methodologies were carried out to show the differences in the binding model and the conformational changes upon five ligands binding. The binding free energies of the guanine riboswitch aptamer with C74U mutation complexes were compared to the binding free energies of the adenine riboswitch (AR) aptamer complexes. The calculated results are in agreement with the experimental data. The differences for the same ligand binding to two different aptamers are related to the electrostatic contribution. Binding dynamical analysis suggests a flexible binding pocket for the pyrimidine ligand in comparison with the purine ligand. The 18 μs of MD simulations in total indicate that both ligand-unbound and ligand-bound aptamers transfer their conformation between open and closed states. The ligand binding obviously affects the conformational change. The conformational states of the aptamer are associated with the distance between the mass center of two key nucleotides (U51 and A52) and the mass center of the other two key nucleotides (C74 and C75). The results suggest that the dynamical character of the binding pocket would affect its biofunction. To design new ligands of the adenine riboswitch, it is recommended to consider the binding affinities of the ligand and the conformational change of the ligand binding pocket.

Applications of water molecules for analysis of macromolecule properties

Water molecules maintain proteins' structures, functions, stabilities and dynamics. They can occupy certain positions or pass quickly via a protein's interior. Regardless of their behaviour, water molecules can be used for the analysis of proteins' structural features and biochemical properties. Here, we present a list of several software programs that use the information provided by water molecules to: i) analyse protein structures and provide the optimal positions of water molecules for protein hydration, ii) identify high-occupancy water sites in order to analyse ligand binding modes, and iii) detect and describe tunnels and cavities. The analysis of water molecules' distribution and trajectories sheds a light on proteins' interactions with small molecules, on the dynamics of tunnels and cavities, on protein composition and also on the functionality, transportation network and location of functionally relevant residues. Finally, the correct placement of water molecules in protein crystal structures can significantly improve the reliability of molecular dynamics simulations.

The Role of the Ω-Loop in Regulation of the Catalytic Activity of TEM-Type β-Lactamases

Bacterial resistance to β-lactams, the most commonly used class of antibiotics, poses a global challenge. This resistance is caused by the production of bacterial enzymes that are termed β-lactamases (βLs). The evolution of serine-class A β-lactamases from penicillin-binding proteins (PBPs) is related to the formation of the Ω-loop at the entrance to the enzyme’s active site. In this loop, the Glu166 residue plays a key role in the two-step catalytic cycle of hydrolysis. This residue in TEM–type β-lactamases, together with Asn170, is involved in the formation of a hydrogen bonding network with a water molecule, leading to the deacylation of the acyl–enzyme complex and the hydrolysis of the β-lactam ring of the antibiotic. The activity exhibited by the Ω-loop is attributed to the positioning of its N-terminal residues near the catalytically important residues of the active site. The structure of the Ω-loop of TEM-type β-lactamases is characterized by low mutability, a stable topology, and structural flexibility. All of the revealed features of the Ω-loop, as well as the mechanisms related to its involvement in catalysis, make it a potential target for novel allosteric inhibitors of β-lactamases.

Dynamical Behavior of β-Lactamases and Penicillin- Binding Proteins in Different Functional States and Its Potential Role in Evolution

β-Lactamases are enzymes produced by bacteria to hydrolyze β-lactam-based antibiotics, and pose serious threat to public health through related antibiotic resistance. Class A β-lactamases are structurally and functionally related to penicillin-binding proteins (PBPs). Despite the extensive studies of the structures, catalytic mechanisms and dynamics of both β-lactamases and PBPs, the potentially different dynamical behaviors of these proteins in different functional states still remain elusive in general. In this study, four evolutionarily related proteins, including TEM-1 and TOHO-1 as class A β-lactamases, PBP-A and DD-transpeptidase as two PBPs, are subjected to molecular dynamics simulations and various analyses to characterize their dynamical behaviors in different functional states. Penicillin G and its ring opening product serve as common ligands for these four proteins of interest. The dynamic analyses of overall structures, the active sites with penicillin G, and three catalytically important residues commonly shared by all four proteins reveal unexpected cross similarities between Class A β-lactamases and PBPs. These findings shed light on both the hidden relations among dynamical behaviors of these proteins and the functional and evolutionary relations among class A β-lactamases and PBPs.

Exploring the Lapse in Druggability: Sequence Analysis, Structural Dynamics and Binding Site Characterization of K-RasG12C Variant, a Feasible Oncotherapeutics Target

BACKGROUND

The difficulty in druggability of K-Ras variant has presented a challenge in the treatment of cancer diseases associated with its dysfunctionality. Despite the identification of different binding sites, limited information exists in the literature about their characteristics. Therefore, identification, crossvalidation and characterization of its druggable sites would aid the design of chemical compounds that will arrest its dysfunctionality related oncogenesis.

OBJECTIVE

This study entails the identification, cross-validation and characterization of K-Ras G12C variant's binding sites for potential druggability, coupled with the elucidation of alterations in 3D conformations and dynamics.

METHOD

Molecular dynamics simulation was carried out on the inactive, the active and the hyperactive K-RasG12Cvariant using the amber software package. The SiteMap software was employed in identifying and characterizing the druggable binding sites while the validation of the binding sites was carried out with the SiteHound and MetaPocket servers.

RESULTS

Four druggable binding sites were identified, validated and characterized based on physicochemical attributes such as size, volume, degree of enclosure or exposure, degree of contact, hydrophobic/hydrophilic character, hydrophobic/hydrophilic balance and hydrogen-bonding features. Conformational studies also revealed that the K-Ras variant exhibited notable structural instability, increased flexibility and a strongly anticorrelated movement compared to the inactive and active wildtype forms.

CONCLUSION

The attributes of the characterized druggable sites will be useful in designing site-specific K-Ras inhibitors for the treatment of K-Ras variant associated cancer diseases.

Machine Learning Classification Model for Functional Binding Modes of TEM-1 β-Lactamase

TEM family of enzymes is one of the most commonly encountered β-lactamases groups with different catalytic capabilities against various antibiotics. Despite the studies investigating the catalytic mechanism of TEM β-lactamases, the binding modes of these enzymes against ligands in different functional catalytic states have been largely overlooked. But the binding modes may play a critical role in the function and even the evolution of these proteins. In this work, a newly developed machine learning analysis approach to the recognition of protein dynamics states was applied to compare the binding modes of TEM-1 β-lactamase with regard to penicillin in different catalytic states. While conventional analysis methods, including principal components analysis (PCA), could not differentiate TEM-1 in different binding modes, the application of a machine learning method led to excellent classification models differentiating these states. It was also revealed that both reactant/product states and apo/product states are more differentiable than the apo/reactant states. The feature importance generated by the training procedure of the machine learning model was utilized to evaluate the contribution from residues at active sites and in different secondary structures. Key active site residues, Ser70 and Ser130, play a critical role in differentiating reactant/product states, while other active site residues are more important for differentiating apo/product states. Overall, this study provides new insights into the different dynamical function states of TEM-1 and may open a new venue for β-lactamases functional and evolutional studies in general.

Influence of substrates and inhibitors on the structure of Klebsiella pneumoniae carbapenemase-2

The hydrolysis of last resort carbapenem antibiotics by Klebsiella pneumoniae carbapenemase-2 (KPC-2) presents a significant danger to global health. Combined with horizontal gene transfer, the emergence KPC-2 threatens to quickly expand carbapenemase activity to ever increasing numbers of pathogens. Our understanding of KPC-2 has greatly increased over the past decade thanks, in great part, to 20 crystal structures solved by groups around the world. These include apo KPC-2 structures, along with structures featuring a library of 10 different inhibitors representing diverse structural and functional classes. Herein we focus on cataloging the available KPC-2 structures and presenting a discussion of key aspects of each structure and important relationships between structures. Although the available structures do not provide information on dynamic motions with KPC-2, and the family of structures indicates small conformational changes across a wide array of bound inhibitors, substrates, and products, the structures provide a strong foundation for additional studies in the coming years to discover new KPC-2 inhibitors.

Impact statement

The work herein is important to the field as it provides a clear and succinct accounting of available KPC-2 structures. The work advances the field by collecting and analyzing differences and similarities across the available structures. This work features new analyses and interpretations of the existing structures which will impact the field in a positive way by making structural insights more widely available among the beta-lactamase community.

Probing Gallate-Mediated Selectivity and High-Affinity Binding of Epigallocatechin Gallate: a Way-Forward in the Design of Selective Inhibitors for Anti-apoptotic Bcl-2 Proteins

Selective inhibition is a key focus in the design of chemotherapeutic compounds that can abrogate the oncogenic activities of anti-apoptotic Bcl-2 proteins. Although recent efforts have led to the development of highly selective BH3 mimetics, setbacks such as toxicities have limited their use in cancer therapy. Epigallocatechingallate (EGCG) has been widely reported to selectively inhibit Bcl-2 and Bcl-xL compared to other green tea phenols due to its gallate group. Herein, we investigate the interaction dynamics of EGCG at the hydrophobic grooves of Bcl-2 and Bcl-xL and the consequential effects on their BH4 domains. Arg143 and Asp108 (Bcl-2), and Glu96 and Tyr195 (Bcl-xL) formed high-affinity hydrogen interactions with the gallate group while non-gallate groups of EGCG formed weak interactions. EGCG-bound proteins showed systemic perturbations of BH4 domains coupled with the burial of crucial surface-exposed residues such as Lys17 (Bcl-2) and Asp11 (Bcl-xL); hence, a distortion of non-canonical domain interactions. Interactions of gallate group of EGCG with key hydrophobic groove residues underlie EGCG selectivity while concurrent BH4 domain perturbations potentiate EGCG inhibitory activities. Findings will aid the optimization and design of selective inhibitors that could suppress anti-apoptotic activities of Bcl2-family proteins with minimal toxicities.

Spectroscopic analysis and docking simulation on the recognition and binding of TEM-1 β-lactamase with β-lactam antibiotics

The interaction between TEM-1 β-lactamase and antibiotics is very important in the hydrolysis of antibiotics. In the present study, the recognition and binding of TEM-1 β-lactamase with three β-lactam antibiotics, including penicillin G, cefalexin and cefoxitin, was investigated by fluorescence and ultraviolet-visible absorption spectra in combination with molecular docking in the temperature range of 278-288 K and under simulated physiological conditions. The results demonstrated that the fluorescence emissions of TEM-1 β-lactamase were extinguished by static quenching and the energy of TEM-1 β-lactamase was transferred in a non-radioactive manner. The binding of TEM-1 β-lactamase with the three antibiotics was a spontaneously exothermic process, with binding constants of 1.41×10⁷, 7.81×10⁶ and 5.43×10⁴ at 278 K. Furthermore, binding was driven by enthalpy change and the binding forces between them were mainly hydrogen bonding and Van der Waals forces. A TEM-1 β-lactamase only bound with one antibiotic at a time and the binding capacity between them was closely associated with the functional groups and flexibility in the antibiotics. In addition, a conformational change occurred in the TEM-1 β-lactamases when they bound with the three antibiotics and TEM-1 β-lactamase-antibiotic complexes were formed. The present study provided an insight into the recognition and binding of TEM-1 β-lactamase with β-lactam antibiotics, which may be helpful for designing a novel substrate for TEM-1 β-lactamase and developing novel antibiotics that are resistant to the enzyme.

Predicting allosteric mutants that increase activity of a major antibiotic resistance enzyme

Allosteric mutations increasing k_cat in a beta lactamase act by changing conformational ensembles of active-site residues identified by machine learning.

The CTX-M family of beta lactamases mediate broad-spectrum antibiotic resistance and are present in the majority of drug-resistant Gram-negative bacterial infections worldwide. Allosteric mutations that increase catalytic rates of these drug resistance enzymes have been identified in clinical isolates but are challenging to predict prospectively. We have used molecular dynamics simulations to predict allosteric mutants increasing CTX-M9 drug resistance, experimentally testing top mutants using multiple antibiotics. Purified enzymes show an increase in catalytic rate and efficiency, while mutant crystal structures show no detectable changes from wild-type CTX-M9. We hypothesize that increased drug resistance results from changes in the conformational ensemble of an acyl intermediate in hydrolysis. Machine-learning analyses on the three top mutants identify changes to the binding-pocket conformational ensemble by which these allosteric mutations transmit their effect. These findings show how molecular simulation can predict how allosteric mutations alter active-site conformational equilibria to increase catalytic rates and thus resistance against common clinically used antibiotics.

Binding of Solvent Molecules to a Protein Surface in Binary Mixtures Follows a Competitive Langmuir Model

The binding of solvent molecules to a protein surface was modeled by molecular dynamics simulations of of Candida antarctica (C. antarctica) lipase B in binary mixtures of water, methanol, and toluene. Two models were analyzed: a competitive Langmuir model which assumes identical solvent binding sites with a different affinity toward water (KWat), methanol (KMet), and toluene (KTol) and a competitive Langmuir model with an additional interaction between free water and already bound water (KWatWat). The numbers of protein-bound molecules of both components of a binary mixture were determined for different compositions as a function of their thermodynamic activities in the bulk phase, and the binding constants were simultaneously fitted to the six binding curves (two components of three different mixtures). For both Langmuir models, the values of KWat, KMet, and KTol were highly correlated. The highest binding affinity was found for methanol, which was almost 4-fold higher than the binding affinities of water and toluene (KMet ≫ KWat ≈ KTol). Binding of water was dominated by the water-water interaction (KWatWat). Even for the three protein surface patches of highest water affinity, the binding affinity of methanol was 2-fold higher than water and 8-fold higher than toluene (KMet > KWat > KTol). The Langmuir model provides insights into the protein destabilizing mechanism of methanol which has a high binding affinity toward the protein surface. Thus, destabilizing solvents compete with intraprotein interactions and disrupt the tertiary structure. In contrast, benign solvents such as water or toluene have a low affinity toward the protein surface. Water is a special solvent: only few water molecules bind directly to the protein; most water molecules bind to already bound water molecules thus forming water patches. A quantitative mechanistic model of protein-solvent interactions that includes competition and miscibility of the components contributes a robust basis for solvent and protein engineering.

Variability of the Cyclin-Dependent Kinase 2 Flexibility Without Significant Change in the Initial Conformation of the Protein or Its Environment; a Computational Study

BACKGROUND

Protein flexibility, which has been referred as a dynamic behavior has various roles in proteins' functions. Furthermore, for some developed tools in bioinformatics, such as protein-protein docking software, considering the protein flexibility, causes a higher degree of accuracy. Through undertaking the present work, we have accomplished the quantification plus analysis of the variations in the human Cyclin Dependent Kinase 2 (hCDK2) protein flexibility without affecting a significant change in its initial environment or the protein per se.

OBJECTIVES

The main goal of the present research was to calculate variations in the flexibility for each residue of the hCDK2, analysis of their flexibility variations through clustering, and to investigate the functional aspects of the residues with high flexibility variations.

MATERIALS AND METHODS

Using Gromacs package (version 4.5.4), three independent molecular dynamics (MD) simulations of the hCDK2 protein (PDB ID: 1HCL) was accomplished with no significant changes in their initial environments, structures, or conformations, followed by Root Mean Square Fluctuations (RMSF) calculation of these MD trajectories. The amount of variations in these three curves of RMSF was calculated using two formulas.

RESULTS

More than 50% of the variation in the flexibility (the distance between the maximum and the minimum amount of the RMSF) was found at the region of Val-154. As well, there are other major flexibility fluctuations in other residues. These residues were mostly positioned in the vicinity of the functional residues. The subsequent works were done, as followed by clustering all hCDK2 residues into four groups considering the amount of their variability with respect to flexibility and their position in the RMSF curves.

CONCLUSIONS

This work has introduced a new class of flexibility aspect of the proteins' residues. It could also help designing and engineering proteins, with introducing a new dynamic aspect of hCDK2, and accordingly, for the other similar globular proteins. In addition, it could provide a better computational calculation of the protein flexibility, which is, especially important in the comparative studies of the proteins' flexibility.

Flexibility Correlation between Active Site Regions Is Conserved across Four AmpC β-Lactamase Enzymes

β-lactamases are bacterial enzymes that confer resistance to β-lactam antibiotics, such as penicillins and cephalosporins. There are four classes of β-lactamase enzymes, each with characteristic sequence and structure properties. Enzymes from class A are the most common and have been well characterized across the family; however, less is known about how physicochemical properties vary across the C and D families. In this report, we compare the dynamical properties of four AmpC (class C) β-lactamases using our distance constraint model (DCM). The DCM reliably predicts thermodynamic and mechanical properties in an integrated way. As a consequence, quantitative stability/flexibility relationships (QSFR) can be determined and compared across the whole family. The DCM calculates a large number of QSFR metrics. Perhaps the most useful is the flexibility index (FI), which quantifies flexibility along the enzyme backbone. As typically observed in other systems, FI is well conserved across the four AmpC enzymes. Cooperativity correlation (CC), which quantifies intramolecular couplings within structure, is rarely conserved across protein families; however, it is in AmpC. In particular, the bulk of each structure is composed of a large rigid cluster, punctuated by three flexibly correlated regions located at the active site. These regions include several catalytic residues and the Ω-loop. This evolutionary conservation combined with active their site location strongly suggests that these coupled dynamical modes are important for proper functioning of the enzyme.

Maintenance of native-like protein dynamics may not be required for engineering functional proteins

Proteins are dynamic systems, and understanding dynamics is critical for fully understanding protein function. Therefore, the question of whether laboratory engineering has an impact on protein dynamics is of general interest. Here, we demonstrate that two homologous, naturally evolved enzymes with high degrees of structural and functional conservation also exhibit conserved dynamics. Their similar set of slow timescale dynamics is highly restricted, consistent with evolutionary conservation of a functionally important feature. However, we also show that dynamics of a laboratory-engineered chimeric enzyme obtained by recombination of the two homologs exhibits striking difference on the millisecond timescale, despite function and high-resolution crystal structure (1.05 Å) being conserved. The laboratory-engineered chimera is thus functionally tolerant to modified dynamics on the timescale of catalytic turnover. Tolerance to dynamic variation implies that maintenance of native-like protein dynamics may not be required when engineering functional proteins.

Variations within class-A β-lactamase physiochemical properties reflect evolutionary and environmental patterns, but not antibiotic specificity

The bacterial enzyme β-lactamase hydrolyzes the β-lactam ring of penicillin and chemically related antibiotics, rendering them ineffective. Due to rampant antibiotic overuse, the enzyme is evolving new resistance activities at an alarming rate. Related, the enzyme's global physiochemical properties exhibit various amounts of conservation and variability across the family. To that end, we characterize the extent of property conservation within twelve different class-A β-lactamases, and conclusively establish that the systematic variations therein parallel their evolutionary history. Large and systematic differences within electrostatic potential maps and pairwise residue-to-residue couplings are observed across the protein, which robustly reflect phylogenetic outgroups. Other properties are more conserved (such as residue pKa values, electrostatic networks, and backbone flexibility), yet they also have systematic variations that parallel the phylogeny in a statistically significant way. Similarly, the above properties also parallel the environmental condition of the bacteria they are from in a statistically significant way. However, it is interesting and surprising that the only one of the global properties (protein charge) parallels the functional specificity patterns; meaning antibiotic resistance activities are not significantly constraining the global physiochemical properties. Rather, extended spectrum activities can emerge from the background of nearly any set of electrostatic and dynamic properties.

Chimeric β-lactamases: global conservation of parental function and fast time-scale dynamics with increased slow motions

Enzyme engineering has been facilitated by recombination of close homologues, followed by functional screening. In one such effort, chimeras of two class-A β-lactamases – TEM-1 and PSE-4 – were created according to structure-guided protein recombination and selected for their capacity to promote bacterial proliferation in the presence of ampicillin (Voigt et al., Nat. Struct. Biol. 2002 9:553). To provide a more detailed assessment of the effects of protein recombination on the structure and function of the resulting chimeric enzymes, we characterized a series of functional TEM-1/PSE-4 chimeras possessing between 17 and 92 substitutions relative to TEM-1 β-lactamase. Circular dichroism and thermal scanning fluorimetry revealed that the chimeras were generally well folded. Despite harbouring important sequence variation relative to either of the two ‘parental’ β-lactamases, the chimeric β-lactamases displayed substrate recognition spectra and reactivity similar to their most closely-related parent. To gain further insight into the changes induced by chimerization, the chimera with 17 substitutions was investigated by NMR spin relaxation. While high order was conserved on the ps-ns timescale, a hallmark of class A β-lactamases, evidence of additional slow motions on the µs-ms timescale was extracted from model-free calculations. This is consistent with the greater number of resonances that could not be assigned in this chimera relative to the parental β-lactamases, and is consistent with this well-folded and functional chimeric β-lactamase displaying increased slow time-scale motions.

Molecular dynamics of class A β-lactamases-effects of substrate binding

The effects of substrate binding on class A β-lactamase dynamics were studied using molecular dynamics simulations of two model enzymes; 40 100-ns trajectories of the free and substrate-bound forms of TEM-1 (with benzylpenicillin) and PSE-4 (with carbenicillin) were recorded (totaling 4.0 μs). Substrates were parameterized with the CHARMM General Force Field. In both enzymes, the Ω loop exhibits a marked flexibility increase upon substrate binding, supporting the hypothesis of substrate gating. However, specific interactions that are formed or broken in the Ω loop upon binding differ between the two enzymes: dynamics are conserved, but not specific interactions. Substrate binding also has a global structuring effect on TEM-1, but not on PSE-4. Changes in TEM-1's normal modes show long-range effects of substrate binding on enzyme dynamics. Hydrogen bonds observed in the active site are mostly preserved upon substrate binding, and new, transient interactions are also formed. Agreement between NMR relaxation parameters and our theoretical results highlights the dynamic duality of class A β-lactamases: enzymes that are highly structured on the ps-ns timescale, with important flexibility on the μs-ms timescale in regions such as the Ω loop.

Interfacial water molecules in SH3 interactions: a revised paradigm for polyproline recognition

In spite of its biomedical relevance, polyproline recognition is still not fully understood. The disagreement between the current description of SH3 (Src homology 3) complexes and their thermodynamic behaviour calls for a revision of the SH3-binding paradigm. Recently, Abl-SH3 was demonstrated to recognize its ligands by a dual binding mechanism involving a robust network of water-mediated hydrogen bonds that complements the canonical hydrophobic interactions. The systematic analysis of the SH3 structural database in the present study reveals that this dual binding mode is universal to SH3 domains. Tightly bound buried-interfacial water molecules were found in all SH3 complexes studied mediating the interaction between the peptide ligand and the domain. Moreover, structural waters were also identified in a high percentage of the free SH3 domains. A detailed analysis of the pattern of water-mediated interactions enabled the identification of conserved hydration sites in the polyproline-recognition region and the establishment of relationships between hydration profiles and the sequence of both ligands and SH3 domains. Water-mediated interactions were also systematically observed in WW (protein–protein interaction domain containing two conserved tryptophan residues), UEV (ubiquitin-conjugating enzyme E2 variant) and EVH-1 [Ena/VASP (vasodilator-stimulated phosphoprotein) homology 1] structures. The results of the present study clearly indicate that the current description of proline-rich sequence recognition by protein–protein interaction modules is incomplete and insufficient for a correct understanding of these systems. A new binding paradigm is required that includes interfacial water molecules as relevant elements in polyproline recognition.

The evolution of cefotaximase activity in the TEM β-lactamase

The development of a molecular-level understanding of drug resistance through β-lactamase is critical not only in designing newer-generation antibacterial agents but also in providing insight into the evolutionary mechanisms of enzymes in general. In the present study, we have evaluated the effect of four drug resistance mutations (A42G, E104K, G238S, and M182T) on the cefotaximase activity of the TEM-1 β-lactamase. Using computational methods, including docking and molecular mechanics calculations, we have been able to correctly identify the relative order of catalytic activities associated with these four single point mutants. Further analyses suggest that the changes in catalytic efficiency for mutant enzymes are correlated to structural changes within the binding site. Based on the energetic and structural analyses of the wild-type and mutant enzymes, structural rearrangement is suggested as a mechanism of evolution of drug resistance through TEM β-lactamase. The present study not only provides molecular-level insight into the effect of four drug resistance mutations on the structure and function of the TEM β-lactamase but also establishes a foundation for a future molecular-level analysis of complete evolutionary trajectory for this class of enzymes.