Department of Chemistry and Biochemistry

Candida parapsilosis is known to cause severe and persistent outbreaks in clinical settings. Patients infected with multidrug-resistant C. parapsilosis (MDR Cp) isolates were identified in a large Turkish hospital from 2017-2020. We subsequently identified three additional patients infected with MDR Cp isolates in 2022 from the same hospital and two echinocandin-resistant (ECR) isolates from a single patient in another hospital. The increasing number of MDR and ECR isolates contradicts the general principle that the severe fitness cost associated with these phenotypes could prevent their dominance in clinical settings. Here, we employed a multidimensional approach to systematically assess the fitness costs of MDR and ECR C. parapsilosis isolates. Whole-genome sequencing revealed a novel MDR genotype infecting two patients in 2022. Despite severe in vitro defects, the levels and tolerances of the biofilms of our ECR and MDR isolates were generally comparable to those of susceptible wild-type isolates. Surprisingly, the MDR and ECR isolates showed major alterations in their cell wall components, and some of the MDR isolates consistently displayed increased tolerance to the fungicidal activities of primary human neutrophils and were more immunoevasive during exposure to primary human macrophages. Our systemic infection mouse model showed that MDR and ECR C. parapsilosis isolates had comparable fungal burden in most organs relative to susceptible isolates. Overall, we observed a notable increase in the genotypic diversity and frequency of MDR isolates and identified MDR and ECR isolates potentially capable of causing persistent outbreaks in the future.

The standard redox potentials of metal nanoparticles are important for understanding their chemical properties. Traditionally, these redox potentials are measured by using voltammetry. Although voltammetry is fast and cost-effective, loading or landing the nanoparticles on electrodes alters their electrochemical properties, posing a challenge for accurately determining their intrinsic redox potentials. Here, a contactless method was developed utilizing chemical assays and the Nernst equation to measure the standard reduction potentials of gold nanoparticles in their colloidal state. To showcase the versatility and accuracy of this Nernstian approach, the reduction potentials were measured for a size range of 5.0-73 nm, revealing their scaling law and dependence on the nanoparticle surface energy.

Dynamin 1 (Dyn1) GTPase, a principal driver of membrane fission during synaptic endocytosis, self-assembles into short mechanoactive helices cleaving the necks of endocytic vesicles. While structural information about Dyn1 helix is abundant, little is known about the nanoscale dynamics of the helical scaffolding at the moment of fission, complicating mechanistic understanding of Dyn1 action. To address the role of the helix dynamics in fission, we used High-Speed Atomic Force Microscopy (HS-AFM) and fluorescence microscopy to track and compare the spatiotemporal characteristics of the helices formed by wild-type Dyn1 and its K44A mutant impaired in GTP hydrolysis on minimal lipid membrane templates. In the absence of nucleotide, membrane-bound WTDyn1 and K44ADyn1 self-assembled into tubular protein scaffolding of similar diameter encaging the lipid bilayer. In both cases, the GTP addition caused scaffold constriction coupled with formation of 20 to 30 nm nanogaps in the protein coverage. While both proteins reached scaffold diameters characteristic for membrane superconstriction causing fission, the fission was detected only with WTDyn1. We associated the fission activity with the dynamic evolution of the nanogaps: K44ADyn1 gaps were static, while WTDyn1 gaps actively evolved via repetitive nonaxisymmetric constriction-bending deformations caused by localized GTP hydrolysis. Modeling of the deformations implicated filament twist as an additional deformation mode which combines with superconstriction to facilitate membrane fission. Our results thus show that the dynamics of the Dyn1 helical scaffold goes beyond radial constriction and involves nonaxisymmetric deformations, where filament twist emerges as a critical driver of membrane fission.

Abstract: Time-dependent Hartree–Fock (TDHF) is one of the fundamental post-Hartree–Fock (HF) methods to describe excited states. In its Tamm-Dancoff form, equivalent to Configuration Interaction Singles, it is still widely used and particularly applicable to big molecules where more accurate methods may be unfeasibly expensive. However, it is rarely implemented in real space, mostly because of the expensive nature of the exact-exchange potential in real space. Compared to widely used Gaussian-type orbitals (GTO) basis sets, real space often offers easier implementation of equations and more systematic convergence of Rydberg states, as well as favorable scaling, effective domain parallelization, flexible boundary conditions, and ability to treat model systems. We implemented TDHF in the Octopus real-space code as a step toward linear-response hybrid time-dependent density-functional theory (TDDFT), other post-HF methods, and ensemble density-functional theory methods involving exact exchange. Calculation of HF’s non-local exact exchange is very expensive in real space. We overcome this limitation with Octopus’ implementation of Adaptively Compressed Exchange, and find the appropriate mixing scheme and starting point to complete the ground-state calculation in a practical amount of time, and thus enable TDHF. We compared our results to those from GTOs on a set of small molecules and confirmed close agreement of results, though with larger deviations than in the case of semi-local TDDFT. We find that convergence of TDHF demands a finer real-space grid than semi-local TDDFT. We also present the subtleties in benchmarking a real-space calculation against GTOs, relating to Rydberg and vacuum states.

Compared to ground-state electronic structure optimizations, accurate simulations of molecular real-time electron dynamics are usually much more difficult to perform. To simulate electron dynamics, the time-dependent density matrix renormalization group (TDDMRG) has been shown to offer an attractive compromise between accuracy and cost. However, many simulation parameters significantly affect the quality and efficiency of a TDDMRG simulation. So far, it is unclear whether common wisdom from ground-state DMRG carries over to the TDDMRG, and a guideline on how to choose these parameters is missing. Here, in order to establish such a guideline, we investigate the convergence behavior of the main TDDMRG simulation parameters, such as time integrator, the choice of orbitals, and the choice of matrix-product-state representation for complex-valued nonsinglet states. In addition, we propose a method to select orbitals that are tailored to optimize the dynamics. Lastly, we showcase the TDDMRG by applying it to charge migration ionization dynamics in furfural, where we reveal a rapid conversion from an ionized state with a σ character to one with a π character within less than a femtosecond.

The conformational ensemble and function of intrinsically disordered proteins (IDPs) are sensitive to their solution environment. The inherent malleability of disordered proteins, combined with the exposure of their residues, accounts for this sensitivity. One context in which IDPs play important roles that are concomitant with massive changes to the intracellular environment is during desiccation (extreme drying). The ability of organisms to survive desiccation has long been linked to the accumulation of high levels of cosolutes such as trehalose or sucrose as well as the enrichment of IDPs, such as late embryogenesis abundant (LEA) proteins or cytoplasmic abundant heat-soluble (CAHS) proteins. Despite knowing that IDPs play important roles and are co-enriched alongside endogenous, species-specific cosolutes during desiccation, little is known mechanistically about how IDP-cosolute interactions influence desiccation tolerance. Here, we test the notion that the protective function of desiccation-related IDPs is enhanced through conformational changes induced by endogenous cosolutes. We find that desiccation-related IDPs derived from four different organisms spanning two LEA protein families and the CAHS protein family synergize best with endogenous cosolutes during drying to promote desiccation protection. Yet the structural parameters of protective IDPs do not correlate with synergy for either CAHS or LEA proteins. We further demonstrate that for CAHS, but not LEA proteins, synergy is related to self-assembly and the formation of a gel. Our results suggest that functional synergy between IDPs and endogenous cosolutes is a convergent desiccation protection strategy seen among different IDP families and organisms, yet the mechanisms underlying this synergy differ between IDP families.

Background:Chlamydia trachomatis is the most prevalent bacterial sexually transmitted pathogen in humans worldwide. Since chlamydial infection is largely asymptomatic with the potential for serious complications, a preventative vaccine is likely the most viable long-term answer to this public health threat. Cell-free protein synthesis (CFPS) utilizes the cellular protein manufacturing machinery decoupled from the requirement for maintaining cellular viability, offering the potential for flexible, rapid, and decentralized production of recombinant protein vaccine antigens. Methods: Here, we use CFPS to produce the full-length putative chlamydial type three secretion system (T3SS) needle-tip protein, CT584, for evaluation as a vaccine antigen in mouse models. High-speed atomic force microscopy (HS-AFM) (RIBM, Tsukuba, Japan) imaging and computer simulations confirm that CFPS-produced CT584 retains a native-like structure prior to immunization. Female mice were primed with CT584 adjuvanted with CpG-1826 intranasally (i.n.) or CpG-1826 + Montanide ISA 720 intramuscularly (i.m.), followed four weeks later by an i.m. boost before respiratory challenge with 104 inclusion forming units (IFU) of Chlamydia muridarum. Results: Immunization with CT584 generated robust antibody responses but weak cell-mediated immunity and failed to protect against i.n. challenge as demonstrated by body weight loss, increased lung weights, and the presence of high numbers of IFUs in the lungs. Conclusion: While CT584 was not a protective vaccine candidate, the speed and flexibility with which CFPS can be used to produce other potential chlamydial antigens make it an attractive technique for antigen production.

Nanotube porins form transmembrane nanomaterial-derived scaffolds that mimic the geometry and functionality of biological membrane channels. We report synthesis, transport properties, and osmotic energy harvesting performance of another member of the nanotube porin family: boron nitride nanotube porins (BNNTPs). Cryo-transmission electron microscopy imaging, liposome transport assays, and DNA translocation experiments show that BNNTPs reconstitute into lipid membranes to form functional channels of ~2-nm diameter. Ion transport studies reveal ion conductance characteristics of individual BNNTPs, which show an unusual C^1/4 scaling with ion concentration and pronounced pH sensitivity. Reversal potential measurements indicate that BNNTPs have strong cation selectivity at neutral pH, attributable to the high negative charge on the channel. BNNTPs also deliver very large power density up to 12 kW/m² in the osmotic gradient transport experiments at neutral pH, surpassing that of other BNNT-based devices by two orders of magnitude under similar conditions. Our results suggest that BNNTPs are a promising platform for mass transport and osmotic power generation.

Abstract: Ensemble Density Functional Theory (EDFT) is a generalization of ground-state Density Functional Theory (GS DFT), which is based on an exact formal theory of finite collections of a system's ground and excited states. EDFT in various forms has been shown to improve the accuracy of calculated energy level differences in isolated model systems, atoms, and molecules, but it is not yet clear how EDFT could be used to calculate band gaps for periodic systems. We extend the application of EDFT toward periodic systems by estimating the thermodynamic limit with increasingly large finite one-dimensional "particle in a box" systems, which approach the uniform electron gas (UEG). Using ensemble-generalized Hartree and Local Spin Density Approximation (LSDA) exchange-correlation functionals, we find that corrections go to zero in the infinite limit, as expected for a metallic system. However, there is a correction to the effective mass, with results comparable to other calculations on 1D, 2D, and 3D UEGs, which indicates promise for non-trivial results from EDFT on periodic systems.

Spin-flip (SF) methods applied to excited-state approaches like the Bethe-Salpeter equation allow access to the excitation energies of open-shell systems, such as molecules and defects in solids. The eigenstates of these solutions, however, are generally not eigenstates of the spin operator Formula presented. Even for simple cases where the excitation vector is expected to be, for example, a triplet state, the value of ⟨Formula presented⟩ may be found to differ from 2.00; this difference is called ‘spin contamination’. The expectation values ⟨Formula presented⟩ must be computed for each excitation vector, to assist with the characterization of the particular excitation and to determine the amount of spin contamination of the state. Our aim is to provide for the first time in the SF methods literature a comprehensive resource on the derivation of the formulas for ⟨Formula presented⟩ as well as its computational implementation. After a brief discussion of the theory of the SF Bethe-Salpeter equation (BSE) and some examples further illustrating the need for calculating ⟨Formula presented⟩ , we present the derivation for the general equation for computing ⟨Formula presented⟩ with the eigenvectors from an SF-BSE calculation, how it is implemented in a Python script, and timing information on how this calculation scales with the size of the SF-BSE Hamiltonian.