Indicators used by clinicians to assess pain in the brain injured

INTRODUCTION. The assessment of pain in critically ill brain-injured patients is challenging for health professionals. In addition to be unable to self-report, the confused and stereotyped behaviors of these patients are likely to alter their ''normal'' pain responses. Therefore, the pain indicators observed in the general critically ill population may not be appropriate. OBJECTIVES. To identify behavioral and physiological indicators used by clinicians to assess pain in critically ill brain-injured patients who are unable to self-report. METHODS.Amixed-method design was used with the first step being the combination of the results of an integrative literature review with the results of nominal groups of 12 nurses and four physicians. The second step involved a web-based survey to establish content validity. Fourteen experts (clinicians and academics) from three French speaking European countries rated the relevance of each indicator. A content validity index (CVI) was computed for each indicator (I-CVI) and for each category (S-CVI). RESULTS. The first step generated 52 indicators. These indicators were classified into six categories: facial expressions, position/movement, muscle tension, vocalization, compliance with ventilator, and physiological indicators. In the second step, the agreement between raters was high with an Intraclass Correlation Coefficient of 0.88 (95% CI 0.83-0.92). The I-CVIs ranged from 0.07 to 1. Indicators with an I-CVI below 0.5 (n = 12) were not retained, resulting in a final list of 30 indicators. The CVI for this final list was 0.75 with categories ranging from 0.67 (compliance with ventilation) to 0.87 (vocalization). CONCLUSIONS. This process identified specific pain indicators for critically ill braininjured patients. Further evaluation is in progress to test the validity and relevance of these indicators in the clinical setting.

New combined serum creatinine and cystatin C quadratic formula for GFR assessment in children.

BACKGROUND AND OBJECTIVES: The estimated GFR (eGFR) is important in clinical practice. To find the best formula for eGFR, this study assessed the best model of correlation between sinistrin clearance (iGFR) and the solely or combined cystatin C (CysC)- and serum creatinine (SCreat)-derived models. It also evaluated the accuracy of the combined Schwartz formula across all GFR levels. DESIGN, SETTING, PARTICIPANTS, & MEASUREMENTS: Two hundred thirty-eight iGFRs performed between January 2012 and April 2013 for 238 children were analyzed. Regression techniques were used to fit the different equations used for eGFR (i.e., logarithmic, inverse, linear, and quadratic). The performance of each model was evaluated using the Cohen κ correlation coefficient and the percentage reaching 30% accuracy was calculated. RESULTS: The best model of correlation between iGFRs and CysC is linear; however, it presents a low κ coefficient (0.24) and is far below the Kidney Disease Outcomes Quality Initiative targets to be validated, with only 84% of eGFRs reaching accuracy of 30%. SCreat and iGFRs showed the best correlation in a fitted quadratic model with a κ coefficient of 0.53 and 93% accuracy. Adding CysC significantly (P<0.001) increased the κ coefficient to 0.56 and the quadratic model accuracy to 97%. Therefore, a combined SCreat and CysC quadratic formula was derived and internally validated using the cross-validation technique. This quadratic formula significantly outperformed the combined Schwartz formula, which was biased for an iGFR≥91 ml/min per 1.73 m(2). CONCLUSIONS: This study allowed deriving a new combined SCreat and CysC quadratic formula that could replace the combined Schwartz formula, which is accurate only for children with moderate chronic kidney disease.

A filtering method to correct time-lapse 3D ERT data and improve imaging of natural aquifer dynamics

We have developed a processing methodology that allows crosshole ERT (electrical resistivity tomography) monitoring data to be used to derive temporal fluctuations of groundwater electrical resistivity and thereby characterize the dynamics of groundwater in a gravel aquifer as it is infiltrated by river water. Temporal variations of the raw ERT apparent-resistivity data were mainly sensitive to the resistivity (salinity), temperature and height of the groundwater, with the relative contributions of these effects depending on the time and the electrode configuration. To resolve the changes in groundwater resistivity, we first expressed fluctuations of temperature-detrended apparent-resistivity data as linear superpositions of (i) time series of riverwater-resistivity variations convolved with suitable filter functions and (ii) linear and quadratic representations of river-water-height variations multiplied by appropriate sensitivity factors; river-water height was determined to be a reliable proxy for groundwater height. Individual filter functions and sensitivity factors were obtained for each electrode configuration via deconvolution using a one month calibration period and then the predicted contributions related to changes in water height were removed prior to inversion of the temperature-detrended apparent-resistivity data. Applications of the filter functions and sensitivity factors accurately predicted the apparent-resistivity variations (the correlation coefficient was 0.98). Furthermore, the filtered ERT monitoring data and resultant time-lapse resistivity models correlated closely with independently measured groundwater electrical resistivity monitoring data and only weakly with the groundwater-height fluctuations. The inversion results based on the filtered ERT data also showed significantly less inversion artefacts than the raw data inversions. We observed resistivity increases of up to 10% and the arrival time peaks in the time-lapse resistivity models matched those in the groundwater resistivity monitoring data.

No use for waist-for-height ratio in addition to body mass index to identify children with elevated blood pressure.

Abstract Background. In children, waist-for-height ratio (WHtR) has been proposed to identify subjects at higher risk of cardiovascular diseases. The utility of WHtR to identify children with elevated blood pressure (BP) is unclear. Design. Cross-sectional population-based study of schoolchildren. Methods. Weight, height, waist circumference and BP were measured in all sixth-grade schoolchildren of the canton de Vaud (Switzerland) in 2005/06. WHtR was computed as waist [cm]/height [cm]. Elevated BP was defined according to sex-, age- and height-specific US reference data. The area under the receiver operating characteristic curve (AUC) statistic was computed to compare the ability of body mass index (BMI) z-score and WHtR, alone or in combination, to identify children with elevated BP. Results. 5207 children participated (76% response) [2621 boys, 2586 girls; mean (± SD) age, 12.3 ± 0.5 years; range: 10.1-14.9]. The prevalence of elevated BP was 11%. Mean WHtR was 0.44 ± 0.05 (range: 0.29- 0.77) and 11% had high WHtR (> 0.5). BMI z-score and WHtR were strongly correlated (Spearman correlation coefficient r = 0.76). Both indices were positively associated with elevated BP. AUCs for elevated BP was relatively low for BMI z-score (0.62) or for WHtR (0.62), and was not substantially improved when both indices were considered together (0.63). Conclusions. The ability of BMI z-score or WHtR to identify children aged 10-14 with elevated BP was weak. Adding WHtR did not confer additional discriminative power to BMI alone. These findings do not support the measurement of WHtR in addition to BMI to identify children with elevated BP.

Relationships between grain yield and accumulation of biomass, nitrogen and phosphorus in common bean cultivars

Shoot biomass is considered a relevant component for crop yield, but relationships between biological productivity and grain yield in legume crops are usually difficult to establish. Two field experiments were carried out to investigate the relationships between grain yield, biomass production and N and P accumulation at reproductive stages of common bean (Phaseolus vulgaris) cultivars. Nine and 18 cultivars were grown on 16 m² plots in 1998 and 1999, respectively, with four replications. Crop biomass was sampled at four growth stages (flowering R6, pod setting R7, beginning of pod filling R8, and mid-pod filling R8.5), grain yield was measured at maturity, and N and P concentrations were determined in plant tissues. In both years, bean cultivars differed in grain yield, in root mass at R6 and R7 stages, and in shoot mass at R6 and R8.5, whereas at R7 and R8 differences in shoot mass were significant in 1998 only. In both years, grain yield did not correlate with shoot mass at R6 and R7 and with root mass at R6. Grain yield correlated with shoot mass at R8 in 1999 but not in 1998, with shoot mass at R8.5 and with root mass at R7 in both years. Path coefficient analysis indicated that shoot mass at R8.5 had a direct effect on grain yield in both years, that root mass at R7 had a direct effect on grain yield in 1998, and that in 1999 the amounts of N and P in shoots at R8.5 had indirect effects on grain yield via shoot mass at R8.5. A combined analysis of both experiments revealed that biomass accumulation, N and P in shoots at R6 and R7 as well as root mass at R6 were similar in both years. In 1998 however bean accumulated more root mass at R7 and more biomass and N and P in shoots at R8 and R8.5, resulting in a 57 % higher grain yield in 1998. This indicates that grain yield of different common bean cultivars is not intrinsically associated with vegetative vigor at flowering and that mechanisms during pod filling can strongly influence the final crop yield. The establishment of a profuse root system during pod setting, associated with the continuous N and P acquisition during early pod filling, seems to be relevant for higher grain yields of common bean.

Critical evaluation of outcome scales assessment of lateral ankle ligament reconstruction.

BACKGROUND: Outcome following foot and ankle surgery can be assessed by disease- and region-specific scores. Many scoring systems exist, making comparison among studies difficult. The present study focused on outcome measures for a common foot and ankle abnormality and compared the results obtained by 2 disease-specific and 2 body region-specific scores. METHODS: We reviewed 41 patients who underwent lateral ankle ligament reconstruction. Four outcome scales were administered simultaneously: the Cumberland Ankle Instability Tool (CAIT) and the Chronic Ankle Instability Scale (CAIS), which are disease specific, and the American Orthopedic Foot & Ankle Society (AOFAS) hindfoot scale and the Foot and Ankle Ability Measure (FAAM), which are both body region-specific. The degree of correlation between scores was assessed by Pearson's correlation coefficient. Nonparametric tests, the Kruskal-Wallis and the Mann-Whitney test for pairwise comparison of the scores, were performed. RESULTS: A significant difference (P < .005) was observed between the CAIS and the AOFAS score (P = .0002), between the CAIS and the FAAM 1 (P = .0001), and between the CAIT and the AOFAS score (P = .0003). CONCLUSIONS: This study compared the performances of 4 disease- and body region-specific scoring systems. We demonstrated a correlation between the 4 administered scoring systems and notable differences between the results given by each of them. Disease-specific scores appeared more accurate than body region-specific scores. A strong correlation between the AOFAS score and the other scales was observed. The FAAM seemed a good compromise because it offered the possibility to evaluate the patient according to his or her own functional demand. CLINICAL RELEVANCE: The present study contributes to the development of more critical and accurate outcome assesment methods in foot and ankle surgery.

Análise de trilha no estudo dos efeitos de atributos físicos e matéria orgânica sobre a compressibilidade e resistência à penetração de um Argissolo cultivado com cana-de-açúcar

Uma das formas de avaliar a condição estrutural do solo é por meio da pressão de pré-compactação e da resistência do solo à penetração. Considerando que o conhecimento da associação entre duas variáveis é de grande importância na ciência do solo, o objetivo deste trabalho foi estudar o desdobramento dos efeitos diretos e indiretos, por meio da análise de correlação de trilha ou "path coefficient analysis", de atributos físicos e teor de matéria orgânica sobre a pressão de pré-compactação média (σpm) e resistência à penetração média em laboratório (RPLm) de um Argissolo Amarelo Distrocoeso cultivado com cana-de-açúcar nos Tabuleiros Costeiros do Estado de Alagoas. O trabalho foi desenvolvido na Fazenda Progresso, de propriedade da Usina Coruripe (AL), onde foram selecionadas cinco áreas cultivadas com cana-de-açúcar com diferentes tempos de uso e dois tipos de cultivo: 4 anos com cultivo de inverno (4 anos I); 14 anos com cultivo de inverno (14 anos I); 14 anos com cultivo de verão (14 anos V); 30 anos com cultivo de inverno (30 anos I); 30 anos com cultivo de verão (30 anos V); e uma área de mata (mata), tipo floresta subperenifólia, que foi considerada como testemunha-padrão. As avaliações foram realizadas nas profundidades de 0 a 0,20, 0,20 a 0,40 e 0,40 a 0,60 m, correspondentes aos horizontes A (mata) ou Ap (áreas cultivadas), AB e Bt, respectivamente, para as quais foram realizadas análises de correlação de "trilha" ("path coefficient analysis"), para determinação do grau de associação entre os atributos físicos e o teor de matéria orgânica do solo com a σpm e a RPLm (variáveis dependentes). Nas três profundidades, o grau de compactação e a densidade do solo foram os atributos que mostraram maior relação causa e efeito com a pressão de pré-compactação e resistência à penetração do Argissolo estudado. Nos horizontes A (mata) ou Ap (áreas cultivadas), o teor de matéria orgânica e a estabilidade de agregados em água apresentaram efeito direto sobre as variações da pressão de pré-compactação e resistência à penetração; a mata apresentou maior teor de matéria orgânica, maior porcentagem de agregados estáveis e maior compressibilidade, em relação às áreas cultivadas.

Influence of liming on residual soil respiration and chemical properties in a tropical no-tillage system

Because of the climate changes occurring across the planet, especially global warming, the different forms of agricultural soil use have attracted researchers´ attention. Changes in soil management may influence soil respiration and, consequently, C sequestration. The objectives of this study were to evaluate the long-term influence of liming on soil respiration and correlate it with soil chemical properties after two years of liming in a no-tillage system. A randomized complete block design was used with six replications. The experimental treatments consisted of four lime rates and a control treatment without lime. Two years after liming, soil CO2 emission was measured and the soil sampled (layers 0-5, 5-10, 10-20, and 20-30 cm). The P, Ca2+ e Mg2+ soil contents and pH and base saturation were determined. CO2 emission from soil limed at the recommended rate was 24.1 % higher, and at twice the recommended rate, 47.4 % higher than from unlimed soil. Liming improved the chemical properties, and the linear increase in soil respiration rate correlated positively with the P, Ca2+ and Mg2+ soil contents, pH and base saturation, and negatively with H + Al and Al3+ contents. The correlation coefficient between soil respiration rate and chemical properties was highest in the 10-20 cm layer.

Prescription footwear for severe injuries of foot and ankle: effect on regularity and symmetry of the gait assessed by trunk accelerometry.

After foot and/or ankle fracture, the restoration of optimal gait symmetry is one of the criteria of recovery. Orthotic insoles and orthopaedic shoes improve gait symmetry and regularity by controlling joint motion and improving alignment. The aim of the present study was to assess the effect of prescription footwear on gait quality by using accelerometers attached to the lower back. Sixteen adult patients with persistent disability after ankle and/or foot fractures performed two 30-s walking trials with and without prescription footwear (insoles and stabilizing shoes). Sixteen control subjects were also tested for comparison. The autocorrelation function was computed from the acceleration signal and the first two dominant periods were assessed (d1 and d2). Two parameters were used: (1) Stride Regularity (SR) which expresses the similarity between strides over time (d2), and (2) Stride Symmetry (SS) a ratio (d1/d2) which expresses the left/right similarity of gait independently of repeatability in the successive movements of each limb. In control subjects, SR and SS were 0.86+/-0.05 (correlation coefficient) and 81+/-10%, respectively. In the patient group, the effect of footwear was significant (SR: 0.88+/-0.06 vs. 0.90+/-0.05, SS: 38+/-23% vs. 46+/-27%). Pain was also significantly reduced (-34%). By using a rapid and low-cost method, we objectively quantified gait quality improvement after footwear intervention, concomitant to pain reduction. Substantial inter-patient variability in the footwear outcome was observed. In conclusion, we believe that trunk accelerometry can be a useful tool in the field of gait rehabilitation.

Can one or two high doses of oral vitamin D3 correct insufficiency in a non-supplemented rheumatologic population?

We evaluated the effectiveness of supplementation with high dose of oral vitamin D3 to correct vitamin D insufficiency. We have shown that one or two oral bolus of 300,000 IU of vitamin D3 can correct vitamin D insufficiency in 50% of patients and that the patients who benefited more from supplementation were those with the lowest baseline levels. INTRODUCTION: Adherence with daily oral supplements of vitamin D3 is suboptimal. We evaluated the effectiveness of a single high dose of oral vitamin D3 (300,000 IU) to correct vitamin D insufficiency in a rheumatologic population. METHODS: Over 1 month, 292 patients had levels of 25-OH vitamin D determined. Results were classified as: deficiency <10 ng/ml, insufficiency ≥10 to 30 ng/ml, and normal ≥30 ng/ml. We added a category using the IOM recommended cut-off of 20 ng/ml. Patients with deficient or normal levels were excluded, as well as patients already supplemented with vitamin D3. Selected patients (141) with vitamin D insufficiency (18.5 ng/ml (10.2-29.1) received a prescription for 300,000 IU of oral vitamin D3 and were asked to return after 3 (M3) and 6 months (M6). Patients still insufficient at M3 received a second prescription for 300,000 IU of oral vitamin D3. Relation between changes in 25-OH vitamin D between M3 and M0 and baseline values were assessed. RESULTS: Patients (124) had a blood test at M3. Two (2%) had deficiency (8.1 ng/ml (7.5-8.7)) and 50 (40%) normal results (36.7 ng/ml (30.5-5.5)). Seventy-two (58%) were insufficient (23.6 ng/ml (13.8-29.8)) and received a second prescription for 300,000 IU of oral vitamin D3. Of the 50/124 patients who had normal results at M3 and did not receive a second prescription, 36 (72%) had a test at M6. Seventeen (47%) had normal results (34.8 ng/ml (30.3-42.8)) and 19 (53%) were insufficient (25.6 ng/ml (15.2-29.9)). Of the 72/124 patients who receive a second prescription, 54 (75%) had a test at M6. Twenty-eight (52%) had insufficiency (23.2 ng/ml (12.8-28.7)) and 26 (48%) had normal results (33.8 ng/ml (30.0-43.7)). At M3, 84% patients achieved a 25-OH vitamin D level >20 ng/ml. The lowest the baseline value, the highest the change after 3 months (negative relation with a correlation coefficient r = -0.3, p = 0.0007). CONCLUSIONS: We have shown that one or two oral bolus of 300,000 IU of vitamin D3 can correct vitamin D insufficiency in 50% of patients.

Pedotransfer functions to estimate water retention parameters of soils in northeastern Brazil

Pedotransfer functions (PTF) were developed to estimate the parameters (α, n, θr and θs) of the van Genuchten model (1980) to describe soil water retention curves. The data came from various sources, mainly from studies conducted by universities in Northeast Brazil, by the Brazilian Agricultural Research Corporation (Embrapa) and by a corporation for the development of the São Francisco and Parnaíba river basins (Codevasf), totaling 786 retention curves, which were divided into two data sets: 85 % for the development of PTFs, and 15 % for testing and validation, considered independent data. Aside from the development of general PTFs for all soils together, specific PTFs were developed for the soil classes Ultisols, Oxisols, Entisols, and Alfisols by multiple regression techniques, using a stepwise procedure (forward and backward) to select the best predictors. Two types of PTFs were developed: the first included all predictors (soil density, proportions of sand, silt, clay, and organic matter), and the second only the proportions of sand, silt and clay. The evaluation of adequacy of the PTFs was based on the correlation coefficient (R) and Willmott index (d). To evaluate the PTF for the moisture content at specific pressure heads, we used the root mean square error (RMSE). The PTF-predicted retention curve is relatively poor, except for the residual water content. The inclusion of organic matter as a PTF predictor improved the prediction of parameter a of van Genuchten. The performance of soil-class-specific PTFs was not better than of the general PTF. Except for the water content of saturated soil estimated by particle size distribution, the tested models for water content prediction at specific pressure heads proved satisfactory. Predictions of water content at pressure heads more negative than -0.6 m, using a PTF considering particle size distribution, are only slightly lower than those obtained by PTFs including bulk density and organic matter content.

Test-retest reliability of thermal quantitative sensory testing on two sites within the L5 dermatome of the lumbar spine and lower extremity.

INTRODUCTION: Quantitative sensory testing (QST) is widely used in human research to investigate the integrity of the sensory function in patients with pain of neuropathic origin, or other causes such as low back pain. Reliability of QST has been evaluated on both sides of the face, hands and feet as well as on the trunk (Th3-L3). In order to apply these tests on other body-parts such as the lower lumbar spine, it is important first to establish reliability on healthy individuals. The aim of this study was to investigate intra-rater reliability of thermal QST in healthy adults, on two sites within the L5 dermatome of the lumbar spine and lower extremity. METHODS: Test-retest reliability of thermal QST was determined at the L5-level of the lumbar spine and in the same dermatome on the lower extremity in 30 healthy persons under 40 years of age. Results were analyzed using descriptive statistics and intraclass correlation coefficient (ICC). Values were compared to normative data, using Z-transformation. RESULTS: Mean intraindividual differences were small for cold and warm detection thresholds but larger for pain thresholds. ICC values showed excellent reliability for warm detection and heat pain threshold, good-to-excellent reliability for cold pain threshold and fair-to-excellent reliability for cold detection threshold. ICC had large ranges of confidence interval (95%). CONCLUSION: In healthy adults, thermal QST on the lumbar spine and lower extremity demonstrated fair-to-excellent test-retest reliability.

Reproducibility of twitch mouth pressure, sniff nasal inspiratory pressure, and maximal inspiratory pressure.

Twitch mouth pressure (Pmo,tw) during magnetic phrenic nerve stimulation and sniff nasal inspiratory pressure (SNIP) were recently proposed as alternative noninvasive methods for assessing inspiratory muscle strength. This study aimed to compare their reproducibility with maximal inspiratory pressure (MIP) in normal subjects. Ten healthy subjects were studied at functional residual capacity in semirecumbent position. Cervical magnetic phrenic nerve stimulation was performed during gentle expiration against an occlusion incorporating a small leak. Constancy of stimulation was controlled by recording diaphragmatic electromyogram. Within and between-session reproducibility of pressure were studied for Pmo,tw, SNIP, and MIP. The subjects were studied during a session of 10 manoeuvres repeated after 1 day and 1 month. The mean values were 16 cmH2O for Pmo,tw, 118 cmH2O for SNIP, and 115 cmH2O for MIP. For the three tests, the within subject variation was small in relation to between-subject variation, with the intraclass correlation coefficient ranging 0.79-0.90 for Pmo,tw, 0.85-0.92 for SNIP, and 0.88-0.92 for MIP. At 1 day interval, the coefficient of repeatability (2 SD of differences) was 3.6 cmH2O for Pmo,tw, 32 cmH2O for SNIP and 28 cmH2O for MIP. At 1 month interval, the coefficient of repeatability was 5.8 cmH2O for Pmo,tw, 23 cmH2O for SNIP and 21 cmH2O for MIP. We conclude that the within session reproducibility of the new tests twitch mouth pressure and sniff nasal inspiratory pressure is sufficient to be clinically useful. For sniff nasal inspiratory pressure, the between session reproducibility established after 1 day was maintained after 1 month. For twitch mouth pressure, the between session reproducibility declined slightly after 1 month. These characteristics should be considered when using these methods to follow an individual patient over time.

Secular changes of spinal canal dimensions in Western Switzerland: a narrowing epidemic?

STUDY DESIGN: Computed tomography-based anatomical study. OBJECTIVE: To study the secular changes in lumbar spinal canal dimensions. SUMMARY OF BACKGROUND DATA: Development of symptomatic lumbar spinal stenosis, among other factors, is related to the dimensions of the bony canal. The canal reaches its adult size early on in life. Several factors, including protein intake, may influence its final dimensions. As with increases in human stature from improvements of socioeconomic conditions, we hypothesized that adult bony canal size has also grown larger in recent generations. METHODS: This study analyzes computed tomographic reconstructions from 184 subjects performed for either trauma (n = 81) or abdominal pathologies (n = 103) and born either between 1940 and 1949 (n = 88) or 1970 and 1979 (n = 96). The cross-sectional area of the bony canal was digitally measured at the level of the pedicle (i.e., at a level not influenced by degenerative changes) for each lumbar vertebra. Intra- and interobserver reliability was assessed. RESULTS: Intra- and interobserver measurement reliability were excellent (interclass correlation coefficient = 0.87) and good (interclass correlation coefficient = 0.61), respectively. Contrary to our hypothesis, the 1940-1949 generation patient group exhibited larger lumbar canals at all levels as compared with the 1970-1979 group. Statistically this difference was highly significant (P < 0.001) and particularly pronounced in the trauma subgroup. CONCLUSION: Given that human stature evolution has stabilized and adult height is established during the first 2 years of long bone growth, it is possible that antenatal factors are responsible for this surprising finding. Maternal smoking and age may be possible explanations. This finding may have significant implications. An increasing number of patients may emerge with lumbar spinal stenosis as degenerative changes develop, putting a strain on health resources. Further studies in different population groups and countries will be important to further confirm this trend. LEVEL OF EVIDENCE: 3.

Protein-ligand binding free energy estimation using molecular mechanics and continuum electrostatics. Application to HIV-1 protease inhibitors.

A method is proposed for the estimation of absolute binding free energy of interaction between proteins and ligands. Conformational sampling of the protein-ligand complex is performed by molecular dynamics (MD) in vacuo and the solvent effect is calculated a posteriori by solving the Poisson or the Poisson-Boltzmann equation for selected frames of the trajectory. The binding free energy is written as a linear combination of the buried surface upon complexation, SASbur, the electrostatic interaction energy between the ligand and the protein, Eelec, and the difference of the solvation free energies of the complex and the isolated ligand and protein, deltaGsolv. The method uses the buried surface upon complexation to account for the non-polar contribution to the binding free energy because it is less sensitive to the details of the structure than the van der Waals interaction energy. The parameters of the method are developed for a training set of 16 HIV-1 protease-inhibitor complexes of known 3D structure. A correlation coefficient of 0.91 was obtained with an unsigned mean error of 0.8 kcal/mol. When applied to a set of 25 HIV-1 protease-inhibitor complexes of unknown 3D structures, the method provides a satisfactory correlation between the calculated binding free energy and the experimental pIC5o without reparametrization.