Accelerating self-consistent field convergence with the augmented Roothaan-Hall energy function.

Based on Pulay's direct inversion iterative subspace (DIIS) approach, we present a method to accelerate self-consistent field (SCF) convergence. In this method, the quadratic augmented Roothaan-Hall (ARH) energy function, proposed recently by Høst and co-workers [J. Chem. Phys. 129, 124106 (2008)], is used as the object of minimization for obtaining the linear coefficients of Fock matrices within DIIS. This differs from the traditional DIIS of Pulay, which uses an object function derived from the commutator of the density and Fock matrices. Our results show that the present algorithm, abbreviated ADIIS, is more robust and efficient than the energy-DIIS (EDIIS) approach. In particular, several examples demonstrate that the combination of ADIIS and DIIS ("ADIIS+DIIS") is highly reliable and efficient in accelerating SCF convergence.

Quality control autophagy: a joint effort of ubiquitin, protein deacetylase and actin cytoskeleton.

Autophagy has been predominantly studied as a nonselective self-digestion process that recycles macromolecules and produces energy in response to starvation. However, autophagy independent of nutrient status has long been known to exist. Recent evidence suggests that this form of autophagy enforces intracellular quality control by selectively disposing of aberrant protein aggregates and damaged organelles--common denominators in various forms of neurodegenerative diseases. By definition, this form of autophagy, termed quality-control (QC) autophagy, must be different from nutrient-regulated autophagy in substrate selectivity, regulation and function. We have recently identified the ubiquitin-binding deacetylase, HDAC6, as a key component that establishes QC. HDAC6 is not required for autophagy activation per se; rather, it is recruited to ubiquitinated autophagic substrates where it stimulates autophagosome-lysosome fusion by promoting F-actin remodeling in a cortactin-dependent manner. Remarkably, HDAC6 and cortactin are dispensable for starvation-induced autophagy. These findings reveal that autophagosomes associated with QC are molecularly and biochemically distinct from those associated with starvation autophagy, thereby providing a new molecular framework to understand the emerging complexity of autophagy and therapeutic potential of this unique machinery.

Electrostatic Energetics of Bacillus subtilis Ribonuclease P Protein Determined by Nuclear Magnetic Resonance-Based Histidine pKa Measurements.

The pKa values of ionizable groups in proteins report the free energy of site-specific proton binding and provide a direct means of studying pH-dependent stability. We measured histidine pKa values (H3, H22, and H105) in the unfolded (U), intermediate (I), and sulfate-bound folded (F) states of RNase P protein, using an efficient and accurate nuclear magnetic resonance-monitored titration approach that utilizes internal reference compounds and a parametric fitting method. The three histidines in the sulfate-bound folded protein have pKa values depressed by 0.21 ± 0.01, 0.49 ± 0.01, and 1.00 ± 0.01 units, respectively, relative to that of the model compound N-acetyl-l-histidine methylamide. In the unliganded and unfolded protein, the pKa values are depressed relative to that of the model compound by 0.73 ± 0.02, 0.45 ± 0.02, and 0.68 ± 0.02 units, respectively. Above pH 5.5, H22 displays a separate resonance, which we have assigned to I, whose apparent pKa value is depressed by 1.03 ± 0.25 units, which is ∼0.5 units more than in either U or F. The depressed pKa values we observe are consistent with repulsive interactions between protonated histidine side chains and the net positive charge of the protein. However, the pKa differences between F and U are small for all three histidines, and they have little ionic strength dependence in F. Taken together, these observations suggest that unfavorable electrostatics alone do not account for the fact that RNase P protein is intrinsically unfolded in the absence of ligand. Multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.

Protein and Drug Design Algorithms Using Improved Biophysical Modeling

This thesis focuses on the development of algorithms that will allow protein design calculations to incorporate more realistic modeling assumptions. Protein design algorithms search large sequence spaces for protein sequences that are biologically and medically useful. Better modeling could improve the chance of success in designs and expand the range of problems to which these algorithms are applied. I have developed algorithms to improve modeling of backbone flexibility (DEEPer) and of more extensive continuous flexibility in general (EPIC and LUTE). I’ve also developed algorithms to perform multistate designs, which account for effects like specificity, with provable guarantees of accuracy (COMETS), and to accommodate a wider range of energy functions in design (EPIC and LUTE).

Calcium/Calmodulin-Dependent Protein Kinase Kinase 2 (CaMKK2) Regulates Dendritic Cells and Myeloid Derived Suppressor Cells Development in the Lymphoma Microenvironment

Calcium (Ca2+) is a known important second messenger. Calcium/Calmodulin (CaM) dependent protein kinase kinase 2 (CaMKK2) is a crucial kinase in the calcium signaling cascade. Activated by Ca2+/CaM, CaMKK2 can phosphorylate other CaM kinases and AMP-activated protein kinase (AMPK) to regulate cell differentiation, energy balance, metabolism and inflammation. Outside of the brain, CaMKK2 can only be detected in hematopoietic stem cells and progenitors, and in the subsets of mature myeloid cells. CaMKK2 has been noted to facilitate tumor cell proliferation in prostate cancer, breast cancer, and hepatic cancer. However, whethter CaMKK2 impacts the tumor microenvironment especially in hematopoietic malignancies remains unknown. Due to the relevance of myeloid cells in tumor growth, we hypothesized that CaMKK2 has a critical role in the tumor microenvironment, and tested this hyopothesis in murine models of hematological and solid cancer malignancies.

We found that CaMKK2 ablation in the host suppressed the growth of E.G7 murine lymphoma, Vk*Myc myeloma and E0771 mammary cancer. The selective ablation of CaMKK2 in myeloid cells was sufficient to restrain tumor growth, of which could be reversed by CD8 cell depletion. In the lymphoma microenvironment, ablating CaMKK2 generated less myeloid-derived suppressor cells (MDSCs) in vitro and in vivo. Mechanistically, CaMKK2 deficient dendritic cells showed higher Major Histocompatibility Class II (MHC II) and costimulatory factor expression, higher chemokine and IL-12 secretion when stimulated by LPS, and have higher potent in stimulating T-cell activation. AMPK, an anti-inflammatory kinase, was found as the relevant downstream target of CaMKK2 in dendritic cells. Treatment with CaMKK2 selective inhibitor STO-609 efficiently suppressed E.G7 and E0771 tumor growth, and reshaped the tumor microenvironment by attracting more immunogenic myeloid cells and infiltrated T cells.

In conclusion, we demonstrate that CaMKK2 expressed in myeloid cells is an important checkpoint in tumor microenvironment. Ablating CaMKK2 suppresses lymphoma growth by promoting myeloid cells development thereby decreasing MDSCs while enhancing the anti-tumor immune response. CaMKK2 inhibition is an innovative strategy for cancer therapy through reprogramming the tumor microenvironment.

Electrostatic Contributions to the Thermodynamics of Ribonuclease P Protein Folding

Electrostatic interactions are of fundamental importance in determining the structure and stability of macromolecules. For example, charge-charge interactions modulate the folding and binding of proteins and influence protein solubility. Electrostatic interactions are highly variable and can be both favorable and unfavorable. The ability to quantify these interactions is challenging but vital to understanding the detailed balance and major roles that they have in different proteins and biological processes. Measuring pKa values of ionizable groups provides a sensitive method for experimentally probing the electrostatic properties of a protein.

pKa values report the free energy of site-specific proton binding and provide a direct means of studying protein folding and pH-dependent stability. Using a combination of NMR, circular dichroism, and fluorescence spectroscopy along with singular value decomposition, we investigated the contributions of electrostatic interactions to the thermodynamic stability and folding of the protein subunit of Bacillus subtilis ribonuclease P, P protein. Taken together, the results suggest that unfavorable electrostatics alone do not account for the fact that P protein is intrinsically unfolded in the absence of ligand because the pKa differences observed between the folded and unfolded state are small. Presumably, multiple factors encoded in the P protein sequence account for its IUP property, which may play an important role in its function.