Cross-sectional associations between questionnaire-measured physical activity and tissue doppler indices of left ventricular diastolic function.

BACKGROUND: The prevalence of left ventricular (LV) diastolic dysfunction has been increasing over the past decade, and to date, effective pharmacotherapies that enhance LV diastolic function have not yet been identified. Though some data has demonstrated the beneficial effects of exercise training on LV diastolic function, little is known about the adaptations of diastolic function to daily physical activity (PA). Accordingly, our study aimed to investigate the impact of daily PA on tissue Doppler indices of LV diastolic function. METHODS: A total of 432 participants were enrolled for clinically indicated echocardiography from July 2019 to July 2020 at Peking University People's Hospital. Participants aged ≥ 18 years were included if they had stable PA in the past six months and normal LV systolic function. A questionnaire was used to collect demographic characteristics, medical history, and daily PA. According to PA Guidelines for Americans, we identified these participants into low-intensity PA (LPA) group and moderate-high-intensity PA (MHPA) group. Propensity score matching (PSM) was performed to match potential confounding factors between the two groups. The clinical characteristics and echocardiographic parameters between LPA group and MHPA group were compared using student's t-test, Mann-Whitney U test, and chi-square test as appropriate. RESULTS: After matching potential confounding factors using PSM with a 1:3 matching ratio, our final analysis included 86 cases in the MHPA group and 214 cases in the LPA group. All demographic characteristics and comorbidities were statistically similar between the two groups. Compared to the LPA group, the MHPA group showed higher septal e' (7.9 ± 2.9 cm/s versus 7.2 ± 2.6 cm/s, P = 0.047). Other echocardiographic parameters associated with LV diastolic function concerning lateral e' and average E/e', also trended towards improved LV diastolic function in the MHPA group, but failed to reach statistical significance. CONCLUSIONS: Our study demonstrated that moderate-high-intensity daily PA was associated with improved septal e', suggesting that moderate-high-intensity PA could potentially ameliorate LV diastolic dysfunction.

The Effects of Vitamin D on Markers of Glucose and Obesity in Postmenopausal Women: A Meta-analysis of Randomized Controlled Trials.

PURPOSE: The effect of vitamin D effect on glucose markers and obesity in postmenopausal women remains controversial. The current literature contains little information on vitamin D dosage and duration for optimal efficacy in postmenopausal women. This meta-analysis was undertaken to assess the impact of vitamin D on glucose markers and obesity in postmenopausal women. METHODS: A number of databases were used dated up to January 5, 2023, with no language restrictions (PubMed/MEDLINE, Web of Science, EMBASE, and Scopus). Treatment response from baseline was estimated from the mean within-group analysis, and SDs were used to calculate the treatment response. FINDINGS: Nine eligible articles with 12 comparisons qualified for the final quantitative analysis. An overall decrease was noted in fasting blood glucose (weighted mean difference [WMD], -3.56 mg/dL; 95% CI, -5.49 to -1.64; P < 0.001), homeostatic model assessment for insulin resistance (WMD, -1.168 mm; 95% CI, -2.001 to -0.33; P = 0.006), insulin (WMD, -2.26 units; 95% CI, -4.35 to -0.18; P = 0.033), and glycosylated hemoglobin (WMD, -0.41%; 95% CI, -0.54 to -0.29; P < 0.001) after vitamin D administration in postmenopausal women. In subgroup analyses, a notable decrease in fasting blood glucose was detected when the intervention course was ˃6 months and dosage ≤1000 IU/d (WMD, -3.48 mg/dL). The present study showed that vitamin D was not associated with body mass index, body weight, or waist circumference in postmenopausal women. IMPLICATIONS: Vitamin D is beneficial for glucose markers but not obesity in postmenopausal women. An individualized dosage regimen of vitamin D should be followed depending on the clinical outcome target of postmenopausal women.

Influence of statin use on prognosis of patients with renal cell cancer: a meta-analysis.

Background: Statin may confer anticancer efficacy, while the studies evaluating the influence of statin on survival of patients with renal cell cancer (RCC) yielded inconsistent results. A systematic review and meta-analysis was performed to investigate the association between statin use and survival of patients with RCC. Materials and Methods: Cohort studies were identified by search of PubMed, Embase, and Web of Science databases according to the objective of the meta-analysis. A random-effect model incorporating the possible between-study heterogeneity was used for meta-analysis. Subgroup analyses according to study characteristics were also performed. Results: Seventeen cohort studies involving 42528 patients with RCC were available for the meta-analysis. Results showed that statin use was associated with a better overall survival (OS, hazard ratio [HR]: 0.73, 95% confidence interval [CI]: 0.65 to 0.84, p < 0.001; I2 = 40%), progression progression-free survival (PFS, HR: 0.82, 95% CI: 0.68 to 0.98, p = 0.03; I2 = 52%), and cancer-specific survival (CSS, HR: 0.76, 95% CI: 0.59 to 0.99, p = 0.04; I2 = 38%). Besides, for the outcome of OS and PFS, subgroup analyses showed similar results in patients with surgical and non-surgical anticancer treatments, and in patients with stage I-III and stage IV RCC (p values for subgroup difference all > 0.05). Conclusions: Statin use may be associated with improved survival outcomes in patients with RCC. Although prospective clinical studies should be considered to validate these results, these findings suggest that statins may be potential adjuvant therapy for patients with RCC.

Hirudo Lyophilized Powder Ameliorates Renal Injury in Diabetic Rats by Suppressing Oxidative Stress and Inflammation.

As diabetic nephropathy (DN) is one of the most common and destructive microvascular complications of diabetes mellitus, the goal of this study, therefore, was to investigate the renal protective effect and latent mechanisms of Hirudo lyophilized powder on diabetic rats. In this study, all rats were randomly assigned into the control group and diabetic group. The rats of diabetic group were injected with low-dose STZ (35 mg/kg) intraperitoneal plus high-fat diet to induce diabetes. Then, the successful diabetic model rats were weighed and randomly assigned into four groups: (1) diabetic model group (DM group); (2) Hirudo lyophilized powder 0.3 g/kg treatment group (SL group); (3) Hirudo lyophilized powder 0.6 g/kg treatment group (SM group); (4) Hirudo lyophilized powder 1.2 g/kg treatment group (SH group). Their fasting blood glucoses (FBG) were measured every 4 weeks. After treatment with Hirudo lyophilized powder at a corresponding dose once a day for 16 weeks, their metabolic and biochemical as well as oxidative stress parameters were tested, and the kidney weight (KW)/body weight (BW) was calculated. The renal tissues were used for histological, mRNA, and protein expression analysis. The results showed that Hirudo lyophilized powder could protect against the structural damages and functional changes of diabetic renal tissue by inhibiting oxidative stress, inflammation, and fibrosis. Furthermore, it was found in the further research that inhibiting the NOX4 expression and JAK2/STAT1/STAT3 pathway activation might be the underlying mechanisms. Collectively, Hirudo lyophilized powder might be a promising therapeutic agent for the treatment of DN.

Clinical imaging characteristics of inpatients with coronavirus disease-2019 in Heilongjiang Province, China: a retrospective study.

OBJECTIVE: To investigate the clinical, laboratory, and radiological characteristics of patients with coronavirus disease-2019 (COVID-19) in Heilongjiang Province. RESULTS: Patients in the ICU group were older and their incidence of cardiovascular disease was higher than those in the non-ICU group. Lymphocyte levels were lower and neutrophil and D-dimer levels were higher in the ICU than that in the non-ICU group. Compared to the non-ICU group, the incidence of pulmonary consolidation and ground-glass opacity with consolidation was significantly higher in the ICU group, all lung lobes were more likely to be involved, with higher number of lung lobes and areas surrounding the bronchi. Of the 59 patients with COVID-19 in this group, 15 received mechanical ventilation. All intubated patients involved lung lobes, and a large number of lesions were observed in the area around the bronchial vessels. CONCLUSION: Significant differences were observed in clinical symptoms, laboratory tests, and computed tomography features between the ICU and non-ICU groups. METHODS: A total of 59 patients with COVID-19, comprising 44 patients in the intensive care unit (ICU) and 15 in the non-ICU, were retrospectively analyzed. Characteristics of the two groups of patients were compared.

Qingshen Buyang Formula Attenuates Renal Fibrosis in 5/6 Nephrectomized Rats via Inhibiting EMT and Wnt/ß-Catenin Pathway.

As renal fibrosis significantly contributes to various kinds of chronic kidney diseases, this study aimed to investigate the renal protective effects of Qingshen Buyang Formula against renal fibrosis on 5/6 nephrectomized rats, and its underlying mechanisms were explored. A total of 24 male Sprague-Dawley rats were randomly divided into sham operation group (Sham group), 5/6 nephrectomy group (5/6Nx group), and Qingshen Buyang Formula treatment group (QBF group). The intervention was intragastric administration for 12 weeks. In the end, the blood samples were collected to test renal functional parameters, urine proteins were measured, and the left kidneys were removed for histological studies, as well as mRNA and protein expression analysis. The results showed that Qingshen Buyang Formula significantly decreased BUN, Scr, and proteinuria in 5/6Nx rats. Meanwhile, it ameliorated the kidney injury and fibrosis, exemplified by the depressed expression of collagen I and fibronectin (FN), which are the main components of ECM. Furthermore, the process of EMT inhibited the Wnt/ß-catenin signaling pathway related genes, such as Wnt4, TCF4, ß-catenin, and p-GSK3ß. Collectively, the Qingshen Buyang Formula can improve renal function and attenuate renal fibrosis, and its underlying mechanisms may be related with inhibiting EMT and Wnt/ß-catenin signaling pathway.

Potent and Selective Inhibition of Plasma Membrane Monoamine Transporter by HIV Protease Inhibitors.

Plasma membrane monoamine transporter (PMAT) is a major uptake-2 monoamine transporter that shares extensive substrate and inhibitor overlap with organic cation transporters 1-3 (OCT1-3). Currently, there are no PMAT-specific inhibitors available that can be used in in vitro and in vivo studies to differentiate between PMAT and OCT activities. In this study, we showed that IDT307 (4-(4-(dimethylamino)phenyl)-1-methylpyridinium iodide), a fluorescent analog of 1-methyl-4-phenylpyridinium (MPP+), is a transportable substrate for PMAT and that IDT307-based fluorescence assay can be used to rapidly identify and characterize PMAT inhibitors. Using the fluorescent substrate-based assays, we analyzed the interactions of eight human immunodeficiency virus (HIV) protease inhibitors (PIs) with human PMAT and OCT1-3 in human embryonic kidney 293 (HEK293) cells stably transfected with individual transporters. Our data revealed that PMAT and OCTs exhibit distinct sensitivity and inhibition patterns toward HIV PIs. PMAT is most sensitive to PI inhibition whereas OCT2 and OCT3 are resistant. OCT1 showed an intermediate sensitivity and a distinct inhibition profile from PMAT. Importantly, lopinavir is a potent PMAT inhibitor and exhibited >120 fold selectivity toward PMAT (IC50 = 1.4 ± 0.2 µM) over OCT1 (IC50 = 174 ± 40 µM). Lopinavir has no inhibitory effect on OCT2 or OCT3 at maximal tested concentrations. Lopinavir also exhibited no or much weaker interactions with uptake-1 monoamine transporters. Together, our results reveal that PMAT and OCTs have distinct specificity exemplified by their differential interaction with HIV PIs. Further, we demonstrate that lopinavir can be used as a selective PMAT inhibitor to differentiate PMAT-mediated monoamine and organic cation transport from those mediated by OCT1-3.

PubChem applications in drug discovery: a bibliometric analysis.

A bibliometric analysis of PubChem applications is presented by reviewing 1132 research articles. The massive volume of chemical structure and bioactivity data in PubChem and its online services have been used globally in various fields including chemical biology, medicinal chemistry and informatics research. PubChem supports drug discovery in many aspects such as lead identification and optimization, compound-target profiling, polypharmacology studies and unknown chemical identity elucidation. PubChem has also become a valuable resource for developing secondary databases, informatics tools and web services. The growing PubChem resource with its public availability offers support and great opportunities for the interrogation of pharmacological mechanisms and the genetic basis of diseases, which are vital for drug innovation and repurposing.

Activation of the constitutive androstane receptor inhibits gluconeogenesis without affecting lipogenesis or fatty acid synthesis in human hepatocytes.

OBJECTIVE: Accumulating evidence suggests that activation of mouse constitutive androstane receptor (mCAR) alleviates type 2 diabetes and obesity by inhibiting hepatic gluconeogenesis, lipogenesis, and fatty acid synthesis. However, the role of human (h) CAR in energy metabolism is largely unknown. The present study aims to investigate the effects of selective hCAR activators on hepatic energy metabolism in human primary hepatocytes (HPH). METHODS: Ligand-based structure-activity models were used for virtual screening of the Specs database (www.specs.net) followed by biological validation in cell-based luciferase assays. The effects of two novel hCAR activators (UM104 and UM145) on hepatic energy metabolism were evaluated in HPH. RESULTS: Real-time PCR and Western blotting analyses reveal that activation of hCAR by UM104 and UM145 significantly repressed the expression of glucose-6-phosphatase and phosphoenolpyruvate carboxykinase, two pivotal gluconeogenic enzymes, while exerting negligible effects on the expression of genes associated with lipogenesis and fatty acid synthesis. Functional experiments show that UM104 and UM145 markedly inhibit hepatic synthesis of glucose but not triglycerides in HPH. In contrast, activation of mCAR by 1,4-bis[2-(3,5-dichloropyridyloxy)]benzene, a selective mCAR activator, repressed the expression of genes associated with gluconeogenesis, lipogenesis, and fatty acid synthesis in mouse primary hepatocytes, which were consistent with previous observations in mouse model in vivo. CONCLUSION: Our findings uncover an important species difference between hCAR and mCAR in hepatic energy metabolism, where hCAR selectively inhibits gluconeogenesis without suppressing fatty acid synthesis. IMPLICATIONS: Such species selectivity should be considered when exploring CAR as a potential therapeutic target for metabolic disorders.

Pathway analysis for drug repositioning based on public database mining.

Sixteen FDA-approved drugs were investigated to elucidate their mechanisms of action (MOAs) and clinical functions by pathway analysis based on retrieved drug targets interacting with or affected by the investigated drugs. Protein and gene targets and associated pathways were obtained by data-mining of public databases including the MMDB, PubChem BioAssay, GEO DataSets, and the BioSystems databases. Entrez E-Utilities were applied, and in-house Ruby scripts were developed for data retrieval and pathway analysis to identify and evaluate relevant pathways common to the retrieved drug targets. Pathways pertinent to clinical uses or MOAs were obtained for most drugs. Interestingly, some drugs identified pathways responsible for other diseases than their current therapeutic uses, and these pathways were verified retrospectively by in vitro tests, in vivo tests, or clinical trials. The pathway enrichment analysis based on drug target information from public databases could provide a novel approach for elucidating drug MOAs and repositioning, therefore benefiting the discovery of new therapeutic treatments for diseases.

Pharmacophore and 3D-QSAR characterization of 6-arylquinazolin-4-amines as Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) inhibitors.

Cdc2-like kinase 4 (Clk4) and dual specificity tyrosine-phosphorylation-regulated kinase 1A (Dyrk1A) are protein kinases that are promising targets for treatment of diseases caused by abnormal gene splicing. 6-Arylquinazolin-4-amines have been recently identified as potent Clk4 and Dyrk1A inhibitors. In order to understand the structure-activity correlation of these analogs, we have applied ligand-based pharmacophore and 3D-QSAR modeling combined with structure-based homology modeling and docking. The high R(2) and Q(2) (0.88 and 0.79 for Clk4, 0.85 and 0.82 for Dyrk1A, respectively) based on validation with training and test set compounds suggested that the generated 3D-QSAR models are reliable in predicting novel ligand activities against Clk4 and Dyrk1A. The binding mode identified through docking ligands to the ATP binding domain of Clk4 was consistent with the structural properties and energy field contour maps characterized by pharmacophore and 3D-QSAR models and gave valuable insights into the structure-activity profile of 6-arylquinazolin-4-amine analogs. The obtained 3D-QSAR and pharmacophore models in combination with the binding mode between inhibitor and residues of Clk4 will be helpful for future lead compound identification and optimization to design potent and selective Clk4 and Dyrk1A inhibitors.

Identification of novel breast cancer resistance protein (BCRP) inhibitors by virtual screening.

Breast cancer resistance protein (BCRP; ABCG2) is an efflux transporter that plays an important role in multidrug resistance to antineoplastic drugs. The identification of drugs as BCRP inhibitors could aid in designing better therapeutic strategies for cancer treatment and will be critical for identifying potential drug-drug interactions. In the present study, we applied ligand-based virtual screening combined with experimental testing for the identification of novel drugs that can possibly interact with BCRP. Bayesian and pharmacophore models generated with known BCRP inhibitors were validated with an external test set. The resulting models were applied to predict new potential drug candidates from a database with more than 2000 FDA-approved drugs. Thirty-three drugs were tested in vitro for their inhibitory effects on BCRP-mediated transport of [(3)H]-mitoxantrone in MCF-7/AdrVp cells. Nineteen drugs were identified with significant inhibitory effect on BCRP transport function. The combined strategy of computational and experimental approaches in this paper has suggested potential drug candidates and thus represents an effective tool for rational identification of modulators of other proteins.

Identification of novel activators of constitutive androstane receptor from FDA-approved drugs by integrated computational and biological approaches.

PURPOSE: The constitutive androstane receptor (CAR, NR1I3) is a xenobiotic sensor governing the transcription of numerous hepatic genes associated with drug metabolism and clearance. Recent evidence suggests that CAR also modulates energy homeostasis and cancer development. Thus, identification of novel human (h) CAR activators is of both clinical importance and scientific interest. METHODS: Docking and ligand-based structure-activity models were used for virtual screening of a database containing over 2000 FDA-approved drugs. Identified lead compounds were evaluated in cell-based reporter assays to determine hCAR activation. Potential activators were further tested in human primary hepatocytes (HPHs) for the expression of the prototypical hCAR target gene CYP2B6. RESULTS: Nineteen lead compounds with optimal modeling parameters were selected for biological evaluation. Seven of the 19 leads exhibited moderate to potent activation of hCAR. Five out of the seven compounds translocated hCAR from the cytoplasm to the nucleus of HPHs in a concentration-dependent manner. These compounds also induce the expression of CYP2B6 in HPHs with rank-order of efficacies closely resembling that of hCAR activation. CONCLUSION: These results indicate that our strategically integrated approaches are effective in the identification of novel hCAR modulators, which may function as valuable research tools or potential therapeutic molecules.

Molecular analysis and structure-activity relationship modeling of the substrate/inhibitor interaction site of plasma membrane monoamine transporter.

Plasma membrane monoamine transporter (PMAT) is a new polyspecific transporter that interacts with a wide range of structurally diverse organic cations. To map the physicochemical descriptors of cationic compounds that allow interaction with PMAT, we systematically analyzed the interactions between PMAT and three series of structural analogs of known organic cation substrates including phenylalkylamines, n-tetraalkylammonium (n-TAA) compounds, and ß-carbolines. Our results showed that phenylalkylamines with a distance between the aromatic ring and the positively charged amine nitrogen atom of ∼6.4 Å confer optimal interactions with PMAT, whereas studies with n-TAA compounds revealed an excellent correlation between IC(50) values and hydrophobicity. The five ß-carbolines that we tested, which possess a pyridinium-like structure and are structurally related to the neurotoxin 1-methyl-4-phenylpyridinium, inhibited PMAT with high affinity (IC(50) values of 39.1-65.5 µM). Cytotoxicity analysis further showed that cells expressing PMAT are 14- to 15-fold more sensitive to harmalan and norharmanium, suggesting that these two ß-carbolines are also transportable substrates of PMAT. We then used computer-aided modeling to generate qualitative and quantitative three-dimensional pharmacophore models on the basis of 23 previously reported and currently identified PMAT inhibitors and noninhibitors. These models are characterized by a hydrogen bond donor and two to three hydrophobic features with distances between the hydrogen bond donor and hydrophobic features ranging between 5.20 and 7.02 Å. The consistency between the mapping results and observed PMAT affinity of a set of test compounds indicates that the models performed well in inhibitor prediction and could be useful for future virtual screening of new PMAT inhibitors.

Active site gating and substrate specificity of butyrylcholinesterase and acetylcholinesterase: insights from molecular dynamics simulations.

Butyrylcholinesterase (BChE) and acetylcholinesterase (AChE) are highly homologous proteins with distinct substrate preferences. In this study we compared the active sites of monomers and tetramers of human BChE and human AChE after performing molecular dynamics (MD) simulations in water-solvated systems. By comparing the conformational dynamics of gating residues of AChE and BChE, we found that the gating mechanisms of the main door of AChE and BChE are responsible for their different substrate specificities. Our simulation of the tetramers of AChE and BChE indicates that both enzymes could have two dysfunctional active sites due to their restricted accessibility to substrates. The further study on catalytic mechanisms of multiple forms of AChE and BChE would benefit from our comparison of the active sites of the monomers and tetramers of both enzymes.

Identification and validation of novel human pregnane X receptor activators among prescribed drugs via ligand-based virtual screening.

Human pregnane X receptor (hPXR) plays a key role in regulating metabolism and clearance of endogenous and exogenous substances. Identification of novel hPXR activators among commercial drugs may aid in avoiding drug-drug interactions during coadministration. We applied ligand-based computational approaches for virtual screening of a commonly prescribed drug database (SCUT). Bayesian classification models were generated with a training set comprising 177 compounds using Fingerprints and 117 structural descriptors. A cell-based luciferase reporter assay was used for evaluation of chemical-mediated hPXR activation in HepG2 cells. All compounds were tested at 10 µM concentration with rifampicin and dimethyl sulfoxide as positive and negative controls, respectively. The Bayesian models showed specificity and overall prediction accuracy up to 0.92 and 0.69 for test set compounds. Screening the SCUT database with this model retrieved 105 hits and 17 compounds from the top 25 hits were chosen for in vitro testing. The reporter assay confirmed that nine drugs, i.e., fluticasone, nimodipine, nisoldipine, beclomethasone, finasteride, flunisolide, megestrol, secobarbital, and aminoglutethimide, were previously unidentified hPXR activators. Thus, the present study demonstrates that novel hPXR activators can be efficiently identified among U.S. Food and Drug Administration-approved and commonly prescribed drugs, which should lead to detection and prevention of potential drug-drug interactions.

Reaction pathway and free energy profile for prechemical reaction step of human butyrylcholinesterase-catalyzed hydrolysis of (-)-cocaine by combined targeted molecular dynamics and potential of mean force simulations.

Combined targeted molecular dynamics (TMD) and potential of mean force (PMF) simulations have been carried out to uncover the detailed pathway and determine the corresponding free energy profile for the structural transformation from the nonprereactive butyrylcholinesterase (BChE)-(-)-cocaine binding to the prereactive BChE-(-)-cocaine binding associated with the (-)-cocaine rotation in the binding pocket of BChE. It has been shown that the structural transformation involves two transition states (TS1(rot) and TS2(rot)). TS1(rot) is mainly associated with the deformation of the nonprereactive complex, whereas TS2(rot) is mainly associated with the formation of the prereactive complex. It has also been demonstrated that the A328W/Y332G mutation significantly reduces the steric hindrance for (-)-cocaine rotation in the binding pocket of BChE and, thus, decreases the free energy barrier for the structural transformation from the nonprereactive binding to the prereactive binding. The calculated relative free energy barriers are all consistent with available experimental kinetic data. The new mechanistic insights obtained and the novel computational protocol tested in this study should be valuable for future computational design of high-activity mutants of BChE. The general computational strategy and approach based on the combined TMD and PMF simulations may be also valuable in computational studies of detailed pathways and free energy profiles for other similar mechanistic problems involving ligand rotation or another type of structural transformation in the binding pocket of a protein.

Structural requirements of the ASBT by 3D-QSAR analysis using aminopyridine conjugates of chenodeoxycholic acid.

The human apical sodium-dependent bile acid transporter (ASBT) is a validated drug target and can be employed to increase oral bioavailability of various drug conjugates. The aim of the present study was to investigate the chemical space around the 24-position of bile acids that influences both inhibition and uptake by the transporter. A series of 27 aminopyridine and aminophenol conjugates of glutamyl-chenodeoxycholate were synthesized and their ASBT inhibition and transport kinetics (parametrized as K(i), K(t), and J(max)) measured using stably transfected ASBT-MDCK cells. All conjugates were potent ASBT inhibitors. Monoanionic conjugates exhibited higher inhibition potency than neutral conjugates. However, neutral conjugates and chloro-substituted monoanionic conjugates were not substrates, or at least not apparent substrates. Kinetic analysis of substrates indicated that similar values for K(i) and K(t) implicate substrate binding to ASBT as the rate-limiting step. Using 3D-QSAR, four inhibition models and one transport efficiency model were developed. Steric fields dominated in CoMFA models, whereas hydrophobic fields dominated CoMSIA models. The inhibition models showed that a hydrophobic or bulky substitute on the 2 or 6 position of a 3-aminopyridine ring enhanced activity, while a hydrophobic group on the 5 position was detrimental. Overall, steric and hydrophobic features around the 24 position of the sterol nucleus strongly influenced bile acid conjugate interaction with ASBT. The relative location of the pyridine nitrogen and substituent groups also modulated binding.

Free energy perturbation simulation on transition states and high-activity mutants of human butyrylcholinesterase for (-)-cocaine hydrolysis.

A unified computational approach based on free energy perturbation (FEP) simulations of transition states has been employed to calculate the mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state for (-)-cocaine hydrolysis catalyzed by the currently most promising series of mutants of human butyrylcholinesterase (BChE) that contain the A199S/A328W/Y332G mutations. The FEP simulations were followed by Michaelis-Menten kinetics analysis determining the individual k(cat) and K(M) values missing for the A199S/F227A/A328W/Y332G mutant in this series. The calculated mutation-caused shifts of the free energy change from the free enzyme to the rate-determining transition state are in good agreement with the experimental kinetic data, demonstrating that the unified computational approach based on the FEP simulations of the transition states may be valuable for future computational design of new BChE mutants with a further improved catalytic efficiency against (-)-cocaine.

Model of human butyrylcholinesterase tetramer by homology modeling and dynamics simulation.

A mutant of human butyrylcholinesterase (BChE) with high activity against cocaine would be highly promising as a drug for therapeutic treatment of cocaine abuse and overdose. It is desirable to design a recombinant BChE mutant with a long half-life in human circulation. Studies showed that BChE subunits can be assembled by a peptide containing the proline-rich attachment domain (PRAD) to form a stable tetramer. The models of BChE tetramer complexed with PRAD with various sequences have been constructed, in the present study, on the basis of homology modeling and molecular dynamics simulation of explicit water-solvated systems. The 3D models enable us to understand how the BChE subunits are arranged in the tetramer and how the tetramerization domain of BChE is associated with PRAD to form a stable tetramer of human BChE. It has been shown that the six conserved hydrophobic residues located on the C-terminal of BChE are responsible for the key electrostatic and hydrophobic interactions between the tetramerization domain of BChE and PRAD. The simulated tetramer structures suggest that mutation of three residues, i.e., Phe547, Met554, and Phe561, to other hydrophobic residues may be beneficial for increasing the binding between the tetramerization domain of BChE and PRAD. Thus, the detailed structural insights obtained from this study may be valuable for rational design of a recombinant BChE tetramer with a longer residence time in circulation.