The electrostatic profile of consecutive Cβ atoms applied to protein structure quality assessment.

The structure of a protein provides insight into its physiological interactions with other components of the cellular soup. Methods that predict putative structures from sequences typically yield multiple, closely-ranked possibilities. A critical component in the process is the model quality assessing program (MQAP), which selects the best candidate from this pool of structures. Here, we present a novel MQAP based on the physical properties of sidechain atoms. We propose a method for assessing the quality of protein structures based on the electrostatic potential difference (EPD) of Cβ atoms in consecutive residues. We demonstrate that the EPDs of Cβ atoms on consecutive residues provide unique signatures of the amino acid types. The EPD of Cβ atoms are learnt from a set of 1000 non-homologous protein structures with a resolution cuto of 1.6 Å obtained from the PISCES database. Based on the Boltzmann hypothesis that lower energy conformations are proportionately sampled more, and on Annsen's thermodynamic hypothesis that the native structure of a protein is the minimum free energy state, we hypothesize that the deviation of observed EPD values from the mean values obtained in the learning phase is minimized in the native structure. We achieved an average specificity of 0.91, 0.94 and 0.93 on hg_structal, 4state_reduced and ig_structal decoy sets, respectively, taken from the Decoys `R' Us database. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134

Numerical investigation into discontinuity-induced bifurcations in an aeroelastic system with coupled non-smooth nonlinearities

peer reviewedThe present study focuses on investigating the bifurcation characteristics of a pitch–plunge aeroelastic system possessing coupled non-smooth nonlinearities, both in structural and aerodynamic fronts. To this end, a freeplay nonlinearity is considered in the stiffness of the pitch degree-of-freedom. The effects of dynamic stall arising due to large instantaneous angles-of-attackare incorporated using the semi-empirical Leishman–Beddoes aerodynamic model. A systematic response analysis is carried out to discern the bifurcation characteristics of the aeroelastic system considering the airspeed as the system parameter. At low airspeeds, a series of dynamical transitions, including aperiodic responses, occur predominantly due to the structural freeplay nonlinearity while the flow remains attached to the surface of the wing. However, beyond a critical value of airspeed, the system response is dominated by high amplitude pitch-dominated limit-cycle oscillations, which can be attributed to stall flutter. It is demonstrated that the freeplay gap plays a key role in combining the effects of structural and aerodynamic nonlinearities. At higher values of the freeplay gap, interesting discontinuity-induced bifurcation scenarios, such as grazing and boundary equilibrium bifurcations arise due to coupled nonlinear interactions, which can significantly impact the safety of the aeroelastic system

Protein structure quality assessment based on the distance profiles of consecutive backbone Cα atoms

Predicting the three dimensional native state structure of a protein from its primary sequence is an unsolved grand challenge in molecular biology. Two main computational approaches have evolved to obtain the structure from the protein sequence - ab initio/de novo methods and template-based modeling - both of which typically generate multiple possible native state structures. Model quality assessment programs (MQAP) validate these predicted structures in order to identify the correct native state structure. Here, we propose a MQAP for assessing the quality of protein structures based on the distances of consecutive Cα atoms. We hypothesize that the root-mean-square deviation of the distance of consecutive Cα (RDCC) atoms from the ideal value of 3.8 Å, derived from a statistical analysis of high quality protein structures (top100H database), is minimized in native structures. Based on tests with the top100H set, we propose a RDCC cutoff value of 0.012 Å, above which a structure can be filtered out as a non-native structure. We applied the RDCC discriminator on decoy sets from the Decoys 'R' Us database to show that the native structures in all decoy sets tested have RDCC below the 0.012 Å cutoff. While most decoy sets were either indistinguishable using this discriminator or had very few violations, all the decoy structures in the fisa decoy set were discriminated by applying the RDCC criterion. This highlights the physical non-viability of the fisa decoy set, and possible issues in benchmarking other methods using this set. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134

The electrostatic profile of consecutive Cβ atoms applied to protein structure quality assessment

The structure of a protein provides insight into its physiological interactions with other components of the cellular soup. Methods that predict putative structures from sequences typically yield multiple, closely-ranked possibilities. A critical component in the process is the model quality assessing program (MQAP), which selects the best candidate from this pool of structures. Here, we present a novel MQAP based on the physical properties of sidechain atoms. We propose a method for assessing the quality of protein structures based on the electrostatic potential difference (EPD) of Cβ atoms in consecutive residues. We demonstrate that the EPDs of Cβ atoms on consecutive residues provide unique signatures of the amino acid types. The EPD of Cβ atoms are learnt from a set of 1000 non-homologous protein structures with a resolution cuto of 1.6 Å obtained from the PISCES database. Based on the Boltzmann hypothesis that lower energy conformations are proportionately sampled more, and on Annsen's thermodynamic hypothesis that the native structure of a protein is the minimum free energy state, we hypothesize that the deviation of observed EPD values from the mean values obtained in the learning phase is minimized in the native structure. We achieved an average specificity of 0.91, 0.94 and 0.93 on hg_structal, 4state_reduced and ig_structal decoy sets, respectively, taken from the Decoys 'R' Us database. The source code and manual is made available at https://github.com/sanchak/mqap and permanently available on 10.5281/zenodo.7134

The dipeptidyl peptidase IV inhibitors vildagliptin and K-579 inhibit a phospholipase C: a case of promiscuous scaffolds in proteins.

The long term side effects of any newly introduced drug is a subject of intense research, and often raging controversies. One such example is the dipeptidyl peptidase-IV (DPP4) inhibitor used for treating type 2 diabetes, which is inconclusively implicated in increased susceptibility to acute pancreatitis. Previously, based on a computational analysis of the spatial and electrostatic properties of active site residues, we have demonstrated that phosphoinositide-specific phospholipase C (PI-PLC) from Bacillus cereus is a prolyl peptidase using in vivo experiments. In the current work, we first report the inhibition of the native activity of PI-PLC by two DPP4 inhibitors - vildagliptin (LAF-237) and K-579. While vildagliptin inhibited PI-PLC at micromolar concentrations, K-579 was a potent inhibitor even at nanomolar concentrations. Subsequently, we queried a comprehensive, non-redundant set of 5000 human proteins (50% similarity cutoff) with known structures using serine protease (SPASE) motifs derived from trypsin and DPP4. A pancreatic lipase and a gastric lipase are among the proteins that are identified as proteins having promiscuous SPASE scaffolds that could interact with DPP4 inhibitors. The presence of such scaffolds in human lipases is expected since they share the same catalytic mechanism with PI-PLC. However our methodology also detects other proteins, often with a completely different enzymatic mechanism, that have significantly congruent domains with the SPASE motifs. The reported elevated levels of serum lipase, although contested, could be rationalized by inhibition of lipases reported here. In an effort to further our understanding of the spatial and electrostatic basis of DPP4 inhibitors, we have also done a comprehensive analysis of all 76 known DPP4 structures liganded to inhibitors till date. Also, the methodology presented here can be easily adopted for other drugs, and provide the first line of filtering in the identification of pathways that might be inadvertently affected due to promiscuous scaffolds in proteins

The dipeptidyl peptidase IV inhibitors vildagliptin and K-579 inhibit a phospholipase C: a case of promiscuous scaffolds in proteins

The long term side effects of any newly introduced drug is a subject of intense research, and often raging controversies. One such example is the dipeptidyl peptidase-IV (DPP4) inhibitor used for treating type 2 diabetes, which is inconclusively implicated in increased susceptibility to acute pancreatitis. Previously, based on a computational analysis of the spatial and electrostatic properties of active site residues, we have demonstrated that phosphoinositide-specific phospholipase C (PI-PLC) from Bacillus cereus is a prolyl peptidase using in vivo experiments. In the current work, we first report the inhibition of the native activity of PI-PLC by two DPP4 inhibitors - vildagliptin (LAF-237) and K-579. While vildagliptin inhibited PI-PLC at micromolar concentrations, K-579 was a potent inhibitor even at nanomolar concentrations. Subsequently, we queried a comprehensive, non-redundant set of 5000 human proteins (50% similarity cutoff) with known structures using serine protease (SPASE) motifs derived from trypsin and DPP4. A pancreatic lipase and a gastric lipase are among the proteins that are identified as proteins having promiscuous SPASE scaffolds that could interact with DPP4 inhibitors. The presence of such scaffolds in human lipases is expected since they share the same catalytic mechanism with PI-PLC. However our methodology also detects other proteins, often with a completely different enzymatic mechanism, that have significantly congruent domains with the SPASE motifs. The reported elevated levels of serum lipase, although contested, could be rationalized by inhibition of lipases reported here. In an effort to further our understanding of the spatial and electrostatic basis of DPP4 inhibitors, we have also done a comprehensive analysis of all 76 known DPP4 structures liganded to inhibitors till date. Also, the methodology presented here can be easily adopted for other drugs, and provide the first line of filtering in the identification of pathways that might be inadvertently affected due to promiscuous scaffolds in proteins