Search

We present efficient analytical gradients of property-based diabatic states and couplings using a Lagrangian formalism. Unlike previous formulations, the method achieves a computational scaling that is independent of the number of adiabatic states used to construct the diabats. The approach is generalizable to other property-based diabatization schemes and electronic structure methods as long as analytical energy gradients are available and integral derivatives with the property operator can be formed. We also introduce a scheme to phase and reorder diabats to ensure their continuity between molecular configurations. We demonstrate this for the specific case of Boys diabatic states obtained from state-averaged complete active space self-consistent field electronic structure calculations with GPU acceleration in the TeraChem package. The method is used to test the Condon approximation for the hole transfer in an explicitly solvated model DNA oligomer.

The X-ray absorption spectrum (XAS) of the hydrated electron (e−(aq)) has been simulated using time- dependent density functional theory with a quantum mechanics/molecular mechanics description. A unique XAS peak at 533 eV is observed with an energy and intensity in quantitative agreement with recent time-resolved experiments, allowing its definitive assignment as arising from water O1s transitions to the singly occupied molecular orbital (SOMO) in which the excess electron resides. The transitions acquire oscillator strength due to the SOMO comprising an admixture of a cavity-localized orbital and water 4a1 and 2b2 antibonding orbitals. The mixing of antibonding orbitals has implications for the strength of couplings between e−(aq) and intramolecular modes of water.

This article considers a range of early twentieth-century projects in which rural space was made subject to an ensemble of institutional forms and practices grounded in emergent urban paradigms. The projects in rural reconstruction I consider here seized on village space and the minds and souls of villagers as their terrain of operation. Rural reconstruction entailed the production of intimate knowledge of the rural population, the development of affective modes of engagement with them, and investments of governmental power not only in state institutions but in more abstract concepts like the ‘rural community’. Many projects were largely propagandistic; their concrete effects were minimal. But rural reconstruction brought urban logics, objects, and institutional forms directly into village milieus in ways that continue to shape how India’s urban and rural areas are conceptualised and operate in relation to one another. If urbanism in India today is deeply enmeshed in agrarian processes, then this article attempts to reverse the gaze by asking how expert knowledge directed towards managing rural change in twentieth-century India depended crucially on emergent urban paradigms.

A polarizable embedding scheme is presented that accounts for differential solvation of ground and excited states in QM/MM simulations. The polarization and dispersion interactions between the quantum-mechanical (QM) and molecular-mechanical (MM) regions are described by the direct reaction field (DRF) Hamiltonian, while the Pauli repulsion between explicitly treated QM electrons and the implicit electron density around MM atoms is modeled with effective core potentials. A single Hamiltonian is used for all electronic states, so that Born-Oppenheimer states belonging to the same geometry are orthogonal and state crossings are well-defined. The method is implemented in TeraChem, where it is combined with multiple electronic structure methods, including Hartree-Fock, Configuration Interaction Singles, and Complete Active Space Self-Consistent Field. In contrast with older implementations of the DRF method, integrals of the polarization operators are evaluated exactly. Expressions for ingredients needed to construct analytical gradients and non-adiabatic coupling vectors are derived and tested by optimizing a conical intersection between two excited states in the presence of a polarizable solvent shell. The method is applied to estimating the solvent shifts of absorption energies of a series of donor-acceptor dyes having low-lying charge transfer states. Even for a non-polar solvent such as n-hexane the inclusion of its static polarizability leads to non-negligible shifts that improve the agreement with the positions of experimental absorption maxima measured in solution.

A hybrid QM/MM study with the integral-exact direct reaction field (IEDRF) polarizable embedding scheme developed in the companion paper is performed on the purple bacteria reaction center of Rhodobacter sphaeroides in order to investigate how polarization induced in the protein matrix helps to stabilize the primary charge-transfer state relative to the exciton states of the special pair. The protein environment is represented by point charges and induced dipoles that are coupled to the QM region via the integral-exact direct reaction field Hamiltonian, which can simultaneously describe differential solvation of multiple electronic states of different polarities. Treating the special pair, $P=P_L P_M$, and one of the bacteriopheophytins ($H_L$ or $H_M$) quantum-mechanically, we compare excitation energies computed at the $\omega$PBEh time-dependent density functional theory level for charge-transfer states along the active and inactive branches. Thermal fluctuations on the electronic-state energies are included by extracting snapshots from a molecular dynamics trajectory. With IEDRF embedding, the reaction field induced in the protein matrix stabilizes long-range charge-transfer (CT) states by over 1.0 eV, shifting them below the exciton states. The protein environment favors charge separation along the active $L$ branch. With only electrostatic embedding, the CT states are found $>$0.5 eV above the exciton states, and asymmetry between the branches is diminished. The polarization in the protein is largely dictated by the secondary structure, with induced dipoles pointing along the axes of $\alpha$-helices. The relaxation of dipoles on the excited state provides similar stabilization of both branches, however the ground state polarization, captured by IEDRF, screens $P^+ H_L^-$ more than $P^+ H_M^-$ and is thus key to directing the CT towards the active branch. The different ground-state dielectric environment of the two branches is confirmed by computing relative local fields at key carbonyl probes. With IEDRF embedding, we find the local fields to be in near quantitative agreement with the interpretations of recent vibrational Stark-effect experiments, while the agreement is diminished for electrostatically embedded QM/MM calculations. These results caution against QM/MM calculations with electrostatic embedding alone, because CT states in a rigid, dense protein matrix are mostly stabilized by the fast electronic polarization of the surrounding medium.

The photoisomerization of isolated ethylene (ethene) was observed in real time from the Franck-Condon region in the

Background:Carbapenem-resistant Enterobacterales (CRE) are common causes of healthcare-associated infections and are often multidrug resistant with limited therapeutic options. Additionally, CRE can spread within and between healthcare facilities, amplifying potential harms.Objective:To better understand the burden, risk factors, and source of acquisition of carbapenemase genes in clinical Escherichia coli and Klebsiella spp isolates from patients in Washington to guide prevention efforts.Design:Multicenter prospective surveillance study.Methods:Escherichia coli and Klebsiella spp isolates meeting the Washington state CRE surveillance case definition were solicited from clinical laboratories and tested at Washington Public Health Laboratories using polymerase chain reaction (PCR) for the 5 most common carbapenemase genes: blaKPC, blaNDM, blaIMP, blaVIM, and blaOXA-48. Case patients positive by PCR were investigated by the public health department.Results:From October 2012 through December 2017, 363 carbapenem-resistant E. coli and Klebsiella spp isolates were tested. Overall, 45 of 115 carbapenem-resistant K. pneumoniae (39%), 1 of 8 K. oxytoca (12.5%), and 28 of 239 carbapenem-resistant E. coli (11.7%) were carbapenemase positive. Of 74 carbapenemase-positive isolates, blaKPC was most common (47%), followed by blaNDM (30%), blaOXA-48 (22%), and blaIMP (1%). Although all cases had healthcare exposure, blaKPC acquisition was associated with US health care, whereas non-blaKPC acquisition was associated with international health care or travel.Conclusions:We report that blaKPC, the most prevalent carbapenemase in the United States, accounts for nearly half of carbapenemase cases in Washington state and that most KPC-cases are likely acquired through in-state health care.

Embedding theory is a powerful computational chemistry approach to exploring the electronic structure and dynamics of complex systems, with QM/MM being the prime example. A challenge arises when trying to apply embedding methodology to systems with diffusible particles, e.g. solvents, if some of them must be included in the QM region, for example in the description of solvent-supported electronic states or reactions involving proton transfer or charge-transfer-to-solvent: without a special treatment, inter-diffusion of QM and MM particles will lead eventually to a loss of QM/MM separation. We have developed a new method called Flexible Boundary Layer using Exchange (FlexiBLE) that solves the problem by adding a biasing potential to the system that closely maintains QM/MM separation. The method rigorously preserves ensemble averages by leveraging their invariance to exchange of identical particles. With a careful choice of the biasing potential, and the use of a tree algorithm to include only important QM and MM exchanges, we find the method has an MM-forcefield-like computational cost and thus adds negligible overhead to a QM/MM simulation. Furthermore, we show that molecular dynamics with the FlexiBLE bias conserves total energy and remarkably, dynamical quantities in the QM region are unaffected by the applied bias. FlexiBLE thus widens the range of chemistry that can be studied with embedding theory.

The FlexiBLE embedding method introduced in the preceding companion paper [Z. Shen and W. J. Glover, J. Chem. Phys. X, X (2021)] is applied to explore the structure and dynamics of the aqueous solvated electron at an all-electron density functional theory QM/MM level. Compared to a one-electron mixed quantum/classical description, we find the dynamics of the many-electron model of the hydrated electron exhibits enhanced coupling to water OH stretch modes. Natural Bond Orbital analysis reveals this coupling is due to significant population of water OH σ* orbitals, reaching 20%. Based on this, we develop a minimal frontier orbital picture of the hydrated electron involving a cavity orbital and important coupling to 4-5 coordinating OH σ* orbitals. Implications for the interpretation of the spectroscopy of this interesting species are discussed.

There has been a great deal of recent controversy over the structure of the hydrated electron and whether it occupies a cavity or contains a significant number of interior waters (noncavity). The questions we address in this work are, from a free energy perspective, how different are these proposed structures? Do the different structures all lie along a single continuum, or are there significant differences (i.e., free energy barriers) between them? To address these questions, we have performed a series of one-electron calculations using umbrella sampling with quantum biased molecular dynamics along a coordinate that directly reflects the number of water molecules in the hydrated electron's interior. We verify that a standard cavity model of the hydrated electron behaves essentially as a hard sphere: the model is dominated by repulsion at short range such that water is expelled from a local volume around the electron, leading to a water solvation shell like that of a pseudohalide ion. The repulsion is much larger than thermal energies near room temperature, explaining why such models exhibit properties with little temperature dependence. On the other hand, our calculations reveal that a noncavity model is highly fluxional, meaning that thermal motions cause the number of interior waters to fluctuate from effectively zero (i.e., a cavity-type electron) to potentially above the bulk water density. The energetic contributions in the noncavity model are still repulsive in the sense that they favor cavity formation, so the fluctuations in structure are driven largely by entropy: the entropic cost for expelling water from a region of space is large enough that some water is still driven into the electron's interior. As the temperature is lowered and entropy becomes less important, the noncavity electron's structure is predicted to become more cavity-like, consistent with the observed temperature dependence of the hydrated electron's properties. Thus, we argue that although the specific noncavity model we study overestimates the preponderance of fluctuations involving interior water molecules, with appropriate refinements to correctly capture the true average number of interior waters and molar solvation volume, a fluxional model likely makes the most sense for understanding the various experimental properties of the hydrated electron.

Summary Two-dimensional electronic spectroscopy (2D ES) probes the energies of chromophores and their coupling, and therefore, the technique is now widely used for studying excitation-energy-transfer mechanisms. The ubiquity and importance of Forster resonance energy transfer (FRET) demand a thorough study of its signatures in 2D ES. Here, using principles from noisy-light spectroscopy to account for the spontaneous transfer event, we describe a model of the signals that arise from FRET in 2D ES. This model yields three theoretical results regarding the energy-transfer peak, which are consistent with our laboratory measurements. A second unanticipated finding is that the solvation processes that give rise to the dynamic Stokes shift of the acceptor signal also induce a relatively slower redshift of the energy-transfer signal, revealing a fundamental difference in solvation relaxation between photoexcited and FRET-excited acceptor molecules. The signatures derived from the model serve as a benchmark for ongoing work on energy-transfer mechanisms. Here, Turner and colleagues develop a noisy-light-inspired model of FRET in 2D ES to account for the vibrations that activate it. Crucially, they reveal several unreported features that distinguish FRET from other signals. These signatures are fundamental and immediately applicable to ongoing work on photosynthetic complexes, solar nanomaterials, and other renewable energy focal points. Future interdisciplinary work ranging from studies of energy transfer in organisms to medically applicable bio-imaging methodologies will take advantage of these spectral and temporal signatures.

We introduce Active Orbital Preservation for Multiconfigurational Self-Consistent Field (AOP-MCSCF), an automated approach to improving the consistency of active space orbitals over multiple molecular configurations. Our approach is based on maximum overlap with a reference set of active space orbitals taken from a single geometry of a chromophore in the gas phase and can be used to automatically preserve the appropriate orbitals of the chromophore across multiple thermally sampled configurations, even when the chromophore is solvated by quantum-mechanically treated water molecules. In particular, using the singular value decomposition of a Molecular Orbital (MO) overlap matrix between the system and reference, we rotate the MOs of the system to align with the reference active space orbitals and use the resulting rotated orbitals as an initial guess to a MCSCF calculation. We demonstrate the approach on aqueous p-hydroxybenzylidene-imidazolinone (HBI) and find that AOP-MCSCF converges to the “correct” orbitals for over 90% of 3000 thermally sampled configurations. In addition, we compute the linear absorption spectrum and find excellent agreement with new experimental measurements up to 5.4 eV (230 nm). We show that electrostatic contributions to the solvation energy of HBI largely explain the observed state-dependent solvatochromism.

We explore the convergence of the many-body expansion for a prototypical hole-transfer reaction between Zn(0) and Zn(I) in a condensed-phase environment. Poor convergence of state energies is seen when the adiabatic representation is used, which can be understood from the fragment single-point calculations at low orders of the many-body expansion incorrectly localizing charges compared to the full system, thus leading to qualitative errors in the electronic structure of the adiabatic states between fragments. Using a charge-localized representation of the electronic Hamiltonian, we introduce a diabatic many-body expansion method with quantitative accuracy for ground- and excited-state potential energy surfaces of a charge-transfer reaction. Combining with a multiconfigurational self-consistent field affords a fragmentation approach that scales quadratically with system size while retaining chemical accuracy (

Carbapenem-resistant Enterobacterales (CRE) are common causes of healthcare-associated infections and are often multidrug resistant with limited therapeutic options. Additionally, CRE can spread within and between healthcare facilities, amplifying potential harms.To better understand the burden, risk factors, and source of acquisition of carbapenemase genes in clinical Escherichia coli and Klebsiella spp isolates from patients in Washington to guide prevention efforts.Multicenter prospective surveillance study.Escherichia coli and Klebsiella spp isolates meeting the Washington state CRE surveillance case definition were solicited from clinical laboratories and tested at Washington Public Health Laboratories using polymerase chain reaction (PCR) for the 5 most common carbapenemase genes: blaKPC, blaNDM, blaIMP, blaVIM, and blaOXA-48. Case patients positive by PCR were investigated by the public health department.From October 2012 through December 2017, 363 carbapenem-resistant E. coli and Klebsiella spp isolates were tested. Overall, 45 of 115 carbapenem-resistant K. pneumoniae (39%), 1 of 8 K. oxytoca (12.5%), and 28 of 239 carbapenem-resistant E. coli (11.7%) were carbapenemase positive. Of 74 carbapenemase-positive isolates, blaKPC was most common (47%), followed by blaNDM (30%), blaOXA-48 (22%), and blaIMP (1%). Although all cases had healthcare exposure, blaKPC acquisition was associated with US health care, whereas non-blaKPC acquisition was associated with international health care or travel.We report that blaKPC, the most prevalent carbapenemase in the United States, accounts for nearly half of carbapenemase cases in Washington state and that most KPC-cases are likely acquired through in-state health care.

Recent interest in polarizable embedding methods for electronic excited states has so far been focused on optical absorption and emission spectra calculations. To explore the suitability of these methods for excited-state reactions, we constructed a simple molecular system with an electronic crossing coupled to a polarizable species: the triatomic LiFBe. We found that current polarizable QM/MM methods inadequately describe the potential energy surfaces in this system, particularly close to the electronic crossing, so we developed a new polarizable embedding method called dynamically weighted polarizable QM/MM. The new method reproduces the potential energy surfaces of LiFBe from full-system multireference configuration interaction singles and doubles calculations with near-quantitative accuracy.

A fundamental understanding of the excited-state properties of molecular crystals is of central importance for their optoelectronics applications. In this study, we developed the electrostatically embedded generalized molecular fractionation (EE-GMF) method for the quantitative characterization of the excited-state properties of locally excited molecular clusters. The accuracy of the EE-GMF method is systematically assessed for oligoacene crystals. Our result demonstrates that the EE-GMF method is capable of providing the lowest vertical singlet (S

The hydrated electron, e–(aq), has often served as a model system to understand the influence of condensed-phase environments on electronic structure and dynamics. Despite over 50 years of study, however, the basic structure of e–(aq) is still the subject of controversy. In particular, the structure of e–(aq) was long assumed to be an electron localized within a solvent cavity, in a manner similar to halide solvation. Recently, however, we suggested that e–(aq) occupies a region of enhanced water density with little or no discernible cavity. The potential we developed was only subtly different from those that give rise to a cavity solvation motif, which suggests that the driving forces for noncavity solvation involve subtle electron-water attractive interactions at close distances. This leads to the question of how dispersion interactions are treated in simulations of the hydrated electron. Most dispersion potentials are ad hoc or are not designed to account for the type of close-contact electron-water ov...

We report a 26-year-old man with Psychrobacter sanguinis cellulitis of a wound sustained during ocean fishing in Washington, USA, in 2017. Psychrobacter spp. are opportunistic pathogens found in a wide range of environments. Clinicians should be aware of Psychrobacter spp. and perform 16S rRNA sequencing if this pathogen is suspected.

The ultrafast excited state dynamics of the smallest polyene, trans-1,3-butadiene, were studied by femtosecond time-resolved photoelectron-photoion coincidence (TRPEPICO) spectroscopy. The evolution of the excited state wavepacket, created by pumping the bright 1Bu (ππ*) electronic state at its origin of 216 nm, is projected via one- and two-photon ionization at 267 nm onto several ionization continua. The results are interpreted in terms of Koopmans’ correlations and Franck-Condon factors for the excited and cationic states involved. The known predissociative character of the cation excited states is utilized to assign photoelectron bands to specific continua using TRPEPICO spectroscopy. This permits us to report the direct observation of the famously elusive S1(21Ag) dark electronic state during the internal conversion of trans 1,3-butadiene. Our phenomenological analysis permits the spectroscopic determination of several important time constants. We report the overall decay lifetimes of the 11Bu and 21Ag states and observe the re-appearance of the hot ground state molecule. We argue that the apparent dephasing time of the S2(11Bu) state, which leads to the extreme breadth of the absorption spectrum, is principally due to large amplitude torsional motion on the 1Bu surface in conjunction with strong non-adiabatic couplings via conical intersections, whereupon nuclear wavepacket revivals to the initial Franck-Condon region become effectively impossible. In Paper II [W. J. Glover et al., J. Chem. Phys. 148, 164303 (2018)], ab initio multiple spawning is used for on-the-fly computations of the excited state non-adiabatic wavepacket dynamics and their associated TRPEPICO observables, allowing for direct comparisons of experiment with theory.

The excited state non-adiabatic dynamics of the smallest polyene, trans 1,3-butadiene (BD), has long been the subject of controversy due to its strong coupling, ultrafast time scales and the difficulties that theory faces in describing the relevant electronic states in a balanced fashion. Here we apply Ab Initio Multiple Spawning (AIMS) using state-averaged complete active space multistate second order perturbation theory [SA-3-CAS(4/4)-MSPT2] which describes both static and dynamic electron correlation effects, providing a balanced description of both the initially prepared bright 11Bu (ππ*) state and non-adiabatically coupled dark 21Ag state of BD. Importantly, AIMS allows for on-the-fly calculations of experimental observables. We validate our approach by directly simulating the time resolved photoelectron-photoion coincidence spectroscopy results presented in Paper I [A. E. Boguslavskiy et al., J. Chem. Phys. 148, 164302 (2018)], demonstrating excellent agreement with experiment. Our simulations reveal that the initial excitation to the 11Bu state rapidly evolves via wavepacket dynamics that follow both bright- and dark-state pathways as well as mixtures of these. In order to test the sensitivity of the AIMS results to the relative ordering of states, we considered two hypothetical scenarios biased toward either the bright 1Bu or the dark 21Ag state. In contrast with AIMS/SA-3-CAS(4/4)-MSPT2 simulations, neither of these scenarios yields favorable agreement with experiment. Thus, we conclude that the excited state non-adiabatic dynamics in BD involves both of these ultrafast pathways.

Quantitative DNA detection of cytomegalovirus (CMV) and BK virus (BKV) is critical in the management of transplant patients. Quantitative laboratory-developed procedures for CMV and BKV have been described in which much of the processing is automated, resulting in rapid, reproducible, and high-throughput testing of transplant patients. To increase the efficiency of such assays, the performance and stability of four commercial preassembled frozen fast qPCR master mixes (Roche FastStart Universal Probe Master Mix with Rox, Bio-Rad SsoFast Probes Supermix with Rox, Life Technologies TaqMan FastAdvanced Master Mix, and Life Technologies Fast Universal PCR Master Mix), in combination with in-house designed primers and probes, was evaluated using controls and standards from standard CMV and BK assays. A subsequent parallel evaluation using patient samples was performed comparing the performance of freshly prepared assay mixes versus aliquoted frozen master mixes made with two of the fast qPCR mixes (Life Technologies TaqMan FastAdvanced Master Mix, and Bio-Rad SsoFast Probes Supermix with Rox), chosen based on their performance and compatibility with existing PCR cycling conditions. The data demonstrate that the frozen master mixes retain excellent performance over a period of at least 10 weeks. During the parallel testing using clinical specimens, no difference in quantitative results was observed between the preassembled frozen master mixes and freshly prepared master mixes. Preassembled fast real-time qPCR frozen master mixes perform well and represent an additional strategy laboratories can implement to reduce assay preparation times, and to minimize technical errors and effort necessary to perform clinical PCR.

Wohlfahrtiimonas species bacteria were isolated from the bloodstream of a patient with septicemia and wound myiasis. Environmental investigations identified a Wohlfahrtiimonas sp. among insects in the Americas and in a previously undescribed vector, the green bottle fly (Lucilia sericata). The isolates possibly represent a new species within the genus Wohlfahrtiimonas.

We use one-electron non-adiabatic mixed quantum/classical simulations to explore the temperature dependence of both the ground-state structure and the excited-state relaxation dynamics of the hydrated electron. We compare the results for both the traditional cavity picture and a more recent non-cavity model of the hydrated electron and make definite predictions for distinguishing between the different possible structural models in future experiments. We find that the traditional cavity model shows no temperature-dependent change in structure at constant density, leading to a predicted resonance Raman spectrum that is essentially temperature-independent. In contrast, the non-cavity model predicts a blue-shift in the hydrated electron's resonance Raman O-H stretch with increasing temperature. The lack of a temperature-dependent ground-state structural change of the cavity model also leads to a prediction of little change with temperature of both the excited-state lifetime and hot ground-state cooling time of the hydrated electron following photoexcitation. This is in sharp contrast to the predictions of the non-cavity model, where both the excited-state lifetime and hot ground-state cooling time are expected to decrease significantly with increasing temperature. These simulation-based predictions should be directly testable by the results of future time-resolved photoelectron spectroscopy experiments. Finally, the temperature-dependent differences in predicted excited-state lifetime and hot ground-state cooling time of the two models also lead to different predicted pump-probe transient absorption spectroscopy of the hydrated electron as a function of temperature. We perform such experiments and describe them in Paper II [E. P. Farr et al., J. Chem. Phys. 147, 074504 (2017)], and find changes in the excited-state lifetime and hot ground-state cooling time with temperature that match well with the predictions of the non-cavity model. In particular, the experiments reveal stimulated emission from the excited state with an amplitude and lifetime that decreases with increasing temperature, a result in contrast to the lack of stimulated emission predicted by the cavity model but in good agreement with the non-cavity model. Overall, until ab initio calculations describing the non-adiabatic excited-state dynamics of an excess electron with hundreds of water molecules at a variety of temperatures become computationally feasible, the simulations presented here provide a definitive route for connecting the predictions of cavity and non-cavity models of the hydrated electron with future experiments.

We recently introduced a Dynamically Weighted Complete Active Space Self-Consistent Field (DW-CASSCF) electronic structure for excited-state dynamics. In this Communication, we reformulate analytical gradients at this level of theory using a Lagrangian approach, thereby reducing the required number of coupled-perturbed CASSCF calculations to one per state gradient. In addition, we derive and implement derivative couplings at the DW-CASSCF level for the first time. We demonstrate the new formulation of DW-CASSCF gradients by optimizing a conical intersection for the p-hydroxybenzylidene-imidazolinone anion, the green fluorescent protein chromophore, to shed light on its observed radiationless decay dynamics in the ultraviolet region.

Human health risks associated with exposure to algal and cyanobacterial toxins (phycotoxins) have been largely concerned with aquatic habitats. People inhabiting desert environments may be exposed to phycotoxins present in terrestrial environments, where cyanobacterial crusts dominate. Seafood comprises a significant portion of the human diet in desert environments proximal to an ocean or sea. Consequently, in addition to terrestrial exposure to cyanotoxins, the potential exists that seafood may be an important exposure route for cyanotoxins in desert regions. Understanding the possible risk of exposure from seafood will help create cyanotoxin health guidelines for people living in environments that rely on seafood. Commonly-consumed local seafood products destined for human consumption were purchased from a fish market in Doha, Qatar. Organs were excised, extracted, and analyzed for the neurotoxic amino acid β-N-methylamino-L-alanine (BMAA) and the isomers 2,4-diaminobutyric acid (DAB) and N-2(aminoethyl)glycine (AEG). The presence and concentration of neurotoxic amino acids were investigated in organisms from various trophic levels to examine the potential for biomagnification. Although BMAA and isomers were detected in marine microbial mats, as well as in marine plankton net trawls associated with diatoms and dinoflagellates, in seafood, only AEG and DAB were present at low concentrations in various trophic levels. The findings of this study suggest that exposure to neurotoxic amino acids through seafood in the Arabian Gulf may be minor, yet the presence of BMAA in phytoplankton confirms the need for further monitoring of marine waters and seafood to protect human health.