Building polarization into protein-inhibitor binding dynamics in rational drug design for rheumatoid arthritis.

Standard force field-based simulations to accomplish structure-based evaluations of lead molecules is a powerful tool. Combining protein fragmentation into tractable sub-systems with continuum solvation method is envisaged to enable quantum mechanics-based electronic structure calculations of macromolecules in their realistic environment. This along with incorporation of many-body polarization effect in molecular dynamics simulations may augment an accurate description of electrostatics of protein-inhibitor systems for effective drug design. Rheumatoid arthritis (RA) is a complex autoimmune disorder plagued by the ceiling effect of current targeted therapies, encouraging identification of new druggable targets and corresponding drug design to tackle the refractory form of disease. In this study, polarization-inclusive force field approach has been used to model protein solvation and ligand binding for 'Mitogen-activated protein kinase' (MAP3K8), a regulatory node of notable pharmacological relevance in RA synovial biology. For MAP3K8 inhibitors belonging to different scaffold series, the calculations illustrated differential electrostatic contribution to their relative binding affinities and successfully explained examples from available structure-activity relationship studies. Results from this study exemplified i) the advantage of this approach in reliably ranking inhibitors having close nanomolar range activities for the same target; and ii) its prospective application in lead molecule identification aiding drug discovery efforts in RA.Communicated by Ramaswamy H. Sarma.

Genetic variations of TLRs and their association with HPV/EBV, co-infection along with nicotine exposure in the development of premalignant/malignant lesions of the oral cavity in Indian population.

BACKGROUND: Despite being most preventable malignancies associated with smoked and smokeless tobacco products, squamous cell carcinoma of oral cavity is one of the most common malignancy in India. The aim of the present study was to evaluate the role of TLRs in oral pre-cancerous, cancerous cases and their genotypic correlation with HPV/EBV, co-infection & lifestyle habits in Indian population. METHODS: The present study was conducted on 300 subjects (100 OSCC, 50 pre-cancer & 150 controls). The amplification of TLRs gene and HPV/EBV co-infection was assessed by Nested PCR, PCR-RFLP and further confirmation by direct sequencing. RESULTS: The TLR 9(-1486 T/C), revealed that the TT vs. CT + CC genotype had a ˜5-fold increased risk for the development of pre-cancerous lesions as compared to controls (p = 0.0001). Further analysis showed that the risk of cancer was extremely pronounced in HPV/EBV, co-infection (p = 0.0141), implicating the possible interaction between TLR 9(-1486T/C) genotype and HPV infection in increasing cancer/pre-cancer risk. The 'G' allele of TLR 4(+896A/G) was also a higher risk of developing pre-cancerous lesions with 4.5 fold and statistically significant (p = 0.0001). The genotypic association of TLR 9(-1486T/C) in OSMF cases showed ˜8 fold increased risk and TLR 4(+896A/G) showed fourteen fold higher risk for leukoplakia (p < 0.0001, OR = 14.000). CONCLUSION: Genetic polymorphism of TLR 9(-1486 T/C) and TLR 4(+896A/G) may influence the effects of HPV/EBV, co-infection and play the significant role in development of the disease. The significance of these TLRs seemed to be enhanced by tobacco chewing and smoking habits also, which act as an important etiological risk factor for OSCC.

Genetic variations of IL-10: Identification of novel variations and evaluation of the impact of the SNPs/haplotype in the promoter region with the progression of Oral Squamous Cell Carcinoma in Indian population.

The correlation of interleukin 10 (IL-10) with the outbreak and progression of cancer has been well established as it contributes to tumor immune evasion. Convincing number of evidences has been accumulated to reflect the critical correlation between IL-10 polymorphism and tumorogenesis. Several polymorphic sites at promoter regions have been reported to be associated with cancer susceptibility. The purpose of this study was to examine the effect of modulated genotypes in the promoter region of IL-10 gene with life-style habits in oral squamous cell carcinoma (OSCC) in the Indian population. A total of 300 subjects (100 OSCC, 50 precancer and 150 healthy controls) were recruited in this study. The IL-10 promoter region was amplified in 14 overlapping fragments by PCR and further screened through the high throughput technique of denaturing high-performance liquid chromatography (dHPLC) followed by sequencing. We identified three novel variations at positions (-924, -1045 & -1066); we also found some known SNPs (-592C/A, -657G/A, -851G/A, -819C/T, -1082A/G). The identified novel variations were submitted to the NCBI Gene Bank (accession numbers KT153594, KT291742 and KT291743). We also noticed a significant association of polymorphisms (-592C/A, -819C/T and -1082A/G) individually as well as in combination (haplotypes) along with lifestyle habits for the risk of oral carcinoma (p<0.0001). We have reported three novel SNPs in the Indian population for the first time, and these SNPs may be associated with OSCC. Besides, we showed the first evidence of IL-10 haplotypes, i.e., CCG and CTG, may act as a biomarker for early detection of oral pre-cancerous/cancerous lesions or treatment management of oral carcinoma.

Oral Screening for Pre-cancerous Lesions Among Areca-nut Chewing Population from Rural India.

PURPOSE: To detect early oral premalignant lesions (OPLs) in a rural population chewing tobacco-free areca nut preparations, determine their awareness level of oral cancer and educate them about maintaining good oral health. MATERIALS AND METHODS: A total of 2175 18- to 65-year-old areca nut chewers (male:female ratio 2.5:1), without a history of consuming tobacco in any form, from the villages of two districts of the West Bengal state of India were screened clinically through oral examination for suspected OPLs. A pre-designed questionnaire was employed to record demographic data, information on tobacco-free areca-nut chewing habit and knowledge about oral diseases. Education on oral health was provided through distribution of printed leaflets, display of banner/posters and a public-announcement system. RESULTS: Chewing areca nut in the form of betel quid was more popular (90.7%) than chewing areca nut alone (9%) or tobacco-free packaged areca nut preparation sold as 'pan masala' (0.3%). OPLs were detected in 7.3% of the subjects, more among the males. An increasing incidence of OPLs could be observed with an increase in age as well as with duration and frequency of areca-nut chewing, while decreasing incidence was observed with an increase in educational level. Oral submucous fibrosis showed the highest prevalence (2.7%) among the various OPLs detected. CONCLUSIONS: Tobacco-free areca-nut chewing is an independent risk factor for the development of OPL and a large rural population still practices such high risk behaviour. In rural areas with limited health care resources, screening by visual oral examination involving minimum cost may prove useful to reduce oral cancer mortality.

Bhageerath: an energy based web enabled computer software suite for limiting the search space of tertiary structures of small globular proteins.

We describe here an energy based computer software suite for narrowing down the search space of tertiary structures of small globular proteins. The protocol comprises eight different computational modules that form an automated pipeline. It combines physics based potentials with biophysical filters to arrive at 10 plausible candidate structures starting from sequence and secondary structure information. The methodology has been validated here on 50 small globular proteins consisting of 2-3 helices and strands with known tertiary structures. For each of these proteins, a structure within 3-6 A RMSD (root mean square deviation) of the native has been obtained in the 10 lowest energy structures. The protocol has been web enabled and is accessible at http://www.scfbio-iitd.res.in/bhageerath.

Protein structure evaluation using an all-atom energy based empirical scoring function.

Arriving at the native conformation of a polypeptide chain characterized by minimum most free energy is a problem of long standing interest in protein structure prediction endeavors. Owing to the computational requirements in developing free energy estimates, scoring functions--energy based or statistical--have received considerable renewed attention in recent years for distinguishing native structures of proteins from non-native like structures. Several cleverly designed decoy sets, CASP (Critical Assessment of Techniques for Protein Structure Prediction) structures and homology based internet accessible three dimensional model builders are now available for validating the scoring functions. We describe here an all-atom energy based empirical scoring function and examine its performance on a wide series of publicly available decoys. Barring two protein sequences where native structure is ranked second and seventh, native is identified as the lowest energy structure in 67 protein sequences from among 61,659 decoys belonging to 12 different decoy sets. We further illustrate a potential application of the scoring function in bracketing native-like structures of two small mixed alpha/beta globular proteins starting from sequence and secondary structural information. The scoring function has been web enabled at www.scfbio-iitd.res.in/utility/proteomics/energy.jsp.

A computational pathway for bracketing native-like structures fo small alpha helical globular proteins.

Impressive advances in the applications of bioinformatics for protein structure prediction coupled with growing structural databases on one hand and the insurmountable time-scale problem with ab initio computational methods on the other continue to raise doubts whether a computational solution to the protein folding problem--categorized as an NP-hard problem--is within reach in the near future. Combining some specially designed biophysical filters and vector algebra tools with ab initio methods, we present here a promising computational pathway for bracketing native-like structures of small alpha helical globular proteins departing from secondary structural information. The automated protocol is initiated by generating multiple structures around the loops between secondary structural elements. A set of knowledge-based biophysical filters namely persistence length and radius of gyration, developed and calibrated on approximately 1000 globular proteins, is introduced to screen the trial structures to filter out improbable candidates for the native and reduce the size of the library of probable structures. The ensemble so generated encompasses a few structures with native-like topology. Monte Carlo optimizations of the loop dihedrals are then carried out to remove steric clashes. The resultant structures are energy minimized and ranked according to a scoring function tested previously on a series of decoy sets vis-a-vis their corresponding natives. We find that the 100 lowest energy structures culled from the ensemble of energy optimized trial structures comprise at least a few to within 3-5 angstroms of the native. Thus the formidable "needle in a haystack" problem is narrowed down to finding an optimal solution amongst a computationally tractable number of alternatives. Encouraging results obtained on twelve small alpha helical globular proteins with the above outlined pathway are presented and discussed.