The Effective Fragment Molecular Orbital Method for Fragments Connected by Covalent Bonds

We extend the effective fragment molecular orbital method (EFMO) into treating fragments connected by covalent bonds. The accuracy of EFMO is compared to FMO and conventional ab initio electronic structure methods for polypeptides including proteins. Errors in energy for RHF and MP2 are within 2 kcal/mol for neutral polypeptides and 6 kcal/mol for charged polypeptides similar to FMO but obtained two to five times faster. For proteins, the errors are also within a few kcal/mol of the FMO results. We developed both the RHF and MP2 gradient for EFMO. Compared to ab initio, the EFMO optimized structures had an RMSD of 0.40 and 0.44 {\AA} for RHF and MP2, respectively.Comment: Revised manuscrip

Clinical Features and Long-Term Follow-Up of Quasi-Moyamoya Disease in Children.

Background: Inherited or acquired disorders and conditions may present in conjunction with moyamoya disease. This condition is known as quasi-moyamoya disease. Methods: A retrospective review of 69 moyamoya disease patients treated for the past 20 years identified 7 patients with quasi-moyamoya disease and 5 of them were pediatric patients. Results: The mean age at initial diagnosis was 6.4 years (range 5-9). Associated disorders were: craniosynostosis, dwarfism with coarctation of aorta, Proteus syndrome, and cranial irradiation for brain tumor. Their clinical type included cerebral ischemia in 3 patients, cerebral bleeding with ischemia in 1 and epilepsy in 1. The 3 patients with cerebral ischemia underwent bypass surgery and their ischemia was improved. One patient died of brain tumor recurrence and the activities of daily living in the remaining patients were affected by mental retardation. Conclusions: The clinical course and radiological finding of quasi-moyamoya disease are diverse because of associated disorders, distinguishing definite moyamoya disease

HER2 testing in breast cancer: Opportunities and challenges

Human epidermal growth factor receptor 2 (HER2) is overexpressed in 15-25% of breast cancers, usually as a result of HER2 gene amplification. Positive HER2 status is considered to be an adverse prognostic factor. Recognition of the role of HER2 in breast cancer growth has led to the development of anti-HER2 directed therapy, with the humanized monoclonal antibody trastuzumab (Herceptin (R)) having been approved for the therapy of HER2-positive metastatic breast cancer. Clinical studies have further suggested that HER2 status can provide important information regarding success or failure of certain hormonal therapies or chemotherapies. As a result of these developments, there has been increasing demand to perform HER2 testing on current and archived breast cancer specimens. This article reviews the molecular background of HER2 function, activation and inhibition as well as current opinions concerning its role in chemosensitivity and interaction with estrogen receptor biology. The different tissue-based assays used to detect HER2 amplification and overexpression are discussed with respect to their advantages and disadvantages, when to test (at initial diagnosis or pre-treatment), where to test (locally or centralized) and the need for quality assurance to ensure accurate and valid testing results

A model for transition of 5 '-nuclease domain of DNA polymerase I from inert to active modes

Bacteria contain DNA polymerase I (PolI), a single polypeptide chain consisting of similar to 930 residues, possessing DNA-dependent DNA polymerase, 3'-5' proofreading and 5'-3' exonuclease (also known as flap endonuclease) activities. PolI is particularly important in the processing of Okazaki fragments generated during lagging strand replication and must ultimately produce a double-stranded substrate with a nick suitable for DNA ligase to seal. PolI's activities must be highly coordinated both temporally and spatially otherwise uncontrolled 5'-nuclease activity could attack a nick and produce extended gaps leading to potentially lethal double-strand breaks. To investigate the mechanism of how PolI efficiently produces these nicks, we present theoretical studies on the dynamics of two possible scenarios or models. In one the flap DNA substrate can transit from the polymerase active site to the 5'-nuclease active site, with the relative position of the two active sites being kept fixed; while the other is that the 5'-nuclease domain can transit from the inactive mode, with the 5'-nuclease active site distant from the cleavage site on the DNA substrate, to the active mode, where the active site and substrate cleavage site are juxtaposed. The theoretical results based on the former scenario are inconsistent with the available experimental data that indicated that the majority of 5'-nucleolytic processing events are carried out by the same PolI molecule that has just extended the upstream primer terminus. By contrast, the theoretical results on the latter model, which is constructed based on available structural studies, are consistent with the experimental data. We thus conclude that the latter model rather than the former one is reasonable to describe the cooperation of the PolI's polymerase and 5'-3' exonuclease activities. Moreover, predicted results for the latter model are presented

Order-of-magnitude speedup for steady states and traveling waves via Stokes preconditioning in Channelflow and Openpipeflow

Steady states and traveling waves play a fundamental role in understanding hydrodynamic problems. Even when unstable, these states provide the bifurcation-theoretic explanation for the origin of the observed states. In turbulent wall-bounded shear flows, these states have been hypothesized to be saddle points organizing the trajectories within a chaotic attractor. These states must be computed with Newton's method or one of its generalizations, since time-integration cannot converge to unstable equilibria. The bottleneck is the solution of linear systems involving the Jacobian of the Navier-Stokes or Boussinesq equations. Originally such computations were carried out by constructing and directly inverting the Jacobian, but this is unfeasible for the matrices arising from three-dimensional hydrodynamic configurations in large domains. A popular method is to seek states that are invariant under numerical time integration. Surprisingly, equilibria may also be found by seeking flows that are invariant under a single very large Backwards-Euler Forwards-Euler timestep. We show that this method, called Stokes preconditioning, is 10 to 50 times faster at computing steady states in plane Couette flow and traveling waves in pipe flow. Moreover, it can be carried out using Channelflow (by Gibson) and Openpipeflow (by Willis) without any changes to these popular spectral codes. We explain the convergence rate as a function of the integration period and Reynolds number by computing the full spectra of the operators corresponding to the Jacobians of both methods.Comment: in Computational Modelling of Bifurcations and Instabilities in Fluid Dynamics, ed. Alexander Gelfgat (Springer, 2018

Biological Stoichiometry in Human Cancer

A growing tumor in the body can be considered a complex ecological and evolutionary system. A new eco-evolutionary hypothesis (the "Growth Rate Hypothesis", GRH) proposes that tumors have elevated phosphorus (P) demands due to increased allocation to P-rich nucleic acids, especially ribosomal RNA, to meet the protein synthesis demands of accelerated proliferation.We determined the elemental (C, N, P) and nucleic acid contents of paired malignant and normal tissues from colon, lung, liver, or kidney for 121 patients. Consistent with the GRH, lung and colon tumors were significantly higher (by approximately two-fold) in P content (fraction of dry weight) and RNA content and lower in nitrogen (N):P ratio than paired normal tissue, and P in RNA contributed a significantly larger fraction of total biomass P in malignant relative to normal tissues. Furthermore, patient-specific differences for %P between malignant and normal tissues were positively correlated with such differences for %RNA, both for the overall data and within three of the four organ sites. However, significant differences in %P and %RNA between malignant and normal tissues were not seen in liver and kidney and, overall, RNA contributed only approximately 11% of total tissue P content.Data for lung and colon tumors provide support for the GRH in human cancer. The two-fold amplification of P content in colon and lung tumors may set the stage for potential P-limitation of their proliferation, as such differences often do for rapidly growing biota in ecosystems. However, data for kidney and liver do not support the GRH. To account for these conflicting observations, we suggest that local environments in some organs select for neoplastic cells bearing mutations increasing cell division rate ("r-selected," as in colon and lung) while conditions elsewhere may select for reduced mortality rate ("K-selected," as in liver and kidney)